Thermodynamic implications of the gravitationally induced particle creation scenario

A rigorous thermodynamic analysis has been done as regards 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 Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document}. Assuming a perfect fluid equation of state p=(γ-1)ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p=(\gamma -1)\rho $$\end{document} with 23≤γ≤2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{2}{3} \le \gamma \le 2$$\end{document}, 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 Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} explicitly. Moreover, a lower bound for the fluid temperature Tf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_f$$\end{document} has also been found which is given by Tf≥83γ2-12γ-1H2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_f \ge 8\left( \frac{\frac{3\gamma }{2}-1}{\frac{2}{\gamma }-1}\right) H^2$$\end{document}. It has been shown that the GSL is satisfied for Γ3H≤1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\Gamma }{3H} \le 1$$\end{document}. Further, when Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} is constant, thermodynamic equilibrium is always possible for 12<Γ3H<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{2}<\frac{\Gamma }{3H} < 1$$\end{document}, while for Γ3H≤min12,2γ-23γ-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\Gamma }{3H} \le \text {min}\left\{ \frac{1}{2},\frac{2\gamma -2}{3\gamma -2} \right\} $$\end{document} and Γ3H≥1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\Gamma }{3H} \ge 1$$\end{document}, equilibrium can never be attained. Thermodynamic arguments also lead us to believe that during the radiation phase, Γ≤H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma \le H$$\end{document}. When Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} is not a constant, thermodynamic equilibrium holds if H¨≥274γ2H31-Γ3H2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ddot{H} \ge \frac{27}{4}\gamma ^2 H^3 \left( 1-\frac{\Gamma }{3H}\right) ^2$$\end{document}, however, such a condition is by no means necessary for the attainment of equilibrium.


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" component 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 a e-mail: subhajit1729@gmail.com literature but observational data from various sources such as supernovae type Ia (SNe Ia), cosmic microwave background (CMB), and baryon acoustic 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-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 suggested-the gravitationally induced particle creation mechanism. Schrödinger [1] pioneered the microscopic description of such a mechanism which was further developed by Parker et al. based on quantum field theory in curved spacetimes [2][3][4][5][6]. Prigogine et al. [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] cannot 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, Sect. 3 along with subsections A, B, and C is dedicated to a detailed thermodynamic analysis of the process, while Sect. 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 ( 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 1 In this manuscript, without any loss of generality, we have assumed that the physical constants, namely c, G,h, and κ B , as well as 8π , are unity.ρ The non-conservation of the total number N of particles in an open thermodynamic system produces an equation given bẏ 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-27] of the fluid and hence non-equilibrium thermodynamics comes into the picture.
Using Eqs. (3) and (4), and the Gibbs 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 non-conserved particle number. Equation (2) now reduces tȯ The deceleration parameter q takes the form and the effective EoS parameter for this model (denoted by w e f f ) becomes which represents the quintessence era for < 3H and the phantom era for > 3H , while = 3H corresponds to a cosmological constant, owing to the fact that w eff = −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 final 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.

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]) 2 The idea of incorporating the GSL in cosmology was first developed by Brustein [28]. This second law is based on the conjecture that causal boundaries and not only event horizons have geometric entropies proportional to their area.
Again, using the expressions of T A and S A given in Eqs. (11) and (12)

, the expression T A d S A becomes
where we have used Eq. (10). From the above analysis, we find that the first law holds at the apparent horizon whenever = 0, or loosely speaking, whenever 3H .

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 Gibbs 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.
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: (20) The differential of the total entropy can then be evaluated as It can easily be 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 eff = −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 Equation (23) implies that we can also impose a lower bound on the fluid temperature T f as For the 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 in Fig. 1.

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 1, 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 where we have substituted the value of˙ evaluated aṡ It is evident from Eq. (26)

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 drawn 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 the 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 the 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.