Palatini formulation of f(R, T) gravity theory, and its cosmological implications

We consider the Palatini formulation of f(R, T) gravity theory, in which a non-minimal coupling between the Ricci scalar and the trace of the energy-momentum tensor is introduced, by considering the metric and the affine connection as independent field variables. The field equations and the equations of motion for massive test particles are derived, and we show that the independent connection can be expressed as the Levi-Civita connection of an auxiliary, energy-momentum trace dependent metric, related to the physical metric by a conformal transformation. Similar to the metric case, the field equations impose the non-conservation of the energy-momentum tensor. We obtain the explicit form of the equations of motion for massive test particles in the case of a perfect fluid, and the expression of the extra force, which is identical to the one obtained in the metric case. The thermodynamic interpretation of the theory is also briefly discussed. We investigate in detail the cosmological implications of the theory, and we obtain the generalized Friedmann equations of the f(R, T) gravity in the Palatini formulation. Cosmological models with Lagrangians of the type f=R-α2/R+g(T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=R-\alpha ^2/R+g(T)$$\end{document} and f=R+α2R2+g(T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=R+\alpha ^2R^2+g(T)$$\end{document} are investigated. These models lead to evolution equations whose solutions describe accelerating Universes at late times.


Introduction
The observational discovery of the recent acceleration of the Universe [1][2][3][4][5] has raised the fundamental theoretical problem if general relativity, in its standard formulation, can fully account for all the observed phenomena at both galactic and extra-galactic scales. The simplest theoretical explanation for the observed cosmological dynamics consists in slightly modifying the Einstein field equations, by adding to it a cosmological constant [6]. Together with the assumption of the existence of another mysterious component of the Universe, called dark matter [7,8], assumed to be cold and pressureless, the Einstein gravitational field equations with the cosmological constant included, can give an excellent fit to all observed data, thus leading to the formulation of the standard cosmological paradigm of our present days, called the Cold Dark Matter ( CDM) model. However, despite its apparent simplicity and naturalness, the introduction of the cosmological constant raises a number of important theoretical and observational question for which no convincing answers have been provided so far. The CDM model can fit the observational data at a high level of precision, it is a very simple theoretical approach, it is easy to use in practice, but up to now no fundamental theory can explain it. Why is the cosmological constant so small? Why is it so fine-tuned? And why did the Universe begin to accelerate only recently? And, after all, would a cosmological constant really be necessary to explain all observations?
From a theoretical point of view two possible answers to the questions raised by the observation of the recent acceleration of the Universe can be formulated. The first, called the dark energy approach, assumes that the Universe is filled by a mysterious and unknown component, called dark energy [9][10][11][12], which is fully responsible for the acceleration of the Universe, as well as for its mass-energy balance. The cosmological constant corresponds to a particular phase of the dynamical dark energy (the ground state of a potential, say), and the recent de Sitter phase may prove to be just an attractor of the dynamical system describing the cosmological evolution. A second approach, the dark gravity approach, assumes the alternative possibility that at large scales the gravitational force may have a very different behavior as compared to the one suggested by standard general relativity. In the general relativistic description of gravity, the starting point is the Hilbert-Einstein action, which can be written down as S = R/2κ 2 + L m √ −g d 4 x, where R is the Ricci scalar, κ is the gravitational coupling constant, and L m is the matter Lagrangian, respectively. Hence in dark gravity type theories for a full understanding of the gravitational interaction a generalization of the Hilbert-Einstein action is necessary.
There are (at least) two possibilities to construct dark gravity theories. The first is based on the modification of the geometric part of the Hilbert-Einstein Lagrangian only. An example of such an approach is the f (R) gravity theory, introduced in [13,14], and in which the geometric part of the action is generalized so that it becomes an arbitrary function f (R) of the Ricci scalar. Hence in f (R) gravity the total Hilbert-Einstein action can be written as S = f (R)/2κ 2 + L m √ −g d 4 x. The recent cosmological observations can be satisfactorily explained in the f (R) theory, and a solution of the dark matter problem, interpreted as a geometric effect in the framework of the theory, can also be obtained [15]. For reviews and in depth discussions of f (R) and other modified gravity theories see [16][17][18][19][20][21][22][23][24]. A second avenue for the construction of the dark gravity theories consists in looking for maximal extensions of the Hilbert-Einstein action, in which the matter Lagrangian L m plays an equally important role as the Ricci scalar. Hence in this more general approach one modifies both the geometric and the matter terms in the Hilbert-Einstein action, thus allowing a coupling between matter and geometry [25,39]. The first possibility for such a coupling is to replace the gravitational action by an arbitrary function of the Ricci scalar and the matter Lagrangian L m , thus obtaining the so-called f (R, L m ) class of modified gravity theories [26]. This class of theories has the potential of explaining the recent acceleration of the Universe without the need of the cosmological constant, and can give some new insights into the dark matter problem, and on the nature of the gravitational motion. The cosmological and physical implications of this theory have been intensively investigated [27][28][29][30][31][32][33][34][35][36][37]. For a review of f (R, L m ) type theories see [38].
A second extension of the Hilbert-Einstein action can be obtained by assuming that the gravitational field couples to the trace T of the energy-momentum tensor of the matter. This assumption leads to the f (R, T ) class of gravitational theories [39]. f (R, T ) theory may give some hints for the existence of an effective classical description of the quantum properties of gravity. As pointed out in [40], by using a non-perturbative approach for the quantization of the metric, proposed [41][42][43], as a consequence of the quantum fluctuations of the metric, a particular type of f (R, T ) gravity naturally emerges, with the Lagrangian given by constant. This interesting theoretical result suggests that a deep connection may exist between the quantum field theoretical description of gravity, which naturally involves parti-cle production in the gravitational field, and the corresponding effective classical description of the f (R, T ) gravity theory [44]. The astrophysical and cosmological implications of f (R, T ) gravity theory were investigated in .
Einstein's theory of general relativity can be obtained by starting from two different theoretical approaches, called the metric and the Palatini formalism, respectively, the latter being introduced by Albert Einstein [67][68][69]. In the Palatini variational approach one takes as independent field variables not only the ten components g μν of the metric tensor, but also the components of the affine connection α βμ , without assuming, a priori, the form of the dependence of the connection on the metric tensor, and its derivatives [70,71]. When applied to the Hilbert-Einstein action, these two approaches lead to the same gravitational field equations. Moreover, the Palatini formalism also provides the explicit form of the symmetric connection as determined by the derivatives of the metric tensor. However, in f (R) modified gravity, as well as in other modified theories of gravity, this does not happen anymore. In fact it turns out that the gravitational field equations obtained by using the metric approach are generally different from those obtained by using the Palatini variation [70,71]. An important difference is related to the order of the field equations. The metric formulation usually leads to higher-order derivative field equations, while in the Palatini formalism the obtained gravitational field equations are always second-order partial differential equations. A number of new algebraic relations also do appear in the Palatini variational formulation, which describe the subtle relation between the matter fields and the affine connection, which can be determined from a set of equations that couples it not only to the metric, but also to the matter fields. The astrophysical and cosmological implications of the Palatini formulation of f (R) gravity have also been intensively investigated [72][73][74][75][76].
Based on a hybrid combination of the metric and Palatini mathematical formalisms, an extension of the f (R) gravity theory was proposed in [77] and was used to construct a new type of gravitational Lagrangian [77,78]. A simple example of such an hybrid metric-Palatini theory can be constructed by adopting for the gravitational Lagrangian the expression is the Palatini scalar curvature. A similar formalism that interpolates between the metric and Palatini regimes was proposed in [79,80] for the study of f (R) type theories. This approach is called C-theory. A generalization of the hybrid metric-Palatini gravity was introduced in [81].
Despite the intensive investigations of the theoretical and observational aspects of the modified gravity theories with geometry-matter coupling, their Palatini formulation and properties have attracted considerably less attention. The Palatini formulation of the linear f (R, L m ) gravity was introduced in [82], where the field equations and the equa-tions of motion for massive test particles were derived. The independent connection can be expressed as the Levi-Civita connection of an auxiliary, matter Lagrangian dependent metric, which is related to the physical metric by means of a conformal transformation. Similar to the metric case, the field equations impose the non-conservation of the energy-momentum tensor. The study of Palatini formulation of f (R, T ) gravity was initiated in [83]. Analogously to its metric counterpart, the field equations impose of the f (R, T ) gravity in the Palatini formulation implies the nonconservation of the energy-momentum tensor, which leads to non-geodesic motion, and to the appearance of an extra force.
It is the purpose of the present paper to derive the gravitational field equations of the generalized f (R, T ) type gravity models, with non-minimal coupling between matter, described by the trace of the energy-momentum tensor, and geometry, characterized by the Ricci scalar By taking separately two independent variations of the gravitational action with respect to the metric and the connection, respectively, we obtain the field equations and the connection associated to the Ricci tensor, which, due to the coupling between the trace of the energy-momentum tensor and the geometry, is also a function of T . The metric that defines the new independent connection is conformally related to the initial spacetime metric, with the conformal factor given by a function of the trace of the energy-momentum tensor, and of the Ricci scalar. After the conformal factor is obtained, the gravitational field equations can be written down easily in both metrics. Similar to the case of the metric f (R, T ) gravity, after taking the divergence of the gravitational field equations we obtain the important result that the energy-momentum tensor of the matter is not conserved. Similar to the metric case, the motion of the particles is not geodesic, and due to the matter-geometry coupling, an extra force arises. However, this force has the same expression as in the metric case, and therefore no new physics is expected to arise during the motion of massive test particles in the Palatini formulation of f (R, T ) gravity. As the next step in our analysis we investigate in detail the cosmological implications of the Palatini formulation of the f (R, T ) gravity theory. We obtain the generalized Friedmann equations, which explicitly contain the extra terms generated by the coupling between the trace of the energy-momentum tensor and geometry. The general properties of the cosmological evolution are obtained, including the behavior of the deceleration parameter, of the effective energy density and pressure, and of the parameter of the equation of state of the dark energy. Cosmological models with Lagrangians of the type L = R − α 2 /R + g(T ) and L = R + α 2 R 2 + g(T ) are considered in detail, and it is shown that these models lead to evolution equations whose solutions tend to a de Sitter type Universe at late times.
The present paper is organized as follows. After a brief review of the metric formalism, the field equations of f (R, T ) gravity theory are obtained by using the Palatini formalism of gravitational theories in Sect. 1. The energy and momentum balance equations are obtained, after taking the divergence of the energy-momentum tensor, in Sect. 3. The thermodynamical interpretation of the theory is also briefly discussed. The cosmological implications of the Palatini f (R, T ) theory is investigated in Sect. 4. We discuss and conclude our results in Sect. 5. The details of the derivation of the field equations in the metric formalism are given in Appendix A, while the divergence of the matter energymomentum tensor is derived in Appendix B. The explicit computations of the various geometric quantities for the Friedmann-Robertson-Walker geometry are presented in Appendix C.

Palatini formulation of f (R, T ) gravity
In the present section, after a brief review of the metric formulation of f (R, T ) gravity theory, we derive the field equations of the theory by using the Palatini formalism.

The metric formalism
The f (R, T ) gravity theory is described by the action [39] where g ≡ det g μν , f is an arbitrary function of the Ricci scalar R = R(g) and of the trace T = g μν T μν of the matter energy-momentum tensor T μν ; the matter Lagrangian L m is assumed to be independent of ∂ λ g μν . T μν is generally obtained as [89] T To describe the variation of the energy-momentum tensor with respect to the metric, we also introduce the tensor μν , defined as μν ≡ g αβ δT αβ δg μν .
1 Throughout this article we use the natural system of units with c = G = 1. For the metric tensor we adopt the signature convention (−, +, +, +).
For a perfect fluid characterized by its energy density ρ and isotropic pressure P only, the energy-momentum tensor is given by where the four-velocity u μ satisfies the normalization condition u μ u μ = − 1. In the comoving frame its components are u μ = (− 1, 0, 0, 0), and in this frame the components of the energy-momentum tensor become T μ ν = (− ρ, P, P, P). The components of the affine connection are defined to be which are currently regarded as functions of the metric. In the following we assume that the connection is symmetric, that is, ρ μν = ρ νμ . Varying Eq. (1) with respect to g μν , we obtain where From the condition δS = 0 we obtain the field equations of f (R(g), T ) gravity theory as (for the computational details see Appendix A), 2.2 Palatini formulation of f (R, T ) gravity

Field equations from metric variation
An alternative formulation of gravitational theories can be obtained within the Palatini formalism, which consists in taking separately in the gravitational action two independent variations, with respect to the metric and the connection, respectively. The action is formally identical to the metric one, but the Riemann tensor and the Ricci tensor are constructed with the independent symmetric connection˜ . Hence in the Palatini formulation the gravitational action of f (R, T ) gravity is given by In Eq. (11) the Ricci scalar is defined as with the Ricci tensorR μν ˜ expressed only in terms of the Palatini connection˜ , with the connection coefficients λ μν determined self-consistently through the independent variation of the gravitational field action Eq. (11), and not constructed directly from the metric by using the usual Levi-Civita definition. In the following we defineR μν ˜ with the help of the yet undetermined Palatini connection as The matter Lagrangian L m (g, ψ) is assumed to be a function of the metric tensor g and of the physical fields ψ only.
We vary now the gravitational action (11) with respect to the metric tensor g μν , under the assumption δR μν ˜ = 0, that is, by keeping the connection constant. As a result we immediately obtain the field equations By contracting the above equation with g μν we obtain for the Ricci scalar (12) the expression where = μ μ . In terms of the Einstein tensor the Palatini field equations can be written as

The Palatini connection
We vary now the gravitational action (11) with respect to the connection˜ , by keeping the metric constant, so that According to the Palatini identity [82], we have where∇ λ is the covariant derivative associated with˜ . Hence the variation of the action (11) with respect to˜ leads to where we have denoted We integrate now by parts to obtain The first term in Eq. (22) is a total derivative, and thus after transforming it into a surface integral it vanishes. Therefore the variation of the action with respect to the connection˜ becomes Equation (23) can be significantly simplified by taking into account that for α = λ the equation is identically zero. Hence for the case λ = α, we find Equation (24) shows that the connection˜ is compatible with a conformal metricg μν , conformally related to the initial metric g μν by means of the relations with the conformal factor F defined as Moreover, we have whereg ≡ det g μν . Thus Eq. (24) gives the geometric interpretation of the Palatini connection˜ as the Levi-Civita connection corresponding to the metricg μν , conformally related to g μν ,

Conformal geometry and g frame field equations
Since the metrics g μν andg μν are conformally related, the connection˜ can be expressed in terms of the Levi-Civita connection , whose components are given in Eq. (5), as In terms of the tensor R μν (g), constructed from the metric by using the Levi-Civita connection (5), the Ricci tensorR μν in the conformally transformed metric is given by [84,85] SinceR ≡g μνR μν = F −1 R g,˜ andG μν ≡R μν − g μνR /2 = G μν g,˜ , the Ricci scalar (12) and the Einstein tensor in Eq. (17) can be obtained in the conformally related frames [84,85]: and respectively, with all the covariant derivatives and algebraic operations performed with the help of the metric g μν .
By using the expression of the Einstein tensor as given by Eq. (32) in Eq. (17), we obtain finally the f (R, T ) field equations of the Palatini formulation expressed solely in the g frame: The Ricci scalar can be obtained:

Field equations in theg frame
In the previous subsection we have obtained the gravitational field equations in the Palatini formulation of f (R, T ) gravity as expressed in terms of the metric tensor g μν . This approach usually involves higher-order differential equations for the physical and geometrical quantities. An alternative approach, which keeps the order of the differential equations of the model not higher than two would be to solve first the gravitational field equations in the conformal metricg μν , and to recover the metric g μν with the use of the conformal transformation (25). To obtain the field equations in the conformal frame we follow the procedure introduced in [83]. First we multiply the Palatini field equation (14) with whereR λ ν is now a function of the metricg only. As for the energy-momentum tensor, we have, respectively, and That is, the mixed components T μ ν of the energy-momentum tensor are conformally invariant. Additionally, we have By contracting Eq. (35) by taking ν = λ we obtain the expression of the Ricci scalar in the conformal frame, thus Therefore in the conformal frame the full set of the Einstein equations can be written with These equations determine the conformal metricg as a function of the thermodynamic parameters that enter in the definition of the matter energy-momentum tensor.

The Newtonian limit
To investigate the Palatini f (R, T ) gravity under the weak field, slow motion and static approximation, namely the Newtonian limit, we assume the metric to be a Minkowski metric plus a perturbation, given by Hence g μν = η μν − γ μν and the relation g αλ g λβ = δ α β still holds. In this context, we also assume the conformal metric (25) to be nearly flat, so that whereγ μν 1 is of the same order as γ μν . Hence F ≈ 1. From Eqs. (41) and (42) we obtain If we take F = e 2W and expand it to F = 1 + 2W , then the above equation Let us now consider the Palatini field equation (14) under the Newtonian limit. First we obtain the g-frame Ricci tensor: Omitting all higher-order terms with respect to O(γ μν ) and taking into account the gauge then some algebra gives the expression of the Ricci tensor (30) as With the use of the above equations, and preserving only the first-order terms, the field equation (14) becomes For perfect fluids, T μ ν = diag(−ρ, P, P, P) and (61)); besides, in the Newtonian limit, P → 0 and ∂ 0 → 0. Hence we can obtain immediately the generalized Poisson equation in the Palatini formulation of f (R, T ) gravity: The same result can also be found in [83].

Violation of the equivalence principle
An interesting feature of the modified gravities, including their Palatini extensions, is the violation of the equivalence principle. In Palatini f (R) gravity this problem was discussed in [86]. In the following we will generalize some results from the case of f (R) theory to the Palatini f (R, T ) gravity. In the conformal frame the field equations of the Palatini f (R, T ) gravity are given by Eq. (40), In the weak field limit we can represent the gravitational Lagrangian as where , are constants, K R(g,˜ ) is an arbitrary function of the argument R(g,˜ ), while g is an arbitrary function of the trace of the matter energy-momentum tensor.
In the limit of small , , with , → 0, it follows that f − R(g,˜ )F ≈ 0. In addition, by neglecting the matter energy-momentum tensor T μ ν for weak sources, it follows that Eq. (40) becomesG μ ν ≈ 0, which leads tog μν ≈ η μν , or, equivalently, The above equation tells us that, similar to the metric and the Palatini formulation of f (R) gravity [86], in the Palatini formulation of f (R, T ) gravity it is impossible to recover the flat Minkowski metric even in local frames with external gravitational fields screened. This result violates the basic postulate of general relativity according to which in locally free falling frames the non-gravitational laws of physics are those of special relativity [86]. Since this postulate assumes that the Einstein equivalence principle holds, it follows that, similar to Palatini f (R) theory, in the Palatini formulation of f (R, T ) gravity the equivalence principle does not hold exactly. The deviation of the current metric with respect to the Minkowski metric is In order to give a quantitative estimate for the deviation of the f (R, T ) metric from the Minkowski metric we will consider the cosmological case, to be discussed in detail in Sect. 4. We assume that the gravitational action takes the form and w is the parameter of the matter equation of state. By estimating all quantities at the present time t = t 0 , then when α, β → 0, and w = 0, a (t 0 ) = 1, the deviation (∂ K /∂ R) from the Minkowki metric can be obtained: where H 0 is the present-day value of the Hubble function, and we have assumed for the present-day density of the Universe the critical value.
Recently the first results obtained by the MICROSCOPE satellite, whose aims are to constrain the weak equivalence principle in outer space by determining the Eötvös parameter η, have been published [87]. The Eötvös parameter is defined as the normalized difference of accelerations between two bodies i and j, located in the same gravitational field. The MICROSCOPE determinations give for η the value η = (− 1 ± 27) × 10 −15 at a 2σ confidence level [87]. These results allow one to constrain possible sources of violation of the weak equivalence principle, like, for example, the existence of light or massive scalar fields with coupling to matter weaker than the gravitational coupling [88]. For a massive scalar field of mass smaller than 10 −12 eV, the coupling is constrained as |α C | < 10 −11 , if the scalar field couples to the baryon number, and to |α C | < 10 −12 if the scalar field couples to the difference between the baryon and the lepton numbers, respectively. We expect a similar order of magnitude for the coupling between matter and geometry in both the metric and the Palatini formulations of f (R, T ) gravity.

Energy and momentum balance equations
An interesting and important consequence of modified gravity theories with geometry-matter coupling is the nonconservation of the matter energy-momentum tensor. This property of the theory has a number of far reaching physical implications, and it may represent the main link between the interpretation of f (R, T ) theory as an effective classical description of the quantum theory of gravity. In the present section we obtain the general expression of the divergence of the energy-momentum tensor in f (R, T ) gravity theory, and, by using it, we obtain the energy-momentum balance equations, which describe the energy transfer processes from geometry to matter, and the deviations from the geodesic motion, respectively.

The divergence of the matter energy-momentum tensor
We begin our analysis by calculating first the divergence of the Einstein tensor in the Palatini frame. Since ∇ μ G μν (g) = 0, Eq. (32) yields where to simplify the calculation we have taken F = e 2w , and we have used the mathematical identities [84] and respectively, where φ (x μ ) is an arbitrary scalar. We have also used the relation g αβ R μβαν = −g αβ R βμαν = −R μν in the last step of Eq. (53). Note that the above two identities are valid in both the metric and the Palatini formulations [84]. Substituting Eq. (30) into Eq. (52), we find The covariant divergence of the field equation (17) yields (note that G μν g,˜ =G μν and g μν R g,˜ =g μνR ) where we have used Eq. (55). One can check now by comparison with Eq. (10) that the above expression gives the same result as in the metric formulation, except for the functional form of the function f , From its definition with the use of Eq. (2), the tensor μν for perfect fluids can be obtained: where the relation δg αβ δg μν = −g αμ g βν (59) can be derived from the relations g αβ δg αβ + g αμ g βν δg μν = g αβ δg αβ + g μν δg μν = 0, ⇒ δg αβ = −g αμ g βν δg μν .

Balance equations in the conformal frame
In the following for notational simplicity we define first a mixed-component vector field, Since∇ μG μ ν ≡ 0, from Eq. (40) we at once get the conservation equations in the conformal frame: where (15), and we have used Eq. (28). Taking into account that F is a scalar and∇ μ F = ∇ μ F = ∂ μ F, then or equivalently where we have used the relation 8π T − (T + ) f T + f = F R g,˜ − f on the left hand side. Hence one can easily see that the equations above are exactly the same as the energy balance equation (57).

Thermodynamic interpretation of f (R, T ) gravity theories
For the sake of completeness we briefly present the thermodynamic interpretation of f (R, T ) gravity theories, as discussed in [46]. The non-conservation of the matter energymomentum tensor strongly suggests that, due to the mattergeometry coupling, a particle creation processes may take place during the cosmological evolution. This phenomenon is also specific to quantum field theories in curved spacetimes, as pointed out in [90][91][92], and it is a consequence of a time varying gravitational field. Hence, f (R, T ) theory, which also involves particle creation, may lead to the possibility of a semiclassical description of quantum field theoretical processes in gravitational fields.

Particle and entropy fluxes, and the creation pressure
Particle creation implies that the covariant divergence of the basic equilibrium quantities, including the particle and entropy fluxes, as well as of the energy-momentum tensor, are now different from zero. Consequently, all the balance equilibrium equations must be modified to include particle creation [93][94][95]. In the presence of gravitationally generated matter, the balance equation for the particle flux N μ ≡ nu μ , where n is the particle number density, becomes where is the particle production rate, which can be neglected in the case that H . The entropy flux vector is defined to be S μ ≡ su μ = nσ u μ , where s is the entropy density, and σ is the entropy per particle. The divergence of the entropy flux gives If we consider a specific σ which is constant, then that is, the variation of the entropy is entirely due to (adiabatic gravitational) particle creation processes. Since s > 0, from the above equation it follows that the particle creation rate must satisfy the condition ≥ 0, that is, gravitational fields can generate particles, but the inverse process is prohibited. The energy-momentum tensor of a fluid in the presence of particle creation must also be modified to take into account the second law of thermodynamics, so that [96] T μν = T μν eq + T μν , where T μν eq denotes the equilibrium component (4), and T μν is the correction due to particle creation. Due to the isotropy and homogeneity of space-time, the extra contribution to the equilibrium energy-momentum tensor must be represented by a scalar process. Generally one can write where P c is the dynamic creation pressure that describes phenomenologically the thermodynamic effect of particle creation in a macroscopic system. In a covariant representation we have [96] T μν = −P c h μν = −P c g μν + u μ u ν , which immediately gives u μ ∇ ν T μν = 3H P c . Therefore in the presence of particle creation the total energy balance equation u μ ∇ ν T μν = 0, which follows from Eq. (75), immediately giveṡ The thermodynamic quantities must also satisfy the Gibbs law, which can be formulated as [94] nT d s where T is the thermodynamic temperature of the system.

Thermodynamic quantities in f (R, T ) gravity
After some simple algebraic manipulations the energy balance Eq. (93) can be reformulated aṡ where the creation pressure P c is defined as where we have denoted w = P/ρ. Then the generalized balance equation (93) can be derived from the divergence of the total energy momentum tensor T μν , defined as On the other hand under the assumption of adiabatic particle production, withσ = 0, the Gibbs law giveṡ which together with the energy balance equation gives immediately the relation between the particle creation rate and the creation pressure as In the framework of f (R, T ) gravity theory we obtain for the particle creation rate the general expression Hence the condition ≥ 0 imposes a strong constraint on the physical parameters of the theory. In the case of pressureless dust, with P = 0, w = 0, under the assumption f T > 0, we obtain the following general cosmological constraint that must be satisfied by the function f T for all times: The divergence of the entropy flux vector can be reformulated in terms of the creation pressure as Finally, we consider the temperature evolution in a system with particle creation. In order to fully determine the time behavior of a relativistic fluid we must add two equations of state for the density and pressure, which have the general form ρ = ρ(n, T ) and p = p(n, T ), respectively. Then we obtaiṅ By using the energy and particle balance equations we find With the use of the thermodynamic identity [96] T Equation (89) yields the temperature evolution of a relativistic fluid in the presence of matter creation: In the particular case (∂ P/∂ρ) n = w = constant, we obtain for the temperature-particle number dependence the simple expression T ∼ n w .

The case w = −1
Based on the homogeneous and isotropic Friedmann-Robertson-Walker metric, and on the energy conservation equationρ + 3H (1 + w)ρ = 0, in [97], general cosmological thermodynamic properties with an arbitrary, varying equation-of-state parameter w(a), where a is the scale factor, were discussed. The w = − 1-crossing problem of w was explicitly pointed out, and the behaviors of the quantities (ρ(a), μ(a), T (a), etc.) at/near w = −1 were discussed. As a result of this study it was concluded that all cosmological quantities must be regular and well defined for all values of w(a), and indeed they are [97]. In the present thermodynamical approach we have assumed that matter is created in an ordinary form, and therefore all our previous results are valid for w ≥ 0. However, the thermodynamic approach and the interpretation of f (R, T ) gravity can be extended to the case w < 0. In the following we will consider this problem, and we show that our results are valid, in the sense of regularity and of being well defined, even in the case of w = −1.
In particular, we concentrate on the temperature evolution equation, which still holds even if w = P/ρ = −1. The demonstration is as follows. First we consider the perfect-fluid energymomentum balance equation, When w = −1, Eq. (93) becomeṡ where P c is the matter creation pressure. Under the assumption of adiabatic particle production, withσ = 0, the Gibbs law giveṡ where σ is the entropy per particle. That is, from the above two equations, we obtaiṅ Since ρ = ρ (n, T ), we havė Combining the above equation and the thermodynamic identity [96] T it immediately follows that Eq. (92) still holds even for negative values of w = − 1. If w = −1 = constant, we can find from Eq. (92) that nT is a constant, or T ∼ 1/n. This relation shows that for very low density "dark energy" particles their thermodynamic temperature is extremely high, while high particle number (density) systems have a very low temperature. In the limit n → ∞, the temperature of the "dark energy" made system tends to zero.

Cosmology of Palatini f (R, T ) gravity
In the present section we investigate the cosmological implications of the Palatini formulation of f (R, T ) gravity. We assume that the Universe is flat, homogeneous and isotropic, with the metric given in comoving coordinates by the Friedmann-Robertson-Walker metric, where a(t) is the scale factor. We also introduce the Hubble function, defined as H =ȧ/a. We assume that the matter content of the Universe consists of a perfect fluid that can be characterized by two thermodynamic parameters, the energy density ρ, and the pressure P, respectively. As for the relations of the geometric quantities in the g andg frames, their detailed computation is presented in Appendix C.

Generalized Friedmann equations in Palatini f (R, T ) gravity
Substituting the expression of the perfect-fluid energymomentum tensor as well as Eqs. (61), (62) and (C11) into the Palatini field Eq. (33), from the 00 component we obtain the first modified Friedmann equation as Similarly, with the help of Eqs. (61), (62) and (C12), we can obtain the second modified Friedmann equation from the 'ii' components of Eq. (33), Note that the first modified Friedmann equation (100) can be written more compactly as Finally, we substitute Eqs. (61) and (C10) into Eq. (34) and obtain the trace equation By eliminating the term 3H 2 between the two generalized Friedmann equations we obtain the evolution equation for the Hubble function as given by The f R g,˜ , T → f R g,˜ limits of the two modified Friedmann equations can be given by and If we go one step further and take f R g,˜ → R g,˜ , then f = R(g) = 6 Ḣ + 2H 2 and the above two equations reduce to and That is, and Hence the first modified Palatini f (R, T ) Friedmann equation (100) reduces to the ordinary second Friedmann equation (109) when f (R, T ) → R.

The energy balance equation
With the help of Eq. (C2), we can directly work out the covariant divergence of the energy-momentum tensor (4) as and Substituting the above two equations and Eqs. (61) and (62) into the already known Palatini f (R, T ) energy balance equation (93), we obtaiṅ When f (R, T ) → f (R), the energy balance equation reduces to the ordinary conservation equation,

Deceleration parameter and equation of state of the Universe
An important cosmological parameter, indicating the accelerating/decelerating nature of the cosmological dynamics is the deceleration parameter q, defined to be Using Eqs. (100) and (104), we immediately obtain For a vacuum Universe with ρ = P = 0, the condition for accelerated expansion, q < 0, reduces tö The deceleration parameter can also be defined, by analogy with the standard general relativistic cosmology, in terms of the effective parameter w eff of the equation of state of the Universe as giving By using the above definition we obtain for the effective parameter of the equation of state of the Universe the expression with the general solution given by where ρ 0 is an integration constant, and we have denoted Equation (104) becomes In the limit of large times the last, exponential term in the above equation becomes negligibly small, and hence we can approximate the solution of Eq. (125) by where F 0 is an arbitrary constant of integration, and without any loss of generality we have taken the second integration constant as zero. Then the first generalized Friedmann equation (100) gives the Lagrangian function of the model as On the other hand the trace of Eq. (103) gives (note that in the following R = R(g)) or Equations (127) and (129) give a parametric representation of f as a function of R, with t taken as a parameter. Once the function t = t (R) is obtained from Eq. (129), by substituting it in Eq. (127) we can find the explicit dependence of f on R.

Comparison with the metric f (R, T ) cosmology
Using the Friedmann-Robertson-Walker metric, and the field equation (9), we can similarly derive the two Friedmann equations in the metric formulation. The cosmological equations in metric f (R, T ) gravity are different from their counterparts in the Palatini formalism, due to the presence of some dynamical terms related to f R , and they are given by Besides, the deceleration parameter q = −Ḣ /H 2 − 1 can be obtained: For ρ = P = 0, an accelerated Universe with q < 0 requires thaẗ The condition for accelerated expansion in the metric f (R, T ) gravity is very different from the similar condition given by Eq. (118), in the Palatini formulation of the theory. The presence of the extra term, 3Ḟ 2 /2F, in Eq. (118) may have a significant effect on the transition from the decelerating to the accelerating phase. In the Palatini formulation the moment of the transition to the accelerated expansion with q ≤ 0 is determined by the equation with f R = f R R g,˜ , T , while in the metric f (R, T ) gravity the same condition is given by the much simpler expression with f = f (R(g), T ). It is interesting to note that the conditions for the transition to an accelerated expansion in the vacuum case are independent in both approaches on f T . However, the Palatini formulation of f (R, T ) gravity allows for a much richer cosmological dynamics, as compared to the metric formulation.
As for the energy-momentum balance equation and the non-geodesic equation of motion of massive test particles: since they are all derived from the divergence of the energymomentum tensor, and since ∇ μ T μ ν is the same in the two formalisms, independently of the functional form of f , the energy balance equations and the equations of motion have the same functional form in both approaches.

Specific cosmological models in the Palatini f (R, T ) gravity
In the present section we will investigate two specific cosmological models in the framework of the Palatini formulation of f (R, T ) gravity. We will assume that the action of the gravitational field has the general form where β is a constant, and k R g,˜ and g(T ) are arbitrary functions of the Ricci scalar and the trace of the matter energy-momentum tensor. For the function g(T ) for simplicity we will assume a simple linear dependence on T , so that g(T ) = T . As for the function k R g,˜ , we will consider two cases, corresponding to the Starobinsky model k(R) ∼ R 2 g,˜ [98], and the 1/R g,˜ case, respectively.
We consider a Palatini f (R, T ) model, specified by a functional form of f (R, T ) given as where α, β are constants, g(T ) is a function that depends on T solely, and for simplicity we set g(T ) = T . For this Lagrangian we immediately find Moreover, we assume that the cosmological fluid satisfies a linear barotropic equation of state of the form P = wρ, w = constant. Consequently, from the trace of (15) of the Palatini field equations we first obtain thus When β → 0, β 0 → 1 − 3w, and R g,˜ → −8π T , respectively. Substituting Eq. (139) back into the expressions of f (R, T ) and F, we obtain and Substitution of the expression of f T into the balance equation (114) giveṡ where we have denoted When β → 0, then β 1 = 3(1 + w). Substituting Eq. (143) into Eq. (142), we find Taking into account the limit αρ → 0, we have the series expansion where we have denoted Equation (143) can be immediately integrated to give Hence Eq. (147) becomes a first-order differential equation, with the general solution given by In the limit of large times we obtain a(t) ∼ t 2/β 1 , and H (t) = (2/β 1 ) (1/t), respectively. The deceleration parameter in this model is given by q = β 1 /2 − 1, and once the coefficient β 1 satisfies the condition β 1 /2 < 1, the Universe will experience an accelerated evolution. For an arbitrary w, the condition for a power law type accelerated expansion is β < − (1 + 3w) /4w, a condition which shows that, for w > 0, β must take negative values.
Now let us consider the following f (R, T ) gravity Palatini type model: which immediately gives F = 1 + α 2 R −2 g,˜ /3. From the trace of (15) of the Palatini field equations we obtain where we have denoted The algebraic equation (152) has two distinct solutions. However, only one of them can be adopted as the physical solution, more exactly the one which in the limit f → R, would give R = −8π T , which is the trace of the Einstein field equation.
Hence when ≤ 0, the physical solution of Eq. (152) is [99] R g,˜ = − In the following we will study this cosmological model under the approximation α | | [99], which allows us to find the evolution of the Universe at later times when R g,˜ is relatively small. Under the adopted approximation we can expand R g,˜ to the order of O (| |), and we obtain For simplicity, we set again g(T ) = T . Using the approximations (155), (156) and (157), respectively, we can easily obtain the Palatini f (R, T ) field equations for this model: When f (R, T ) → f (R), β → 0, and the above field equations reduce to the field equations considered within the Palatini formulation of f (R) gravity, considered in [99], 2 : With the equation of state w = P/ρ, and since T = −ρ + 3P for a perfect fluid, we have = −8πβ 0 ρ; besides, we already know from Eq. (143) thatρ = −β 1 Hρ. Thus, similar to Eq. (145), we have In the first-order approximation in ρ/α we obtain Hence in the present model a cosmological constant α/12 is automatically generated, due to the 1/R modification of the gravitational Lagrangian. Hence for ρ → 0, the Universe will end in a de Sitter phase, with H = H 0 = √ α/12 = constant. For ρ = 0, we have ρ = ρ 0 a −β 1 , and the evolution of the scale factor is determined by the equatioṅ where we have denoted β 4 = β 0 (5 − β 1 )/6 + (3 + 4β)w + 1 π , with the general solution given by where we have used the initial condition a (t 0 ) = a 0 . For the Hubble function we obtain while the deceleration parameter of this model is given by where sech t = 1/cosh t. In the limit t → ∞, q → −1, and hence the Universe ends in a de Sitter type accelerating phase, independent of the matter equation of state.
a. The matter dominated phase The matter dominated phase corresponds to the choice w = 0 in the matter equation of state, that is, to a Universe filled with pressureless baryonic matter. In order to investigate the behavior of the cosmological model during matter domination, we consider the series expansion of the cosmological parameters. Thus we obtain These equations describe the main cosmological parameters during the matter dominated era. The expansion is decelerating, and, depending on the model parameters, the deceleration parameter can have a large range of positive values. The transition to the accelerating phase occurs at a time interval t tr , which in the first-order approximation is found as Since t tr must be greater than t 0 , it follows that in order for the model to admit a matter dominated era followed by a transition to an accelerated phase, the model parameters must satisfy the condition 6 (β 1 − 2) β 4 ρ 0 − αa β 1 0 > 0, or, equivalently, a condition that can easily be satisfied by appropriately choosing the free parameters in the gravitational action.

Discussions and final remarks
In the present paper we have considered in the framework of the Palatini formalism the gravitational field equations for the modified gravity f (R, T ) theory, implying a geometrymatter coupling, with the trace of the energy-momentum tensor included as a field variable in the gravitational action. We have derived the field equations by independently varying the metric and the connection in the f (R, T ) type gravitational action, and we have formulated them in both the initial metric frame and the conformal one, in which the independent connection can be expressed as the Levi-Civita connection of an auxiliary, energy-momentum trace dependent metric, which is related to the physical metric by a conformal transformation. Similar to the metric case [46] the energymomentum tensor of the matter is not conserved, and the energy and momentum balance equations take the same form as in the metric theory. Generally, Palatini type theories have a number of special properties that make them especially attractive for analyzing strong gravity phenomena, like, for example, the dynamics of the early Universe or stellar collapse processes [100][101][102][103][104][105][106][107]. The coupling of the trace of the energy-momentum tensor with the curvature scalar generates some extra terms in the gravitational field equations, which strongly depend on the possible functional forms for the geometry-energy momentum trace coupling. If, for example, it would be possible to generate through the geometryenergy momentum trace coupling some repulsive forces, then one could obtain cosmological models that are non-singular at extremely high densities and high geometric curvatures, or one could even construct models for non-singular collapsing stars as viable alternatives for the black hole paradigm.
To obtain such repulsive gravitational forces, in modified gravity theories with geometry-matter coupling no new degrees of freedom in the matter side (exotic sources) or in the gravitational side are required in the total action. In these models the extra force is simply induced by the coupling between matter and geometry. Our present results show that Palatini type theories might play an important role in the phenomenology of gravity at both high densities (energies), as well as in the very low density limit. On the other hand, in the variational process the assumption of the independence between metric and connection is essential to obtain secondorder differential equations for the metric tensor. It is thus possible to assume that at large/small scales the effective descriptions of the gravitational forces, going beyond standard general relativity, could come from the Palatini formulation of gravity theories.
In the Palatini type formulation of f (R, T ) gravity, the equation of motion of massive particles is non-geodesic, and in three dimensions and in the Newtonian limit, Eq. (66) can be formally represented as an ordinary vector equation in three dimensions of the form a = a N + a f where a represents the total acceleration of the particle, a N denotes the Newtonian gravitational acceleration, while a f is the acceleration due to the presence of the extra force induced by the coupling between geometry and matter. This shows that one observational or experimental possibility of testing the effects of the coupling between geometry and the trace of the energy-momentum tensor could be in the physical domain of extremely small particle accelerations, with values of the order of 10 −10 m/s 2 . Such an acceleration could explain the observed behavior of test particles rotating around galaxies, which is usually explained by postulating the existence of dark matter. However, as a possible astrophysical application of the gravitational field equations derived with the Palatini formalism one may consider an alternative view to the dark matter problem, in which the mass discrepancy in galaxies and clusters of galaxies as well as the galactic rotation curves are explained by the existence of a non-minimal coupling between matter and geometry.
We have also briefly investigated the intriguing feature of the non-conservation of the energy-momentum tensor of the matter in f (R, T ) gravity theory by interpreting it in the framework of the thermodynamics of open systems. We have interpreted this effect as describing phenomenologically the particle production in the cosmological fluid filling the Universe, with the extra terms induced by the nonminimal coupling between R and T assumed to describe particle creation processes, with the gravitational field acting as a source for particles. We have explicitly obtained the particle creation rates, the entropy flux, the creation pressure and entropy generation rate in a covariant form, as functions of the Lagrangian density f (R, T ) of the theory, and of its derivatives, respectively. On the other hand it is natural to assume that such particle production processes are of the same nature as the similar processes that appear in the framework of quantum field theory in curved space-times. A static gravitational field does not produce particles. But a time dependent gravitational field can generate new particles. This interesting analogy between gravitational theories with geometry-matter coupling and quantum field theory in curved space-times may open the possibility of an effective classical description of quantum gravity on small geometric scales.
As a cosmological application of the Palatini formalism of f (R, T ) theory we have briefly considered two classes of cosmological models, corresponding to two choices of the gravitational Lagrange density f (R, T ) = k(R) + g(T ), with k(R) ∼ R 2 and k(R) ∼ 1/R, respectively. In both cases we have assumed for the function g(T ) the simple form g(T ) ∼ T . We have explicitly shown that both models δ ∇ μ g νλ = δ ∇ ν g λμ = δ ∇ λ g μν = 0. A variation δg μν of the metric tensor then leads to δT = δ g αβ T αβ = T αβ δg αβ δg μν + μν δg μν = T μν + μν δg μν , μν ≡ g αβ δT αβ δg μν (A3) and respectively, where the above relation is called the Palatini identity. Then for the variation of the metric we obtain Assuming that the variation of δg μν vanishes at infinity, then and respectively.