Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Stochastic Chaos and Markov Blankets
Next Article in Special Issue
Inhomogeneous and Radiating Composite Fluids
Previous Article in Journal
Learning in Convolutional Neural Networks Accelerated by Transfer Entropy
Previous Article in Special Issue
Phase Transition in Modified Newtonian Dynamics (MONDian) Self-Gravitating Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hyperbolically Symmetric Versions of Lemaitre–Tolman–Bondi Spacetimes

1
Instituto Universitario de Física Fundamental y Matemáticas, Universidad de Salamanca, 37007 Salamanca, Spain
2
Escuela de Física, Facultad de Ciencias, Universidad Central de Venezuela, Caracas 1050, Venezuela
3
Departamento de Matemáticas Aplicada and Instituto Universitario de Física Fundamental y Matemáticas, Universidad de Salamanca, 37007 Salamanca, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Entropy 2021, 23(9), 1219; https://doi.org/10.3390/e23091219
Submission received: 30 August 2021 / Revised: 13 September 2021 / Accepted: 14 September 2021 / Published: 16 September 2021
(This article belongs to the Special Issue Complexity of Self-Gravitating Systems)

Abstract

:
We study fluid distributions endowed with hyperbolic symmetry, which share many common features with Lemaitre–Tolman–Bondi (LTB) solutions (e.g., they are geodesic, shearing, and nonconformally flat, and the energy density is inhomogeneous). As such, they may be considered as hyperbolic symmetric versions of LTB, with spherical symmetry replaced by hyperbolic symmetry. We start by considering pure dust models, and afterwards, we extend our analysis to dissipative models with anisotropic pressure. In the former case, the complexity factor is necessarily nonvanishing, whereas in the latter cases, models with a vanishing complexity factor are found. The remarkable fact is that all solutions satisfying the vanishing complexity factor condition are necessarily nondissipative and satisfy the stiff equation of state.
PACS:
04.40.-b; 04.20.-q; 04.40.Dg; 04.40.Nr

1. Introduction

In a recent paper [1], we presented a general study on the dynamics of hyperbolic symmetric fluids (DHSFs). Our main motivation behind such an endeavor (but not the only one) was to describe the dynamic regime preceding the final equilibrium state of static hyperbolic symmetric fluids described in [2], which in its turn could be used to model the source of the hyperbolic symmetric black hole described in [3,4] where the region interior to the horizon is described by the line element:
d s 2 = 2 M R 1 d t 2 + d R 2 2 M R 1 + R 2 d Ω 2 , d Ω 2 = d θ 2 + sinh 2 θ d ϕ 2 .
The rationale behind such a proposal stems from the well-known fact that any transformation that maintains the static form of the Schwarzschild metric (in the whole spacetime) is unable to remove the coordinate singularity in the line element [5]. In other words, the regular extension of the Schwarzschild metric to the whole spacetime (including the region inner to the horizon) may be achieved, but at the price of admitting a nonstatic spacetime inside the horizon [6,7].
Since any dynamic regime should eventually end in an equilibrium final state, it would be desirable to have a static solution over the whole spacetime.
Thus, the model proposed in [3] describes the space time as consisting of two four-dimensional manifolds, one described by the usual Schwarzschild metric on the exterior side of the horizon and a second one in the interior of it, described by (1).
The metric (1) is a static solution admitting the four Killing vectors:
K ( 0 ) = t ,
and:
K ( 2 ) = cos ϕ θ + coth θ sin ϕ ϕ , K ( 1 ) = ϕ , K ( 3 ) = sin ϕ θ + coth θ cos ϕ ϕ .
The above Killing vectors (3) define the hyperbolic symmetry. Solutions to the Einstein equations endowed with this type of symmetry have been the subject of research in different contexts (see [8,9,10,11,12,13,14,15,16] and the references therein).
Besides the general properties of DHSFs analyzed in [1], some exact solutions were found. Particular attention is owed to two nondissipative solutions, which could be regarded as the hyperbolic symmetric versions of the Friedman–Robertson–Walker spacetime, since they are shear-free, nondissipative, and conformally flat. All solutions presented in [1] satisfy the condition of the vanishing complexity factor and evolve in the quasi-homologous regime. This last condition, in the nondissipative case, implies that the fluid is shear-free, thereby excluding the possibility to obtain a hyperbolic symmetric version of the LTB spacetimes under the above-mentioned conditions.
Due to the huge relevance of LTB spacetimes, we devote this work to study in some detail its possible hyperbolic symmetric versions. To do that, we must abandon the condition of quasi-homologous evolution.
It is worth recalling that LTB dust models [17,18,19] are among the most appealing solutions to the Einstein equations. They describe the spherically symmetric distribution of inhomogeneous nondissipative dust (see [20,21] for a detailed description of these spacetimes). Although LTB spacetimes are usually associated with an inhomogeneous dust source, it is known that the most general source compatible with LTB spacetimes is an anisotropic fluid [20,22].
LTB spacetimes have been invoked as cosmological models (see [23,24,25] and the references therein), in the study of gravitational collapse, when dealing with the problem of cosmic censorship [26,27,28,29,30,31], and in quantum gravity [32,33].
The apparent accelerated expansion of the universe, as inferred from some observations of type Ia supernovae, has renewed the interest in LTB spacetimes. Indeed, even though there is general consensus to invoke dark energy as a source of antigravity for understanding the cosmic acceleration, a growing number of researchers are now considering that inhomogeneities may account for the observed cosmic acceleration, without resorting to dark energy (see [34,35,36,37,38,39] and the references therein).
In this work, we present several models that could be considered as hyperbolic symmetric versions of LTB spacetimes, in the sense that they are geodesic, shearing, and nonconformally flat and the energy density is inhomogeneous, but with the spherical symmetry replaced by the hyperbolic symmetry.
We consider both nondissipative and dissipative models. The general approach used for reaching our goal was already outlined in [1]; however, for the sake of completeness, we present the basic steps in the following sections. The last section is devoted to the discussion of the obtained results.

