Next Article in Journal
Torsion at Different Scales: From Materials to the Universe
Next Article in Special Issue
Efficient Computation of Null Affine Parameters
Previous Article in Journal
Tidal Resonance: A Factor Worth Considering in the Orbital Evolution of Heartbeat Stars
Previous Article in Special Issue
Schrödinger Symmetry in Gravitational Mini-Superspaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spherically Symmetric Configurations in Unimodular Gravity

by
Júlio César Fabris
1,2,*,
Mahamadou Hamani Daouda
3 and
Hermano Velten
4
1
Núcleo Cosmo-Ufes & Departamento de Física, Universidade Federal do Espírito Santo, Av. Fernando Ferrari, 514, Goiabeiras, Vitória 29060-900, Brazil
2
Moscow Engineering Physics Institute, National Research Nuclear University MEPhI, Kashirskoe Shosse 31, Moscow 115409, Russia
3
Département de Physique, Université de Niamey, Niamey 10662, BP, Niger
4
Departamento de Física, Universidade Federal de Ouro Preto (UFOP), Campus Universitário Morro do Cruzeiro, Ouro Preto 35402-136, Brazil
*
Author to whom correspondence should be addressed.
Universe 2023, 9(12), 515; https://doi.org/10.3390/universe9120515
Submission received: 6 November 2023 / Revised: 4 December 2023 / Accepted: 12 December 2023 / Published: 14 December 2023
(This article belongs to the Special Issue Universe: Feature Papers 2023—Gravitation)

Abstract

:
Unimodular gravity (UG) is often deemed comparable to General Relativity (GR) in many respects, despite the theory exhibiting invariance under a more limited set of diffeomorphic transformations. The discussion we propose in this work relies on the criteria for establishing the equivalence between these two formulations, specifically exploring UG’s application to static and spherically symmetric configurations with the energy-momentum tensor originating from either a scalar field or an electromagnetic field. We find that the equivalence between UG and GR might be disrupted when scrutinizing the stability of solutions at a perturbative level.

1. Introduction

General Relativity (GR) is the modern theory of gravitational interaction. The gravitational phenomena are considered as the structure of the space-time itself induced dynamically by matter. GR is considered a very successful theory: all local tests confirm the predictions of GR with high precision. At cosmological scales, it leads to the Standard Cosmological Model (SCM) which also addresses consistently all available observations, from scales of galaxies up to the larger structures of the universe. It accounts also successfully for the different phases of the evolution of the universe, including the primordial phases at least to the primordial nucleosynthesis scales. From this point of view, the SCM, based on GR, is an almost perfect model to describe the entire evolution of the universe.
However, seen from a different perspective, the SCM is at least problematic. To account for the observations at the different scales, it demands the introduction of two until now undetected components in the matter/energy content of the universe [1,2]. The dynamics of galaxies and cluster of galaxies, and even the formation of such structures, asks for an additional pressureless component, dubbed dark matter, which manifests only indirectly. Moreover, to explain the present accelerated phase of the universe, the CMB spectrum and to obtain an age of the universe consistent with the age of globular clusters, the SCM asks for another component, with negative pressure, which does not agglomerate, dubbed dark energy.
Dark energy is now frequently associated with vacuum energy as predicted by quantum field theory. However, its observed value seems not consistent with the theoretical predictions by dozens of orders of magnitude [3,4,5]. There are many proposals to cope with this problem. One of them is to replace it with a self-interacting scalar field, called quintessence [6]. However, it must be explained, in the quintessence program, why the vacuum energy must be exactly or, at least, nearly zero. Therefore, the vacuum energy must, somehow, degravitate [7,8,9]. There are many mechanisms to implement such degravitating mechanism, but until now such proposals are, in some sense, in construction. For a general overview of the dark energy problem, see Ref. [10].
One interesting approach to the cosmological constant problem described above is through the unimodular gravity (UG) class of theories [11,12,13,14] where the determinant of the metric, g, is fixed. Originally, g = 1 , but other possibilities can be explored, see next section. UG leads to traceless gravitational equations. The energy-momentum tensor does not conserve necessarily anymore, since UG is not invariant by the full diffeomorphism group, but by a more restricted structure called transverse diffeomorphism [15]: While the general diffeomorphism is characterized by a coordinate transformation x μ x μ + ξ μ , x μ being an arbitrary vector, the transverse diffeomorphism requires that this function satisfies the condition ξ ; μ μ = 0 , as it will be seen later. If the conservation of the energy-momentum tensor is imposed, GR is recovered with a cosmological term that appears as an integration constant. If the conservation of the energy-momentum tensor is not imposed, as we will prove below, a class of dynamical vacuum theories is obtained, implying an interaction of the matter sector with the decaying cosmological term.
The present motivation for the UG deviates from the original one, which is closer to the choice of a coordinate system in GR. From the modern perspective, UG has two main motivations. The first one is related to the cosmological constant problem: if we introduce a cosmological constant term in the matter content, it disappears from the UG gravitational equations. It may reappear if the energy-momentum tensor is supposed to conserve, but the possible connection with the quantum vacuum energy density is lost and the usual cosmological constant problem is alleviated since it remains only to explain the observed value, but without a necessary connection with the vacuum energy. UG may have some interesting properties at a quantum level related to the volume fixing constraint, see Ref. [14] and references therein.
However, there are many discussions, from the original formulation of UG to its present applications, if the theory is equivalent to GR with a constant or dynamical cosmological term. In previous works we have explored the distinction between GR and UG mainly at a perturbative level in the cosmological context, see [16] and references therein (see also [17]). Here we will extend such studies to the static, spherically symmetric configurations. In UG, with the imposition of the conservation of the energy-momentum tensor, the static, spherically symmetric solutions are identical to those of GR but now contain a cosmological constant. The non-conservation of the energy-momentum tensor, on the other hand, can be mapped in the GR structure with a dynamical cosmological term. Indeed, the non-conservation of the energy-momentum tensor is allowed in this context, leading to a new formulation of the UG theory. It is also worth mentioning that such a non-conservation mechanism appears in many other situations. For a review, see Ref. [18]. All these aspects are discussed in the next section. In Section 3 the general equations for a static, spherically symmetric configuration are settled out. Some examples of interacting models, resulting from the non-conservation of the energy-momentum tensor, will be shown in Section 4, both in the presence of an electromagnetic field as well as of a self-interacting scalar field in Section 5. For the latter case, we perform, in Section 6, a perturbative analysis aiming to show how the usual results of GR change in the unimodular context. In particular, the unimodular condition on the determinant of the metric implies vanishing perturbations at the linear level. The results obtained are discussed in Section 7.

