Thermodynamic implications of the gravitationally induced particle creation scenario

A rigorous thermodynamic analysis has been done at the apparent horizon of a spatially flat Friedmann-Lemaitre-Robertson-Walker universe for the gravitationally induced particle creation scenario with constant specific entropy and an arbitrary particle creation rate $\Gamma$. Assuming a perfect fluid equation of state $p=(\gamma -1)\rho$ with $\frac{2}{3} \leq \gamma \leq 2$, the first law, the generalized second law (GSL), and thermodynamic equilibrium have been studied and an expression for the total entropy (i.e., horizon entropy plus fluid entropy) has been obtained which does not contain $\Gamma$ explicitly. Moreover, a lower bound for the fluid temperature $T_f$ has also been found which is given by $T_f \geq 8\left(\frac{\frac{3\gamma}{2}-1}{\frac{2}{\gamma}-1}\right)H^2$. It has been shown that the GSL is satisfied for $\frac{\Gamma}{3H} \leq 1$. Further, when $\Gamma$ is constant, thermodynamic equilibrium is always possible for $\frac{1}{2}<\frac{\Gamma}{3H}<1$, while for $\frac{\Gamma}{3H} \leq \text{min}\left\lbrace \frac{1}{2},\frac{2\gamma -2}{3\gamma -2} \right\rbrace$ and $\frac{\Gamma}{3H} \geq 1$, equilibrium can never be attained. Thermodynamic arguments also lead us to believe that during the radiation phase, $\Gamma \leq H$. When $\Gamma$ is not a constant, thermodynamic equilibrium holds if $\ddot{H} \geq \frac{27}{4}\gamma ^2 H^3 \left(1-\frac{\Gamma}{3H}\right)^2$, however, such a condition is by no means necessary for the attainment of equilibrium.

given by T f ≥ 8 and Γ 3H ≥ 1, equilibrium can never be attained. Thermodynamic arguments also lead us to believe that during the radiation phase, Γ ≤ H. When Γ is not a constant, thermodynamic equilibrium holds ifḦ ≥ 27 4 γ 2 H 3 1 − Γ

INTRODUCTION
There have been several attempts to incorporate the present stage of cosmic acceleration into standard cosmology, the most notably being the introduction of an "exotic" matter termed dark energy (DE) which is believed to have a huge negative pressure, however, its nature and origin is still a mystery despite extensive research over the past one and a half decades. Several DE models have been proposed in the literature but observational data from various sources such as Supernovae Type Ia (SNe Ia), Cosmic Microwave Background (CMB), and Baryon Accoustic Oscillations (BAO) have established that the cosmological constant Λ is the most viable candidate among them. The cosmic concordance ΛCDM model in which the Universe is believed to contain a cosmological constant Λ associated with DE, and cold (i.e., pressureless) dark matter (abbreviated CDM) fits rather well the current astronomical data.
Nevertheless, there are severe drawbacks corresponding to a finite but incredibly small value of Λ such as the fine-tuning problem which leads to a discrepancy of 50 to 120 orders of magnitude with respect to its observed value which is about 3 × 10 −11 eV 4 . Then there is the coincidence problem which is related to the question of "why are the energy densities of pressureless matter and DE of the same order precisely at the present epoch although they evolve so differently with expansion?" Several models such as decaying vacuum models, interacting scalar field descriptions of DE, and a single fluid model with an antifriction dynamics have been proposed with a view to alleviate such problems. Moreover, in order to solve the flatness and horizon problems, an inflationary stage for the very early universe was introduced but this again gave rise to several new problems, like the initial conditions, the graceful exit, and multiverse problems.
Other attempts to explain the late time accelerating stage are modified gravity models, inhomogeneous cosmological models etc. but each one of them comes with several problems that are yet to be settled. Because of these said difficulties in various cosmological models, another well known proposal has been suggestedthe gravitationally induced particle creation mechanism. Schrodinger [1] pioneered the microscopic description of such a mechanism which was further developed by Parker and others based on quantum field theory in curved spacetimes [2][3][4][5][6]. Prigogine and collaborators [7] provided a macroscopic description of particle creation mechanism induced by the gravitational field. A covariant description was later proposed [8,9] and the physical difference between particle creation and bulk viscosity was clarified [10]. The process of particle creation is classically described by introducing a backreaction term in the Einstein field equations whose negative pressure may provide a self-sustained mechanism of cosmic acceleration. Indeed, many phenomenological particle creation models have been proposed in the literature [11][12][13][14][15][16]. It has also been shown that phenomenological particle production [17][18][19][20] can not only incorporate the late time cosmic acceleration but also provide a viable alternative to the concordance ΛCDM model.
Despite rigorous investigation of various aspects of particle creation mechanism, its thermodynamic implications have never been explored. Such a study has been undertaken in this paper and the essence of this work is that the particle creation rate has been considered arbitrary, not a phenomenological one. The conclusions drawn from the present analysis are valid for any expression of the creation rate, constant or otherwise. The paper is organized as follows. Section 2 contains a brief discussion of the gravitationally induced adiabatic particle creation scenario, Section 3 along with Subsections A, B, and C are dedicated to detailed thermodynamic analysis of the process, while Section 4 provides a short discussion and possible scope for future work.

