Gravitational induced particle production through a nonminimal curvature-matter coupling

We consider the possibility of a gravitationally induced particle production through the mechanism of a nonminimal curvature-matter coupling. An interesting feature of this gravitational theory is that the divergence of the energy-momentum tensor is nonzero. As a first step in our study we reformulate the model in terms of an equivalent scalar-tensor theory, with two arbitrary potentials. By using the formalism of open thermodynamic systems, we interpret the energy balance equations in this gravitational theory from a thermodynamic point of view, as describing irreversible matter creation processes. The particle number creation rates, the creation pressure, and the entropy production rates are 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 are briefly investigated, and it is shown that the late-time cosmic acceleration may be due to particle creation processes. Furthermore, it is also shown that due to the curvature--matter coupling, during the cosmological evolution a large amount of comoving entropy is also produced.


I. INTRODUCTION
The simplest explanation of the late-time accelerated expansion of the Universe [1,2] is to invoke a cosmological constant, Λ, which can be associated to the vacuum energy [3]. Despite the excellent fit to observational data, the presence of Λ suffers from two serious drawbacks, namely, the cosmological constant problem and the coincidence problem. However, it has recently been argued that the vacuum energy is not constant but decays into other particle constituents [4]. Indeed, phenomenological models, with a variable cosmological constant [5], have been proposed to address the above problems. For instance, a simple and thermodynamically consistent cosmology with a phenomenological model of quantum creation of radiation due to vacuum decay was presented in [6], where the thermodynamics and Einstein's equations lead to an equation in which H is determined by the particle number N . It was shown that the evolution equation for H has a remarkably simple exact solution, in which a non-adiabatic inflationary era exits smoothly to the radiation era, without a reheating transition. In [7], a new accelerating flat model without dark energy that is fully dominated by cold dark matter (CDM) was investigated. It was shown that the number of CDM par- * Electronic address: t.harko@ucl.ac.uk † Electronic address: fslobo@fc.ul.pt ‡ Electronic address: jpmimoso@fc.ul.pt § Electronic address: diego.pavon@uab.es ticles is not conserved and the present accelerating stage is a consequence of the negative pressure describing the irreversible process of gravitational particle creation.
In [4], the correspondence between cosmological models powered by a decaying vacuum energy density and gravitationally induced particle production was explored. Although being physically different, it was shown that under certain conditions both classes of cosmologies can exhibit the same dynamical and thermodynamical behavior. By using current type Ia supernovae data, recent estimates of the cosmic microwave background shift parameter and baryon acoustic oscillations measurements, the authors performed a statistical analysis to test the observational viability of the models and the best-fit of the free parameters was also obtained. Furthermore, the particle production cosmologies (and the associated decaying Λ(t)-models) were modelled in the framework of field theory by a phenomenological scalar field model.
In this context, a new cosmic scenario with gravitationally induced particle creation was proposed [8], where the Universe evolves from an early to a late time de Sitter era, with the recent accelerating phase driven only by the negative creation pressure associated with the cold dark matter component. The model can be interpreted as an attempt to reduce the so-called cosmic sector (dark matter plus dark energy) and relate the two cosmic accelerating phases (early and late time de Sitter expansions). A detailed thermodynamic analysis including possible quantum corrections was also carried out. For a very wide range of the free parameters, it was found that the model presents the expected behavior of an ordinary macroscopic system in the sense that it approaches ther-modynamic equilibrium in the long run (i.e., as it nears the second de Sitter phase). Moreover, an upper bound was found for the Gibbons-Hawking temperature of the primordial de Sitter phase [9]. Finally, when confronted with the recent observational data, the current 'quasi'-de Sitter era, as predicted by the model, it was verified to pass the cosmic background tests very comfortably.
In this work, we consider an alternative mechanism for the gravitational particle production, namely, through a nonminimal curvature-matter coupling, in modified theories of gravity. A general property of these theories [10][11][12][13][14][15], is the non-conservation of the energy-momentum tensor (for a recent review of modified gravity models with curvature-matter coupling see [16]). Thus, the coupling between the matter and the higher derivative curvature terms may be interpreted as an exchange of energy and momentum between both. This latter mechanism induces a gravitational particle production. We note that the generalized energy balance equations in these gravitational theories have been interpreted from a thermodynamic point of view as describing irreversible matter creation processes in [17]. Thus, the coupling between matter and geometry generates an irreversible energy flow from the gravitational field to newly created matter constituents, with the second law of thermodynamics requiring that the geometric curvature transforms into matter.
Here we extend and refine the analysis initiated in [17] by investigating in detail the thermodynamic interpretation of the curvature-matter gravitational coupling for the so-called linear version of the f (R, L m ) gravity theory, where R is the Ricci scalar and L m is the matter Lagrangian, which we denote as Lf (R, L m ) theory, with the gravitational Lagrangian given by f 1 (R)/2 + [1 + λf 2 (R)] L m , where f 1 (R) and f 2 (R) are arbitrary functions of R, and λ is a coupling constant. We note that in the Lf (R, L m ) theory, the gravitational action is linear in the matter Lagrangian L m , and not in the Ricci scalar R. As a first step in our study we introduce the equivalent scalar-tensor description of the theory [18], in which the action is equivalent to a Brans-Dicke type theory, with a single scalar field ψ, a vanishing Brans-Dicke parameter ω, and a coupling U (ψ) between the scalar field and matter.
By using the formalism of open thermodynamic systems [19][20][21][22], we interpret the energy balance equation of the theory as describing a matter creation process. Indeed, the irreversible thermodynamics of open systems, and its implications for cosmology have been extensively analyzed [23]. Here, we obtain the equivalent particle number creation rates, the creation pressure and the entropy production rates as functions of the scalar field, of the two scalar potentials, and of the matter Lagrangian, respectively. The temperature evolution of the newly created particles is also obtained. Due to the curvaturematter coupling, during the cosmological evolution a large amount of comoving entropy could be produced. The cosmological implications of the theory in its scalar-tensor representation are also investigated.
The present paper is organized as follows. In Section II, we present the action and the field equations of the modified gravity model with a linear curvature-matter coupling in both their standard and scalar-tensor representations. The thermodynamic interpretation of the theory is developed in Section III, where the particle creation rates, the creation pressure and the entropy production are analyzed in detail. In Section IV, the cosmological implications of the theory in its scalar-tensor representation are considered. In Section V, the behavior of the entropy of the Universe in the Lf (R, L m ) gravity theory, with the horizon entropy included, is analyzed. We discuss and conclude our results in Section VI.

