Next Article in Journal
A General Class of Differential Hemivariational Inequalities Systems in Reflexive Banach Spaces
Next Article in Special Issue
Dynamic Analytical Solution of a Charged Dilaton Black Hole
Previous Article in Journal
A Novel Hybrid Approach for Optimal Placement of Non-Dispatchable Distributed Generations in Radial Distribution System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Can We Prescribe the Physical Parameters of Multiple Black Holes?

Wigner RCP, Konkoly Thege Miklós út 29-33, H-1121 Budapest, Hungary
Mathematics 2021, 9(24), 3170; https://doi.org/10.3390/math9243170
Submission received: 6 November 2021 / Revised: 29 November 2021 / Accepted: 6 December 2021 / Published: 9 December 2021
(This article belongs to the Special Issue Mathematics of Black Holes)

Abstract

:
The parabolic-hyperbolic form of the constraints and superposed Kerr-Schild black holes have already been used to provide a radically new initialization of binary black hole configurations. The method generalizes straightforwardly to multiple black hole systems. This paper is to verify that each of the global Arnowitt-Deser-Misner quantities of the constructed multiple black hole initial data can always be prescribed, as desired, in advance of solving the constraints. These global charges are shown to be uniquely determined by the physical parameters of the involved individual Kerr-Schild black holes.

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) C 2 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 t = c o n s t 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 h i j . 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.

2. Preliminaries

Initial data relevant for the vacuum Einstein’s equations is comprised of a Riemannian metric h i j and a symmetric tensor field K i j . Both of these fields are assumed to be given on a three-dimensional manifold Σ . They are not arbitrary as they have to satisfy the constraints which read as (see, e.g., [3])
( 3 ) R + K j j 2 K i j K i j = 0
D j K j i D i K j j = 0 ,
where ( 3 ) R and D i denote the scalar curvature and the covariant derivative operator associated with h i j , respectively.

2.1. The Parabolic-Hyperbolic Form of Constraints

The essential steps in deriving the parabolic-hyperbolic form of the constraints are as follows: Assume, for simplicity, that there exists a smooth function ρ : Σ R such that the ρ = c o n s t surfaces (denoted also by S ρ ) provide a foliation of Σ . We assume that the S ρ surfaces are homologous to each other and the transversal one form D i ρ to these level surfaces does not vanish on Σ . The unite normal n ^ i to the S ρ surfaces is given then as
n ^ i = h e f ( D e ρ ) ( D f ρ ) 1 2 D i ρ .
Choosing then a vector field ρ i on Σ such that ρ i i ρ = 1 , and considering its parallel and orthogonal parts we get
ρ i = N ^ n ^ i + N ^ i ,
where N ^ and N ^ i stand for the lapse and shift of ρ i , respectively, and n ^ i = h i j n ^ j .
Analogous decomposition of the metric h i j and the symmetric tensor field K i j gives
h i j = γ ^ i j + n ^ i n ^ j , and K i j = κ n ^ i n ^ j + n ^ i k j + n ^ j k i + K i j ,
where γ ^ i j = γ ^ e i γ ^ f j h e f is the induced metric on the ρ = c o n s t level surfaces, and
κ = n ^ k n ^ l K k l , k i = γ ^ k i n ^ l K k l , K i j = γ ^ k i γ ^ l j K k l ,
where γ ^ k i denotes the projection operator
γ ^ k i = h k i n ^ k n ^ i .
It is also rewarding to introduce the trace and trace-free part of K i j defined as
K l l = γ ^ k l K k l and K i j = K i j 1 2 γ ^ i j K l l .
Using the new variables N ^ , N ^ i , γ ^ i j , κ , k i , K i j , K l l the constraints can be seen to be equivalent to the parabolic-hyperbolic system for N ^ , k i and K l l [5]
K [ ( ρ N ^ ) N ^ l ( D ^ l N ^ ) ] = N ^ 2 ( D ^ l D ^ l N ^ ) + A N ^ + B N ^ 3 ,
L n ^ k i 1 2 D ^ i ( K l l ) D ^ i κ + D ^ l K l i + N ^ K k i + [ κ 1 2 ( K l l ) ] n ^ ˙ i n ^ · l K l i = 0
L n ^ ( K l l ) D ^ l k l N ^ K [ κ 1 2 ( K l l ) ] + N ^ K k l K k l + 2 n ^ ˙ l k l = 0 .
In (9)–(11), D ^ i denotes the covariant derivative operator associated with γ ^ i j and
K i j = 1 2 L ρ γ ^ i j D ^ ( i N ^ j ) ,
K = 1 2 γ ^ i j L ρ γ ^ i j D ^ j N ^ j ,
A = ( ρ K ) N ^ l ( D ^ l K ) + 1 2 [ K 2 + K k l K k l ]
B = 1 2 R ^ + 2 κ ( K l l ) + 1 2 ( K l l ) 2 2 k l k l K k l K k l
n ^ k = n ^ l D l n ^ k = D ^ k ( ln N ^ ) .
The variables N ^ i , γ ^ i j , κ and K i j are unconstrained whence they are freely specifiable throughout Σ . The well-posedness of the coupled system (9)–(11) is guaranteed if (9) is uniformly parabolic. It was shown in [5] that this happens in those subregions of Σ where K is either strictly positive or negative. Note also that, as K depends on the freely specifiable fields γ ^ i j and N ^ i exclusively, the sign of K (at least locally) is adjustable according to the needs of specific problems to be solved [5,6].
Note that, in addition to the freely specifiable variables, on one of the ρ = c o n s t level surfaces initial data has also to be chosen for the constrained variables [5]. Once this has been done, i.e., smooth data has been chosen for N ^ , k i and K l l , then a unique smooth solution exists to (9)–(11) in the domain of dependence of that ρ = c o n s t level surface. This, in general, can be seen to be global as the hyperbolic part of the system, see (10) and (11), is linear in k i and K l l . It is also important that the fields h i j and K i j that can be reconstructed from such a solution, and from the freely specifiable fields, do satisfy the Hamiltonian and momentum constraints (1) and (2).