2. Field Equations

The Einstein–Hilbert action, in the presence of the cosmological term and the matter Lagrangian,
S = d 4 x g { R 16 π G + 2 Λ + L m } ,
implies in the following field equations:
R μ ν 1 2 g μ ν R = 8 π G T μ ν + g μ ν Λ .
The application of the Bianchi identities leads to the energy-momentum tensor T μ ν conservation:
T μ ν ; μ = 0 .
The conservation laws related to the energy-momentum tensor can be alternatively deduced from the invariance of the Einstein–Hilbert Lagrangian by diffeomorphic transformations [19].
In order to obtain the UG equations, we introduce a constraint in the action via a Lagrange Multiplier χ and an external field ξ [17]:
S = d 4 x { g R χ ( g ξ ) } + d 4 x g L m .
The presence of the external field allows one to use a suitable coordinate system according to the problem under analysis, for example, the usual spherical coordinates or the quasi-global coordinates employed in spherical symmetric space-time.
The final field equations for this case are
R μ ν 1 4 g μ ν R = 8 π G ( T μ ν 1 4 g μ ν T ) ,
R ; ν 4 = 8 π G ( T μ ν ; μ 1 4 T ; ν ) .
The above Equation (6) is obtained by using the Bianchi identities in (5).
As highlighted in Ref. [19], it is important to note that in UG, the conservation of the energy-momentum tensor cannot be derived through the conventional diffeomorphism invariance because the theory only exhibits invariance with respect to a limited set of diffeomorphisms, referred to as transverse diffeomorphisms. The latter implies that the energy-momentum divergence tensor is equal to the gradient of a (undetermined) scalar function:
T ν μ ; μ = Θ ; ν ,
On one hand, it is entirely permissible to set the gradient of Θ to zero. If this is conducted, we recover (2), with Λ appearing as an integration constant. On the other hand, one can also choose,
Θ = R 4 + 2 π G T .
As it will be seen later, if the energy-momentum tensor is not conserved ( Θ 0 ), the final system of equations is underdetermined, and a supplementary ansatz must be introduced to close consistently the system. In practice, this means that the functional form of the relation (8) depends on the ansatz to be adopted. Due to this, from now on, we will identify Θ Λ , in order to keep in contact with the usual notation in the literature for decaying vacuum models. If Λ is constant, as already stressed, we return to the GR equations in the presence of a cosmological constant. But, if Λ is a function of the space-time coordinates, we end up with the following set of equations,
R μ ν 1 4 g μ ν R = 8 π G { T μ ν 1 4 g μ ν T } ,
T ν μ ; μ = Λ ; ν .
To proceed with a dynamical Λ , we must have some guiding physical perspectives. Equations (9) and (10) represent generally the idea of a decaying vacuum. A decaying vacuum may be an effective mechanism to begin with a huge initial value for the vacuum energy density, reaching later a small value as observed today. To compute this decaying process, a deep comprehension of the dynamics of quantum fields in a gravitational field is required. Perhaps this may require a full quantum gravity theory. In the absence of complete quantum gravity theory, we can only try to guess, by using some reasonable assumptions, how the vacuum can decay and that is the usual approach in the literature, for example, the case of running gravitational coupling, for example, Refs. [20,21,22], or scale dependent gravity, see Ref. [23].
This is equivalent (up to the restriction in the diffeomorphic class of transformation) to the GR equations in the presence of a dynamical cosmological term:
R μ ν 1 2 g μ ν R = 8 π G T μ ν + g μ ν Λ ,
T ν μ ; μ = Λ ; ν ,
provided that Λ is identified with Θ as given by (8). Hence, the non-conservation of the energy-momentum tensor allows us to map the UG theory into GR equipped with a dynamical cosmological term, implying an interacting-like model in the GR context.
The mapping described above seems to indicate an equivalence of UG and GR. In fact, sometimes two or more apparently different theories may describe the same phenomena (in occurrence, gravity) but in a completely equivalent way. One example, is the so-called trinity of gravity involving Riemannian geometry, torsion and non-metricity, see Refs. [24,25,26]. Are the relations displayed above showing that UG is GR disguised? To answer this question, we must show that the equivalence is complete. As discussed below, and also in Ref. [16], there are indications that the equivalence between UG and GR is broken at the perturbative level.
It is also convenient, for reasons that will become clear later on in the work, to write down the UG equations in a more compact form such as
E μ ν = 8 π G τ μ ν ,
with the definitions,
E μ ν = R μ ν 1 4 g μ ν R = G μ ν + 1 4 g μ ν R ,
τ μ ν = T μ ν 1 4 g μ ν T .
We will call E μ ν the unimodular gravitational tensor and τ μ ν the unimodular energy-momentum tensor.

3. Equations for a Symmetric and Static Configuration