II. NONMINIMAL CURVATURE-MATTER COUPLING
A. General formalism The action of f (R) gravity can be generalized with the introduction of a linear nonminimal coupling between matter and curvature. The corresponding gravitational theory, a particular case of the general f (R, L m ) theory [12], and which we denote as Lf (R, L m ), has the action given by [10], (1) where the factors f i (R) (with i = 1, 2) are arbitrary functions of the Ricci scalar R. The coupling constant λ determines the strength of the interaction between f 2 (R) and the matter Lagrangian. L m is the matter Lagrangian density, which is a function of the metric g µν and of the matter fields Φ. Now, varying the action with respect to the metric g µν provides the following field equations: where F i (R) = f ′ i (R), i = 1, 2, and the prime represents a derivative with respect to the scalar curvature R. We have defined for notational simplicity. The matter energy-momentum tensor is defined as An important property of any gravitational theory is its stability with respect to local perturbations.
In the standard f (R) gravity, a fatal instability (the Dolgov-Kawasaki instability) appears once the condition f ′′ (R) < 0 [24]. The instability develops on time scales of the order of 10 26 s. The stability properties of the gravitational models described by the action given by Eq. (1) were studied in [25], and it turns out that the corresponding stability criterion is f ′′ 1 (R) + 2λf ′′ 2 (R) > 0. In order to obtain the stability condition one expands the parameters of the model as the sum of a background field with constant curvature, and a small perturbation, so that With the use of the linearized field equations we obtain the stability condition ǫφ ′′ (R) + 2λf ′′ 2 (R) > 0 [25], generalizing the stability condition f ′′ (R) = ǫφ ′′ (R) > 0, found in f (R) gravity [24]. From a physical point of view the reason for stability is once the stability conditions are satisfied the effective mass m ef f of the dynamical degree of freedom associated to the small perturbation R 1 of the background curvature is non-negative.
A general property of these nonminimal curvaturematter coupling theories is the non-conservation of the energy-momentum tensor. This can be easily verified by taking into account the covariant derivative of the field Eq. (2), the Bianchi identities, ∇ µ G µν = 0, and the following identity, which then implies the following relationship: Note that in the absence of the coupling, λ = 0, one obtains the conservation of the energy-momentum tensor [27], which can also be verified through the diffeomorphism invariance of the matter part of the action. The conservation of the energy-momentum tensor also follows from Eq. (7), if f 2 (R) is a constant or the matter Lagrangian is not an explicit function of the metric. In order to test the motion in our model, we consider for the energy-momentum tensor of matter a perfect fluid where ρ is the total energy density and p, the pressure, respectively. The four-velocity, U µ , satisfies the conditions U µ U µ = 1 and ∇ ν U µ U µ = 0. We also introduce the projection operator h µλ = g µλ − U µ U λ from which one obtains h µλ U µ = 0. From Eq. (7), we deduce the equation of motion for a fluid element: where the extra force is given by: An intriguing feature is that the extra force depends on the form of the Lagrangian density. Note that considering the Lagrangian density L m = −p, where p is the pressure, the contribution of the nonminimal curvaturematter vanishes [28]. It has been argued that this is not the unique choice for the matter Lagrangian density and that more natural forms for L m , such as L m = ρ, do not imply the vanishing of the extra-force. Indeed, in the presence of the nonminimal coupling, they give rise to two distinct theories with different predictions [29,30].
B. Scalar-tensor representation of the linear curvature-matter coupling As has been shown in [31][32][33][34], f (R) gravity is equivalent to a scalar-tensor theory. In this context, the equivalence between the modified gravity models to a linear curvature-matter coupling and scalar-tensor gravity models was also established in [18]. More specifically, it was shown that the action given by Eq. (1) is equivalent with a two-potential scalar-tensor Brans-Dicke type theory, with a single scalar field, a vanishing Brans-Dicke parameter ω, and an unusual coupling of the second potential U (ψ) of the theory to matter.
As a first step in the scalar-tensor formulation of the theory we introduce a new field φ, and reformulate the action (1) as (11) Next, we introduce the second field ψ(φ) ≡ f ′ 1 (φ) (with a prime denoting a differentiation with respect to φ), and thus we obtain for the action the final expression where the two potentials V (ψ) and U (ψ) of the theory are defined as and respectively. The function φ(ψ) must be obtained by inverting ψ(φ) ≡ f ′ 1 (φ). The actions (1) and (12) are equivalent when f ′′ 1 (R) = 0 [18], similarly to the case of pure f (R) gravity [31][32][33][34].
The action given by Eq. (12) can be written, via a conformal rescaling g µν = exp (−αφ/2)ĝ µν , α = constant, as a four-dimensional dilaton gravity whose action, in the "Einstein frame", has the form [25] In the Einstein frame representation of the modified gravity with linear coupling between matter and geometry the extra force is due to the coupling between the matter Lagrangian and the Brans-Dicke-like scalar φ. On the other hand it is important to mention that in the Lf (R, L m ) theory there is no scaling of units with some powers of the conformal factor of the conformal transformation [25]. For this reason it is impossible to reduce Lf (R, L m ) gravity to a standard scalar-tensor theory, or to find a string gravity equivalent.
By varying the action (12) with respect to g µν provides the gravitational field equations of the scalar-tensor theory as while the variation of the action with respect to the field ψ gives the relation The contraction of the field equation Eq. (16) yields the scalar relation where T = T µ µ is the trace of the energy-momentum tensor. By combining Eqs. (17) and (18) we obtain the field equation of the field ψ as By eliminating the term ψ, and taking into account Eq. (19) the field equations (16) provides Now, taking the covariant divergence of the field equation (16), with the use of Eq. (6), we obtain first Then, by eliminating R/2 with the help of Eq. (17) we obtain for the divergence of the energy-momentum tensor Equation (22) allows the formulation of the energy and momentum balance equations in the scalar-tensor representation of the modified theory of gravity with a linear coupling between matter and geometry. By assuming that the energy-momentum tensor has the perfect fluid form given by Eq. (8), then Eq. (22) can be written in the equivalent form By multiplying Eq. (23) with U ν we obtain the energy balance equation in the scalar-tensor representation of the linear curvature-matter coupling given bẏ where we have introduced the Hubble function H = (1/3)∇ µ U µ , and we have denoted˙= U µ ∇ µ = d/ds, respectively, where ds is the line element corresponding to the metric g µν , ds 2 = g µν dx µ dx ν . After acting on Eq. (23) with the projection operator h νλ , provides the momentum balance equation for a perfect fluid as where the extra-force is given by (26) It is important to note that in the scalar-tensor representation of modified gravity with a linear curvaturematter coupling, the extra-force acting on test fluids is non-zero independently of the choice for the matter Lagrangian.

