Thermodynamics of Gravitationally Induced Particle Creation Scenario in DGP Braneworld

In this paper, we discuss the thermodynamical analysis for gravitationally induced particle creation scenario in the framework of DGP braneworld model. For this purpose, we consider apparent horizon as the boundary of the universe. We take three types of entropy such as Bakenstein entropy, logarithmic corrected entropy and power law corrected entropy with ordinary creation rate $\Gamma$. We analyze the first law and generalized second law of thermodynamics analytically for these entropies which hold under some constraints. The behavior of total entropy in each case is also discussed which implies the validity of generalized second law of thermodynamics. Also, we check the thermodynamical equilibrium condition for two phases of creation rate, that is constant and variable $\Gamma$ and found its validity in all cases of entropy.


Introduction
It is believed that the universe undergoes an accelerated expansion due to mysterious form of force called dark energy (DE) which was firstly predicted by two independent teams of cosmologists [1,2]. Both of them used distant type Ia supernova as standard candles to measure the expansion of the universe. This discovery was unexpected, because before this invention, cosmologists just think that the expansion of the universe would be decelerating because of the gravitational attraction of the matter in the universe. In the accelerated expansion of the universe, DE plays major role but its nature is still unknown. The simplest candidate of DE is the cosmological constant, but its composition and mechanism are unknown. More generally, the detail of its equation of state (EoS) and relationship with the standard model of particle physics continue to be investigated both through observationally and theoretically [3]. In order to explain this cosmic acceleration, various DE models and modified theories of gravity have been developed such as f (R) theory, f (T ) theory, Brans-Dicke theory, dynamical Chern-Simons modified gravity, DGP Braneworld model, etc.
The concept of thermodynamics in cosmological system originates through black hole physics. It was suggested [4] that the temperature of Hawking radiations emitting from black holes is proportional to their corresponding surface gravity on the event horizon. Jacobson [5] found a relation between thermodynamics and the Einstein field equations. He derived this relation on the basis of entropy-horizon area proportionality relation along with first law of thermodynamics (also called Clausius relation) dQ = T dS, where dQ, T and dS indicate the exchange in energy, temperature and entropy change for a given system. It was shown that the field equations for any spherically symmetric spacetime can be expressed as T dS = dE + P dV (E, P and V represent the internal energy, pressure and volume of the spherical system) for any horizon [6].
The generalized second law of thermodynamics (GSLT) has been studied extensively in the scenario of expanding behavior of the universe. The GSLT states that the entropy of matter inside the horizon plus entropy of the horizon remains positive and increases with the passage of time [7]. In order to discuss GSLT, horizon entropy of the universe can be taken as one quarter of its horizon area [8] or power law corrected [9] or logarithmic corrected [10] forms. Many people have explored the validity of GSLT of different systems including interaction of two fluid components like DE and dark matter [11], as well as interaction of three components of fluid [12] in the FRW universe by using simple horizon entropy of the universe. The thermodynamical analysis widely performed in modified theories of gravity [13].
The gravitationally induced particle creation is another well-known mechanism which was firstly introduced by Schrodinger [14] on microscopic level. This mechanism was further extended by Parker et al. towards quantum field theory in curved spacetimes [15]- [19]. The macroscopic description of particle creation mechanism induced by gravitational field was presented by Prigogine et al. [20]. Later on, the covariant description of this mechanism was developed [21,22] as well as the physical difference between particle creation and bulk viscosity was also given [23]. The particle creation process can be described with the inclusion of backreaction term in the Einstein field equations whose negative pressure may help in explaining the cosmic acceleration. In this way, various phenomenological models of particle creation have been presented [24]- [29]. It is also shown that phenomenological particle production [30]- [33] help in explaining the cosmic acceleration and paved the alternative way to the concordance ΛCDM model.
The fact that we resides in a three-dimensional space embedded in an extra-dimensional world and five-dimensional models in which universe would be a hypersurface has attain a great attention. The four-dimensional Einsteins equations projected onto the brane have been explored by Shiromizu et al. [34]. The approaches which made on the basis of brane-world in the early-time cosmology favor a particular model of cosmic evolution featured by quadratic relations between the energy density and the Hubble parameter, dubbed quadratic cosmology [35].
Recently, by assuming the gravitationally induced particle scenario with constant specific entropy and arbitrary particle creation rate (Γ), thermodynamics on the apparent horizon for FRW universe has been discussed [36]. They have investigated the first law, GSLT and thermodynamical equilibrium by assuming the EoS for perfect fluid and put forward various constraints on Γ for which thermodynamical laws hold. Our aim is to discuss the thermodynamical analysis on the apparent horizon for gravitationally induced particle creation scenario with ordinary creation rate Γ by assuming entropies (Bakenstein entropy or usual entropy, logarithmic corrected entropy and power law corrected entropy) in DGP braneworld model. The scheme of the paper is as follows: In the next section, we will present the basic equations of DGP brane-world, particle creation rate and cosmological parameters. Section 3, 4, 5 contain the discussion of thermodynamic quantities as well as its laws corresponding to usual, logarithmic and power law corrected entropies, respectively. The last section comprises of concluding remarks on our results.

