Cosmological implications of modified gravity induced by quantum metric fluctuations

We investigate the cosmological implications of modified gravities induced by the quantum fluctuations of the gravitational metric. If the metric can be decomposed as the sum of the classical and of a fluctuating part, of quantum origin, then the corresponding Einstein quantum gravity generates at the classical level modified gravity models with a non-minimal coupling between geometry and matter. As a first step in our study, after assuming that the expectation value of the quantum correction can be generally expressed in terms of an arbitrary second order tensor constructed from the metric and from the thermodynamic quantities characterizing the matter content of the Universe, we derive the (classical) gravitational field equations in their general form. We analyze in detail the cosmological models obtained by assuming that the quantum correction tensor is given by the coupling of a scalar field and of a scalar function to the metric tensor, and by a term proportional to the matter energy-momentum tensor. For each considered model we obtain the gravitational field equations, and the generalized Friedmann equations for the case of a flat homogeneous and isotropic geometry. In some of these models the divergence of the matter energy-momentum tensor is non-zero, indicating a process of matter creation, which corresponds to an irreversible energy flow from the gravitational field to the matter fluid, and which is direct consequence of the non-minimal curvature-matter coupling. The cosmological evolution equations of these modified gravity models induced by the quantum fluctuations of the metric are investigated in detail by using both analytical and numerical methods, and it is shown that a large variety of cosmological models can be constructed, which, depending on the numerical values of the model parameters, can exhibit both accelerating and decelerating behaviors.