III. GRAVITATIONALLY INDUCED PARTICLE CREATION
In the present Section, we analyze the physical interpretation of the curvature-matter coupling in the scalartensor representation by adopting the point of view of the thermodynamics of open systems, in which matter creation irreversible processes may take place at a cosmological scale [19]- [22]. As we have already seen in the previous Section, the energy conservation equation of the curvature-matter coupling, given by Eq. (24), contains, as compared to the standard adiabatic conservation equation, an extra term, which can be interpreted in the framework of the open thermodynamic systems as an irreversible matter creation rate. According to irreversible thermodynamics, matter creation also represents an entropy source, generating an entropy flux, and thus leading, in the presence of the curvature-matter coupling, to a modification in the temperature evolution.
In the following, we investigate only the case in which all the non-diagonal components of the energymomentum tensor of the matter are equal to zero, so that T µν = 0, µ = ν. Generally, the energy-momentum tensor for a viscous dissipative fluid in the presence of heat conduction is given by T µν = (ρ + p + Π) − (p + Π)g µν + q µ u ν + q ν u µ + π µν , where Π is the bulk viscous pressure, q ν is the heat flux, and π µν is the tensor of viscous dissipation [26]. q µ and π µν must satisfy the conditions q µ u µ = 0 and π µν u ν = π µ µ = 0, respectively. In the following we neglect the viscous effects in the cosmological fluid, thus assuming Π ≡ 0 and π µν ≡ 0. In a comoving reference frame with u µ = (1, 0), the form of the heat flux vector is fixed by the normalization condition as q µ = (0, q). Therefore in a comoving reference frame all components of the form u µ q ν of the energy-momentum tensor are identically equal to zero, and T µν is a diagonal tensor. From the point of view of the thermodynamics of the irreversible processes, this condition implies the impossibility of heat transfer in the considered gravitational system. In particular, this condition is always satisfied in homogeneous and isotropic cosmological models, described by the Friedmann-Robertson-Walker geometry, since in these models the condition T 0i ≡ 0, i = 1, 2, 3 must always hold.

