1. Introduction
The accelerated cosmic expansion indicates the presence of a negative-pressure component to the total energy density of the universe today. This component is the source of the accelerated expansion and, either material fluid or geometry, is known as dark energy. The observational evidence and theoretical consistency of dark energy is well supported; for a review see Refs. [1-3] and references therein. Dark energy nowadays is a fundamental element of the standard cosmological model, however the physical mechanism behind it still a mystery. The simplest and best-known model of dark energy is the energy of the vacuum represented as the cosmological constant added to the Einstein field equations. The main feature of vacuum energy is that its energy density is constant in time and is spatially smooth. The cosmological constant with cold dark matter are the key elements of the standard cosmological model or ΛCDM.
In cosmology, a perfect fluid description is adequate to model the known cosmic material components (e.g. photons, baryons, neutrinos) and also dark matter (even though we don’t have a consolidated microscopic theory of dark matter, when modeled as dust, a pressureless perfect fluid, there is theoretical and observational consistency). Dark energy, like the other main cosmological components, as a first approach is modeled in the framework of perfect fluids in a homogeneous and isotropic expanding universe. The perfect fluid approach applied to describe the cosmic components fit very well the cosmological observations at background and linear perturbation level; and also for some components has an underlying microscopic theory that support it.
To obtain an accelerated cosmic expansion the material content of the universe must violate the strong energy condition, ∑ i (ρ i + 3p i ) < 0, where ρ i and p i are the energy density and the pressure of each component respectively. The total pressure must be negative. Since baryonic and dark matter are pressureless and radiation pressure is ρ r /3 > 0, then there must be a source of negative pressure, this source is the dark energy. In the standard approach, the dark energy is considered as a fluid with negative pressure and barotropic equation of state, p = wρ (with w constant), and value w = −1. The dark energy equation of state p = −ρ implies that the energy density of dark energy is a constant, i.e., has no dynamics.
Recent observational reconstructions of the dark energy equation of state shows that the parameter w could be dynamical, i.e., depends on time (scale factor) [4, 5] allowing to cross the phantom divide line (w < −1), and there are even studies that suggest a slight preference to the cosmological phantom regime w < −1 [6-8], but it is not clear if it corresponds properly to a phantom fluid or to an additional physical effect which gives an effective phantom dark energy equation of state.
Dynamical dark energy, w = w(α), together with the assumption of dark energy as a perfect fluid brings up some thermodynamic problems such as the positiveness of the entropy, temperature, and chemical potential implies that w ≥ −1, which is in direct conflict with phantom dark energy [9,10].
One way to avoid some thermodynamic problems is to suppose that the dark energy is a fluid with bulk viscosity breaking the perfect fluid hypothesis. A fluid with bulk viscosity means that dissipative processes occur, which in cosmic fluids has the advantage of allowing the violation of the dominant energy condition (ρ + p < 0) without the dark energy necessarily being phantom [11]. In this case we have an effective pressure given by p eff = p + Π, where p = wρ is the barotropic pressure and Π < 0 is the viscous pressure. The possibility of explaining the accelerated expansion of the universe at the late future as an effect of the effective negative pressure due to bulk viscosity in the cosmic fluids was first considered in [12,13]. Viscous matter or radiation cosmologies can be mapped into the phantom dark energy scenario [14,15] and also a viscous fluid is able to produce a Little Rip cosmology as a purely viscous effect [16].
A perfect fluid in equilibrium generates no entropy and no frictional heating because its dynamic is reversible and without dissipation. We know that real fluids behave irreversibly, and if we want to model dissipative processes we require a relativistic theory of dissipative fluids. A classical irreversible thermodynamics was first extended from Newtonian to relativistic fluids by Eckart in 1940 [17]. In the Eckart theory, the effective pressure of the cosmic fluid is modeled as Π = −3ξH, where ξ is a function and H the Hubble parameter. The Eckart theory has the problem that dissipative perturbations propagate at infinite speeds; this noncausal feature is its main limitation and therefore this approach could be useful to find insight from toy models but not as a realistic theory. Nevertheless, Eckart theory has been used widely to model bulk viscosity in dark matter and dark energy models, for example, interacting viscous dark matter and dark energy [18, 19], the possibility of crossing the phantom divide line [20], the magnitude of the viscosity to achieve this crossing using cosmological data [21], Big Rip singularities for various forms of the equation of state parameter and the bulk viscosity [22], and unified dark fluid cosmologies [23-26]. The causal extension of the Eckart theory is the so-called Israel-Stewart theory [27,28]. This approach presents a better description than Eckart theory, including a casual description of dissipative processes associated to small deviations of equilibrium. The Israel-Stewart theory converges to the Eckart theory, when the collision time-scale in the transport equation of fluid is zero, i.e., when the bulk viscous model is noncausal and unstable. An interesting review about viscous cosmology can be found in [29] and for a dynamical analysis for a bulk viscosity dark matter model in the full Israel-Stewart formalism see [30]. A pioneer work about dissipative processes in cosmology is [31] and a highly recommended summary about relativistic fluid dynamics, dissipative relativistic fluids, applications to cosmology and astrophysics and bulk viscous perturbations is [32]. Bulk viscous dark energy has been studied in several contexts: as phantom dark energy [15,33,34]; dark energy with bulk viscosity observational constraints were studied in Refs. [35-38]; the present acceleration of the universe as effect of bulk viscosity of self-interacting scalar field [39];
An important problem in cosmology is that most of the dark energy models are able to adjust fairly well the observational data, this degeneracy hinders the tests and selection of more well-grounded dark energy models. The data shows only an accelerated cosmic expansion, but does not reveal the intrinsic nature of the source that causes this acceleration. Despite unknown, dark energy should be consistent with the known laws of physics; for that reason the dark energy when modeled as an exotic fluid must satisfy the bounds imposed by the laws of thermodynamics. In this vein, dark energy thermodynamics has been studied by several authors: a remarkable theoretical treatment in the context of perfect fluids can be found in Ref. [9]; thermodynamic properties of dark energy with varying equation of state parameter w(α) [10, 40, 41]; and, thermodynamic properties of dark energy as a self-interacting complex scalar field [42]. Studies about the thermodynamics of dark energy with viscosity: the validity of the generalized second law of thermodynamics in a non-flat universe in the presence of viscous dark energy [43]; some thermodynamic aspects of an alternative to the standard expression for bulk viscosity [44]; the thermodynamic stability analysis of non-interacting diffusive cosmic fluids with barotropic equation of state [45]; and, thermodynamics of viscous dark energy in the braneworld context [46,47]. Following this path, in this work we investigate some general thermodynamic aspects of dark energy modeled as a perfect fluid with constant equation of state parameter w = w 0, dynamical dark energy w = w(α), and bulk viscosity dark energy in the Eckart and Israel-Stewart theories.
This paper is organized as follows: in Sec. 2 is presented the general cosmological and thermodynamic considerations used all along the paper. The paper is divided into two parts, the first one, dedicated to dark energy as a perfect fluid, first, Sec. 3 summarizes the main results of the thermodynamics of dark energy with constant equation of state parameter, the same treatment but with dynamical dark energy is presented in Sec. 4. In second part, first in Sec. 5 are established the main ideas of a dissipative fluid in a flat, homogeneous and isotropic universe used in this work, then in Sec. 6 is presented the dark energy thermodynamic analysis in the Eckart framework to then in Sec. 7 in the Israel-Stewart. Finally, in Sec. 8 we presented the conclusions.
2. General considerations
In this section we present the general cosmological and thermodynamic considerations used all along the paper. Let us consider a Friedman-Lemaître-Robertson-Walker universe (FLRW), i.e., a homogeneous, isotropic spacetime, and the flat universe case. We are considering the late future time universe in which the dominant component is the dark energy, ∑ i ρ i ≈ ρ. In what follows we assume only one fluid as the main component of the universe, which experiments a dissipative process during cosmic evolution. With these considerations the Friedmann equations are
where the dot denotes derivatives with respect to the cosmic time,
In a FLRW cosmology the energy density conservation equation is
where ρ is the energy density and p is the pressure.
The thermodynamic assumptions are: the physical three dimensional volume of the universe at a given time can be expressed in terms of the scale factor V = V 0 α 3(t) (where V 0 is the three dimensional volume at the present time), the internal energy of the a cosmological fluid is U = ρV , and the first law of thermodynamics is expressed like
where T, S, µ and N are the temperature, entropy, chemical potential and the number of particles respectively.
The temperature will be assumed as a function of the number of particle density, n = N/V , and the energy density, ρ, therefore T = T(n,ρ). The last assumption gives the following useful relation [31,32]
These are the main general work hypothesis of this work.
3. Dark energy as a perfect fluid with w = constant
In this section we will present the main thermodynamic properties of dark energy as a barotropic perfect fluid with constant equation of state parameter. This section is a summary of the previous work [9].
The main considerations in this section are: a barotropic equation of state p = wρ (with the condition w < −1/3 to ensure the accelerated expansion) and the conservation of the number of particle current, n
α
= nu
α
(u
α
is the 4-velocity), which for a perfect fluid this is a conserved quantity
where w is a constant. The last two Eqs. (6) and (7) in the energy density conservation Eq. (3) and the first law of thermodynamics (4) gives:
Now we can calculate some relevant thermodynamic quantities.
3.1. Temperature and energy density
Assuming that the energy density is a function of the temperature and volume, ρ = ρ(T,V ), then
combining the last equation with the perfect fluid energy density conservation Eq. (8) we have
using now the relation of the temperature (5) and the equation of state (6), we obtain a relation for the temperature
Integrating the last equation we have
The last expression relates temperature, energy density and internal energy in the expected way for a perfect fluid, the temperature is directly proportional to the internal energy.
To calculate the energy density in terms of the scale factor we must integrate directly the perfect fluid energy density conservation Eq. (8) to obtain the well-known expression
This is the usual scale factor power-law energy density. Let us remind that in the ΛCDM case, for w = −1, which gives ρ = ρ 0, the energy density of the dark energy is a constant. When, −1 < w < −1/3, we have a dark energy that dilutes as the scale factor grows (but slower than dust or radiation). And, for the phantom case, w < −1, we have ρ = ρ 0 α r where r is a positive constant, in this case the energy density grows as the scale factor increase, and consequently the temperature and internal energy too.
To calculate the energy density in terms of the temperature we use (13) and (8)
where
3.2. Entropy
Using the Gibbs relation (9) and the perfect fluid energy density conservation Eq. (8) we can easily calculate a relation for the entropy
The energy density conservation equation and first law of the thermodynamics for a perfect fluid with w constant in FLRW imply that the entropy is a constant, i.e. is adiabatic. There is no entropy produced within the system (no friction, viscous dissipation, etc.) by the work done due to the cosmic expansion. Notice that this result is independent of the value of w.
Lets calculate the entropy, in the perfect fluid case the entropy of a cosmic fluid can be calculated from the well-known Euler relation:
we can easily calculate the entropy taking into account the relation U/U 0 = T/T 0 from (14), with this we have
For a perfect fluid the entropy S and the number of particles N are constants, consequently the ratio µ/T must be constant. We can define µ/T = µ 0 /T 0 and by the second law of thermodynamics, S ≥ 0, then we have the following relation
We have three options: First µ 0 = 0, in this case we have w ≥ −1 which is the phantom divide line restriction, that is, the phantom-like behavior of the dark energy fluid is forbidden. For µ 0 > 0, then w has a minimal value that does not reach the cosmological constant point w > −1 and again a phantom regime is forbidden. Finally, if the chemical potential is negative µ 0 < 1 then values of w < −1 are allowed, i.e, phantom dark energy. In Ref. [10] the authors examine these cases for several dark energy models.
4. Dark energy as a perfect fluid with w = w(α)
In this section we will study dynamical dark energy as a perfect fluid with the same background hypothesis of the previous section; the conservation of particle Eq. (7), the perfect fluid energy density conservation Eq. (8) and the first law of thermodynamics (9), hold in this case. But with the difference that now the equation of state has the form
The core idea of dynamical dark energy is that the equation of state parameter is a variable quantity, w(α), since it is constructed from two variable quantities.
4.1. Temperature and energy density
For a dynamical dark energy the relation U/U 0 = T/T 0 = (ρ/ρ 0)α 3 (14) still holds. When we consider p = w(α)ρ the equation (16) cannot be integrated immediately, it is necessary to know the explicit dependence with the scale factor of the equation of state parameter.
To calculate the temperature in terms of the energy density we may solve the equation of T(n,ρ) (5) in this particular case
whose solution using the method of characteristics yields
in which F is an arbitrary function. Notice that, for a the particular case of w = const., from equation (15) we have ρα 3(1+w) = ρ 0 = const. and from (7) we have nα 3 = n 0 = const.; then, ρ (1/w+1) /n = const. which imply that F(ρ 1/w+1 /n) = const. recovering the previous result ρ = ρ 0 T (1+w)/w of the equation (17).
To calculate the energy density in terms of the scale factor we integrate the energy density conservation equation (3).
Again, it is necessary to know the explicit dependence with the scale factor of the equation of state parameter in order to obtain a closed ρ(α).
4.2. Entropy
Using the first law of thermodynamics Eq. (4), the energy density conservation Eq. (3) and the equation of state (22) we can easily calculate
There is no entropy production in a universe filled with perfect fluid dark energy despite having a dynamical equation of state, w(α). This is an interesting result.
If we use the Euler relation U = TS − pV + µN, the temperature relation (24), the equation of state (27) and N = nV = n 0 V 0, then we have
If the chemical potential is null, µ = 0, or if it is of the form µ = µ 0 E, then there is no dynamical dark energy, w(α) = w = const. ≥ −1. A positive chemical potential, µ > 0, results in a w(α) > −1 avoiding the cosmological constant case. Finally, for µ < 0 and |µ| < (S 0 T/V 0 n 0), again we have w(α) > −1; and for |µ| > (S 0 T/V 0 n 0), we have the phantom case w(α) < −1 [10]. This model allows the phantom dark energy with positive definite temperature and entropy in certain cases of negative chemical potential.
5. Dissipative dark energy
A common feature of many of dark energy models is the assumption that it can be modeled as a perfect fluid. As stated before, a perfect fluid generates no entropy and no frictional heating because its dynamics is reversible and without dissipation. This approach works quite well in standard cosmology, but the observational reconstructions of the dark energy equation of state suggest the possibility of phantom dark energy which is incompatible with the perfect fluid hypothesis at thermodynamic level. Real fluids behave irreversibly, for this reason, if we want to maintain the fluid hypothesis to dark energy and ensure its thermodynamic compatibility, it is interesting to consider dissipative processes like bulk viscosity. Qualitatively, the bulk viscosity can be interpreted as a macroscopic consequence coming from the frictional effects within the fluid.
For a dissipative fluid, the particle 4-current will be taken to be of the same form as
where Π is the bulk viscous pressure. We are considering in the last Eqs. (32) that the thermodynamical system is in the near equilibrium regime, if the fluid is out of equilibrium as a result of dissipative effects then there is no unique average 4-velocity spoiling the formalism [31]. In the presented case, the near equilibrium condition occurs when the local equilibrium pressure p is dominant over the viscous pressure (which encodes deviations from equilibrium), i.e., when the condition
is satisfied. This condition is the same for both non causal and causal approaches. For notation, from now we will drop the bar of the thermodynamic variables and write the pressure as p eff = p+Π, where p is the usual barotropic pressure p = wρ. The energy density conservation equation in this case is
From the first law of thermodynamics (4) and remembering that for the energy U = ρV , volume V = V 0 α 3 and number of particle density n = N/V , can be derived an expression for the entropy in terms of the bulk viscous pressure
The above equation sets a strong general constraint in this framework: due to n and T are positive, and H in an expanding universe is positive too to ensure a non-negative entropy production (second law of thermodynamics) the bulk viscous pressure necessarily must be Π ≤ 0.
6. Eckart bulk viscosity dark energy model
The simplest approach to treat bulk viscosity in cosmology is the Eckart theory [17]. This approach has some important limitations, the main one is that it is a noncausal approach to dissipative processes.
In this framework the bulk viscosity pressure is given by
where the bulk viscosity depends on the function ξ and the Hubble parameter. This formalism has been widely used at background level. Parameter constraints from observational data of the Eckart bulk viscosity dark energy can be found in Refs. [35, 36]; in both works the main results are that the Bayesian evidence show that the Eckart viscous dark parameters are small and statistically indistinguishable from ΛCDM with the current data and the effective EoS is slightly phantom. First, in Ref. [35] the authors study three Eckart viscous dark energy models:
From (34), the energy density conservation equation is
First, we can sketch some general features without determining the functional form of ξ. From the Friedmann Eqs. (1) and (2) we have
Note that for the case without viscosity ξ = 0 and w = −1 the last equation boils down into
where C is an integration constant. Integrating again, for w ≠ −1, we obtain an expression for the scale factor:
where D is a another integration constant.
6.1. Temperature and energy density
As in the perfect fluid case, to calculate the temperature in terms of the energy density we have to solve the Eq. (5) whose solution is (26)
using the Eqs. (3) and (7) in
Then the Eq. (41) reduces to (16) whose solution has already w been calculated,
In order to calculate specific expressions for the energy density and temperature lets analyze two particular popular proposals for ξ.
Model 1 ξ(t) = ξ
0 = const.: this is simplest Eckart bulk viscosity model. Solving (38) for model 1 we have
Integrating again with the condition α(t 0) = 1 we obtain the scale factor
from which can be calculated the energy density in terms of the scale factor
with the constant W = ξ 0 /(H 0(1 + w)). The near equilibrium condition (33) for ξ(t) = ξ 0 implies |ξ 0| ≪ H 0, and consequently for −1 < w < −1/3 imply that W is small. Notice that the energy density has the form ρ = c 1 + c 2 α −3(1+w) + c 2 α −3(1+w)/2 , where ci are constants.
The first term behaves as the cosmological constant, the second as a perfect fluid in ΛCDM and the last as a slow diluting fluid corresponding to the bulk viscosity effect.
Model 2 ξ = (β/√3)ρ 1/2 : in Ref. [49] was first assumed that the viscosity has a power-law dependence upon the energy density, ξ = αρ s , where α > 0 and s are constant parameters. In model 2 is a great difficulty to obtain easily manageable solutions to the main equations, only some particular results have been found. In the special case where the bulk viscosity coefficient takes the form ξ(ρ) ∝ ρ 1/2 , a Big Rip singularity solution was obtained in this formalism for late future times FLRW flat universe filled with only one barotropic fluid with bulk viscosity [48]. For s = 1/2 and α = β/√3 in ξ = αρ s , the Eq. (38) simplifies into
which can be easily integrated to calculate the Hubble function and the scale factor
With this, we can calculate the energy density
The energy density is a power-law of the scale factor. It is important to notice that for w = −1 dark energy the exponent is positive and consequently is an increasing energy density and temperature. The near equilibrium condition (33) in this model ξ(t) = (β/√3)ρ 1/2 implies β ≪ |w|. For values w ≠ −1, some combinations of the parameters can result into different evolutions of the energy density. For a decreasing energy density (diluted by the cosmic expansion like occurs with ordinary matter and radiation) it is needed negative powers, i.e., β < w + 1 must hold; observations show that w + 1 is close to zero (and positive or negative), so for a decreasing energy density the conditions w +1 > 0 and β ≪ 1 are necessary. Phantom dark energy, w +1 < 0, or relatively big values of β, always yields β > w + 1, the increasing energy density.
6.2. Entropy
We calculate the entropy using the Eq. (35), substituting the Eckart the bulk viscosity pressure, Π = −3ξH, we have
As in the previous subsection we will treat two cases, first we calculate the thermodynamic variables n and T, substitute in Ref. (50) and then integrate.
Model 1 ξ(t) = ξ 0 = const.: first note that a growing entropy in time, dS/dt > 0, immediately imply that the bulk viscosity parameter is a positive constant ξ 0 > 0. Considering only this option the Hubble function is given by Eq. (43).
The particle density n = n 0 α −3 can be easily calculated using (44)
And the temperature
Substituting these two previous expressions in Eq. (50) we have a differential equation for the entropy that depends only on time
A quick examination of the last equation gives a constriction to the equation of state parameter w, in order to ensure the second law of thermodynamics, it must satisfy w > −1, i.e. the phantom case is forbidden for a barotropic fluid. Another possibility is the to allow unusual assumptions like a negative temperature
Equation (53) can be integrated easily if we assume w = constant,
where S 0 is the entropy at the present time. The entropy in this particular model (54) grows exponentially.
Model 2 ξ = (β/√3)ρ 1/2 : Repeating the same procedure of the previous case, calculating the particle density, temperature and substituting in Ref. (50); first we obtain the differential equation for the entropy
Let’s quickly examine Eq. (55) before solving it. First, clearly for β = 0 (no bulk viscosity) then dS/dt = 0, i.e. the entropy is a constant in time as expected for a perfect fluid. For small values of the bulk viscosity function parameter, 0 < β ≪ 1, we have δ ≈ 0, which implies that dS/dt ∝ (1 + [3H 0 /2][1 + w − β][t − t 0])−1, the condition of increasing entropy dS/dt > 0 requires that w ≥ −1.
Integrating we obtain the expression of the entropy
again, S 0 is the entropy at the present time. The entropy in model 2 is a power-law. An increasing entropy condition imposes δ > 0, from which we have two cases: for w + 1 > 0 (quintessence) it is required β < 1 + w and as we know the equation of state is close to the cosmological constant case, i.e. 0 < w + 1 ≪ 1 then β has to be very small; the second case is w+1 < 0 (phantom) which implies 1+w < β which it is easily fulfilled.
7. Israel-Stewart bulk viscosity dark energy model
The Israel-Stewart theory provides a better description than the Eckart theory. It is a causal and stable theory of thermal phenomena in the presence of gravitational fields. This theory besides solving the non-causal problem of the Eckart theory, enrich the framework including new features like the entropy has terms of second order in the dissipative variables and incorporates transient phenomena on the scale of the mean free path/time, outside the quasi-stationary regime of the classical theory. In the Israel-Stewart theory we have the same Friedmann equations with an equation for the causal evolution of the bulk viscous pressure given by
The last equation is known as the transport equation of the viscous pressure Π. Where τ is the relaxation time for bulk viscous effects (in the limit τ = 0 the theory is noncausal), ξ the bulk viscosity coefficient and T the temperature.
It is clear that the Israel-Stewart theory is much more complex than the Eckart theory. In addition to having to propose or infer a function for ξ and for the relaxation time τ, we must know the temperature T and solve a differential equation involving these physical quantities. Let’s see how we can make a proposal and look under which conditions an analytical solution of Eq. (58) can be obtained.
First for the temperature, following [15, 31, 32], we assume that the temperature depends only on the particle number an energy density, T = T(n,ρ), with this
The conservation equations of energy density and particle density in the Israel-Stewart theory are
substituting both expressions in Eq. (60) and using (5) we have
If w is not a function of ρ, then (∂p/∂ρ) n = w (which w is not necessarily a constant). Taking this into account this, it is easy to prove that the expression for the temperature
is solution of (63). Which is the same as for the perfect fluid case and Eckart bulk viscosity.
For ξ, we assume a power-law for the bulk viscosity in terms of the energy density of the fluid [49],
where s is a constant arbitrary parameter and ξ 0 a positive constant.
And finally we follow [48] in which a simple relation between τ and ξ is proposed
Equations (65) and (66) are the two main hypotheses assumed in order to find a solution for the transport equation of the viscous pressure (58).
Combining the Friedmann Eqs. (1), (2), with the pressure in this context p eff = p + Π = wρ + Π, it is easy to obtain
With the Eq. (67), a differential equation for the Hubble parameter can be constructed substituting (64), (65) and (66) into the transport Eq. (58), then differentiate (67) and equal both expressions considering the energy density conservation Eq. (61). The resulting equation for the Hubble parameter after these considerations is [48]:
There is not a general analytical solution to the last equation, it may be solved numerically for specific cases of the parameters of the model. In Ref. [15] the authors solved this equation for several values of s numerically. The authors found that for s ≠ 1/2 (or s ≤ −1/2) there is no phantom solution. However, for the particular case of s = 1/2, i.e. ξ = ξ 0 √ρ, the Eq. (68) has a particular analytical solution. Under this assumption the equation reduces to
Rewriting the last equation we have
Notice that clearly if the bulk viscosity parameter is zero ξ 0 = 0, then (from the transport Eq. (58)) the bulk viscosity pressure is null Π = 0, and we recover the standard Hubble equation
For the cosmological constant case, w = −1, Eq. (70) boils down into the well-known de Sitter case
Henceforth we will assume ξ 0 ≠ 0 and w + 1 ≠ 0 to to exclude the standard case. For mathematical convenience, we define
then the Eq. (70) transforms into
To sketch a solution, the following Ansatz is proposed for the Hubble function in the late future universe [15,48,50,51]:
where t s is a finite time in the future at which the Big Rip occurs and A is a positive constant in order to describe an expanding universe. For s ≠ 1/2 it can be seen by direct inspection that a more general Ansatz with the form, H(t) = A(t s − t) p with p < 0, does not reduce the differential equation to an algebraic polynomial in A [15].
Substituting the Ansatz (75) into Eq. (74) it is obtained a quadratic equation for the constant A
whose solution is
Integrating (75) we obtain the scale factor as a function of time, [15,50]
where t 0 is the present time. With this it can be easily calculated the number of particles
It is clear that at the time t = t s the size of the universe becomes infinite and the number density of particles goes to zero.
Finally, with the Ansatz (75) we calculate the bulk viscous pressure:
Notice that the bulk viscous pressure is divergent at the Big Rip time.
7.1. Temperature
To calculate the expression of the temperature as a function of the scale factor we use the temperature as a function of the energy density Eq. (64), and the equation of the Hubble parameter (75) and the scale factor (78):
with T 0 the temperature today. It is a power-law in terms of the scale factor. For the quintessence case, −1 < w < −1/3, the power is always negative, giving a temperature that decrease as the universe expands; whereas for phantom dark energy the power is positive giving an increasing temperature. Once again, the characteristic behavior where the temperature increases as a function of the scale factor, becomes present.
7.2. Entropy
The entropy change can be evaluated from (35),
To calculate the entropy we integrate (83)
where
and since A > 0 then, the last inequality is satisfied for values 0 < w < 1/2. This is an interesting result, in the IsraelStewart framework, for the particular case ξ = ξ 0 ρ 1/2 ; the natural conditions of positive entropy (η + 1 < 0) and expanding universe (A > 0) imply that the equation of state parameter must be positive, w > 0. This can be interpreted as the fluid is not a “barotropic dark energy” (w < −1/3), but, it is important to take in mind that we have discussed the thermodynamic properties of a Big Rip solution H = A(t s − t)−1 in the framework of the Israel-Stewart theory and despite it gives w > 0, the cosmic expansion is accelerating due to the global effect of the bulk viscosity in which w = p/ρ < −1/3. The solution is for a cosmological scenario with barotropic fluid p = wρ and the constriction 0 < w < 1/2, but, globally behaves like a phantom fluid, due to the bulk viscosity provides the sufficient negative pressure, allowing to cross the phantom divide line.
8. Conclusions
We have discussed in the present work a general treatment for dark energy thermodynamics in several scenarios for the late future time: by considering it as a perfect fluid with constant and variable equation of state parameter w, and as a dissipative fluid with bulk viscosity both in the Eckark and Israel-Stewart frameworks. General equations have been derived for the dark energy temperature, energy density and entropy in a flat, homogeneous and isotropic universe. From the results of the entropy evolution, some theoretical thermodynamic constraints are imposed for the w parameter in order to satisfy a positive the entropy and grow rate.
We first recall some main results of dark energy as a perfect fluid with constant w; previously presented in a clear way in the work [9]. In this context the temperature is directly proportional to the internal energy T ∝ U, the energy density is a power-law of the temperature ρ ∝ T
(1+w)/w
, the entropy is constant in time (adiabatic)
In the second part of this work, the dark energy is modeled as a dissipative fluid with bulk viscosity where the (effective) pressure is assumed as the sum of the barotropic and the bulk viscous pressures p eff = p + Π. This is the core of this work. Within the Eckart framework the viscous pressure is given by Π = −3Hξ(t) and two cases were studied: ξ(t) = ξ 0 = const. and ξ = (β/√3)ρ 1/2 . Both cases have the relation T ∝ ρ w/(1+w) . In the first case, the entropy grows exponentially and puts the condition w > −1 to preserve the entropy and temperature positive, or, if w > −1 then T < 0. Although the possibility of negative temperature for phantom energy has been studied [40], it seems a cumbersome hypothesis and therefore quintessence is favored. In the second case, the entropy grows as a power-law, and quite an interesting fact in this case is that the phantom regime is allowed under not so strict conditions. The next step is to use the causal Israel-Stewart theory. We have solved the bulk viscous pressure transport equation assuming ξ = ξ 0 ρ 1/2 and τ = ξ/ρ which for the late future time universe gives rise to Big Rip cosmological solutions, H = A(t s −t)−1 [15,48,50]. Under these particular assumptions the temperature evolves as T = T 0 α 2w/A(w+1) , that is, in the quintessence regime the temperature decreases as the universe expands, whereas the temperature increases in the phantom regime. This is a characteristic of phantom and Big Rip solutions, here the energy density increases with the cosmic expansion as well the internal energy and consequently the temperature too. Calculating the entropy we obtain a power-law in time and imposing the conditions A > 0 and S > 0, an important constraint 0 < w < 1/2 is derived. The solution implies a cosmological scenario of a barotropic fluid with 0 < w < 1/2 which behaves like a dark fluid (even phantom) driven by the bulk viscosity. So, in the Israel-Stewart formalism, some solutions behave as dark energy and even allow to cross the phantom divide line with an effective equation of state lesser than −1 and constant in time without evoking an exotic barotropic fluid with w < −1. An important extension to the study done is to obtain other cosmological solutions additional to H = A(t s − t)−1 and constrain the model parameters with observational data to track the most plausible model.
Despite its intrinsic nature is not well understood yet, reconstructions of the dark energy equation of state from observational data seem to suggest a time-dependent w and even the crossing of the phantom divide line w = −1 from above to below is not only possible but could indeed be a condition for a successful description of observations. For this reason it is important to examine the phantom regime carefully. It has important problems and if modeled as a perfect fluid, it unfeasible from the point of view of classical thermodynamics. Hence, new perspectives beyond the perfect fluid approach in standard cosmology have to be considered in order to have a better understanding of the phantom regime, such as fluids with bulk viscosity that reduce the kinetic pressure of the cosmological fluid and provide a richer thermodynamic framework. In fact, one could actually dive deeper and seek for the foundations of dark energy thermodynamics from statistical physics (and possibly quantum mechanics) to establish basic general principles for a dark energy theory, or if it is the case, give solid arguments in favor of alternatives like modified gravity avoiding the hypothesis of dark energy as a substance; all this important but long beyond the scope of this work.
Finally, dark energy thermodynamics and bulk viscous dark energy are topics that have been studied from several different approaches for years, in this work we summarize particular results previously obtained by other authors on both subjects to make a broader and unified study, extending its scope in order to search for insights about the nature of the dark energy from a general theory that is classical thermodynamics and more general fluids than the perfect fluid.