# The van der Waals fluid and its role in cosmology

## Abstract

We analyze the properties of a generic cosmological fluid described by the van der Waals equation of state. Exact solutions for the energy density evolution are found as implicit functions of the scale factor for a flat Friedmann-Robertson-Walker space-time. The possible values of the free parameter in the van der Waals equation are selected a posteriori, in accordance with asymptotic behaviors that are physically relevant. The stability of the model against small perturbations is studied through the hydrodynamic perturbations of the fluid for the relevant cases. It is found that a van der Waals fluid seems appropriate to describe non-eternal inflationary scenarios.

## 1 Introduction

With the advent of the inflationary cosmology [1, 2] and the relatively recent discovery of the accelerated expansion of the Universe at late times [3, 4], the cosmological components enlarged and exotic equations of state were needed to model the acceleration of the cosmic expansion [5, 6, 7]. The Chaplygin gas type equation of state is frequently studied in the literature and used to describe standard matter and dark energy as a single fluid [8, 9, 10, 11, 12, 13, 14]. Such equations of state can be generated by a standard scalar field with a specific potential [8, 15] or a tachyon scalar field [16], which can be also applied in various other contexts [17, 18, 19].

In the concept of standard matter and dark energy as a single fluid, the van der Waals equation of state has received some attention and was investigated in the literature [20, 21, 22]. The van der Waals fluid as a cosmological constituent in a variety of other contexts was also investigated [23, 24, 25, 26, 27]. Due to the non-linearities of such equations of state, the cosmological dynamics is usually solved numerically and asymptotic analytic solutions can be obtained just in some special cases. The main feature explored in these models concerns the ability to reproduce the accelerated and matter-dominated ages with a single component. Further, depending on the characteristics of the model in question, the properties of the early and/or the late Universe may be possible to describe.

In this work, we do not impose any desired behavior on the van der Waals fluid. Instead we look for general analytical solutions of the energy density in a flat Friedmann-Robertson-Walker space-time, and so select the physically relevant cases. Once the suitable cases are selected, the corresponding stability is tested through the analysis of the hydrodynamic perturbations.

## 2 Basic equations

Let us consider a cosmological fluid that obeys the van der Waals equation of state, which is written in the following form [23]

(1) |

where and are the pressure and energy density of the fluid, respectively. The parameter is a non-negative constant. This non-linear equation of state presents some peculiar qualitative features. It has a vertical asymptote in , where the pressure diverges, and its form approaches the form when is sufficiently small while the pressure is persistently negative when is large enough. Further, the pressure can eventually switch from a positive value to a negative one and vice versa. Some of these features are very compelling for cosmology and justify a more detailed investigation of the van der Waals fluid in this context.

In a flat Friedmann-Robertson-Walker (F-R-W) background, assuming that the perfect fluid conservation equation holds, we have

(2) |

with the dot denoting time derivative and the scale factor.

From equations (1) and (2), one obtains a differential equation for and , namely

(3) |

which can be rewritten in the form

(4) |

This differential equation furnishes the necessary information to determine the energy density evolution of the fluid.

By assuming that general relativity governs the large scale dynamics, the energy density of the van der Waals fluid has an effect on the cosmic evolution through the Friedmann equation

(5) |

where stands for the Hubble parameter.

## 3 Background solutions

With the above equations in hand, we can now try to obtain the functional behavior of the energy density of the van der Waals fluid in a flat F-R-W background. Equation (4) is separable in terms of functions of and , generating the following differential expression

(6) |

The corresponding integration furnishes

where is a constant of integration. The integral in the right hand side of the above expression contains a polynomial of the form , with , and , whose integration for gives

(8) |

Here is a constant and this solution is valid for .

In the case for , the solution is

(9) |

which is valid only for .

Finally, the solution corresponding to the case reads

(10) |

which holds for .

All the above solutions have asymptotic roots, i.e. they allow values of that satisfy , meaning that the equation of state considered here is compatible with what is expected for the Universe at early times. As we can realize, expressions (8), (9) and (10) involve terms quadratic in and one can expect two roots for each of them. However, there is only one root common to all solutions and another is a root only of (8) and (9). Let us see this in more details as follows.

Expressions (8), (9) and (10) satisfy

(11) |

So is an asymptotic root for all the above solutions. But notice that for , solution (8), we have that (11) is complex and must be discarded for physical purposes. Then the value given by (11) is physically valid for and is the asymptotic initial condition from which the Universe expands as decreases. Here it is important to point out that such a property means that a van der Waals fluid does not have a divergent energy density when , i.e. it does not present a big bang singularity. Interestingly enough, is the result from equation (1) when evaluated at , namely,

(12) |

This result establishes that for a sufficiently small the van der Waals dynamics for , encompassed by (9) and (10), is close to the one of the cosmological constant. Hence, the equation of state holds in the beginning of the expansion and we have a de Sitter-like inflation for . One more important point to stress is that once we take initial conditions around (11), the energy density will never take values near the vertical asymptote, , given that and always will decrease with the increase of .

The second asymptotic root is

(13) |