Introduction
Modified gravity theories may provide an attractive alternative to the standard explanations of the present day observations that have shaken the well-established foundations of theoretical physics. Astronomical observations have confirmed that our Universe does not according to standard general relativity, as derived from the Hilbert-Einstein action, S = −R/2κ 2 + L m √ −gd 4 x, where R is the Ricci scalar, κ is the gravitational coupling constant, and L m is the matter Lagrangian, respectively. Extremely successful on the Solar System scale, somehow unexpectedly, general relativity faces, on a fundamental theoretical level, two important challenges, the dark energy and the dark matter problems, respectively. Several high precision astronomical observations, with the initial goal of improving the numerical values of the basic cosmological parameters by using the properties of the distant Type Ia Supernovae, have provided the result that the Universe underwent recently a transition to an accelerating, de Sitter type phase [1][2][3][4][5]. The necessity of explaining the late time acceleration lead to the formulation of a new paradigm in theoretical physics and cosmology, which postulates that the explanation of the late time acceleration is the existence of a mysterious component of the Universe, called dark energy (DE), which can describe the late time dynamics of the Universe [6,7], and can explain all the observed features of the recent (and future) cosmological evolution. However, in order to close the matter-energy balance of the Universe, a second, and equally mysterious component, called Dark Matter, is required. Dark Matter, assumed to be a nonbaryonic and non-relativistic (cold) component of the Universe, is necessary for explaining the dynamics of the hydrogen clouds rotating around galaxies, and having flat rotation curves, as well as the virial mass discrepancy in clusters of galaxies [8,9]. The direct detection/observation of the dark matter is extremely difficult due to the fact that it interacts only gravitationally with the baryonic matter. After many decades of intensive observational and experimental efforts there is no direct evidence as regards the particle nature of the dark matter.
One of the best theoretical descriptions that fits almost perfectly the observational data is based on the simplest theoretical extension of general relativity, which includes in the gravitational field equations the cosmological constant [10,11]. Based on this theoretical formalism the basic paradigm of modern cosmology has been formulated as the CDM-Cold Dark Matter Model, in which the dark energy is nothing but the simple constant introduced almost one hundred years ago by Einstein. Even that the CDM model fits the data well, it raises some fundamental theoretical questions about the possibility of explaining it. There is no theoretical explanation for the physical/geometrical nature of , and, moreover, general relativity cannot give any hints on why it is so small, and why it is so fine tuned [10,11]. Therefore, the possibility that dark energy can be explained as an intrinsic property of a generalized gravity theory, going beyond general relativity, and its Hilbert-Einstein gravitational action, cannot be rejected a priori. In this context a large number of modified gravity models, all trying to extend and generalize the standard Einsteinian theory of gravity, have been proposed. Historically, in going beyond Einstein gravity, the first, and most natural step, was to extend the geometric part of the Hilbert-Einstein action. One of the first attempts in this direction is represented by the f (R) gravity theory, in which the gravitational action is generalized to be an arbitrary function of the Ricci scalar, so that S = 1 2κ 2 f (R) √ −gd 4 x + L m √ −gd 4 x [12][13][14][15][16][17][18]. However, this as well as several other modifications of the Hilbert-Einstein action focused only on the geometric part of the gravitational action, by explicitly postulating that the matter Lagrangian plays a subordinate and passive role only, as compared to the geometry [19]. From a technical point of view such an approach implies a minimal coupling between matter and geometry. But a fundamental theoretical principle, forbidding an arbitrary coupling between matter and geometry, has not been formulated yet, and perhaps it may simply not exist. On the other hand, if general matter-geometry couplings are introduced, a large number of theoretical gravitational models, with extremely interesting physical and cosmological properties, can easily be constructed. The first of the modified gravity theory with arbitrary geometry-matter coupling was the f (R, L m ) modified grav-ity theory [20][21][22][23], in which the total gravitational action takes the form S = 1 2κ 2 f (R, L m ) √ −gd 4 x. In this kind of theories matter is essentially indistinguishable from geometry, and plays an active role in generating the geometrical properties of the spacetime. A different geometrymatter coupling is introduced in the f (R, T ) [24,25] gravity theory, where matter and geometry are coupled via T , the trace of the energy-momentum tensor. The gravitational action of the f (R, T ) theory is given by S = f (R, T )/2κ 2 + L m √ −gd 4 x. A recent review of the generalized f (R, L m ) and f (R, T ) type gravitational theories with non-minimal curvature-matter couplings can be found in [26]. Several other gravitational theories involving geometry-matter couplings have also been proposed, and extensively studied, like, for example, the Weyl-Cartan-Weitzenböck (WCW) gravity theory [27], hybrid metric-Palatini f (R, R) gravity [28][29][30], where R is the Ricci scalar formed from a connection independent of the metric, f R, T, R μν T μν gravity theory, where R μν is the Ricci tensor, and T μν the matter energy-momentum tensor, respectively [31,32], or f (T , T ) gravity [33], in which a coupling between the torsion scalarT , essentially a geometric quantity, and the trace T of the matter energy-momentum tensor is introduced. Gravitational models with higher derivative matter fields were investigated in [34].
One of the interesting (and intriguing) properties of the gravitational theories with geometry-matter coupling is the non-conservation of the matter energy-momentum tensor, whose four-divergence is usually different of zero, ∇ μ T μν = 0. This property can be interpreted, from a thermodynamic point of view, by using the formalism of open thermodynamic systems [25]. Hence one can assume that the generalized conservation equations in these gravitational theories describe irreversible matter creation processes. Thus the non-conservation of the energy-momentum tensor describes an irreversible energy flow from the gravitational field to the newly created matter constituents, with the second law of thermodynamics requiring that spacetime transforms into matter. In [25] the equivalent particle number creation rates, the creation pressure and the entropy production rates were obtained for both f (R, L m ) and f (R, T ) gravity theories. The temperature evolution laws of the newly created particles was also obtained, and studied. Due to the non-conservation of the energy-momentum tensor, which is a direct consequence of the geometry-matter coupling, during the cosmological evolution of the Universe a large amount of comoving entropy could also be produced.
The prediction of the production of particles from the cosmological vacuum is one of the remarkable results of the quantum field theory in curved spacetimes [35][36][37][38][39]. Particles 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 [39]. Therefore, the presence of particle creation processes in both quantum theories of gravity and modified gravity theories with geometry-matter coupling may suggest that a deep connection between these two, apparently very different physical theories, may exist. And, interestingly enough, such a connection has been found in [40], where it was pointed out that by using a non-perturbative approach for the quantization of the metric, proposed in [41][42][43], a particular type of f (R, T ) gravity, with Lagrangian given by L = where α is a constant, naturally emerges as a result of the quantum fluctuations of the metric. This result suggests that an equivalent microscopic quantum description of the matter creation processes in f (R, T ) or f (R, L m ) gravity is possible, and such a description could shed some light on the physical mechanisms leading to particle generation via gravity and matter geometry coupling. Such mechanisms do indeed exist, and they can be understood, at least qualitatively, in the framework of some quantum/semi-classical gravity models. It is the goal of the present paper to further investigate the cosmological implications of modified gravities induced by the quantum fluctuations of the gravitational metric, as initiated in [40][41][42][43]. As a starting point we assume that a general quantum metric can be decomposed as the sum of the classical and of a fluctuating part, the latter being of quantum (or stochastic) origin. If such a decomposition is possible, the corresponding Einstein quantum gravity generates at the classical level modified gravity models with a non-minimal interaction between geometry and matter, as previously considered in [20,24,31]. After assuming that the expectation value of the quantum correction can be generally expressed in terms of an arbitrary second order tensor K μν , which can be constructed from the metric and from the thermodynamic quantities characterizing the matter content of the Universe, we derive from the first order quantum gravitational action the (classical) gravitational field equations in their general form. We analyze in detail the cosmological models obtained from the quantum fluctuations of the metric in tqo cases. First we assume that the quantum correction tensor K μν is given by the coupling of a scalar field and of a scalar function to the metric tensor, respectively. As a second case we consider that K μν is given by a term proportional to the matter energymomentum tensor. The first choice gives a particular version of the f (R, T ) gravity model [24], while the second choice corresponds to specific case of the modified gravity theory of the form f R, T, R μν T μν , T μν T μν [31]. For each considered model we obtain the gravitational field equations, and the generalized Friedmann equations for the case of a flat homogeneous and isotropic geometry. In some of these models the divergence of the matter energy-momentum ten-sor is non-zero, indicating a process of matter creation. From a physical point of view a non-zero divergence of the energymomentum tensor can be interpreted as corresponding to an irreversible energy flow from the gravitational field to the matter fluid. Such an irreversible thermodynamic process is the direct consequence of the non-minimal curvature-matter coupling, induced in the present case by the quantum fluctuations of the metric [25,44,45]. The cosmological evolution equations of these modified gravity models induced by the quantum fluctuations of the metric are investigated in detail by using both analytical and numerical methods. As a result of this analysis we show that a large variety of cosmological models can be constructed. Depending on the numerical values of the model parameters, these cosmological models can exhibit both late time accelerating, or decelerating behaviors.
The present paper is organized as follows. In Sect. 2 we derive the general set of field equations induced by the quantum fluctuations of the metric. The relation between this approach and the standard semi-classical formulation of quantum gravity is also briefly discussed. In Sect. 3 we investigate in detail the cosmological implications of the quantum fluctuations induced modified gravity models with the fluctuation tensor proportional to the metric. Two cases are considered, in which the fluctuation couples to the metric via a scalar field, and a scalar function, respectively. Modified gravity models induced by quantum metric fluctuations proportional to the energy-momentum tensor are investigate in Sect. 4. Finally, we discuss and conclude our results in Sect. 5. The details of the derivation of the gravitational field equations for an arbitrary metric fluctuation tensor and for a fluctuation tensor proportional to the matter energy-momentum tensor are presented in Appendices A and B, respectively. In the present paper we use a system of units with c = 1.