GRAVITATIONALLY INDUCED PARTICLE CREATION MECHANISM: A BRIEF DISCUSSION
Let us consider a spatially flat, homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) universe with matter content endowed with the mechanism of particle creation. The dynamics of such a model is governed by the Friedman equations given by 1 In the above equations, ρ and p are the energy density and thermostatic pressure of the cosmic fluid respectively and they are related by the equation of state (EoS) p = (γ − 1)ρ with 2 3 ≤ γ ≤ 2, H =ȧ (t) a(t) is the Hubble parameter [a(t) is the scale factor of the Universe], and Π is the creation pressure related to the gravitationally induced process of particle creation. The lower bound on γ ensures that the perfect fluid does not become exotic, or equivalently, the strong energy condition remains valid. As a consequence, the energy conservation law gets reduced toρ Now, the non-conservation of the total number N of particles in an open thermodynamic system produces an equation given byṅ In Eqs. (3) and (4), Θ is the fluid expansion scalar which turns out to be 3H in our case, Γ denotes the rate of change of the number of particles (N = na 3 ) in a comoving volume a 3 , and n is the number density of particles. So, a positive Γ implies production of particles while a negative Γ indicates particle annihilation. Further, a non-zero Γ produces an effective bulk viscous pressure [21][22][23][24][25][26][27] of the fluid and hence non-equilibrium thermodynamics comes into the picture. Using Eqs. (3) and (4), and the Gibb's relation we can obtain an equation relating the creation pressure Π and the creation rate Γ, which can be expressed as under the customary assumption that the specific entropy s (in other words, the entropy per particle) is constant, i.e., the process is adiabatic (or isentropic). Thus, a dissipative fluid is equivalent to a perfect fluid with a nonconserved particle number. Eq. (2) now reduces tȯ The deceleration parameter q takes the form and the effective EoS parameter for this model (denoted by w ef f ) becomes which represents quintessence era for Γ < 3H and phantom era for Γ > 3H, while Γ = 3H corresponds to a cosmological constant, owing to the fact that w ef f = −1.

THERMODYNAMIC ANALYSIS
In the following subsections, we shall study the first law, the generalized second law (GSL) 2 , and thermodynamic equilibrium for an arbitrary particle creation rate Γ. We shall consider an apparent horizon as our thermodynamic boundary, since, unlike the event horizon, a cosmic apparent horizon always exists and it coincides with the event horizon in the case of a last de Sitter expansion. Moreover, in a flat FLRW universe, the apparent horizon coincides with the Hubble horizon H −1 . So, the apparent horizon can be considered to be located at R A = 1 H and its first order derivative with respect to the cosmic time t can be evaluated aṡ The (Bekenstein) entropy and (Hawking) temperature of the apparent horizon are given by and respectively.

A. First law
The first law of thermodynamics at the horizon is governed by the Clausius relation The differential dE A of the amount of energy crossing the apparent horizon can be evaluated as (see Eq. (27) of Ref. [29]) Again, using the expressions of T A and S A given in Eqs. (11) and (12), the expression T A dS A becomes where we have used relation (10). From the above analysis, we find that the first law holds at the apparent horizon whenever Γ = 0, or loosely speaking, whenever Γ ≪ 3H.

B. Generalized second law: An expression for total entropy
According to thermodynamics, the equilibrium configuration of an isolated macroscopic physical system should be the maximum entropy state, consistent with the constraints imposed on the system. Thus if S is the total entropy of the system, the following conditions should hold -(a) dS ≥ 0 [i.e., the entropy function cannot decrease (the second law of thermodynamics)], and (b) d 2 S < 0 [i.e., the entropy function attains a maximum (thermodynamic equilibrium)]. In our context, the Universe bounded by an apparent horizon and filled with some cosmic fluid forms an isolated macroscopic physical system for which the above inequalities can be generalized as respectively, where S f is the entropy of the cosmic fluid contained within the horizon. The inequality (i) is sometimes called the GSL.