2.2. Asymptotic Flattness and the ADM Quantities

Before turning to the specific class of solutions to the parabolic-hyperbolic system (9)–(11), it is rewarding to have a glance at the generic notion of asymptotic flatness, along with the conditions ensuring the existence of well-defined ADM charges.
Our model is based on the use of the superposed Kerr-Schild metric, whence, the singularities of the multiple black hole system will be arranged to be located in a finite ball B in R 3 . Therefore it suffices to assume the existence of a single asymptotically flat end. (See Section 3 below for further specifications.)
The initial data set ( Σ , h i j , K i j ) is called strictly (Note that in Section 4—likewise in many other analogous investigations (see, e.g., [14])—weaker fall-off conditions could also be used. Nevertheless, hereafter, for convenience, and for definiteness, the above recalled stronger fall-off conditions will be applied.) asymptotically Euclidean [4], to order , if in exterior to B i.e., in Σ \ B , admissible asymptotically flat coordinates x i = ( x , y , z ) exist such that
| α h i j δ i j ( x ) | = O ( | x | | α | 1 ) , | β K i j ( x ) | = O ( | x | | β | 2 ) ,
hold, where | x | = x 1 2 + x 2 2 + x 3 2 , and where α , with multi index α = α 1 + α 2 + α 3 , stands for the composition of partial derivative operators x 1 α 1 x 2 α 2 x 3 α 3 , and for the multi indices α and β , for some value of , the inequalities | α | + 1 and | β | hold. For the arguments applied in this paper 1 will suffice. In most cases the operators x i will also be abbreviated as i .
It is known that conditions in (17) can only guarantee the existence and finiteness of the four-momentum, and to have, in addition, a well-defined center of mass and angular momentum the so-called Regge-Teitelboim asymptotic parity conditions need to be used which, in admissible coordinates, read as
| α h i j ( x ) h i j ( x ) | = O ( | x | | α | 2 ) , | β K i j ( x ) K i j ( x ) | = O ( | x | | β | 3 ) .
Assuming that both the asymptotic flatness and the Regge-Teitelboim conditions hold, the ADM mass, center of mass, linear and angular momenta are given by the flux integrals [4]
M A D M = 1 16 π i h i j j h i i n j d S
M A D M d i = 1 16 π x i k h k j j h k k h k j δ k i h k k δ i j n j d S
P i A D M = 1 8 π K i j h k j K l l n j d S
J i A D M = 1 8 π K k j h k j K l l ϵ i l k x l n j d S ,
where the symbol is meant to denote limits of integrals over spheres while their radii tend to infinity, whereas n i and d S denote the outward pointing unit normal and the volume element of the individual spheres of the sequences, respectively. Note that in (22) ϵ i j k x j stands for the components of the three rotational Killing vector fields, defined with respect to the applied admissible asymptotically Euclidean coordinates x i .

2.3. Superposed Kerr-Schild Black Holes