Modified gravity from quantum metric fluctuations
In the present section we will briefly review the fluctuating metric approach to quantum gravity, we will discuss its relation with standard semi-classical gravity, and we will point out the quantum mechanical origins of the modified gravity models with geometry-matter coupling. Moreover, we derive the general set of field equations induced by the quantum fluctuations of the metric for an arbitrary form of the tensor K μν .

Modified gravity as the semi-classical approximation of quantum gravity
In the standard quantum mechanics physical (or geometrical) quantities must be represented by operators. Therefore, in a full non-perturbative quantum approach the Einstein gravi-tational field equations must take an operator form, given by [41][42][43] As shown in [41][42][43], in order to extract meaningful physical information from the Einstein operator equations one must average it over all possible products of the metric operatorŝ g (x 1 ) . . .ĝ (x n ), and thus to solve the infinite set of equations for the Green functionsĜ μν . In the above equations |Q > is quantum state that might not be the ordinary vacuum state. These equations cannot be solved analytically, and hence we have to use some approximations [41][42][43]. In [41] it was suggested to decompose the metric operator into the sum of an average metric g μν , and a fluctuating part δĝ μν , according tô Assuming that where K μν is a classical tensor quantity, and ignoring higher order fluctuations, the gravitational Lagrangian will be modified into [41] L = L g ĝ μν + √ −gL m ĝ μν ≈ L g + δL g δg μν δĝ μν where κ 2 = 8π G/c 4 , and where we have defined the matter energy-momentum tensor as Even that the above formalism starts from a full quantum approach of gravity, after performing the decomposition of the metric we are still considering semi-classical theories. In this paper we will consider several functional forms of the so called quantum perturbation tensor K μν , we will obtain the field equations of the corresponding gravity theory, and we will investigate their cosmological implications, respectively.
The gravitational field equations corresponding to the first order corrected quantum Lagrangian (4) are given, in a general form, by where K = g μν K μν and A αβ δ K αβ = δg μν (γ αβ μν A αβ ). Here A αβ is either R αβ or T αβ , and γ αβ μν is an algebraic tensor, an operator, or the combination of them. The detailed derivation of Eq. (6) is presented in Appendix A.
It is interesting to compare the formalism based on the quantum metric fluctuation proposal to the standard semiclassical gravity approach, which also leads to particle production via the non-conservation of the matter energymomentum tensor [45]. Semi-classical gravity is constructed from the basic assumption that the gravitational field is still classical, while the classical matter (bosonic) fields φ are taken as quantized. The direct coupling of the quantized matter fields to the classical gravitational fields is performed via the replacement of the quantum energy momentum tensor T μν by its expectation value T μν , obtained with respect to some quantum state . Therefore the effective semi-classical Einstein equation can be postulated as [46], From Eq. (7) it follows immediately that the classical energy-momentum tensor T μν of the gravitating system is obtained from its quantum counterpart through the definition |T μν | = T μν . The semi-classical Einstein equations (7) can be derived from the variational principle [47] where x is the standard general relativistic classical action of the gravitational field, while the quantum part of the action is given by In Eq. (9)Ĥ is the Hamiltonian operator of the gravitating system, and α is a Lagrange multiplier. The variation of Eq. (8) with respect to the wave function leads to the normalization condition for the quantum wave function | = 1, to the Schödinger equation for , and to the semi-classical Einstein equations (7), respectively. Hence, in this simple phenomenological approach to semi-classical gravity, the Bianchi identities still require the conservation of the effective energy-momentum tensor, ∇ μ |T μν | = 0. A more general set of semi-classical Einstein equations, having many similarities with the modified gravity models, can be derived by introducing a new coupling between the quantum fields and the classical curvature scalar of the spacetime. One such model was proposed in [47], where the total action containing the geometry-quantum matter coupling term was assumed to be of the form where F and f are arbitrary functions, and . Such a geometry-matter coupling term modifies the Hamiltonian H (t) in the Schrödinger equation (10) into [47] where N is the lapse function, while ξ i are intrinsic spacetime coordinates, chosen in such a way that the normal vector to a space-like surface is time-like on the entire spacetime manifold. The scalar function ς can be obtained: ς = det ς rs , where ς rs is the metric induced on a space-like surface σ (t).
The surface σ (t) globally slices the spacetime manifold into space-like surfaces. Then, by taking into account the effect of the geometry-quantum matter coupling, the effective semiclassical Einstein equations become [47] In the semi-classical gravitational model introduced through Eq. (13), the matter energy-momentum tensor is not conserved anymore, since ∇ μ T μν = 0. Thus, describes an effective particle generation process, in which there is a quantum-mechanically induced energy transfer from spacetime to matter. Equation (13) also gives an effective semiclassical description of the quantum processes in a gravitational field, which are intrinsically related to the matter and energy non-conservation.
It is interesting to compare Eqs. (6) and (13), both based on some assumptions on the quantum nature of gravity. While Eq. (6) is derived through a first order approximation to quantum gravity, Eq. (13) postulates the existence of a quantum coupling between geometry and matter. While the coupling in Eq. (13) is introduced via a scalar function, the quantum effects are introduced in Eq. (6) through the fluctuations of the metric, having a tensor algebraic structure. However, in both models, the matter energy-momentum tensor is generally not conserved, indicating the possibility of the energy transfer between geometry and matter. Therefore, the physical origin of the modified gravity theories with geometrymatter coupling, which "automatically" leads to matter creation processes, may be traced back to the semi-classical approximation of the quantum field theory in a Riemannian curved spacetime geometry.