A. Matter creation rates and the creation pressure
We assume that the cosmological metric is given by the flat isotropic and homogeneous Friedmann-Robertson-Walker (FRW) metric, where a(t) is the scale factor, describing the expansion of the Universe. In this geometry the cosmological matter is comoving with the cosmological expansion, and therefore the four velocity of the cosmological fluid is U µ = (1, 0, 0, 0), while the Hubble function takes the form H =ȧ/a, since U µ ∇ µ =˙= d/dt. To investigate the thermodynamical implications at the cosmological scale with a curvature-matter coupling we consider that the Universe contains N particles in a volume V , with an energy density ρ and a thermodynamic pressure p, respectively. For such a cosmological system, the second law of thermodynamics, in its most general form, is given by [20] where dQ is the heat received by the system during time dt, and n = N/V is the particle number density, respectively. Due to our choice of the geometry of the Universe, and of the cosmological principle, only adiabatic transformations, defined by the condition dQ = 0, are possible. Therefore in the following we ignore proper heat transfer processes in the Universe. However, as one can see from Eq. (28), under the assumption of adiabatic transformations the second law of thermodynamics contains the term [(ρ+p)/n]d na 3 /dt, which explicitly takes into account the time variation of the cosmological particles in a given volume V . Hence, in the irreversible thermodynamics description of open systems, even for adiabatic transformations dQ = 0, one can consider the "heat" (internal energy), received/lost by the system, and which is due to the change in the particle number n. For modified gravity with a curvature-matter coupling the change in the particle number is due to the transfer of energy from gravity to matter. Thus, via matter creation, gravity acts as a source of internal energy, and of entropy. For adiabatic transformations dQ/dt = 0, we can reformulate Eq. (28) in an equivalent form aṡ Therefore, from the point of view of the thermodynamics of open systems, Eq. (24), giving the energy balance in the presence of a curvature-matter coupling, can be interpreted as describing particle creation in an homogeneous and isotropic geometry, with the time variation of the particle number density obtained aṡ where the particle creation rate Γ is a non-negative quantity defined as (31) Therefore, the energy conservation equation can be reformulated in the alternative forṁ As proven initially in [20], for adiabatic transformations Eq. (28), describing irreversible particle creation in an open thermodynamic system, can be rewritten as an effective energy conservation equation, or, in an equivalent form, as, where we have introduced a new thermodynamic quantity, p c , denoted the creation pressure and defined as [20] Therefore in modified gravity with a linear curvaturematter coupling the creation pressure is given by Note that from Eq. (7), the coupling between the matter and the higher derivative curvature terms may be interpreted as an exchange of energy and momentum between both. In the standard formulation of the linear curvature-matter coupling, by taking into account a FRW background, and from Eq. (7) we obtain the energy balance equation aṡ Hence, by considering that the mechanism for gravitational particle production is through the nonminimal curvature-matter coupling, so that comparing Eqs. (35) and (37), we have for the creation pressure with the requirement that p c be negative. We will only consider the case of L m = −p, i.e., α = −ω, in order to have a non-vanishing creation pressure.