The Kerr solution [15] is known to take the Kerr-Schild form given as
g α β = η α β + 2 H α β ,
where
H = r 3 M r 4 + a 2 z 2 and α = 1 , r x + a y r 2 + a 2 , r y a x r 2 + a 2 , z r ,
and where the Boyer-Lindquist radial coordinate r is related to the spatial part of the inertial coordinates x α = ( t , x , y , z ) —which are asymptotically flat admissible coordinates—via the implicit relation
x 2 + y 2 r 2 + a 2 + z 2 r 2 = 1 .
It is well-known that generic displaced, boosted and spinning black holes can be produced by performing suitable Poincaré transformations on a Kerr black hole. It is also important that the Kerr-Schild metric is form-invariant under these transformations. In particular, even if a Lorentz transformation x α = Λ α β x β is performed the metric g α β in the new coordinates will retain the Kerr-Schild form g α β = η α β + 2 H α β , where
H ( x α ) = H [ Λ α β ] 1 x β and β ( x ε ) = Λ α β α [ Λ ε φ ] 1 x φ .
As boosts and spatial rotations are special Lorentz transformations, it is straightforward to construct models of moving black holes with preferably oriented speed and spin by performing suitable combinations of boosts and rotations on a Kerr black hole that is in rest and suitably oriented with respect to an auxiliary Minkowski background. Note also that displacement of these boosted and spinning black holes may be represented by a straightforward change in the argument of H and α whence all the displaced, boosted and spinning individual black holes may be produced by applying the indicated transformations.
Now we are almost ready to combine the parabolic-hyperbolic form of the constraints with superposed Kerr-Schild black holes. As indicated in the introduction solving the constraints in their parabolic-hyperbolic form (9)–(11) requires specifications of the unconstrained variables N ^ i , γ ^ i j , κ and K i j everywhere on Σ , and, in addition, an initialization of constrained variables N ^ , k i and K l l on one of the ρ = c o n s t level surfaces.
In virtue of the results in [6] it is rewarding to start with the auxiliary metric
g α β = η α β + I = 1 N 2 H [ I ] α [ I ] β [ I ]
yielded by superposing the contributions of individual black holes represented by the ( H [ I ] , α [ I ] ) , I = 1 , , N , pairs. This metric retains much of the algebraic simplicity of (23). For instance, the vector fields α [ I ] remain null with respect to the background Minkowski metric, and they satisfy the geodesic equation β [ I ] β α [ I ] = 0 , the functions H [ I ] satisfy the background wave equation η α β α β H [ I ] = 0 (see, e.g., [16]).
From now on, to distinguish the real physical quantities from the auxiliary ones deduced from superimposed Kerr-Schild form (27), the latter will be labeled by the “pre upper index” ( A ) . For instance, ( A ) h i j will stand for the three-metric
( A ) h i j = δ i j + I = 1 N 2 H [ I ] i [ I ] j [ I ]
induced by the superposed Kerr-Schild form (27) on t = c o n s t hypersurfaces, where t is the time coordinate of the inertial system x α = ( t , x , y , z ) . Clearly, such a t = c o n s t hypersurface may be assumed to be a Kerr-Schild time slice for each of the individual black holes, which implies that topologically it is simply the complement of the individual “ring” singularities in R 3 as it was indicated in Section 2.2.
Note that the metric (27) is not a solution to Einstein’s equations yet it is asymptotically flat. This, in particular, means that the integrability of the scalar curvature of ( A ) h i j (necessary to have, for instance, well-defined ADM mass [17]) is guaranteed, and that conditions in (17) and (18) hold for the metric in (27), and, in turn, the superposed Kerr-Schild black holes can always be assigned with well-defined ADM charges. These quantities can be determined either by evaluating the flux integrals given in (19)–(22) or by taking into account the Poincaré transformations associated with displacements, boosts and rotations performed on the individual black holes. In either way we get the remarkably simple relations
M A D M = I = 1 N γ [ I ] M [ I ] ,
M A D M d = I = 1 N γ [ I ] M [ I ] d [ I ]
P A D M = I = 1 N γ [ I ] M [ I ] v [ I ] ,
J A D M = I = 1 N γ [ I ] M [ I ] d [ I ] × v [ I ] + M [ I ] a [ I ] s [ I ] ,
where v [ I ] and d [ I ] are the velocity and position vectors of the individual black holes with respect to the background inertial frame, and s [ I ] denotes the unit vectors pointing to directions of the individual spin vectors.

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 R 3 (see Figure 1) with boundary comprised by six squares each with edges 2 A .
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 z = 0 plane in R 3 , and a foliation of Σ by z = c o n s t level surfaces will be applied. By an argument, analogous to the one applied in [6], the principal coefficient K of the parabolic Equation (9) can be shown to vanish on the z = 0 plane dividing Σ into two disjoint subsets. Accordingly, (9) is uniformly parabolic on the disjunct subregions of Σ located above and below of the z = 0 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 z = ± A 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 N ^ i , γ ^ i j , κ and K i j to coincide with the auxiliary fields ( A ) N ^ i , ( A ) γ ^ i j , ( A ) κ and ( A ) K i j , 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 N ^ , K l l and k i 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 N ^ , K l l and k i are expected to satisfy
| α N ^ 1 ( x ) | = O ( | x | | α | 1 ) , | β K l l ( x ) | = O ( | x | | β | 2 ) , | β k i ( x ) | = O ( | x | | β | 2 ) ,
and the parity conditions
| α N ^ ( x ) N ^ ( x ) | = O ( | x | | α | 2 )
| β K l l ( x ) K l l ( x ) | = O ( | x | | β | 3 )
| β k i ( x ) k i ( x ) | = O ( | x | | β | 3 ) .
Since all the ADM quantities are well-defined for the superposed Kerr-Schild configurations, the auxiliary fields ( A ) N ^ , ( A ) K l l and ( A ) k i do also satisfy conditions analogous to (33)–(36).
Consider now the asymptotic expansions of the constrained fields N ^ , K l l and k i and those of the corresponding auxiliary fields ( A ) N ^ , ( A ) K l l and ( A ) k i given in terms of various powers of 1 / | x | n , with integer n > 0 . 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 N ^ , K l l and k i associated with the pertinent asymptotically flat solution, the leading order terms in their asymptotic expansions are assumed to be equal to those of ( A ) N ^ , ( A ) K l l and ( A ) k i , respectively. Obviously, all the higher order contributions may—and, in general, do indeed—differ from each other. Nevertheless, as the fields N ^ , K l l , k i and ( A ) N ^ , ( A ) K l l , ( A ) k i 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
| α Δ [ N ^ ] ( x ) | = | α N ^ ( A ) N ^ ( x ) | = O ( | x | | α | 2 )
| β Δ [ K l l ] ( x ) | = | β K l l ( A ) K l l ( x ) | = O ( | x | | β | 3 )
| β Δ [ k i ] ( x ) | = | β k i ( A ) k i ( x ) | = O ( | x | | β | 3 )
and
| α Δ [ N ^ ] ( x ) Δ [ N ^ ] ( x ) | = O ( | x | | α | 3 )
| β Δ [ K l l ] ( x ) Δ [ K l l ] ( x ) | = O ( | x | | β | 4 )
| β Δ [ k i ] ( x ) Δ [ k i ] ( x ) | = O ( | x | | β | 4 ) .