2. Fluid Distribution, Kinematical Variables, and Basic Equations

We consider hyperbolic symmetric distributions of a geodesic fluid, which may be bounded (or not) from outside by a surface Σ e . As we already know (see [1] for a detailed discussion on this point), hyperbolic symmetric fluids cannot fill the central region, and therefore, such a region should be described either by an empty vacuole or by a fluid distribution not endowed with hyperbolic symmetry. Thus, our fluid is also bounded from inside by a surface Σ i .
Choosing comoving coordinates, the general metric can be written as:
d s 2 = d t 2 + B 2 d r 2 + R 2 ( d θ 2 + sinh 2 θ d ϕ 2 ) ,
where B and R are assumed positive and due to the symmetry defined by (3) are functions of t and r. We number the coordinates x 0 = t , x 1 = r , x 2 = θ , and x 3 = ϕ .
The general energy momentum tensor T α β of the fluid distribution may be written as:
T α β = ( μ + P ) V α V β + P g α β + ( P r P ) χ α χ β + q α V β + V α q β ,
where μ is the energy density, P r the radial pressure, P the tangential pressure, and q α the heat flux; these physical variables, due to the symmetry defined by (3), are functions of t and r. Furthermore, V α and χ α denote the four-velocity of the fluid and a unit four-vector along the radial direction, respectively; they satisfy:
V α V α = 1 , V α q α = 0 , χ α χ α = 1 , χ α V α = 0 .
Since we consider comoving observers, we have:
V α = δ 0 α , q α = q B 1 δ 1 α , χ α = B 1 δ 1 α .

2.1. Einstein Equations

For (4) and (5), the Einstein equations:
G α β = 8 π T α β ,
read:
8 π μ = 2 B ˙ B + R ˙ R R ˙ R 1 B 2 2 R R + R R 2 2 B B R R + B R 2 ,
4 π q = 1 B R ˙ R B ˙ B R R ,
8 π P r = 2 R ¨ R + R ˙ R 2 + R B R 2 + 1 R 2 ,
8 π P = B ¨ B + R ¨ R + B ˙ B R ˙ R + 1 B 2 R R B B R R ,
where dots and primes denote derivatives with respect to t and r, respectively. It is worth stressing the difference between these equations and those corresponding to the spherically symmetric LTB case.

2.2. Kinematical Variables and the Mass Function

The expansion Θ is given by:
Θ = V α ; α = B ˙ B + 2 R ˙ R ,
and for the shear, we have (remember that the four-acceleration and the vorticity vanish):
σ α β = V ( α ; β ) 1 3 Θ h α β ,
where h α β = g α β + V α V β .
The nonvanishing components of (14) are:
σ 11 = 2 3 B 2 σ , σ 22 = σ 33 sinh 2 θ = 1 3 R 2 σ ,
with:
σ α β σ α β = 2 3 σ 2 ,
being:
σ = B ˙ B R ˙ R ,
σ α β may also be written as:
σ α β = σ χ α χ β 1 3 h α β .
Next, the mass function m ( t , r ) introduced by Misner and Sharp [40] (see also [41]) is given by:
m = R 2 R 232 3 = R 2 R ˙ 2 + R B 2 + 1 ,
where the components R 232 3 of the Riemann tensor are calculated with (4).
Defining as usual the “areal” velocity U of the fluid as the variation of R with respect to proper time, i.e.,
U = R ˙ ,
then since U < 1 , it follows at once from (19) that m is a positive defined quantity.
With the above, we can express (19) as:
E R B = U 2 + 2 m R 1 1 / 2 .
From (19) and the field equations, we obtain:
m ˙ = 4 π R 2 P r U + q E ,
and:
m = 4 π R R 2 μ + q U E .
The integration of (23) produces:
m = 0 r 4 π R 2 μ + q U E R d r ,
whose partial integration yields:
3 m R 3 = 4 π μ + 4 π R 3 0 r R 3 μ 3 q U B R d r .
Then, it follows from (24) that μ is necessarily negative, if we assume the condition R > 0 to avoid shell crossing and remind that m > 0 .
Furthermore, it follows from (24) that whenever the energy density is regular, then m r 3 as r tends to zero. However, in this same limit, U 0 and R r , implying because of (21) that the central region cannot be filled with our fluid distribution. Among the many possible scenarios, we shall assume here that the center is surrounded by a vacuum cavity. However, this is just one of the possible choices, which even if having implications on specific models, does not affect the general properties of the fluids endowed with hyperbolic symmetry.

2.3. The Exterior Spacetime and Junction Conditions

In the case of bounded configurations, we assume that outside Σ e , we have the hyperbolic symmetric version of the Vaidya spacetime, described by:
d s 2 = 2 M ( v ) ρ 1 d v 2 2 d ρ d v + ρ 2 ( d θ 2 + sinh 2 θ d ϕ 2 ) ,
where M ( v ) denotes the total mass and v is the retarded time.
Now, from the continuity of the first differential form, it follows (see [42] for details),
d t = Σ e d v 2 M ( v ) ρ 1 ,
R = Σ e ρ ( v ) ,
and:
d v d t 2 = Σ e 2 M ( v ) ρ 1 + 2 d ρ d v ,
whereas the continuity of the second differential form produces:
m ( t , r ) = Σ e M ( v ) ,
and:
P r = Σ e q ,
where = Σ e means that both sides of the equation are evaluated on Σ e .
The corresponding junction conditions on Σ i are:
m ( t , r ) = Σ i 0 ,
and:
P r = Σ i 0 .
When either of the above conditions cannot be satisfied, we have to admit the presence of thin shells.

