Next Article in Journal
Quantum Gravity Phenomenology from the Thermodynamics of Spacetime
Next Article in Special Issue
Optical Features of AdS Black Holes in the Novel 4D Einstein-Gauss-Bonnet Gravity Coupled to Nonlinear Electrodynamics
Previous Article in Journal
Modern Approaches to Non-Perturbative QCD and Other Confining Gauge Theories
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Note on Singularity Avoidance in Fourth-Order Gravity

DIME, Sez. Metodi e Modelli Matematici, Università di Genova, Via all’Opera Pia 15, 16145 Genova, Italy
Universe 2022, 8(1), 51; https://doi.org/10.3390/universe8010051
Submission received: 5 December 2021 / Revised: 6 January 2022 / Accepted: 11 January 2022 / Published: 13 January 2022
(This article belongs to the Special Issue Higher-Derivative Theories of Gravity)

Abstract

:
We consider the fourth-order differential theory of gravitation to treat the problem of singularity avoidance: studying the short-distance behaviour in the case of black-holes and the big-bang we are going to see a way to attack the issue from a general perspective.

1. Introduction

With the detection of gravitational waves, each original experimental prediction of Einstein gravity has now been settled. In fact, if dark matter is indeed a form of matter, and not a gravitational effect, there are no observational open issues left in modern gravity. Nevertheless, from a purely theoretical perspective, there are yet two problems that need fixing, and which are, more or less, connected: one is solving the nature of singularity formation, which seems to be an occurrence that is unavoidable, in light of the Hawking-Penrose theorem; the other, to have gravity made into a renormalizable theory, hence to fit, with the standard model of particles, into one single framework.
Because renormalizability needs all fields with a kinetic term of mass dimension 4 and because for gravitation the mass dimension is assigned to be equal to 2 for curvature terms, Einsteinian gravitation cannot meet the necessary requirements. However, the solution is simple: raise to 2 the number of curvatures appearing in the Lagrangian.
This strategy has led to a variety of extensions of Einstein gravity, of which the f ( R ) types are just the most famous to have arisen in recent times (for a general overview, we refer the reader to [1,2,3,4,5] and references therein).
This type of gravitational extension may meet also the criteria to solve the problem of singularity formation [6].
However, by not stopping at two curvatures, they tend to include in the Lagrangian terms that have too high a differential order for renormalizability. Renormalizability is instead granted for theories displaying a conformal invariance [7,8,9]. Yet, when torsion is included, it is not difficult to see that this theory is not continuous, in the sense that conformal gravity with torsion taken in its torsionless limit does not give the conformal gravity that we would have without having torsion in the first place [10].
This problem is in fact very general, affecting virtually everyone of the higher-order theories of gravitation [11].
As a matter of fact, a thorough application of the principle requiring continuity for the torsionless limit of any torsion gravity imposes severe restrictions to the possible forms that the gravitational Lagrangian can have [12].
Remarkably, by following this principle the Lagrangian found in [12] is the one that grants at most 2 curvatures and the renormalizable kinetic term for torsion [13,14].
So a question naturally arises, asking whether this Lagrangian can also provide some solution to the singularity problem. The answer is partly positive, because a renormalizable torsion theory, where torsion is an axial-vector field coupled to the spin of spinors, ensures the conditions for which the Hawking-Penrose inequality is violated and singularities no longer unavoidable [15] (and as a matter of fact, even the Higgs field can do the same [16]), so long as black holes are considered. Singularities that involve the big bang instead are more delicate because symmetry arguments may be used to demonstrate that torsion does not impact the energy density for the gravitational field equations [17] (for the Higgs the argument is that at the energy scales of the big bang no spontaneous symmetry breaking has occurred yet). Thus the question, could we have a way to avoid singularity formation, or at least its inevitability, even for early cosmological scenarios?

2. Fourth-Order Differential Torsion-Gravity