In this section, we will write down the general expressions for a symmetric and static configuration. In the Appendix A the corresponding expressions with a time dependence will be derived, which are necessary to perform the perturbative analysis to be described later.
Let us consider the metric,
d s 2 = e 2 γ d t 2 e 2 α d r 2 e 2 β d Ω 2 .
The non-vanishing Christoffel symbols are the following.
Γ 10 0 = γ , Γ 00 1 = e 2 ( γ α ) γ ,
Γ 11 1 = α , Γ 22 1 = e 2 ( β α ) β , Γ 33 1 = e 2 ( β α ) β sin 2 θ ,
Γ 12 2 = Γ 12 3 = β , Γ 33 2 = sin θ cos θ , Γ 23 3 = cot θ .
Primes denote derivatives with respect to r.
Also, the non-vanishing components of the Ricci tensor and the Ricci scalar are the following.
R 00 = e 2 ( γ α ) [ γ + γ ( γ + 2 β α ) ] ,
R 11 = γ 2 β + γ ( α γ ) + 2 β ( α β ) ,
R 22 = 1 e 2 ( β α ) [ β + β ( γ + 2 β α ) ] ,
R 33 = R 22 sin 2 θ ,
R = 2 e 2 β + 2 e 2 α [ γ + 2 β + 3 β 2 + γ ( γ + 2 β α ) 2 α β ] .
Consequently, the non-vanishing components of the unimodular gravitational tensor defined in (14) are the following:
E 00 = e 2 ( γ α ) [ γ 2 β 3 2 β 2 + γ 2 ( γ + 2 β α ) + β α ] + e 2 ( γ β ) 2 ,
E 11 = γ 2 β β 2 2 γ 2 ( γ α ) + β ( α + γ ) e 2 ( α β ) 2 ,
E 22 = 1 2 + e 2 ( β α ) 2 [ γ β 2 + γ ( γ α ) ] ,
E 33 = E 22 sin 2 θ .
The left-hand side of the UG field equations for the symmetric and static configuration has been set up. The next step is to characterize the source field. In the next couple of sections, the electromagnetic field and a scalar field will be considered as sources of the gravitational field.

4. The Electromagnetic Field

In the case of an electromagnetic field as the source of the energy-momentum tensor one has,
8 π G T μ ν E M = 2 { F μ ρ F ν ρ 1 4 g μ ν F ρ σ F ρ σ } .
It is worth mentioning that it has zero trace:
T E M = 0 .
Equations (5) and (6) become,
R μ ν 1 4 g μ ν R = 2 { F μ ρ F ν ρ 1 4 g μ ν F ρ σ F ρ σ } ,
F μ ρ ; μ F ν ρ = R ; ν 8 .
Remark that, contrary to GR, the traceless character of the energy-momentum tensor does not imply R = 0 , unless the Maxwell equations are obeyed.
Imposing the spherical symmetry, the only non-vanishing component is F 01 = E E ( r ) . Then, the equations are:
γ 2 β 3 2 β 2 + γ 2 ( γ + 2 β α ) + β α + e 2 ( α β ) 2 = e 2 γ + 4 α E 2 ,
γ 2 β β 2 2 + γ 2 ( α + 2 β γ ) + β α e 2 ( α β ) 2 = e 2 γ + 4 α E 2 ,
1 2 [ γ β 2 + γ ( γ α ) ] + e 2 ( α β ) 2 = e 2 γ + 4 α E 2 ,
( E 2 ) + 2 ( α + γ + 2 β ) E 2 = e 2 ( α + γ ) 4 R ,
with R given by (24).
Until now, no coordinate condition has been imposed. Adding (33) and (34), we obtain,
β ( α + γ ) β β 2 = 0 .
The use of the quasi-global coordinates, with α = γ , leads to,
β = log r .
As in the usual Reissner–Nordström (RN) solution in GR, there is a center at r = 0 . Equation (8), with T = 0 and with the identification of Θ with Λ implies in R = 4 Λ . The equations of motion reduce to,
γ + 2 γ 2 1 r 2 + e 2 γ r 2 = 2 e 2 γ E 2 ,
( E 2 ) + 4 E 2 r = Λ ,
In order to proceed further, we must impose a condition. This is a crucial step in working with UG as already stressed in [16]. One possibility is to fix the R = 4 Λ constant. This leads to the Reisnner–Nordström–de Sitter (RNdS) solution. In fact, this implies to recover the conservation law F μ ν ; μ = 0 . If Λ = 0 , we re-obtain the RN solution. If Λ > 0 , the RNdS solution is obtained, and if Λ < 0 , the Reisnner–Nordström-(Anti) de Sitter (RNAdS) solution is recovered, as it will be seen below. On the other hand, there are also also other possibilities that are to be explored since Λ can be nonconstant, covering the possibility of a dynamical cosmological term.
Three cases will be considered, namely a constant and two dynamical cosmological terms, corresponding to either the usual or the modified conservation laws.

4.1. Constant Cosmological Term

If Λ = constant,
E = Q r 2 ,
after identifying an integration constant with the total charge Q. The Coulomb law is recovered, as in the RN solution.
Using the quasi-global coordinate condition and the solution for the electric field E, Equation (33) becomes:
e 2 γ ( γ + 2 γ 2 ) e 2 γ r 2 = 1 r 2 + 2 Q 2 r 4 .
Defining A = e 2 γ , the equation takes the form,
A 2 A r 2 = 2 r 2 + 4 Q 2 r 4
This is a second-order, linear, non-homogeneous differential equation whose solution is
A = 1 + C 1 r + Q 2 r 2 + C 2 3 r 2 ,
C 1 , 2 being integration constants. Inserting this solution into the condition R = Λ , it is satisfied provided C 2 = Λ , while C 1 remains arbitrary, being fixed by using the Newtonian limit.
The final solution is given by,
A = 1 2 G M r + Q 2 r 2 Λ 3 r 2 .
This is the RNdS solution. It coincides with the static and spherically symmetric solution in GR with an electromagnetic field and a cosmological constant. This could be expected from the beginning since UG (satisfying the usual conservation laws) leads to the same field equations as GR with a cosmological term, with the only (but important, as we will see later) difference that UG is restricted to transverse diffeomorphism instead of the full diffeomorphism.