The cosmological model
For cosmological applications we adopt the Friedmann-Robertson-Walker metric, where a(t) is the scale factor. The components of the Ricci tensor for this metric are We define the Hubble function as H =ȧ/a. As an indicator of the possible accelerated expansion we consider the deceleration parameter q, defined as Negative values of q indicate accelerating evolution, while positive ones correspond to decelerating expansion. In order to perform the comparison between the observational and theoretical results, instead of the time variable t we introduce the redshift z, defined as where we have adopted for the scale factor a (z) the normalization a(t 0 ) = 1, where t 0 is the present age of the Universe. Hence as a function of the redshift the time derivative operator can be expressed as 3 Modified gravity from quantum perturbation of the metric proportional to the classical metric, As a first example of cosmological evolution in modified gravity models induced by the quantum fluctuations of the metric we consider the simple case in which the expectation value of the quantum fluctuation tensor is proportional to the classical metric [40,41], where α(x) is an arbitrary function of the spacetime coor- [40,41], respectively. This condition suggest that we have an additional part of the metric, due to the quantum perturbation effects, which is proportional to the classical one. In the following we will consider two distinct cases, by assuming first than α(x) is a scalar field, with a specific selfinteraction potential. As a second model we will assume that α(x) is a simple scalar function.

Scalar field-metric coupling
If α(x) is a scalar field, we add an additional Lagrangian as the matter source into the general quantum perturbed Lagrangian (4). Hence we obtain the Lagrangian of our model as By taking the variation of the Lagrangian with respect to the field α, the Euler-Lagrange equation gives the generalized Klein-Gordon equation for the scalar field, which has the form By taking the variation of Eq. (21) with respect to the metric tensor we obtain the Einstein gravitational field equations as where After contraction of the Einstein field equations we obtain (25) and thus we can reformulate the gravitational field equations as The divergence of the matter energy-momentum tensor can be obtained from Eq. (23), and is given by In the following we will adopt for the matter energymomentum tensor the perfect fluid form where ρ is the matter energy density, p is the thermodynamic pressure, and u μ is the matter four velocity, satisfying the normalization condition u μ u μ = 1. For this choice of the energy-momentum tensor we have To obtain the above equation we have adopted for the matter Lagrangian the representation L m = p. For the scalar θ we obtain With the use of Eq. (25) the Klein-Gordon equation now reads