We start by recalling the results of [12]. In this work it was shown we get continuity, thus the torsionless limit of torsion gravity gives pure gravity, only if the Lagrangian is restricted to a very special type: this Lagrangian is
L = 1 4 ( α W ν ν W α ) ( α W ν ν W α ) Y R α μ R α μ Z R 2 + 1 2 M 2 W ν W ν R + + i ψ ¯ γ μ μ ψ X ψ ¯ γ μ π ψ W μ m ψ ¯ ψ
where W α is the axial-vector field known as torsion, R α μ and R are the Ricci tensor and scalar, ψ and ψ ¯ the pair of adjoint spinor fields. The metric is given by g μ ν and the connection Λ μ ν ρ is used to compute the covariant derivatives, with γ μ being the Clifford matrices ( π is normally denoted by a gamma with an index 5, but because in the space-time this index has no meaning we will employ the definition with no index). Finally, m and M are the mass of the spinor and torsion with X being the torsion-spinor coupling constant. Notice that the case Y = Z = 0 reduces the Lagrangian to its least-order derivative form, but in general this Lagrangian has 2 curvatures and so it has fourth-order differential character in the metric. Remark also that the term R should be multiplied by a constant with mass dimension 2 and that is the Newton constant, although we have normalized it to unity and therefore it does not explicitly show within the Lagrangian. Because the metric is dimensionless, every curvature counts for a mass dimension 2 so that the general kinetic term is mass dimension 4 with strong consequences for the requirement of renormalizability as discussed in references [13,14].
Varying the torsion axial-vector, the metric tensor and the spinor field, we get respectively the field equations
ρ ( W ) ρ μ + M 2 W μ = X ψ ¯ γ μ π ψ
for torsion with
Y 2 R μ ν + 1 2 ( 4 Z + Y ) 2 R g μ ν ( 2 Z + Y ) μ ν R + + 2 Y R μ ρ ν σ R ρ σ 1 2 Y R α ρ R α ρ g μ ν + 2 Z R R μ ν 1 2 Z R 2 g μ ν + R μ ν 1 2 g μ ν R = = 1 2 [ 1 4 ( W ) 2 g μ ν ( W ) ν α ( W ) μ μ α + M 2 ( W μ W ν 1 2 W 2 g μ ν ) + + i 4 ( ψ ¯ γ μ ν ψ ν ψ ¯ γ μ ψ + ψ ¯ γ ν μ ψ μ ψ ¯ γ ν ψ ) 1 2 X ( W ν ψ ¯ γ μ π ψ + W μ ψ ¯ γ ν π ψ ) ]
for gravity and
i γ μ μ ψ X W σ γ σ π ψ m ψ = 0
for the spinor field. By taking the divergence of (2) and contracting (3) we get the constraints
M 2 μ W μ = 2 X m i ψ ¯ π ψ
and
4 ( 3 Z + Y ) 2 R 2 R = M 2 W 2 + m ψ ¯ ψ
which will be important in the following. Again, the case Y = Z = 0 reduces all to the least-order derivative form we would have in the usual Einsteinian second-order gravity, although in general the kinetic term ensures fourth-order character to gravitation. Because each curvature has the mass dimension equal to 2, single-curvature kinetic terms of second-order gravity scale as l 2 and double-curvature kinetic terms of fourth-order gravity scale as l 4 with the consequence that Einstein gravity is more relevant in the case of large l while its fourth-order extension is dominant for small l in general. Consequently, the usual theory of gravitation is recovered in the infrared whereas some new physics is expected in the ultraviolet as we will see next.

High Energy and Averaged Spins

