A family of scalar-Einstein and vacuum Brans-Dicke axisymmetric solutions
Bertrand Chauvineau

TL;DR
This paper introduces a novel method linking static axially symmetric vacuum solutions in general relativity to a class of scalar-Einstein solutions with specific mass and angular momentum, revealing new Kerr-like solutions in Brans-Dicke theory.
Contribution
It establishes a one-to-one correspondence between vacuum GR solutions and scalar-Einstein solutions, and constructs new stationary axially symmetric solutions in Brans-Dicke theory.
Findings
A new method for solution correspondence is proposed.
A specific asymptotically flat solution with a ringlike naked singularity is analyzed.
The solution reduces to Kerr when scalar charge vanishes but not to Fisher-Janis-Newman-Winicour for zero rotation.
Abstract
A new method is proposed, that establishes a one to one correspondance between the whole set of static axially symmetric vacuum GR solutions and a specific class of stationary axially symmetric scalar-Einstein solutions having a given mass and a given angular momentum. The method explicitly takes advantage of the Kerr metric Ricci flatness. This also results in a class of stationary axially symmetric vacuum, ie Kerrlike, Brans-Dicke solutions. A particular solution, that is asymptotically flat, is more closely considered. It converges to Kerr for a vanishing scalar charge, but fails to converge to the Fisher-Janis-Newman-Winicour solution for a vanishing "rotation parameter". This solution exhibits a naked singularity having a ringlike structure.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
A family of scalar-Einstein and vacuum Brans-Dicke axisymmetric
solutions
Bertrand Chauvineau
(Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange, France
(e-mail : [email protected]))
Abstract
A new method is proposed, that establishes a one to one correspondance between the whole set of static axially symmetric vacuum GR solutions and a specific class of stationary axially symmetric scalar-Einstein () solutions having a given mass and a given angular momentum. The method explicitly takes advantage of the Kerr metric Ricci flatness. This also results in a class of stationary axially symmetric vacuum, ie Kerrlike, Brans-Dicke solutions. A particular solution, that is asymptotically flat, is more closely considered. It converges to Kerr for a vanishing scalar charge, but fails to converge to the Fisher-Janis-Newman-Winicour solution for a vanishing ”rotation parameter”. This solution exhibits a naked singularity having a ringlike structure.
I- Introduction
Obtaining exact solutions of a relativistic gravity theory is far from being an easy task, due to the highly non-linear character of the involved field equations. The most known solutions are Schwarzschild, Kerr, and Friedmann-Lemaitre-Robertson-Walker (FLRW), that solve the general relativity (GR) field equation in vacuum for the two first ones. These three solutions (or family of solutions for FLRW) are of obvious usefulness in the astronomical framework. Some other exact GR solutions are known, that are also very useful for astronomical purposes [1][2][3]. On the other hand, an incredibly huge number of solutions, the usefulness of which is not obvious a priori, have been obtained so far [4]. However, despite the non immediate relevance of a solution, its usefulness should not be underestimated. Indeed, any solution may reveal some unexpected property of the considered theory. It may also happen that a solution only later turns out to be of genuine astronomical relevance: the emergence of the black hole (BH) concept from the Schwarzschild solution is probably the most obvious exemple.
Since many attemps to quantify gravity, and/or to unify gravity with other interactions, return a classical gravitationnal sector having not a GR, but a scalar-tensor (ST) structure [5][6][7], Brans-Dicke (BD) and ST gravity theories are considered as valuable alternatives to GR, despite the fact that the latter successfully passes solar system tests up to now [8]. (Let us also point out that many ST theories are driven to mimic the GR behaviour as a consequence of the cosmic expansion [9][10].) Thence the interest in ST theories, and specifically in BD, that just involves a constant parameter instead of an arbitrary function . In this context, looking for exact BD/ST solutions is particularly appropriate to point out relevant qualitative features with respect to GR. For instance, the spherical Brans class I solution [11][12] generally exhibits a naked singularity (NakS) or wormhole structure [13] (see also [14] for pionering works on the detailled significance of the Class I, II, III and IV Brans solutions), such features being absent from the Schwarzschild GR solution. Besides, it has been recently shown that a particle orbiting a large Brans class I solution results in an observed (by a far observer) unbound orbital frequency, depending on how much the solution is scalarized [15]. One may then suspect striking qualitative differences in extreme mass ratio binaries BD/ST gravitational radiation, with respect to GR, since GR orbital frequencies can not exceed the innermost circular orbital value. Nevertheless, let us remind that under some conditions (mainly regularity, asymptotical flatness and finite area horizon), vacuum and stationary BH like solutions are the same in BD/ST as in GR [16][17].
It is known from long that, by the means of a conformal transformation, any vacuum BD solution is associated to a massless scalar filled GR solution [18]. Thence seeking vacuum BD solutions can be reformulated as a scalar-GR problem. From the Hawking’s theorem [16], it is clear that any stationary axisymmetric (SAS) vacuum BD, but non Kerr, solution should exhibit a NakS structure, unless exhibiting some peculiar feature that allow it to evade the theorem (like being not asymptotically flat).
A method has been proposed by [19], that allows to generate a BD vacuum SAS solution from a GR vacuum SAS seed one, provided one is able to solve a given non linear PDE system. A class of electrovacuum BD solutions has been obtained by [20], that generalises the Majumdar-Papapetrou GR solution. The same authors later proposed a method that allows to generate SAS vacuum BD solutions from SAS vacuum GR ones (and also from static axisymmetric vacuum BD ones) [21], but the solutions then obtained are generally not asymptotically flat. The SAS, but only one coordinate dependent, case has been solved by [22]. The case of ST SAS solutions has been considered by [23], and explicit solutions are given in a few case of very specific ST (not BD) theories. Along the same lines as [19], it is claimed in [24] that any SAS vacuum BD solution can be obtained by nonlinearly combining any SAS vacuum GR solution and any vacuum solution of the Weyl class. A Kerr-type BD solution has been obtained in [25], that was built in such a way that it reduces to the Kerr solution in the limit. Inspired by [20] and [22] works, a BD BH-like solution is proposed by [26], but this solution is not asymptotically flat. Motivated by particle collisions near a Kerr-like BD BH, another non asymptotically flat solution is obtained in [27], that is also derived using the [20] method. Starting from an unusual formulation of the Lewis metric, a two parameters extension of the method initiated in [21] allowed [28] to derive a BD version of the Ernst equation. The method is applied to some examples, but here again, the obtained spacetimes are generally not asymptotically flat. The matter filled case has been considered by [29] in self interacting (ie with a potential ) BD, but only static axisymmetric spacetimes were considered. Very recently, a way to generate new solutions starting from known ones has been proposed by [30]. The technique makes use of a BD symmetry in the traceless case ().
It is worth reminding that a scalar-metric was proposed by [31] as an SAS vacuum BD solution. Its conformally related scalar-GR version is given in [32]. This metric, or its [32] form, has been used in several papers to characterize physics in a NakS (versus BH) field [32][33][34][35], and also to suggest rotating antiscalar solutions as alternatives to Kerr BHs [36]. However, the [31] scalar-metric is actually not a vacuum BD solution, as it is explicitly shown in [37]. The reason is that the authors of [31] applied, without justification, to a BD spherical vacuum solution (the Brans Class I) the Newman-Janis (NJ) algorithm, the authors (NJ) having just noticed it to be an unexpected way to recover Kerr from Schwarzschild [38]. Let us stress that, even in the GR framework, determining the ability of the NJ algorithm (or some equivalent reformulations) to generate new solutions from a seed one is not an easy task [39]. The NJ algorithm, as well as some modified versions, nevertheless received continuous interest, not only as a way to suggest new (generally non perfect fluid filled) solutions, but also in the context of other theories, like supergravity. See for instance [40] for a recent review. Related to the NJ finding [38], let us also mention that another (but fully justified) way to derive Kerr from Schwarzschild was obtained in [41].
In this paper, a new one to one correspondance between static axially symmetric vacuum GR solutions and SAS massless scalar GR solutions, with metrics having some prior form, is established. This prior form is inspired from the metric found in [27], refered to as the SB solution in the following. It is defined as a modification of the Kerr metric, with given parameters, by inserting two unknown metric functions in a suitable manner. Thence these functions are demanded to be such that the scalar-Einstein field equations are satisfied. Thence, it uses Kerr as a seed in some sense, but in a way that differs from the generating techniques previously reviewed. The definition of some well suited coordinates then allows to establish the correspondance. Let us stress that, unlike the NJ algorithm, the method is completely justified since it involves an explicit solving of the relevant field equations. The SB solution is recovered as a special case, and other explicit solutions are built, one of them being asymptotically flat. In some sense, the introduction of the coordinates cures the problem of solving the non linear PDE system obtained in [19], for the considered prior metric form.
Outline of the paper
The prior form of the metric and the related index notations are defined in section II. The first half of the Einstein equations is considered in II-1. The Klein-Gordon (KG) equation is used in II-2, that suggests the definition of the coordinates. In II-3, the second half of the Einstein equations is considered. The correspondance with the general static axisymmetric vacuum GR case is explicited in II-4. A particular asymptotically flat solution is then considered in section III. The links with the BD theory is discussed in section IV, while section V is dedicated to a brief conclusion.
II- The considered set of axisymmetric metrics
We consider in this paper metrics having the form
[TABLE]
[TABLE]
in Boyer-Lindquist coordinates, and where and are dependent functions. One has introduced the usual quantities
[TABLE]
Besides the usual index convention , let us also make the convention
[TABLE]
while the indexes take the four spacetime values. The ordering in (0.2) explicits the block diagonal structure of Kerr’s metric. The metric (0.1a-0.1b) is also block diagonal, each block being ”conformally” related to the Kerr corresponding one, but with different ”conformal” factors. The requirement (0.1a-0.1b) implicitly means imposing two prior relations between the four metric functions describing the general form of an SAS metric (see for instance eq. (1) of [42]).
The metric (0.1a-0.1b) is required to solve the scalar-Einstein equation
[TABLE]
where is an dependent scalar field. The special case solves (0.4) for , since the Kerr metric (0.2) is Ricci flat.
The SB solution [27] corresponds to and . Of course it should be recovered as a special solution of (0.4) with the corresponding scalar, that reads . This will be checked later.
From the dependence of the considered scalar-metric, one has . The non zero connexion components then read
[TABLE]
[TABLE]
[TABLE]
For convenience, let us introduce the following notation
[TABLE]
II-1- The Einstein equation components
From (0.4), one has , that writes
[TABLE]
Using (0.5a), (0.5b), but also (Kerr’s metric being Ricci flat) and (from the Ricci identity on Kerr’s metric), (0.9) yields an equation that is only dependent
[TABLE]
Eliminating the bracket term thanks to the contraction by , one obtains (since )
[TABLE]
Using the Kerr metric (0.2) and the connexion components (0.5b), (0.11) yields
[TABLE]
This equation has to be satisfied by both , and . Writing out the two equations for and , one obtains two homogeneous equations on and . One then shows that having a non trivial solution requires . Thence, is necessarily constant for . One can then set by redefining and the units.
The fact that has two straigtforward consequences: (1) the equatorial () circular orbits and their linear planar stability (both do not involve ) are obtained solving the same equations as the Kerr’s case, and (2) the horizon and ergosphere are the ”same” as Kerr’s
[TABLE]
[TABLE]
To be precise, let us point out that these claims just concern the functions that describe these orbits and surfaces in terms of and of the parameters. Indeed, the metric components enter their geometric and relative properties. For instance, the geometric radial distance between two circular orbits having circumferences and depends on .
Since , the metric (0.1a-0.1b) achieves the form
[TABLE]
II-2- The Klein-Gordon equation
To pursue the integration, it would be sufficient to solve the components of (0.4), since the KG equation is a direct consequence of (0.4). It is nevertheless useful to write out the KG equation
[TABLE]
Using (0.15), it reads
[TABLE]
This form suggests defining the alternative radial coordinate
[TABLE]
that yields
[TABLE]
[TABLE]
For convenience, one has introduced the following notation
[TABLE]
One sees that this form returns two obvious solutions
[TABLE]
[TABLE]
In terms of the initial Boyer-Lindquist coordinates, these coordinates read
[TABLE]
[TABLE]
that turn out to be useful in the calculations to do. The derivation operators transform as
[TABLE]
Reinserting in (0.20) returns, after a lengthy calculation, and using (0.25), the KG equation in the nice form
[TABLE]
This form points out two other obvious solutions, besides (0.22a) and (0.22b)
[TABLE]
It may be worth spoting that if solves (0.27), so do its successive derivatives with respect to . (Remark that .) Along these lines, new solutions can also be obtained by integration with respect to . For instance
[TABLE]
See the appendix for a quicker demonstration of (0.27), using the coordinates.
II-3- The Einstein equation components
From (0.4), one has , that writes
[TABLE]
Using (0.5a), (0.5c) and (0.7), with , but also (Kerr’s metric is Ricci flat) and (from the Ricci identity on Kerr’s metric), (0.30) yields
[TABLE]
Expliciting the , and components yields, using the radial coordinate and (0.19a)
[TABLE]
[TABLE]
See the appendix for a quicker demonstration of these results, using the coordinates from the start.
II-4- Generating solutions
From the previous sections, a solution can be obtained by (1) solving first the KG equation (0.27), and then (2) integrating the system (0.33a)-(0.33b). There is then a direct correspondance with the issue of seeking a static axisymmetric solution of vacuum GR. Indeed, such a metric can be written [2]
[TABLE]
where the metric functions and have to solve
[TABLE]
Considering the four first solutions of the KG equation obtained in II-2, one obtains
[TABLE]
From the previously pointed out correspondance, the SB solution is associated to the Minkowski spacetime (in some non cartesian coordinates), while the Newtonian one is associated to the Curzon-Chazy spacetime [2].
It is worth also having a look on the fifth solution of the KG equation obtained in II-2. Integrating for yields
[TABLE]
ie, in terms of coordinates, using (0.23)
[TABLE]
This, with (0.15), suggests that the spacetime is not asymptotically flat. From the previous correspondance, the solution is associated to the Gautreau-Hoffman spacetime [2].
The solution is obviously not asymptotically flat. Indeed, (0.4) shows that it is even not Ricci flat at infinity, since and do not vanish at infinity.
III- The Newtonian solution
Since it is asymptotically flat, let us have a closer look on the Newtonian solution. Its metric (0.36d) explicitly reads
[TABLE]
The and metric tensor components being identical to Kerr, the and constants have the same ADM mass and angular momentum meanings. As already mentionned, the Kerr horizon (0.13) and ergosphere (0.14) surfaces are recovered, since they just depend on the metric components, which are the same as Kerr.
The scalar associated to (0.39) reads, in terms of coordinates
[TABLE]
The explicit dependence of shows that the (0.39)-(0.40) solution is not the BD-Kerr solution obtained in [19] (equations (24)-(26) of [19]).
III-1- Singularities
One knows that the Kerr metric is singularity free outside its (external) horizon (0.13), and that the horizon itself is a regular surface. However, from the Hawking theorem [16], this should not be true for the spacetime (0.39) since while (1) it is asymptotically flat and, (2) its horizon surface is finite. Indeed, the scalar curvature reads, from (0.4)
[TABLE]
which, from (0.36d), yields
[TABLE]
One sees that (0.42) diverges whatever the way , ie the way (and only in this case for ). The points having and are then NakS points. Thence, the horizon is not regular everywhere, since its equator , and the equator only (considering points having ), is scalar curvature singular.
However, the fact that the scalar curvature does not diverge on the points having is not sufficient to prove that none of these points is singular. A quantity that is often regarded as a reliable singularity indicator is the Kretschmann scalar
[TABLE]
where is the Riemann-Christoffel curvature tensor. From (0.5a)-(0.5c), one finds that
[TABLE]
where is the Kerr’s metric Riemann-Christoffel curvature tensor, and
[TABLE]
If , which implies from (0.7), is finite since the Kerr’s metric is regular in the considered region. Any divergence of can then only appear from the and , ie from the and , quantities entering the terms in (0.44). It is then easy to see from (0.36d) that a divergence of can only occur at points where . The set of Kretschmann singularity points is then included in the scalar curvature points set. This strongly suggests that the horizon is regular at points not belonging to the equatorial NakS ring 111Extending to the region, one sees on (0.42) that is a scalar curvature singularity of the considered spacetime (it is a singularity, but not a scalar curvature one, in Kerr’s spacetime). The spacetime has then two ring singularities: one is naked, the second being hidden (the latter being the counterpart of the usual Kerr’s, in some sense)..
Let us stress that it was recently found by [43] that the (not asymptotically flat) SB solution, ie (0.36a), also exhibits such a ringlike NakS structure.
III-2- The subcase (static case)
Let us now specify the Newtonian solution (0.36d) to the case. The vacuum Kerr’s solution returns the (vacuum) spherical Schwarzschild solution for . In the non vacuum, but massless scalar filled, case, a spherical solution is known, often named the JNW metric, refering to the Janis-Newman-Winicour 1968’s paper [44]. It is worth to point out that this solution was in fact earlier discovered by Fisher in his 1948 paper [45], but using an areal radial coordinate. It seems then fair to use the FJNW acronym when refering to this solution, and so will I do in this paper. The FJNW solution includes Schwarzschild as a limit (non scalarized) case. It is then naturel to suspect the existence of a massless scalar filled GR solution, that would depend on a ”rotation parameter” , and that would return (1) the Kerr spacetime for a vanishing scalar, and (2) the FJNW solution for .
The solution (0.36d) indeed returns the Kerr spacetime for a vanishing scalar, ie for . On the other hand, making in (0.39) does not return the FJNW solution, but
[TABLE]
This static but non spherical solution is known from long, see equations (17)-(18) of Penney’s paper [46]. From (0.40), the scalar field reads
[TABLE]
The metric (0.46) has an horizon, that reads
[TABLE]
and whose equator’s points , and only these points, are singular, as it can be directly checked from the expression entering (0.42), that reads for
[TABLE]
This agrees with the fact that the metric, while static since , is not spherical because of the presence of the ( dependent) exponential in front of the bracket in (0.46). Incidentally, this also shows that (0.46) can not be FJNW in deguise. The ring character of the NakS, stuck on the horizon, is obviously a property inherited from the general case.
The fact that the FJNW metric does not appear as a subcase of (0.46) means that (0.39) can not be interpreted as a ”rotating version of FJNW”. It also strongly suggests the existence of asymptotically flat Kerr like solutions of (0.4), that do not fulfil (0.1a). Indeed, recovering FJNW as a rotationless limit case, with its properties ([15]), is incompatible with the conclusion obtained in II-1.
III-3- Orbits in the equatorial plane of the static solution
As already mentioned, the presence of the scalar field affects neither the existence condition of equatorial circular orbits nor their linear stability (using the radial coordinate defined by the form of the metric (0.15)). The situation is very similar to the case of a scalarized version of the -metric reported in [47]. However, this does not mean that the physics is unaffected by the scalar, even in the equatorial plane. The aim of this subsection is to illustrate this point. For convenience, we specify to the static case, that makes all the calculations easily tractable.
In the plane, the metric (0.46) simplifies into
[TABLE]
One has already pointed out that the exponential term in impacts the distance between close circular orbits. The effect is mostly important near the NakS, where the distance between two close orbits of radii and goes to zero, because of the vanishing of the exponential factor. Related to this, the (coordinate) time needed for a photon to radially propagate from to an external observer at coordinate , that reads
[TABLE]
converges for , unlike what happens in the Schwarzschild case. On the other hand, the convergence does not occur for a radial polar photon ( or ), since the argument of the exponential term in (0.46) cancels then. More generally, the convergence does not occur for any photon crossing the Penney’s sphere at any non equatorial point, since then , so that the exponential goes to when . However, let us also spot that despite the convergence of (0.51), the far observed behaviour of a clock at rest close to the Penney’s sphere is frozen, even in the equatorial plane, since for any , the case included. The consequence is that in the equatorial plane, while circular close to the NakS (non geodesic) motions are frozen (the areal NakS circumference being ), radial infallings are not.
Let us point out that the same observations, ie (1) (coordinate) time convergence for a radial photon propagating from any NakS point, but (2) frozen rest clock behaviour near any NakS point, also occurs in the FJNW spherical metric [44][45]. This is obvious from the isotropic form of this metric [15]
[TABLE]
where ( corresponding to the GR Schwarzschild metric)222Let us point out that in the FJNW case, the fact that close to the NakS rest clocks are frozen does not mean that circular motions are frozen, since the NakS circumference is zero. Indeed, close to NakS circular geodesics (that exist for ) return a divergent far observed frequency [15].. At least for these two metrics, the time convergence for a photon reaching the horizon only concerns orbits approaching points belonging to the NakS location.
IV- Brans-Dicke vacuum solutions
It is well-known [48] that in a four-dimensional spacetime, the conformal transformation
[TABLE]
yields, for any constant (supposed to be )
[TABLE]
where is any positive scalar function and
[TABLE]
This means that the vacuum BD action of the BD gravitational field identifies with the GR action, with gravitational field , but filled by the (matter source) massless scalar . Thence, any solution of (0.4) is conformally associated to a vacuum BD solution [18]. This ()BD solution reads
[TABLE]
Experiments strongly constrain BD/ST theories to satisfy , where is the BD parameter, or the present value of in the ST case [8]. In such circumstances, vacuum BD gravity, and also to some extent ST gravity, is asymptotically equivalent to massless scalar filled GR [15][28][49]. In other words, the solutions obtained solving (0.27) and (0.33a)-(0.33b) can directly serve as vacuum BD/ST solutions in the large case, without having to explicitly consider the conformal correspondance (0.56a)-(0.56b). Linked to this, let us remind that for any scalar function chosen ”independently on ”, (0.56a) yields
[TABLE]
Despite that goes to a constant value, the term entering the full BD equation does then not vanish in the large limit, but results in a contribution. This is coherent with the fact that a vacuum BD solution does not reduce to a GR vacuum solution (that would have been Kerr in the SAS case) in the ” limit” [49].
Let us remark that while the massive geodesics are not the same in the scalar-Einstein and in the corresponding BD solutions, they are asymptotically identical in the limit, since the conformal factor is constant in this limit.
V- Conclusion and outlook
One has established a non trivial one to one correspondance that allows to built an SAS massless scalar-GR solution from any static axisymmetric vacuum GR solution. While this in turn also results in a way to obtain an SAS vacuum BD solution, one has pointed out that the so obtained massless scalar-GR solutions can also directly serve as a (large ) vacuum BD/ST solutions for practical purposes.
It was claimed in [25] that a Kerr-like BD solution should: (1) only depend on the three parameters, (2) go to Kerr spacetime in the limit, and (3) return Schwarzschild’s spacetime for . From the results obtained in this paper, it seems that the truth is by far more complex. The solutions explicitly presented in this paper can be interpreted as limit of BD solutions for , but they differ from Kerr spacetime. The particular case (0.46) of the Newtonian solution has but differs from Schwarzschild. The fact that the Newtonian solution (0.39), for instance, depends on a parameter , besides , shows that the ()BD solution built from it using (0.56a)-(0.56b) results in a family of Kerr-like solutions that depends on more than the three parameters.
Let us mention that a lot of authors explicitly or implicitly still suppose that (2) is true (see for instance [26] or [27]). Such a presupposition results in a drastic impoverishment of the (large ) ST potential predictions. Among these the drastically different behaviour of far observed orbital frequencies in the spherical case with respect to GR [15]. It would be of interest to know whether such a different from GR behaviour could also occur in the rotating case. This requires exploring more general SAS solutions of (0.4) with , which means ruling out the restrictive hypothesis (0.1a).
Acknowledgments
I would like to thank here an anonymous referee for several very interesting historical references, especially concerning the often named JNW (but accordingly here renamed FJNW) scalar-Einstein solution.
Appendix: calculations using the coordinates
Using (0.18) and (0.23), one obtains, in the coordinates
[TABLE]
where
[TABLE]
From (0.23) and the trigonometric identities, one could explicitly obtain in terms of , but this is not needed here. Indeed, since , the Dalembertian operator reads
[TABLE]
from which one directly obtains (0.27).
Writing now (0.58) in the form
[TABLE]
(the version of which being Kerr in coordinates), the connexion components achieve a form like (0.5a)-(0.5c) with , but with instead of . Thence, (0.31) is replaced by
[TABLE]
with . Writing out the component directly returns (0.33b). The and components read
[TABLE]
References
[1] A. Krasiński, Inhomogeneous Cosmological Models, Cambridge University Press (1997).
[2] J. B. Griffiths, J. Podolský, Exact Space-Times in Einstein’s General Relativity, (Cambridge University Press, 2009).
[3] K. Bolejko, A. Krasiński, C. Hellaby, M.-N. Célérier, Structures in the Universe by exact methods (Cambridge University Press, 2010).
[4] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, E. Herlt, Exact solutions of Einstein’s field equations (Cambridge University Press, 2003).
[5] Y. Fujii, K. Maeda, The scalar-tensor theory of gravitation (Cambridge University Press, 2003).
[6] V. Faraoni, Cosmology in scalar-tensor gravity (Kluwer Academic Publishers, 2004).
[7] S. Capozziello, V. Faraoni, Beyond Einstein’s gravity, Fundamental Theories of Physics, volume 170, Springer (2011).
[8] C. M. Will, *The confrontation between general relativity and experiments *in www.livingreviews.org/Irr-2014-4 (living reviews in relativity).
[9] T. Damour, K. Nordtvedt, Phys. Rev. Lett. 70, 2217 (1993).
[10] T. Damour, K. Nordtvedt, Phys. Rev. D 48, 3436 (1993).
[11] C. Brans and R. H. Dicke, Phys. Rev. 124, 925 (1961).
[12] C. H. Brans, Phys. Rev. 125, 2194 (1962).
[13] V. Faraoni, F. Hammad, S. D. Belknap-Keet, Phys. Rev. D 94, 104019 (2016).
[14] K. A. Bronnikov, Acta Phys. Pol. B4, 251 (1973).
[15] B. Chauvineau, Gen. Rel. Grav. 49, 143 (2017).
[16] S. W. Hawking, Commun. Math. Phys. 25, 167 (1972).
[17] J. D. Bekenstein and A. Meisels, Phys. Rev. D 18, 4378 (1978).
[18] R. H. Dicke, Phys. Rev. 125, 2163 (1962).
[19] C. B. G. McIntosh, Commun. Math. Phys. 37, 335 (1974).
[20] R. N. Tiwari, B. K. Nayak, Phys. Rev. D 14, 2502 (1976).
[21] B. K. Nayak, R. N. Tiwari, J. Math. Phys. 18, 289 (1977).
[22] T. Singh, L. N. Rai, J. Math. Phys. 22, 136 (1981).
[23] T. Singh, T. Singh, Astrophys. Sp. Sc. 100, 309 (1984).
[24] A. García N. Bretón, I. Hauser, Int. J. Theor. Phys. 27, 635 (1988).
[25] G.-W. Ma, Int. J. Theor. Phys. 34, 2331 (1995).
[26] H. Kim, Phys. Rev. D 60, 024001 (1999).
[27] J. Sultana, B. Bose, Phys. Rev. D 92, 104022 (2015).
[28] P. Kirezli, Ö. Delice, Phys. Rev. D 92, 104045 (2015).
[29] M. Sharif, R. Manzoor, Commun. Theor. Phys. 68, 39 (2017).
[30] V. Faraoni, D. K. Çiftci, S. D. Belknap-Keet, Phys. Rev. D 97, 064004 (2018).
[31] K. D. Krori, D. R. Bhattacharjee, J. Math. Phys. 23, 637 (1982).
[32] K. K. Nandi, P. M. Alsing, J. C. Evans, T. B. Nayak, Phys. Rev. D 63, 084027 (2001).
[33] G. N. Gyulchev, S. S. Yazadjiev, Phys. Rev. D 78, 083004 (2008).
[34] Z. Kovács, T. Harko, Phys. Rev. D 82, 124047 (2010).
[35] T. Karmakar, T. Sarkar, Gen. Rel. Grav. 50, 85 (2018).
[36] M. A. Makukov, E. G. Mychelkin, Phys. Rev. D 98, 064050 (2018).
[37] B. Chauvineau, Phys. Rev. D 98, 088501 (2018).
[38] E. T. Newman, A. I. Janis, J. Math. Phys. 6, 915 (1965).
[39] H. Erbin, Gen. Rel. Grav. 47, 19 (2015).
[40] H. Erbin, Universe 3, 19 (2017).
[41] G. Clément, Phys. Rev. D 57, 4885 (1998).
[42] H. Levy, W. J. Robinson, Proc. Camb. Phil. Soc. 60, 279 (1964).
[43] Z. Kovács, T. Harko, S. Shahidi, Phys. Rev. D 98, 088502 (2018).
[44] A. I. Janis, E. T. Newman, J. Winicour, Phys. Rev. Lett. 20, 878 (1968).
[45] I. Z. Fisher, Z. Exp. Teor. Fiz. 18, 636 (1948); see also arXiv:gr-qc/9911008 for an English translation.
[46] R. Penney, Phys. Rev. 174, 1578 (1968).
[47] B. Turimov, B. Ahmedov, M. Kološ, Z. Stuchlík, arXiv:gr-qc/1810.01460 (2018).
[48] R. V. Wagoner, Phys. Rev. D 1, 3209 (1970).
[49] B. Chauvineau, Gen. Rel. Grav. 39, 297 (2007).