4.2. Varying Cosmological Term

For a varying cosmological term, it is necessary to impose an ansatz on the behavior of the function Λ . This is also true in GR when the cosmological term is dynamical. Since a static and spherically symmetric configuration is considered, the cosmological term must be a function on the coordinate r only: Λ Λ ( r ) .
Let us restrict ourselves again to the condition R = 4 Λ . Using the previous results and also identifying β = ln r , α = γ and A = e 2 γ , then:
A + 4 A r + 2 A r 2 = 2 r 2 4 Λ ( r ) .
The solution for the homogenous equation is,
A h = C 1 r + C 2 r 2 .
To obtain the inhomogeneous solution, we write,
A = f r 2 ,
obtaining,
f = 2 4 r 2 Λ ( r ) .
with a solution that depends on r:
f = r 2 4 [ r r 2 Λ ( r ) d r ] d r .
We will consider two different configurations for the function Λ ( r ) , corresponding to two distinct behaviors both asymptotically as well as at the center ( r = 0 ).

4.2.1. Case A

First, it imposes a power law behavior for Λ ( r ) ,
Λ ( r ) = Λ 0 + Λ 1 r p ,
with Λ 0 , 1 constants.
The final solution is given by the following expressions.
  • p 4 :
    A = 1 2 G M r + Q 2 r 2 Λ 0 3 r 2 4 Λ 1 ( p + 3 ) ( p + 4 ) r p + 2 ,
    E 2 = Q 2 r 4 p p + 4 Λ 1 r p ;
  • p = 4 :
    A = 1 2 G M r + Q 2 r 2 Λ 0 3 r 2 16 Λ 1 9 r { 3 ( ln r ) 2 + ln r } ,
    E 2 = Q 2 r 4 + 4 Λ 1 r 4 ln r .
The case p = 4 is clearly pathological since the electric field becomes imaginary near r = 0 when Λ 1 > 0 or for large r if Λ 1 < 0 . For p 4 a change of sign of E 2 can be avoided by choosing Λ 1 > 0 for 4 < p < 0 , or Λ 1 < 0 for p < 4 or p > 0 . The values p = 0 , 3 correspond to the cases already included in the constants C 1 and C 2 of the homogenous solution.
The solution p 0 , with the required conditions to avoid an imaginary electric field, contains either multiple horizon black holes, with a singularity at r = 0 , or naked singularities similar to the de Sitter–Reissner–Nordström (dSRN) solution in GR but the metric functions, in the UG case, may present a different shape but without change the general properties. These solutions are asymptotically non-flat except if Λ 0 = 0 and p > 2 . The corresponding equations in GR equipped with a cosmological term with the same functional dependence, using the same symmetries, lead to the same solution as it can be explicitly verified.

4.2.2. Case B

We will exploit now the functional form,
Λ ( r ) = Λ 0 + Λ 1 ( r 2 + a 2 ) 2 .
If Λ 0 = 0 , this functional form represents an asymptotically constant cosmological term near the origin, which becomes zero at infinity. Following the same steps as the previous case, the final form of the metric function is:
A = 1 2 G M r + Q 2 r 2 Λ 0 3 r 2 2 Λ 1 r 2 [ r a arctan r a ln ( 1 + r 2 a 2 ) ] ,
E 2 = Q 2 r 4 + Λ 1 { 1 ( r 2 + a 2 ) 2 + 2 [ 1 r 2 a 2 1 r 4 ln ( 1 + r 2 a 2 ) ] } .
Again, the same solution is obtained in the GR with a varying cosmological term given by (56). There are multiple horizons and naked singularities, as in the previous case. No change of sign in the E 2 term can be assured by imposing Λ 1 > 0 .
In all the cases discussed above, the presence of the cosmological term, constant or not, introduces new features in the solutions with respect to the usual Reissner–Nordström (RN) solution but does not remove the singularity at r = 0 . The behavior of this class of solution is depicted in Figure 1 for some parameter values.

5. Scalar Field