2.4. Weyl Tensor

The Weyl tensor is defined through the Riemann tensor R α β μ ρ , the Ricci tensor R α β , and the curvature scalar R as:
C α β μ ρ = R α β μ ρ 1 2 R β ρ g α μ + 1 2 R α β δ μ ρ 1 2 R α μ δ β ρ
+ 1 2 R μ ρ g α β + 1 6 R ( δ β ρ g α μ g α β δ μ ρ ) .
In our case, the magnetic part of the Weyl tensor vanishes, whereas its electric part, defined by:
E α β = C α μ β ν V μ V ν ,
has the following nonvanishing components:
E 11 = 2 3 B 2 E , E 22 = 1 3 R 2 E , E 33 = E 22 sinh 2 θ ,
where:
E = 1 2 R ¨ R B ¨ B R ˙ R B ˙ B R ˙ R + 1 2 B 2 R R + B B + R R R R + 1 2 R 2 .
Observe that we may also write E α β as:
E α β = E ( χ α χ β 1 3 h α β ) .
Finally, using (9), (11) and (12) with (19) and (37), we obtain:
3 m R 3 = 4 π μ + 4 π ( P r P ) + E .

3. Structure Scalars and Complexity Factor

Some of the models exhibited below are obtained from the conditions imposed on a scalar function that appears in a natural way in the orthogonal splitting of the Riemann tensor (see [43] for details) and that is identified as the complexity factor.
Thus, let us introduce the tensor Y α β (which is an element of that splitting [44,45,46,47]), defined by:
Y α β = R α γ β δ V γ V δ .
Tensor Y α β may be expressed in terms of two scalar functions Y T , Y T F (structure scalars) as:
Y α β = 1 3 Y T h α β + Y T F ( χ α χ β 1 3 h α β ) ,
where:
Y T = 4 π ( μ + 3 P r 2 Π ) , Y T F = E 4 π Π ,
with Π = P r P .
Combining (39) with (25) and (42) produces:
Y T F = 8 π Π + 4 π R 3 0 r R 3 μ 3 q U B R d r .
The complexity factor is a scalar function intended to measure the complexity of a given fluid distribution (see [48,49,50] for details). For static hyperbolic symmetric fluids (as well as for spherically symmetric ones), the complexity factor is identified with the scalar function Y T F defined, in the dynamic case, by Equations (42) and (43) (see [1]). The main reason behind such a proposal resides, on the one hand, in the basic assumption that one of the less complex systems corresponds to a homogeneous (in the energy density) fluid distribution with isotropic pressure. Thus, any variable measuring complexity should vanish for this specific case. On the other hand, the scalar function Y T F contains contributions from the energy density inhomogeneity and the local pressure anisotropy, combined in a very specific way, which (in the static case) vanishes for the homogeneous and locally isotropic fluid distribution. Furthermore, this scalar measures the departure of the value of the Tolman mass for the homogeneous and isotropic fluid, produced by the energy density inhomogeneity and the pressure anisotropy.
It is worth mentioning that the complexity factor so defined not only vanishes for the simple configuration mentioned above, but also may vanish when the terms appearing in its definition cancel each other. Thus, vanishing complexity may correspond to very different systems.
In the time-dependent case, we face two different problems: on the one hand, we have to generalize the concept of the complexity of the structure of the fluid distribution to time-dependent dissipative fluids, and on the other hand, we also have to evaluate the complexity of the mode of evolution. Following the strategy outlined in [49], the complexity factor for the dissipative case of the fluid distribution is assumed to be the function Y T F , as in the static case, which now includes the dissipative variables. With respect to the complexity of the mode of evolution, let us recall that in the past, the homologous and quasi-homologous conditions have been used to characterize the simplest mode of evolution. However, we know that in the nondissipative case, the homologous and the quasi-homologous conditions imply the vanishing of the shear (see Equation (59) in [1]), and therefore, we shall not adopt such restrictions here.

4. Hyperbolically Symmetric Lemaitre–Tolman–Bondi Metric: The Nondissipative Dust Case