Basic Equations
A most particular version was proposed by Dvali et al. [37], in which the four dimensional FRW universe is enclosed in a five dimensional Minkoski bulk with infinite size. The gravitational laws were obtained by adding an Einstein-Hilbert term to the action of brane computed with the brane's intrinsic curvature. The presence of such a term in the action is generically induced by quantum corrections coming from the bulk gravity and its coupling with matter living on the brane and must be included for a large class of theories for self-consistency [38,39]. Here, we consider 3-brane embedded in a 5D space-time with an intrinsic curvature term included in the brane whose action can be written as where L m is the brane curvature term, given by and κ 2 = 8πG (5) = M −3 (5) , µ 2 = 8πG (4) = M −2 (4) . The Eq.(1) represents the Einstein-Hilbert action in five dimensions for a five-dimensional metricg AB (bulk metric) of Ricci scalarR. Similarly, Eq.(2) indicates the Einstein-Hilbert action for the induced metricg cd on the brane with R appeared as its scalar curvature. From Eq.(1), we can get modified Friedmann equation as [40] where H =˙a a is Hubble parameter with a(t) is the scale factor. Also, ρ = ρ m + ρ D , the subscripts m and D represent the energy densities corresponding to dark matter and DE respectively, r c is the crossover length which represents the scale that has length away from which gravity starts opening into the bulk [40]. Moreover it is the distance scale follow the comparison among 4D and 5D effects of gravity and can be written as [40] where M 2 5 stands for the 5D Planck mass and M 2 p is the 4D Planck mass. For the spatially flat DGP braneworld (k = 0), Eq.(3) reduces to There exist two different branches for the DGP model depending on the sign of ǫ. These are as follows: • For ǫ = +1, there is a de Sitter solution for Eq.(5) with constant Hubble parameter, i.e., H = 1 rc ⇒ a(t) = a 0 e t rc in the absence of any kind of energy or matter field on the brane (i.e., ρ = 0). However, this branch faces some problems like ghost instabilities [41].
• For ǫ = −1, the accelerated expansion of the universe can only be explained through the inclusion of DE component in the DGP scenario.
We consider the latter case in the present work. The equation of continuity for this model will becomeρ where Π is a particle creation pressure which represents the gravitationally induced process of particle creation and Θ = 3H is the fluid expansion. Differentiating Eq.(5) and replacing the value ofρ using Eq.(6), we obtaiṅ The respective EoS for this model is giving by p = γ − 1 ρ with 2 3 ≤ γ ≤ 2. The non-conservation of the total rate of change of number of particles, N = na 3 with comoving volume a 3 and n is the number density of particle production in an open thermodynamical system yieldṡ where Γ is a particle creation rate has negative and positive phases. Negative Γ represents the particle destruction and positive Γ describes the elimination of particles. Furthermore, a non-zero Γ produces effective bulk viscous pressure [42]- [48]. Now using the Eqs. (6), (8) and Gibbs relation, we get An equation related to the creation pressure Π and the creation rate Γ has the form Under traditional assumption that the specific entropy of each particle is constant, i.e., the process is adiabatic or isentropic. This implies a dissipative fluid is similar to a perfect fluid with a non-conserved particle number. To discuss cosmological parameters, we insert Eqs.(10) and p = (γ − 1)ρ in Eq. (7), to obtainḢ The deceleration parameter q can be written as The effective EoS parameter for this model turns out to be This parameter has ability to explain the different phases of the universe on the basis of Γ, i.e., if Γ < 3H, then we have quintessence era of the universe (ω eff < −1), if Γ > 3H then effective EoS parameter represents the phantom era of the universe (ω eff > −1) while for Γ = 3H, effective EoS parameter exhbits the cosmological constant (ω eff = −1).
In the following, we will discuss first and second thermodynamical laws in the presence of particle creation rate Γ on the apparent horizon.