B. Entropy and temperature evolution
According to the basic principles of the thermodynamics of open systems the entropy change consists of two components: the entropy flow term d e S, and the entropy creation term d i S. The total entropy S of an open thermodynamic system can be represented as [19,20] where by definition d i S > 0. Both the entropy flow and the entropy production can be obtained from the total differential of the entropy given by [20], where T is the temperature of the open thermodynamic system,s = S/a 3 is the entropy per unit volume, and µ is the chemical potential, defined as where h = ρ + p is the enthalpy of the system. In the case of a closed thermodynamic system and for adiabatic transformations we have dS = 0 and d i S = 0. However, in the presence of a curvature-matter coupling, leading to effective matter creation, there is a non-zero contribution to the total entropy. For a homogeneous and isotropic Universe the entropy flow term d e S vanishes, so that d e S = 0. On the other hand matter creation also represents a source for entropy creation, and the time variation of the corresponding entropy is obtained as [20] T Equation (42) gives the time variation of the entropy as so that the entropy increase due to particle production yields the expression where S 0 = S(0) is a constant. With the use of Eq. (43), we obtain for the entropy creation in the scalar-tensor representation of the linear coupling between matter and geometry the following equation The entropy flux four-vector S µ is defined as [21] S µ = nσU µ , where σ = S/N is the specific entropy per particle. S µ must satisfy the second law of thermodynamics, which imposes the constraint ∇ µ S µ ≥ 0 on its four-dimensional divergence. The Gibbs relation [21], together with the definition of the chemical potential µ of the open thermodynamic system, yields where we have used the relation which immediately follows from Eq. (29). With the use of Eq. (30) we obtain the entropy production rate due to the particle creation processes given by A general thermodynamic system is described by two fundamental thermodynamic variables, the particle number density n, and the temperature T , respectively. If the system is in an equilibrium state, the energy density ρ and the thermodynamic pressure p are obtained, in terms of n and T , from the equilibrium equations of state of the matter, ρ = ρ(n, T ), p = p(n, T ).
Therefore the energy conservation equation (32) can be obtained in the following general form By using the general thermodynamic relation [21] ∂ρ it follows that the temperature evolution of the newly created particles due to the curvature-matter coupling is given by the expressioṅ where the speed of sound c s is defined as c 2 s = ∂p/∂ρ. If the geometrically created matter satifies a barotropic equation of state of the form p = (γ − 1) ρ, 1 ≤ γ ≤ 2, the temperature evolution follows the simple equation C. Bulk-viscosity description of matter creation processes with a curvature-matter coupling An alternative physical interpretation of particle creation processes in cosmology was suggested by Zeldovich [35], and later on by Murphy [36] and Hu [37]. According to this interpretation, the viscosity of the cosmological fluid represents a phenomenological description of the effect of the creation of particles by the non-stationary gravitational field of the expanding universe. Therefore, from a physical point of view, a non-vanishing particle production rate is equivalent to the introduction of a bulk viscous pressure in the energy-momentum tensor of the cosmological fluid. From a quantum mechanical point of view, such a viscous pressure can also be related to the viscosity of the vacuum [35][36][37]. This physical interpretation follows from the simple circumstance that any source term in the energy balance equation of a general relativistic fluid may be formally rewritten in terms of an effective bulk viscosity [26].
The energy-momentum tensor of a general relativistic fluid with bulk viscosity as the only dissipative process can be written as [26] T µν = (ρ + p + Π) U µ U ν − (p + Π) g µν , where Π is the bulk viscous pressure. The particle flow vector N µ is defined as N µ = nU µ . In the framework of causal thermodynamics the entropy flow vector S µ takes the form [38] where τ is the relaxation time, and ξ is the coefficient of bulk viscosity. In Eq. (58), we have limited ourselves to considering only second-order deviations from equilibrium. In the case of homogeneous and isotropic geometries, in the presence of bulk viscous dissipative phenomena, the energy conservation equation is obtained aṡ By comparing Eq. (59), giving the energy conservation equation for a cosmological fluid in the presence of bulk viscosity, with Eq. (34), which includes in the energy balance the creation of particles from the gravitational field due to the curvature-matter coupling, it follows that these two equations are equivalent if (60) Therefore particle creation can be indeed described from a phenomenological point of view by introducing an effective bulk viscous pressure in the energy-momentum tensor of the cosmological fluid. Hence it follows that the causal bulk viscous pressure Π acts as a creation pressure.
Hence, it would be interesting to investigate matter creation processes in modified gravity with a curvaturematter coupling from the point of view of bulk viscous thermodynamic processes. Hence we shall consider in the following that there is a change in the number of particles, due to matter creation processes, with bulk viscous pressure playing the role of the creation pressure. We introduce a simple toy model in which the newly created particles obey, as a function of the particle number density n, an equation of state of the form where ρ 0 , n 0 and γ are constants, we have denoted k = ρ 0 /n γ 0 , and 1 ≤ γ ≤ 2, respectively. Using Eq. (61), then Eq. (59) takes the form of a particle balance equation, where is the particle production rate, proportional to the bulk viscous pressure. Combining the equation of state Eq. (61) with the Gibbs relation T ds = d(ρ/n) + pd(1/n) we obtains =s 0 = constant, that is, particles are created with constant entropy density.
However, there is a major difference between the particle creation irreversible processes in open thermodynamic systems and bulk viscous processes, and this difference is related to the expression for entropy production rate. While the entropy production rate associated to particle creation is given by [21] in the presence of bulk viscous dissipative processes the entropy production rate can be obtained as [26] ∇ (65) In the particle creation model in open thermodynamics systems the entropy production rate is proportional to the creation pressure, while in the viscous dissipative processes thermodynamic interpretation ∇ µ S µ is quadratic in the creation pressure, ∇ µ S µ ∝ p 2 c /ξT , and, moreover, involves a new dynamical variable, the bulk viscosity coefficient.