We start our search of hyperbolic symmetric versions of LTB (HSLTB) by considering the simplest case, i.e., we assume nondissipative geodesic dust. Extensions to the dissipative, anisotropic case are discussed in the next section, along the lines developed in [51].
Under the above-mentioned conditions, we find after integration of (10):
B ( t , r ) = R k ( r ) 1 1 / 2 ,
where k is an arbitrary function of r.
Then, from (19) and (44), it follows:
R ˙ 2 = 2 m R + k ( r ) .
Equation (45) implies k ( r ) > 2 m R . Thus, unlike the spherically symmetric LTB spacetime, we now have only one case k ( r ) > 0 .
The solution to (45) may be written as:
R = m k ( cosh η + 1 ) , m k 3 / 2 ( sinh η + η ) = t t 0 ( r ) ,
where t 0 ( r ) is an integration function of r.
Thus, for the line element, we have:
d s 2 = d t 2 + ( R ) 2 k ( r ) 1 d r 2 + R 2 ( d θ 2 + sinh 2 θ d ϕ 2 ) .
In order to prescribe an explicit model, we have to provide the three functions k ( r ) , m ( r ) , and t 0 ( r ) . However, since (47) is invariant under transformations of the form r = r ( r ˜ ) , we only need two functions of r.
Assuming m 0 = m k = c o n s t a n t and t 0 ( r ) = c o n s t a n t , the expressions for Θ and σ read:
Θ = k m 0 sinh η sinh η + η + cosh η + 2 sinh η + 1 ( cosh η + 1 ) 2 ,
σ = k m 0 sinh η sinh η + η + cosh η sinh η + 1 ( cosh η + 1 ) 2 ,
from which it is clear that the expansion is always positive.
Since, as we have already mentioned, our fluid distribution cannot reach the central region, then we do not need to consider any regularity conditions there.
The only nontrivial conservation law in this case reads:
μ ˙ + μ Θ = 0 ,
or:
μ ˙ + μ B ˙ B + 2 R ˙ R = 0 ,
producing:
μ = h ( r ) B R 2 ,
or, using (44):
μ = 3 h ( r ) k ( r ) 1 1 / 2 ( R 3 ) ,
where h ( r ) is a function of integration, which due to the fact that the energy density is negative, must be necessarily negative.
Scalar Y T F for (47) reads:
Y T F = R ¨ R R ¨ R .
As is evident from (43), since we are considering nondissipative inhomogeneous dust, the complexity factor Y T F cannot vanish. However, this situation may change in the dissipative, anisotropic case, as we see in the models exhibited below. On the other hand, since in the nondissipative case, the quasi-homologous condition implies the vanishing of the shear (see Equation (59) in [1]), we have to abandon such a restriction for our models.

5. Dissipative Case

We now consider the possibility that the system radiates, and the pressure is nonvanishing and may be anisotropic. To do so, following the approach presented in [51], let us assume:
B ( t , r ) = R K ( t , r ) 1 ,
then integrating (10), we find:
K ( t , r ) 1 = 4 π q R d t + C ( r ) 2 ,
since in the nondissipative case, (55) becomes (44), then C ( r ) = k ( r ) 1 ) .
Thus, the line element reads:
d s 2 = d t 2 + ( R ) 2 d r 2 4 π q R d t + C ( r ) 2 + R 2 ( d θ 2 + sinh 2 θ d ϕ 2 ) .
Since we are considering dissipative systems, we need a transport equation. For simplicity, we adopt here the transport equation ensuing from the so-called “truncated” theory [52]; it reads:
τ h α β V γ q β ; γ + q α = κ h α β T , β ,
whose only nonvanishing independent component becomes:
τ q ˙ + q = κ B T ,
where κ and τ denote the thermal conductivity and the relaxation time, respectively.
In order to obtain specific models, we need to impose additional conditions. A first family of models is obtained from conditions on the complexity factor, while a second family is obtained by a specific restriction on the function B, particularly suitable for describing situations where a cavity surrounding the central region appears.

5.1. Models Obtained upon Conditions on the Complexity Factor

In our case, the complexity factor Y T F may be written as:
Y T F = R ¨ R R ¨ R + K ¨ 2 ( K 1 ) + K ˙ K 1 R ˙ R 3 4 K ˙ K 1 .
In order to obtain the models, we shall assume that the above structure scalar has the same form as in the nondissipative case, implying:
K ¨ 2 ( K 1 ) + K ˙ K 1 R ˙ R 3 4 K ˙ K 1 = 0 .
The integration of (61) produces:
R K ˙ ( K 1 ) 3 4 = C 1 ( r ) ,
where C 1 is an integration function. Integrating the above equation, we obtain:
K 1 = 4 [ C 1 2 d t ( R ) 2 + C 2 ( r ) ] 2 ,
where C 2 is another integration function.
Combining (56) and (63), it follows that:
4 π q R d t + C ( r ) = 2 C 1 2 d t ( R ) 2 + C 2 ( r ) ,
or:
C 2 ( r ) = 2 4 π q R d t + C ( r ) + C 1 2 d t ( R ) 2 .
Using (62), it follows that C 2 ( r ) = 0 ; thus:
2 4 π q R d t + C ( r ) + C 1 2 d t ( R ) 2 = 0 ,
implying that (64) may be written as:
2 π q = 1 R R C 1 ( r ) d t ( R ) 2 2 .
Let us first try to obtain models of dissipative dust satisfying (63).
Then, using (55) and (61), we obtain:
B ˙ B = R ˙ R K ˙ 2 ( K 1 ) B ¨ B = R ¨ R .
Feeding back (68) in (11) and (12) produces:
K = 2 R ¨ R + R ˙ 2 ,
K = 2 R R R ¨ R + R ¨ R + R ˙ R R ˙ R C 1 2 K 1 2 ( R ) 2 .
Next, taking the r derivative of (69) and replacing it in (), we find:
R ˙ R C 1 2 K 1 2 ( R ) 2 = 0 ,
from which it follows at once that there are no radiating dust solutions in this case. Therefore, in the following subsection, we relax the dust condition, and we consider models of radiating anisotropic fluids.

P = 0 , P r 0