The energy-momentum tensor for a self-interacting scalar field is,
T μ ν = ϵ ( ϕ ; μ ϕ ; ν 1 2 g μ ν ϕ ; ρ ϕ ; ρ ) + g μ ν V ( ϕ ) .
The ordinary scalar field is denoted by ϵ = + 1 and the phantom scalar field by ϵ = 1 . In GR, in four dimensions, black holes exist only for the phantom case [27,28].
Inserting the expression for the energy-momentum tensor (59) in the UG equations one obtains,
R μ ν 1 4 g μ ν R = ϵ ( ϕ ; μ ϕ ; ν 1 4 g μ ν ϕ ; ρ ϕ ; ρ ) ,
R ; ν 4 = ϵ ( ϕ ; ν ϕ + ϕ ; ρ ϕ ; ν ; ρ 2 ) .
One distinguishing feature of the above equations is the absence of the potential (or, as before, the cosmological term, which is the particular case of a constant potential): it naturally disappears due to the traceless structure of the UG equations. Equation (61) can be written as,
( R 4 ϵ 4 ϕ ρ ϕ ; ρ ) ; ν = ϵ ϕ ; ν ϕ .
Identifying,
R 4 ϵ 4 ϕ ρ ϕ ; ρ = V ( ϕ ) ,
Equations (60) and (61) take the following form,
R μ ν 1 2 g μ ν R = ϵ ( ϕ ; μ ϕ ; ν 1 2 g μ ν ϕ ; ρ ϕ ; ρ ) + g μ ν V ( ϕ ) ,
ϕ = ϵ V ϕ ( ϕ ) .
In this way, we recover the GR equations equipped with a self-interacting scalar field.
In a static and spherically symmetric configuration, the UG field Equation (60) read,
γ 2 β 3 2 β 2 + γ 2 ( γ + 2 β α ) + β α + e 2 ( α β ) 2 = ϵ ϕ 2 4 ,
γ 2 β β 2 2 γ 2 ( γ α ) + β ( α + γ ) e 2 ( α β ) 2 = ϵ 3 4 ϕ 2 ,
1 2 [ γ β 2 + γ ( γ α ) ] + e 2 ( α β ) 2 = ϵ ϕ 2 4 .
Combining these equations, we have the following relations:
γ β 2 β 2 + γ ( γ + β α ) + α β + e 2 ( α β ) = 0 ,
β β + β ( γ + α ) = ϕ 2 2 .
Remark that in the UG equations, there is no potential, even if it appears in the energy-momentum tensor. Moreover, there are three metric functions (which can be reduced to two functions by gauging the radial coordinate) and the scalar field to be determined, and just two independent equations, (69) and (70). Hence an ansatz must be introduced. From the conservation law, we have the relation,
R + e α ϕ 2 = 4 V ( ϕ ) ,
where V ( ϕ ) is a function to be determined. We have slightly changed the notation (V instead of Λ ) to identify the unknown function with the potential. With this identification, the UG equations become identical to the GR equations with a potential. In GR, the potential must be chosen. In UG a functional form for the scalar field (or for one of the metric functions) must be chosen. There is a correspondence between the choice of the functional form of the scalar field and the choice of the potential in GR.
Two possible examples are the following.
  • If the scalar field is chosen such that,
    ϕ = ϵ C 2 k P ,
    P = 1 2 k ρ ,
    we find,
    d s 2 = P a d t 2 P a d ρ 2 P 1 a ρ 2 d Ω 2 ,
    a 2 = 1 ϵ C 2 k 2 .
    Here, ρ denotes the radial coordinate in the quasi-global coordinate system.
    Using (71) we find V = 0 . This solution represents a black hole only if ϵ = 1 . This solution has been determined in the GR context in Ref. [29].
  • The regular black hole determined in Ref. [28], is also the solution in the UG case, without a potential. Imposing that the scalar field is given by,
    ψ = ϕ 2 = arctan ρ b ,
    the metric is then given, in the quasi-global coordinates, by,
    d s 2 = A d t 2 d ρ 2 A r 2 ( ρ ) d Ω 2 ,
    A = 1 + c b 2 r 2 + ρ 0 b 3 ( b ρ + r 2 arctan ρ b ) .
    In these expressions b, c and ρ 0 are constants. Using relation (71) the potential in GR context is given by,
    V c b 3 ( 3 2 cos 2 ψ ρ 0 b 3 3 sin ψ cos ψ + ψ ( 3 2 cos 2 ψ ) ] ,
    the same used in the GR context.

6. Remarks on the Birkhoff Theorem and the Stability of the Solutions

Initially, we present a straightforward argument to demonstrate that, in the cases of both electromagnetism and scalar fields, the Birkhoff theorem holds the same significance as it does in GR. For electric-charged static solutions in GR, the Birkhoff theorem is valid. The same occurs for the corresponding solution in the UG. This can be seen by supposing radial time-dependent configurations. The argumentation follows the same reasoning used in GR, see for example [30]. From the expression presented in the Appendix, and considering only a radial electric field, the 0 1 component of the field equations, using the Schwarzschild coordinate system, with β = ln r , implies that α must be time-independent since the right-hand side of the equation is zero for a pure radial electric field. Combining equations 0 0 and 1 1 , it comes out that α = γ . Hence, all metric functions are time-independent.
For the scalar field case, the Birkhoff theorem is not valid because the right-hand side contains a term of the time ϕ ˙ ϕ which forbids considering the metric function α as time independent, as it happens in the GR case. The Birkhoff theorem is verified only if the scalar field is static [31].
In the two examples discussed in the previous section, having the scalar field as the source of the geometry, and considering the GR context, the solutions are unstable, except for the regular solution in the very special case where the minimum of the areal function coincides with the horizon [32,33]. However, this result can change in the UG context since the unimodular condition implies new relations for the perturbed functions that are absent in GR.
We will illustrate the special features of the perturbative analysis considering the case of black holes with a scalar field. Only radial perturbations will be considered. In the GR context, this is enough to conclude about the instability of the solution [32]. We will show that in the UG, if we try to follow the same procedure as in GR, the perturbations at first order are strictly zero due to the unimodular condition.
The unimodular condition implies,
g = det g μ ν = e α + γ + 2 β = ξ .
Since the function ξ is fixed, the unimodular condition leads, at linear perturbative order,
δ α + δ γ + 2 δ β = 0 .
There is still the freedom to impose a coordinate condition due to the diffeomorphic (even if transverse) invariance. The choice δ β = 0 is related to the gauge-invariant variables [32]. In fact, in GR this choice is equivalent to the full employment of gauge invariant variables and we will suppose that this property is also valid in UG. Hence, we end up with the conditions,
δ α = δ γ , δ β = 0 .
We write down the perturbations in a generic way as
δ f ( x , t ) = f ( x ) e i ω t .
The perturbed equations, under the conditions above, are the following:
δ γ + 4 γ δ γ { ω 2 . e 4 γ + 2 e 2 ( γ + β ) } δ γ = ϕ δ ϕ ,
δ γ + 4 γ δ γ { ω 2 . e 4 γ + 2 e 2 ( γ + β ) } δ γ = 3 ϕ δ ϕ ,
δ γ + 4 γ δ γ { ω 2 . e 4 γ + 2 e 2 ( γ + β ) } δ γ = ϕ δ ϕ ,
2 β δ γ = ϕ δ ϕ .
It is clear that the equations are consistent only in the trivial case: δ ϕ = δ γ = 0 . Hence, it is not possible to obtain information on the stability of the solution, at least at the linear level and follow a procedure close to that used in GR. This indicates a distinguishing feature of unimodular gravity in comparison with GR.
However, some cautions are necessary. In cosmology, the Newtonian gauge is equivalent to the gauge-invariant formalism in GR, in the sense that, at least in the absence of anisotropic stress, the final equations are the same [34]. In Ref. [17] it has been argued that, in cosmology and in the UG context, the gauge-invariant formalism can be applied while the Newtonian gauge is not. Does the same happen here, when static, spherically symmetric perturbations are considered? In our point of view, this asks for a separate analysis, especially if we take into account the point of view presented in Ref. [35]. We must remember, on the other hand, that the same solution may appear in different theories, but with different perturbative behavior, see an example in Ref. [36] where the same wormhole solution has been analyzed at background and perturbative levels in three different theoretical contexts.