Cosmological applications
In the following we will restrict our analysis to the case of the flat Friedmann-Robertson-Walker metrics, and hence we set K = 0 in the gravitational field equations. The Friedmann and the Klein-Gordon equations describing our generalized gravity model obtained from a fluctuating metric take the form where a. The energy conservation equation By multiplying Eq. (32) with a 3 , and taking the time derivative of its both sides, we obtain With the use of Eq. (33) we obtain Equation (37) can be rewritten in the equivalent form The same conservation equation can be obtained directly from the field equations (27). By taking into account the explicit form of the matter energy-momentum tensor as given by Eq. (28), and that we have (40) and respectively. Hence Eq. (28) takes the form By taking into account Eq. (81), after eliminating α from the above equation, we reobtain again the conservation equation Eq. (38).

b. The dimensionless form of the cosmological evolution equation
In order to simplify the mathematical formulation of the cosmological model we rescale first the field α and its potential V as α → κα, and V → κ 2 V , respectively. Next, we introduce a set of dimensionless variables (τ, h, r, P, v), defined as where H 0 is the present day value of the Hubble function. Then the cosmological evolution equation take the dimensionless form In order to close the system of Eqs. (44)-46) we must specify the equation of state of the matter P = P(r ), and the functional form of the self-interaction potential of the scalar field v. In the following we will restrict our analysis to the case of the dust, with P = 0. Then, by denoting u = dα/dτ , from Eq. (44) we obtain for the dimensionless energy density r the expression Then, by introducing the redshift z as an independent variable, it follows that the cosmological evolution is described by the following system of equations: For the deceleration parameter we obtain the expression For the scalar field self-interaction potential we adopt a Higgs type form, so that where μ 2 > 0 and λ > 0 are constants. In the following we will consider two cases. In the first case we assume that the Universe dominated by the quantum fluctuations of the metric evolves in the minimum of the Higgs potential, so that ∂v/∂α = 0, implying α = ± μ 2 /λ, and v(α) = μ 4 /4λ = constant. Secondly, we will investigate the evolution of the Universe in the presence of a time varying "full" Higgs type scalar field self-interaction potential (52). Once the form of the potential is fixed, the system of differential equations Eqs.  Fig. 4, shows a complex behavior. The Universe starts its evolution at a redshift z = 1 in a decelerating phase, with q ≈ 1. The Universe begins to accelerates, and it enters in a marginally accelerating phase, with q = 0, at a redshift of around z ≈ 0.3. The variation of the deceleration parameter strongly depends on the numerical values of v 0 , and, depending on this numerical value, the evolution of the Universe at the present time can either be de Sitter, with q = −1, or have higher values of q, of the order of q ≈ −0.5. The evolution of the cosmological and physical parameters of a Universe filled with a scalar field with Higgs potential (52), coupled to the fluctuating quantum metric, are presented in Figs. 5, 6, 7, and 8. To numerically integrate the gravitational field equations (48)-(50) in the redshift range z ∈ [0, 2] we have adopted the initial conditions α(0) = 0.1 and u(0) = 0.1, respectively. We have fixed the value of the coefficient λ in the Higgs potential as λ = 150, and we have varied the mass μ 2 > 0 of the Higgs field.

d. Cosmological evolution in the presence of the Higgs potential
The Hubble function, plotted in Fig. 5, is an increasing function of the redshift, indicating an expansionary evolution. However, h has a complex behavior, which is strongly dependent on the numerical value of μ 2 . The variation with z of the Higgs potential is represented in Fig. 6. The potential has a damped harmonic oscillator type behavior, with the amplitude of the oscillations decreasing while we are approaching the present day moment, with z = 0. The same damped oscillator type pattern can be observed in the redshift evolution of the scalar field α, depicted in Fig. 7. An oscillatory behavior is also characteristic for the variation with respect to z of the deceleration parameter q, presented in Fig. 8. The behavior is strongly dependent on the numerical values of μ 2 , so that at z = 2 the Universe can be in either a decelerating (q.0), or in an accelerating phase, with q < 0. There is an alternation of the accelerating and decelerating epochs, but, (almost) independently of the values of μ 2 , the Universe enters a very rapid accelerating phase at z ≈ 0.10, which in a very short cosmological time interval decreases the deceleration parameter from q ≈ 1.5 to q ≈ −1. Hence the de Sitter solution is also an attractor of this model.