Thermodynamical Analysis with Usual Entropy
For flat FRW universe, Hubble parameter coincides with the apparent horizon as R A = 1 H . Differentiating the apparent horizon with respect to time, we getṘ The Bekenstein entropy and Hawking temperature of the apparent horizon are given by (8π = G = 1) where A = 4πR 2 A . The first law of thermodynamics at the horizon can be obtained through the Clausius relation as For the sake of convenance, we consider Ξ = T A dS A + dE A . The differential dE A is the amount of energy crossing the apparent horizon can be evaluated as [49] − From Eq.(15), the differential of surface entropy at apparent horizon yields which leads to Thus, Ξ turns out to be From this relation, it can be seen that first law of thermodynamics holds (i.e., Ξ → 0) at the apparent horizon for Γ = 3H 1 − 2H− ǫ rc 2H . Next, we will discuss the GSLT and thermodynamical equilibrium of a system containing perfect fluid distribution bounded by apparent horizon in DGG brane-world scenario. For GSLT, the total entropy of the system can not be decrease, i.e., d(S A +S f ) ≥ 0. In this relation, S A and S f appear as the entropy at apparent horizon and the entropy of cosmic fluid enclosed within the horizon, respectively. The Gibbs equation is given by where T f is the temperature of the cosmic fluid and E f is the energy of the fluid (E f = ρV ). The evolution equation for fluid temperature having constant entropy can be described as [50] T where T 0 is the constant of integration. The differential of the fluid entropy can be obtained by using the Eq. (21) as follows Using Eqs. (18) and (24), we get the rate of change of total entropy aṡ where S T = S A + S f . We discuss the validity of GSLT on the basis of Γ such that • Γ < 3H: The GSLT holds if the following constraint satisfies. This shows that the GSLT satisfies in the qunitessence era of the evolving universe.
• Γ > 3H: For this case, we have the constraint which implies the GSLT holds in phantom era of the universe.
Replacing dt to dḢ H and integrating Eq.(25), we get The plot between total entropy and parameter γ is shown in Figure 1 for three values of T by setting constant values as H = 67, r c = 1 67 , ǫ = −1. We observe that S T ≥ 0 for all the values of T which leads to the validity of GSLT.

Thermal Equilibrium Scenario
Further, we will discuss the thermal equilibrium scenario in the present case. For thermodynamical equilibrium, the entropy function attains a maximum value and satisfies the condition d 2 S T = d 2 (S A + S f ) < 0. For this purpose, we consider two cases of particle creation rate (Γ).

Case 1: Γ = constant
Firstly, we consider particle creation rate Γ as a constant. Under this scenario, differentiating Eq.(25) w.r.t time, it results the following second order differential equation where λ = 1 + 4H . The plot betweenS T versus γ for three values of T with constant values of H = 67, r c = 1 67 . ǫ = −1, q = −0.53 as shown in Figure 2. One can observe that the thermodynamical equilibrium condition holds for all values of T with specific ranges of γ. For example, for T = 10 2 , thermodynamical equilibrium holds for the range 1.4 < γ ≤ 1.8 and does not obey for 0.6 ≤ γ ≤ 1.4. For T = 10 2.3 , thermal equilibrium holds for the range 1.3 < γ ≤ 1.8 and does not showing the validity for 0.6 ≤ γ ≤ 1.3. However, for T = 10 2.5 , thermodynamic equilibrium condition holds for the range 1.2 < γ ≤ 1.8 and disobey for the range 0.6 ≤ γ < 1.2.

Case 2: Γ = Γ(t)
Here we take Γ as variable parameter, i.e., Γ = Γ(t), for which Eq.(25) becomes The plot ofS T versus γ for three values of T as shown in Figure 3

Logarithmic Corrected Entropy
Quantum gravity allows the logarithmic corrections in the presence of thermal equilibrium fluctuations and quantum fluctuations [51]- [57]. The loga- rithmic entropy corrections can be defined as where α and β are constants whose values are still under consideration. The differential form of above equation leads to which gives In view of this entropy, the quantity Ξ takes the form It can be observed from Eq.(32) that the first law of thermodynamics holds To discuss the GSLT for logarithmic corrected entropy of horizon, we obtain the total entropy by using Eqs. (24) and (30) as followṡ The GSLT will hold under these constraints.
• For the case Γ < 3H, the GSLT satisfy in the quintessence region of the universe if the following constraint holds • Γ > 3H For this case, we obtain the following constraint which implies the GSLT holds in phantom era of the universe.
The expression of total entropy in the form of Hubble parameter is given by The plot of total entropy S T versus γ with respect to three values of T is shown in Figure 4 with constant values as α = −2, β = −0.00001, L p = 1. It is observed that the total entropy is positive, i.e, S T > 0 which leads to the validity of GSLT for all values of T .