7. Conclusions

Unimodular gravity (UG) is one of the first alternatives to General Relativity (GR). It is a geometric theory that is invariant with respect to a restricted diffeomorphic class of transformations, the transverse diffeomorphism, due to the imposition of a constraint on the determinant of the metric. In UG the usual conservation of the energy-momentum tensor is not assured: The conservation of the energy-momentum tensor is a choice. If it is imposed, UG becomes in principle equivalent to GR with a cosmological term. However, the restriction on the determinant of the metric may lead to some important new features at the perturbative level. We have shown here that if the conservation of the energy-moment tensor is relaxed, UG becomes equivalent to GR with a dynamical cosmological term, with still the same important difference due to the UG constraint which can manifested at the perturbative level.
We have discussed, in this context, the static and spherically symmetric solutions in UG. For the vacuum configuration, the Schwarzschild solution is also verified in UG. The same occurs with the Reissner–Nordström solution, but only if the energy-momentum tensor is conserved. If not, the dynamical cosmological term induces new features, but it does not prevent the appearance of the singularity at r = 0 . Similar features appear in the case when a scalar field appears as the main source. In this case, the potential term, representing the self-interaction of the scalar field, disappears in the UG context and an ansatz must be imposed in order to close the set of equations. This mounts, in the GR context, to choose a given potential for the scalar field. For a discussion of the UG in static, spherical configurations but focusing on compact objects, see Ref. [37].
We have shown that the Birkhoff theorem follows the same features as in GR, being satisfied for a charged solution, being possibly violated for a dynamical scalar field. The linear radial perturbations have been analyzed when a scalar field is present. Once more, GR black hole solutions are generically unstable in the latter case. In UG, using the gauge invariant approach employed in GR and restricting to radial perturbations, the condition on the determinant of the metric leads to vanishing perturbations at linear order, and possibly also for higher order. As already discussed in the cosmological context, this result seems to point out a breaking of the equivalence of UG and GR at the perturbative level. There are other viewpoints on the implementation of the UG constraints in performing a perturbative analysis, see for example Ref. [35]. However, the results reported here indicate that a direct application of the procedures used in GR combined with the unimodular constraint may lead to conclusions different from those obtained in GR.
In fact, the main open problem related to the results obtained above refers to the perturbative analysis. As pointed out in previous works (see Ref. [16] and references therein) the main difference between GR and UG is related to the perturbative analysis. While the UG gravity can be mapped into GR with a cosmological constant, when the energy-momentum tensor is conserved, or in a dynamical vacuum model, when the energy-momentum tensor is not conserved, previous investigations indicate that such mapping seems not to be complete when perturbations are considered even at the perturbative level. We have given an example here, applying directly the relations coming from the unimodular constraints besides the usual coordinate condition. However, there are claims (for example, [35]) that this issue is not closed and, for example, all coordinate conditions can be used in UG in spite of the unimodular constraint in the determinant of the metric. We intend to present in the future a more concrete and detailed analysis of this problem.

Author Contributions

Conceptualization, J.C.F. and H.V.; methodology, J.C.F., H.V. and M.H.D.; validation, M.H.D.; investigation, J.C.F., H.V. and M.H.D.; writing—original draft preparation, J.C.F., H.V. and M.H.D. All authors have read and agreed to the published version of the manuscript.

Funding

We thank CNPq, FAPES and FAPEMIG for partial financial support.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We thank K.A. Bronnikov for enlighting discussions on some aspects of the problem treated in this work and L.F. de Oliveira Guimarães for his remarks on the text.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Spherically Symmetric Non Static Metric