4. The Determination of the ADM Quantities

Now we are in the position to compare the ADM quantities of the superposed Kerr-Schild metric with those of the corresponding multiple black hole initial data. In particular, we shall show that conditions (37)–(42), along with the choices we made for the other auxiliary variables, guarantee that the two sets of ADM quantities are pairwise equal to each other. In the following subsections a case by case verification of this claim will be provided.

4.1. The ADM Mass

Consider first the ADM mass. Start by replacing the flux integral applied in (19) by a slightly different flux integral
M A D M = 1 16 π i h i j j h i i n j d C ,
where the symbol denotes the limit of integrals over the boundary of a sequence of co-centered cubes while the length of their edges tend to infinity, whereas n i and d C denote the outward pointing unit normal vector and the volume element on the squares bounding individual cubes in this sequence.
At the first glance it may not be obvious that the flux integrals over concentric spheres can be replaced by flux integrals over boundaries of co-centered cubical regions. Note, however, that to any individual member of these cubes there always exists a minimal radius sphere that contains the cube, and a maximal radius sphere that is contained by the cube. Clearly, either the minimal or maximal radius of spheres are applied to construct a sequence, the flux integrals defined with respect to them tend to the ADM mass. Thereby, the flux integrals evaluated on the boundaries of the cubes have to tend to the ADM mass as well. Accordingly, the limits of the integrals in (19) and (44) have to be equal to each other.
Returning to the main line of the argument, note that our aim here is to show that the difference Δ [ M A D M ] = M A D M ( A ) M A D M between the ADM mass of the physical solution and that of the superposed Kerr-Schild black holes is zero. To see that this is indeed the case, note first that by virtue of (43)
Δ [ M A D M ] = 1 16 π i ( h i j ( A ) h i j ) j ( h i i ( A ) h i i ) n j d C .
To evaluate the integrands we need to determine first the involved derivatives. In doing so note that the difference h i j ( A ) h i j , in virtue of the relations
γ ^ i j = ( A ) γ ^ i j , n ^ i = N ^ δ i z and ( A ) n ^ i = ( A ) N ^ δ i z ,
reads as
h i j ( A ) h i j = n ^ i n ^ j ( A ) n ^ i ( A ) n ^ j = δ i z δ j z ( N ^ 2 ( A ) N ^ 2 ) = δ i z δ j z Δ [ N ^ 2 ] .
(46) implies then that
k ( h i j ( A ) h i j ) = δ i z δ j z k ( Δ [ N ^ 2 ] ) ,
and also that
i ( h i j ( A ) h i j ) = δ j z z ( Δ [ N ^ 2 ] ) and j ( h i i ( A ) h i i ) = j ( Δ [ N ^ 2 ] ) .
By combining (43), (44) and (48) we get then that
Δ [ M A D M ] = 1 16 π δ j z z ( Δ [ N ^ 2 ] ) j ( Δ [ N ^ 2 ] ) n j d C = 1 16 π x ( Δ [ N ^ 2 ] ) n ± x + y ( Δ [ N ^ 2 ] ) n ± y d C ,
where n ± x and n ± y denote the outward pointing unit normal vectors to the squares in the x = ± A and y = ± A plains, bounding the cubical region on Figure 1. Since n x = n + x and n y = n + y , the two terms given explicitly in the integrands stand indeed for four terms.
By applying then the replacements y x υ and z x ζ in the first term of (49) evaluated on the squares in the x = ± A plains, and also the replacements x y ξ and z y ζ in the second term of (49) evaluated on the squares in the y = ± A plains we get
Δ [ M A D M ] = 1 16 π 1 , 1 1 , 1 lim x x 2 x ( Δ [ N ^ 2 ] ) n ± x d υ d ζ + 1 , 1 1 , 1 lim y y 2 y ( Δ [ N ^ 2 ] ) n ± y d ξ d ζ ,
where n ± x and n ± y are scalars taking the values n ± x = ± 1 and n ± y = ± 1 on the squares in the x = ± A and y = ± A plains, respectively. Note that the replacements applied in the above integral transformations are analogous to the ones used in case of sequences of spheres, where the integrals formally are given over a unit sphere with angular coordinates θ and ϕ ranging through their usual intervals. In both cases with the help of these integral transformations the limits of integrals can be evaluated by inspecting the limits of the yielded integrands.
Using then
Δ [ N ^ 2 ] = ( Δ [ N ^ ] ) 2 + 2 ( A ) N ^ Δ [ N ^ ] ,
along with (37), applied for Δ [ N ^ ] , and (34), applied for ( A ) N ^ , respectively, we get that each of the terms x i Δ [ N ^ ] , x i 2 ( x i Δ [ N ^ ] ) and x i ( x i ( A ) N ^ ) appearing in the integrand of (50) is at most of order O ( | x | 1 ) . Here the relations x i = ( x i / | x | ) | x | = x i | x | and | x i | 1 were also used. This, in turn, implies that the limits exist and they all vanish which verifies that Δ [ M A D M ] = 0 .

