Quantifying closeness between black hole spacetimes: a superspace approach
Abstract
The set of all metrics that can be placed on a given manifold defines an infinite-dimensional ‘superspace’ that can itself be imbued with the structure of a Riemannian manifold. Geodesic distances between points on measure how close two different metrics over are to one another. Restricting our attention to only those metrics that describe physical black holes, these distances may therefore be thought of as measuring the level of geometric similarity between different black hole structures. This allows for a systematic quantification of the extent to which a black hole, possibly arising as an exact solution to a theory of gravity extending general relativity in some way, might be ‘non-Kerr’. In this paper, a detailed construction of a superspace for stationary black holes with an arbitrary number of hairs is carried out. As an example application, we are able to strengthen a recent claim made by Konoplya and Zhidenko about which deviation parameters describing a hypothetical, non-Schwarzschild black hole are likely to be most relevant for astrophysical observables.
I Introduction
The study of geometrodynamics championed by Wheeler wheeler1 ; wheeler2 ; wheeler3 aims to explore, in some precise sense, the configuration space of general relativity (GR). Remarkably, the collection of all metrics that can be placed on some manifold admits a rich geometric structure: the set of metrics that can be placed over can itself be viewed as an (infinite-dimensional) Riemannian manifold gil2 ; clark10 ; ebin ; gil1 . Considering a subspace of which includes only those metrics which are Einstein defines Wheeler’s superspace fishcer ; edwards ; guil09 , and provides a means to quantify the relationship between different spacetime structures in a purely geometrical setting. In particular, geodesics on superspace can be used to ‘measure’ a distance between distinct metrics. More ambitiously, understanding the configuration space of GR may provide a stepping stone to the building of a quantum path integral over spacetime histories hawk ; eucqg .
In GR, all stationary, vacuum, and asymptotically flat black holes (BHs) are locally isometric to the Kerr solution kerr63 ; mars1 . This result is related to the no-hair theorems chand1 , and implies that an appropriately defined configuration space of BHs consists only of various Kerr metrics with different mass and spin parameters (see also Ref. schconfig ). However, if GR is in some sense incomplete and provides but an inexact description for the gravitational field in the strong-field regime, it is possible that more general BHs, described by some finite set of generalised charges or ‘hairs’ (e.g. exact1 ; exact2 ; exact3 ; exact4 ), may exist in Nature. In this paper, we note that this vector of hairs can be used to introduce a coordinate basis over a generalised BH superspace . Some technical considerations necessary to build the superspace, such as the extraction of a single Riemannian manifold from a collection of BH spacetimes, are resolved in detail. Distances on this space, essentially quantifying the difference between two different hairy BHs, can then be used for benchmarking ‘non-Kerrness’; see also Ref. bhag18 .
In general, the nonlinear and (often) higher-than-second order nature of non-GR field equations renders the task of finding exact solutions challenging, and only a handful of physically relevant ones are known (e.g. exact1 ; exact2 ; exact3 ; exact4 ). To address the dearth of exact solutions, various approaches to constructing metrics representing generic BHs in a theory-agnostic manner have been developed nkerr1 ; nkerr2 ; nkerr3 ; contfrac . These metrics are designed to represent parameterised departures from a Kerr description in some particular way, such as including deformations that still preserve the Killing tensor symmetry nkerr5 ; nkerr4 . Recently, Konoplya and Zhidenko (KZ hereafter) kz20 derived a metric which contains “the only parameters that matter”, in the sense that they found including higher-order terms beyond those in their parameterised metric affected electromagnetic observables very little. As an application of the formalism developed here, we further validate the claims made in the aforementioned study by showing that including additional, non-KZ parameters, leads to relatively small curvature differentials on the BH superspace.
This paper is organised as follows. In Section 2 we introduce the concept of a relativistic superspace, and describe how it may be applied to the case of BH geometries. Section 3 then details the necessary technical aspects of the construction (Secs. 3.1 and 3.2), and walks the reader through the various steps for the simple case of static, spherically symmetric BHs (Sec. 3.3). Section 4 is devoted to some worked examples, namely the GR case of the Kerr kerr63 metric (Sec. 4.1) and the theory-agnostic KZ kz20 solution (Sec. 4.2). Some general remarks about the applicability, and limitations, of the approach as a measure for BH closeness are detailed in Section 5. Some discussion is given in Section 6.
II Superspaces
The collection of all metrics definable over a Riemannian manifold can be shown to admit enough structure that it itself defines an (infinite-dimensional) Riemannian manifold ebin . Points on are Riemannian metrics on , i.e. each corresponds to a symmetric, positive-definite -tensor over . We take the base-space to be 3-dimensional (though a generalisation to the higher-dimensional case is straightforward), as later it will be identified with the leaves of a spacelike foliation of a 4-dimensional BH spacetime.
A metric can be introduced over , in the -topology clark10 , as ebin ; gil1
(1) |
where and are tangent vectors to the space of metrics at the ‘reference point’ , provided that the integral converges (see Sec. 3). It is important to note that other metrics exist on , all of which fall under the class of the so-called -metrics osmo which include (1) and the DeWitt gil2 metric as special cases. However, the choice (1) is invariant under the action of the diffeomorphism group on gil1 , and is therefore ‘canonical’ in some appropriate sense. Nevertheless, using a different metric (e.g. DeWitt gil2 ), does not qualitatively change the picture discussed herein.
As it stands, the metric (1) is defined over the entire manifold . In this work we are exclusively interested in those metrics that could describe BHs, so we restrict our consideration to only those and which are tangent vectors to some appropriately defined (see below) space of black hole metrics, . The submanifold inherits a metric from its parent space via pullback, which we also call through a slight abuse of notation. In theories other than GR, it is difficult to define in total generality since, depending on the validity of the no-hair theorem chand1 , there may be an arbitrarily large (but finite) number of parameters (‘hairs’) which characterise the BH (cf. Refs. pappas15 ; suvmel16 ). Nevertheless, suppose that a BH can be described by macroscopic hairs . The parameters can be used to define a natural coordinate basis for the -dimensional submanifold . In vacuum GR, Kerr uniqueness mars1 implies that is a two-dimensional space, spanned by appropriate mass (, say) and spin (, say) vectors. Electric chand1 or more general Yang-Mills ymbh1 ; ymbh2 charges giving rise to additional hairs may also be relevant for electrovacuum BHs in GR.
With respect to the basis, the tensor components of (1) read suv20
(2) |
where . Assuming a (Levi-Civita) connection, all relevant geometric quantities on can be defined from (2), including the Christoffel symbols . The distance between two metrics and , described by parameter vectors and , respectively, is then given by the length of a geodesic connecting these points, viz.
(3) |
for affine parameter , where and , and satisfies the geodesic equation,
(4) |
From the Christoffel symbols , we can also calculate the Ricci curvature
(5) |
whose contraction, , contains additional information about the structure of ; see Sec. 4.
III Stationary, asymptotically flat black holes and 3+1 decompositions
In GR and its extended variants, dynamics traditionally take place on a Lorentzian spacetime rather than on a Riemannian manifold, i.e. the metric is not positive-definite. To construct the space we must first extract a Riemannian manifold from the 4-dimensional pseudo-Riemannian manifold . The standard approach to achieving this is through a Arnowitt-Deser-Misner (ADM) split admsplit ; for a globally hyperbolic spacetime, a foliation of by spacelike (Cauchy) hypersurfaces parameterised by some global time coordinate is always possible ger70 . For stationary spacetimes, the underlying leaves of the foliation are necessarily time-independent in some appropriate sense and we need only consider one such , which will ultimately play the role of our (see Sec. 3.1).
Although we will not present a detailed account of the rich theory of ADM decompositions (which can be found in Refs. gour06 ; gour07 , for example), we note that the 4-dimensional line element, in adapted coordinates where is such that each spatial is constant along its field lines, can be written as
(6) |
where and are known as the lapse function and shift vectors, respectively, and defines a Riemannian metric on each . The metric is just the restriction of to , and will serve as the natural Riemannian metric uniquely associated to any given BH.
III.1 The Riemannian submanifold
A BH is defined by the existence of a horizon chand1 , the interior of which we do not include as part of the spacetime . For a given (stationary) BH geometry, the spatial slice can be identified as , where is a ball of radius centered around some origin for horizon radius ; i.e., one exorcises the interior. This construction suggests a natural candidate for for any given BH spacetime.
However, the volume of is invariably a function of the hairs (e.g. for the Schwarzschild metric of mass , ), which complicates the procedure of choosing a single from a family of BH spacetimes. One way to resolve this problem involves exploiting the diffeomorphism invariance and introducing a ‘horizon-adapted’ radial coordinate (say), so that the horizon is always located at and all choices agree. Physically speaking, however, it is reasonable to consider only BHs with finite horizon area, so that there always exists some largest such , , among any given set of BHs. For the Kerr case, this would correspond to setting a (arbitrarily large) maximum mass that any given BH could attain, and defining . Spatial slices of spacetimes containing smaller BHs can still be defined on this particular in a consistent way through a restriction mapping. For concreteness, we adopt the latter approach here. A schematic representation of how is defined from some ‘maximal’ BH spacetime is shown in Figure 1 (see also Figure 7.1 in Ref. gour07 ).