Thermal Equilibrium Scenario
Now we will discuss the thermodynamic equilibrium by assuming two cases for particle creation rate Γ as follows: In this way, the second order differential equation can be obtained from Eq.(33) for Γ as a constanẗ where The plot betweenS T and γ for three values of T by fixing the constant values q = −0.53, α = −2, β = −0.000000001 and L p = 1 as shown in Figure 5. It can be seen that thermodynamical equilibrium is obeying the condition S T < 0 for all values of T which leads to the thermal equilibrium.

Case 2: Γ = Γ(t)
Taking Γ as a function of t, Eq.(33) yields Figure 6 reperesents the plot betweenS T and γ for three values of T for variable Γ for same constant values. Figure 6 indicates that the trajectories ofS T

Power Law Corrected Entropy
The power-law correction to the entropy-area law comes from association of the wave-function of the scalar field between the ground state and the exited state [58]- [60]. The correction term is also more significant for higher excitations. It is important to note that the correction term decreases faster with A and hence in the semi-classical limit (large area) the entropy-area law is recovered. The power entropy can b given as where δ is dimensionless constant and r c is the crossover scale. From Eq.(37) the differential of surface entropy at horizon can be expressed as which gives In this entropy correction, we have 1 The plot of total entropy (S T ) is shown in Figure 7 versus γ for three values of T by assigning the constant values as L = 1, H = 67, 1 67 , α = −2, and ǫ = −1. We observe that total entropy is positive, i.e, (S T ) > 0 which leads to the validity of GSLT for all values of T .

Thermal Equilibrium Scenario
For thermodynamic equilibrium, we assume two cases of Γ, i.e., Γ is constant and Γ is variable.

Case 1: Γ = constant
For constant Γ, the second order differential equation of total entropy leads toS where The graphical behavior ofS T against γ is shown in Figure 8 for three values of T by keeping the constant values as H = 67, r c = 1 67 , ǫ = −1, q = −0.53, δ = −2, L = 1. It is found that the trajectories of this plot satisfy the conditionS T < 0 for all values of T which leads to the thermodynamical equilibrium of the system.

Case 2: Γ = Γ(t)
For variable Γ, second order differential equation of total entropy turns out to bë The plot betweenS T and γ for three values of T by setting the constant values as H = 67, r c = 1 67 , ǫ = −1, q = −0.53, δ = −2, L = 1 as shown in Figure  9. We observe that thermodynamical equilibrium satisfying the condition S T < 0 for all values of T which leads to thermodynamical disequilibrium of the system.