A dynamical spherically symmetric metric, admitting radial oscillations, is given by,
d s 2 = e 2 γ ( t , r ) d t 2 e 2 α ( t , r ) d r 2 e 2 β ( t , r ) d Ω 2 .
The non-vanishing Christoffel symbols are the following.
Γ 00 0 = γ ˙ , Γ 10 0 = γ , Γ 11 0 = e 2 ( α γ ) α ˙ ,
Γ 22 0 = e 2 ( β γ ) β ˙ , Γ 33 0 = Γ 22 0 sin 2 θ ,
Γ 00 1 = e 2 ( γ α ) γ , Γ 01 1 = α ˙ , Γ 11 1 = α ,
Γ 22 1 = e 2 ( β α ) β , Γ 33 1 = Γ 22 1 sin 2 θ ,
Γ 02 2 = Γ 03 3 = β ˙ , Γ 12 2 = Γ 12 3 = β ,
Γ 33 2 = sin θ cos θ , Γ 23 3 = cot θ .
The dots mean derivatives with respect to t and the primes with respect to r.
The non-vanishing components of the Ricci tensor and the Ricci scalar are the following.
R 00 = α ¨ 2 β ¨ + γ ˙ ( α ˙ + 2 β ˙ ) α ˙ 2 2 β ˙ 2 + e 2 ( γ α ) [ γ + γ ( γ + 2 β α ) ] ,
R 11 = e 2 ( α γ ) { α ¨ + α ˙ ( α ˙ γ ˙ + 2 β ˙ ) } γ 2 β + γ ( α γ ) + 2 β ( α β ) ,
R 22 = 1 + e 2 ( β γ ) [ β ¨ + β ˙ ( α ˙ + 2 β ˙ γ ˙ ) ] e 2 ( β α ) [ β + β ( γ + 2 β α ) ] ,
R 33 = R 22 sin 2 θ ,
R 01 = 2 { β ˙ + β ˙ ( γ β ) + α ˙ β } ,
R = 2 e 2 β + 2 e 2 α [ γ + 2 β + 3 β 2 + γ ( γ + 2 β α ) 2 α β ] 2 e 2 γ [ α ¨ + 2 β ¨ + 3 β ˙ 2 + α ˙ ( α ˙ + 2 β ˙ γ ˙ ) 2 γ ˙ β ˙ ]
The non-vanishing components of the unimodular gravitational tensor,
E μ ν = R μ ν 1 4 g μ ν R ,
are the following:
E 00 = α ¨ 2 β ¨ + γ ˙ 2 ( α ˙ + 2 β ˙ ) α ˙ 2 2 β ˙ 2 2 + α ˙ β ˙ + e 2 ( γ α ) [ γ 2 β 3 2 β 2 + γ 2 ( γ + 2 β α ) + β α ] + e 2 ( γ β ) 2 ,
E 11 = e 2 ( α γ ) [ α ¨ 2 β ¨ 3 2 β ˙ 2 + α ˙ 2 ( α ˙ + 2 β ˙ γ ˙ ) + γ ˙ β ˙ ] γ 2 β β 2 2 γ 2 ( γ α ) + β ( α + γ ) e 2 ( α β ) 2 ,
E 22 = 1 2 e 2 ( β γ ) 2 [ α ¨ β ˙ 2 + α ˙ ( α ˙ γ ˙ ) ] + e 2 ( β α ) 2 [ γ β 2 + γ ( γ α ) ] ,
E 01 = 2 [ β ˙ + β ˙ ( γ β ) + α ˙ β ] ,
E 33 = E 22 sin 2 θ .
For the static case, the above expressions reduce to,
E 00 = e 2 ( γ α ) [ γ 2 β 3 2 β 2 + γ 2 ( γ + 2 β α ) + β α ] + e 2 ( γ β ) 2 ,
E 11 = γ 2 β β 2 2 γ 2 ( γ α ) + β ( α + γ ) e 2 ( α β ) 2 ,
E 22 = 1 2 + e 2 ( β α ) 2 [ γ β 2 + γ ( γ α ) ] ,
E 33 = E 22 sin 2 θ .