4.2. The Center of Mass

Rephrasing (Section 4.2), by using integrals over cubical domains, we get that the center of mass can be given by the flux integral
M A D M d l = 1 16 π x l k h k j j h k k h k j δ k l h k k δ l j n j d C .
Taking now into account the consequences of (46) we get
Δ [ h k j ] = δ k z δ j z Δ [ N ^ 2 ] and Δ [ h k k ] = Δ [ N ^ 2 ] .
A straightforward calculation, consisting of steps analogous to the ones applied in the previous subsection, yields then
Δ [ M A D M d l ] = 1 16 π { x l k ( h k j ( A ) h k j ) j ( h k k ( A ) h k k ) . ( h k j ( A ) h k j ) δ k l ( h k k ( A ) h k k ) δ l j } n j d C = 1 16 π x l k ( Δ [ h k j ] ) j ( Δ [ h k k ] ) Δ [ h k j ] δ k l Δ [ h k k ] δ l j n j d C = 1 16 π x l x ( Δ [ N ^ 2 ] ) n ± x + y ( Δ [ N ^ 2 ] ) n ± y + δ l x n ± x + δ l y n ± y Δ [ N ^ 2 ] d C = 1 16 π x l x ( Δ [ N ^ 2 ] ) δ l x Δ [ N ^ 2 ] n ± x + x l y ( Δ [ N ^ 2 ] ) δ l y Δ [ N ^ 2 ] n ± y d C .
Accordingly, for the x-component of the deviation Δ [ M A D M d l ]
Δ [ M A D M d x ] = 1 16 π x x ( Δ [ N ^ 2 ] ) Δ [ N ^ 2 ] n ± x + x y ( Δ [ N ^ 2 ] ) ] n ± y d C = 1 16 π 1 , 1 1 , 1 lim x x 3 x ( Δ [ N ^ 2 ] ) x 2 Δ [ N ^ 2 ] n ± x d υ d ζ + 1 , 1 1 , 1 lim y y 3 y ( Δ [ N ^ 2 ] ) n ± y ξ d ξ d ζ
hold.
Taking into account (51), along with the relations n + x = n x = 1 and n + y = n y = 1 , in virtue of (34), (37) and the Regge-Teitelboim condition (40), it follows that either of the terms x 2 Δ [ N ^ 2 ] , x 3 ( x Δ [ N ^ 2 ] ) , y 3 ( y Δ [ N ^ 2 ] ) and x i 2 ( x i ( A ) N ^ ) in the integrands, is at most of order O ( | x | 1 ) . This, in turn, implies, as above, that the pertinent limits exist and they vanish as we intended to show.
In virtue of the relations,
Δ [ M A D M d y ] = 1 16 π 1 , 1 1 , 1 lim x x 3 x ( Δ [ N ^ 2 ] ) n ± x υ d υ d ζ + 1 , 1 1 , 1 lim y y 3 y ( Δ [ N ^ 2 ] ) y 2 Δ [ N ^ 2 ] n ± y d ξ d ζ
and
Δ [ M A D M d z ] = 1 16 π 1 , 1 1 , 1 lim x x 3 x ( Δ [ N ^ 2 ] ) n ± x ζ d υ d ζ + 1 , 1 1 , 1 lim y y 3 y ( Δ [ N ^ 2 ] ) n ± y ζ d ξ d ζ ,
arguments, analogous to the one applied above, can be used to show the vanishing of the y- and z-components of the deviation Δ [ M A D M d l ] .