Let us consider models with vanishing tangential pressure.
Replacing (55) in (12), the following expression for P is found:
8 π P = B ¨ B + R ¨ R + R ˙ R R ˙ R K ˙ 2 ( K 1 ) + K 2 R R .
In order to obtain a model, let us choose:
R ˙ R R ˙ R = K 2 R R ,
implying:
K 1 = R ˙ 2 .
Then, from the condition P = 0 , we obtain from (72):
B ¨ = 0 , B = b 1 ( r ) t + b 2 ( r ) ,
where b 1 and b 2 are two arbitrary functions.
Using (55), we find for R:
R B R ˙ = 0 , R = Φ a 1 ( r ) t + a 2 ( r ) ,
where Φ is an arbitrary function of its argument: and
a 1 ( r ) = e b 1 ( r ) d r , a 2 ( r ) = b 2 ( r ) e b 1 ( r ) d r d r .
Then, the physical variables read:
8 π μ = 1 R 2 R ˙ R K ˙ ( K 1 ) ,
4 π q = K ˙ 2 R K 1 ,
8 π P r = R ˙ R K ˙ K 1 + 1 R 2 .
To further specify the model, let us choose b 1 ( r ) and b 2 ( r ) as:
b 1 ( r ) = β 1 r + β 2 , b 2 ( r ) = ( r + β 2 ) α ,
implying:
a 1 ( r ) = ( r + β 2 ) β 1 , a 2 ( r ) = ( r + β 2 ) α + β 1 + 1 α + β 1 + 1 ,
and:
B = β 1 t r + β 2 + ( r + β 2 ) α ,
R = ( a 1 t + a 2 ) n ,
where β 1 , β 2 , α , and n are arbitrary constants.
Thus, the physical and kinematical variables for this model read:
8 π μ = 1 ( a 1 t + a 2 ) 2 n 2 ( n 1 ) n a 1 2 ( a 1 t + a 2 ) 2 ,
4 π q = ( n 1 ) n a 1 2 ( a 1 t + a 2 ) 2 ,
8 π P r = 1 ( a 1 t + a 2 ) 2 n 2 ( n 1 ) n a 1 2 ( a 1 t + a 2 ) 2 ,
m = ( a 1 t + a 2 ) n 2 ,
Θ = b 1 b 1 t + b 2 + 2 n a 1 a 1 t + a 2 ,
σ = b 1 b 1 t + b 2 n a 1 a 1 t + a 2 .
A simple calculation of the complexity factor ( Y T F ) for this model produces:
Y T F = ( n 1 ) n a 1 2 ( a 1 t + a 2 ) 2 .
It is worth noticing that it has exactly the same expression as q as given by (). Therefore, any solution of this family satisfying the vanishing complexity factor is necessarily nondissipative. On the other hand, Y T F is zero if n = 0 and/or a 1 = 0 and/or n = 1 . The first two conditions are ruled out at once from (82) and (85). Thus, the solution of this family with a vanishing complexity factor is characterized by n = 1 , which using (85) and (), produces:
P r = μ .
The above is the stiff equation of state originally considered by Zeldovich (see [48]).

5.2. Models with B = 1

We next assume B = 1 in order to obtain some additional analytical models. As discussed in [53], such a condition is particularly suitable for describing fluid distributions whose center is surrounded by an empty cavity, a scenario we expect for the kind of fluid distributions we are dealing with in this work.
The corresponding Einstein equations may be written as:
8 π μ = 1 R 2 2 R R R R 2 + R ˙ 2 R 2 ,
4 π q = R ˙ R ,
8 π P r = 1 R 2 + R R 2 R ˙ R 2 + 2 R ¨ R ,
8 π P = R R R ¨ R .

5.2.1. Nondissipative Case

Let us first consider the nondissipative case ( q = 0 ). In this case, it follows at once from (94) that R is a separable function, i.e., it takes the form:
R = R 1 ( t ) + R 2 ( r ) ,
where R 1 and R 2 are arbitrary functions of their arguments.
Using (97) in (44), it follows at once that:
R 2 = k ( r ) 1 .
In order to exhibit an exact solution, let us further assume P = 0 . Using this condition in (96) produces:
R 1 ( t ) = a t 2 + b 1 t + c 1 , R 2 ( r ) = a r 2 + b 2 r + c 2 ,
where a , b 1 , c 1 , b 2 , c 2 are arbitrary constants.
The physical and kinematical variables for this model are:
8 π μ = 1 α 2 1 4 a α β 2 + γ 2 ,
8 π P r = 1 α 2 1 4 a α + β 2 γ 2 ,
P r + μ = a π α ,
Θ = 2 γ α ,
σ = γ α ,
m = α 2 β 2 γ 2 + 1 ,
where:
α a ( t 2 + r 2 ) + b 1 t + b 2 r + c 1 + c 2 ; β 2 a r + b 2 , γ 2 a t + b 1 .
For this model, the expression for Y T F reads:
Y T F = 2 a α .
Therefore, the vanishing complexity factor implies a = 0 , which using (100) and (101) produces:
P r = μ .
Thus, the solution of this family with the vanishing complexity factor condition is also characterized by the stiff equation of state.

5.2.2. Dissipative Case

Let us now consider the dissipative case ( q 0 ). If we impose the condition P = 0 , then we obtain the equation R ¨ = R , whose general solution is of the form:
R ( t , r ) = c 1 Ψ ( t + r ) + c 2 Φ ( t r ) ,
where c 1 , c 2 are arbitrary constants and Ψ , Φ arbitrary functions of their arguments.
As an example, let us choose:
R ( t , r ) = c sin a ( t r ) ,
where a , c are constants. Then, for the kinematical and physical variables, we obtain:
8 π μ = 2 a 2 1 c 2 sin 2 [ a ( t r ) ] ,
4 π q = a 2 ,
8 π P r = 2 a 2 + 1 c 2 sin 2 [ a ( t r ) ] ,
Θ = 2 a cot [ a ( t r ) ] ,
σ = a cot [ a ( t r ) ] ,
m = c sin [ a ( t r ) ] 2 ,
For this case, the temperature T ( t , r ) , calculated from (59), reads:
T ( t , r ) = a 2 r 4 π κ + T 0 ( t ) ,
whereas the expression for Y T F is:
Y T F = a 2 ,
implying because of (110) that no solution of this family has a vanishing complexity factor.
Finally, as an alternative model, we may assume:
R = a ( t r ) n ,
where a , n are constants.
The ensuing physical and kinematical variables are in this case:
8 π μ = 1 a 2 t r 2 n 2 n n 1 t r 2 ,
4 π q = n n 1 t r 2 ,
8 π P r = 1 a 2 t r 2 n 2 n n 1 t r 2 ,
Θ = 2 n t r ,
σ = n t r ,
m = a t r n 2 .
Using (59), the expression for the temperature becomes:
T ( t , r ) = n ( n 1 ) 4 π κ ( t r ) n ( n 1 ) τ 4 π κ ( t r ) 2 + T 0 ( t ) .
The complexity factor for this family of solutions reads:
Y T F = n ( n 1 ) ( t r ) 2 .
The above scalar may vanish only if n = 0 and/or n = 1 . The first possibility has to be ruled out from a simple inspection of (120), and therefore, the vanishing complexity factor conditions require n = 1 , implying because of (121) that the fluid is nondissipative and because of (120) and (122) that the fluid satisfies the stiff equation of state P r = μ .