The Gibb's equation can be rewritten in the form
where T f is the temperature of the cosmic fluid respectively, and E f = ρV is the energy of the fluid. Now, the assumption of a constant specific entropy leads us to an evolution equation for the fluid temperature given by (see the second relation in Eq. (35) of Ref. [30]) Noting from Eq.
On integration, we obtain, where T 0 is the constant of integration. Note that Γ does not appear explicitly in the equation. From Eq. (17), the differential of the fluid entropy can be obtained in the following form: The differential of the total entropy can then be evaluated as It can be easily seen from the previous equation that GSL holds if Γ ≤ 3H, or equivalently, if Γ 3H ≤ 1. Therefore, the GSL is not consistent with the phantom fluid. Furthermore, S T is a constant of motion when Γ = 3H, i.e., when w ef f = −1, a cosmological constant.
Another remarkable fact is that Eq. (21) gives us an opportunity [by replacing dt by dḢ H and using Eq. (7)] to derive an expression for the total entropy in terms of the Hubble parameter H in the form The essence of Eq. (22) lies in the fact that the particle creation rate Γ does not occur explicitly in the equation.
Requiring that the total entropy be always positive, we can, in principle, obtain a lower bound on T 0 given by Eq. (23) implies that we can also impose a lower bound on the fluid temperature T f as For radiation era (i.e., γ = 4 3 ) and matter dominated era (i.e., γ = 1), the lower bounds on T f become T f ≥ 16H 2 and T f ≥ 4H 2 respectively. Using Maple software, the total entropy S T has been plotted against γ for H = 67 and T 0 = 10 5 , and presented in Figure 1.

C. Thermodynamic equilibrium
Case I: Γ is constant -If the particle creation rate Γ is assumed to be constant, then the second order differential of the total entropy can be obtained from Eq. (21) as In Table I, we have explored relevant subintervals of Γ in order to test for the validity of thermodynamic equilibrium. From the table, it is evident that thermodynamic equilibrium holds unconditionally for 1 2 < Γ 3H < 1, while it never holds for Γ 3H ≤ min 1 2 , 2γ−2 3γ−2 and Γ 3H ≥ 1. Thus, thermodynamic equilibrium in this case is inconsistent with the cosmological constant as well as the phantom fluid.
From different observational sources, it has been well established that the radiation phase was followed by a matter dominated era which eventually transited to a second de Sitter phase. Accordingly, it can be expected that in the radiation dominated era, the entropy increased but the thermodynamic equilibrium was not achieved [31]. If this were not true, the Universe would have attained a state of maximum entropy and would have stayed in it forever unless acted upon by some "external agent." However, it is a well known fact [6] that the production of particles was suppressed during the radiation phase, so in this model, there would be no external agent to remove the system from thermodynamic equilibrium. Therefore, our present analysis leads us to conclude that during the radiation phase, if Γ is constant, then Γ 3H ≤ 1 3 , or equivalently, Γ ≤ H.
where we have substituted the value ofΓ evaluated aṡ It is evident from Eq. (26) that it is quite difficult to perform an analysis similar to the one that we have done in the previous case. The only definite conclusion which can be made here is that the thermodynamic equilibrium holds ifḦ ≥ 27

DISCUSSION AND FUTURE WORK
This paper dealt with a rigorous thermodynamic analysis at the apparent horizon of a spatially flat FLRW universe for the gravitationally induced particle creation scenario with constant specific entropy and an arbitrary particle creation rate Γ. Assuming a perfect fluid EoS p = (γ − 1)ρ with 2 3 ≤ γ ≤ 2 (the lower bound ensures that the strong energy condition remains valid), the first law, the GSL, and thermodynamic equilibrium have been studied and the following results have been found-• The first law holds at the apparent horizon either for a zero particle creation rate or, loosely speaking, when the creation rate is infinitesimally small as compared to 3H.
• The GSL holds if Γ ≤ 3H, or equivalently, if Γ 3H ≤ 1, which implies that the GSL is not consistent with the phantom fluid.
• For a constant particle creation rate, thermodynamic equilibrium always holds for 1 2 < Γ 3H < 1, while it never holds for Γ 3H ≤ min 1 2 , 2γ−2 3γ−2 and Γ 3H ≥ 1. Thus, thermodynamic equilibrium in this case is inconsistent with the cosmological constant as well as the phantom fluid.
• When Γ is not constant, the only definite conclusion which can be made is that the thermodynamic equilibrium holds ifḦ ≥ 27 4 γ 2 H 3 1 − Γ 3H 2 , however, such a condition is by no means necessary for the attainment of equilibrium.
An expression for the total entropy with no explicit dependence on Γ has also been found. Such an expression suggests that for Γ = 3H, i.e., a cosmological constant, the total entropy is a constant of motion. Further, imposing the condition that the total entropy is always positive, a lower bound on the fluid temperature T f has been obtained. It is evident that T f ≥ 16H 2 and T f ≥ 4H 2 for radiation and matter dominated eras respectively. Thermodynamic arguments also lead us to believe that if Γ is a constant, then Γ ≤ H during the radiation phase.
For future work, thermodynamics of the particle creation scenario at any arbitrary horizon can be investigated. The present thermodynamic analysis can also help to constrain various parameters of phenomenological particle creation rates that have been considered in recent literature [18,19,[31][32][33][34][35][36][37][38]. Further, attempts to incorporate matter creation in inhomogeneous cosmological models can be made and its thermodynamic implications can be studied.