References

  1. Brandenberger, R. Perspectives on the Dark Sector. arXiv 2022, arXiv:2211.11273. [Google Scholar]
  2. Perivolaropoulos, L.; Skara, F. Challenges for ΛCDM: An update. arXiv 2022, arXiv:2105.05208. [Google Scholar] [CrossRef]
  3. Weinberg, S. The cosmological constant problem. Rev. Mod. Phys. 1989, 61, 1. [Google Scholar] [CrossRef]
  4. Sahni, V.; Starobinsky, A.A. The Case for a positive cosmological Lambda term. Int. J. Mod. Phys. D 2000, 9, 373–443. [Google Scholar] [CrossRef]
  5. Padmanabhan, T. Cosmological constant—The weight ofthe vacuum. Phys. Rep. 2003, 380, 235–320. [Google Scholar] [CrossRef]
  6. Zlatev, I.; Wang, L.-M.; Steinhardt, P.J. Quintessence, Cosmic Coincidence, and the Cosmological Constant. Phys. Rev. Lett. 1999, 82, 896. [Google Scholar] [CrossRef]
  7. Arkani-Hamed, N.; Dimopoulos, S.; Dvali, G.; Gabadadze, G. Non-local modification of gravity and the cosmological constant problem. arXiv 2002, arXiv:hep-th/0209227. [Google Scholar]
  8. Dvali, G.; Hofmann, S.; Khoury, J. Degravitation of the cosmological constant and graviton width. Phys. Rev. D 2007, 76, 084006. [Google Scholar] [CrossRef]
  9. de Rham, C.; Hofmann, S.; Khoury, J.; Tolley, A.J. Cascading gravity and degravitation. J. Cosmol. Astropart. Phys. 2008, 802, 11. [Google Scholar] [CrossRef]
  10. Amendola, L.; Tsujikawa, S. Dark Energy: Theory and Observations; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  11. Einstein, A. Spielen Gravitationsfelder im Aufber der materiellen Elementarteilchen eine wessentliche Rolle? (Sizungsberichte der Preussisschen Akad. d. Wissenschaften, 1919). English translation Do gravitational fields play an essential part in the structure of the elementary particles of matter? In The Principle of Relativity; Volume 7: The Berlin Years: Writings, 1918–1921; Dover Publications: New York, NY, USA, 1952. [Google Scholar]
  12. Alvarez, E.; Anero, J. Unimodular cosmological models. arXiv 2021, arXiv:2109.08077. [Google Scholar]
  13. de Brito, G.P.; Melichev, O.; Percacci, R.; Pereira, A.D. Can quantum fluctuations differentiate between standard and unimodular gravity? J. High Energy Phys. 2021, 12, 90. [Google Scholar] [CrossRef]
  14. Carballo-Rubio, R.; Garay, L.J.; García-Moreno, G. Unimodular gravity vs general relativity: A status report. Class. Quantum Grav. 2022, 39, 243001. [Google Scholar] [CrossRef]
  15. Lopez-Villarejo, J.J. Transverse Diff gravity is to scalar-tensor as unimodular gravity is to General Relativity. J. Cosmol. Astropart. Phys. 2011, 11, 2. [Google Scholar] [CrossRef]
  16. Fabris, J.C.; Alvarenga, M.H.; Velten, H. Using Cosmological Perturbation Theory to Distinguish between GR and Unimodular Gravity. Symmetry 2023, 15, 1392. [Google Scholar] [CrossRef]
  17. Gao, C.; Brandenberger, R.H.; Cai, Y.; Chen, P. Cosmological perturbations in unimodular gravity. J. Cosmol. Astropart. Phys. 2014, 9, 21. [Google Scholar] [CrossRef]
  18. Velten, H.; Caramês, T.R.P. To conserve, or not to conserve: A review of nonconservative theories of gravity. Universe 2021, 7, 38. [Google Scholar] [CrossRef]
  19. Wald, R.M. General Relativity; Chigago University Press: Chicago, IL, USA, 1984. [Google Scholar]
  20. Peracaula, J.S. The dynamics of vacuum, gravity and matter: Implications on the fundamental constants. arXiv 2023, arXiv:2308.13349. [Google Scholar]
  21. Mavromatos, N.E.; Peracaula, J.S.; Gómez-Valent, A. String-inspired running-vacuum cosmology, quantum corrections and the current cosmological tensions. arXiv 2023, arXiv:2307.13130. [Google Scholar]
  22. Gómez-Valent, A.; Mavromatos, N.E.; Peracaula, J.S. Stringy Running Vacuum Model and current Tensions in Cosmology. arXiv 2023, arXiv:2305.15774. [Google Scholar] [CrossRef]
  23. Panotopoulos, G.; Rincon, A. Quasinormal spectra of scale-dependent Schwarzschild-de Sitter black holes. arXiv 2021, arXiv:2011.02860. [Google Scholar] [CrossRef]
  24. Heisenberg, L. A systematic approach to generalisations of General Relativity and their cosmological implication. Phys. Rep. 2019, 796, 1–113. [Google Scholar] [CrossRef]
  25. Jiménez, J.B.; Heisenberg, L.; Koivisto, T.S. The geometrical trinity of gravity. Universe 2019, 5, 173. [Google Scholar] [CrossRef]
  26. Heisenberg, L. Review on f(Q) Gravity. arXiv 2023, arXiv:2309.15958. [Google Scholar]
  27. Xanthopoulos, B.C.; Zannias, T. Einstein gravity coupled to a massless scalar field in arbitrary spacetime dimensions. Phys. Rev. D 1989, 40, 2564. [Google Scholar] [CrossRef] [PubMed]
  28. Bronnikov, K.; Fabris, J.C. Regular Phantom Black Holes. Phys. Rev. Lett. 2006, 96, 251101. [Google Scholar] [CrossRef] [PubMed]
  29. Bronnikov, K.A.; Chernakova, M.S.; Fabris, J.C.; Pinto-Neto, N.; Rodrigues, M.E. Cold black holes and conformal continuations. Int. J. Mod. Phys. D 2008, 17, 25. [Google Scholar] [CrossRef]
  30. Weinberg, S. Gravitation and Cosmology; Wiley: New York, NY, USA, 1972. [Google Scholar]
  31. Bronnikov, K.A.; Kovalchuk, M.A. On a generalisation of Birkhoff’s theorem. J. Phys. A 1980, 13, 187. [Google Scholar] [CrossRef]
  32. Bronnikov, K.A.; Fabris, J.C.; Zhidenko, A. On the stability of scalar-vacuum space-times. Eur. Phys. J. C 2011, 71, 1791. [Google Scholar] [CrossRef]
  33. Bronnikov, K.A.; Konoplya, R.A.; Zhidenko, A. Instabilities of wormholes and regular black holes supported by a phantom scalar field. Phys. Rev. D 2012, 86, 024028. [Google Scholar] [CrossRef]
  34. Mukhanov, V.F.; Feldman, H.A.; Brandenberger, R.H. Theory of cosmological perturbations. Phys. Rep. 1992, 215, 203–333. [Google Scholar] [CrossRef]
  35. Bengochea, G.R.; Leon, G.; Perez, A.; Sudarsky, D. A clarification on prevailing misconceptions in unimodular gravity. J. Cosmol. Astropart. Phys. 2023, 2023, 011. [Google Scholar] [CrossRef]
  36. Bronnikov, K.A.; Barcellos, V.A.G.; de Carvalho, L.P.; Fabris, J.C. The simplest wormhole in Rastall and k-essence theories. Eur. Phys. J. C 2021, 81, 395. [Google Scholar] [CrossRef]
  37. Astorga-Moreno, J.A.; Chagoy, J.; Flores-Urbinad, J.C.; García-Aspeitia, M.A. Compact objects in unimodular gravity. J. Cosmol. Astropart. Phys. 2019, 9, 5. [Google Scholar] [CrossRef]
Figure 1. Solution of the case B using Λ 0 = 0 , Λ 1 = 1 / 2 , with m = 1 , Q = 0.5 (Dashed), and Q = 1.5 (Dotted). The continuous line is the RN solution with Λ 0 = 0 , Λ 1 = 0 , with m = 1 and Q = 0.5 .
Figure 1. Solution of the case B using Λ 0 = 0 , Λ 1 = 1 / 2 , with m = 1 , Q = 0.5 (Dashed), and Q = 1.5 (Dotted). The continuous line is the RN solution with Λ 0 = 0 , Λ 1 = 0 , with m = 1 and Q = 0.5 .
Universe 09 00515 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fabris, J.C.; Daouda, M.H.; Velten, H. Spherically Symmetric Configurations in Unimodular Gravity. Universe 2023, 9, 515. https://doi.org/10.3390/universe9120515

AMA Style

Fabris JC, Daouda MH, Velten H. Spherically Symmetric Configurations in Unimodular Gravity. Universe. 2023; 9(12):515. https://doi.org/10.3390/universe9120515

Chicago/Turabian Style

Fabris, Júlio César, Mahamadou Hamani Daouda, and Hermano Velten. 2023. "Spherically Symmetric Configurations in Unimodular Gravity" Universe 9, no. 12: 515. https://doi.org/10.3390/universe9120515

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