A Heuristic Approach to the Far-Future State of a Universe Dominated by Phantom Energy ()
1. Introduction
The analysis of this paper starts with a brief review of the standard arguments that lead to the introduction of dark energy. The framework is a homogeneous, isotropic universe, which is spatially flat (
). Its spacetime can be described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric1
(1)
with a(t) the scale factor. Solving the Einstein equations
for this type of metric leads to the well-known Friedmann equations
(2)
(3)
with
the energy density and P the pressure of the respective component. P and
follow an equation of state
(4)
with w the equation of state parameter. Combining (2) and (3) results in
(5)
Simultaneously, the relation for the energy-momentum conservation
leads to the continuity equation
(6)
The observation of the accelerating expansion of the universe has been one of the most remarkable discoveries in Astrophysics and Cosmology [1] [2]. This expansion cannot be explained by the current forms of matter and radiation known so far. Looking at Equation (3), a form of energy with negative pressure and
should exist, so that
. This substance, called dark energy, has the opposite effect (anti-gravitational) to that of gravity, accounts for almost 3/4 of the energy content of the universe and despite intense research since its discovery, it is yet unknown what this form of energy is [3].
One interesting case to consider is that of dark energy as a fluid with equation of state
. Then the dark energy is called phantom energy. In this case, Equations (5) and (6) lead to
and
, i.e., the Hubble rate and the energy density increase with time. These models violate the so-called Null Energy Condition (NEC), which states that for any light-like vector
, with
, the following condition for the energy-momentum tensor should hold
. In the case of a spatially-flat FLRW universe, this equivalently leads to the fact that
, which is violated by a fluid with equation of state with
.
On the one hand, it is exactly the violation of NEC that has made a large part of the academic community not consider the models of dark energy with equation of state parameter
viable. This is mainly based on the grounds that the violation of NEC is proven to lead to instabilities (in the form of “ghosts” and “tachyons”) for a large class of models [4].
On the other hand, it has been proven that these instabilities can successfully be controlled in the context of effective field theories, despite the violation of NEC [5]. In this work, it will be assumed that these instabilities are under control. Additionally, although there is the dominant belief in the scientific community that the dark energy ought to be the energy of the vacuum in the form of a cosmological constant
or equivalently a substance with
, it is still not clear from the data whether the dark energy fluid has w greater, less or equal to −1 [6] [7]. Only future experiments have the potential to distinguish
from percent-level deviations.
Therefore, the idea that dark energy could be phantom energy is still an open possibility that could be discovered in the future. It is definitely worth exploring the implications of this idea into more depth. Over the last years, there has been extensive analysis in this direction. The most interesting realisation is that if dark energy is indeed phantom energy, then it could lead to the so-called “Big Rip” scenario, where eventually every part of the universe could be “ripped” apart in a finite amount of time due to the super-exponential expansion of the underlying space2 [8] [9] [10]. However, infinities usually indicate an incompleteness of the respective theory in some specific limits. Therefore, a reasonable question to ask is whether the Big Rip singularity can be avoided in any way possible in the far-future of the universe.
2. Far-Future State of the Universe
2.1. Gravitational Schwinger Pair-Production in the Far-Future of the Universe
To answer the question posed in the previous section, one has to think of what could possibly act as a counterpart to the super-exponential expansion of the universe caused by the phantom energy3. In this work, the attention is turned to a phenomenon called the Schwinger effect. This phenomenon, first derived by Julian Schwinger in 1951, is the production of a particle-antiparticle pair out of the vacuum in the presence of a strong electric field [28]. The spectrum of the produced pairs is given by the formula [29] [30]
(7)
with
the transverse momentum,
the time component of the spacetime 4-vector,
, m the mass, qe the charge of each of the produced particles of the pair and E the electric field.
It is easy to observe that one should create a very strong electric field, i.e.
, in order for the exponential factor in Equation (7) not to suppress the pair production from the vacuum. Such large electric fields are difficult to be produced in the laboratories which is why the Schwinger effect has not yet been observed in Nature. However, one may still try to think how this effect could potentially play a role on a cosmological scale. To make this connection, the horizon temperature
for the case of de Sitter spacetime is defined [31] [32] [33].
(8)
where
is the acceleration of gravity at the surface of de Sitter horizon and is proportional to the Hubble rate H in the case of de Sitter spacetime. Connecting the relation (8) with the Unruh temperature of a uniformly accelerated observer [34], one can find the following relation between the gravitational acceleration
of a particle of mass m, charge qe and transverse momentum
with the electric field E [35] [36] 4
(9)
The combination of the Equations (8) and (9) results in
(10)
One may calculate the contribution of the gravitational Schwinger effect to the energy density. According to the definition of the energy density, this is as following
(11)
Then using the relations (7) and (10), (11) becomes
(12)
Assuming light-produced particles (
) (12) becomes
(13)
with
the respective coefficient. This is a radiation term, which is expected given the light-produced particles and whose energy density scales like
according to the Stefan-Boltzmann Law. However, in the case of the gravitational production of particles, the coefficient does not have to be the same as in the Stefan-Boltzmann Law.5
Considering the contributions of the phantom energy and the Schwinger effect, the Friedmann Equation (2) becomes
(14)
Simultaneously, considering (13) and that
, the continuity Equation (6) becomes
(15)
where it is assumed that
is a constant. Also
is a constant function of
and
6.
The analysis above has been done in the context of an almost de Sitter space, although
. This is a reasonable assumption because, even if
, still it should be close to −1, according to the observations, as discussed in paragraph (1). Additionally, there is no concrete and universally accepted definition of the surface gravity for general curved spacetimes [37], so (8) could not be used for the general case of
. Also for our model we assume that there no catastrophic instabilities, which can be justified in the context of [5], as mentioned in the paragraph (1).
It is noted that the phantom fluid does not evolve independently from the Schwinger pair production, so the continuity Equation (6) does not hold independently for each component. One can see that the Schwinger effect contributes with a term proportional to H4 and H5 in (14) and (15) respectively. These terms become important to higher values of H, therefore later in the evolution of the universe if it is dominated by a phantom substance. Even though at later times the solutions of (14) and (15) start deviating from the de Sitter case, one expects gravitational particle production to occur whenever there is an event horizon (as with Hawking radiation in the case of backholes [31]). Therefore, one would expect at least a radiation term proportional to H4 and a source term proportional to H5 in the Friedmann and continuity equations as in (14) and (15), independent of whether the spacetime is described by a metric close to de Sitter or not.
The goal now is to solve the Equations (14) and (15) simultaneously. First of all, the Equation (14) has two solutions in H2
(16)
In (16), only the solution with the relative minus sign is kept, since this is the one that reduces to
for small H, where the H4 term coming from the Schwinger effect in (14) is negligible7. Additionally, this solution has an upper value for
, which is
. This leads to a maximum value of H2, which is
8. Setting also
, the solution with the relative minus sign in Equation (16) becomes
(17)
Using the Equation (17), the continuity Equation (15) becomes
(18)
Integrating the Equation (18) and using the fact that
results in
(19)
Putting (19) into (16) leads to
(20)
From Figure 1 and Figure 2 one can see that the energy density and Hubble rate do not reach infinite values in a finite amount of time, as it would happen
Figure 1. The plot of
as a function of
.
as
.
Figure 2. The plot of
as a function of
.
as
.
in the case of a Big Rip scenario. On the contrary, they reach the maximum values
and
respectively as
and therefore the Big Rip singularity is avoided. This happens thanks to the ever increasing rate of Schwinger pair-production that an observer inside a causal horizon of radius
would observe as the Hubble rate increases.
2.2. Horizon Entropy in the Far-Future of the Universe
It is also important to have a qualitative understanding of the evolution of the horizon entropy in the scenario discussed in this paper. Its calculation is quite straightforward. Using the analogy between the thermodynamics of black holes and cosmological horizons, the Generalised Second Law (GSL) for black holes [38] [39] [40] extended to de Sitter horizons is [32] [33]
(21)
(22)
where
is the horizon entropy,
is the entropy outside the horizon,
is the area of the de Sitter horizon and
is the Hubble radius9.
The relations (21) and (22) have been proven to hold true also in the case of accelerated horizons and are independent of whether the horizon area
increases or decreases, as long as the rate of increase of the entropy outside the horizon
outweighs the rate of decrease of the horizon entropy
[41].
In the case analysed in this paper, where the universe is dominated by the phantom fluid, the horizon entropy decreases because of the increase of the Hubble rate, until it reaches asympotically a minimum constant value
as
which is
(23)
One observes that the entropy does not reach a zero value thanks to the particle-pair production because of the Schwinger effect. When
reaches its minimum value
in the far-future, the universe enters into a de Sitter phase of constant and high Hubble rate
and energy density
. The fact that the horizon entropy decreases is not a problem as long as the causal patch of an observer is not an isolated system and the rate of increase of the entropy outside the horizon outweighs the rate of decrease of the horizon entropy, as discussed above10.
Using (17) and (23) results in the following evolution of the horizon entropy (22) (Figure 3)
(24)
Figure 3. The plot of
as a function of
.
as
.
Figure 4. The plots of the phantom energy density
and the radiation/matter energy density
as a function of time t in two cycles, each of period T, according to the model described in this work. One can see how in each period the phantom energy density eventually dominates over radiation/matter. It never reaches infinity, but rather a constant and high value
, thanks to the backreaction from the gravitational Schwinger pair-production. Then the universe passes smoothly to an inflationary phase with
until its decay refills the universe with radiation and matter. Therefore a continuous cycle can be assumed where superacceleration, followed by inflation, followed by reheating, followed by matter/radiation domination, followed by phantom substance domination takes place periodically in the history of the universe. In this graph, the same scaling for matter and radiation is assumed for simplicity. In general, each cycle can have a different value of period T, which can be defined by the specific characteristics of each cycle of the universe.
2.3. Cyclic/Periodic Model of the Universe?
The combination of the results of paragraphs (2.1) and (2.2) leads to the following conclusion: If the universe is dominated by a phantom fluid with
and if the gravitational Schwinger pair-production takes place, the Big Rip singularity could be avoided in the far-future. In particular, because the
and
are solutions to the Equations (14) and (15), the energy density and the Hubble rate would remain constant after they reach these values in the far-future, according to the analysis in this paper, and the universe could start a phase of de Sitter inflation. Then, by taking one more step of speculation, the phantom substance could decay through an unknown hypothetical mechanism, which would lead to the reheating of the universe and to a new Big Bang; thus the beginning of a new cycle of the universe. Of course, the specific characteristics of the decay of the phantom fluid to matter and radiation would have to be understood and this is beyond the scope of this paper11. This idea can be seen graphically in Figure 4.
3. Conclusion
In this work, the behaviour of a universe dominated by phantom energy with a generic equation of state parameter
is analysed. Working solely in the context of general relativity and assuming no instabilities caused by the violation of NEC, it is found that the Big Rip singularity could potentially be avoided in the far-future because of the gravitational Schwinger pair-production. The universe would reach a high but constant Hubble rate
and energy density
, passing to a de Sitter inflationary phase in a finite amount of time
. As a final step, it is assumed that the phantom substance could decay to matter/radiation and reheat the universe until the first would dominate again, leading thus to a cyclic/periodic model of the universe. This situation could in theory be repeated an infinite amount of times, unless some other process stops it or changes it. This cyclic/periodic model of the universe is highly speculative and the specific dynamics of the phantom substance decay to radiation and matter are left to be investigated in future works. Even more importantly, it should be clarified from future observations whether dark energy is a phantom fluid, a cosmological constant or something completely different and unexpected.
Acknowledgements
The author would like to thank Hitoshi Murayama and Yasunori Nomura for useful discussions, and especially Simone Ferraro for the useful guidance, valuable inputs and his support throughout this work.
NOTES
1In this paper, the author works in Planck units
.
2In particular, the energy density and thus the Hubble rate reach infinite values in a finite amount of time.
3In this paper the author solely works in the context of general relativity, without assuming f(R) theories of modified gravity (for a review on this topic, see [11] [12]). In this context, there have been various analysis of how the Big Rip or other type of singularities could be avoided in the far-future. Here, some of the analyses tackling of these topics are mentioned [13] - [27].
4Given that one may equate the gravitational acceleration kh with the acceleration of the particle.
5From now on, the gravitational Schwinger effect is simply referred to as Schwinger effect.
6For convenience, one may set
.
7The other solution with the plus sign reduces to
in the same limit.
8
as
, which is the Big Rip case if the Schwinger particle production would be negligible. Because this analysis is done only in the context of general relativity, without considering modifications of gravity in the UV regime, one should require that
, and
, which equivalently means that
.
9“Outside” for the case of the black holes is “inside” for the case of de Sitter space, and in general for the case of cosmological horizons.
10An idea has been recently introduced: the central dogma about cosmological horizons [42]. This idea is an extension of the central dogma about black holes [43] to cosmological horizons and considers that every causal patch is supposed to be an isolated system. If this conjecture would hold true, the presence of a phantom substance would violate the second law of thermodynamics, which is one of the most sacred laws of Physics, therefore the whole discussion in this paper would have to be abandoned. However, this dogma is based on holographic arguments, which silently imply the NEC, which by definition is violated in the case of a phantom fluid. This discussion is also done in [44], but for bouncing models. Therefore, the whole discussion in this paper would be ruled out just by bias and not by any independent arguments.
11An idea similar to the time-periodic (cyclic) universe, discussed in this work, has already been suggested in [5], but there in the context of effective field theories. However, in this paper the author uses additional entropic arguments and a mechanism of stopping smoothly the ever increasing energy density and Hubble rate through the gravitational Schwinger effect, before the decay of the phantom fluid to radiation/matter in the next cycle.