IV. COSMOLOGICAL APPLICATIONS
In the present Section, we consider several cosmological applications of the scalar-tensor formulation of modified gravity with a linear curvature-matter coupling, as its interpretation as a particle creation theory. For a homogeneous and isotropic geometry the gravitational field equations (17)-(20) take the form and respectively. Equation (30) can be integrated to give the particle number time variation as where n 0 is an arbitrary constant of integration, while, by assuming a barotropic equation of state of the form p = p (ρ) we obtain for the density evolution where ρ 0 is an integration constant.
A. Cosmological models satisfying the condition The linear curvature-matter coupling depends on the two arbitrary (and independent) potentials U (ψ) and V (ψ). As a simple toy model, in the following, we assume that the two potentials U and V are related to the matter Lagrangian via the relation Moreover, we restrict our analysis to the case of dust, with negligible thermodynamic pressure p = 0. With the choice of Eq. (71), the particle creation rate, given by Eq. (31), takes the form Therefore, in this approach the matter creation rate is determined by the function U (ψ), which describes the coupling between the matter Lagrangian and the scalar field, its derivative, and the time variation of the scalar field ψ only. With this choice, the variations of the particle number and of the matter energy density is given by With the assumption of Eq. (71) on the potentials, the gravitational field equations (66)-(68) take the form respectively. In the following, we consider only de Sitter type accelerating solutions of the system given by Eqs. (74)-(76), with H = H 0 = constant and a(t) = exp (H 0 t). Then we obtain first while the evolution equation for the scalar field ψ, which can be obtained from Eqs. (73)-(75), is given bÿ with the general solution given by where we have used the initial conditions ψ (t 0 ) = ψ 0 , and ψ ′ (t 0 ) = ψ 01 , respectively. Since for the dust fluid the matter Lagrangian is L m = ρ, Eq. (77) yields or, equivalently, The above equation determines the time variation of the potential U (ψ) as where U 0 is an arbitrary constant of integration. The time variation of the particle creation rate is obtained as where we have assumed ψ 01 < 0. The particle creation rate is non-negative for all times t, and is a monotonically increasing function of time for all t ≥ t 0 , and its initial value is given by Γ (t 0 ) = 12H 2 0 |ψ 01 |/ρ 0 ≥ 0. Therefore we obtain for the time variation of the matter energy density the equation (84) For the time variation of the potential V (t) we obtain the evolution equation which gives the potential V in an integral form as where V 0 is an arbitrary constant of integration. The creation pressure, defined as p c = −ρΓ/(3H), is given by Finally, for the comoving entropy of the de Sitter type expanding Universe in the linear curvature-matter coupling theory we obtain The entropy is a monotonically increasing function of time, with the propertyṠ(t) ≥ 0, ∀t ≥ t 0 . In the first order approximation, and for small times we obtain for the entropy the following expression Therefore, the curvature-matter coupling allows the production of a large amount of entropy during a de Sitter type evolutionary phase of the Universe.
B. de Sitter type expansionary models with constant matter creation rate In the following, in the scalar-tensor representation of Lf (R, L m ) gravity, we consider a second simple cosmological toy model by assuming that the cosmological expansion of a dust Universe, with p = 0, is accelerating with a = exp (H 0 t), where H 0 = constant, and the particle creation rate is a constant during the entire accelerating phase, and it is given by Γ = Γ 0 = 3H 0 = constant. Moreover, we take the matter Lagrangian as L m = ρ. Then from Eq. (70) it follows immediately that the matter density of the Universe is also a constant, Then the constancy of the matter creation rate, given by Eq. (31) imposes the following condition on the potentials U and V , while the field equation Eq. (68) yieldṡ From the field equations Eqs. (66) and (67) we obtain the evolution equation for ψ as From Eqs. (91) and (92) we obtain By taking the time derivative of the above equation, and by eliminatingψ with the help of Eq. (93), it follows that U satisfies the equation with the general solution given by where we have used the initial conditions U (t 0 ) = U 0 and U ′ (t 0 ) = U 01 , respectively. The time dependence of the potential V is obtained from Eq. (91) in the form With the use of the explicit time dependence of U in Eq. (93), we obtain for ψ(t) the following expression where we have used the initial conditions ψ (t 0 ) = ψ 0 anḋ ψ (t 0 ) = ψ 01 , respectively. Due to the cosmological particle production the entropy of the Universe increases as However, the specific entropy s = S/V remains a constant during the cosmological evolution, s = s 0 = constant.