Scalar function coupling to the metric
As a second example of modified gravity model induced by a quantum perturbation tensor of the form (19) we assume that alpha is a "pure" function of the coordinates, and there is no specific physical field associated to it. Then, by using the least action principle, from Eq. (4) we obtain the gravitational field equations as where we have denoted θ μν = g αβ δT αβ /δg μν , and respectively. As usual, T = T μν g μν is the trace of the energymomentum tensor, and we have denoted θ = θ μν g μν . After contraction Eq. (53) reads Combining Eqs. (53) and (55) we obtain the field equations in the form Assuming that the Lagrangian density of the matter L m depends only on the metric tensor but not on its derivatives, we obtain and respectively. Different kinds of matter Lagrangians may lead to different θ μν , and therefore to different gravitational theories. By taking into account the mathematical identity ∇ ν G μν = 0, after taking the covariant divergence of Eq. (53) we obtain Now it is easy to check that the divergence of the matter energy-momentum tensor is In the following discussion, we consider a perfect fluid described by ρ, the matter energy density, and by p, the thermodynamic pressure. In the comoving frame with u μ = (1, 0, 0, 0) the components of the energy-momentum tensor have the components T μ ν = diag (ρ, − p, − p, − p). Moreover, we will consider that matter obeys an equation of state of the form p = ωρ, where ω = constant, and 0 ≤ ω ≤ 1. Then the generalized Friedmann equations that follow from the gravitational field equations Eq. (53) are or, equivalently, Equation (62) can be alternatively written as Hencė For the deceleration parameter we obtain The presence of an accelerated expansion requiresä >, which imposes on the function α(t) the condition

Conservative models -the ∇ ν T μν = 0 case
In standard general relativity theory, the energy-momentum tensor is conserved (∇ ν T μν = 0), while in modified theories of gravity the classical-defined energy-momentum tensor is not always conserved. However, in the modified gravity theory induced by the quantum fluctuations of the metric, with expectation value of the fluctuations proportional to the metric tensor, thanks to the addition of the function α, we can maintain the conservation of the energy and momentum by constraining the new function. In the following we assume for the matter Lagrangian the form L m = p, and, since the second derivatives of L m with respect to the metric tensor are zero, we obtain for the tensor θ μν the expression Demanding that ∇ ν T μν = 0 we obtain Multiplying by u μ both sides of Eq. (69) we obtaiṅ givinġ For the linear barotropic equation of state with p = ωρ, for ω = 1 we obtaiṅ which gives the density as a function of α as where ρ 0 is an arbitrary constant of integration. On the other hand the conservation of the energy-momentum tensor gives the equation which give for the matter density the standard relation where ρ 0 is an arbitrary constant of integration. From Eqs. (73) and (75) we obtain the scale factor dependence of the function α as ω) . The the first Friedmann equation from which we obtain In the limit of small values of the scale factor α 0 (ω − 1)a 3n << 2, we obtain q ≈ (1 + 3ω)/2 > 0, indicating a decelerating expansion during the early stages of the cosmological evolution. For α 0 (ω − 1)a 3n >> 2, and for very large values of a, q ≈ [3(ω − n) + 1] /2 = 4(ω+2)/(1−ω). Since generally for any realistic cosmological matter equation of state ω < 1, it follows that in both small and large time limits the time evolution of the Universe is decelerating.

Non-conservative cosmological models with
By taking into account that T = ρ − 3 p, from Eq. (60) we immediately obtain or, equivalently, After multiplying Eq. (81) by the four-velocity vector u μ , defined in the comoving reference frame, we obtain By taking into account the mathematical identity we immediately finḋ For a general linear barotropic equation of state of the for p = ω(t)ρ, Eq. (84) becomes For ω = constant we have After taking the derivative of Eq. (61) with respect to the time, and after substitutingḢ with the use of Eq. (65), we obtain for the time deriative of the density the equatioṅ Substituting Eqa. (87) into (86) gives the equation Due to the first generalized Friedmann equation (61), Eq. (88) is identically satisfied during the cosmological evolution. However, a second solution of the field equation can be obtained by also imposing the condition that the first term in Eq. (88) also vanishes identically.
a. The case α = 1 − 1−ω 6(1+ω) 1 ln(a/a 0 ) By requiring that the first term in the left-hand side of Eq. (88) also vanishes, we obtain for the scalar function α(t) the differential equatioṅ from which we obtain where a 0 is an arbitrary constant of integration. By substituting this expression of α into Eq. (63), we obtain the following second order differential equation describing the time evolution of the scale factor: By integration we first obtain where ζ is an arbitrary integration constant. Hence for the time variation of the scale factor we obtain where t 0 is an arbitrary constant of integration. For the time variation of the Hubble function we obtain while the deceleration parameter q of this model is given by In the limit of large times t → ∞, we have q → −1, and therefore the Universe ends in an (approximately de Sitter) accelerating phase. For the time variation of the function α we find In the limit of large times α(t) tends to a constant, lim t→∞ α(t) = 1.