4.3. The Linear Momentum

Consider now the linear momentum determined by the flux integral
P i A D M = 1 8 π K i j h k j K l l n j d C .
In verifying that the individual components of Δ [ P i A D M ] vanish, respectively, it is rewarding to rephrase first the term K i j h i j K l l in (58) in terms of the new variables we introduced. In doing so we get, in virtue of (5)–(8), that
K i j = κ n ^ i n ^ j + n ^ i k j + n ^ j k i + K i j + 1 2 γ ^ i j K l l ,
and
K l l = κ + K l l .
It is straightforward to see then that
K i j h i j K l l = n ^ i k j + n ^ j k i + K i j + 1 2 γ ^ i j K l l γ ^ i j κ + ( γ ^ i j + n ^ i n ^ j ) K l l ,
and, by applying (45), that
Δ [ K i j h i j K l l ] = δ i z Δ [ N ^ k j ] + δ j z Δ [ N ^ k i ] 1 2 γ ^ i j Δ [ K l l ] δ i z δ j z Δ [ N ^ 2 K l l ] .
It follows then that for the individual components of Δ [ P i A D M ] the relations
Δ [ P x A D M ] = 1 8 π 1 2 γ ^ x x Δ [ K l l ] n ± x 1 2 γ ^ x y Δ [ K l l ] n ± y + Δ [ N ^ k x ] n ± z d C
Δ [ P y A D M ] = 1 8 π 1 2 γ ^ y x Δ [ K l l ] n ± x 1 2 γ ^ y y Δ [ K l l ] n ± y + Δ [ N ^ k y ] n ± z d C
Δ [ P z A D M ] = 1 8 π Δ [ N ^ k x ] n ± x + Δ [ N ^ k y ] n ± y Δ [ N ^ 2 K l l ] n ± z d C
hold, where k i n ^ i = 0 and γ ^ i j n ^ i = 0 had also been used.
In verifying that each of the components of Δ [ P A D M ] vanish—besides the replacements y x υ z x ζ , and x y ξ z y ζ on the x = ± A and y = ± A surfaces, respectively—the transformations x z ξ and y z υ have also to be performed on the z = ± A surfaces since the integrands there, as opposed to the previous two cases, are not identically zero. By applying the corresponding integral transformations we get from (63) that
Δ [ P x A D M ] = 1 16 π 1 , 1 1 , 1 lim x γ ^ x x n ± x x 2 Δ [ K l l ] d υ d ζ + 1 , 1 1 , 1 lim y γ ^ x y n ± y y 2 Δ [ K l l ] d ξ d ζ 2 1 , 1 1 , 1 lim z n ± z z 2 Δ [ N ^ ] ( A ) k x + ( A ) N ^ z 2 Δ [ k x ] d ξ d υ .
The fall off conditions in (37) and (39), along with the one satisfied, in virtue of (33), by ( A ) k x , imply that all the terms in braces are at most of order O ( | x | 1 ) which, in turn, verifies that Δ [ P x A D M ] vanishes.
As the terms in (64) are similar to those in (63), a completely analogous argument applies to the y-component of Δ [ P i A D M ] .
To show that z-component, Δ [ P z A D M ] , does also vanish note first that (65) can be given as
Δ [ P z A D M ] = 1 8 π 1 , 1 1 , 1 lim x n ± x x 2 Δ [ N ^ ] ( A ) k x + ( A ) N ^ x 2 Δ [ k x ] d υ d ζ + 1 , 1 1 , 1 lim y n ± y y 2 Δ [ N ^ ] ( A ) k y + ( A ) N ^ y 2 Δ [ k y ] d ξ d ζ 1 , 1 1 , 1 lim z n ± z z 2 Δ [ N ^ 2 ] ( A ) K l l + ( A ) N ^ 2 z 2 Δ [ K l l ] d ξ d υ .
The fall off conditions listed in relations (37)–(39), along with those satisfied, in virtue of (33), by ( A ) k x , ( A ) k y and ( A ) K l l , imply again that all the terms in braces are at most of order O ( | x | 1 ) which verifies that Δ [ P z A D M ] vanishes as we desired to show.
Note that, as expected, in verifying that the physical ADM mass and linear momentum are equal to the ADM mass and linear momentum of the superposed Kerr-Schild black holes, no use of the Regge-Teitelboim conditions had to be made.

4.4. The Angular Momentum