V. TOTAL ENTROPY BEHAVIOR IN Lf (R, Lm) GRAVITY WITH PARTICLE CREATION
In the present paper, we have defined the entropy through the particle production rate, given by Eq. (43), as depending on the positive particle creation rate Γ via the relationṠ Therefore, in an ever expanding Universe with particle creation, the matter entropy will increase indefinitely. On the other hand, all natural systems tend to approach a state of thermodynamic equilibrium, implying that the entropy of equilibrium systems never decreases,Ṡ ≥ 0, and that it is concave when approaching the equilibrium state,S ≤ 0. However, in the present, and several other, cosmological models, these fundamental requirements for the behavior of the entropy do not seem to be satisfied. The problem of the validity of the second law of thermodynamics in cosmology was investigated in detail in [9,39], where it was shown that the Universe approaches thermodynamic equilibrium in a de Sitter phase, if one defines the total entropy S tot of the Universe as the entropy of the apparent horizon plus that of matter and radiation inside it. Then it follows that S tot increases, and that it is concave, thus leading to the result that the second law of thermodynamics is still valid for the case of the cosmological expansion. In the following, we investigate the thermodynamic properties of the total entropy in a Universe with matter creation.
In the standard thermodynamic description of physical systems the time parameter t is not a thermodynamic equilibrium variable. Therefore, the variation of the thermodynamic quantities should be considered with respect to some extensive variable. In the following we will adopt, following [9,39], as extensive variable for the cosmological system the proper volume enclosed by the apparent horizon, or, more specifically, its scale factor a. Then the relation between d/dt and d/da is simply In the following we denote by a prime the derivative with respect to the extensive variable a.
We define the total entropy of a FRW Universe with dust as the sum of the entropy of the apparent horizon S ah , proportional to its area, and that of the matter particles within it S m [9,39]. For practical purposes, in the case of the flat FRW model, the total entropy is where we have used the fact that the radius of the apparent horizon is r ah = H −1 [40]. We also introduce an important observational quantity, the deceleration parameter q, defined as Therefore, the variation of the total entropy can be obtained as respectively. The particle number n satisfies Eq. (30), and is rewritten as aH dn da + 3nH = Γn.
From its definition the derivatives of the entropy with respect to the scale factor can be evaluated as respectively.
In terms of the deceleration parameter we can express the variation of the total entropy with respect to a as Therefore the standard thermodynamic requirements S ′ tot ≥ 0 and S ′′ tot ≤ 0 impose the following constraints on the particle creation rate Γ, and its derivative with respect to the scale factor respectively. Since in the Lf (R, L m ) gravity theory the matter creation rate Γ is determined by the coupling functions U and V , the thermodynamic conditions impose some strong constraints on the allowed physical form of these functions. A particularly interesting case is that of the de Sitter evolution of the Universe, with H = H 0 = constant. For this situation the total entropy of the Universe is given by giving The thermodynamic condition of the non-negativity of the total entropy derivative with respect to the scale factor imposes the conditions Γ ≥ 3H 0 and Γ ′ (a) ≤ 7Γ(a) − Γ 2 (a)/H 0 − 12H 0 /a on the particle creation rate Γ. In the particular case Γ = 3H 0 , we obtain S tot = constant, showing that in this case the cosmological evolution is isentropic, with the total entropy being a constant.