Modified gravity from quantum fluctuations proportional to the matter energy-momentum tensor-K μν = αT μν
Many extensions of standard general theory of relativity are based on the assumption that in certain physical situations spacetime and matter may couple to each other [20,24,31].
Hence it is natural to also consider the case in which the average of the quantum fluctuations of the metric is proportional to the matter energy-momentum tensor, K μν = αT μν , where α is a constant. This approach suggests that the quantum perturbations of the spacetime may also be strongly influenced by the presence of the classical matter.

The gravitational field equations
With the choice K μν = αT μν of the classical form of the average of the quantum fluctuations of the metric tensor we obtain for the first order quantum corrected gravitational Lagrangian the expression By varying the gravitational action given by Eq. (97) with respect to the metric tensor g μν it follows that the classical gravitational field equations corresponding to the gravitational action (97) are given by (for the full details of the derivation see Appendix B) where we have canceled the terms containing the second derivatives of L m with respect to the metric tensor, since in most cases of physical interest they vanish. After contraction of Eq. (98) we find Equation (98) can then be rewritten as

The divergence of the energy-momentum tensor
By taking the divergence of the gravitational field equations (98) we obtain first Hence for the divergence of the matter energy-momentum tensor in the modified gravity model induced by the quantum fluctuations of the metric proportional to the matter energymomentum tensor we find The above results show that generally in this class of models the matter energy-momentum tensor is not conserved. The non-conservation of T μν can be related to particle production processes that takes place due to the quantum fluctuations of the spacetime metric.

Cosmological applications
For simplicity in the following we consider a spatially flat spacetime, in which K = 0. By taking into account the intermediate results we obtain first from the 00 component of Eq. (98) 3ä For μ = ν = 1 we obtain Hence the 11 component of Eq. (98) giveṡ From Eqs. (105) and (107) we obtain Therefore the generalized Friedmann equations for this gravity model take the form where we have denoted For the deceleration parameter we obtain the expression For α = 0 we have = κ 2 (ρc 2 + 3 p)/2, = 2, = −κ 2 (ρc 2 + p)/2, = 0, = −6, = 1, and = 3, respectively. Hence in this limit we recover the standard Friedmann equations of general relativity.

Dust cosmological models with p = 0
By assuming that the matter content of the Universe consists of pressureless dust with p = 0, we obtain immediately λ = 1. Then the field equations can be written as a a = − 3α Hρ + κ 2 ρ 6(αρ + 1) and respectively. Thus we obtain the generalized Friedmann equations of the present model as For the deceleration parameter we obtain By multiplying with a 3 both sides of Eq. (128), taking its time derivative, and considering Eq. (129), we obtain the time evolution of the matter density as

de Sitter type evolution of the dust Universe
By assuming that H = H 0 = constant, the cosmological evolution Eqs. (123) and (124) give for the time evolution of the density the first order differential equatioṅ with the general solution given by where t 0 is an arbitrary constant of integration. When αρ >> ln(ρ), the matter energy density linearly increases in time as ρ(t) ∝ t, indicating that the late time de Sitter expansion of the Universe is triggered by (essentially quantum) particle creation processes.

Cosmological evolution of the dust Universe
From Eq. (122) we obtain for the time evolution of the matter density the differential equatioṅ Substitution of this equation into Eq. (124) gives the evolution equation of the Hubble function as By introducing the set of dimensionless variables (τ, h, α 0 , r ), defined as where H 0 is the present day value of the Hubble function, and after changing the independent time variable from τ to the redshift z, Eqs. (129) and (130) take the form The deceleration parameter is given by .
The system of differential equation must be integrated with the initial conditions h(0) = h 0 and r (0) = r 0 , respectively. It is interesting to note that for h(0) = r (0) = 1, the present day deceleration parameter takes the value q(0) = 1/2. Thus, in order to obtain accelerating models we will adopt an initial condition for the matter energy density so that r (0) = 0.67, while for h we adopt the initial condition h(0) = 1. The variation with respect to the redshift z of the dimensionless Hubble function h, of the dimensionless matter energy density r , and of the deceleration parameter q are represented for z ∈ [0, 1] in Figs. 9, 10, and 11, respectively.
The Hubble function, depicted in Fig. 9, is a monotonically increasing function of z, indication an expansionary evolution of the Universe. The numerical values of h in the chosen redshift range show a very mild dependence on the numerical values of the parameter α 0 . The energy density of the matter r , presented in Fig. 10, is also a monotonically increasing function of the redshift, indicating a time decrease of r during the cosmological evolution. The numerical values of r depend strongly on α 0 . Finally, the deceleration parameter q, shown in Fig. 11, also has a strong dependence on α 0 . In the range z ∈ (0.14, 1], the Universe is in a marginally decelerating phase, with q ≈ 0. In this redshift range the evolution is independent on the numerical values of α 0 . For z < 0.14, the Universe enters in an accelerating phase, and, depending on the values of α 0 , a large range of accelerating models can be constructed. The numerical values of q rapidly increase with increasing α 0 , so that for α 0 = 1, q(0) ≈ −1,

Discussions and final remarks
In the present paper we have considered the cosmological properties of some classes of modified gravity models that are obtained from a first order correction of the quantum metric, as proposed in [40,41]. By assuming that the quantum metric can be decomposed into two components, and by substituting the fluctuating part by its (classical) average value K μν , a large class of modified gravity models can be obtained. As a first step in our study we have derived the general Einstein equations corresponding to an arbitrary K μν . An important property of this class of models is the non-conservation of the matter energy-momentum tensor, which can be related to the physical processes of particle creation due to the quantum effects in the curved spacetime. By assuming that K μν ∝ g μν , a particular class of the modified f (R, T ) gravity theory is obtained. We have investigated in detail the cosmological properties of these models, by assuming that the coupling coefficient between the metric and the average value of the quantum fluctuation tensor is a scalar field with a non-vanishing self-interaction potential, and a simple scalar function. The scalar field self-interaction potential was assumed to be of Higgs type [48], which plays a fundamental role in elementary particle physics as describing the generation of mass in the quantum world.
We have investigated two cosmological models, in which the scalar field is in the minimum of the Higgs potential, and the case of the "complete" Higgs potential. In both approaches in the large time limit the Universe enters an accelerating phase, with the accelerating de Sitter solution acting as an attractor for these cosmological models. However, in the case of the "complete" Higgs potential the redshift evolution of the deceleration parameter q indicates at low redshifts an extremely complex, oscillating behavior. Such a dynamics, as well as the corresponding cosmological evolution may play a significant role in the inflationary/post inflationary reheating phase in the history of the Universe. By assuming that the coupling between the metric and the quantum fluctuations is given by a scalar function, two distinct cases of cosmological models can be obtained. By imposing the conservativity of the energy-momentum tensor for a Universe filled with matter obeying a linear barotropic equation of state, a decelerating cosmological model is obtained. Thus a model could be useful to describe the evolution of the high density Universe at large redshifts. An alternative model with matter creation can also be constructed, by imposing a specific equation of evolution for the time evolution of the coupling function α. This model leads to an approximately de Sitter type expansion, with the deceleration parameter tending to minus one in the large time limit.
A second modified gravity model can be obtained by assuming that the average value of the quantum fluctuation tensor is proportional to the matter energy-momentum tensor, K μν ∝ T μν . This choice leads to an extension of the f R, T, R μν T μν gravity theory [31,32], with the quantum corrected action including an extra term T μν T μν . Hence by considering the effects of the quantum fluctuations of the metric proportional to the matter energy-momentum tensor a particular case of a general f R, T, R μν T μν , T μν T μν modified gravity theory is obtained. The numerical analysis of the cosmological evolution equations for this model show that, depending on the numerical values of the coupling constant α, for a Universe filled with dust matter a large variety of cosmological behaviors can be obtained at low redshifts, with the deceleration parameter varying between a constant (approximately) zero value on a large redshift range, and a de Sitter phase reached at z = 0.
Quantum gravity represents the greatest challenge present day theoretical physics faces. Since no exact solutions for this problem are known, resorting to some approximate methods for studying quantum effects in gravity seems to be the best way to follow. A promising path may be represented by the inclusion of some tensor fluctuating terms in the metric, whose quantum mechanical origin can be well understood. Interestingly enough, such an approach leads to classical gravity models involving geometry-matter coupling, as well as non-conservative matter energy-momentum tensors, and, consequently, to particle creation processes. Hence even the study of the gravitational models with first order quantum corrections can lead to a better understanding of the physical foundations of the modified gravity models with geometry-matter coupling. In the present paper we have investigated some of the cosmological implications of these models, and we have developed some basic tools that could be used to further investigate the quantum effects in gravity, + δT μν T μν = δg μν T α(μ T α ν) − T g μν L m + T μν +2L m T μν + 4T αβ ∂ 2 L m ∂g μν ∂g αβ .
And thus for the total variation of the action we obtain . (B8) The condition δL = 0 gives immediately the field equations (98).