Universe consisting of diffusive dark fluids: thermodynamics and stability analysis

The present work deals with homogeneous and isotropic FLRW model of the Universe having a system of non-interacting diffusive cosmic fluids with barotropic equation of state (constant or variable equation of state parameter). Due to diffusive nature of the cosmic fluids, the divergence of the energy momentum tensor is chosen to be proportional to the diffusive current. The thermodynamic stability analysis of individual fluids is done and the stability conditions are expressed as restrictions on the equation of state parameter.


II. INTRODUCTION
The greatest challenge of standard cosmology today is to accommodate the recent observational predictions [1][2][3][4][5][6][7][8][9]. The modern cosmology is facing the challenging issue of explaining the present accelerated expansion of the universe. In the framework of Einstein gravity, cosmologists are speculating some hypothetical exotic matter (known as dark energy (DE) having large negative pressure) to explain this accelerating phase. It is estimated [10] that about 70 percent of the cosmic fluid consists of this unknown DE component. The simplest as well as the common candidate for this dark fluid is the cosmological constant (the zero point energy of the quantum fields). Although a large number of available observational data are in support of this cosmological parameter as a DE candidate but there are severe drawback of it at the interface of Cosmology and Particle Physics : the cosmological constant problem [11,12] and the coincidence problem [13]. As a result there are several alternative proposed dynamical DE models into the picture. These DE models have been studied for the last several years [14], yet the cosmological constant is still the best observationally supported DE candidate. However, these different dynamic DE models cannot be compared from observational view point as these models try to adjust the data seamlessly. As a result, cosmologists have been trying with interacting DE model to have a better understanding of the mechanism of this cosmic acceleration.
From recent past, interacting dark fluid models have been receiving much attention as they can provide small value of the cosmological constant and due to their ability of explaining the cosmic coincidence problem [15][16][17]. Moreover, recent observed data [18][19][20][21][22][23] favor interacting dark fluid models and it is possible to have an estimate of the coupling parameter in the interaction term by various observations [21][22][23][24][25][26][27][28][29][30]. Further,the cosmologists are of the opinion that these interacting dark fluid models may have the solution to the current tensions on σ 8 and the local value of the Hubble constant (H 0 ) [22][23][24][25][26][27][28][29][30][31]. Moreover, it is speculated that cosmological perturbation analysis may be affected by the interaction terms and consequently, the lowest multipoles of the CMB spectrum [32,33] should have an imprint of it.
On the other hand, the unknown nature of DE may have some clue from the thermodynamic laws which are applicable to all types of macroscopic systems and are based on experimental evidence. However, unlike classical mechanics or electromagnetism, thermodynamical analysis can not predict any definite value for observables, it may only give limit on physical processes.
So it is reasonable to believe that thermodynamical study of dark cosmic fluid may indicate some unknown character of it. Investigation in this direction has been initiated recently by Barboza et al [34]. But their stability conditions demand that the DE should have constant (-ve) equation of state and it is not supported by observation.
Subsequently, Chakraborty and collaborators [35][36][37][38]  The change of the energy-density is proportional to the particle density. The resulting nonequilibrium thermodynamics is studied and stability conditions are explicitly determined. The paper is organized as follows: Section III deals with thermodynamical analysis of non-interacting cosmic fluids with constant equation of state . The stability criteria for non-interacting cosmic fluids with variable equation of state has been organized in Section IV . Finally, the paper ends with a brief discussion and the conclusion from the result in dark energy range in Section V .
where ω i ,a constant ,is the barotropic index of the fluid and p i and ρ i are the pressure and energy density of the fluid respectively. So the Friedmann equations for the whole system take the form where H = a −1ȧ is the Hubble parameter , a(t) is called the scale factor of the universe as all physical distance is scaled with same factor " a " due to homogeneity and isotropy of space-time.
Now due to the diffusive nature of the fluids, they do not obey the matter conservation equation (T µν i ;µ = 0), rather they follow, where N ν i is the current of diffusion corresponding to that fluid and k is a constant. But for the whole universe,the total stress-energy tensor is conserved i.e.
In simplest form, the conservation equation(non conservation) takes the form for FLRW universe [39][40][41] as, where, γ i is a constant for a particular fluid but it is different for different fluids with i γ i = 0. Now, in the context of classical thermodynamics,the contribution to the energy of the universe by i-th fluid is given by where, v = v 0 a 3 (t) is the physical volume of universe at given time t and at present time t 0 ,the physical volume is v 0 with a(t 0 ) = 1. According to the 1st law of thermodynamics(which is an energy conservation equation), Integrating equation (6) one can easily find the expression for energy density of the i-th fluid as, where ρ i0 is the value of ρ i at t = t 0 i.e. at the present time .The energy density of i-th fluid, from equation (9) can be arranged as where , taking entropy S i as a function of two independent thermodynamic variables namely volume(v) and temperature (T ) and also considering dS i to be an exact differential, one obtains from equation (8), Now combining equations (10) and (11), a new relation among energy density ,temperature and time can be written as where β is an integration constant. Effectively at present time t = t 0 , Now equations (12) and (13) imply the evolution of energy of the i-th fluid with time and temperature as, This can be termed as the modified equation of state of this fluid where as before E i = ρ i v and A. Determination of thermodynamic derivatives : C p , C v , K S , K T Now, the motivation of the present work is to find the criteria of thermodynamic stability of the universe . So the conditions which may be satisfied by the parameters , indicate the thermodynamic features of the evolution of the expanding universe. For any system the thermodynamic derivatives i.e. the heat capacities,compressibility and isobaric expansibility determine whether the system is stable or not. So it will be very important to analyze these parameters for each fluid in the present context.
From the 1st law of thermodynamics in equation (8), one can determine the two heat capacities : heat capacity at constant volume C v and heat capacity at constant pressure C p by the relations , where H = E + pv, is called the enthalpy of the system. Now the relation between temperature of a system at a constant volume (T a ) and the physical equilibrium temperature is given by (according to [42]) The temperature at a constant volume T a also evolves as So from equations (17) and (18), Now from equations (15)and (19) one obtains, and equations (1) , (16) and (19) yield ( using the definition of enthalpy) Also equation (21) can be rewritten (using 17) as ation of volume can be expressed as (See ref. [34] and [35]) where is familiar as thermal expansibility. The isothermal compressibility (K T ) is given by Similarly the adiabatic compressibility is In isothermal process, Also one has the well established relation amongst heat capacities and compressibilities as [34] C So one can find the expression for isothermal expansibility α from equation (11) as, Now equations (27) and (28) yield Hence one obtains the expressions for compressibilities from equations (29) , (30) and (31) containing the dependence of diffusion parameter γ i as, For thermodynamic stability of any fluid , it must follow the conditions [34,35] namely For the present work, using the expressions of these thermodynamical parameters , the stability conditions are presented in the following table I : Ai ∂ ln T pi and ∂ ln 1+ which is not possible hence unstable. where, one can write equation (34) as A. Derivation of the relation between physical equilibrium temperature (T) and temperature at constant volume (T a ): From 1st law of thermodynamics (i.e. equation (8)), by choosing entropy (S) as a function of two independent variables volume (v) and temperature (T), the condition for dS to be an exact differential gives and Now, from equation (39) one obtains the energy density as a function of temperature as Now, as for an adiabatic system ∆S = 0, so one can determine from equation (41), the relation between T 0 and T a (i.e. evolution of T a ) as Now combining equations (35) , (36) and (42) one obtains , or equivalently . Now as the adiabatic contribution of entropy for the i-th fluid (excluding the diffusive entropy does not vary (i.e. d ln(S i(ad) ) = 0) , so temperature can be written as, where , L is a constant . Now combining (45) and (46) yields with Q = ω i0 L 1+ω i0 , a constant .

B. Thermodynamic derivatives
Similar to the above section , one can determine the heat capacities and compressibilities of the present thermodynamical system as follows: now the solution of equation (6),(or equivalently equation ( 34)) can be expressed as Again from equation (39), one finds Combining the above equations (48) and (49) one can write, Now imposing the initial conditions , one can formulate the modified equation of state in terms of evolution of energy with time and temperature as The above equation (51) can be written according to equation (47) as So using equations (15) , (16), (47) and (51) ,one can easily express the heat capacities as and hence one obtains the relation Again equations (48) and (49) yield which implies along with equation (47) that So from equations (27) , (53) , (55) and (57) one writes and from (28),one obtains and the expressions for K T and K S are given by (using equations (58) , (59) and (60)) and  Table II . Ai ∂ ln T pi    (Table I )and variable (Table II)  only in REGION-II,which is acceptable for expanding universe (Table III).
Further, it is interesting to note that, for constant equation of state parameter, (i.e. Table   I ), the system cannot be thermodynamically stable for −1 < ω i ≤ 0. This result is very much similar to that of [34]. However, the present work is thermodynamically stable in phantom region in contrary with non-diffusive cosmic fluids [34]. Finally, it is found that diffusive fluid with variable equation of state parameter is thermodynamically stable for all possible values of ω i under some restrictions on the diffusion parameter γ i and nature of variation of ω i with respect to temperature.