In the next pages we will be focusing on the problems usually met in the ultraviolet, that is at short distances, and hence high energies. In addition, since the situations we want to study (black holes and big bang) involve very large number of spinors with randomly distributed spin, an average on spin will also be assumed. When the high energy density condition is implemented, whenever short distances are considered, terms with high mass dimension become dominant as compared to all else, and thus such a condition of ultraviolet regimes can be implemented by requiring that only the mass dimension 4 terms are kept in all field equations. After doing this, we lose the mass and all single-curvature terms, remaining with
ρ ( W ) ρ μ = X ψ ¯ γ μ π ψ
for torsion and
Y 2 R μ ν + 1 2 ( 4 Z + Y ) 2 R g μ ν ( 2 Z + Y ) μ ν R + 2 Y R μ ρ ν σ R ρ σ 1 2 Y R α ρ R α ρ g μ ν + + 2 Z R R μ ν 1 2 Z R 2 g μ ν = 1 2 [ 1 4 ( W ) 2 g μ ν ( W ) ν α ( W ) μ μ α + + i 4 ( ψ ¯ γ μ ν ψ ν ψ ¯ γ μ ψ + ψ ¯ γ ν μ ψ μ ψ ¯ γ ν ψ ) 1 2 X ( W ν ψ ¯ γ μ π ψ + W μ ψ ¯ γ ν π ψ ) ]
for gravity together with
i γ μ ( μ ψ + i X W μ π ψ ) = 0
for the spinor field. The constraints are
i ψ ¯ π ψ = 0
and
2 R = 0
as easy to see. Notice that condition i ψ ¯ π ψ = 0 is equivalent to the vanishing of the Yvon-Takabayashi angle, and this is expected in ultra-relativistic limits such as those obtained in massless conditions [18]. Condition 2 R = 0 will prove to be important later on. Averaging spins has the same effect of requiring that all instances of the spin axial-vector ψ ¯ γ μ π ψ be neglected, which leaves
ρ ( W ) ρ μ = 0
for torsion with
Y 2 R μ ν ( 2 Z + Y ) μ ν R + 2 Y R μ ρ ν σ R ρ σ 1 2 Y R α ρ R α ρ g μ ν + + 2 Z R R μ ν 1 2 Z R 2 g μ ν = 1 2 [ 1 4 ( W ) 2 g μ ν ( W ) ν α ( W ) μ μ α ]
for gravity and no more occurrences of spinors. This fact is understood by recalling that with no Yvon-Takabayashi angle nor spin axial-vector, we have that P ν = m u ν [18], with the spinor contribution in the energy turning into a mass dimension 3 term that, again in high energy density conditions, will become negligible. These last equations, the last of which now implying the subsidiary condition
2 R = 0
in general, are the basis for the study of those symmetric situations that are at the foundation of modern cosmology, and that is the big bang and black holes. In detail, in the following we will consider two situations of symmetry given when the space-time is isotropic, a first of pure isotropy and a second where isotropy is accompanied by homogeneity. The first case, corresponding to stationary spherical symmetry, is the best suited to represent black holes (unless we want to consider the rotating case where axial symmetry should be used instead), and the second case, corresponding to maximal symmetry, is best suited to represent the scenarios of cosmological evolution.

3. Isotropic Spaces

The treatment of singularity formation revolves around ways to evaluate whether the Hawking-Penrose dominant energy condition R μ ν u μ u ν > 0 is respected. For Einstein gravity the field equations are given in terms of the Ricci tensor R μ ν and so this condition can be examined rather straightforwardly, but in higher-order extensions the field equations have an altogether different form and therefore this direct study cannot generally be done. One alternative is that of focusing on a particular situation, described by a specific symmetry, finding exact solutions, and see that in such a case the dominant energy condition comes to be violated. This method is not general, as it is fully and strictly tied to the particular situation one wants to examine. Hence, one would not be able to prove that the dominant energy condition is always violated, but solely that it can be violated. Nevertheless, for now this is the most we can hope for. So, in the following, we shall focus on two particular (although well known) situations.

3.1. Black Holes