Consider, finally, the angular momentum determined by the flux integral
J i A D M = 1 8 π K k l h k l K m m ϵ i j k x j n l d C .
Thus we have
Δ [ J i A D M ] = J i A D M ( A ) J i A D M = 1 8 π Δ [ J ˜ i l ] n l d C ,
where Δ [ J ˜ i l ] = Δ K k l h k l K m m ϵ i j k x j .
By a direct calculation, consisting of steps analogous to the ones applied in the previous subsection in evaluating Δ [ K i j h i j K l l ] in (62), the relations
Δ [ J ˜ x l ] = y Δ [ N ^ k l ] δ l z Δ [ N ^ 2 K m m ] z δ l z Δ [ N ^ k y ] 1 2 γ ^ l y Δ [ K m m ]
Δ [ J ˜ y l ] = z δ l z Δ [ N ^ k x ] 1 2 γ ^ l x Δ [ K m m ] x Δ [ N ^ k l ] δ l z Δ [ N ^ 2 K m m ]
Δ [ J ˜ z l ] = x δ l z Δ [ N ^ k y ] 1 2 γ ^ l y Δ [ K m m ] y δ l z Δ [ N ^ k x ] 1 2 γ ^ l x Δ [ K m m ]
can be seen to hold.
It is then straightforward to verify that
Δ [ J x A D M ] = 1 8 π y Δ [ N ^ k l ] + 1 2 z γ ^ x y Δ [ K m m ] n ± x + y Δ [ N ^ k y ] + 1 2 z γ ^ y y Δ [ K m m ] n ± y y Δ [ N ^ 2 K m m ] + z Δ [ N ^ k y ] n ± z d C
from which—by applying the integral transformations already used several times in the previous subsections, on the individual x i = ± A plains—we get
Δ [ J x A D M ] = 1 8 π 1 , 1 1 , 1 lim x x 3 υ Δ [ N ^ k l ] + 1 2 ζ γ ^ x y Δ [ K m m ] n ± x d υ d ζ + 1 , 1 1 , 1 lim y y 3 Δ [ N ^ k y ] + 1 2 ζ γ ^ y y Δ [ K m m ] n ± y d ξ d ζ 1 , 1 1 , 1 lim z z 3 υ Δ [ N ^ 2 K m m ] + Δ [ N ^ k y ] n ± z d ξ d υ .
Taking now into account (51), along with the relations n + x = n x = 1 and n + y = n y = 1 , and also the boundedness of the components of γ ^ i j and that of the coordinates ξ , υ , ζ , in virtue of (34)–(39) and the Regge-Teitelboim condition (40)–(42), it follows that either of the terms involved by the integrands, is at most of order O ( | x | 1 ) which implies the vanishing of Δ [ J x A D M ] .
In virtue of the use of analogous terms in Δ [ J ˜ x l ] and Δ [ J ˜ y l ] a completely analogous argument can be seen to apply to the y-component of Δ [ J i A D M ] .
Finally, Δ [ J z A D M ] can be evaluated as
Δ [ J z A D M ] = 1 16 π 1 , 1 1 , 1 lim x x 3 γ ^ x y υ γ ^ x x Δ [ K m m ] n ± x d υ d ζ + 1 , 1 1 , 1 lim y y 3 ξ γ ^ y y γ ^ y x Δ [ K m m ] n ± y d ξ d ζ 2 1 , 1 1 , 1 lim z z 3 ξ Δ [ N ^ k y ] υ Δ [ N ^ k x ] n ± z d ξ d υ .
Then, by making use of the fall off and Regge-Teitelboim conditions, (37)–(42), and by an argument that has already been applied several times above, the vanishing of Δ [ J z A D M ] can be inferred.
Putting all the results of the previous subsections together we get the desired verification of Theorem 1.

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 z = 0 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 h i j and K i j . 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.

Funding