which holds only for (8) and (9) - solution (10) diverges for . This root also produces a complex for , but even for , solution (9), it is not physically motivated as an asymptotic initial condition because it leads to limiting equations of state that are not viable.

Now consider the dynamics in the range , were the initial condition (11) holds and the energy density dilutes as the Universe expands. Then, from solutions (9) and (10) it is straightforward to see that when , i.e. is large, we have

(14) |

This is in accordance with (1), from which results in the referred limit.

As follows, we will analyze the cases of physical interest: and .

1) Case

It corresponds to the equation of state , see equation (1), which through (2) furnishes the differential expression

(15) |

The respective integration renders

(16) |

where is a constant of integration.

By taking the Friedmann equation, , and applying (16), we have the differential expression bellow

(17) |

where the change of variable was performed.

Integration of (17) gives the implicit time evolution of the scale factor, namely,

(18) |

with being a constant of integration.

The above expression has the asymptotic behaviors

(19) |

for a sufficiently small , and

(20) |

corresponding to a large . Hence, this case describes an inflation followed by a matter-dominated era type dynamics.

2) Case

The range of validity of solution (10) comprises , from which we have

(21) |

From the general result (11), this expression has the root , which is the initial energy density and corresponds to the equation of state , as stated by (12). Also from (21), when is sufficiently diluted we have , i.e. a radiation-dominated era, as anticipated in (14). Thus this case describes an inflationary period followed by a radiation-dominated era.

## 4 Hydrodynamic perturbations

In this section we will develop the general equations to evaluate the small scalar perturbations in a flat F-R-W background for a van der Waals fluid. As the background fluid does not present anisotropy, the line element for small scalar perturbations in the longitudinal gauge is [28, 29, 30]

(22) |

where is the conformal time, , and is the gauge-invariant perturbation potential.

By introducing the Mukhanov-Sasaki variables [29, 30], which are defined by

(23) |

the adiabatic perturbations, in the general relativity context, are described by the following compact differential equation

(24) |

where the prime represents derivative with respect to the conformal time and is the speed of sound.

For a plane-wave solution of the form , the short-wavelength perturbations hold the condition and the last term of equation (24) can be neglected, furnishing

(25) |

The corresponding solution, valid for a sufficient slowly varying speed of sound (WKB approximation) reads

(26) |

where is a constant of integration.

For long-wavelength perturbations, when the condition is valid, the spatial derivatives of (24) can be neglected and the following solution holds

(27) |

All the quantities appearing in the asymptotic solutions above can be expressed in a compact form as functions of the energy density . The Friedmann equation furnishes the infinitesimal conformal time interval in the form

(28) |

On the other hand, by using equation of state (1) one directly determines the speed of sound, namely,

(29) |

Also equation (1) allows to write the variable given by as a function of in the form

(30) |

Furthermore, the variable is also equivalent to

(31) |

where is the Hubble parameter in terms of the conformal time. To get (31) from (30) one uses the acceleration equation

(32) |

along with the equation of state (1). From the variable in the form (31) one can write the integral in (27) as

(33) |

through an integration by parts. Using this result in (27) and returning to , with a plane-wave solution of the form , finally we have

(34) |

where (28) was used and is a constant defined by . With the quantities (26), (28), (29), (30) and (34) in hand, we can evaluate the perturbations for any and determine the density contrast. The perturbed dynamics for the cases of interest will be analyzed as follows.

Let us now apply the above results to the case 1, when , and case 2, when . We immediately infer from equation of state of the case 1, , that the speed of sound (29) is always imaginary, which will generate perturbations with increasing amplitude all the time. Hence, the case 1 is unstable against small perturbations. The same occurs for the case 2 when is small, since in this situation and solution (26) will give a combination of real exponentials for short-wavelength perturbations. However, the equation of state for the case 2 rapidly evolves to (radiation) as dilutes with inflation, so reaching a real speed of sound and it is expected that persistent growing instabilities are prevented. Thus from now on we are focused on the perturbations for the case with , which is in fact the most interesting one since it naturally describes inflation followed by a radiation-dominated era.

Substituting (21) into (34), we get the perturbation potential for in the long-wavelength limit as a function only of , namely,

(35) |

where

(36) |

and denotes the initial energy density corresponding to . Here the constants and are set to unity without loss of generality.

With this result in hand, using the following perturbation equation from the perturbed Einstein’s equations

(37) |

where is the gauge-invariant density perturbation, we can now infer the density contrast.

Since we are analyzing the long-wavelength limit, the first term of (37) can be neglected (spatial derivatives). Further, using the Friedmann equation in terms of the conformal time, , equation (37) furnishes

(38) |

which gives the density contrast if is known.

If we derive (34) with respect to conformal time, one has

(39) |

Now we use (31) to eliminate in the last expression, which is then substituted into (38). The corresponding result can be expressed in terms of by means of (28) and (30). So, by considering the density perturbation of the form , one obtains

(40) |

Thus we have an expression that gives the density contrast for as a function only of .

The scale factor (21) and amplitudes of the density contrast (40) and perturbation potential (35) are plotted in the left frame of Fig. 1. By observing the corresponding curves, we can see that the amplitude of the perturbation potential increases slowly with the dilution of , i.e. with increasing . On the other hand, the amplitude of the density contrast increases initially with the dilution of and starts to decrease when is sufficiently diluted, i.e. when is sufficiently large and a radiation-dominated era is expected.

An alternative visualization of the perturbations is shown in the right frame Fig. 1, where the amplitudes of the perturbation potential (35) and density contrast (40) are expressed as functions of the scale factor. From this figure it is clear to see the slow growth of the amplitude of the perturbation potential with the cosmic expansion while the amplitude of the density contrast grows in the beginning and starts to decrease after a certain point of the cosmic expansion.

In view of these results, the van der Waals fluid could describe inflation followed by a radiation era, i.e. a non-eternal inflation which is stable against small perturbations in the limit of the long-wave perturbations, as shown by the above analysis of the hydrodynamic perturbations.

## 5 Concluding summary

In this work, we studied a generic cosmological fluid described by the van der Waals equation of state for a flat F-R-W space-time. General background solutions for the energy density of the fluid were obtained in analytical form. From these solutions it was possible to identify the physically relevant cases. It was found that in general the energy density of a cosmological van der Waals fluid does not diverge as the scale factor goes to zero (no big bang singularity) and behaves like a cosmological constant type component. Two cases of interest in cosmology were found and they correspond to the values 0 and 1/8 for the parameter of the equation of state. For the background solution describes an inflation followed by a matter-dominated era while for it describes an inflation followed by a radiation-dominated era.

The viability of the relevant background solutions is tested through the analysis of the hydrodynamics perturbations. The case for is unstable against small perturbations since the speed of sound is persistently imaginary. On the other hand, the case for is stable against small perturbations in the limit of long-wave perturbations. This entails that the van der Waals model presents some skill in describing a non-eternal inflation, i.e. an initial quasi-cosmological constant dynamics followed by a radiation-dominated era.

## Acknowledgments

M. H. B. C. acknowledges the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and G. M. K. the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Brazil, for financial support.

### Footnotes

- rudijantsch@gmail.com
- marcushbc@gmail.com
- kremer@fisica.ufpr.br

### References

- A. H. Guth, Phys. Rev. D 23 (1981) 347.
- A. Linde, Phys. Lett. B 108 (1982) 389.
- A. G. Riess et al., Astron. J. 116 (1998) 1009.
- S. Perlmutter et al., Astrophys. J. 517 (1999) 565.
- R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80 (1998) 1582.
- P. J. Steinhardt, I. Zlatev and L. Wang, Phys. Rev. Lett. 82 (1999) 896.
- P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75 (2003) 559.
- A. Yu. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B 511 (2001) 265.
- J. C. Fabris, S. V. B. Gonçalves and P. E. de Souza, Gen. Rel. Grav. 34 (2002) 53.
- G. M. Kremer, Gen. Rel. Grav. 35 (2003) 1459.
- J. C. Fabris, H. E. S. Velten, C. Ogouyandjou and J. Tossa, Phys. Lett. B 694 (2011) 289.
- S. Kumar and A. A. Sen, J. Cosmol. Astrop. Phys. 10 (2014) 036.
- P. P. Avelino, K. Bolejko and G. F. Lewis, Phys. Rev. D 89 (2014) 103004.
- E. O. Kahya and B. Pourhassan, Mod. Phys. Lett. 30 (2015) 1550070.
- C. Uggla, Phys. Rev. D 88 (2013) 064040.
- J. C. Fabris, T. C. C. Guio, M. Hamani Daouda and O. F. Piattella, Grav. Cosmol. 17 (2011) 259.
- P. P. Avelino, L. Losano and J. J. Rodrigues, Phys. Lett. B 699 (2011) 10.
- R. K. Jain, P. Chingangbam and L. Sriramkumar, Nucl. Phys. B 852 (2011) 366.
- R. C. de Souza and G. M. Kremer, Phys. Rev. D 89 (2014) 027302.
- S. Capozziello, S. De Martino and M. Falanga, Phys. Lett. A 299 (2002) 494.
- G. M. Kremer, Phys. Rev. D 68 (2003) 123507.
- S. Capozziello, S. Carloni and A. Troisi, Recent Res. Dev. Astron. Astrophys. 1 (2003) 625.
- G. M. Kremer, Gen. Rel. Grav. 36 (2004) 1423.
- S. Capozziello, V. F. Cardone, S. Carloni, S. De Martino, M. Falanga, A. Troisi and M. Bruni, J. Cosmol. Astrop. Phys. 04 (2005) 005.
- V. F. Cardone, C. Tortora, A. Troisi and S. Capozziello, Phys. Rev. D 73 (2006) 043508.
- B. Saha, Int. J. Theor. Phys. 45 (2006) 952.
- M. Khurshudyan, B. Pourhassan and E. O. Kahya, Int. J. Geom. Met. Mod. Phys. 11 (2014) 1450061.
- J. M. Bardeen, Phys. Rev. D 22 (1980) 1882.
- V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215 (1992) 203.
- V. Mukhanov, Physical Foundations of Cosmology (Cambridge University Press, Cambridge, 2005).