Let us begin the study by testing our equations against the case of black holes. This case will be considered not rotating, so to exploit spherical symmetry, and therefore, in polar coordinates, we have that the metric is given by
g t t = A g r r = B g θ θ = r 2 g φ φ = r 2 | sin θ | 2
with A ( r ) and B ( r ) in general. The connection is
Λ t r t = A 2 A Λ t t r = A 2 B Λ r r r = B 2 B Λ θ θ r = r B Λ φ φ r = r B | sin θ | 2 Λ φ φ θ = cos θ sin θ Λ θ r θ = Λ φ r φ = 1 r Λ φ θ φ = cot θ
from which covariant derivatives and curvatures can be computed. In particular, the curvatures are given by
R t r t r t r = A 2 A B A 2 4 A 2 B A B 4 A B 2 R t θ t θ t θ = R t φ t φ t φ = A 2 A B r R r θ r θ r θ = R r φ r φ r φ = B 2 B 2 r R θ φ θ φ θ φ = 1 B r 2 1 r 2
with contraction
R t t t = A 2 A B A 2 4 A 2 B A B 4 A B 2 + A A B r R r r r = A 2 A B A 2 4 A 2 B A B 4 A B 2 B B 2 r R θ θ θ = A 2 A B r B 2 B 2 r + 1 B r 2 1 r 2 R φ φ φ = A 2 A B r B 2 B 2 r + 1 B r 2 1 r 2
and contraction
R = A A B A 2 2 A 2 B A B 2 A B 2 + 2 A A B r 2 B B 2 r + 2 B r 2 2 r 2
from which we calculate all terms in the field equations.
Torsion would have only W t ( r ) and W r ( r ) albeit only W t ( r ) can appear in the curl within the field equations.
However, we notice that precisely because of (12) there can be no dynamical solution for torsion, which can then be taken to vanish in the energy density contribution.
As another consequence of the field equations, (13) this time, we can pick A B = 1 without any loss of generality.
With this wisdom, we can proceed to observe that the constraint 2 R = 0 can be solved with some non-singular R only if R is a constant. To see this we have to consider that the Laplacian is generally expressed by
2 R 1 | g | μ | g | g μ ν ν R
where | g | is the determinant of the metric, and so
2 R = 1 r 2 r 2 B 1 R
for our metric and where the prime stands for derivative with respect to the radial coordinate. So 2 R = 0 can be integrated easily as
R = a B r 2
with a a constant. This yields a divergent behaviour for R in the origin unless B behaves as r b with b > 2 although in such a case the metric would degenerate in the origin and thus the solution would still be singular. Hope for a non-singular behaviour can only be found in having a = 0 and hence R constant. Setting R = 12 k we get therefore
12 k = A + 4 A / r + 2 ( A 1 ) / r 2 = [ r 2 ( A 1 ) ] / r 2
which is now easy to solve. In fact
A = 1 + k r 2
is the only metric still non-singular. This metric is exact solution for the field Equation (13) in this particular case.
The above solution (24) in the limit r 0 becomes flat consequently showing no formation of singularity at all.
A point is worth commenting. The solution (24) represents what we expect to be the solution near the center of the matter distribution, that is black holes with spherical symmetry in conditions of high energy density, although this does not mean that (24) is the full solution. The full solution is a metric that has (24) as the limit behaviour at the center of the matter distribution where large energy densities are found. Because (24) is an exact solution of the approximated field Equation (13), is a limiting case of the exact solution to the complete field Equation (3).
It is considerably difficult to obtain such a most general solution. But for the purposes we have here all we need is to assess its ultraviolet behaviour, that is for large energy density near the center of the material distribution.

3.2. Big Bang

