Ricci Flow Gravity
PMC Physics A volume 1, Article number: 3 (2007)
A theory of gravitation is proposed, modeled after the notion of a Ricci flow. In addition to the metric an independent volume enters as a fundamental geometric structure. Einstein gravity is included as a limiting case. Despite being a scalar-tensor theory the coupling to matter is different from Jordan-Brans-Dicke gravity. In particular there is no adjustable coupling constant. For the solar system the effects of Ricci flow gravity cannot be distinguished from Einstein gravity and therefore it passes all classical tests. However for cosmology significant deviations from standard Einstein cosmology will appear.
PACS Codes: 04.50.+h, 04.20.Cv, 02.40.Hw
A generalization of Einstein's theory of gravity is developed. It has a purely geometric foundation, including in addition to a metric an independent volume. Although related to scalar-tensor theories à la Jordan-Brans-Dicke and to string cosmology, it is nevertheless different: the basic field equations are in the form of Ricci flow equations, generalized to include matter. Einstein's theory is included as the limiting case of no flow. The volume scalar has two interpretations: geometrically, it is responsible for volume-preservation and physically it obeys a mass-zero real scalar wave equation. This is also the main difference with Jordan-Brans-Dicke theories, where the scalar couples to the trace of the energy-momentum tensor. As a consequence, in general the energy-momentum tensor is not anymore "conservative" in the ordinary sense of ∇·T = 0, and point particles do not move anymore on geodesics, having a Newton-Nordström potential. But from Noether's fundamental relation conservation still follows from symmetry. Nevertheless, in "ordinary" solar-system and astrophysical settings, the corrections are negligible: the theory cannot be distinguished from Einstein's and therefore passes all the standard tests. However in a cosmological setting, deviations from standard Einstein gravity are to be expected. This will be the subject of a forthcoming paper.
We will proceed as follows: after this brief introduction, in section 2 the motivations for this kind of extension are discussed. In section 3 a short introduction to volumetric manifolds is given, emphasizing the notion of volume-preservation in section 4. Before establishing the definitive field equations of Ricci flow gravity in section 6, several other choices are discussed in section 5, with emphasis on the main differences with respect to Einstein's and in particular to Jordan-Brans-Dicke theories. Section 7 refers to the physical interpretation of the volume scalar in Ricci flow gravity. Finally, in section 8 the viability of the theory with respect to the standard tests is discussed. The conclusion in section 9 ends this paper.
2 Motivation and Inputs
The present work is principally motivated by the conviction that the notion of "volume" has an existence independent from any metric – in fact, it must be considered to be a pre-metric concept. Curiously, such an independent volume had not been taken into consideration in physical theories until relatively recently. Even in differential geometry it is almost ignored. Perhaps the reason for this neglect is that in most circumstances there is a canonical volume element, based on other geometric structures considered to be more basic. For example, in Riemannian geometry the volume element density is defined in terms of the metric. In particular, the important operation of Hodge dual for differential forms is conventionally based on such a Riemannian volume element.
However, from the gravitational sector of the low-energy limit of string theory (i.e., compactification to dimension n = 4) there comes the suggestive hint (cmp. Garfinkle, Horowitz and Strominger ) that when both a dilaton scalar and a two-form are present, the dilaton scalar enters the expression for the volume element density when defining a "natural" Hodge dual operator.1 This was taken as the starting point to develop a theory of geometric dilaton gravity (Graf ). Although the particular coupling does not exactly correspond to the coupling suggested by string theory, wormholelike solutions were obtained.
Recently a breakthrough on Ricci flow methods was achieved by Perelman [5–7], developing the decisive tools to solve the famous Poincaré conjecture on the topological characterization of the three-sphere. Based on 3-d (compact and positive-definite) Riemann spaces, smoothly deformed by a Ricci flow (RF), the "basic" RF equation was originally defined by
∂ t g ik = -2 R ik ,
where R ik is the Ricci-tensor corresponding to a "time-dependent" three-metric g ik . Also a special class of diffeomorphisms was considered, with vector which is essentially the gradient of a scalar φ in the sense of vi= gik∂ k φ. The so generalized RF equation then becomes
∂ t g ik = -2 (R ik + ∇ i ∇ k φ).
Although such equations have already been studied since the early eighties starting with the seminal works of Hamilton  and DeTurck , an essential insight of Perelman was to recognize that the r.h.s. of this equation2 can be expressed as the gradient of an appropriate functional. This functional involves a "measure" given in terms of an independent volume element. The gradient property allows to apply a series of standard analytical tools. And the introduction of the measure gives an extra flexibility, analogous to a choice of gauge.3
Since the works of Hamilton and DeTurck just mentioned, geometric flows have been applied to a variety of geometric, topological and analytical problems.4 Flow-like equations are also not unfamiliar to physicists, the earliest and most well-known being the renormalization group equations in quantum field theory (for an introduction, see Mitter ), and also the nonlinear σ-model (e.g., Lott , Oliynyk, Suneeta and Woolgar , Tseytlin ). After Ellis  called attention to the cosmological "fitting problem", the usefulness of the Ricci flow to deal with volume-averaged inhomogeneities was immediately recognized and continues to be an active area of research (e.g., Carfora and Marzuoli , Buchert and Carfora , and the recent review by Buchert ). An overview of flow techniques in physics is given in Bakas .
Another motivation came however with the insight, that the basic equations derived from the low-energy limit of string theory can be put into a form suprisingly similar to Ricci flow equations when besides the metric only a dilaton scalar is kept. The main formal difference is the number of dimensions and the signature of the corresponding Riemann spaces: whereas the "classic" RF equations refer to a parameter-dependent truly Riemannian three space evolved by an extrinsic "time" parameter, the reformulated string theory equations refer to a four-dimensional Lorentzian spacetime, which is evolved along the directions of an intrinsic vector field.
Neither the "classic" RF approach nor string theory suggest any hints about the coupling of geometry to external matter fields. Therefore we will spend some time to prepare the field in order to include other external matter. As not only geodesy of the motion of "test particles" will in general be violated, but also "conservation" (in the sense of ∇·T = 0), we will be especially careful to lay a coherent and stringent foundation. The Noether identities will be our main guide. As result we will get Ricci flow gravity (RFG).
3 Volumetrical Manifolds
In Graf  we already introduced the notion of a volume manifold and its specialization for the case a nondegenerate metric exists. Let us briefly recapitulate the main notions. First, we introduced the fundamental concept of a volume structure, which has to be considered as independent from any metric. This is just a non-negative n-form density ω, and makes the manifold a volume manifold. Secondly, we will need of course a metric structure. However, it does not need to be compatible with the volume structure. This difference is encoded by means of the volume scalar φ by ω= ωe-φ,5 where ω := |det g|1/2dx1 ∧ dx2 ∧ ... ∧ dxnis the usual metrical (i.e., Riemannian) volume element density. Furthermore, for the metric derivative along a vector X we have ∇ X ω= -(X·∂φ)ω(the dot denoting a contraction) as a measure of incompatibility.
Such a manifold, endowed both with an independent volume and a metric structure, we will denote by volumetrical manifold. Whereas the manifold is considered to be smooth, both metric and volume element density are allowed to diverge or to be degenerate, when they are not locally smooth.
Already in a volume manifold the Gauss theorem for a vector ξ can be expressed very compactly in terms of differential form densities as
where bΩ is a two-sided hypersurface bounding the n-dimensional region Ω. The scalar factor div ξ in the relation d(ξ·ω) = (div ξ) ωis also better known under the name of divergence of the vector ξ. Evidently the div-operator only depends on the particular choice of ωand not on any metric.
4 Volume Preserving Lie Flows
In a differentiable manifold, the thing coming closest to an autonomous first order differential equation for a "vector" x(t),
is the notion of a Lie equation
£ ξ F = G,
where F and G are geometric objects (e.g., tensors), ξ is some vector field and £ ξ F denotes the Lie derivative of F along ξ. In the simplest case ξ and G are considered as given and F to be determined. However in the applications we have in mind, all elements of the equation will be dynamically determined, G depending nonlinearly on F and its partial derivatives, and even ξ will become dynamical.
In the theory of ordinary differential equations, such systems of first-order equations which guarantee the long-term existence both to the future and the past, are also called flows and can be characterized by the one-parameter Abelian group property of their solutions. As well-known, Lie operators share exactly the same one-parameter Abelian group property (at least locally) by means of the exponential map. We can therefore speak of a Lie flow.
In a volume manifold, a Lie flow with vector ξ is called volume-preserving (or VP)6 if
£ ξ ω= 0, or equivalently, div ξ = 0.
In the rest of this paper we will try to make plausible a particular scalar-tensor extension of Einstein gravity in terms of a volume-preserving Ricci Lie flow in a volumetric manifold.
5 A choice of Scalar-Tensor Field Equations
Assuming that the total Lagrangian (or at least the field equations) can be uniquely split into a pure geometrical part and the "physical" part, we can already draw important conclusions about both the algebraic and the differential properties of the "physical" energy-momentum tensor just from examining the purely geometrical part. Note that whereas Perelman's analysis is "metric-centered", with an auxiliary scalar, in the following physical applications this scalar will play a role at the same conceptual level as the metric. Therefore the "classical" Lagrangian approach is appropriate.
Let us start with the "geometrical" Lagrangian living on a volumetric manifold M,
where ω:= e-φω and (∇φ)2 := gij∂ i φ∂ j φ and λ is a constant parameter. Despite its simple form it not only includes the one used initially by Perelman and in the low-energy limit of string theory (when ignoring the axion and the other moduli fields), but which also is essentially the Jordan-Brans-Dicke Lagrangian.
Defining the volume factor Φ := e-φ, we then have as variational derivatives (up to volume element, g-dualizations of P and a common sign -1)
where G denotes the Einstein tensor G ik := R ik – Rg ik and Δ := ∇2 the d'Alembertian. For the above Lagrangian the Noether identity can be written compactly as
with some tensor to be determined by it. More conventionally,
As this identity must hold for any smooth vector ξ, we get separately
Note that from equation (11) follows conservation in the proper sense, if is a simultaneous Killing vector both of the metric and of the scalar, even if Q ≠ 0.
The following algebraico-differential relations evidently hold:
symmetry: , and
balance: ∇ i (Φ) = Q∂ k Φ.
Up to this point we only made use of identities, but not of any field equations. In particular if we equate P and Q to their corresponding physical quantities, these algebraico-differential relations will be "impressed" on them. In fact, it is not even necessary that they follow from a Lagrangian.
But note that there is a dependency not only on λ but also on the number n of dimensions, making n = 4 and n = 3 (for λ = 1) somewhat special. This is most evident in the relation
where P is the "trace" P := gikP ik . Assuming from now on n = 4, this simplifies to
P + Q = (3 - 2λ)Φ-1ΔΦ.
Let us define the geometrical energy-momentum tensor P ik := 1/2 and more closely examine the corresponding balance relation
The following cases can be distinguished, when equating the geometrical quantities P and Q to their "physical" counterparts P m and Q m :
Pure Einstein, φ = 0
Conformally Einstein, λ = 3/2
"Conservative", Q = 0
VP flow, ΔΦ = 0
Each of these choices will now be discussed individually.
5.1 Pure Einstein
This is just the "compatibility mode", or "Einstein-limit" φ → 0 (if it exists). It is thus a volumetric theory only in the trivial sense of φ = 0.
5.2 Conformally Einstein
The system of equations is underdetermined. This becomes evident by going to the Einstein-frame by means of the conformal transformation , where the scalar field drops out completely.
5.3 "Conservative": Jordan-Brans-Dicke
From Q = 0 there follows "conservation" in the usual sense of ∇ i (Φ) = 0. This assumes the particular relation PΦ = (3 - 2λ) ΔΦ, which could either be postulated or obtained by a specially tailored Lagrangian. In order to have a more familiar looking equation, P ik could be equated to the physical quantity T ik over P ik = Φ-1T ik , so that in fact "conservation" in the sense of ∇ i = 0 would result. This kind of "conservation" was considered as absolutely essential in the closely related scalar-tensor theories of Jordan and Brans-Dicke (in short, JBD theories).7 In fact, this can be achieved as follows: their scalars φ (resp. κ) must be identified with Φ, and we must identify their coupling parameters -ω (resp. ζ) with λ. Moreover λ ≠ 3/2 has to be assumed, otherwise the conformally Einstein theory would result. Then 1/φ (resp. κ) is interpreted as the (variable) gravitational constant. Both for Jordan's and Brans-Dicke's material energy-momentum tensor it is supposed that it does not depend on the scalar φ. However in Jordan's theory it is the product κ2 which is "conserved".8 In particular, for a "dust model" geodesy of (resp. κ2) and conservation of ρ (resp. κ2 ρ) still follow when staying in the original conformal ("Jordan") frame.
5.4 Volume-preserving Flow
When not a conformally Einstein coupling, from equation (14) the condition P + Q = 0 is equivalent to ΔΦ = 0, which in turn is equivalent to volume-preservation Â£ ξ ω = 0. This translates to the scalar condition
P m + Q m = 0
for the corresponding "material" quantities. Let us call such a coupling to matter a volume-preserving material coupling (VPMC), and assume it to hold troughout this section. Then
Evidently, if the trace P m of the energy-momentum tensor vanishes the VPMC is satisfied if we set Q m = 0. Then the standard "conservation" continues to hold. This is the case e.g. for the Maxwell field.
As an important example where P ≠ 0, let us take the ideal fluid model, where the material 6 energy-momentum tensor is given by := Tik= ρuiuk+ p∏ik, and is the projector orthogonal to the trajectory with (normalized) tangent . Its trace is T = 3p - ρ. To satisfy the VPMC, we must set Q m = -T. Specializing to pure dust we get ∇ i (Φρuiuk) = ρgik∂ i Φ. Splitting into tangential and orthogonal parts, we then get the separate equations
Due to the nonvanishing of the r.h.s. of these equations, both "conservation of matter" and geodesy for "test particles" are broken unless φ = const. And due to the particular form of the equation of motion 18b (i.e., being proportional to a gradient) we have in fact got a Newton-Nordström-term.9
Concerning the divergence expression (18a), it can nevertheless be rewritten as a conservation law, ∇ i (Φ1/2ρui) = 0. Therefore, for such a theory with volume-preserving flow, both the equation of motion as well as the "conservation of dust matter" are not anymore the well-known standard expressions from Einstein or Jordan-Brans-Dicke theory. It can be expected that this will have profound consequences in a cosmological setting.
We will continue the discussion of volume-preserving theories in section 6, where we further specialize to the coupling parameter λ = 1.
5.5 Fully dynamic Scalar Field
Here the scalar φ is dynamically determined by a set of field equations obtained via a suitable Lagrangian, and no case of the previously discussed ones fits. This would normally be the "standard" procedere in physics, where not only the Lagrangian is set up as a linear combination of individual Lagrangians, each one describing a different matter model, but in addition possibly introducing some extra "potential terms" containing φ and ∂φ, or even let λ depend on φ. However this will in general prevent a simple geometrical interpretation in terms of a flow, and in particular will lack the crucial VP property. For example, our geometric dilaton gravity (Graf ) belongs to this more general class.
6 Ricci Flow Gravity
The class of volume-preserving volumetric theories can be further refined by requiring the particular value λ = 1 of the coupling, as is common in the low-energy limit of string theory. With this particular value the field equations can be rearranged into an explicit flow-like form and we get the Ricci flow gravity equations (RFG equations)
describing Ricci flow gravity (RF gravity). Here the flow vector is defined in terms of the volume scalar as = -g-1∂φ, 10 and := 8π(T ik – Tg ik ).11 In contrast to the JBD equations, they have a much simpler structure and an immediate geometric character. Through their particular flow-like form, they exhibit a strong dynamical touch: broadly speaking, the rate of change of the metric is driven by the difference of the geometrical and the physical energy momentum tensors. Evidently, when the flow vector can be ignored (e.g., when it vanishes) equations equivalent to Einstein's are obtained. In this sense Einstein gravity is a special case of Ricci flow gravity.
The RFG vacuum equations are equivalent to JBD's vacuum equations with ω = -1. More remarkably is the fact that they are also equivalent to the equations following from the low-energy limit of string theory when the standard dilaton coupling with λ = 1 is chosen and besides the metric only the dilaton scalar is kept. And of course there is a strong resemblance to Perelman's Ricci flow equations which can be made more evident as follows. Consider tentatively on M := M3 × T the vector := -(∂ t + ) and the metric g ik with line element ds2 = dσ2 – dt2. Then the generalized RF equation (2) can be written as , which differs (in content, but not in form) only on the r.h.s. from the corresponding RFG vacuum equation.12
This coincidence of seemingly different approaches could signal a deeper raison d'être.
7 On the physical interpretation of the scalar Φ
The scalar Φ was here interpreted geometrically in the context of a volumetrical manifold as the volume factor. In the theories of JBD the corresponding scalar is essentially interpreted as "gravitational constant" κ – more precisely κ = Φ in Jordan's theory, whereas κ = 1/Φ in Brans-Dicke theory. However this physical interpretation cannot be upheld anymore in a volume-preserving theory like RFG where κ is constant.
Due to the fact that the volume factor Φ of a volume-preserving theory obeys the d'Alembertian wave equation ΔΦ = 0 it must therefore be interpreted as a massless real scalar field. By the tenets of relativistic quantum mechanics this corresponds to a totally uncharged massless bosonic particle.13 The volume-preservation will also be instrumental to guarantee an almost perfect compliance with the standard solar system tests of gravity. This is a fair return for the price we had to pay for giving up the geodesy of "test particles".
8 Compliance of Ricci flow gravity with the Standard Tests
For the standard solar system tests the corresponding generalization of the Schwarzschild metric is needed. The general asymptotically flat and static spherically symmetric vacuum RFG solution with Φ → 1 for r → ∞ can be written as
ds2 = -Yγ-σdt2 + Y-γ-σ(dr2 + Z2dΩ2),
with γ2 + σ2 = 1, and it is assumed that r ≥ r+ > r- ≥ 0.14 Being for r > r+ a vacuum RFG solution, it is also the corresponding general JBD vacuum solution in the "Jordan"-frame. But whereas in JBD gravity the source of the volume factor Φ for a mass point has to be a certain nonzero distribution supported by r = r+, in RF gravity due to ΔΦ = 0 it must be sourceless. This can be shown to hold even for a compactly supported smooth static spherically symmetric energy-momentum tensor as source, if both metric and volume factor are smooth and the manifold is simply-connected. Therefore for RFG σ = 0, whereas for JBD σ = 1/2 (3 + 2ω)-1.
This can also be expressed more conveniently in terms of the "dilaton charge" D, which in the context of the low-energy limit of string theory is defined for a static solution with Killing vector η (normalized to η2 = -1 at infinity) as
where the integral is taken over a closed and externally orientable 2-sphere at spatial infinity.15 For RF gravity the two-form density χ := η·ξ·ωis even closed, dχ = 0, for any stationary solution with Killing vector η so that the above integral only depends on the homology class of the closed externally orientable 2-sphere. In particular it vanishes if this 2-sphere bounds. With the flow vector = σ (r+ – r-)Yγ+σZ-2 ∂ r for the above solution this results in D = σ (r+ - r-). For vanishing dilaton charge the Schwarzschild solution is evidently reobtained after substituting r by r + r-, setting m = (r+ - r-)/2 and assuming m > 0. Thus for the standard solar-system tests the flow vector vanishes and we have full compatibility with Einstein gravity, which passes these tests with ever increasing accuracy (cf. Will ).16
Of course where the flow vector does not vanish, Ricci flow gravity and Einstein gravity will lead to different answers. Using heuristically the term "charge" as introduced above (possibly without stationarity) we note that differently from the "mass charge" m, the "dilaton charge" D can have any sign.17 Therefore the contributions to the total charge of a collection of "charged regions" can still sum up to zero, so as to make the Newton-Nordström-terms of the equation of motion insignificant for sufficiently big distances.18 This should be considered to be in fact the case for "ordinary matter" building up planets, stars and perhaps, galaxies. Significant differences are however to be expected in a cosmological setting, where the "big bang" will affect the behaviour of the volume scalar φ.
Although for the "compliant mode" φ = const evidently it makes no difference if the metric is interpreted in the geometric frame or in the Einstein frame, this is not so in the general case where even the equations of motion for a point particle are modified. We have to chose the particular conformal frame, where the field equations find their "most natural expression". This is the geometrical frame with an independent volume element density.
Motivated by the neglect of the notion of an independent volume and led by the appeal of Perelman's approach to solve the Poincaré conjecture, as well as by the equations following from the low-energy limit of string theory, we developed the equations of Ricci flow gravity as a natural extension of Einstein gravity. The main differences with regard to other scalar-tensor theories were worked out in the framework of volumetric manifolds. The volume-preservation of the flow turned out to be of decisive importance for the theory and allowed it to essentially agree with Einstein's under non-cosmological settings and not too small distances in the case of vanishing "total dilatonic charge".
2the expression in parenthesis in the r.h.s. of equation (2) has its own measure-theoretic meaning and is also known under the name "Bakry-Émery" tensor (cmp. Lott )
3in particular he envisages volume-preserving flows and certain diffeomorphic images thereof
5the factor -1 of φ is purely conventional – here we follow Perelman , in contrast to string theory, where -2 is preferred
6this is a local concept in contrast to the much weaker global definition of Huisken 
8this is suggested by his interpretation of the Kaluza-Klein decomposition
9recall that around 1912–13 Nordström developed a precursor relativistic gravitational theory, where the gravitational potential φ obeys a Minkowskian potential equation, Δφ = 0. This was shown in 1914 by Einstein and Fokker to admit a conformally Minkowskian formulation
10the arbitrary minus-sign is taken in view of cosmological applications
11we use troughout the sign- and units conventions of Misner, Thorne and Wheeler 
13except for a ''dilaton charge''; see the discussion in next section
14in the ''degenerate'' case r+ = r- the metric is locally flat and the volume factor constant
15cmp. Garfinkle, Horowitz and Strominger 
16to compare, for JBD gravity to pass the current tests |ω| > 4 × 104 must be assumed
17this allows the dilaton scalar to act ''repulsively'', as shown in Graf 
18e.g., for a ''multipole charge'' when the distance is much bigger than the individual ''charges''
Garfinkle D, Horowitz GT, Strominger A: Phys Rev D. 1991, 43: 3140-10.1103/PhysRevD.43.3140. [Erratum: Phys Rev D 1992, 45:3888].
Sundrum A: arXiv:hep-th/0312212
Ortín T: Gravity and Strings. 2004, Cambridge: Cambridge University Press
Graf W: Phys Rev D. 2003, 67: 024002-10.1103/PhysRevD.67.024002. [ArXiv:gr-qc/0209002]
Perelman G: arXiv:math.DG/0211159
Perelman G: arXiv:math.DG/0303109
Perelman G: arXiv:math.DG/0307245
Hamilton RS: J Diff Geom. 1982, 17: 255-
DeTurck DM: J Diff Geom. 1983, 18: 157-
Lott J: arXiv:math.DG/0211065
Chow B, Knopf D: The Ricci Flow: An Introduction. 2004, Providence: American Mathematical Society
Topping PM: Lectures on the Ricci flow. 2006, Cambridge: Cambridge University Press
Morgan J, Tian G: arXiv:math.DG/0607607
Mitter P: arXiv:math-ph/0505008
Lott J: Commun Math Phys. 1986, 107: 165-10.1007/BF01206956.
Oliynyk T, Suneeta V, Woolgar E: Nucl Phys B. 2006, 739: 441-10.1016/j.nuclphysb.2006.01.036. [ArXiv:hep-th/0510239]
Tseytlin AA: Phys Rev D. 2007, 75: 064024-10.1103/PhysRevD.75.064024. [ArXiv:hep-th/0612296]
Ellis GFR: Proceedings of the Tenth International Conference on General Relativity and Gravitation. Edited by: Bertotti B, de Felice F, Pasolini A. 1984, Dordrecht: Reidel, 215-
Carfora M, Marzuoli A: Phys Rev Lett. 1984, 53: 2245-10.1103/PhysRevLett.53.2445.
Buchert T, Carfora M: Phys Rev Lett. 2002, 90: 031101-10.1103/PhysRevLett.90.031101.
Buchert T: arXiv:0707.2153 [gr-qc]
Bakas I: arXiv:hep-th/0511057
Brans C: arXiv:gr-qc/0506063
Fujii Y, Maeda KI: The Scalar-Tensor Theory of Gravitation. 2003, Cambridge: Cambridge University Press
Faraoni V: Cosmology in Scalar-Tensor Theory. 2004, Dordrecht: Kluwer
Huisken G: J reine angew Math. 1987, 382: 35-
Jordan P: Schwerkraft und Weltall. 1955, Braunschweig: Vieweg
Weinberg S: Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity. 1972, New York: Wiley
Misner CW, Thorne KS, Wheeler JA: Gravitation. 1973, San Francisco: Freeman
Bleecker DD: J Geom Phys. 1990, 7: 363-10.1016/0393-0440(90)90028-2.
List B: PhD thesis. 2005, Freie Universität Berlin
Will CM: [http://www.livingreviews.org/lrr-2006-3]
I thank the members of the relativity group of the University of Vienna for discussions and useful suggestions, and in particular Peter C. Aichelburg. Partial financial support from the Fundación Federico is also acknowledged.