VI. DISCUSSION AND CONCLUSIONS
In the present paper, we have considered the thermodynamic interpretation of modified theories of gravity with a linear coupling between matter and geometry, which we denote as Lf (R, L m ) gravity. This theory represents a particular class, corresponding to a specific choice of the gravitational Lagrangian, of a very general class of theories, in which the action is an arbitrary function of the Ricci scalar and of the matter Lagrangian. An interesting characteristic of these theories is the nonconservation of the energy-momentum tensor of the matter, indicating that matter and energy fluxes can be generated by the conversion of the geometric curvature, describing the gravitational field, into matter. Hence the presence of matter, and its possible coupling to geometry could modify the cosmological evolution in a way that goes far beyond the standard description of general relativity. The presence of a source term in the energy balance equation can be naturally interpreted in the framework of the thermodynamics of open systems as describing a particle creation process, in which the "geometric energy" of the gravitational field is transferred to "real" matter. During the particle production phase a large amount of entropy is produced.
In order to estimate the effective thermodynamics quantities we have first introduced the equivalent scalartensor representation of the Lf (R, L m ) theory, which can be formulated in terms of a scalar field ψ, with two independent potentials V (ψ) and U (ψ), with the potential U (ψ) coupled to the matter Lagrangian. Using the scalar-tensor representation of the Lf (R, L m ) theory, we have obtained the particle creation rate, the creation pressure and the entropy associated to the gravitational energy transfer to matter. The cosmological implications of the particle creation have also been investigated, by assuming a specific relation between the two potentials. The imposed condition makes the particle creation rate a function of the second scalar potential U (ψ), which directly couples to the matter Lagrangian. The gravitational field equations corresponding to these choices have a de Sitter type accelerating solution, where the cosmic acceleration is triggered by the particle creation process, which generates a negative creation pressure. Thus, it was argued that the negative creation pressure is responsible for the accelerated expansion of the Universe.
Matter creation processes are supposed to play a fundamental role in the quantum field theoretical approaches to gravity, where they naturally appear. It is a standard result of quantum field theory in curved spacetimes that quanta of the minimally-coupled scalar field are created in the expanding Friedmann-Robertson-Walker universe [41]. That's why finding an equivalent microscopic quantum description of the matter creation processes considered in the present paper could shed some light on the physical mechanisms leading to particle generation via gravity and matter geometry coupling. In the following we will briefly point out that such mechanisms do exist, and can be understood, at least qualitatively, in the framework of some semiclassical gravity models.
In semiclassical gravity it is assumed that the gravitational field remains classical, while the classical bosonic fields φ are quantized. In order to couple quantized fields to classical gravitational fields the quantum energy momentum tensorT µν is replaced by its expectation value with respect to some quantum state Ψ, thus leading to the effective semiclassical Einstein equation [42], Hence the classical energy-momentum tensor of the system T µν is defined as Ψ|T µν |Ψ = T µν . The semiclassical equation Eq. (118) can be obtained from the variational principle [43] δ (S g + S ψ ) = 0, where x is the classical action of the gravitational field, and whereĤ is the Hamiltonian operator of the system, and α is a Lagrange multiplier. The variation of Eq. (119) provides the normalization condition for the wave function Ψ|Ψ = 1, the Schödinger equation for the wave function as well as the semiclassical Einstein Eq. (118). In this simple case the Bianchi identities require the conservation of the energy-momentum tensor, ∇ µ Ψ|T µν |Ψ = 0.
A very different set of semiclassical Einstein equations can be obtained by assuming a coupling between the quantum fields and the curvature of the space-time. In the model introduced in [43] the contribution to the total action of the geometry-quantum matter coupling term was assumed to be of the form where F and f are arbitrary functions, and ( f (φ) ) Ψ = Ψ(t)| f [φ(x)] |Ψ(t) . Then, in the presence of such a geometry-matter coupling the Hamiltonian H(t) in the Schrödinger Eq. (121) is modified to [43] where N is the lapse function, ξ i are intrinsic coordinates, such that the normal is everywhere time-like, and γ = det γ rs , where γ rs is the metric induced on a surface σ(t), which gives a global slicing of the space-time into space-like surfaces. The effective semiclassical Einstein equation takes the form [43] In Eq. (124) the matter energy-momentum tensor is not conserved, ∇ µ T µν Ψ = 0. Thus, this equation describes an effective particle production process, and can be interpreted as giving an effective semiclassical description of the quantum processes in a gravitational field. By modifying the classical part of the gravitational action we can recover the field equations Eqs. (2) and (16), respectively, used in the present paper. Therefore the physical origin of the matter creation processes considered in the present paper can be traced back to the semiclassical approximation of the quantum field theory in a Riemannian curved geometry.
An interesting and important question is the physical nature of the particles that could be created via gravitationally induced creation processes. The most natural assumption would be that these particles are dark matter particles. It has been conjectured that dark matter may consist of ultra-light particles with masses of the order of m ≈ 10 −24 eV (see [44] and references therein). From a physical point of view such a particle may represent a pseudo Nambu-Goldstone boson. Axions are other ultralight dark matter candidates, with masses in the range m ≤ 10 −22 eV [45]. Such extremely very low mass particles can be created even in very weak gravitational fields. An alternative description of dark matter is provided by the so-called scalar field dark matter models [46],in which it is assumed that dark matter is a real scalar field, minimally coupled to gravity, with the mass of the scalar field having a very small value of the order of m < 10 21 eV. For zero temperature scalar field dark matter models all particles in the system condense to the same quantum ground state, thus forming a Bose-Einstein condensate. Therefore scalar field dark matter models are equivalent to the Bose-Einstein condensate dark matter models [47]. This implies, from a physical point of view, that in the open irreversible thermodynamic model introduced in the present paper particle creation can take place also in the form of a scalar field. In such a model the evolution of the scalar field dark energy particles, with energy density ρ φ and pressure p φ , and having a particle number density n φ , is governed by an equation of the forṁ where Γ 1 is the particle decay rate, determined by the coupling between matter and geometry. For the energy density and pressure of the scalar field dark matter we can assume the standard form where U (φ) is the scalar field self-interaction potential. The creation pressure corresponding to the scalar field creation processes can be obtained as It is interesting to note that Eq. (125), which describes the creation of a scalar field as a result of the geometrymatter coupling, can be written in an equivalent form asφ where we have denoted Γ φ,φ, U = Γ 1 ρ φ /n φ . Therefore in the scalar field dark matter model a friction term in the scalar field evolution equation Eq. (128) does appear naturally, and in a general form, as a direct consequence of the irreversible thermodynamics of open systems as applied to the dark matter case. Hence scalar field dark matter can be a result of the cosmological particle production due to the geometry-matter coupling in modified gravity theories. For gravitational models with an action given by an arbitrary function of the Ricci scalar, the matter Lagrangian density, a scalar field and a kinetic term constructed from the gradients of the scalar field, respectively, see [48]. The Lf (R, L m ) gravitational theory investigated in the present paper predicts the possibility that matter creation, associated with the curvature-matter coupling, could also occur in the present-day universe, as proposed by Dirac [49] a long time ago. The late expansion of the Universe [1, 2] may be considered as an empirical evidence for matter creation, and a viable alternative to the mysterious dark energy. Presently the existence of some forms of the curvature-matter coupling leading to matter creation processes cannot be fully ruled out by the existing cosmological observations or by astrophysical data. Presumably, the functional forms of the potentials V (ψ) and U (ψ) that completely characterize the Lf (R, L m ) gravitational theory will be provided by fundamental quantum field theoretical models of the gravitational interaction, thus opening the possibility of an in depth comparison of the predictions of the Lf (R, L m ) gravity with cosmological and astrophysical observational data.