Having worked out a rather general way to treat such a singularity problem in the case of black holes, we will move to study the case of the big bang. Now we have both isotropy and homogeneity, and so in polar coordinates
g t t = 1 g r r = A 2 g θ θ = A 2 r 2 g φ φ = A 2 r 2 | sin θ | 2
with A ( t ) in general. The connection is
Λ r r t = A A ˙ Λ θ θ t = A A ˙ r 2 Λ φ φ t = A A ˙ r 2 | sin θ | 2 Λ r t r = Λ θ t θ = Λ φ t φ = A ˙ A Λ θ θ r = r Λ φ φ r = r | sin θ | 2 Λ φ φ θ = cos θ sin θ Λ θ r θ = Λ φ r φ = 1 r Λ φ θ φ = cot θ
from which covariant derivatives and curvatures can be computed. In particular, the curvatures are given by
R t r t r t r = R t θ t θ t θ = R t φ t φ t φ = A ¨ A R r θ r θ r θ = R r φ r φ r φ = R θ φ θ φ θ φ = A ˙ 2 A 2
with contraction
R t t t = 3 A ¨ A R r r r = R θ θ θ = R φ φ φ = 2 A ¨ 2 A + A ˙ 2 A 2
and contraction
R = 6 A ¨ A + A ˙ 2 A 2
from which we calculate all terms in the field equations.
Torsion has only W t ( t ) and therefore no contribution.
With the same reasoning as above we can try to search for the cases in which R is constant. So setting R = 12 k
2 k = A ¨ A + A ˙ 2 A 2 = A 2 ¨ 2 A 2
or equivalently
A 2 ¨ + 4 k A 2 = 0
as a well known ordinary differential equation. With freedom to choose k we may notice that for k = a 2 we obtain
A 2 ¨ 4 a 2 A 2 = 0
as the only field equation admitting non-vanishing solutions and with
A 2 = e 2 a t
as the only expanding solution. So
A = e a t
is the only possible solution that will remain non-singular at all times. Again, field Equation (13) are verified.
The above solution (34) in the limit t 0 becomes flat consequently showing no formation of singularity again, but now the specific type of solution has acquire an additional piece of information. In fact, allowing only solutions that do not degenerate the metric, that is excluding solutions of the type sin ( a t ) or cos ( a t ) , not only we have eliminated the possibility of recurring universes, but we have also fixed the Ricci scalar to be negative. And this will have a very specific effect on the Ricci scalar as well.
In fact the Ricci tensor is R μ ν = 3 a 2 g μ ν and in turn R μ ν u μ u ν < 0 for all curves of tangent vector u α time-like.
This blunt violation of the Hawking-Penrose dominant energy condition is the clearest sign that the singularity formation is not longer a necessary feature of these types of extended theories of the gravitational field in general.
Finally, it may be worth mentioning that the solution (34) also gives rise to an era of inflationary expansion as the one we would expect to have in the early universe.
Notice that in this model such an evolution is obtained without any cosmological constant and this is important since any cosmological constant term has mass dimension zero and so it is negligible at high energy densities.

4. General Consideration

In the previous analysis, we have seen two cases where it was in fact possible to find exact solutions as metrics that were non-singular and asymptotically flat, such that the corresponding space-times violated the dominant energy condition needed for singularity formation. Hence, the formation of singularities was not a necessity in these circumstances. Being restricted to specific situations, the analysis was not general, and thus we could only show the possibility of a solution to the problem, and not its whole applicability. Just the same, the results are intriguing.
Apart from the symmetries of the specific situations at hand, the only other element playing a role was the existence of fourth-order differential gravitational field equations. The found solutions would in fact not be solutions in the least-order derivative Einsteinian gravitation. Nor would they work for higher-order derivative extensions based solely on the Ricci scalar. In fact, in this case the analogous of field Equation (13) would be given by
f ( R ) R μ ν 1 2 f ( R ) g μ ν = 0
and hence
6 a 2 f ( a 2 ) + f ( a 2 ) = 0
which is not generally verified. They can for special situations such as f ( R ) = R 2 but this would merely coincide with the case Y = 0 in the field equations used above.
The method we used is certainly tied to the underlying symmetries of the system, but it is a character of fourth-order theories of gravitation and no other one.

5. Conclusions