6. Discussion and Conclusions

We investigated in some detail all possible solutions of fluids endowed with the hyperbolic symmetry (3), characterized by nonvanishing shear, inhomogeneous energy–density, and vanishing four-acceleration (geodesics). So defined, these solutions are entitled to be considered as hyperbolic symmetric versions of LTB spacetimes.
The first class of solution corresponds to nondissipative dust configurations. Comparing with the spherically symmetric case, we observe that only one family of solutions ( k ( r ) > 0 ) exists, instead of the three families existing in this latter case ( k ( r ) 0 ).
These solutions cannot satisfy the vanishing complexity factor, neither can they evolve in the quasi-homologous regime. On the other hand, the scalar expansion is positive as expected for pure dust submitted to a repulsive gravity.
Next, we analyzed the case of dissipative anisotropic fluids. To do this, we generalized the expression (44) by assuming (55). Different specific models were found from two different conditions. One class of solutions was obtained from a condition imposed on the complexity factor (61). It was shown that in this case, the pressure must be anisotropic. A solution of this type was found assuming further that P = 0 . The subclass of this solution satisfying the vanishing complexity factor is necessarily nondissipative and satisfies the stiff equation of state P r = μ .
The other class of solutions was found under the condition B = 1 . For the nondissipative case, a family of solutions was found under the additional condition P = 0 . In this case as well, the vanishing complexity factor condition implies the stiff equation of state P r = μ . In the dissipative case, two families of solutions were found from different assumptions on the specific form of R. Thus, assuming (110), we found a solution never satisfying the vanishing complexity factor condition, whereas assuming (119), such a condition can be satisfied, implying that the fluid is nondissipative and satisfies the stiff equation of state P r = μ . It is worth noticing that the temperature for the first of the above solutions (117) does not contain terms depending on the relaxation time. In other words, the model behaves as if the fluid is always in the thermal stationary state, a result that becomes intelligible when we observe that the dissipative flux (112) is constant. Instead, for the second family of solutions, the temperature (126) clearly exhibits the effects of transient phenomena (i.e., those depending on τ ).
Finally, we would like to conclude with a general comment: all the models exhibited above were found with the sole purpose of illustrating the richness of solutions endowed with hyperbolic symmetry and sharing the general physical and geometrical properties (excluding the isometry group) characterizing the LTB spacetimes. It is now up to cosmologists and astrophysicists to decide if any of the above models (or any other HSLTB not described in this manuscript) could be of any use in the study of specific scenarios, as for example cosmological models beyond the standard FRW solution [54,55].

Author Contributions

All authors contributed equally to this work. Conceptualization, L.H.; methodology, L.H., A.D.P., and J.O.; software, J.O.; formal analysis, L.H., A.D.P., and J.O.; writing—original draft preparation, L.H. writing—review and editing, L.H., A.D.P., and J.O.; funding acquisition, L.H., J.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministerio de Ciencia, Innovación y Universidades, Grant Number PGC2018-096038-B-I00, and Junta de Castilla y León, Grant Number SA096P20.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; nor in the decision to publish the results.