Comparison and Conclusion
We have studied the thermodynamics on the apparent horizon for gravitationally induced particle creation scenario by assuming entropy corrections (Bakenstein entropy, logarithmic corrected entropy and power law corrected entropy) in DGP braneworld gravity. Considering the perfect fluid EoS p = (γ − 1)ρ, we have analyzed the first law of thermodynamics, GSLT and thermodynamic equilibrium. We have fixed model parameters in the way to obtain the validity of GSLT. However, the results may be changed for other model parameters. The results have been summarized as follows: For Bakenstein entropy, the first law of thermodynamics holds at appar- . The GSLT valid under some constraints in different eras of the evolving universe. In terms of total entropy scenario, GSLT remains valid for all values of T with 2 3 ≤ γ ≤ 2 ( Figure 1). In the presence of usual entropy of horizon and for constant Γ, thermodynamical equilibrium holds for all values of T with specific ranges of γ. For example, thermodynamical equilibrium is satisfying the condition d 2 S T dt 2 < 0 for 1.4 < γ ≤ 1.8, 3 < γ ≤ 1.8, 1.2 < γ ≤ 1.8 and does not showing the validity 0.6 ≤ γ ≤ 1.4, 0.6 ≤ γ ≤ 1.3 and 0.6 ≤ γ < 1.2 with T = 10 2 , T = 10 2.3 and T = 10 2.5 respectively (Figure 2). For variable Γ, Figure 3 indicates that thermodynamical equilibrium holds for all values of T with different ranges of γ such as 1.4 < γ ≤ 1.8, 1.3 < γ ≤ 1.8, 1.2 ≤ γ ≤ 1.8 and does not obey 0.6 ≤ γ ≤ 1.4, 0.6 ≤ γ ≤ 1.3 and 0.6 ≤ γ < 1.2 for T = 10 2 , T = 10 2.3 and T = 10 2.5 respectively.
For logarithmic corrected entropy, it has observed that first law of thermodynamics holds at apparent horizon when Γ = 3H 1− . The GSLT valid under some constraints in different eras of the evolving universe. In terms of total entropy scenario, GSLT is valid for all values of T with 2 3 ≤ γ ≤ 2 ( Figure 4). Thermodynamical equilibrium condition ( d 2 S T dt 2 < 0) satisfied for all values of T with all values of γ for constant as well as variable Γ ( Figure 5 and 6).
For power law entropy corrected, it has been observed that the first law of thermodynamics satisfies under constraint Γ = 3H 1 − and GSLT holds under some conditions on model parameters in quintessence, phantom and cosmological constant phases. Also, GSLT meets for all values of T in the range ( 2 3 ≤ γ ≤ 2) (Figure 7). It is found that, for constant as well as variable Γ, the thermodynamic equilibrium condition d 2 S T dt 2 < 0 obey for all values of T which leads to the thermodynamical equilibrium of the system.
Here we provide some details about past works and compare with underlying work. Harko et al. [61] considered the possibility of a gravitationally induced particle production through the mechanism of a non-minimal curvaturematter coupling. An interesting feature of this gravitational theory is that the divergence of the energy momentum tensor is nonzero. Firstly, they have reformulated the model in terms of an equivalent scalartensor theory, with two arbitrary potentials. By using the formalism of open thermodynamic systems, they have interpreted the energy balance equations in this gravitational theory from a thermodynamic point of view, as describing irre-versible matter creation processes. The particle number creation rates, the creation pressure, and the entropy production rates have explicitly obtained as functions of the scalar field and its potentials, as well as of the matter Lagrangian. The temperature evolution laws of the newly created particles are also obtained. The cosmological implications of the model have briefly investigated, and it is shown that the late-time cosmic acceleration may be due to particle creation processes. Furthermore, it has also shown that due to the curvaturematter coupling, during the cosmological evolution a large amount of comoving entropy is also produced.
Mitra et al. [62] have studied thermodynamics laws by assuming flat FRW universe enveloped by by apparent and event horizon in the framework of RSII brane model and DGP brane scenario. Assuming extended Hawking temperature on the horizon, the unified first law is examined for perfect fluid (with constant equation of state) and Modified Chaplygin Gas model. As a result there is a modification of Bekenstein entropy on the horizons. Further the validity of the generalized second law of thermodynamics and thermodynamical equilibrium are also investigated.
Pan et al. [63] investigated the expansion of the universe powered by the gravitationally induced adiabatic matter creation by developing general creation rate and their dynamical analysis. They also developed dynamical analysis in the the presence of a non-singular universe (without the big bang singularity) with two successive accelerated phases, one at the very early phase of the universe (i.e. inflation), and the other one describes the current accelerating universe, where this early, late accelerated phases are associated with an unstable fixed point (i.e. repeller) and a stable fixed point (attractor), respectively.
Sal and Haro [64] performed a qualitative and thermodynamic study of two models when one takes into account adiabatic particle production. In the first one, there is a constant particle production rate, which leads to solutions depicting the current cosmic acceleration but without inflation. The other one has solutions that unify the early and late time acceleration. These solutions converge asymptotically to the thermal equilibrium.
Recently, by assuming the gravitationally induced particle scenario with constant specific entropy and arbitrary particle creation rate (Γ), thermodynamics on the apparent horizon for FRW universe has been discussed [36]. They have investigated the first law, GSLT and thermodynamical equilibrium by assuming the EoS for perfect fluid and put forward various constraints on Γ for which thermodynamical laws hold. We have extended the work of [36] in the DGP brane-world scenario by assuming usual entropy as well as its entropy corrections (power law as well as logarithmic corrected) in flat FRW universe. We have extracted EoS parameter and obtained its various constraints with respect to quintessence, vacuum and phantom era of the universe. For variable as well as constant particle creation rate (Γ), the first law of thermodynamics, GSLT and thermal equilibrium condition is satisfied in all cases of entropies forms.