1. Introduction
Binary black holes are considered to be the foremost vital sources for the emerging field of gravitational-wave astrophysics. Multiple black hole systems, above obvious curiosities, may also serve as natural gravitational waves sources. Investigation of the dynamics of these systems starts with a careful initialization. This may be done by applying the elliptic method [
1,
2] (see also [
3,
4]) or either of the evolutionary form of the constraints introduced in [
5].
In [
6,
7], the parabolic-hyperbolic formulation of Hamiltonian and momentum constraints, along with superposing individual Kerr-Schild black holes, was applied to construct initial data for binary black hole configurations. The only technical restriction was that each of the initial speeds is parallel to and each of the spin vectors is orthogonal to a plane of some background Euclidean space, respectively. For this class, the existence and uniqueness of (at least)
solutions to the parabolic-hyperbolic form of the constraints is outlined in [
6,
7]. This new method has also been successfully applied to determine initial data numerically for individual and binary black hole systems [
8,
9]. Notably, the very same construction can also be used to initialize multiple black hole systems. There are no restrictions on the masses, speeds, spins, and distances of the individual black holes; thereby, this set contains many physically realistic initial data configurations.
This paper will focus on the Arnowitt-Deser-Misner (ADM) charges of multiple black hole systems. Accordingly, this involves superposing individual Kerr-Schild black holes, each located momentarily on a plane in some background Euclidean space with speeds parallel to and spins orthogonal to the distinguished plane. By adopting constructive elements of the proposal in [
6,
7], we shall choose the free data to the initial-boundary value problem, derived from the parabolic-hyperbolic form of the constraints, using superposed Kerr-Schild black holes. Clearly, if the individual Kerr-Schild black holes are widely separated, the initial data, satisfying the constraints, will only slightly differ from the superposed data induced on a
time-slice. Therefore, it is highly plausible that the craved global solutions exist. Note, however, that the verification of the primary result of the present paper refers only to the specific choice of the free data, and it does not require detailed knowledge of solutions. Therefore, no attempt will be made to deal with the global existence and uniqueness of solutions to the aforementioned initial-boundary value problem. Instead, assuming that appropriate free data has been chosen, we will assume that the global existence of asymptotically flat solutions to the aforementioned initial-boundary value problem to have an analogy to the arguments applied in [
6,
7].
Note that the metric (
27) of superposed Kerr-Schild black holes—though not satisfying Einstein’s equations—is asymptotically flat [
6]. Therefore, it is plausible that solutions to the specific initial-boundary value problem will also be asymptotically flat. The asymptotic form of the metric (
27), along with the geometric assumptions imposed in our construction, guarantee that well-defined ADM mass, center of mass, linear and angular momenta can always be associated with the corresponding multiple black hole system. The main point in this paper is that all the ADM quantities of multiple black hole systems are determined by the rest masses, positions, velocities, and spins of the involved individual black holes.
This result immediately raises the question as to whether any other formulation of the constraints can provide an analogous determination of the ADM quantities. For instance, the method proposed by Bowen and York [
10], in principle, allows for the prescription of the ADM linear and angular momenta by solving the momentum constraint explicitly [
10,
11]. However, to do so, they had to apply a restricted set of basic variables. In particular, to guarantee that the Hamiltonian and momentum constraints decouple, the authors had to assume the vanishing of the mean curvature, in addition to assuming the conformal flatness of the Riemannian metric
. One of the unfavorable consequences of these technical assumptions is that they are known to be so strong that they exclude even the Kerr black hole solution from the outset [
12,
13]. One should also mention here that, within the setup proposed by Bowen and York [
10], there is no way to get an analogous control on the ADM mass or the center of mass.
In the context of the determinacy of the ADM quantities, one should also mention the construction applied in [
4], where, by combining the gluing techniques with Kerr-Schild black holes, an interesting initialization of multiple black hole systems was proposed. Indeed, as Kerr-Schild black holes were applied in [
4] and our proposal also rests upon using these types of black holes, one would expect that analogous determinacy of the ADM quantities applies to both of these approaches. It is, however, not the case, as gluing requires the use of the elliptic method that starts by a conformal rescaling of the basic variables. In turn, gluing gets somewhat implicit, which does not allow—apart from the extreme case with infinitely separated individual black holes—to have complete control on the ADM quantities [
4]. Yet another unfavorable consequence of using the conformal method is that intermediate regions—where the gluing happens—have to be allocated to each of the involved Kerr-Schild black holes. This, however, does not allow to set the initial distances of these black holes to be arbitrary, and, as stated explicitly in [
4], “they must be separated by a distance above a certain threshold”.
It is important to emphasize that our proposal does not impose analogous restrictions on the distances of individual black holes. Yet, it provides an unprecedented complete control on the ADM parameters of multiple black hole systems. As this happens in advance of solving the constraints, an unprecedented fine-tuning of the complete set of ADM parameters of the to-be solutions is possible.
This paper is structured as follows. In
Section 2, first a brief account on the parabolic-hyperbolic form of the constraints is given. This is followed by recalling the notion of asymptotic flatness, the definition of the ADM quantities and the superposition of Kerr-Schild black holes in
Section 2.2 and
Section 2.3. In
Section 3, the choice of the freely specifiable variables and the initial-boundary data, applied in determining multiple black hole initial data, and the pertinent fall-off properties are discussed.
Section 4 is to present our mane result, containing a case by case verification of the statement that the ADM quantities of the superposed Kerr-Schild black holes and the corresponding multiple black hole initial data are pairwise equal to each other. The paper is closed, in
Section 5, by our final remarks.
3. The Initial-Boundary Value Problem
In advance of determining the ADM quantities relevant for asymptotically flat multiple black hole initial data configurations, one has to choose free data for the underlying initial-boundary problem. As a preparation for the asymptotic case first, by a straightforward adaptation of the method applied in [
6], considerations will be restricted to finite cubical domains.
Accordingly, the initial data surface
is chosen to be a cube, centered at the origin in
(see
Figure 1) with boundary comprised by six squares each with edges
.
By choosing the value of
A sufficiently large all the individual black holes will be contained in this cubical domain with suitable margin. The parabolic-hyperbolic system (
9)–(11) has to be solved then as an initial-boundary value problem to which (local) well-posedness is guaranteed (see, e.g., [
18]) in those subregions of
where (
9) is uniformly parabolic.
As in [
6], the ring singularities of the individual black holes are assumed to be located momentarily on the
plane in
, and a foliation of
by
level surfaces will be applied. By an argument, analogous to the one applied in [
6], the principal coefficient
of the parabolic Equation (
9) can be shown to vanish on the
plane dividing
into two disjoint subsets. Accordingly, (
9) is uniformly parabolic on the disjunct subregions of
located above and below of the
surface, and one could also investigate the well-posedness of the parabolic-hyperbolic system (
9)–(11) in these subregions. Nevertheless, as indicated in the introduction, instead of attempting to do so we shall simply assume that the global existence and uniqueness of solutions, along with their proper matching, is guaranteed. Note that the pertinent initial values are supposed to be specified on the horizontal
squares, whereas the boundary values have to be given on the four vertical sides of the cube (see
Figure 1).
The Asymptotic Properties of the Initial-Boundary Data
The model of multiple black hole configurations, as introduced in the previous section, makes use of finite cubical domains. Therefore, to investigate the asymptotic properties of the corresponding initial-boundary data configurations, one has to consider sequences of solutions to the initial-boundary value problem such that the edges tend to infinity. The individual members of such a sequence will differ slightly from global (up to spacelike infinity) solutions to the constraints. Nevertheless, due to the asymptotic flatness of the auxiliary metric (
27), the deviations are expected to be smaller and smaller as the boundary is pushed further and further towards spacelike infinity.
Next we fix the freely specifiable fields
and
to coincide with the auxiliary fields
and
, respectively. The proper fall off property of this part of the data is guaranteed by the fact that (
27) is asymptotically flat. As indicated in the previous sections, the initialization of the constrained variables
and
also happens by utilizing the auxiliary metric (
27).
In proceeding, let us sum up what we already have by hand. Taking into account that the solution is assumed to be asymptotically flat and the Regge-Teitelboim parity conditions also hold, the fields
and
are expected to satisfy
and the parity conditions
Since all the ADM quantities are well-defined for the superposed Kerr-Schild configurations, the auxiliary fields
and
do also satisfy conditions analogous to (
33)–(36).
Consider now the asymptotic expansions of the constrained fields
and
and those of the corresponding auxiliary fields
and
given in terms of various powers of
, with integer
. Referring to what we have inferred concerning sequences of solutions defined on finite cubical domains, it is plausible to assume that as the boundaries of these domains are pushed further and further towards spacelike infinity higher than the leading order terms of the asymptotic expansion play less and less important role. Therefore, as the only sensible asymptotic behavior of the fields
and
associated with the pertinent asymptotically flat solution, the leading order terms in their asymptotic expansions are assumed to be equal to those of
and
, respectively. Obviously, all the higher order contributions may—and, in general, do indeed—differ from each other. Nevertheless, as the fields
,
and
,
are assumed to agree at leading order, respectively, their deviations (These deviations will be denoted by
followed by the pertinent quantities in square brackets, as they are spelled out explicitly in the first two terms of (
37)–(39)) are expected to satisfy the following relations
and
5. Conclusions
Our primary aim was to answer the question raised in the title concerning the freedom we have in specifying the physical parameters of multiple black hole configurations. In doing so, a combination of the parabolic-hyperbolic formulation of constraints and superposed Kerr-Schild black holes were used. We treated only the case of multiple black hole systems where the ring singularities and the speeds of the individual Kerr-Schild black holes were confined to the plane. This, also meant that the spins were required to be aligned or anti-aligned to the z-axis. As there were no further restrictions on the input parameters, a significant number of multiple black hole configurations with immediate physical interest are covered by the investigated set.
The main result of this paper can be formulated as follows:
Theorem 1. Suppose that an asymptotically flat solution to the initial-boundary value problem—deduced from the parabolic-hyperbolic form of the constrains, (9)–(11)—exists such that the free data is chosen, as described in Section 3, by applying the superposed Kerr-Schild metric (27). Then, the ADM mass, center of mass, linear and angular momenta, relevant for the initialization of the corresponding multiple black hole system, can be given, as in Section 4.1, Section 4.2, Section 4.3, Section 4.4, in terms of the rest masses, positions, velocities, and spins of the involved individual Kerr-Schild black holes. In addition, all of these parameters can be prescribed in advance of solving the constraints. Several remarkable features characterize the applied initial data construction. First, as we do not use conformal rescalings, our method retains the physically distinguished nature of and . Second, the input parameters are the rest masses, the sizes, and orientations of the displacements, velocities, and spins of the individual black holes. Note that these are essentially the same as the input parameters of the post-Newtonian (PN) formalism. This provides significant interrelations between the PN and our fully relativistic setups. In particular, physically adequate choices of the orbital parameters could be made using the insights earned within the PN. Notably, as shown in this paper, each of the global ADM charges can also be given in terms of the input parameters, which is unprecedented in other methods to solve the constraints. More strikingly, the ADM mass, center of mass, linear and angular momenta of the binary system can be fixed in advance of solving the constraints.
Despite the advantages discussed above, it is important to keep in mind that there is room for further investigations. There is an obvious interest to generalize the applied initial data construction to the case where the speed and spin vectors of involved individual Kerr-Schild black holes are arbitrarily pointing. Once this is done, there will also be a need to generalize the results covered by this paper.
Another physically important generalization could be to replace asymptotic flatness with more realistic geometric assumptions. For instance, black holes could be placed in the environment of the expanding universe modeled by the standard Friedman-Lemaitre-Robertson-Walker solutions. Similarly, investigations of black hole systems in the McVittie background, studied, e.g., in [
19,
20], could also be of interest. Note that such a replacement will require a more substantial generalization of the construction as even the freely specifiable variables have to be altered significantly.
Both of the indicated problems—which certainly deserve further attention—are left open for future investigations.