In this paper, we have recalled the fourth-order theory of gravitation showing, in high-energy average-spin configuration, that the field equations admit exact solutions of specific symmetries: we have studied cases of isotropy and isotropy and homogeneity. Isotropy, that is spherical symmetry, was best suited to examine the behaviour we might find in black holes, for which we had demonstrated the existence of metrics that were non-singular near the origin of the radial coordinate, whereas isotropy and homogeneity, that is maximal symmetry, was best suited for examining the behaviour one might find at the big bang, for which we have proven the existence of metrics that are non-singular near the origin of time and whose structure provides the Ricci scalar with a negative value imposing the violation of the Hawking-Penrose energy conditions.
Such a procedure pertains only to fourth-order theories of gravitation and to no other theory in general cases.
However, the method is highly dependent on the symmetries of the system. Can this methodology be extended so to include more general situations?

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Manuscript has no associated data in any repository.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Cai, Y.F.; Duplessis, F.; Saridakis, E.N. F(R) nonlinear massive theories of gravity and their cosmological implications. Phys. Rev. D 2014, 90, 064051. [Google Scholar] [CrossRef] [Green Version]
  2. Nashed, G.G.L.; Saridakis, E.N. New rotating black holes in nonlinear Maxwell f(R) gravity. Phys. Rev. D 2020, 102, 124072. [Google Scholar] [CrossRef]
  3. Capozziello, S.; Vignolo, S. The Cauchy problem for f(R)-gravity: An Overview. Int. J. Geom. Meth. Mod. Phys. 2012, 9, 1250006. [Google Scholar] [CrossRef] [Green Version]
  4. Vignolo, S.; Cianci, R.; Carloni, S. On the junction conditions in f(R)-gravity with torsion. Class. Quant. Grav. 2018, 35, 095014. [Google Scholar] [CrossRef] [Green Version]
  5. Vignolo, S. Some Mathematical Aspects of f(R)-Gravity with Torsion: Cauchy Problem and Junction Conditions. Universe 2019, 5, 224. [Google Scholar] [CrossRef] [Green Version]
  6. Cai, Y.F.; Easson, D.A.; Brandenberger, R. Towards a Nonsingular Bouncing Cosmology. JCAP 2012, 8, 20. [Google Scholar] [CrossRef]
  7. Mannheim, P.D. Solution to the ghost problem in fourth order derivative theories. Found. Phys. 2007, 37, 532. [Google Scholar] [CrossRef] [Green Version]
  8. Mannheim, P.D. Comprehensive Solution to the Cosmological Constant, Zero-Point Energy, and Quantum Gravity Problems. Gen. Relativ. Gravit. 2011, 43, 703. [Google Scholar] [CrossRef] [Green Version]
  9. Mannheim, P.D. Making the Case for Conformal Gravity. Found. Phys. 2012, 42, 388. [Google Scholar] [CrossRef] [Green Version]
  10. Fabbri, L. Metric solutions in torsionless gauge for vacuum conformal gravity. J. Math. Phys. 2013, 54, 062501. [Google Scholar] [CrossRef] [Green Version]
  11. Fabbri, L.; Mannheim, P.D. Continuity of the torsionless limit as a selection rule for gravity theories with torsion. Phys. Rev. D 2014, 90, 024042. [Google Scholar] [CrossRef] [Green Version]
  12. Fabbri, L. A discussion on the most general torsion-gravity with electrodynamics for Dirac spinor matter fields. Int. J. Geom. Meth. Mod. Phys. 2015, 12, 1550099. [Google Scholar] [CrossRef] [Green Version]
  13. Stelle, K.S. Classical Gravity with Higher Derivatives. Gen. Relativ. Gravit. 1978, 9, 353. [Google Scholar] [CrossRef]
  14. Stelle, K.S. Renormalization of Higher Derivative Quantum Gravity. Phys. Rev. D 1977, 16, 953. [Google Scholar] [CrossRef]
  15. Fabbri, L. Singularity-free spinors in gravity with propagating torsion. Mod. Phys. Lett. A 2017, 32, 1750221. [Google Scholar] [CrossRef]
  16. Fabbri, L. Black Hole singularity avoidance by the Higgs scalar field. Eur. Phys. J. C 2018, 78, 1028. [Google Scholar] [CrossRef]
  17. Fabbri, L. A simple assessment on inflation. Int. J. Theor. Phys. 2017, 56, 2635. [Google Scholar] [CrossRef] [Green Version]
  18. Fabbri, L. Spinors in Polar Form. Eur. Phys. J. Plus 2021, 136, 354. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fabbri, L. A Note on Singularity Avoidance in Fourth-Order Gravity. Universe 2022, 8, 51. https://doi.org/10.3390/universe8010051

AMA Style

Fabbri L. A Note on Singularity Avoidance in Fourth-Order Gravity. Universe. 2022; 8(1):51. https://doi.org/10.3390/universe8010051

Chicago/Turabian Style

Fabbri, Luca. 2022. "A Note on Singularity Avoidance in Fourth-Order Gravity" Universe 8, no. 1: 51. https://doi.org/10.3390/universe8010051

APA Style

Fabbri, L. (2022). A Note on Singularity Avoidance in Fourth-Order Gravity. Universe, 8(1), 51. https://doi.org/10.3390/universe8010051

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