This research was funded by NKFIH grant K-115434.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author is grateful to Piotr Chruściel, Andor Frenkel, Jan Metzger, Christopher Nerz, Harald Pfeiffer, Bob Wald and Jeff Winicour for helpful comments. Thanks are due to the Albert Einstein Institute in Golm for its kind hospitality. This work was also supported in parts by the NKFIH grant K-115434.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Lichnerowicz, A. L’integration des Equations de la Gravitation Relativiste et le Probleme des n Corps. J. Math. Pures Appl. 1944, 23, 39–63. [Google Scholar]
  2. York, J.W. Role of conformal three-geometry in the dynamics of gravitation. Phys. Rev. Lett. 1972, 28, 1082–1085. [Google Scholar] [CrossRef]
  3. Choquet-Bruhat, Y. General Relativity and Einstein’s Equations; Oxford University Press Inc.: New York, NY, USA, 2009. [Google Scholar]
  4. Chruściel, P.T.; Corvino, J.; Isenberg, J. Constructions of N-body initial data sets in general relativity. Commun. Math. Phys. 2011, 304, 637–647. [Google Scholar] [CrossRef] [Green Version]
  5. Rácz, I. Constrains as evolutionary systems. Class. Quant. Grav. 2016, 33, 015014. [Google Scholar] [CrossRef] [Green Version]
  6. Rácz, I. A simple method of constructing binary black hole initial data. Astron. Rep. 2018, 62, 953–958. [Google Scholar] [CrossRef] [Green Version]
  7. Rácz, I. Supplemental Material. 2016. Available online: https://www.kfki.hu/$\sim$iracz/SM-BH-data.pdf (accessed on 10 November 2018).
  8. Doulis, Y. Construction of high precision numerical single and binary black hole initial data. Phys. Rev. D 2019, 100, 024064. [Google Scholar] [CrossRef] [Green Version]
  9. Nakonieczna, A.; Nakonieczny, L.; Rácz, I. Black hole initial data by numerical integration of the parabolic–hyperbolic form of the constraints. Int. J. Mod. Phys. D 2021, 2150111. [Google Scholar] [CrossRef]
  10. Bowen, J.M.; York, J.W. Time-asymmetric initial data for black holes and black-hole collisions. Phys. Rev. D 1980, 21, 2047–2056. [Google Scholar] [CrossRef]
  11. Baumgarte, T.W.; Shapiro, S.L. Numerical Relativity: Solving Einstein’s Equations on the Computer; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  12. Garat, A.; Price, R.H. Nonexistence of conformally flat slices of the Kerr spacetime. Phys. Rev. D 2000, 61, 124011. [Google Scholar] [CrossRef] [Green Version]
  13. Kroon, J.A.V. Asymptotic expansions of the Cotton? York tensor on slices of stationary spacetimes. Class. Quant. Grav. 2004, 21, 3237–3250. [Google Scholar] [CrossRef] [Green Version]
  14. Chruściel, P.T. A remark on the positive energy theorem. Class. Quant. Grav. 1986, 3, L115–L121. [Google Scholar] [CrossRef]
  15. Kerr, R.P.; Schild, A. Some algebraically degenerate solutions of Einstein’s gravitational field equations. Proc. Symp. Appl. Math. 1965, 17, 199. [Google Scholar]
  16. Rácz, I.; Winicour, J. Black hole initial data without elliptic equation. Phys. Rev. D 2015, 91, 124013. [Google Scholar] [CrossRef] [Green Version]
  17. Bartnik, R. The mass of an asymptotically flat manifold. Commun. Pure Appl. Math. 1986, 39, 661–693. [Google Scholar] [CrossRef] [Green Version]
  18. Kreiss, H.-O.; Lorenz, J. Initial-Boundary Value Problems and the Navier-Stokes Equations; Academic Press: Boston, MA, USA, 1989; Reprinted in SIAM Classic 2004; Available online: https://my.siam.org/Store/Product/viewproduct/?ProductId=786 (accessed on 10 November 2021).
  19. Antoniou, I.; Papadopoulos, D.; Perivolaropoulos, L. Spinning particle orbits around a black hole in an expanding background. Class. Quant. Grav. 2019, 36, 085002. [Google Scholar] [CrossRef] [Green Version]
  20. Nolan, B.C. Local properties and global structure of McVittie spacetimes with non-flat Friedmann–Lemaître–Robertson–Walker backgrounds. Class. Quant. Grav. 2017, 34, 225002. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The initial data surface Σ , with a triple black hole system, is chosen to be the cube centered at the origin in R 3 with edges 2 A . The initial data, to the system (9)–(11), is supposed to be specified on the horizontal squares, at z = ± A , bounding the cube from above and below, whereas boundary values have to be given on the complementary part of the boundary comprised by four vertical squares.
Figure 1. The initial data surface Σ , with a triple black hole system, is chosen to be the cube centered at the origin in R 3 with edges 2 A . The initial data, to the system (9)–(11), is supposed to be specified on the horizontal squares, at z = ± A , bounding the cube from above and below, whereas boundary values have to be given on the complementary part of the boundary comprised by four vertical squares.
Mathematics 09 03170 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rácz, I. Can We Prescribe the Physical Parameters of Multiple Black Holes? Mathematics 2021, 9, 3170. https://doi.org/10.3390/math9243170

AMA Style

Rácz I. Can We Prescribe the Physical Parameters of Multiple Black Holes? Mathematics. 2021; 9(24):3170. https://doi.org/10.3390/math9243170

Chicago/Turabian Style

Rácz, István. 2021. "Can We Prescribe the Physical Parameters of Multiple Black Holes?" Mathematics 9, no. 24: 3170. https://doi.org/10.3390/math9243170

APA Style

Rácz, I. (2021). Can We Prescribe the Physical Parameters of Multiple Black Holes? Mathematics, 9(24), 3170. https://doi.org/10.3390/math9243170

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop