Numerical integration of the field equations in bimetric relativity is necessary to obtain solutions describing realistic systems. Thus, it is crucial to recast the equations as a well-posed problem. In general relativity, under certain assumptions, the covariant BSSN formulation is a strongly hyperbolic formulation of the Einstein equations, hence its Cauchy problem is well-posed. In this paper, we establish the covariant BSSN formulation of the bimetric field equations. It shares many features with the corresponding formulation in general relativity, but there are a few fundamental differences between them. Some of these differences depend on the gauge choice and alter the hyperbolic structure of the system of partial differential equations compared to general relativity. Accordingly, the strong hyperbolicity of the system cannot be claimed yet, under the same assumptions as in general relativity. In the paper, we stress the differences compared with general relativity and state the main issues that should be tackled next, to draw a road map towards numerical bimetric relativity.

We study general properties of static and spherically symmetric bidiagonal black holes in Hassan-Rosen bimetric theory. In particular, we explore the behaviour of the black hole solutions both at the common Killing horizon and at the large radii. The former study leads to a new classification for black holes within the bidiagonal ansatz. The latter study shows that, among the great variety of the black hole solutions, the only solutions converging to Minkowski, Anti-de Sitter and de Sitter spacetimes at large radii are those of General Relativity, i.e., the Schwarzschild, Schwarzschild-Anti-de Sitter and Schwarzschild-de Sitter solutions.