References

  1. Herrera, L.; Di Prisco, A.; Ospino, J. Dynamics of hyperbolic symmetric fluids. Symmetry 2021, 13, 1568. [Google Scholar] [CrossRef]
  2. Herrera, L.; Di Prisco, A.; Ospino, J. Hyperbolically symmetric static fluids: A general study. Phys. Rev. D 2021, 103, 024037. [Google Scholar] [CrossRef]
  3. Herrera, L.; Witten, L. An alternative approach to the static spherically symmetric vacuum global solutions to the Einstein’s equations. Adv. High Energy Phys. 2018, 2018, 3839103. [Google Scholar]
  4. Herrera, L.; Di Prisco, A.; Ospino, J.; Witten, L. Geodesics of the hyperbolic symmetric black hole. Phys. Rev. D 2020, 101, 064071. [Google Scholar] [CrossRef] [Green Version]
  5. Rosen, N. The nature of the Schwarzschild singularity. In Relativity; Carmeli, M., Fickler, S.I., Witten, L., Eds.; Plenum Press: New York, NY, USA, 1970; pp. 229–258. [Google Scholar]
  6. Rindler, W. Relativity. Special, General and Cosmological; Oxford University Press: New York, NY, USA, 2001; pp. 260–261. [Google Scholar]
  7. Caroll, S. Spacetime and Geometry. An Introduction to General Relativity; Addison Wesley: San Francisco, CA, USA, 2004; p. 246. [Google Scholar]
  8. Harrison, B.K. Exact Three–Variable Solutions of the Field Equations of General Relativity. Phys. Rev. 1959, 116, 1285. [Google Scholar] [CrossRef]
  9. Ellis, G. Dynamics of Pressure–Free Matter in General Relativity. J. Math. Phys. 1967, 8, 1171. [Google Scholar] [CrossRef]
  10. Stephani, H.; Kramer, D.; MacCallum, M.; Honselaers, C.; Herlt, E. Exact Solutions to Einsteins Field Equations, 2nd ed.; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  11. Gaudin, M.; Gorini, V.; Kamenshchik, A.; Moschella, U.; Pasquier, V. Gravity of a static massless scalar field and a limiting Schwarzschild-like geometry. Int. J. Mod. Phys. D 2006, 15, 1387–1399. [Google Scholar] [CrossRef] [Green Version]
  12. Rizzi, L.; Cacciatori, S.L.; Gorini, V.; Kamenshchik, A.; Piattella, O.F. Dark matter effects in vacuum spacetime. Phys. Rev. D 2010, 82, 027301. [Google Scholar] [CrossRef] [Green Version]
  13. Lobo, F.S.N.; Mimoso, J.P. Possibility of hyperbolic tunneling. Phys. Rev. D 2010, 82, 044034. [Google Scholar] [CrossRef] [Green Version]
  14. Kamenshchik, A.Y.; Pozdeeva, E.O.; Starobinsky, A.A.; Tronconi, A.; Vardanyan, T.; Venturi, G.; Verno, S.Y. Duality between static spherically or hyperbolic symmetric solutions and cosmological solutions in scalar-tensor gravity. Phys. Rev. D 2018, 98, 124028. [Google Scholar] [CrossRef] [Green Version]
  15. Madler, T. On the affine-null metric formulation of General Relativity. Phys. Rev. D 2019, 99, 104048. [Google Scholar] [CrossRef] [Green Version]
  16. Maciel, A.; Delliou, M.L.; Mimoso, J.P. New perspectives on the TOV equilibrium from a dual null approach. Class. Quantum Gravity 2020, 37, 125005. [Google Scholar] [CrossRef]
  17. Lemaître, G. L’Univers en expansion. Ann. Soc. Sci. Bruxelles 1933, A 53, 51. [Google Scholar]
  18. Tolman, R.C. Effect of Inhomogeneity on Cosmological Models. Proc. Natl. Acad. Sci. USA 1934, 20, 169–176. [Google Scholar] [CrossRef] [Green Version]
  19. Bondi, H. Spherically symmetrical models in general relativity. Mon. Not. R. Astron. Soc. 1947, 107, 410–425. [Google Scholar] [CrossRef] [Green Version]
  20. Krasinski, A. Inhomogeneous Cosmological Models; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
  21. Plebanski, J.; Krasinski, A. An Introduction to General Relativity and Gravitation; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
  22. Sussman, R. Quasilocal variables in spherical symmetry: Numerical applications to dark matter and dark energy sources. Phys. Rev. D 2009, 79, 025009. [Google Scholar] [CrossRef] [Green Version]
  23. Sussman, R. Radial asymptotics of Lemaître-Tolman-Bondi dust models. Gen. Relativ. Gravit. 2010, 42, 2813–2864. [Google Scholar] [CrossRef] [Green Version]
  24. Matravers, D.R.; Humphreys, N.P. Matching Spherical Dust Solutions to Construct Cosmological Models. Gen. Rel. Grav. 2001, 33, 531–552. [Google Scholar] [CrossRef]
  25. Hellaby, C.; Krasinski, A. Alternative methods of describing structure formation in the Lemaitre-Tolman model. Phys. Rev. D 2006, 73, 023518. [Google Scholar] [CrossRef] [Green Version]
  26. Eardley, D.M.; Smarr, L. Time functions in numerical relativity: Marginally bound dust collapse. Phys. Rev. D 1979, 19, 2239. [Google Scholar] [CrossRef]
  27. Waugh, B.; Lake, K. Shell-focusing singularities in spherically symmetric self-similar spacetimes. Phys. Rev. D 1989, 40, 2137. [Google Scholar] [CrossRef]
  28. Joshi, P.S.; Dwivedi, I.H. Naked singularities in spherically symmetric inhomogeneous Tolman-Bondi dust cloud collapse. Phys. Rev. D 1993, 47, 5357. [Google Scholar] [CrossRef] [Green Version]
  29. Joshi, P.S.; Singh, T.P. Role of initial data in the gravitational collapse of inhomogeneous dust. Phys. Rev. D 1995, 51, 6778. [Google Scholar] [CrossRef] [Green Version]
  30. Mimoso, J.; Le Delliou, M.; Mena, F. Separating expansion from contraction in spherically symmetric models with a perfect fluid: Generalization of the Tolman-Oppenheimer-Volkoff condition and application to models with a cosmological constant. Phys. Rev. D 2010, 81, 123514. [Google Scholar] [CrossRef] [Green Version]
  31. Le Delliou, M.; Mena, F.; Mimoso, J. Separating expansion from contraction: Generalized TOV condition, LTB models with pressure and ΛCDM. AIP Conf. Proc. 2010, 1241, 1011–1015. [Google Scholar]
  32. Vaz, C.; Witten, L.; Singh, T.P. Toward a midisuperspace quantization of Lemaître-Tolman-Bondi collapse models. Phys. Rev. D 2001, 63, 104020. [Google Scholar] [CrossRef] [Green Version]
  33. Bojowald, M.; Harada, T.; Tibrewala, R. Lemaitre-Tolman-Bondi collapse from the perspective of loop quantum gravity. Phys. Rev. D 2008, 78, 064057. [Google Scholar] [CrossRef] [Green Version]
  34. Coley, A.A.; Pelavas, N. Averaging spherically symmetric spacetimes in general relativity. Phys. Rev. D 2006, 74, 087301. [Google Scholar] [CrossRef] [Green Version]
  35. Coley, A.A.; Pelavas, N. Averaging in spherically symmetric cosmology. Phys. Rev. D 2007, 75, 043506. [Google Scholar] [CrossRef] [Green Version]
  36. Viaggiu, S. A Lemaitre-Tolman-Friedmann Universe without Dark Energy. arXiv 2009, arXiv:0907.0600v1. Available online: https://arxiv.org/abs/0907.0600v1 (accessed on 3 July 2009).
  37. Sussman, R. Back-reaction and effective acceleration in generic LTB dust models. Class. Quantum Gravity 2011, 28, 235002. [Google Scholar] [CrossRef]
  38. Krasinski, A.; Hellaby, C.; Bolejko, K.; Celerier, M.N. Imitating accelerated expansion of the Universe by matter inhomogeneities: Corrections of some misunderstandings. Gen. Rel. Gravit. 2010, 42, 2453–2475. [Google Scholar] [CrossRef] [Green Version]
  39. Celerier, M.N. Some clarifications about Lemaître–Tolman models of the Universe used to deal with the dark energy problem. Astr. Astrophys. 2012, 543, A71. [Google Scholar] [CrossRef] [Green Version]
  40. Misner, C.; Sharp, D. Relativistic Equations for Adiabatic, Spherically Symmetric Gravitational Collapse. Phys. Rev. 1964, 136, B571. [Google Scholar] [CrossRef]
  41. Cahill, M.; McVittie, G. Spherical Symmetry and Mass–Energy in General Relativity. I. General Theory. J. Math. Phys. 1970, 11, 1382. [Google Scholar] [CrossRef]
  42. Chan, R. Radiating gravitational collapse with shear viscosity. Mon. Not. R. Astron. Soc. 2000, 316, 588–604. [Google Scholar] [CrossRef]
  43. Herrera, L.; Ospino, J.; Di Prisco, A.; Fuenmayor, E.; Troconis, O. Structure and evolution of self-gravitating objects and the orthogonal splitting of the Riemann tensor. Phys. Rev. D 2009, 79, 064025. [Google Scholar] [CrossRef] [Green Version]
  44. Bel, L. Sur la radiation gravitationelle. C. R. Acad. Sci. 1958, 247, 1094–1096. [Google Scholar]
  45. Bel, L. Radiation states and the problem of energy in general relativity. Gen. Relativ. Gravit. 2000, 32, 2047–2078. [Google Scholar] [CrossRef]
  46. Bel, L. Introduction d’un tenseur du quatrieme order. C. R. Acad. Sci. Paris 1959, 248, 1297–1300. [Google Scholar]
  47. García-Parrado Gómez Lobo, A. Dynamical laws of superenergy in general relativity. Class. Quantum Gravity 2008, 25, 015006. [Google Scholar] [CrossRef]
  48. Herrera, L. New definition of complexity for self-gravitating fluid distributions: The spherically symmetric static case. Phys. Rev. D 2018, 97, 044010. [Google Scholar] [CrossRef] [Green Version]
  49. Herrera, L.; Di Prisco, A.; Ospino, J. Definition of complexity for dynamical spherically symmetric dissipative self-gravitating fluid distributions. Phys. Rev. D 2018, 98, 104059. [Google Scholar] [CrossRef] [Green Version]
  50. Herrera, L.; Di Prisco, A.; Ospino, J. Quasi–homologous evolution of self–gravitating systems with vanishing complexity factor. Eur. Phys. J. C 2020, 80, 631. [Google Scholar] [CrossRef]
  51. Herrera, L.; Di Prisco, A.; Ospino, J.; Carot, J. Lemaitre–Tolman–Bondi dust spacetimes: Symmetry properties and some extensions to the dissipative case. Phys. Rev. D 2010, 82, 024021. [Google Scholar] [CrossRef] [Green Version]
  52. Triginer, J.; Pavón, D. Heat transport in an inhomogeneous spherically symmetric universe. Class. Quantum Gravity 1995, 12, 689–698. [Google Scholar] [CrossRef]
  53. Herrera, L.; Le Denmat, G.; Santos, N.O. Cavity evolution in relativistic self–gravitating fluids. Class. Quantum Gravity 2010, 27, 135017. [Google Scholar] [CrossRef] [Green Version]
  54. Kopteva, E.; Bormotova, I.; Churilova, M.; Stuchlik, Z. Accelerated Expansion of the Universe in the Model with Nonuniform Pressure. Astrophys. J. 2019, 887, 98. [Google Scholar] [CrossRef]
  55. Cadoni, M.; Sanna, A.P. Emergence of a cosmological constant in anisotropic fluid cosmology. Int. J. Mod. Phys. A 2021, 36, 2150156. [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

Herrera, L.; Di Prisco, A.; Ospino, J. Hyperbolically Symmetric Versions of Lemaitre–Tolman–Bondi Spacetimes. Entropy 2021, 23, 1219. https://doi.org/10.3390/e23091219

AMA Style

Herrera L, Di Prisco A, Ospino J. Hyperbolically Symmetric Versions of Lemaitre–Tolman–Bondi Spacetimes. Entropy. 2021; 23(9):1219. https://doi.org/10.3390/e23091219

Chicago/Turabian Style

Herrera, Luis, Alicia Di Prisco, and Justo Ospino. 2021. "Hyperbolically Symmetric Versions of Lemaitre–Tolman–Bondi Spacetimes" Entropy 23, no. 9: 1219. https://doi.org/10.3390/e23091219

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