To pick an , one could note that H emission-line measurements of the ultraluminous quasar TON 618 reveal the existence of a BH with mass shem04 , which is the current record holder for heaviest known object. This suggests that would be a suitable choice, for example, though the geometric picture discussed in this paper is largely unaffected by the exact, numerical value of since arbitrary units of mass can be chosen.
III.2 Ensuring a finite superspace metric
Although we have described a means to construct a universal Riemannian manifold from a given set of BH spacetimes, the construction of cannot be completed unless the superspace metric is well defined. This requires the finiteness of the integral within expression (2).
The Boyer-Lindquist coordinates traditionally used to describe BH spacetimes are singular at the horizon chand1 , and does not exist there. To remedy this, we introduce ingoing null coordinates (also known as Eddington-Finkelstein or Kerr coordinates) gour06 which define a principal null congruence (i.e. curves of constant , constant , and constant are null geodesics). An explicit time coordinate [as needed to build (6)] is then reintroduced through . The metric (and hence ), when expressed in these ingoing null coordinates, is finite over the entirety of (see Secs. 3.3 and 4.1 for some explicit constructions).
There is, however, an additional problem at infinity. Even for asymptotically flat spacetimes, the volume of is unbounded, and the integral defining (going from to infinity) within (2) diverges. This issue may be overcome by introducing a conformal factor into the spacetime metric, which decays fast enough with to ensure the finiteness of the integral (2) over . As such, we introduce a positive function through , with the property that . This defines a conformal 3-metric from (6). Importantly, conformal transformations do not disrupt the null congruence detailed above, since null geodesics for conformally related metrics are identical. In practice, the choice is sufficient, though some investigation suggests that other choices do not introduce any qualitative differences.
III.3 Static black holes
For the case of static, spherically symmetric BHs, we carry out the above procedure in full with an arbitrary metric for demonstrative purposes. In this case, a general spacetime line element, in Boyer-Lindquist coordinates , may be written as
(7) |
where the vanishing of defines the location of the event horizon [i.e. is determined through ]. We begin by introducing the null coordinate through the transformation , where the radial term within the parentheses is sometimes called the tortoise function. An explicit time coordinate is then reintroduced through . In these latter coordinates , the line element (7) becomes
(8) |
which is more complicated than (7) because of the off-diagonal components, though is well behaved at the horizon if the product is finite there. The line element on is then found by directly comparing (8) with (6), viz.
(9) |
Finally, a conformal factor is introduced by rescaling , leading to . Putting everything together, the configuration metric components can be read off from (2),
(10) |
Although relatively simple, the components (10) typically require numerical evaluation given some functions and and BH hairs .
We close this section by noting that for static spacetimes, an alternate procedure to define a BH superspace could be through introducing Wick rotations , as is done in Euclidean quantum gravity schemes eucqg . This would remove the need to introduce ADM splits completely, and one could calculate distances using an appropriately restricted version of (1) defined over the (now) Riemannian 4-space with a Wick-rotated metric. Such an approach will be investigated elsewhere.
IV Illustrative examples
Here we present some explicit examples of the mathematical machinery developed in the previous sections. We focus on the case of vacuum GR (i.e. Kerr superspace; Sec 4.1), and on KZ static spacetimes (Sec. 4.2).
IV.1 General relativity: the Kerr metric
As is well-known, the Kerr metric uniquely represents stable, asymptotically flat BHs in vacuum GR mars1 . The Kerr metric has only two hairs, namely the mass and spin , and so the configuration manifold of BHs in GR is rather simple. Expressed using the ingoing null (Kerr) coordinates described in Sec. 3.2 (see Appendix D of Ref. gour06 for how to relate these to the usual Boyer-Lindquist coordinates), the Kerr line element may be written as kerr63
(11) | ||||
where . The (outer) horizon is located at
(12) |
which is maximal for , and we note that the metric (11) is regular there in any case. Because the horizon surface (12) ceases to be real for , the Kerr (BH) configuration space is only defined for . From (6), we can read off the 3-metric components,
(13) |
which can be used to define the configuration metric (2) after a suitable conformal factor has been introduced. The superspace metric on can now be evaluated from (2) where . For example, we find that the ‘-component’ of the metric reads
(14) |
which can be evaluated using standard numerical techniques without much difficulty after picking some value for . Similar expressions can be obtained for the - and -components for the 2-dimensional metric . The methodology presented above can be adapted to more general stationary metrics in a straightforward way, provided that an appropriate null congruence can be found.
In Figure 2 we show a variety of geodesics over the Kerr configuration manifold for and with an overplotted colour scale that shows the scalar curvature , as computed from (5). The curvature is related to the geodesic distance since, in general, the ratio of volumes between a ball of radius centered at some point on and a corresponding ball on the plane is proportional to the curvature at gray79 , i.e.
(15) |
Although the curvature is relatively weak everywhere, we see that, by comparing the green [joining to ] and black [joining to ] paths (for example), geodesic curves deviate from straight lines more noticeably near the extremal boundary . In particular, as can be inferred from expression (15), the superspace distance separating two near-extremal BHs (e.g. between ones with and ) is greater than the distance separating two slowly rotating holes (e.g. between ones with and ) even though the spin-value differences are the same. This agrees with the physical expectation that stronger spacetime distortions should accompany more rapidly rotating objects, and that a geometric distance measure should reflect this.
IV.2 Static, non-Einstein black holes: Konoplya and Zhidenko metric
As discussed in the Introduction, finding exact solutions describing BHs in theories of gravity aiming to extend GR in some way can be a challenging task. However, progress can still be made in understanding how hypothetical, non-Einstein hairs may manifest in astrophysical data by studying the properties of some parameterised, non-Kerr metric.
For example, a very general BH metric can be constructed by simply assuming Laurent expansions for the metric coefficients with some arbitrary coupling parameters nkerr3 , which are to be constrained by observations. For the static case, one might write the metric potential appearing in (7) as a series, (or something similar), where only negative powers are kept to preserve asymptotic flatness. The , which are to be constrained by observation, together with appropriate appearing in a series for the potential , define the BH hairs. However, practical application usually requires that one truncate (or sum) the series, and so the question becomes: which terms should be kept? Recently, KZ kz20 showed that only a handful of terms [see expressions (16) and (17) below] are necessary to explore the quantitative features of these parameterised metrics for a variety of electromagnetic-based applications. Using the procedure outlined in Sec. 3.2, we can build the configuration manifold of KZ metrics to further test the relative importance of each and , though in a purely geometric setting.
In particular, the KZ metric has four free parameters: (related to the BH mass; in the notation of KZ), , , and , which appear within the metric in some particular combination. We can introduce an additional free parameter, (say), and see how this adjusts the structure of : if distances between metrics with large differences are comparatively small, for instance, this would imply that a non-zero has a weaker geometric impact than the KZ parameters. In particular, consider the metric (7) with
(16) |
and
(17) |
which reduces to the Schwarzschild metric in the limit , and the KZ metric kz20 for . We note that the product is finite at so that (9) is well-behaved, and that the absence of terms in the coefficient ensures that the metric respects post-Newtonian constraints will18 . We may now evaluate the metric (10) where the define coordinates on the 5-dimensional space .
In Figure 3 we show several geodesics (solid lines) together with the scalar curvature (colour scale) , computed by contracting expression (5), on a 2-dimensional slice of varying and with fixed , where we set . The curvature is everywhere negative, and is greatest in magnitude for larger values of horizon radius; . This figure shows that the KZ configuration manifold admits, in general, a non-trivial geometric structure, and that bigger BHs with larger induce greater curvatures on the geometric superspace. Using the formula (3), we find that (for example) the geodesic joining the points and (shown by the black curve) has length , while the flat distance .
For comparison, and to investigate the importance of the non-KZ parameter , we show in Figure 4 a variety of geodesics and the scalar curvature on the slice with varying and but fixed . Overall, the respective scalar curvatures behave in a similar manner, though is generally more negative for the varying case, i.e. .
We find that the distance between the points and (the geodesic connecting these points is shown by the black curve in Fig. 4) is . These endpoints are the same for the black curve shown in Fig. 3, though with replaced by . This calculation reveals an important feature of the KZ metric: the distance between differing values is less than the corresponding distance for different values, and therefore is less important, geometrically speaking, as metrics with different values are closer together. In short, when changing as opposed to , the spaces are in a natural, geometric sense, more distinct. Although not shown, similar conclusions can be drawn by comparing other slices (e.g. the - subspace) or other geodesics on the same slice. The superspace approach presented here therefore supports the conclusions given in Ref. kz20 : parameters other than , , , and appear to matter less.
V Some general remarks about the applicability of the superspace approach
It is important to stress that the localised superspace metric (2) is tailored to measure distances between various members of a parameterised family of metrics. Under these circumstances, distances (3) are relatively straightforward to compute and yield physically reasonable results, at least for the Kerr and KZ cases considered in Sec. 4. However, the general idea may break down when considering theories of gravity that admit separate families of BH solutions with mutually exclusive hairs. In this sense, there may not be an obvious parameterisation that interpolates between two disjoint solution branches. For example, there are no slowly rotating BHs that are simultaneously solutions in the Hořova-Lifshitz and Einstein-aether theories horlif , even though the two are related in the infrared limit and admit the Schwarzschild solution as a limiting case. When considering sets of solutions with mutually exclusive hairs, the interpretation of distances on as measuring BH similarity becomes less clear.
The superspace of metrics also cannot account for non-metric hairs, and may thus fail to capture certain facets of ‘closeness’. As an extreme example, consider two scalar-tensor theories that differ only by the dynamics in their respective scalar sectors. These theories may admit BH solutions described by the same metric, but with different scalar field profiles (due to, e.g., spontaneous scalarization dam93 ; don18 ). These two BH families are not to be thought of as identical, and yet the distance (3) between them on for fixed metric parameters would be zero. In the same way, the metric (2) may also not reasonably quantify differences between solutions in bi- or multi-metric theories bimetric .
VI Discussion
Our ability to experimentally validate GR (or otherwise) in the strong-field regime is continually improving, with now over a dozen measurements of gravitational waves from events involving BHs abb . For this reason, there is interest in developing generic, parameterised BH spacetimes nkerr1 ; nkerr2 ; nkerr3 , with the goal of using these to test for a variety of hypothetical deviations from the Kerr metric in a theory-agnostic way. On the other hand, a high level of generality implies that many of the proposed metrics can be mapped to one another nkerr5 ; nkerr4 , though it is not always obvious how this can be achieved in practice. It is therefore useful to have a tool which can quantitatively distinguish between various parameterised metrics, or indeed between non-Kerr metrics arising as exact solutions in beyond-Einstein theories of gravity. A similar approach for neutron star geometries was recently developed in Ref. suv20 .
Borrowing ideas from the study of geometrodynamics wheeler1 ; wheeler2 ; wheeler3 , we show in this paper how one can build a configuration manifold of BHs; points on this non-flat manifold (see Figs. 2–4) correspond to spatial -metrics associated to ADM-split BH spacetimes. Distances on this manifold can be computed by solving the geodesic equations (4). The construction of a geometric distance between two distinct metrics allows for a natural way to quantify the ‘closeness’ of two BH geometries and could be used, for example, as a general measure for ‘non-Kerrness’; see also Ref. bhag18 . For the particular case of the Kerr configuration manifold (Sec. 4.1), which is 2-dimensional and is spanned by the mass and spin , we found that the Ricci curvature, though remaining finite, develops stronger gradients as the extremal limit is approached (see Fig. 2). The superspace approach therefore reaffirms the expectation that greater spacetime distortions accompany more rapidly rotating compact objects, as distances between BHs with greater values are larger than distances between slowly rotating holes.
For static BHs, we further validate the claims made by Konoplya and Zhidenko kz20 about which parameters are most important (that is, those which modulate astrophysical observables the most) when constructing a parameterised BH metric. This is achieved by noting that distances between (generalised) KZ metrics with an additional free parameter and those used in Ref. kz20 are comparatively small (see Sec. 4.2), which means that these spaces are geometrically similar.
Although we have only considered asymptotically flat BHs here, it is likely that the approach can be generalised to the asymptotically de Sitter case without much difficulty, since the conformal factor can be used to tune the convergence of (2) even when grows as . The inclusion of non-stationary black holes, where time can be thought of as an additional hair of sorts, could also be achieved given a clear choice for a single among the family of metrics under consideration. Finally, as discussed in Sec. 5, we note that the approach presented here may not capture all features that might be expected of BH ‘closeness’. Importantly, the localised metric (2) on the superspace requires a parameterised family of BH metrics . It is not obvious how to construct such a parameterisation for theories admitting disjoint branches of solutions with mutually exclusive hairs horlif ; don18 . Purely non-metric charges are also not accounted for in expression (2), as non-metric (e.g. scalar) fields do not appear.
Acknowledgements
We thank Prof. Bill Moran for hospitality shown during a research visit to Melbourne, where some of this work was completed. We are grateful to Sebastian Völkel for discussions. We thank the anonymous referees for providing feedback which helped improve the clarity of this manuscript. This work was supported by the Alexander von Humboldt Foundation.
References
References
- (1) J. A. Wheeler, Einsteins Vision, Springer Verlag, Berlin, 1968
- (2) J. A. Wheeler, in Battelle Recontres, Editors, DeWitt and Wheeler, W. A. Benjamin, Inc., 1968.
- (3) J. A. Wheeler, Annals of Physics 2, 604 (1957).
- (4) B. S. DeWitt, Phys. Rev. 160, 1113 (1967).
- (5) B. Clarke, Calculus of Variations and Partial Differential Equations 39, 533 (2010).
- (6) D. G. Ebin, Proc. Symp. Pure Math. AMS 15, 11 (1970).
- (7) O. Gil-Medrano and P. W. Michor, The Quarterly Journal of Mathematics 42, 183 (1991).
- (8) A. Fischer, The theory of superspace, in Proceedings of the Relativity Conference, edited by M. Carmeli, S. I. Fickler, and L. Witten (Plenum, New York, 1970).
- (9) D. A. Edwards, Studies in topology. Academic Press, 121, (1975).
- (10) D. Giulini, Gen. Rel. Grav. 41, 785 (2009).
- (11) S. W. Hawking, Comm. Math. Phys. 55, 133 (1977).
- (12) G. W. Gibbons and S. W. Hawking, Euclidean Quantum Gravity, (World Scientific, Singapore, 1993).
- (13) R. P. Kerr, Phys. Rev. Lett. 11, 237 (1963).
- (14) M. Mars, Class. Quant. Grav. 17, 3353 (2000).
- (15) S. Chandrasekhar, The Mathematical Theory of Black Holes, (Oxford University Press, Oxford, 1998).
- (16) K. V. Kuchar̆, Phys. Rev. D 50, 3961 (1994).
- (17) D. Grumiller and N. Yunes, Phys. Rev. D 77, 044015 (2008).
- (18) T. Clifton and J. D. Barrow, Phys. Rev. D 72, 103005 (2005).
- (19) C. Charmousis et al., Phys. Rev. D 100, 084020 (2019).
- (20) R. A. Konoplya and A. Zhidenko, Phys. Rev. D 101, 084038 (2020).
- (21) S. Bhagwat et al., Phys. Rev. D 97, 104065 (2018).
- (22) K. Glampedakis and S. Babak, Class. Quant. Grav. 23, 4167 (2006).
- (23) S. Vigeland, N. Yunes, and L. C. Stein, Phys. Rev. D 83, 104027 (2011).
- (24) T. Johannsen and D. Psaltis, Phys. Rev. D 83, 124015 (2011).
- (25) L. Rezzolla and A. Zhidenko, Phys. Rev. D 90, 084009 (2014).
- (26) G. O. Papadopoulos and K. D. Kokkotas, Class. Quant. Grav. 35, 185014 (2018).
- (27) Z. Carson and K. Yagi, Phys. Rev. D 101, 084030 (2020).
- (28) R. A. Konoplya and A. Zhidenko, Phys. Rev. D 101, 124004 (2020).
- (29) O. Pekonen, J. Geom. Phys. 4, 493 (1987).
- (30) G. Pappas and T. P. Sotiriou, Phys. Rev. D 91, 044011 (2015).
- (31) A. G. Suvorov and A. Melatos, Phys. Rev. D 93, 024004 (2016).
- (32) R. Bartnik and J. McKinnon, Phys. Rev. Lett. 61, 141 (1988).
- (33) P. Bizon, Phys. Rev. Lett. 64, 2844 (1990).
- (34) A. G. Suvorov, Class. Quant. Grav. 37, 065013 (2020).
- (35) R. L. Arnowitt, S. D. Deser, and C. W. Misner, Gen. Rel. Grav. 40, 1997 (1962).
- (36) R. Geroch, J. Math. Phys. 11, 437 (1970).
- (37) E. Gourgoulhon and J. L. Jaramillo, Physics Reports 423, 159 (2006).
- (38) E. Gourgoulhon, Lecture Notes in Physics, (Springer Verlag, Berlin, 2012), Vol. 846
- (39) O. Shemmer et al., Astrophys. J. 614, 547 (2004).
- (40) A. Gray and L. Vanhecke, Acta Mathematica 142, 157 (1979).
- (41) C. M. Will, Theory and experiment in gravitational physics, (Cambridge university press, Cambridge, 2018).
- (42) E. Barausse and T. P. Sotiriou, Phys. Rev. D 87, 087504 (2013).
- (43) T. Damour and G. Esposito-Farése, Phys. Rev. D 54, 1474 (1996).
- (44) D. D. Doneva and S. S. Yazadjiev, Phys. Rev. Lett. 120, 131103 (2018).
- (45) E. Babichev and R. Brito, Class. Quant. Grav. 32, 154001 (2015).
- (46) B. P. Abbott et al., Phys. Rev. X 9, 031040 (2019).