Generalized curvature-matter couplings in modified gravity
GGalaxies , , 410-465; doi:10.3390/galaxies2030410 OPEN ACCESS galaxies
ISSN 2075-4434
Review
Generalized Curvature-Matter Couplings in Modified Gravity
Tiberiu Harko and Francisco S.N. Lobo , * Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK;E-Mail: [email protected] Centro de Astronomia e Astrofísica da Universidade de Lisboa, Campo Grande,Edifício C8 1749-016 Lisboa, Portugal * Author to whom correspondence should be addressed; E-Mail: fl[email protected];Tel.: +351-217-500-986; Fax: +351-217-500-977.
Received: 30 May 2014; in revised form: 7 July 2014 / Accepted: 8 July 2014 /Published: 28 July 2014
Abstract:
In this work, we review a plethora of modified theories of gravity withgeneralized curvature-matter couplings. The explicit nonminimal couplings, for instance,between an arbitrary function of the scalar curvature R and the Lagrangian density ofmatter, induces a non-vanishing covariant derivative of the energy-momentum tensor,implying non-geodesic motion and, consequently, leads to the appearance of an extra force.Applied to the cosmological context, these curvature-matter couplings lead to interestingphenomenology, where one can obtain a unified description of the cosmological epochs.We also consider the possibility that the behavior of the galactic flat rotation curves can beexplained in the framework of the curvature-matter coupling models, where the extra terms inthe gravitational field equations modify the equations of motion of test particles and inducea supplementary gravitational interaction. In addition to this, these models are extremelyuseful for describing dark energy-dark matter interactions and for explaining the late-timecosmic acceleration. Keywords: modified gravity; curvature-matter couplings; cosmology; dark matter
1. Introduction
Recent observations of supernovae, together with the Wilkinson Microwave Anisotropy Probe(WMAP) and Sloan Digital Sky Survey (SDSS) data, lead to the remarkable conclusion that our universe a r X i v : . [ g r- q c ] J u l alaxies , is not just expanding, but has begun to accelerate [1–3]. The resolution of this fundamental questionis extremely important for theoretical cosmology, looking beyond the standard theory of gravity. Thestandard model of cosmology has favored the dark energy models as fundamental candidates responsiblefor the cosmic expansion. However, it is clear that these questions involve not only gravity, butalso particle physics. String theory provides a synthesis of these two branches of physics and iswidely believed to be moving towards a viable quantum gravity theory. One of the key predictionsof string theory is the existence of extra spatial dimensions. In the braneworld scenario, motivatedby recent developments in string theory, the observed three-dimensional universe is embedded in ahigher-dimensional spacetime [4]. The new degrees of freedom belong to the gravitational sector andcan be responsible for the late-time cosmic acceleration.Indeed, detailed theoretical and phenomenological analysis of the relation between the effectivefield yielding dark energy, non-canonical Lagrangians and non-linear gauge kinetic functions havebeen extensively explored in the literature. Note that generalizations of the action functional can beapproached in several ways. For instance, prescriptions consist in replacing the linear scalar curvatureterm in the Einstein–Hilbert action by a function of the scalar curvature, f ( R ) , or by more general scalarinvariants of the theory. This class of theories is often termed higher-order gravity theories [5–7]. Inthis context, infra-red modifications of General Relativity have been extensively explored, includingfour-dimensional modifications to the Einstein–Hilbert action, and the consistency of various candidatemodels has been extensively analyzed. All modified gravity theories induce observational signaturesat the post-Newtonian level, which are translated by the parameterized post-Newtonian (PPN) metriccoefficients arising from these extensions of General Relativity (GR). Thus, generalizations of theEinstein–Hilbert Lagrangian, including quadratic Lagrangians that involve second order curvatureinvariants have also been extensively explored [5]. While the latter modified theories of gravity offeran alternative explanation to the standard cosmological model [8], it offers a paradigm for naturefundamentally distinct from dark energy models of cosmic acceleration, even those that perfectly mimicthe same expansion history. It is a fundamental question to understand how one may observationally andtheoretically differentiate these modified theories of gravity from dark energy models. Thus, one shouldtest these models against large-scale structure and lensing, astrophysical and laboratory measurements,as well as laboratory and space-based equivalence principle experiments. These tests from the SolarSystem, the large-scale structure and lensing essentially restrict the range of allowed modified gravitymodels and, thus, offer a window into understanding the perplexing nature of the cosmic acceleration,and of gravity itself.In this review, we consider an interesting possibility, which includes non-minimal couplings betweenthe scalar curvature and the matter Lagrangian density, introduced in [9], and its extensions andgeneralizations. The specific cases of a curvature coupling to dark energy [10–12], the Maxwell field [13]and an explicit curvature-Yang–Mills coupling [14] was also explored in the context of inflation andof the late-time cosmic acceleration. Indeed, in the context of f ( R ) modified theories of gravity, itwas shown that an explicit coupling between an arbitrary function of the scalar curvature R and theLagrangian density of matter induces a non-vanishing covariant derivative of the energy-momentum,implying non-geodesic motion and, consequently, leads to the appearance of an extra force [9]. Thelatter extra force is orthogonal to the four-velocity, and the corresponding acceleration law was obtained alaxies , in the weak field limit. Connections with Modified Newtonian Dynamics (MOND) and with the Pioneeranomaly were further discussed. These curvature-matter coupling theories include the more evolvedgeneralizations of f ( R, L m ) [15] and f ( R, T ) gravities [16]. Amongst other features, these modelsallow for an explicit breaking of the equivalence principle (EP), which is highly constrained by SolarSystem experimental tests [17,18], by imposing a matter-dependent deviation from geodesic motion.Note that the weak equivalence principle is considered one of the pillars of General Relativity,and in fact, even a sizable part of the modified-gravity community considers this principle as trulyfundamental [19]. This fact has to be stressed, because it demonstrates somehow a limitation for thisclass of theories. However, it has been recently reported, from data of the Abell Cluster A586, that theinteraction of dark matter and dark energy does imply the violation of the equivalence principle [20].Notice that the violation of the equivalence principle is also found as a low-energy feature of somecompactified version of higher-dimensional theories. Indeed, as emphasized by Thibault Damour, itis important to note that the EP is not one of the “universal” principles of physics [21]. It is a heuristichypothesis, which was introduced by Einstein in 1907, and used by him to construct his theory of GeneralRelativity. In modern language, the (Einsteinian) Equivalence Principle (EP) consists in assumingthat the only long-range field with gravitational-strength couplings to matter is a massless spin-2 field.Modern unification theories, and notably String Theory, suggest the existence of new fields (in particular,scalar fields: “dilaton” and “moduli”) with gravitational-strength couplings. In most cases, the couplingsof these new fields violate the EP. If the field is long-ranged, these EP violations lead to many observableconsequences, such as the variation of fundamental “constants”, the non-universality of free fall and therelative drift of atomic clocks, amongst others. The best experimental probe of a possible violation of theEP is to compare the free-fall acceleration of different materials. Further tests of this principle remainimportant and relevant for new physics and do indeed strongly restrict the parameters of the consideredtheory [22,23]. However, it is important to note that, in this context, the EP does not in principle rule outthe specific theory.More specifically, f ( R, L m ) gravity further generalizes f ( R ) gravity by assuming that thegravitational Lagrangian is given by an arbitrary function of the Ricci scalar R and of the matterLagrangian L m [15]. This may also be considered a maximal extension of the Einstein–Hilbertaction. Here, we use the term “maximal extension” in a strict mathematical sense, and wedefine it as “an element in an ordered set that is followed by no other”. For instance,if we consider the set { Grav Theor } of gravity theories as given by { Grav Theor } = { standard general relativity , f ( R ) gravity , linear curvature − matter coupling , ..., f ( R, L m ) } , the f ( R, L m ) gravity theory represent the maximal extension of the set of all gravitational theoriesconstructed in a Riemann space and with the action depending on the Ricci scalar and the matterLagrangian only. In [15], the gravitational field equations of the f ( R, L m ) gravity theory in the metricformalism were obtained, as well as the equations of motion for test particles, which follow from thecovariant divergence of the energy-momentum tensor. The equations of motion for test particles canalso be derived from a variational principle in the particular case in which the Lagrangian density of thematter is an arbitrary function of the energy-density of the matter only. In [15], the Newtonian limit ofthe equation of motion was also considered, and a procedure for obtaining the energy-momentum tensorof the matter was presented. In particular, the gravitational field equations and the equations of motion alaxies , for a particular model in which the action of the gravitational field has an exponential dependence onthe standard general relativistic Hilbert–Einstein Lagrange density were also derived. The f ( R, L m ) gravitational theory was further generalized by considering the novel inclusion of a scalar field and akinetic term constructed from the gradients of the scalar field, respectively [24]. Specific models witha nonminimal coupling between the scalar field and the matter Lagrangian were further explored. Weemphasize that these models are extremely useful for describing an interaction between dark energy anddark matter and for explaining the late-time cosmic acceleration.In the context of f ( R, T ) modified theories of gravity, the gravitational Lagrangian is given byan arbitrary function of the Ricci scalar R and of the trace of the energy-momentum tensor T [16].The gravitational field equations in the metric formalism were obtained, and it was shown that thesefield equations explicitly depend on the nature of the matter source. The field equations of severalparticular models, corresponding to some explicit forms of the function f ( R, T ) , were also presented.Furthermore, in [16], the Newtonian limit of the equation of motion was also analyzed, and constraintson the magnitude of the extra-acceleration were obtained by analyzing the perihelion precession of theplanet Mercury in the framework of the present model. An interesting specific case, namely the f ( R, T φ ) model, was also analyzed in detail, where T φ is the trace of the stress-energy tensor of a self-interactingscalar field. The cosmological implications of the model were briefly considered. We refer the readerto [16] for more details.The non-minimally curvature-matter-coupled f ( R, T ) gravitational theory was further generalizedby considering the inclusion of a contraction of the Ricci tensor with the matter energy-momentumtensor [25–27]. We emphasize that examples of such couplings can be found in the Einstein–Born–Infeldtheories [28] when one expands the square root in the Lagrangian. An interesting feature of this theory isthat in considering a traceless energy-momentum tensor, i.e. , T = 0 , the field equations of f ( R, T ) gravity reduces to those of f ( R ) gravity theories, while the presence of the R µν T µν coupling termstill entails a non-minimal coupling to matter. In [25], the Newtonian limit of the f ( R, T, R µν T µν ) gravitational theory was considered, and an explicit expression for the extra-acceleration, which dependson the matter density, was obtained in the small velocity limit for dust particles. The so-calledDolgov–Kawasaki instability [29] was also analysed in detail, and the stability conditions of the modelwith respect to local perturbations was obtained. A particular class of gravitational field equations canbe obtained by imposing the conservation of the energy-momentum tensor [25]. In this context, thecorresponding field equations for the conservative case was also derived by using a Lagrange multipliermethod, from a gravitational action that explicitly contains an independent parameter multiplying thedivergence of the energy-momentum tensor [25]. The cosmological implications of the model wereextensively investigated for both the conservative and non-conservative cases, and several classesof analytical solutions were obtained. In [26], the Friedmann-Lemaître-Robertson-Walker (FRLW)cosmological dynamics for several versions of the f ( R, T, R µν T µν ) gravity theory was also considered.The reconstruction of the above action was explicitly analyzed, including the numerical reconstructionfor the occurrence of the Λ CDM (Lambda-Cold Dark Matter) universe. De Sitter universe solutionsin the presence of non-constant fluids were also presented, and the problem of matter instability wasfurther discussed. alaxies , All of the above gravitational modifications are based on the curvature description of gravity.However, an interesting and rich class of modified gravity arises by modifying the action of the equivalenttorsional formulation of General Relativity. The latter approach has been denoted the “TeleparallelEquivalent of General Relativity” (TEGR) [30–35] and consists of replacing the torsion-less Levi–Civitaconnection by the curvature-less Weitzenböck one and using the vierbein instead of the metric as thefundamental field. In this formulation, the Lagrangian of the theory is constructed by contractions of thetorsion tensor. Thus, in analogy to f ( R ) gravity, if one wishes to modify gravity in this formulation, thesimplest approach would be in generalizing the torsion scalar T to an arbitrary function f ( T ) [36–38].Now, in the context of TEGR models, one may also construct an extension of f ( T ) gravity with theinclusion of a non-minimal torsion-matter coupling in the action [39]. The resulting theory is a novelgravitational modification, since it is different from both f ( T ) gravity, as well as from the non-minimalcurvature-matter-coupled theory. The cosmological application of this new theory proves to be veryinteresting. In particular, an effective dark energy sector was obtained, where the equation-of-stateparameter can be quintessence or phantom-like, or cross the phantom-divide, while for a large range ofthe model parameters, the Universe results in a de Sitter, dark-energy-dominated, accelerating phase.Additionally, early-time inflationary solutions were also obtained, and thus, one can provide a unifieddescription of the cosmological history. We refer the reader to [39] for more details.In addition to the latter, non-minimal torsion-matter coupling, an alternative extension of f ( T ) gravitywas explored by considering a general coupling of the torsion scalar T with the trace of the matterenergy-momentum tensor T [40]. The resulting f ( T , T ) theory is a new modified gravity, since it isdifferent from all of the existing torsion- or curvature-based constructions. Applied to a cosmologicalframework, it also leads to interesting phenomenology. In particular, one can obtain a unified descriptionof the initial inflationary phase, the subsequent non-accelerating, matter-dominated expansion and,then, the transition to late-time accelerating phase. In the far future, the Universe results either ina de Sitter exponential expansion or in eternal power-law accelerated expansions. A similar analysiswas investigated in [41], where using a perturbational approach, the stability of the solutions and, inparticular, of the de Sitter phase was explored. Furthermore, the constraints imposed by the energyconditions were also considered. We refer the reader to [40,41] for more details.The possibility that the behavior of the rotational velocities of test particles gravitating around galaxiescan be explained in the framework of modified gravity models with nonminimal curvature-mattercoupling was considered in [42]. Generally, the dynamics of test particles around galaxies, as wellas the corresponding mass deficit, is explained by postulating the existence of dark matter. The extraterms in the gravitational field equations with curvature-matter coupling modify the equations of motionof test particles and induce a supplementary gravitational interaction. Starting from the variationalprinciple describing the particle motion in the presence of the non-minimal coupling, the expression ofthe tangential velocity of a test particle, moving in the vacuum on a stable circular orbit in a sphericallysymmetric geometry, was derived. The tangential velocity depends on the metric tensor components, aswell as on the coupling function between matter and geometry. The Doppler velocity shifts were alsoobtained in terms of the coupling function. If the tangential velocity profile is known, the coupling termbetween matter and curvature can be obtained explicitly in an analytical form. The functional form of thisfunction was obtained in two cases, for a constant tangential velocity and for an empirical velocity profile alaxies , obtained from astronomical observations, respectively. Therefore, these results open the possibility ofdirectly testing the modified gravity models with a non-minimal curvature-matter coupling by usingdirect astronomical and astrophysical observations at the galactic or extra-galactic scale. These issueswill be reviewed in detail below.It is rather important to discuss the theoretical motivations for these curvature-matter couplingtheories more carefully. For instance, one may wonder if there is a fundamental theory or model thatreduces to one of these theories in some limit. Indeed, for the specific case of the linear nonminimalcurvature-matter coupling, one can show that this theory can be expressed as a scalar-tensor theorywith two scalar fields [43]. These are reminiscent of the extensions of scalar-tensor gravity, whichinclude similar couplings, such as the theories considered by Damour and Polyakov [44], where it wasshown that string-loop modifications of the low-energy matter couplings of the dilaton may providea mechanism for fixing the vacuum expectation value of a massless dilaton. In addition to this, theresults presented in [44] provide a new motivation for trying to improve by several orders of magnitudethe various experimental tests of Einstein’s equivalence principle, such as the universality of free falland the constancy of the fundamental constants, amongst others. In addition to this, note that themodern revival of the Kaluza–Klein theory often leads to the introduction of several scalar fields, i.e. ,the “compacton” [45]. Superstring theory also leads to several scalar fields coupled to the macroscopicdistribution of energy; indeed, the “dilaton” is already present in ten dimensions, and the “compacton”comes from a dimensional reduction [45]. Quantum-motivated, higher order generalizations of Einstein’sGeneral Relativity, under some conditions, are also equivalent to adding several scalar fields to theEinstein–Hilbert action [46]. In addition to this, we note that the nonminimal curvature-matter couplingalso arises from one-loop vacuum-polarization effects in the formulation of quantum electrodynamicsin curved spacetimes [47]. All of the above considerations make it natural to consider modified fieldsof gravity containing several scalar fields, which, in some cases, are equivalent, in the present context,to the presence of a nonminimal curvature-matter couplings. These considerations further motivate theanalysis of the coupling between matter and curvature.In this work, based on published literature, we review these generalized curvature-mattercoupling-modified theories of gravity. The paper is outlined in the following manner: In Section 2,we briefly introduce the linear curvature-matter coupling and present some of its interesting features. InSection 3, we generalize the latter linear curvature-matter coupling by considering the maximal extensionof the Einstein–Hilbert action, which results in the f ( R, L m ) gravitational theory, and in Section 4, weextend the theory with the inclusion of general scalar field and kinetic term dependencies. In Section 5,we consider another extension of general relativity, namely f ( R, T ) modified theories of gravity, wherethe gravitational Lagrangian is given by an arbitrary function of the Ricci scalar R and of the trace ofthe energy-momentum tensor T . In Section 6, we include an explicit invariant Ricci-energy-momentumtensor coupling and explore some of its astrophysical and cosmological phenomenology. This lattertheory is motivated considering a traceless energy-momentum tensor, T = 0 ; the gravitational fieldequations for the f ( R, T ) theory reduce to that of the f ( R ) gravity, and all non-minimal couplings ofgravity to the matter field vanish, while the inclusion of the R µν T µν term still allows a nonminimalcoupling. The possibility of explaining dark matter as a consequence of the curvature-matter coupling alaxies , is considered in Section 7. Finally, in Section 8, we briefly discuss our results and conclude ourreview paper.
2. Gravity Theories with Linear Curvature-Matter Coupling
In the present Section, we review some of the basic features and properties of the simplest classof models involving a nonminimal coupling between matter and geometry, which is linear in thematter Lagrangian. For the sake of comparison of the different approaches for the description of thegravitational interaction, we also briefly introduce the f ( R ) modified gravity theory. f ( R ) Gravity
One of the simplest extensions of standard general relativity, based on the Hilbert–Einstein action, isthe f ( R ) -modified theory of gravity, whose action takes the form [48,49]: S = (cid:90) (cid:20) πG f ( R ) + L m (cid:21) √− g d x (1)where f ( R ) is an arbitrary analytical function of the Ricci scalar R and L m is the Lagrangian densitycorresponding to matter.Varying the action with respect to the metric g µν yields the field equations of f ( R ) gravity: F R µν − f g µν − ( ∇ µ ∇ ν − g µν (cid:3) ) F = 8 πG T µν (2)where we have denoted F = df /dR , and T µν is the standard, minimally coupled, matterenergy-momentum tensor. Note that the covariant derivative of the field equations and of the matterenergy-momentum tensor vanishes for all f ( R ) by means of generalized Bianchi identities [9,50].By contracting the field equations, Equation (2), we obtain the useful relation: (cid:3) F + F R − f = 8 πG T (3)where T is the trace of the energy-momentum tensor and from which one verifies that the Ricci scalar isnow a fully dynamical degree of freedom.By introducing the Legendre transformation { R, f } → { φ, V } defined as: φ ≡ F ( R ) , V ( φ ) ≡ R ( φ ) F − f ( R ( φ )) (4)the field equations of f ( R ) gravity can be reformulated as [51–54]: R µν − g µν R = 8 π Gφ T µν + θ µν (5)where: θ µν = − V ( φ ) g µν + 1 φ ( ∇ µ ∇ ν − g µν (cid:3) ) φ (6)Using these variables, Equation (3) takes the following form: (cid:3) φ + 2 V ( φ ) − φ dVdφ = 8 πG T (7) alaxies , In this scalar-tensor representation, the field equations of f ( R ) gravity can be derived from aBrans–Dicke-type gravitational action, with parameter ω = 0 , given by: S = 116 πG (cid:90) [ φR − V ( φ ) + L m ] √− g d x (8)The only requirement for the f ( R ) model equations to be expressed in the form of a Brans–Dicketheory is that F ( R ) be invertible, that is, R ( F ) exists [55]. This condition is necessary for theconstruction of V ( φ ) .The modification of the standard Einstein–Hilbert action leads to the appearance in the field equationsof an effective gravitational constant G eff = G/φ , which is a function of the curvature. Secondly, anew source term for the gravitational field, given by the tensor θ µν , is also induced. The tensor θ µν is determined by the trace of the energy-momentum tensor via Equation (7), which thus acts as anindependent physical parameter determining the metric of the space-time. The action of the f ( R ) modified gravity theory can be generalized by introducing in the action thelinear nonminimal coupling between matter and geometry. The action of this modified gravity theory isgiven by [9]: S = (cid:90) (cid:26) f ( R ) + [1 + λf ( R )] L m (cid:27) √− g d x (9)where the factors f i ( R ) (with i = 1 , ) are arbitrary functions of the Ricci scalar R and L m is the matterLagrangian density, and the coupling constant λ determines the strength of the interaction between f ( R ) and the matter Lagrangian.Now, varying the action with respect to the metric g µν provides the following field equations: F ( R ) R µν − f ( R ) g µν − ∇ µ ∇ ν F ( R ) + g µν (cid:3) F ( R ) = − λF ( R ) L m R µν +2 λ ( ∇ µ ∇ ν − g µν (cid:3) ) L m F ( R ) + [1 + λf ( R )] T ( m ) µν (10)where F i ( R ) = f (cid:48) i ( R ) , and the prime represents the derivative with respect to the scalar curvature. Thematter energy-momentum tensor is defined as: T ( m ) µν = − √− g δ ( √− g L m ) δ ( g µν ) (11)Throughout this paper, we use the metric formalism. However, the field equations and theequations of motion were derived, for massive test particles in modified theories of gravity for a linearcurvature-matter coupling, using the Palatini formalism in [56].A general property of these nonminimal curvature-matter coupling theories is the non-conservationof the energy-momentum tensor. This can be easily verified by taking into account the covariantderivative of the field Equation (10), the Bianchi identities, ∇ µ G µν = 0 , and the following identity ( (cid:3) ∇ ν − ∇ ν (cid:3) ) F i = R µν ∇ µ F i , which then imply the following relationship: ∇ µ T ( m ) µν = λF λf (cid:2) g µν L m − T ( m ) µν (cid:3) ∇ µ R (12) alaxies , Thus, the coupling between the matter and the higher derivative curvature terms may be interpretedas an exchange of energy and momentum between both. Analogous couplings arise after a conformaltransformation in the context of scalar-tensor theories of gravity and also in string theory. In the absenceof the coupling, one verifies the conservation of the energy-momentum tensor [50], which can also beverified from the diffeomorphism invariance of the matter part of the action. It is also interesting tonote that, from Equation (12), the conservation of the energy-momentum tensor is verified if f ( R ) is aconstant or the matter Lagrangian is not an explicit function of the metric.In order to test the motion in our model, we consider for the energy-momentum tensor of matter aperfect fluid T ( m ) µν = ( ρ + p ) U µ U ν + pg µν , where ρ is the overall energy density and p , the pressure,respectively. The four-velocity, U µ , satisfies the conditions U µ U µ = − and ∇ ν U µ U µ = 0 . We alsointroduce the projection operator h µλ = g µλ + U µ U λ from which one obtains h µλ U µ = 0 .By contracting Equation (12) with the projection operator h µλ , one deduces the following expression: ( ρ + p ) g µλ U ν ∇ ν U µ − ( ∇ ν p )( δ νλ − U ν U λ ) − λF λf ( L m − p ) ( ∇ ν R ) ( δ νλ − U ν U λ ) = 0 (13)Finally, contraction with g αλ gives rise to the equation of motion for a fluid element: DU α ds ≡ dU α ds + Γ αµν U µ U ν = f α (14)where the extra force is given by: f α = 1 ρ + p (cid:34) λF λf ( L m − p ) ∇ ν R + ∇ ν p (cid:35) h αν (15)As one can immediately verify, the extra force f α is orthogonal to the four-velocity of the particle, f α U α = 0 which can be seen directly, from the properties of the projection operator. This is consistentwith the usual interpretation of the force, according to which only the component of the four-force thatis orthogonal to the particle’s four-velocity can influence its trajectory.It has also been shown that an f ( G ) -modified Gauss–Bonnet gravity with a non-minimal coupling tomatter also induces an extra force, which is normal to their four-velocities and, as a result, moves alongnongeodesic world-lines [57].A particularly intriguing feature is that the extra force depends on the form of the Lagrangiandensity. Note that considering the Lagrangian density L m = p , where p is the pressure, the extra-forcevanishes [58]. It has been argued that this is not the unique choice for the matter Lagrangian densityand that more natural forms for L m , such as L m = − ρ , do not imply the vanishing of the extra-force.Indeed, in the presence of nonminimal coupling, they give rise to two distinct theories with differentpredictions [59], and this issue has been further investigated in different contexts [60,61], includingphantom energy [62]. In this context, a matter Lagrangian density as an arbitrary function of theenergy-density of the matter only was explored [63] (this possibility has also been explored in fivedimensions [64] and in the braneworld context [65]). It was also argued that the correspondingenergy-momentum tensor of the matter in modified gravity models with non-minimal coupling is moregeneral than the usual general-relativistic energy-momentum tensor for perfect fluids [66], and it containsa supplementary, equation of state-dependent term, which could be related to the elastic stresses in the alaxies , body or to other forms of internal energy. Therefore, the extra force induced by the coupling betweenmatter and geometry never vanishes as a consequence of the thermodynamic properties of the systemor for a specific choice of the matter Lagrangian, and it is non-zero in the case of a fluid of dustparticles. In the following subsection, we discuss in detail the problem of the matter Lagrangian andof the energy-momentum tensor in modified gravity theories with linear curvature-matter coupling [66]. It is an interesting and novel feature of the f ( R ) gravity with non-minimal curvature-mattercoupling that the matter Lagrangian, and the energy-momentum tensor obtained from it, are notmodel-independent quantities, but they are completely and uniquely determined by the coupling functionbetween matter and geometry. This important result can be obtained by deriving first the equations ofmotion of test particles (or test fluid) in the modified gravity model from a variational principle andthen considering the Newtonian limit of the particle action for a fluid obeying a barotropic equation ofstate [66]. The energy-momentum tensor of the matter obtained in this way is more general than theusual general-relativistic energy-momentum tensor for perfect fluids, and it contains a supplementaryterm that may be related to the elastic stresses in the test fluid or to other sources of internal energy.Since we assume that the matter obeys a barotropic equation of state, the matter Lagrangian can beexpressed either in terms of the density or in terms of the pressure, and in both representations, thephysical description of the system is equivalent. Therefore, the presence (or absence) of the extra-force isindependent of the specific form of the matter Lagrangian, and it never vanishes, except in the particularcase of (un)physical matter systems with zero sound speed. In particular, in the case of dust particles,the extra-force is always non-zero.We define the energy-momentum tensor of the matter as [67]: T µν = − √− g (cid:20) ∂ ( √− gL m ) ∂g µν − ∂∂x λ ∂ ( √− gL m ) ∂ ( ∂g µν /∂x λ ) (cid:21) (16)By assuming that the Lagrangian density L m of the matter depends only on the metric tensorcomponents g µν , and not on its derivatives, we obtain: T µν = L m g µν − ∂L m ∂g µν (17)By taking into account the explicit form of the field equations Equation (10), one obtains for thecovariant divergence of the energy-momentum tensor the equation: ∇ µ T µν = 2 {∇ µ ln [1 + λf ( R )] } ∂L m ∂g µν (18)As a specific example of generalized gravity models with linear curvature-matter coupling,we consider the case in which matter, assumed to be a perfect thermodynamic fluid, obeys a barotropicequation of state, with the thermodynamic pressure p being a function of the rest mass density of thematter (for short: matter density) ρ only, so that p = p ( ρ ) . In this case, the matter Lagrangian density,which, in the general case, could be a function of both density and pressure, L m = L m ( ρ, p ) , or of only alaxies , one of the thermodynamic parameters, becomes an arbitrary function of the density of the matter ρ only,so that L m = L m ( ρ ) . Then, the matter energy-momentum tensor is obtained as [66]: T µν = ρ dL m dρ U µ U ν + (cid:18) L m − ρ dL m dρ (cid:19) g µν (19)where the four-velocity U µ = dx µ /ds satisfies the condition g µν U µ U ν = − . To obtain Equation (19),we have imposed the condition of the conservation of the matter current: ∇ ν ( ρU ν ) = 0 (20)and we have used the relation: δρ = 12 ρ ( g µν − U µ U ν ) δg µν (21)whose proof is given in the Appendix of [66]. With the use of the identity U ν ∇ ν U µ = d x µ /ds +Γ µνα U ν U α , from Equations (12) and (19), we obtain the equation of motion of a massive test particle,or of a test fluid in the modified gravity model with linear coupling between matter and geometry, as: dU µ ds + Γ µνα U ν U α = f µ (22)where: f µ = −∇ ν ln (cid:26) [1 + λf ( R )] dL m ( ρ ) dρ (cid:27) ( U µ U ν − g µν ) (23)The extra-force f µ , generated due to the presence of the coupling between matter and geometry, isperpendicular to the four-velocity, f µ U µ = 0 .The equation of motion Equation (22) can be obtained from the variational principle [66]: δS p = δ (cid:90) L p ds = δ (cid:90) (cid:112) Q (cid:112) g µν U µ U ν ds = 0 (24)where S p and L p = √ Q (cid:112) g µν U µ U ν are the action and the Lagrangian density for test particles (testfluid), respectively, and: (cid:112) Q = [1 + λf ( R )] dL m ( ρ ) dρ (25)To prove this result, we start with the Lagrange equations corresponding to the action Equation (24): dds (cid:18) ∂L p ∂U α (cid:19) − ∂L p ∂x α = 0 (26)Since: ∂L p ∂U α = (cid:112) Qu α (27)and: ∂L p ∂x α = 12 (cid:112) Qg µν,α U µ U ν + 12 Q ,α Q (28)where a comma indicates the derivative with respect to x λ , a straightforward calculation gives theequations of motion of the particle as: d x µ ds + Γ µνα U ν U α + ( U µ U ν − g µν ) ∇ ν ln (cid:112) Q = 0 (29)By simple identification with the equation of motion of the modified gravity model with linearcurvature-matter coupling, given by Equation (22), we obtain the explicit form of √ Q , as givenby Equation (25). When √ Q → , we reobtain the standard general relativistic equation forgeodesic motion. alaxies , ds ≈ (cid:112) φ − (cid:126)v dt ≈ (cid:18) φ − (cid:126)v (cid:19) dt (30)where φ is the Newtonian potential and (cid:126)v is the usual tridimensional velocity of the fluid. By representingthe function √ Q as: (cid:112) Q = dL m ( ρ ) dρ + λf ( R ) dL m ( ρ ) dρ (31)in the first order of approximation, the equations of motion of the fluid can be obtained from thevariational principle: δ (cid:90) (cid:20) dL m ( ρ ) dρ + λf ( R ) dL m ( ρ ) dρ + φ − (cid:126)v (cid:21) dt = 0 (32)and are given by: (cid:126)a = −∇ φ − ∇ dL m ( ρ ) dρ − ∇ U E = (cid:126)a N + (cid:126)a H + (cid:126)a E (33)where (cid:126)a is the total acceleration of the system: (cid:126)a N = −∇ φ (34)is the Newtonian gravitational acceleration, and: (cid:126)a E = −∇ U E = − λ ∇ (cid:20) f ( R ) dL m ( ρ ) dρ (cid:21) (35)is a supplementary acceleration induced due to the curvature-matter coupling. As for the term: (cid:126)a H = −∇ (cid:20) dL m ( ρ ) dρ (cid:21) (36)it has to be identified with the hydrodynamic acceleration term in the perfect fluid Euler equation: (cid:126)a H = −∇ dL m ( ρ ) dρ = −∇ (cid:90) ρρ dpdρ dρρ (37)where ρ , an integration constant, plays the role of a limiting density.2.3.2. The Matter Lagrangian in Modified Gravity Theories with Curvature-Matter CouplingWith the use of Equation (37), the matter Lagrangian in modified gravity theories withcurvature-matter coupling can be obtained by a simple integration as [66]: L m ( ρ ) = ρ [1 + Π ( ρ )] − (cid:90) pp dp (38)where: Π ( ρ ) = (cid:90) pp dpρ (39) alaxies , and we have normalized an arbitrary integration constant to one; p is an integration constant, or alimiting pressure. The corresponding energy-momentum tensor of matter is given by: T µν = { ρ [1 + Φ ( ρ )] + p ( ρ ) } U µ U ν − p ( ρ ) g µν (40)respectively, where: Φ ( ρ ) = (cid:90) ρρ pρ dρ = Π ( ρ ) − p ( ρ ) ρ (41)and with all of the constant terms included in the definition of p . By introducing the energy density ofthe body according to the definition: ε = ρ [1 + Φ ( ρ )] (42)the energy-momentum tensor of a test fluid can be written in the modified gravity models withcurvature-matter coupling in a form similar to the standard general relativistic case: T µν = [ ε ( ρ ) + p ( ρ )] U µ U ν + p ( ρ ) g µν (43)From a physical point of view, Φ ( ρ ) can be interpreted as the elastic (deformation) potential energyof the body, and therefore, Equation (40) corresponds to the energy-momentum tensor of a compressibleelastic isotropic system [68]. The matter Lagrangian can also be written in the simpler form: L m ( ρ ) = ρ Φ ( ρ ) (44)If the pressure does not have a thermodynamic or radiative component, one can take p = 0 . If thepressure is a constant background quantity, independent of the density, so that p = p , then L m ( ρ ) = ρ ,and the matter energy-momentum tensor takes the form corresponding to dust: T µν = ρU µ U ν (45)Since matter is supposed to obey a barotropic equation of state, these results are independent ofthe concrete representation of the matter Lagrangian in terms of the thermodynamic quantities [66].The same results are obtained by assuming L m = L m ( p ) ; due to the equation of state, ρ and p are freelyinterchangeable thermodynamic quantities, and the Lagrangians expressed in terms of ρ and p only arecompletely equivalent. More general situations, in which the density and pressure are functions of theparticle number and temperature, respectively, and the equation of state is given in a parametric form,can be analyzed in a similar way.The forms of the matter Lagrangian and the energy-momentum tensor are strongly dependent on theequation of state of the test fluid. For example, if the barotropic equation of state is linear, p = ( γ − ρ , γ = constant, ≤ γ ≤ , then: L m ( ρ ) = ρ (cid:26) γ − (cid:20) ln (cid:18) ρρ (cid:19) − (cid:21)(cid:27) (46)and Φ ( ρ ) = ( γ −
1) ln ( ρ/ρ ) , respectively. For the case of a polytropic equation of state p = Kρ /n , K, n = constant, we have: L m ( ρ ) = ρ + K (cid:18) n n + 1 − (cid:19) ρ /n (47) alaxies , and Φ ( ρ ) = Knρ /n = np ( ρ ) , respectively, where we have taken for simplicity ρ = p = 0 . For atest fluid satisfying the ideal gas equation of state: p = k B µ ρT (48)where k B is Boltzmann’s constant, T is the temperature and µ is the mean molecular weight, we obtain: L m ( ρ ) = ρ (cid:26) k B Tµ (cid:20) ln (cid:18) ρρ (cid:19) − (cid:21)(cid:27) + p (49)In the case of a physical system satisfying the ideal gas equation of state, the extra accelerationinduced by the presence of the non-minimal coupling between matter and geometry is given by: (cid:126)a E ≈ − λ k B Tµ ∇ (cid:20) f ( R ) ln ρρ (cid:21) (50)and it is proportional to the temperature of the fluid. It is also interesting to note that the limiting densityand pressure ρ and p generate in the energy-momentum tensor some extra constant terms, which maybe interpreted as dark energy.In conclusion, the extra-force induced by the coupling between matter and geometry does not vanishfor any specific choices of the matter Lagrangian. In the case of the dust, with p = 0 , the extra force isgiven by: f µ = −∇ ν ln [1 + λf ( R )] ( U µ U ν − g µν ) (51)and it is independent of the thermodynamic properties of the system, being completely determined bythe geometry, kinematics and coupling. In the limit of small velocities and weak gravitational fields, theextra-acceleration of a dust fluid is given by: (cid:126)a E = − λ ∇ [ f ( R )] (52)The thermodynamic condition for the vanishing of the extra-force is: ∂L m ∂g µν = 12 (cid:18) ∂L m ∂ρ (cid:19) ρ ( g µν − U µ U ν ) = 0 (53)only. If the matter Lagrangian is written as a function of the pressure, then: ∂L m ∂ρ = (cid:18) ∂L m ∂p (cid:19) (cid:18) ∂p∂ρ (cid:19) (54)and for all physical systems satisfying an equation of state (or, equivalently, for all systems with anon-zero sound velocity), the extra-force is non-zero. Therefore, the curvature-matter coupling isintroduced in the generalized gravity models with a curvature-matter coupling in a consistent way. Thecoupling determines all of the physical properties of the system, including the extra-force, the matterLagrangian and the energy-momentum tensor, respectively.Hence, we have shown that in f ( R ) gravity with a non-minimal coupling, the matter Lagrangian andthe corresponding energy-momentum tensor are not model and thermodynamic parameters (independentquantities), but they are completely and uniquely determined by the nature of the coupling between alaxies , matter and geometry, which, in the present model, is given by the function f ( R ) . We have obtainedthis result by deriving first the equations of motion in the modified gravity model with curvature-mattercoupling from a variational principle and, then, by taking the Newtonian limit of the particle action for afluid obeying a barotropic equation of state. The energy-momentum tensor of the matter obtained in thisway is more general than the usual general-relativistic energy-momentum tensor for perfect fluids, and itcontains a supplementary term that may be related to the elastic stresses in the body, or to other sourcesof internal energy. The matter Lagrangian can be expressed either in terms of the density or in terms ofthe pressure, and in both representations, the physical description of the system is equivalent. Therefore,the presence (or absence) of the extra-force is independent of the specific form of the matter Lagrangian,and it never vanishes, except in the case of (un)physical systems with zero sound speed. In particular, inthe case of dust particles, the extra-force is always non-zero. As we have shown in Section 2.1, pure f ( R ) gravity, with action given by Equation (1), is equivalentto a scalar-tensor theory, with action given by Equation (8). The equivalence between the modifiedgravity models with linear geometry-matter coupling was established by Faraoni [69]. In the following,we show that the action Equation (9) is equivalent with a scalar-tensor Brans–Dicke-type theory, with asingle scalar field, a vanishing Brans–Dicke parameter ω and an unusual coupling of the potential U ( ψ ) of the theory to matter.By introducing a new field φ , the action Equation (9) becomes: S (cid:90) d x √− g (cid:26) f ( φ )2 + 12 df dφ ( R − φ ) + [1 + λf ( φ )] L m (cid:27) (55)Next, we further introduce the field ψ ( φ ) ≡ f (cid:48) ( φ ) (with a prime denoting the differentiation withrespect to φ ), and we obtain for the action the expression: S = (cid:90) d x √− g (cid:20) ψR − V ( ψ ) + U ( ψ ) L m (cid:21) (56)where: V ( ψ ) = φ ( ψ ) f (cid:48) [ φ ( ψ )] − f [ φ ( ψ )]2 (57) U ( ψ ) = 1 + λf [ φ ( ψ )] (58) φ ( ψ ) must be obtained by inverting ψ ( φ ) ≡ f (cid:48) ( φ ) . The actions Equations (9) and (56) are equivalentwhen f (cid:48)(cid:48) ( R ) (cid:54) = 0 [69]. We can see this by setting φ = R , and then Equation (56) reduces trivially toEquation (9). On the other hand, the variation of Equation (55) with respect to φ gives: ( R − φ ) f (cid:48)(cid:48) ( φ ) + 2 λf (cid:48) ( φ ) L m = 0 (59)In a vacuum, we have L m = 0 , and Equation (59) gives φ = R whenever f (cid:48)(cid:48) (cid:54) = 0 [51–54].However, in the presence of matter, there seem to be other possibilities, which, however, can be excludedas follows. alaxies , When L m (cid:54) = 0 , the action Equations (9) and (55) are equivalent if: ( R − φ ) f (cid:48)(cid:48) ( φ ) + 2 λf (cid:48)(cid:48) ( φ ) L m (cid:54) = 0 (60)When Equation (59) is satisfied, we have a pathological case. It corresponds to: λf ( φ ) L m = f (cid:48) ( φ )2 ( φ − R ) − f ( φ )2 (61)However, if Equation (61) holds, then the action Equation (55) reduces to the trivial case of purematter without the gravity sector. Then, it follows that in the modified gravity theories with lineargeometry-matter coupling, the actions Equations (9) and (56) are equivalent when f (cid:48)(cid:48) ( R ) (cid:54) = 0 , similar tothe case of pure f ( R ) gravity [51–54]. This linear nonminimal curvature-matter coupling has been extensively explored in a plethora ofcontexts in the literature. For instance, the equations of motion of test bodies for a theory withnonminimal coupling by means of a multipole method [70] was also studied, and it was shown thatthe propagation equations for pole-dipole particles allow for a systematic comparison with the equationsof motion of general relativity and other gravity theories.The consequences that a non-minimal coupling between curvature and matter can have on thedynamics of perfect fluids has also been investigated [71]. It was argued that the presence of a static,axially-symmetric, pressureless fluid does not imply an asymptotically Minkowski space-time, such asin General Relativity. This feature can be attributed to a pressure mimicking mechanism related to thenon-minimal coupling. The case of a spherically symmetric black hole surrounded by fluid matter wasanalyzed, and it was shown that under equilibrium conditions, the total fluid mass is about twice that ofthe black hole [71].The consequences of the curvature-matter coupling on stellar equilibrium and constraints on thenonminimal coupling was considered, where particular attention was paid to the validity of theNewtonian regime and on the boundary and exterior matching conditions [72]. This explicit “anomalous”coupling of the Ricci curvature to matter has also raised the question of curvature instabilities, andin [69,73–76], constraints imposed by the energy condition and the conditions in order to avoidthe notorious Dolgov–Kawasaki instability [29] were obtained. It has also been claimed that thecurvature-matter action leads to a theory of gravity that includes higher order derivatives of the matterfields without introducing more dynamics in the gravity sector and, therefore, cannot be a viable theoryfor gravitation [77,78]. However, we emphasize that the results of [77,78] only apply to the specific caseof f ( R ) = R and not to action Equation (9) in general. For more generic functional forms, the theorypropagates extra degrees of freedom and the conclusions of [77,78] do not apply. The relation betweenthese theories and ordinary scalar-tensor gravity was also analyzed, as well as its implications for theequivalence principle [43].In fact, the theoretical consistency of these nonminimal curvature-matter couplings was studiedusing a scalar field Lagrangian to model the matter content [79]. The conditions that the couplingdoes not introduce ghosts, classical instabilities or superluminal propagation of perturbations were alaxies , derived. These consistency conditions were then employed to rule out or severely restrict the formsof the non-minimal coupling functions [79]. Lagrange–Noether methods have been used to derivethe conservation laws for models in which matter interacts non-minimally with the gravitationalfield [80]. Furthermore, a covariant derivation of the equations of motion for test bodies for a wideclass of gravitational theories with nonminimal coupling was presented in [81], encompassing a generalinteraction via the complete set of nine parity-even curvature invariants. The equations of motionfor spinning test bodies in such theories were explicitly derived by means of Synge’s expansiontechnique. The authors’ findings generalize previous results in the literature and allowed for a directcomparison to the general relativistic equations of motion of pole-dipole test bodies. In [82], the authorsderived multipolar equations of motion for gravitational theories with general nonminimal coupling inspacetimes admitting torsion. Their general findings allow for a systematic testing of whole classes oftheories by means of extended test bodies. One peculiar feature of certain subclasses of nonminimaltheories turns out to be their sensitivity to post-Riemannian spacetime structures, even in experimentswithout microstructured test matter. Weak field constraints have also been studied in detail [83,84], andwormhole solutions were further explored, where it is the higher curvature coupling terms that supportthese exotic geometries [85–87].In a cosmological context, the perturbation equation of matter on subhorizon scales was deduced, andspecific bounds on the theory from weak lensing observations and the primordial nucleosynthesis wereobtained in order to constrain the parameters of the model [88]. It was also shown that a non-minimalcoupling between the scalar curvature and the matter Lagrangian density may account for the acceleratedexpansion of the Universe [89–91] and provide, through mimicking, for a viable unification of darkenergy and dark matter [89]. It was shown that a generalized non-minimal coupling between curvatureand matter is compatible with Starobinsky inflation and leads to a successful process of preheating [92],and the problem of a cosmological constant was further explored [93].The effects of the non-minimal curvature-matter coupling on the evolution of cosmologicalperturbations around a homogeneous and isotropic Universe and, hence, the formation of large-scalestructure have also been analyzed [94]. This framework places constraints on the terms, which arisedue to the coupling with matter and, in particular, on the modification in the growth of matter densityperturbations. Approximate analytical solutions were obtained for the evolution of matter overdensitiesduring the matter dominated era, and it was shown that these favor the presence of a coupling functionthat is compatible with the late-time cosmic acceleration.The observations related to the growth of matter has also shown that there is a small, but finite window,where one can distinguish the non-minimally-coupled f ( R ) models with the concordance Λ CDM [95].The possibility that the behavior of the rotational velocities of test particles gravitating around galaxiescan be explained in the framework of modified gravity models with nonminimal curvature-mattercoupling has also been extensively explored [42,96,97].In fact, the literature is extremely vast, and rather than enumerate all of the features of thesemodels, we refer the reader to [98,99] for a review on the topic of the linear nonminimalcurvature-matter coupling. alaxies , f ( R, L m ) Gravity
In this section, we generalize the f ( R ) -type gravity models by assuming that the gravitationalLagrangian is given by an arbitrary function of the Ricci scalar R and of the matter Lagrangian L m [15]. This consists in a maximal extension of the Hilbert–Einstein action, and the action takesthe following form: S = (cid:90) f ( R, L m ) √− g d x (62)where f ( R, L m ) is an arbitrary function of the Ricci scalar R , and of the Lagrangian densitycorresponding to matter, L m . The energy-momentum tensor of the matter is defined by Equation (11).Thus, by assuming that the Lagrangian density L m of the matter depends only on the metric tensorcomponents g µν , and not on its derivatives, we obtain T µν = g µν L m − ∂L m /∂g µν , which will beuseful below.Now, varying the action with respect to the metric yields the following field equation: f R ( R, L m ) R µν + ( g µν (cid:3) − ∇ µ ∇ ν ) f R ( R, L m ) −
12 [ f ( R, L m ) − f L m ( R, L m ) L m ] g µν = 12 f L m ( R, L m ) T µν (63)If f ( R, L m ) = R/ L m (the Hilbert–Einstein Lagrangian), we recover the standard Einstein fieldequation of general relativity, R µν − (1 / g µν R = T µν . For f ( R, L m ) = f ( R ) + f ( R ) G ( L m ) , where f and f are arbitrary functions of the Ricci scalar and G a function of the matter Lagrangian density,respectively, we reobtain the field equations of the modified gravity with arbitrary curvature-mattercoupling, considered in [63].The contraction of Equation (63) provides the following relation between the Ricci scalar R , thematter Lagrange density L m and the trace T = T µµ of the energy-momentum tensor: f R ( R, L m ) R + 3 (cid:3) f R ( R, L m ) − f ( R, L m ) − f L m ( R, L m ) L m ] = 12 f L m ( R, L m ) T (64)By eliminating the term ∇ µ ∇ µ f R ( R, L m ) between Equations (63) and (64), we obtain another formof the gravitational field equations as: f R ( R, L m ) (cid:18) R µν − Rg µν (cid:19) + 16 [ f ( R, L m ) − f L m ( R, L m ) L m ] g µν =12 f L m ( R, L m ) (cid:18) T µν − T g µν (cid:19) + ∇ µ ∇ ν f R ( R, L m ) (65)By taking the covariant divergence of Equation (63), with the use of the mathematical identity [50]: ∇ µ (cid:20) f R ( R, L m ) R µν − f ( R, L m ) g µν + ( g µν (cid:3) − ∇ µ ∇ ν ) f R ( R, L m ) (cid:21) ≡ (66)we obtain for the divergence of the energy-momentum tensor T µν , the following equation: ∇ µ T µν = ∇ µ ln [ f L m ( R, L m )] { L m g µν − T µν } = 2 ∇ µ ln [ f L m ( R, L m )] ∂L m ∂g µν (67) alaxies , The requirement of the conservation of the energy-momentum tensor of matter, ∇ µ T µν = 0 , gives aneffective functional relation between the matter Lagrangian density and the function f L m ( R, L m ) : ∇ µ ln [ f L m ( R, L m )] ∂L m ∂g µν = 0 (68)Thus, once the matter Lagrangian density is known, by an appropriate choice of thefunction f ( R, L m ) , one can construct, at least in principle, conservative models with arbitrarycurvature-matter dependence.Now, assuming that the matter Lagrangian is a function of the rest mass density ρ of the matter only,from Equation (67), we obtain explicitly the equation of motion of the test particles in the f ( R, L m ) gravity model as: D x µ ds = U ν ∇ ν U µ = d x µ ds + Γ µνλ U ν U λ = f µ (69)where the world-line parameter s is taken as the proper time, U µ = dx µ /ds is the four-velocity of theparticle, Γ νσβ are the Christoffel symbols associated with the metric and the extra-force f µ is defined as: f µ = −∇ ν ln (cid:20) f L m ( R, L m ) dL m ( ρ ) dρ (cid:21) ( U µ U ν − g µν ) (70)Note that as in the linear curvature-matter coupling, the extra-force f µ , generated by thecurvature-matter coupling, is perpendicular to the four-velocity, f µ U µ = 0 . Due to the presence of theextra-force f µ , the motion of the test particles in modified theories of gravity with an arbitrary couplingbetween matter and curvature is non-geodesic. From the relation U µ ∇ ν U µ ≡ , it follows that the force f µ is always perpendicular to the velocity, so that U µ f µ = 0 .Building on this analysis, the geodesic deviation equation, describing the relative accelerations ofnearby particles, and the Raychaudhury equation, giving the evolution of the kinematical quantitiesassociated with deformations (expansion, shear and rotation) were extensively considered in thisframework of modified theories of gravity with an arbitrary curvature-matter coupling, by taking intoaccount the effects of the extra force [100,101]. As a physical application of the geodesic deviationequation, the modifications of the tidal forces due to the supplementary curvature-matter coupling wereobtained in the weak field approximation. The tidal motion of test particles is directly influenced notonly by the gradient of the extra force, which is basically determined by the gradient of the Ricci scalar,but also by an explicit coupling between the velocity and the Riemann curvature tensor. As a specificexample, the expression of the Roche limit (the orbital distance at which a satellite will begin to betidally torn apart by the body it is orbiting) was also obtained for this class of models. These aspects willbe presented below.Furthermore, the energy conditions and cosmological applications were explored [102].In [103], the Wheeler–DeWitt equation of f ( R, L m ) gravity was analysed in a flat FRW(Friedmann-Robertson-Walker) universe, which is the first step of the study of quantum f ( R, L m ) cosmology. In the minisuperspace spanned by the FRW scale factor and the Ricci scalar, the equivalenceof the reduced action was examined, and the canonical quantization of the f ( R, L m ) model wasundertaken and the corresponding Wheeler–DeWitt equation derived. The introduction of the invariantcontractions of the Ricci and Riemann tensors were further considered, and applications to black holeand wormhole physics were analyzed [104]. alaxies , f ( R, L m ) Gravity
One of the basic predictions of the modified gravity theories with a curvature-matter coupling is theexistence of an extra force, which makes the motion of the test particles non-geodesic. The existenceof this force can also be tested at the level of the Solar System, by estimating its effects on the orbitalparameters of the motion of the planets around the Sun. The impact on the planetary motion of the extraforce can be obtained in a simple way by using the properties of the Runge–Lenz vector, defined as: (cid:126)A = (cid:126)v × (cid:126)L − α(cid:126)e r (71)where (cid:126)v is the velocity of the planet of mass m relative to the Sun, with mass M (cid:12) , (cid:126)r = r(cid:126)e r the two-bodyposition vector, (cid:126)p = µ(cid:126)v the relative momentum, µ = mM (cid:12) / ( m + M (cid:12) ) the reduced mass: (cid:126)L = (cid:126)r × (cid:126)p = µr ˙ θ(cid:126)k (72)is the angular momentum, and α = GmM (cid:12) , respectively (see [63] and the references therein). For anelliptical orbit of eccentricity e , major semi-axis a and period T , the equation of the orbit is given by ( L /µα ) r − = 1 + e cos θ . The Runge–Lenz vector can be expressed as: (cid:126)A = (cid:32) (cid:126)L µr − α (cid:33) (cid:126)e r − ˙ rL(cid:126)e θ (73)and its derivative with respect to the polar angle θ is given by d (cid:126)A/dθ = r [ dV ( r ) /dr − α/r ] (cid:126)e θ ,where V ( r ) is the potential of the central force [63].The gravitational potential term acting on a planet consists of the Post–Newtonian potential: V P N ( r ) = − αr − α mr (74)plus the gravitational term induced by the general coupling between matter and geometry. Thus, we have: d (cid:126)Adθ = r (cid:20) α mr + m(cid:126)a E ( r ) (cid:21) (cid:126)e θ (75)where we have also assumed that µ ≈ m . Then, we obtain the change of ∆ φ of the perihelion with achange of θ of π as: ∆ φ = (cid:18) αe (cid:19) (cid:90) π (cid:12)(cid:12)(cid:12)(cid:12)(cid:12) ˙ (cid:126)L × d (cid:126)Adθ (cid:12)(cid:12)(cid:12)(cid:12)(cid:12) dθ (76)which can be explicitly calculated as: ∆ φ = 24 π (cid:16) aT (cid:17) − e + L π me (1 − e ) / ( a/T ) (cid:90) π a E (cid:2) L (1 + e cos θ ) − /mα (cid:3) (1 + e cos θ ) cos θdθ (77)where we have used the relation α/L = 2 π ( a/T ) / √ − e . The first term of Equation (77) givesthe expression of the standard general relativistic precession of the perihelion of the planets, while thesecond term gives the contribution to the perihelion precession due to the presence of the extra forcegenerated by the coupling between matter and curvature. alaxies , As an example of the application of Equation (77), we consider the simple case for whichthe extra force may be considered as a constant, a E ≈ constant, an approximation that couldbe valid for small regions of space-time. This case may also correspond to a MOND-typeacceleration a E ≈ √ a a N = √ GM (cid:12) a /r , where a is a constant acceleration, which was proposedphenomenologically as a dynamical model for dark matter [105,106]. In the Newtonian limit, theextra-acceleration generated by the curvature-matter coupling can be expressed in a similar form [9].With the use of Equation (77), one finds for the perihelion precession of a planet in the Solar Systemthe expression: ∆ φ = 6 πGM (cid:12) a (1 − e ) + 2 πa √ − e GM (cid:12) a E (78)where we have also used Kepler’s third law, T = 4 π a /GM (cid:12) . For the planet Mercury a = 57 . × cm, and e = 0 . , respectively, while M (cid:12) = 1 . × g [19]. With theuse of these numerical values, the first term in Equation (78) gives the standard general relativistic valuefor the precession angle, (∆ φ ) GR = 42 . arcsec per century. On the other hand, the observed valueof the precession of the perihelion of Mercury is (∆ φ ) obs = 43 . ± . arcsec per century [19].Therefore, the difference (∆ φ ) E = (∆ φ ) obs − (∆ φ ) GR = 0 . arcsec per century can be attributed toother physical effects. Hence, the observational constraint on the perihelion precession of the planetMercury requires that the value of the constant acceleration a E induced by the curvature-matter couplingat the scale of the Solar System must satisfy the condition a E ≤ . × − cm/s . This value of a E ,obtained from the high precision Solar System observations, is somewhat smaller than the value of theextra acceleration a ≈ − cm/s , necessary to explain the “dark matter” properties, as well as thePioneer anomaly [9]. However, this value does not rule out the possibility of the presence of some extracurvature-matter coupling-type gravitational effects, acting at both the Solar System and galactic levels,since the assumption of a constant extra acceleration may not be correct on larger astronomical scales. f ( R, L m ) Gravity
In the present subsection, we present the equation of the geodesic deviation and the Raychaudhuryequation in f ( R, L m ) gravity, which explicitly contain the effects of the curvature-matter coupling andof the extra force. Some of the physical implications of the geodesic deviation equation, namely, theproblem of the tidal forces in f ( R, L m ) gravity is also considered, and the generalization of the Rochelimit will be considered in the next subsection. In our presentation, we closely follow [100,101].Consider a one-parameter congruence of curves x µ ( s ; λ ) , so that for each λ = λ = constant, x µ ( s, λ ) satisfies Equation (69). We suppose the parametrization to be smooth, and hence, we canintroduce the tangent vector fields along the trajectories of the particles as U µ = ∂x µ ( s ; λ ) /∂s and n µ = ∂x µ ( s ; λ ) /∂λ , respectively. We also introduce the four-vector: η µ = (cid:20) ∂x µ ( s ; λ ) ∂λ (cid:21) δλ ≡ n µ δλ (79)joining points on infinitely close geodesics, corresponding to parameter values λ and λ + δλ , which havethe same value of s . By taking into account Equation (69), we obtain the geodesic deviation equation alaxies , (Jacobi equation), giving the second order derivative with respect to the parameter s of the deviationvector η µ as [67,107]: D η µ ds = R µναβ η α U β U ν + η α ∇ α f µ (80)In the case f µ ≡ , we reobtain the standard Jacobi equation, corresponding to the geodesic motionof test particles in standard general relativity. The interest in the deviation vector η µ derives from thefact that if x µ ( s ) = x µ ( s ; λ ) is a solution of Equation (80), then to first order x µ ( s ) = x µ ( s ) + η µ is asolution of the geodesic equation, as well, since x µ ( s ; λ ) ≈ x µ ( s ; λ )+ n µ ( s ; λ ) δλ ≈ x µ ( s ; λ + δλ ) .By taking into account the explicit form of the extra force given by Equation (70), in f ( R, L m ) -modified theories of gravity with a curvature-matter coupling, the geodesic deviation equationcan be written as: D η µ ds = R µναβ η α U β U ν + η α ∇ α (cid:26) ∇ ν ln (cid:20) f L m ( R, L m ) dL m ( ρ ) dρ (cid:21) ( L m g µν − T µν ) (cid:27) (81)Explicitly, the geodesic deviation equation becomes: D η µ ds = R µναβ η α U β U ν + η α (cid:26) ∇ α ∇ ν ln (cid:20) f L m ( R, L m ) dL m ( ρ ) dρ (cid:21)(cid:27) ( L m g µν − T µν )+ η α ∇ ν ln (cid:20) f L m ( R, L m ) dL m ( ρ ) dρ (cid:21) ( g µν ∇ α L m − ∇ α T µν ) (82)In the presence of an extra force, the Raychaudhury equation is obtained as [107]: ˙ θ + 13 θ + (cid:0) σ − ω (cid:1) = ∇ µ f µ + R µν U µ U ν (83)where θ = ∇ ν U ν is the expansion of the congruence of particles, σ = σ µν σ µν and ω = ω µν ω µν , respectively.With the explicit use of the field equation Equation (63) and the expression of the extra force given byEquation (70), in f ( R, L m ) -modified theories of gravity with an arbitrary coupling between curvatureand matter, the Raychaudhury equation assumes the following generalised form: ˙ θ = − θ − (cid:0) σ − ω (cid:1) + Λ ( R, L m ) + ∇ µ (cid:26) ∇ ν ln (cid:20) f L m ( R, L m ) dL m ( ρ ) dρ (cid:21) ( L m g µν − T µν ) (cid:27) + 1 f R ( R, L m ) U µ U ν ∇ µ ∇ ν f R ( R, L m ) + Φ ( R, L m ) (cid:18) T µν U µ U ν − T (cid:19) (84)where we have denoted: Λ (
R, L m ) = 2 f R ( R, L m ) R − f ( R, L m ) + f L m ( R, L m ) L m f R ( R, L m ) (85)and: Φ (
R, L m ) = f L m ( R, L m ) f R ( R, L m ) (86)respectively. alaxies , Tides are common astrophysical phenomena, and they are due to the presence of a gradient of thegravitational force field induced by a mass above an extended body or a system of particles. In the SolarSystem, tidal perturbations act on compact celestial bodies, such as planets, moons and comets. On largerscales than the Solar System, as, for example, in a galactic or cosmological context, one can observetidal deformations or disruptions of a stellar cluster by a galaxy or in galaxy encounters [108–110].In the relativistic theories of gravitation, as well as in Newtonian gravity, a local system of coordinatescan be chosen, which is inertial, except for the presence of the tidal forces. In strong gravitationalfields, relativistic tidal effects can lead to important physical phenomena, such as the emission of tidalgravitational waves [108–110].In the following, we denote by a prime the reference frame in which all of the Christoffel symbolsvanish. In such a system, one can always take the deviation vector component η (cid:48) = 0 , which meansthat the particle accelerations are compared at equal times. η (cid:48) i is then the displacement of the particlefrom the origin [111]. Moreover, in the static/stationary case, in which the metric, the Ricci scalar andthe thermodynamic parameters of the matter do not depend on time, we have f (cid:48) = 0 . With the use ofthe equation of motion, this condition implies U (cid:48) = constant = 1 . Therefore, with these assumptions,in a static or stationary space-time, the equation of the geodesic deviation (the Jacobi equation) takes theform [100]: d η (cid:48) i dt (cid:48) = R (cid:48) i l η (cid:48) l + R (cid:48) ijlm η (cid:48) l U (cid:48) j U (cid:48) m + η (cid:48) l ∂f (cid:48) i ∂x (cid:48) l (87)Equation (87) can be reformulated as: F (cid:48) i = d η (cid:48) i dt (cid:48) = K ij η (cid:48) j (88)where F (cid:48) i is the tidal force, and we have also introduced the generalized tidal matrix K il [108–110],which is defined as: K ij = R (cid:48) i j + R (cid:48) ikjm U (cid:48) k U (cid:48) m + ∂f (cid:48) i ∂x (cid:48) j (89)The tidal force has the property: ∂F (cid:48) i ∂η (cid:48) j = K ij (90)and its divergence is given by: ∂F (cid:48) i ∂η (cid:48) i = K (91)where the trace K of the tidal matrix is: K = K ii = R (cid:48) + R (cid:48) km U (cid:48) k U (cid:48) m + ∂f (cid:48) j ∂x (cid:48) j (92) alaxies , With the use of the gravitational field equations Equation (63), we can express K as: K = Λ ( R (cid:48) , L (cid:48) m ) η + 1 f R ( R (cid:48) , L m ) ∂ ∂t (cid:48) f R ( R (cid:48) , L (cid:48) m ) + Φ ( R (cid:48) , L (cid:48) m ) (cid:18) T (cid:48) − T (cid:48) η (cid:19) +Λ ( R (cid:48) , L (cid:48) m ) η km U (cid:48) k U (cid:48) m + 1 f R ( R (cid:48) , L (cid:48) m ) U (cid:48) k U (cid:48) m ∂ ∂x (cid:48) k ∂x (cid:48) m f R ( R (cid:48) , L (cid:48) m )+ φ ( R (cid:48) , L (cid:48) m ) (cid:18) T (cid:48) km − T η km (cid:19) U (cid:48) k U (cid:48) m + ∂∂x (cid:48) k (cid:26) ∂∂x (cid:48) m ln (cid:20) f L m ( R (cid:48) , L (cid:48) m ) dL (cid:48) m ( ρ ) dρ (cid:21) (cid:0) L (cid:48) m η km − T (cid:48) km (cid:1)(cid:27) (93)Then, since in the Newtonian limit, one can omit all time derivatives, we obtain for R i l theexpression [111]: R i l = ∂ φ∂x i ∂x l (94)For simplicity, in the following, we omit the primes for the geometrical and physical quantities in theNewtonian approximation. Therefore, in f ( R, L m ) -modified gravity with a curvature-matter coupling,we obtain for the tidal acceleration of the test particles the expression: d η i dt = F i ≈ ∂ φ∂x i ∂x l η l + R ijlm η l V j V m + η l ∂f i ∂x l (95)where V j and V m are the Newtonian three-dimensional velocities. In the Newtonian approximation, inmodified theories of gravity with a curvature-matter coupling, the tidal force tensor is defined as: ∂F i ∂η l = ∂ φ∂x i ∂x l + R ijlm V j V m + ∂f i ∂x l (96)and its trace gives the generalized Poisson equation: ∂F i ∂η i = ∆ φ + R jm V j V m + ∂f i ∂x i (97)In Newtonian gravity: χ il = − ∂ φ∂x i ∂x l (98)represents the Newtonian tidal tensor [111]. In the Newtonian approximation, the spherical potentialof a given particle with mass M is φ ( r ) = − M/ πr . By choosing a frame of reference, so that the x -axis passes through the particle’s position, corresponding to the radial spherical coordinate, that is, ( x = r, y = 0 , z = 0) , the Newtonian tidal tensor is diagonal and has the only non-zero components: χ ii = diag (cid:18) M πr , − M πr , − M πr (cid:19) (99)The Newtonian tidal force (cid:126)F t can be written as F tx = 2 M ∆ x/ πr , F ty = − GM ∆ y/ πr and F tz = − M ∆ z/ πr , respectively [111]. These results can be used to derive the generalization of theRoche limit in modified gravity with an arbitrary coupling between matter and curvature.The Roche limit is the closest distance r Roche that a celestial object with mass m , radius R m anddensity ρ m , held together only by its own gravity, can come to a massive body of mass M , radius R M alaxies , and density ρ M , respectively, without being pulled apart by the massive object’s tidal (gravitational)force [108–110]. For simplicity, we will consider M (cid:29) m , so that the center of mass of the systemnearly coincides with the geometrical center of the mass M .The elementary Newtonian theory of this process is as follows. Consider a small mass ∆ m located atthe surface of the small object of mass m . There are two forces acting on ∆ m , the gravitational attractionof the mass m , given by: F G = m ∆ m πR m (100)and the tidal force exerted by the massive object, which is given by: F tr = M ∆ mR πr (101)where r is the distance between the centers of the two celestial bodies. The Roche limit is reached atthe distance r = r Roche , when the gravitational force and the tidal force exactly balance each other, F G = F tr , thus giving [108–110]: r Roche = R m (cid:18) Mm (cid:19) / = 2 / R M (cid:18) ρ M ρ m (cid:19) / (102)In modified gravity with a curvature-matter coupling, the equilibrium between gravitational and tidalforces occurs at a distance r Roche given by the equation: (cid:18) M πr Roche + R rjrm V j V m + ∂f r ∂r (cid:19) R m = m πR m + f r (103)where f r is the radial component of the extra-force, which modifies the Newtonian gravitational force,and the curvature tensor R rjrm (no summation upon r ) must be evaluated in the coordinate system inwhich the Newtonian tidal tensor is diagonal. Hence, we obtain the generalized Roche limit in thepresence of arbitrary curvature-matter coupling as [100]: r Roche ≈ R m (cid:18) Mm (cid:19) / (cid:20) πR m m (cid:18) R rjrm V j V m + ∂f r ∂r (cid:19) − πR m m f r (cid:21) (104)where we have assumed that the gravitational effects due to the coupling between matter and curvatureare small as compared to the Newtonian ones. Relativistic corrections to the Newtonian tidalaccelerations caused by a massive rotating source, such as, for example, the Earth, could be determinedexperimentally, at least in principle, thus leading to the possibility of testing relativistic theories ofgravitation by measuring such effects in a laboratory.
4. Extended f ( R, L m ) Gravity with the Generalized Scalar Field and Kinetic Term Dependencies
We generalize f ( R, L m ) gravity, outlined in the previous section, by considering a novel gravitationalmodel with an action given by an arbitrary function of the Ricci scalar, the matter Lagrangian density, alaxies , a scalar field and a kinetic term constructed from the gradients of the scalar field, respectively [24].The action is given by: S = (cid:90) f ( R, L m , φ, g µν ∇ µ φ ∇ ν φ ) √− g d x (105)where √− g is the determinant of the metric tensor g µν , and f ( R, L m , φ, g µν ∇ µ φ ∇ ν φ ) is an arbitraryfunction of the Ricci scalar R , the matter Lagrangian density, L m , a scalar field φ and the gradientsconstructed from the scalar field, respectively. The only restriction on the function f is to be an analyticalfunction of R , L m , φ , and of the scalar field kinetic energy, respectively, that is, f must possess a Taylorseries expansion about any point.Note that one may motivate the introduction of the action Equation (105) by the extensive interestin the literature between couplings between the scalar and the matter fields [112–126]. Indeed, suchcouplings generically appear in Kaluza–Klein theories with compactified dimensions [127] or in thelow energy effective limit of string theories [112,113,115–117], where the dilaton has been proposedto be a good candidate for quintessence [112] or the inflaton [113]. The theory represented by actionEquation (105) offers a generalization of the above theories in a single theoretical framework.Here, we introduce first the “reduced” energy-momentum tensor τ µν of the matter, which is definedas [67], τ µν = − √− g δ ( √− gL m ) δg µν = g µν L m − δL m δg µν (106)In addition to this, we assume that the scalar field φ is independent of the metric, i.e. , δφ/δg µν ≡ ,and in the following, we define, for simplicity, ( ∇ φ ) = g µν ∇ µ φ ∇ ν φ .By varying the action S with respect to the metric tensor and to φ , we obtain the followingfield equation: f R R µν + (cid:0) g µν ∇ λ ∇ λ − ∇ µ ∇ ν (cid:1) f R −
12 ( f − f L m L m ) g µν = 12 f L m τ µν − f ( ∇ φ ) ∇ µ φ ∇ ν φ (107)and the evolution equation for the scalar field: (cid:3) ( ∇ φ ) φ = 12 f φ (108)respectively. The subscript of f denotes a partial derivative with respect to the arguments, i.e. , f R = ∂f /∂R , f L m = ∂f /∂L m , f ( ∇ φ ) = ∂f /∂ ( ∇ φ ) , f φ = ∂f /∂φ , and: (cid:3) ( ∇ φ ) = 1 √− g ∂∂x µ (cid:20) f ( ∇ φ ) √− gg µν ∂∂x ν (cid:21) (109)is the generalized d’Alembert operator of f ( R, L m , φ, ∇ µ φ ∇ µ φ ) gravity.The contraction of Equation (107) provides the following relation between the Ricci scalar R , thematter Lagrangian density L m , the derivatives of the scalar field and the trace τ = τ µµ of the “reduced”energy-momentum tensor: f R R + 3 ∇ µ ∇ µ f R − f − f L m L m ) = 12 f L m τ − f ( ∇ φ ) ∇ µ φ ∇ µ φ (110)By taking the covariant divergence of Equation (107), we obtain for the covariant divergence of the“reduced” energy-momentum tensor the following expression: ∇ σ ( f L m τ µσ ) = 12 ( L m ∇ µ f L m − f φ ∇ µ φ ) + f ( ∇ φ ) ∇ µ φ ∇ σ ∇ σ φ + ∇ µ φ ∇ σ φ ∇ σ f ( ∇ φ ) (111) alaxies , This relationship was deduced by taking into account the following mathematical identities: ∇ µ R µν = 12 ∇ ν R , ( ∇ ν (cid:3) − (cid:3) ∇ ν ) f R = − ( ∇ µ f R ) R µν (112)and considering torsion-free space-times, such that [ ∇ σ ∇ (cid:15) − ∇ (cid:15) ∇ σ ] ψ = 0 , where ψ is any scalar-field.Now, using Equation (109), we get: ∇ σ ( f L m τ µσ ) = L m ∇ µ f L m (113)For φ ≡ , Equation (107) reduces to the field equations of the f ( R, L m ) model considered inthe previous section [15]. For φ (cid:54) = 0 , one recovers the conservation equations for either GeneralRelativity and Brans–Dicke-like scalar-tensor theories (with and without scalar/matter coupling [128]).For instance, the total Lagrangian of the simplest matter-scalar field-gravitational field theory, with ascalar field kinetic term and a self-interacting potential V ( φ ) , corresponds to the choice: f = R L m + λ g µν ∇ µ φ ∇ ν φ + V ( φ ) (114)where λ is a constant. The corresponding field equations can be immediately obtained fromEquation (107) as: R µν − Rg µν = τ µν − λ ∇ µ φ ∇ ν φ + (cid:20) λ g αβ ∇ α φ ∇ β φ + V ( φ ) (cid:21) g µν (115)The scalar field satisfies the evolution equation: √− g ∂∂x µ (cid:20) √− gg µν ∂φ∂x ν (cid:21) = 1 λ dV ( φ ) dφ (116)while the energy-momentum tensor obeys the conservation equation ∇ σ τ µσ = 0 . As an example of the application of the formalism developed in the previous section, we consider asimple phenomenological model, in which a scalar field is non-minimally coupled to pressureless matterwith rest mass density ρ . For the action of the system, we consider: S = (cid:90) (cid:20) R − F ( φ ) ρ + λg µν ∇ µ φ ∇ ν φ (cid:21) √− gd x (117)where F ( φ ) is an arbitrary function of the scalar field that couples non-minimally to ordinary matter.The field equations for this model are given by: R µν − Rg µν = F ( φ ) ρU µ U ν + 2 λ (cid:20) ∇ µ φ ∇ ν φ − g µν ∇ α φ ∇ α φ (cid:21) (118)where U µ is the four-velocity of the matter fluid. The scalar field satisfies the evolution equation: (cid:3) φ = − λ dF ( φ ) dφ ρ (119) alaxies , where (cid:3) is the usual d’Alembert operator defined in a curved space. The total energy-momentum tensorof the scalar field-matter system is given by: T µν = F ( φ ) ρU µ U ν + 2 λ (cid:20) ∇ µ φ ∇ ν φ − g µν ∇ α φ ∇ α φ (cid:21) (120)The Bianchi identities imply that ∇ ν T µν = 0 , and in the following, we assume that the mass densitycurrent is conserved, i.e. , ∇ ν ( ρU ν ) = 0 . Using the latter condition and the mathematical identity givenby [ ∇ σ ∇ (cid:15) − ∇ (cid:15) ∇ σ ] ψ = 0 , we have: F ( φ ) ρ ( U ν ∇ ν U µ ) + ρU µ U ν dF ( φ ) dφ ∇ ν φ + 2 λ ( ∇ µ φ ) (cid:3) φ = 0 (121)Then, by eliminating the term (cid:3) φ with the help of Equation (119), we obtain: U ν ∇ ν U µ + (cid:20) ddφ ln F ( φ ) (cid:21) ( U µ U ν ∇ ν φ − ∇ µ φ ) = 0 (122)Using the identity U ν ∇ ν U µ ≡ d x µ ds + Γ µαβ u α u β , where Γ µαβ are the Christoffel symbols correspondingto the metric, the equation of motion of the test particles non-minimally coupled to an arbitrary scalarfield takes the form: d x µ ds + Γ µαβ U α U β + (cid:20) ddφ ln F ( φ ) (cid:21) ( U µ U ν ∇ ν φ − ∇ µ φ ) = 0 (123)A particular model can be obtained by assuming that F ( φ ) is given by a linear function: F ( φ ) = Λ + 12 (cid:20) − φ (cid:21) (124)where Λ is a constant. Then, the equation of motion becomes: d x µ ds + Γ µαβ U α U β + ( U µ U ν − g µν ) ∇ ν ln (cid:20) − φ (cid:21) = 0 (125)In order to simplify the field equations, we adopt for λ the value λ = − (Λ − / . Then,Equation (119), determining the scalar field, takes the simple form (cid:3) φ = ρ .The gravitational field equations take the form: R µν − g µν R = Λ + 12 T µν (126)with the total energy-momentum tensor given by: T µν = (cid:20) − φ (cid:21) ρU µ U ν − Λ − (cid:20) ∇ µ φ ∇ ν φ − g µν ∇ α φ ∇ α φ (cid:21) (127)For Λ = 1 , we reobtain the general relativistic model for dust. Other possible choices of the function F ( φ ) , such as F ( φ ) = exp( φ ) , can be discussed in a similar way.A more general model can be obtained by adopting for the matter Lagrangian the generalexpression [66,68,128]: L m = − (cid:20) ρ + ρ (cid:90) dp ( ρ ) ρ − p ( ρ ) (cid:21) (128) alaxies , where ρ is the rest-mass energy density and p is the thermodynamic pressure, which, by assumption,satisfies a barotropic equation of state, p = p ( ρ ) . By assuming that the matter Lagrangian doesnot depend on the derivatives of the metric and that the particle matter fluid current is conserved( ∇ ν ( ρu ν ) = 0 ), the Lagrangian given by Equation (128) is the unique matter Lagrangian that can beconstructed from the thermodynamic parameters of the fluid, and it is valid for all gravitational theoriessatisfying the two previously mentioned conditions [128].The gravitational field equations and the equation describing the matter-scalar field coupling aregiven by: R µν − g µν R = F ( φ ) (cid:15) U µ U ν − pg µν + λQ µν (129) (cid:3) φ = 12 λ dF ( φ ) dφ (cid:15) (130)where Q µν = ∇ µ φ ∇ ν φ − ∇ λ φ ∇ λ φg µν and where the total energy density is (cid:15) = ρ + ρ (cid:82) dp/ρ − p [68,128]. With the use of the conservation equation ∇ ν ( ρU ν ) = 0 , one obtainsthe equation of motion of massive test particles as: d x µ ds + Γ µαβ U α U β + ( U µ U ν − g µν ) ∇ ν ln (cid:20) (cid:90) dpρ (cid:21) = 0 (131)Models with scalar field-matter coupling were considered in the frameworkof the Brans–Dicke theory [124], with the action of the model given by S = (cid:82) [ φR/ ω/φ ) ∇ µ φ ∇ µ φ + F ( φ ) L m ] √− gd x . Such models can give rise to a late-timeaccelerated expansion of the Universe for very high values of the Brans–Dicke parameter ω .Other models with interacting scalar field and matter have been considered in [120]. We emphasize thatthe gravitational theory considered in this work generalizes all of the above models.In conclusion, the general formalism outlined in this work can be extremely useful in a variety ofscenarios, such as, in describing the interaction between dark energy, modeled as a scalar field, anddark matter, or ordinary matter (neutrinos), with or without pressure, matter-scalar field interactions ininflation, as well as in the study of the interactions of the scalar field (representing dark matter and/or darkenergy) and the electromagnetic component in the very early Universe. Moreover, they can provide arealistic description of the late expansion of the Universe, where a possible interaction between ordinarymatter and dark energy cannot be excluded a priori . f ( R, T ) Gravity
In this section, we consider another extension of standard General Relativity, namely, f ( R, T ) -modified theories of gravity, where the gravitational Lagrangian is given by an arbitrary functionof the Ricci scalar R and of the trace of the energy-momentum tensor T [16]. The dependence from T may be induced by exotic imperfect fluids or quantum effects (conformal anomaly). The action takes thefollowing form: S = 116 π (cid:90) f ( R, T ) √− g d x + (cid:90) L m √− g d x (132) alaxies , where f ( R, T ) is an arbitrary function of the Ricci scalar, R , and of the trace T of the energy-momentumtensor of the matter, T µν . For the matter content, we assume that it consists of a fluid that can becharacterized by two thermodynamic parameters only, the energy density and the pressure, respectively.With the use of Equation (17), it follows that the trace of the energy-momentum tensor T can beexpressed as a function of the matter Lagrangian as: T = g µν T µν = 4 L m − g µν δL m δg µν (133)Hence, f ( R, T ) gravity theories can be interpreted as extensions of the f ( R, L m ) -type gravitytheories, with the gravitational action depending not only on the matter Lagrangian only, but also onits variation with respect to the metric.By varying the action S of the gravitational field with respect to the metric tensor components g µν ,we obtain the field equations of the f ( R, T ) gravity model as: f R ( R, T ) R µν − f ( R, T ) g µν + ( g µν (cid:3) − ∇ µ ∇ ν ) f R ( R, T )= 8 πT µν − f T ( R, T ) T µν − f T ( R, T ) Θ µν (134)We have defined the variation of T with respect to the metric tensor as: δ (cid:0) g αβ T αβ (cid:1) δg µν = T µν + Θ µν (135)where: Θ µν ≡ g αβ δT αβ δg µν (136)Note that when f ( R, T ) ≡ f ( R ) , from Equation (134), we obtain the field equations of f ( R ) gravity.Contracting Equation (134) gives the following relation between the Ricci scalar R and the trace T ofthe energy-momentum tensor: f R ( R, T ) R + 3 (cid:3) f R ( R, T ) − f ( R, T ) = 8 πT − f T ( R, T ) T − f T ( R, T ) Θ (137)where we have denoted
Θ = Θ µµ .By eliminating the term (cid:3) f R ( R, T ) between Equations (134) and (137), the gravitational fieldequations can be written in the form: f R ( R, T ) (cid:18) R µν − Rg µν (cid:19) + 16 f ( R, T ) g µν = 8 π (cid:18) T µν − T g µν (cid:19) − f T ( R, T ) (cid:18) T µν − T g µν (cid:19) − f T ( R, T ) (cid:18) Θ µν −
13 Θ g µν (cid:19) + ∇ µ ∇ ν f R ( R, T ) (138)Taking into account the covariant divergence of Equation (134), with the use of the followingmathematical identity [50]: ∇ µ (cid:20) f R ( R, T ) R µν − f ( R, T ) g µν + ( g µν (cid:3) − ∇ µ ∇ ν ) f R ( R, T ) (cid:21) ≡ − g µν f T ( R, T ) ∇ µ T (139) alaxies , we obtain, for the divergence of the energy-momentum tensor T µν , the equation: ∇ µ T µν = f T ( R, T )8 π − f T ( R, T ) (cid:20) ( T µν + Θ µν ) ∇ µ ln f T ( R, T ) + ∇ µ Θ µν − g µν ∇ µ T (cid:21) (140)Once the matter Lagrangian is known, the calculation of the tensor Θ µν yields (we refer the readerto [16] for more details): Θ µν = − T µν + g µν L m − g αβ ∂ L m ∂g µν ∂g αβ (141)More specifically, in the case of the electromagnetic field, the matter Lagrangian is given by: L m = − π F αβ F γσ g αγ g βσ (142)where F αβ is the electromagnetic field tensor, we obtain Θ µν = − T µν . For a massless scalar field φ withLagrangian L m = g αβ ∇ α φ ∇ β φ , we obtain Θ µν = − T µν + (1 / T g µν .The problem of the perfect fluids, described by an energy density ρ , pressure p and four-velocity U µ ,is more subtle, since there is no unique definition of the matter Lagrangian. However, in the present studywe assume that the energy-momentum tensor of the matter is given by T µν = ( ρ + p ) U µ U ν + pg µν , andthe matter Lagrangian can be taken as L m = p . The four-velocity U µ satisfies the conditions U µ U µ = − and U µ ∇ ν U µ = 0 , respectively. Then, with the use of Equation (141), we obtain, for the variation of theenergy-momentum of a perfect fluid, the expression: Θ µν = − T µν + pg µν (143)The choice L m = − ρ leads to a different form for Θ µν , Θ µν = − T µν + ρg µν . With this form of thematter Lagrangian, the gravitational field equations take a different form. The nature of the qualitativeand quantitative differences between the two different versions of the f ( R, T ) gravity, corresponding todifferent matter Lagrangians, is still an open question. In the present Section, we consider a particular class of f ( R, T ) -modified gravity models, obtainedby explicitly specifying the functional form of f . Generally, the field equations also depend, through thetensor Θ µν , on the physical nature of the matter field. Hence, in the case of f ( R, T ) gravity, dependingon the nature of the matter source, for each choice of f , we can obtain several theoretical models,corresponding to different matter models [129].As a specific case of a f ( R, T ) model, consider a correction to the Einstein–Hilbert action given by f ( R, T ) = R + 2 f ( T ) , where f ( T ) is an arbitrary function of the trace of the energy-momentum tensorof matter. The gravitational field equations immediately follow from Equation (134), and is given by: R µν − Rg µν = 8 πT µν − f (cid:48) ( T ) T µν − f (cid:48) ( T )Θ µν + f ( T ) g µν (144)where the prime denotes a derivative with respect to the argument.Consider a perfect fluid, so that Θ µν = − T µν − pg µν , and the field equations become: R µν − Rg µν = 8 πT µν + 2 f (cid:48) ( T ) T µν + [2 pf (cid:48) ( T ) + f ( T )] g µν (145) alaxies , In the case of dust with p = 0 , the gravitational field equations are given by: R µν − Rg µν = 8 πT µν + 2 f (cid:48) ( T ) T µν + f ( T ) g µν (146)These field equations were proposed in [130] to solve the cosmological constant problem.The simplest cosmological model can be obtained by assuming a dust universe ( p = 0 , T = − ρ )and by choosing the function f ( T ) , so that f ( T ) = λT , where λ is a constant. By assuming that themetric of the universe is given by the flat Robertson–Walker metric: ds = − dt + a ( t ) (cid:0) dx + dy + dz (cid:1) (147)the gravitational field equations are given by: a a = (8 π + 3 λ ) ρ (148) aa + ˙ a a = λρ (149)respectively. Thus, this f ( R, T ) gravity model is equivalent to a cosmological model with an effectivedark energy density Λ eff ∝ H , where H = ˙ a/a is the Hubble function [130]. It is also interestingto note that, generally, for this choice of f ( R, T ) , the gravitational coupling becomes an effective andtime-dependent coupling, of the form G eff = G ± f (cid:48) ( T ) . Thus, the term f ( T ) in the gravitationalaction modifies the gravitational interaction between matter and curvature, replacing G by a runninggravitational coupling parameter.The field equations reduce to a single equation for H : H + 3 8 π + 2 λ π + 3 λ H = 0 (150)with the general solution given by: H ( t ) = 2 (8 π + 3 λ )3 (8 π + 2 λ ) 1 t (151)The scale factor evolves according to a ( t ) = t α , with α = 2 (8 π + 3 λ ) / π + 2 λ ) . The f ( R, T ) gravitational model analysed above has been given a great amount of recent attention.We will briefly outline a few cosmological applications. In [131], the cosmological reconstruction of f ( R, T ) gravity describing matter-dominated and accelerated phases was analysed. Special attention waspaid to the specific case of f ( R, T ) = f ( R ) + f ( T ) . The use of an auxiliary scalar field was consideredwith two known examples for the scale factor corresponding to an expanding Universe. In the firstexample, where ordinary matter is usually neglected for obtaining the unification of matter-dominatedand accelerated phases with f ( R ) gravity, it was shown that this unification can be obtained in thepresence of ordinary matter. In the second example, as in f ( R ) gravity, the model of f ( R, T ) gravitywith a transition of the matter-dominated phase to the acceleration phase was obtained.In [132], it was shown that the dust fluid reproduces Λ CDM, he phantom-non-phantom era andphantom cosmology. Furthermore, different cosmological models were reconstructed, including the alaxies , Chaplygin gas and the scalar field with some specific forms of f ( R, T ) . The numerical simulation for theHubble parameter shows good agreement with the baryon acoustic oscillations (BAO) observational datafor low redshifts z < . In [133], through a numerical reconstruction, it was shown that specific f ( R, T ) models are able to reproduce the same expansion history generated, in standard General Relativity, bydark matter and holographic dark energy. It was further shown that these theories are able to reproducethe four known types of future finite-time singularities [134]. In [135], a non-equilibrium picture ofthermodynamics was discussed at the apparent horizon of an FRW universe in f ( R, T ) gravity, wherethe validity of the first and second law of thermodynamics in this scenario were checked. It was shownthat the Friedmann equations can be expressed in the form of the first law of thermodynamics and thatthe second law of thermodynamics holds both in the phantom and non-phantom phases.The energy conditions have also been extensively explored in f ( R, T ) gravity, for instance, in anFRW universe with perfect fluid [136], and in the context of exact power-law solutions [137], wherecertain constraints were found that have to be satisfied to ensure that power law solutions may be stableand match the bounds prescribed by the energy conditions. In [138], it was also shown that the energyconditions were satisfied for specific models. Furthermore, an analysis of the perturbations and stabilitiesof de Sitter solutions and power-law solutions was performed, and it was shown that for those models inwhich the energy conditions are satisfied, de Sitter solutions and power-law solutions may be stable.Solutions of the Gödel universe in the framework of f ( R, T ) -modified theories of gravity were alsoobtained [139]. Algorithms were derived for constructing five-dimensional Kaluza–Klein cosmologicalspace-times in the presence of a perfect fluid source in f ( R, T ) gravity [140]. Spatially homogeneousand anisotropic Bianchi type-V cosmological models in a scalar-tensor theory of gravitation [141] werealso explored. The Bianchi type I model with perfect fluid as a matter content in f ( R, T ) gravity wasanalysed, and the physical and kinematical properties of the model were also discussed [142]. In [143],the specific case of the conservation of the energy-momentum tensor was considered, and cosmologicalsolutions were obtained for a homogeneous and isotropic model of the Universe.In [144], the evolution of scalar cosmological perturbations in the metric formalism were analysed.According to restrictions on the background evolution, a specific model within these theories wasassumed in order to guarantee the standard continuity equation. Using a completely general procedure,the complete set of differential equations for the matter density perturbations was found. In the case ofsub-Hubble modes, the density contrast evolution reduces to a second-order equation, and it was shownthat for well-motivated f ( R, T ) Lagrangians, the quasistatic approximation provide different resultsfrom the ones derived in the frame of the concordance Λ CDM model, which constrains severely theviability of such theories. In [145], the Ricci and modified Ricci dark energy models were considered inthe context of f ( R, T ) gravity, and it was found that specific models can reproduce the expansion historyof the Universe in accordance with the present observational data. The Dolgov–Kawasaki stabilitycondition were also obtained for the specific reconstructed f ( R, T ) functions. alaxies , f ( R, T , R µν T µν ) Gravity
For the specific case of a traceless energy-momentum tensor, T = 0 , for instance, when theelectromagnetic field is involved, the gravitational field equations for the f ( R, T ) theory [16] reduceto that of the field equations for f ( R ) gravity, and all non-minimal couplings of gravity to the matterfield vanish. This motivates a further generalization of f ( R, T ) gravity that consists of including anexplicit first order coupling between the matter energy-momentum T µν and the Ricci tensor [25,26].In contrast to f ( R, T ) gravity, for T = 0 , this extra coupling still has a non-minimal coupling to theelectromagnetic field via the R µν T µν coupling term in the action, which is non-zero in general. As inthe previous section, we consider the matter content as consisting of a perfect fluid characterized by itsenergy density and thermodynamic pressure only.The action, is given by: S = 116 πG (cid:90) d x √− gf ( R, T, R µν T µν ) + (cid:90) d x √− gL m (152)The only requirement imposed on the function f ( R, T, R µν T µν ) is that it is an arbitrary analyticalfunction in all arguments. The gravitational field equations take the following form: ( f R − f RT L m ) G µν + (cid:20) (cid:3) f R + 12 Rf R − f + f T L m + 12 ∇ α ∇ β (cid:0) f RT T αβ (cid:1)(cid:21) g µν −∇ µ ∇ ν f R + 12 (cid:3) ( f RT T µν ) + 2 f RT R α ( µ T αν ) − ∇ α ∇ ( µ (cid:2) T αν ) f RT (cid:3) − (cid:18) f T + 12 f RT R + 8 πG (cid:19) T µν − (cid:0) f T g αβ + f RT R αβ (cid:1) ∂ L m ∂g µν ∂g αβ = 0 (153)The trace of the gravitational field equation, Equation (153), is obtained as: (cid:3) f R + 12 (cid:3) ( f RT T ) + ∇ α ∇ β (cid:0) f RT T αβ (cid:1) + Rf R − T f T − RT f RT + 2 R αβ T αβ f RT + Rf RT L m + 4 f T L m − f − πGT − g µν (cid:0) g αβ f T + R αβ f RT (cid:1) ∂ L m ∂g µν ∂g αβ = 0 (154)The second derivative of the matter Lagrangian with respect to the metric is non-zero if the matterLagrangian is of second or higher order in the metric. Thus, for a perfect fluid with L m = − ρ , or a scalarfield with L m = − ∂ µ φ∂ µ φ/ , this term can be dropped. However, for instance, considering the Maxwellfield, we have L m = − F µν F µν / ; this term results in ∂ L m /∂g µν ∂g αβ = − F µα F νβ / , thus giving anon-zero contribution to the field equations.In analogy with the standard Einstein field equation, one can write the field Equation (153) as: G µν = 8 πG eff T µν − Λ eff g µν + T eff µν (155)where we have defined the effective gravitational coupling G eff , the effective cosmological constant Λ eff and an effective energy-momentum tensor T eff µν as: G eff = G + π (cid:0) f T + f RT R − (cid:3) f RT (cid:1) f R − f RT L m (156) alaxies , Λ eff = 2 (cid:3) f R + Rf R − f + 2 f T L m + ∇ α ∇ β ( f RT T αβ )2( f R − f RT L m ) (157)and: T eff µν = 1 f R − f RT L m (cid:40) ∇ µ ∇ ν f R − ∇ α f RT ∇ α T µν − f RT (cid:3) T µν − f RT R α ( µ T αν ) + ∇ α ∇ ( µ (cid:2) T αν ) f RT (cid:3) + 2 (cid:0) f T g αβ + f RT R αβ (cid:1) ∂ L m ∂g µν ∂g αβ (cid:41) (158)respectively. Note that, in general, G eff and Λ eff are not constants and depend on the specificmodel considered. f ( R, T, R µν T µν ) Gravity Theory
The covariant divergence of the energy-momentum tensor can be obtained by taking the divergenceof the gravitational field equation, Equation (153), which takes the following form: ∇ µ T µν = 2(1 + Rf T R + 2 f T ) (cid:40) ∇ µ ( f RT R σµ T σν ) + ∇ ν ( L m f T ) − (cid:18) f RT R ρσ + f T g ρσ (cid:19) ∇ ν T ρσ − G µν ∇ µ ( f RT L m ) −
12 [ ∇ µ ( Rf RT ) + 2 ∇ µ f T ] T µν (cid:41) (159)where we have assumed that ∂ L m /∂g µν ∂g αβ = 0 and have used the mathematical identities: ∇ µ (cid:18) f R R µν + (cid:3) f R g µν − f g µν − ∇ µ ∇ ν f R (cid:19) = − (cid:20) f T ∇ ν T + f RT ∇ ν ( R ρσ T ρσ ) (cid:21) (160) T µτ ; δ [; ρ ; σ ] = T µτ ; α R αδρσ + T ατ ; δ R αµρσ + T µα ; δ R ατρσ (161)and [ (cid:3) , ∇ ν ] T = R µν ∇ µ T , respectively.In order to find the equation of motion for a massive test particle, as in the previous sections,we consider the energy-momentum tensor of a perfect fluid. Following the procedure outlined above, weobtain the equation of motion for a massive test particle, considering the matter Lagrangian L m = p , as: d x λ ds + Γ λµν U µ U ν = f λ (162)where the extra force acting on the test particles is given by: f λ = 1 ρ + p (cid:104) ( f T + Rf RT ) ∇ ν ρ − (1 + 3 f T ) ∇ ν p − ( ρ + p ) f RT R σρ ( ∇ ν h σρ − ∇ ρ h σν ) − f RT R σρ h σρ ∇ ν ( ρ + p ) (cid:105) h λν f T + Rf RT (163)Contrary to the nonminimal coupling presented in [9] and as can be seen from the above equations,the extra force does not vanish, even with the Lagrangian L m = p .The extra-force is perpendicular to the four-velocity, satisfying the relation f µ U µ = 0 . In the absenceof any coupling between matter and geometry, with f T = f RT = 0 , the extra-force takes the usual alaxies , form of the standard general relativistic fluid motion, i.e. , f λ = − h λν ∇ ν p/ ( ρ + p ) . In the case of f ( R, T, R µν T µν ) gravity, there is an explicit dependence of the extra-force on the Ricci tensor R σρ ,which makes the deviation from the geodesic motion more important for regions with strong curvatures. f = R + αR µν T µν Let us now consider some examples of cosmological solutions of the theory. In order to obtainexplicit results and as the first step, one has to fix the functional form of the function f ( R, T, R µν T µν ) .We analyze the evolution and dynamics of the Universe, assuming that the Universe is isotropic andhomogeneous, with the matter content described by the energy density ρ and thermodynamic pressure p with the matter Lagrangian chosen as L m = − ρ . Consider the Friedmann–Lemaitre–Robertson–Walker(FLRW) metric, the Hubble parameter H = ˙ a/a and the deceleration parameter q , defined as: q = ddt H − (164)Note that if q < , the expansion of the Universe is accelerating, while positive values of q > describe decelerating evolutions.For simplicity, consider the case in which the interaction between matter and geometry takes placeonly via the coupling between the energy-momentum and Ricci tensors, i.e. , f = R + αR µν T µν (165)In order to pass the Solar System and the other astrophysical tests, the correction term inEquation (165) must be small, which implies that α is a small parameter. This simple case serves as anexample to show the main differences of the present theory with f ( R, T ) gravity, considered before [16].The gravitational field equations for this form of f are given by: G µν + α (cid:34) R σ ( µ T σν ) − R ρσ T ρσ g µν − RT µν − (cid:0) ∇ σ ∇ ( ν T σµ ) − (cid:3) T µν − ∇ α ∇ β T αβ g µν (cid:1) − G µν L m − R αβ ∂ L m ∂g µν ∂g αβ (cid:35) − πGT µν = 0 (166)For the case of the FLRW metric, the modified Friedmann equations are: H = κ − αρ ρ + 32 α − αρ H ( ˙ ρ − ˙ p ) (167)and: H + 3 H = 2 α αp H ˙ ρ − κp αp + 12 α αp ( ¨ ρ − ¨ p ) (168)respectively, where we have denoted κ = 8 πG for simplicity. When α = 0 , we recover the standardFriedmann equations. To remove the under determinacy of the field equations, we impose an equationof state for the cosmological matter, p = p ( ρ ) . A standard form of the cosmological matter equation ofstate is p = ωρ , where ω = constant and ≤ ω ≤ . alaxies , α is small, so that αρ (cid:28) and αp (cid:28) ,respectively. In the high-energy limit, ρ = p , then Equations (167) and (168) take the approximate form: H = κρ (169) H + 3 H = − κρ + 2 αH ˙ ρ (170)The time evolution of the Hubble parameter is described by the equation: (cid:18) − ακ (cid:19) H + 3 H = 0 (171)and, hence, for this model, the evolution of the Hubble parameter is given by: H ( t ) = (cid:112) ( C + 3 κt ) − ακ + C + 3 κt α (172)where C > is an arbitrary integration constant. The scale factor of the Universe is given by: a ( t ) = C exp (cid:20) ( C +3 κt ) √ ( C +3 κt ) − ακ +9 κt +6 κC t ακ (cid:21) (cid:113)(cid:112) ( C + 3 κt ) − ακ + C + 3 κt (173)where C > is an arbitrary integration constant. In order to have a physical solution, the integrationconstant C must satisfy the constraint C ≥ √ κα .The values of the integration constant can be determined from the condition H (0) = H , and a (0) = a , where H and a are the initial values of the Hubble parameter and of the scale factorof the Universe, respectively. This condition immediately provides for C the following value C = (6 αH + κ ) /H . For the integration constant C , we obtain: C = a (cid:118)(cid:117)(cid:117)(cid:116)(cid:115) ( κ − αH ) H + 6 αH + κH exp − (cid:114) ( κ − αH ) H (6 αH + κ )72 αH κ (174)In the small time limit, the scale factor is approximated by: a ( t ) ≈ a (cid:18) κ H α t (cid:19) (175)The deceleration parameter is obtained as: q ( t ) = − αH κ (cid:114) ( αH +3 H κt + κ ) H − ακ (cid:34) αH + H (cid:114) ( αH +3 H κt + κ ) H − ακ + 3 H κt + κ (cid:35) − (176) alaxies , and it can be represented in a form of a power series as: q ( t ) ≈ − − αH κ − αH + 6 H κ ( κ − αH ) t (177)For small values of time, if αH (cid:28) κ , q ≈ − , and the Universe starts its expansion from ade Sitter-like phase, entering, after a finite time interval, into a decelerating phase. On the other hand,if κ > αH , q < − , and the non-singular Universe experiences an initial super-accelerating phase.6.3.2. The Case of DustNext, we consider the case of low density cosmological matter, with p = 0 . Moreover, we assumeagain that the condition αρ (cid:28) holds. Then, the gravitational field equations, Equations (167) and (168),corresponding to a FLRW universe, take the approximate form: H = κρ + 32 αH ˙ ρ (178) H + 3 H = 2 αH ˙ ρ + 12 α ¨ ρ (179)First, we consider the matter dominated phase of the model, in which the non-accelerating expansionof the Universe can be described by a power law form of the scale factor, so that a = t m , m = constant and H = m/t , respectively. The deceleration parameter is given by q = 1 /m − . Therefore,Equation (178) gives, for the time evolution of the density, the equation: αm t ˙ ρ + κρ − m t = 0 (180)with the general solution given by: ρ ( t ) = e − κt α (cid:20) ρ αe κt α + Ei (cid:16) t κ α (cid:17) − Ei (cid:16) t κ α (cid:17)(cid:21) α (181)where Ei ( z ) = − (cid:82) ∞− z e − t dt/t is the exponential integral function, and we have used the initialcondition ρ ( t ) = ρ . By substituting the expressions of the density and of the Hubble parameterinto Equation (179), in the first order of approximation, we obtain the following constraint on m , (9 m − m + 1) / t + O ( t ) ≈ , which is (approximately) satisfied if m is given by the algebraicequation m − m + 1 = 0 , having the solutions m = 1 and m = 1 / , respectively. The decelerationparameters corresponding to these solutions are q = 0 and q = 8 , respectively. Since a value of thedeceleration parameter of the order of q = 8 seems to be ruled out by the observations, the physicalsolution has a scale factor a = t , and q = 0 . The cosmological solutions with zero value of thedeceleration parameter are called marginally accelerating, and they describe the pre-accelerating phaseof the cosmic expansion.Now, we look for a de Sitter-type solution of the field equations for the pressureless matter,Equations (178) and (179), by taking H = H = constant . Then, it follows that, in order to havean accelerated expansion, the matter density must satisfy the equation: ¨ ρ − H ˙ ρ + 2 κα ρ = 0 (182) alaxies , with the general solution given by: ρ ( t ) = e H ( t − t ) (cid:40) √ α (2 ρ − H ρ ) (cid:112) αH − κ sinh (cid:34) (cid:112) αH − κ √ α ( t − t ) (cid:35) + cosh (cid:34) (cid:112) αH − κ √ α ( t − t ) (cid:35) (cid:41) (183)where we have used the initial conditions ρ ( t ) = ρ and ˙ ρ ( t ) = ρ , respectively. Therefore, in thepresence of a non-trivial curvature-matter coupling, once the evolution of the matter density is given byEquation (183), the time evolution of the Universe is of the de Sitter-type.
7. Dark Matter As a Curvature-Matter Coupling Effect
According to Newton’s gravitation theory, the rotation of hydrogen clouds around galaxies shouldshow a Keplerian decrease with distance r of the orbital rotational speed v tg outside the luminousbaryonic matter, v tg ∝ M ( r ) /r , where M ( r ) is the dynamical mass of the galaxy. However, observationsshow that the rotation curves are rather flat [146–149]. The rotational velocities increase near the centerof the galaxy, and then, they remain nearly constant at a value of v tg ∞ ∼ M ( r ) ≈ rv tg ∞ /G [146–149]. Consequently, even at largedistances, where very little luminous matter can be detected, the mass within a distance r from the centerof the galaxy increases linearly with the radial distance r .This unusual behavior of the galactic rotation curves is usually explained by postulating the existenceof some dark (invisible) matter, distributed in a spherical halo around the galaxies. The rotationcurves are obtained by measuring the frequency shifts z of the 21-cm radiation emission from theneutral hydrogen gas clouds. Usually, astronomers report the resulting z in terms of a velocity field v tg [146–148]. In the standard cosmological model, dark matter is assumed to be a cold, pressurelessmedium. Many possible candidates for dark matter have been proposed, the most popular ones beingthe weakly interacting massive particles (WIMP) (for a review of the particle physics aspects of darkmatter, see [150]). The interaction cross-section with normal baryonic matter is expected to be small,but non-zero, and we expect that dark matter particles can be detected directly. However, after more than20 years of intense experimental and observational effort, presently, no convincing non-gravitationalevidence for dark matter exists.In the present section, we will briefly analyze, following [42], the possibility that the dynamicbehavior of hydrogen clouds rotating around the galactic center can be explained due to the presenceof a curvature-matter coupling. We will restrict our analysis to the case of modified gravity modelswith linear curvature-matter coupling. The extra-terms in the gravitational field equations modify theequations of motion of test particles and induce a supplementary gravitational interaction, which canaccount for the observed behavior of the galactic rotation curves.Note that the modified gravity approach has been explored, as an alternative to dark matter,in different contexts [151–154], such as in f ( R ) gravity and the recently proposed hybrid metric-Palatinigravity [155–158]. alaxies , In the following, we assume that gas clouds behave like massive test particles, moving in a staticand spherically-symmetric space-time outside the galactic baryonic mass distribution. In the galacticspace-time, we consider two observers O E and O ∞ , with four-velocities u µE and u µ ∞ , respectively.Observer O E is the light emitter ( i.e. , the gas clouds placed at a point P E of the space-time), and O ∞ represents the detector at point P ∞ , located far from the emitter, at “spatial infinity” [159].Without any loss of generality, we can assume that the gas clouds move in the galactic plane θ = π/ ,so that u µE = (cid:16) ˙ t, ˙ r, , ˙ φ (cid:17) E , where the dot stands for derivation with respect to the affine parameter s .On the other hand, we suppose that the detector is static ( i.e. , O ∞ ’s four-velocity is tangent to the staticKilling field ∂/∂t ), and in the chosen coordinate system, its four-velocity is u µ ∞ = (cid:0) ˙ t, , , (cid:1) ∞ [159].In this section, we use the (+ , − , − , − ) signature.The static spherically symmetric metric outside the galactic baryonic mass distribution is given by: ds = e ν ( r ) dt − e λ ( r ) dr − r (cid:0) dθ + sin θdφ (cid:1) (184)where the metric coefficients ν and λ are functions of the radial coordinate r only. The motion of a testparticle in the gravitational field in modified gravity with a linear curvature-matter coupling is describedby the Lagrangian given by Equation (24): L = Q (cid:34) e ν ( r ) (cid:18) dtds (cid:19) − e λ ( r ) (cid:18) drds (cid:19) − r (cid:18) d Ω ds (cid:19) (cid:35) (185)where d Ω = dθ + sin θdφ , and the function Q is given by: Q = [1 + ζf ( R )] dL m ( ρ ) dρ (186)where ζ is a constant.For θ = π/ , d Ω = dφ . From the Lagrange equations, it follows that we have two constants ofmotion, the energy E : E = Qe ν ( r ) ˙ t (187)and the angular momentum l , given by: l = Qr ˙ φ (188)where a dot denotes the derivative with respect to the affine parameter s . The condition of thenormalization of the four-velocity U µ U µ = 1 gives: e ν ( r ) ˙ t − e λ ( r ) ˙ r − r ˙ φ (189)from which, with the use of the constants of motion, we obtain: E = Q e ν + λ ˙ r + e ν (cid:18) l r + Q (cid:19) (190)Equation (190) shows that the radial motion of the particles in modified gravity with linearmatter-geometry coupling is the same as that of a particle in ordinary Newtonian mechanics, with alaxies , velocity ˙ r , position-dependent mass m eff = 2 Q e ν + λ and energy E , respectively, moving in theeffective potential: V eff ( r ) = e ν ( r ) (cid:18) l r + Q (cid:19) (191)The conditions for stable circular orbits ∂V eff /∂r = 0 and ˙ r = 0 determine the energy and theangular momentum of the particle as: l = 12 r Q ( ν (cid:48) Q + 2 Q (cid:48) )1 − rν (cid:48) / (192)and: E = e ν Q ( rQ (cid:48) + Q )1 − rν (cid:48) / (193)respectively.The line element given by Equation (184) can be rewritten in terms of the spatial components of thevelocity, normalized with the speed of light, measured by an inertial observer far from the source, as: ds = dt (cid:0) − v (cid:1) (194)where: v = e − ν (cid:34) e λ (cid:18) drdt (cid:19) + r (cid:18) d Ω dt (cid:19) (cid:35) (195)For a stable circular orbit dr/dt = 0 , and the tangential velocity of the test particle is obtained as: v tg = e − ν r (cid:18) d Ω dt (cid:19) (196)In terms of the conserved quantities E and l , the angular velocity is given, for θ = π/ , by: v tg = e ν r l E (197)With the use of Equations (192) and (193), we obtain: v tg = 12 r ( ν (cid:48) Q + 2 Q (cid:48) ) rQ (cid:48) + Q (198)Thus, the rotational velocity of the test body in modified gravity with linear coupling between matterand geometry is determined by the metric coefficient exp ( ν ) , and by the function Q and its derivativewith respect to the radial coordinate r . In the standard general relativistic limit ζ = 0 , Q = 1 , and weobtain v tg = rν (cid:48) / . alaxies , The velocities of the hydrogen clouds rotating around galactic centers are determined from the redand blue shifts of the radiation emitted by the gas moving on circular orbits on both sides of the centralregion. The radiation travels on null geodesics with tangent k µ . We may restrict without any loss ofgenerality k µ to lie in the equatorial plane θ = π/ and evaluate the frequency shift for a light signalemitted from O E , in motion, in a circular orbit, and detected by O ∞ at infinity. The frequency shiftassociated with the emission and detection of the light signal is given by: z = 1 − ω E ω ∞ (199)where ω I = − k µ u µI , and the index I refers to emission ( I = E ) or detection ( I = ∞ ) at thecorresponding space-time point [159,160]. Two frequency shifts, corresponding to the frequency shiftsof a receding or approaching gas cloud, respectively, are associated with light propagation in the sameand opposite direction of motion of the emitter, respectively. In terms of the tetrads e (0) = e − ν/ ∂/∂t , e (1) = e − λ/ ∂/∂r , e (2) = r − ∂/∂θ , e (3) = ( r sin θ ) − ∂/∂φ , the frequency shifts take the form [159]: z ± = 1 − e [ ν ∞ − ν ( r )] / (1 ∓ v ) Γ (200)where v = (cid:104)(cid:80) i =1 (cid:0) u ( i ) /u (0) (cid:1) (cid:105) / , with u ( i ) the components of the particle’s four velocity along thetetrad ( i.e. , the velocity measured by an Eulerian observer, whose world-line is tangent to the staticKilling field). Γ = (1 − v ) − / is the usual Lorentz factor, and exp ( ν ∞ ) is the value of exp [ ν (( r ))] for r → ∞ . In the case of circular orbits in the θ = π/ plane, we obtain: z ± = 1 − e [ ν ∞ − ν ( r )] / ∓ (cid:112) r ( ν (cid:48) Q + 2 Q (cid:48) ) / rQ (cid:48) + Q ) (cid:112) − r ( ν (cid:48) Q + 2 Q (cid:48) ) / rQ (cid:48) + Q ) (201)It is convenient to define two other quantities, z D = ( z + − z − ) / and z A = ( z + + z − ) / ,respectively [159]. In the modified gravity model with a linear curvature-matter coupling, the redshiftfactors are given by: z D ( r ) = e [ ν ∞ − ν ( r )] / (cid:112) r ( ν (cid:48) Q + 2 Q (cid:48) ) / rQ (cid:48) + Q ) (cid:112) − r ( ν (cid:48) Q + 2 Q (cid:48) ) / rQ (cid:48) + Q ) (202)and: z A ( r ) = 1 − e [ ν ∞ − ν ( r )] / (cid:112) − r ( ν (cid:48) Q + 2 Q (cid:48) ) / rQ (cid:48) + Q ) (203)respectively, which can be easily related to the galactic observations [159]. z A and z D satisfy the relation: ( z A − − z D = exp [2 ( ν ∞ − ν ( r ))] (204)and thus, in principle, by assuming that the metric tensor component exp [ ν (( r ))] is known, Q and Q (cid:48) can be obtained directly from the astrophysical observations. Hence, the observations of the red- andblue-shifts of the radiation emitted by hydrogen clouds rotating around the galactic center could providea direct observational test of the galactic geometry and, implicitly, of the modified gravity models withlinear coupling between matter and geometry. alaxies , The tangential velocities v tg of gas clouds, considered as massive test particles moving around thecenters of the galaxies, cannot be measured directly. Instead, they can be inferred from the redshift z ∞ of the radiation observed at spatial infinity, for which: z ∞ = exp [( ν ∞ − ν ) /
2] (1 ± v tg ) / (cid:113) − v tg (205)Since the velocities of the gas clouds are non-relativistic, with v tg ≤ (4 / × − , we observe v tg ≈ z ∞ (as the first term of a geometric series), with the consequence that the lapse function exp ( ν ) necessarilytends at infinity to unity, i.e. , e ν ≈ e ν ∞ / (cid:0) − v tg (cid:1) ≈ e ν ∞ → . The observations show that at distanceslarge enough from the galactic center, v tg ≈ constant.In the following, we use this observational constraint to reconstruct the coupling term betweencurvature and matter in the “dark matter”-dominated region, far away from the baryonic matterdistribution. By assuming that v tg = constant, Equation (198) can be written as: v tg rQ ddr ( rQ ) = ν (cid:48) Q (cid:48) Q (206)and can be immediately integrated to give: Q ( r ) = (cid:18) rr (cid:19) v tg / ( − v tg ) exp (cid:34) − ν (cid:0) − v tg (cid:1) (cid:35) (207)where r is an arbitrary constant of integration. Since hydrogen clouds are a pressureless dust ( p = 0 )that can be characterized by their density ρ only, the Lagrangian of the matter (gas cloud) is given by L mat ( ρ ) = ρ . Therefore, from Equations (186) and (207), we obtain: ζf ( R ) = (cid:18) rr (cid:19) v tg / ( − v tg ) exp (cid:34) − ν (cid:0) − v tg (cid:1) (cid:35) − ≈ v tg (cid:0) − v tg (cid:1) ln rr − ν (cid:0) − v tg (cid:1) − v tg (cid:0) − v tg (cid:1) ν ln rr (208)For an arbitrary velocity profile v tg = v tg ( r ) , the general solution of Equation (206) is given by: (cid:112) Q ( r ) = 1 + ζf ( R ) = (cid:112) Q exp (cid:20) (cid:90) v tg ( r ) /r − ν (cid:48) / − v tg ( r ) dr (cid:21) (209)where Q is an arbitrary constant of integration.For v tg , we assume the simple empirical dark halo rotational velocity law [161]: v tg = v x a + x (210)where x = r/r opt , r opt is the optical radius containing 83 % of the galactic luminosity. The parameters a , the ratio of the halo core radius and r opt and the terminal velocity v are all functions of thegalactic luminosity L . For spiral galaxies a = 1 . L/L ∗ ) / and v = v opt (1 − β ∗ ) (1 + a ) , where v opt = v tg ( r opt ) , and β ∗ = 0 .
72 + 0 .
44 log ( L/L ∗ ) , with L ∗ = 10 . L (cid:12) . alaxies , One can assume that the coupling between the neutral hydrogen clouds and the geometry is small, ζf ( R ) L mat (cid:28) , and consequently, the galactic geometry is not significantly modified in the vacuumoutside the baryonic mass distribution with mass M B , corresponding to L mat ≈ . Moreover, weassume, for simplicity, that outside the baryonic matter distribution, the galactic metric is given bythe Schwarzschild metric (which is also a solution of the vacuum field equations of the f ( R ) -modifiedgravity [162] and still gives the dominant contribution to the total metric, even if the curvature-mattercoupling is not small), written as: e ν = e − λ = 1 − R x (211)where R = GM B /r opt , from Equation (209), we obtain: ζf ( R ) = exp (cid:34) α × arctanh (cid:32) (cid:112) − v a x (cid:33)(cid:35) × (cid:18) − R x (cid:19) β [(1 − v ) x + a ] γ Q − / x / (212)where: α = − aR v (cid:112) − v [ a + 4 (1 − v ) R ] (213) β = a − R a + 4 (1 − v ) R ] (214)and: γ = − v [ a + 6 (1 − v ) R ]4 (1 − v ) [ a + 4 (1 − v ) R ] (215)respectively. Thus, the geometric part of the curvature-matter coupling can be completely reconstructedfrom the observational data on the galactic rotation curves.As one can see from Equation (198), in the limit of large r , when ν (cid:48) → (in the case of theSchwarzschild metric ν (cid:48) ≈ M B /r ), the tangential velocity of test particles at infinity is given by: v tg ∞ = rr + Q/Q (cid:48) (216)which, due to the presence of the curvature-matter coupling, does not decay to zero at large distancesfrom the galactic center, a behavior that is perfectly consistent with the observational data and that isusually attributed to the existence of dark matter.By using the simple observational fact of the constancy of the galactic rotation curves, thecurvature-matter coupling function can be completely reconstructed, without any supplementaryassumption. If, for simplicity, we consider again that the metric in the vacuum outside the galaxy canbe approximated by the Schwarzschild metric, with exp ( ν ) = 1 − GM B /r , where M B is the massof the baryonic matter of the galaxy, then in the limit of large r , we have ν → . Therefore,fromEquation (208), we obtain: ζ lim r →∞ f ( R ) ≈ v tg (cid:0) − v tg (cid:1) ln rr (217)If the galactic rotation velocity profiles and the galactic metric are known, the coupling function can bereconstructed exactly over the entire mass distribution of the galaxy. alaxies , One can formally associate an approximate “dark matter” mass profile M DM ( r ) with the tangentialvelocity profile, which is determined by the non-minimal curvature-matter coupling and is given by: M DM ( r ) ≈ G r ( ν (cid:48) Q + 2 Q (cid:48) ) rQ (cid:48) + Q (218)The corresponding “dark matter” density profile ρ DM ( r ) can be obtained as: ρ DM ( r ) = 14 πr dMdr = 14 πG × (cid:20) ν (cid:48) Q + 2 Q (cid:48) r ( rQ (cid:48) + Q ) + ν (cid:48)(cid:48) Q + νQ (cid:48) + 2 Q (cid:48)(cid:48) rQ (cid:48) + Q ) − ( ν (cid:48) Q + 2 Q (cid:48) ) ( rQ (cid:48)(cid:48) + 2 Q (cid:48) )2 ( rQ (cid:48) + Q ) (cid:21) (219) Other observational constraints on M DM and ρ DM can be obtained from the study of the galacticstellar populations. We assume that each galaxy consists of a single, pressure-supported stellarpopulation that is in dynamic equilibrium and traces an underlying gravitational potential resulting fromthe non-minimal curvature-matter coupling. In spherical symmetry, the equivalent mass profile inducedby the geometry-matter coupling (the mass profile of the “dark matter” halo) relates to the moments ofthe stellar distribution function via the Jeans equation [149]: ddr (cid:2) ρ s (cid:10) v r (cid:11)(cid:3) + 2 ρ s ( r ) βr = − Gρ s M DM ( r ) r (220)where ρ s ( r ) , (cid:104) v r (cid:105) and β ( r ) = 1 − (cid:104) v θ (cid:105) / (cid:104) v r (cid:105) describe the three-dimensional density, radial velocitydispersion and orbital anisotropy of the stellar component, where (cid:104) v θ (cid:105) is the tangential velocitydispersion. By assuming that the anisotropy is a constant, the Jeans equation has the solution [163]: ρ s (cid:10) v r (cid:11) = Gr − β (cid:90) ∞ r s − β ) ρ s ( s ) M DM ( s ) ds (221)With the use of Equation (218), we obtain for the stellar velocity dispersion equation: ρ s (cid:10) v r (cid:11) ≈ r − β (cid:90) ∞ r s − β ) ρ s ( s ) ν (cid:48) ( s ) Q ( s ) + 2 Q (cid:48) ( s ) sQ (cid:48) ( s ) + Q ( s ) ds (222)After projection along the line of sight, the “dark matter” mass profile can be related to two observableprofiles, the projected stellar density I ( R ) , and to the stellar velocity dispersion σ p ( R ) , according to therelation [149]: σ P ( R ) = 2 I ( R ) (cid:90) ∞ R (cid:18) − β R r (cid:19) ρ s (cid:104) v r (cid:105) r √ r − R dr (223)Given a projected stellar density model I ( R ) , one recovers the three-dimensional stellar densityfrom [149]: ρ s ( r ) = − (1 /π ) (cid:90) ∞ r (cid:18) dIdR (cid:19) (cid:0) R − r (cid:1) − / dR (224)Therefore, once the stellar density profile I ( R ) , the stellar velocity dispersion (cid:104) v r (cid:105) and the galacticmetric are known, with the use of the integral Equation (222), one can obtain the explicit form ofthe curvature-matter coupling function Q and the equivalent mass profile induced by the non-minimal alaxies , coupling between matter and curvature. The simplest analytic projected density profile is the Plummerprofile [149], given by: I ( R ) = L (cid:0) πr half (cid:1) − (cid:0) R /r half (cid:1) − (225)where L is the total luminosity and r half is the projected half-light radius (the radius of the cylinder thatencloses half of the total luminosity). An important physical requirement for the circular orbits of the test particles moving around galaxiesis that they must be stable. Let r be a circular orbit and consider a perturbation of it of the form r = r + δ , where δ (cid:28) r [160]. Taking expansions of V eff ( r ) , exp ( ν + λ ) and Q ( r ) about r = r ,it follows from Equation (190) that: ¨ δ + 12 Q ( r ) e ν ( r )+ λ ( r ) V (cid:48)(cid:48) eff ( r ) δ = 0 (226)The condition for the stability of the simple circular orbits requires V (cid:48)(cid:48) eff ( r ) > [160]. This givesfor the coupling function Q the constraint: d Qdr (cid:12)(cid:12)(cid:12)(cid:12) r = r < (cid:20) ν (cid:48)(cid:48) (cid:18) l r + Q (cid:19) + ν (cid:48) (cid:18) − l r + 2 QQ (cid:48) (cid:19) − l r (cid:21)(cid:12)(cid:12)(cid:12)(cid:12) r = r (227)a condition that must be satisfied at any point r of the galactic space-time.From the observational, as well as from the theoretical point of view, an important problem is toestimate an upper bound for the cutoff of the constancy of the tangential velocities. If in the large r limit, the coupling function satisfies the condition Q (cid:48) /Q → , then Q/Q (cid:48) → ∞ , and in this limit, thetangential velocity decays to zero. If the exact functional form of Q is known, the value of r at whichthe rotational velocity becomes zero can be accurately estimated.
8. Discussion and Conclusions
In the present paper, we have reviewed the cosmological and astrophysical applications of generalizedcurvature-matter couplings in f ( R ) -type gravity models, and we have extensively analysed some of theirapplications. Specific models were explored and presented in detail. The gravitational field equationsin the metric formalism, in the presence of a nonminimal coupling between curvature and matter,were presented, as well as the equations of motion for test particles, which follow from the covariantdivergence of the energy-momentum tensor. Generally, the motion is non-geodesic and takes place inthe presence of an extra force orthogonal to the four-velocity. The Newtonian limit of the equation ofmotion was also described, and a procedure for obtaining the energy-momentum tensor of the matterin the framework of these gravity models was presented. On the other hand, it was shown that thegravitational field equations are equivalent to the effective Einstein equations of the f ( R ) model in emptyspacetime, but differ from them, as well as from standard General Relativity, in the presence of matter.Therefore, the predictions of these gravity theories could lead to some major differences, as compared tothe predictions of standard General Relativity or to its extensions that ignore the role of matter, in several alaxies , problems of current interest, such as cosmology, dark matter, dark energy, gravitational collapse or thegeneration of gravitational waves.We have also reviewed, based on the already existing literature, the implications of thecurvature-matter coupling models on the galactic dynamics, and we have shown that the behavior ofthe neutral hydrogen gas clouds outside of the galaxies and their flat rotation curves can be explainedin terms of a non-minimal coupling between matter and curvature. We have also shown that in thecurvature-matter coupling “dark matter” model, all of the relevant physical quantities, including the“dark mass” associated with the coupling, and which plays the role of dark matter, its correspondingdensity profile, as well as the curvature-matter coupling function can be expressed in terms of observableparameters: the tangential velocity, the baryonic (luminous) mass and the Doppler frequency shiftsof test particles moving around the galaxy. Therefore, this opens the possibility of directly testingthe modified gravity models with non-minimal coupling between matter and geometry by using directastronomical and astrophysical observations at the galactic or extra-galactic scale. Since the observationson the galactic rotation curves are obtained from the Doppler frequency shifts, we have generalized theexpression of the frequency shifts by including the effect of the curvature-matter coupling. Thus, at leastin principle, the coupling function can be obtained directly from astronomical observations. We havealso presented two other possibilities for testing these classes of gravity theories, at the level of theSolar System, by using the perihelion precession of Mercury, and at astrophysical scales, from possibleobservations of the gravitational tidal effects.As future avenues of research, one should aim to characterize as much as possible the phenomenologypredicted by these theories with curvature-matter couplings in order to find constraints arising fromobservations. The study of these phenomena may also provide some specific signatures and effects,which could distinguish and discriminate between the various theories of modified gravity. We alsopropose to use a background metric to analyse the dynamic system for specific curvature-matter couplingmodels and to use the data of SNIa (Type Ia supernovae), BAO and CMB shift parameter to obtainrestrictions for the respective models and to explore in detail the analysis of structure formation. Acknowledgments
Francisco S.N. Lobo acknowledges financial support of the Fundação para a Ciência e Tecnologiathrough an Investigador FCT (Fundação para a Ciência e Tecnologia) Research contract, withReference IF/00859/2012, funded by FCT/MCTES (Fundação para a Ciência e Tecnologia eMinisterio da Ciencia, Tecnologia e Ensino Superior) (Portugal) and grants CERN/FP/123618/2011 andEXPL/FIS-AST/1608/2013.
Author Contributions
Both authors have contributed to the writing of this manuscript through their joint published workover the past few years.
Conflicts of Interest
The authors declare no conflict of interest. alaxies ,
1. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.;Fabbro, S.; Goobar, A.; Groom, D.E.; et al . Measurements of Ω and Λ from 42 High-RedshiftSupernovae. Astrophys. J. , , 565–586.2. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiattia, A.; Diercks, A.; Garnavich, P.M.;Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al . Observational Evidence fromSupernovae for an Accelerating Universe and a Cosmological Constant. Astron. J. , ,1009–1038.3. Riess, A.G.; Strolger, L.-G.; Tonry, J.; Casertano, S.; Ferguson, H.C.; Mobasher, B.; Challis, P.;Filippenko, A.V.; Jha, S.; Li, W.; et al . Type Ia Supernova Discoveries at z > 1 from the HubbleSpace Telescope: Evidence for Past Deceleration and Constraints on Dark Energy Evolution. Astrophys. J. , , 665–687.4. Maartens, R.; Koyama, K. Brane-world gravity. Living Rev. Relativ. , , 7.5. Lobo, F.S.N. The Dark side of gravity: Modified theories of gravity. ArXiv E-Prints , ,arXiv:0807.1640.6. Nojiri, S.; Odintsov, S.D. Unified cosmic history in modified gravity: From F ( R ) theory toLorentz non-invariant models. Phys. Rep. , , 59–144.7. Nojiri, S.; Odintsov, S.D.; Sasaki, M. Gauss-Bonnet dark energy. Phys. Rev. D , ,doi:10.1103/PhysRevD.71.123509.8. Carroll, S.M.; Duvvuri, V.; Trodden, M.; Turner, M.S. Is cosmic speed-up due to newgravitational physics? Phys. Rev. D , , doi:10.1103/PhysRevD.70.043528.9. Bertolami, O.; Boehmer, C.G.; Harko, T.; Lobo, F.S.N. Extra force in f ( R ) modified theories ofgravity. Phys. Rev. D , , doi:10.1103/PhysRevD.75.104016.10. Nojiri, S.; Odintsov, S.D. Gravity assisted dark energy dominance and cosmic acceleration. Phys.Lett. B , , 137–142.11. Nojiri, S.; Odintsov, S.D. Dark energy and cosmic speed-up from consistent modified gravity. PoS WC , , 024.12. Allemandi, G.; Borowiec, A.; Francaviglia, M.; Odintsov, S.D. Dark energy dominanceand cosmic acceleration in first order formalism. Phys. Rev. D , , doi:10.1103/PhysRevD.72.063505.13. Bamba, K.; Odintsov, S.D. Inflation and late-time cosmic acceleration in non-minimalMaxwell- F ( R ) gravity and the generation of large-scalemagnetic fields. J. Cosmol. Astropart.Phys. , , doi:10.1088/1475-7516/2008/04/024.14. Bamba, K.; Nojiri, S.; Odintsov, S.D. Inflationary cosmology and the late-time acceleratedexpansion of the universe in non-minimal Yang-Mills- F ( R ) gravity and non-minimalvector- F ( R ) gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.77.123532.15. Harko, T.; Lobo, F.S.N. f ( R, L m ) gravity. Eur. Phys. J. C , , 373–379.16. Harko, T.; Lobo, F.S.N.; Nojiri, S.; Odintsov, S.D. f ( R, T ) gravity. Phys. Rev. D , ,doi:10.1103/PhysRevD.84.024020.17. Faraoni, V. Cosmology in Scalar-Tensor Gravity ; Springer: Berlin, Germany, 2004. alaxies ,
18. Bertolami, O.; Paramos, J.; Turyshev, S. General theory of relativity: Will it survive the nextdecade? In
Lasers, Clocks, and Drag-Free Control ; Springer: Berlin, Germany, 2006; pp. 27–67.19. Will, C.M. The Confrontation between General Relativity and Experiment.
ArXiv E-Prints , ,arXiv:gr-qc/0510072.20. Bertolami, O.; Pedro, F.G.; le Delliou, M. Dark Energy-Dark Matter Interaction and the Violationof the Equivalence Principle from the Abell Cluster A586. Phys. Lett. B , , 165–169.21. Damour, T. Questioning the equivalence principle. ArXiv E-Prints , , arXiv:gr-qc/0109063.22. Damour, T. Testing the equivalence principle: Why and how? Class. Quantum Gravity , ,doi:10.1088/0264-9381/13/11A/005.23. Damour, T.; Donoghue, J.F. Equivalence Principle Violations and Couplings of a Light Dilaton. Phys. Rev. D , , doi:10.1103/PhysRevD.82.084033.24. Harko, T.; Lobo, F.S.N.; Minazzoli, O. Extended f ( R, L m ) gravity with generalized scalar fieldand kinetic term dependences. Phys. Rev. D , , doi:10.1103/PhysRevD.87.047501.25. Haghani, Z.; Harko, T.; Lobo, F.S.N.; Sepangi, H.R.; Shahidi, S. Further matters in space-timegeometry: f ( R, T, R µν T µν ) gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.88.044023.26. Odintsov, S.D.; Sáez-Gómez, D. f ( R, T, R µ νT µ ν ) gravity phenomenology and Λ CDM universe.
Phys. Lett. B , , 437–444.27. Haghani, Z.; Harko, T.; Sepangi, H.R.; Shahidi, S. Matter may matter. ArXiv E-Prints , ,arXiv:1405.3771.28. Deser, S.; Gibbons, G.W. Born-Infeld-Einstein actions? Class. Quantum Gravity , ,doi:10.1088/0264-9381/15/5/001.29. Dolgov, A.D.; Kawasaki, M. Can modified gravity explain accelerated cosmic expansion? Phys.Lett. B , , 1–4.30. Unzicker, A.; Case, T. Translation of Einstein’s attempt of a unified field theory withteleparallelism. ArXiv E-Prints , , arXiv:physics/0503046.31. Möller, C. Conservation laws and absolute parallelism in general relativity. Math. Phys. Skr.Danske Vid. Selsk. , , 3–50.32. Pellegrini, C.; Plebanski, J. Tetrad fields and gravitational fields. Math. Phys. Skr. Danske Vid.Selsk. , , 1–39.33. Hayashi, K.; Shirafuji, T. New general relativity. Phys. Rev. D , , doi:10.1103/PhysRevD.19.3524.34. Maluf, J.W. The teleparallel equivalent of general relativity. Ann. Phys. , , 339–357.35. Aldrovandi, R.; Pereira, J.G. Teleparallel Gravity: An Introduction ; Springer: Dordrecht,The Netherlands, 2013.36. Ferraro, R.; Fiorini, F. Modified teleparallel gravity: Inflation without an inflaton.
Phys. Rev. D , , doi:10.1103/PhysRevD.75.084031.37. Bengochea, G.R.; Ferraro, R. Dark torsion as the cosmic speed-up. Phys. Rev. D , ,124019. doi:10.1103/PhysRevD.79.124019.38. Linder, E.V. Einstein’s other Gravity and the Acceleration of the Universe. Phys. Rev. D , , doi:10.1103/PhysRevD.81.127301 . alaxies ,
39. Harko, T.; Lobo, F.S.N.; Otalora, G.; Saridakis, E.N. Non-minimal torsion-matter couplingextension of f ( T ) gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.89.124036.40. Harko, T.; Lobo, F.S.N.; Otalora, G.; Saridakis, E.N. f ( T, T ) gravity and cosmology. ArXiv E-Prints , , arXiv:1405.0519.41. Kiani, F.; Nozari, K. Energy conditions in F ( T, Θ) gravity and compatibility with a stable deSitter solution. Phys. Lett. B , , 554–561.42. Harko, T. Galactic rotation curves in modified gravity with non-minimal coupling between matterand geometry. Phys. Rev. D , , doi:10.1103/PhysRevD.81.08405043. Sotiriou, T.P.; Faraoni, V. Modified gravity with R -matter couplings and (non-)geodesic motion. Class. Quantum Gravity , , doi:10.1088/0264-9381/25/20/205002.44. Damour, T.; Polyakov, A.M. The String dilaton and a least coupling principle. Nucl. Phys. B , , 532–558.45. Damour, T.; Esposito-Farese, G. Tensor-multi-scalar theories of gravitation. Class. QuantumGravity , , doi:10.1088/0264-9381/9/9/015.46. Gottlober, S.; Schmidt, H.J.; Starobinsky, A.A. Sixth-order gravity and conformaltransformations. Class. Quantum Gravity , , doi:10.1088/0264-9381/7/5/018.47. Drummond, I.T.; Hathrell, S.J. QED Vacuum Polarization in a Background Gravitational Fieldand Its Effect on the Velocity of Photons. Phys. Rev. D , , doi:10.1103/PhysRevD.22.343.48. Buchdahl, H.A. Non-linear Lagrangians and cosmological theory. Mon. Not. R. Astron. Soc. , , 1–8.49. Barrow, J.D.; Ottewill, A.C. The stability of general relativistic cosmological theory. J. Phys. AMath. Gen. , , 2757–2776.50. Koivisto, T. A note on covariant conservation of energy-momentum in modified gravities. Class.Quantum Gravity , , doi:10.1088/0264-9381/23/12/N01.51. Teyssandier, P.; Tourrenc, Ph. The Cauchy problem for the R + R theories of gravity withouttorsion. J. Math. Phys. , , doi:10.1063/1.525659.52. Whitt, B. Fourth-order gravity as general relativity plus matter. Phys. Lett. B , , 176–178.53. Wands, D. Extended gravity theories and the Einstein–Hilbert action. Class. Quantum Gravity , , doi:10.1088/0264-9381/11/1/025.54. Faraoni, V. De Sitter space and the equivalence between f ( R ) and scalar-tensor gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.75.067302.55. Olmo, G.J. Limit to general relativity in f ( R ) theories of gravity. Phys. Rev. D , ,doi:10.1103/PhysRevD.75.023511.56. Harko, T.; Koivisto, T.S.; Lobo, F.S.N. Palatini formulation of modified gravity witha non-minimal curvature-matter coupling. Mod. Phys. Lett. A , , doi:10.1142/S0217732311035869.57. Mohseni, M. Non-geodesic motion in f(G) gravity with non-minimal coupling. Phys. Lett. B , , 89–92.58. Bertolami, O.; Lobo, F.S.N.; Páramos, J. Nonminimum coupling of perfect fluids to curvature. Phys. Rev. D , , doi:10.1103/PhysRevD.78.064036. alaxies ,
59. Faraoni, V. Lagrangian description of perfect fluids and modified gravity with an extra force.
Phys. Rev. D , , doi:10.1103/PhysRevD.80.124040.60. Bertolami, O.; Páramos, J. Homogeneous spherically symmetric bodies with a nonminimalcoupling between curvature and matter: The choice of the Lagrangian density for matter. ArXiv E-Prints , , arXiv:1306.1177.61. Minazzoli, O. Conservation laws in theories with universal gravity/matter coupling. Phys. Rev. D , , doi:10.1103/PhysRevD.88.027506.62. Bisabr, Y. Non-minimal Gravitational Coupling of Phantom and Big Rip Singularity. Gen. Relativ. Gravit. , , 1559–1566.63. Harko, T. Modified gravity with arbitrary coupling between matter and geometry. Phys. Lett. B , , 376–379.64. Wu, Y.-B.; Zhao, Y.-Y.; Lu, J.-W.; Zhang, X.; Zhang, C.-Y.; Qiao, J.-W. Five-dimensionalgeneralized f ( R ) gravity with curvature-matter coupling. Eur. Phys. J. C , , doi:10.1140/epjc/s10052-014-2791-9.65. Olmo, G.J.; Rubiera-Garcia, D. Brane-world and loop cosmology from a gravity-matter couplingperspective. ArXiv E-Prints , , arXiv:1405.7184.66. Harko, T. The matter Lagrangian and the energy-momentum tensor in modified gravity withnonminimal coupling between matter and geometry. Phys. Rev. D , , doi:10.1103/PhysRevD.81.044021.67. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Fields ; Butterworth-Heinemann: Oxford,UK, 1998.68. Fock, V.
The Theory of Space, Time and Gravitation ; Pergamon Press: London, UK, 1959.69. Faraoni, V. Viability criterion for modified gravity with an extra force.
Phys. Rev. D , ,doi:10.1103/PhysRevD.76.127501.70. Puetzfeld, D.; Obukhov, Y.N. Motion of test bodies in theories with nonminimal coupling. Phys.Rev. D , , doi:10.1103/PhysRevD.78.121501.71. Bertolami, O.; Martins, A. On the dynamics of perfect fluids in non-minimally coupled gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.85.024012.72. Bertolami, O.; Paramos, J. Do f ( R ) theories matter? Phys. Rev. D , 77, doi:10.1103/PhysRevD.77.084018.73. Bertolami, O.; Sequeira, M.C. Energy Conditions and Stability in f ( R ) theories of gravity withnon-minimal coupling to matter. Phys. Rev. D , , doi:10.1103/PhysRevD.79.104010.74. Wang, J.; Wu, Y.-B.; Guo, Y.-X.; Yang, W.-Q.; Wang, L. Energy conditions and stability ingeneralized f ( R ) gravity with arbitrary coupling between matter and geometry. Phys. Lett. B , , 133–138.75. Wang, J.; Wu, Y.-B.; Guo, Y.-X.; Qi, F.; Zhao, Y.-Y.; Sun, X.-Y. Conditions and instability in f ( R ) gravity with non-minimal coupling between matter and geometry. Eur. Phys. J. C , , 541–546.76. Wang, J.; Wu, Y.-B.; Guo, Y.-X.; Yang, W.-Q.; Wang, L. Energy conditions and stability ingeneralized f ( R ) gravity with arbitrary coupling between matter and geometry. Phys. Lett. B , , 133–138. alaxies ,
77. Sotiriou, T.P. The viability of theories with matter coupled to the Ricci scalar.
Phys. Lett. B , , 225–228.78. Bertolami, O.; Paramos, J. On the non-trivial gravitational coupling to matter. Class. QuantumGravity , , doi:10.1088/0264-9381/25/24/245017.79. Tamanini, N.; Koivisto, T.S. Consistency of non-minimally coupled f ( R ) gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.88.064052.80. Obukhov, Y.N.; Puetzfeld, D. Conservation laws in gravitational theories with generalnonminimal coupling. Phys. Rev. D , , doi:10.1103/PhysRevD.87.081502.81. Puetzfeld, D.; Obukhov, Y.N. Covariant equations of motion for test bodies in gravitationaltheories with general nonminimal coupling. Phys. Rev. D , , doi:10.1103/PhysRevD.87.044045.82. Puetzfeld, D.; Obukhov, Y.N. Equations of motion in gravity theories with nonminimal coupling:A loophole to detect torsion macroscopically? Phys. Rev. D , , doi:10.1103/PhysRevD.88.064025.83. Bertolami, O.; March, R.; Páramos, J. Solar System constraints to nonminimally coupled gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.88.064019.84. Castel-Branco, N.; Páramos, J.; March, R. Perturbation of the metric around a sphericalbody from a nonminimal coupling between matter and curvature. ArXiv E-Prints , ,arXiv:1403.7251.85. Garcia, N.M.; Lobo, F.S.N. Wormhole geometries supported by a nonminimal curvature-mattercoupling. Phys. Rev. D , , doi:10.1103/PhysRevD.82.104018.86. Garcia, N.M.; Lobo, F.S.N. Nonminimal curvature-matter coupled wormholes with mattersatisfying the null energy condition. Class. Quantum Gravity , , doi:10.1088/0264-9381/28/8/085018.87. Bertolami, O.; Ferreira, R.Z. Traversable wormholes and time machines in nonminimally coupledcurvature-matter f ( R ) theories. Phys. Rev. D , , doi:10.1103/PhysRevD.85.104050.88. Nesseris, S. Matter density perturbations in modified gravity models with arbitrary couplingbetween matter and geometry. Phys. Rev. D , , doi:10.1103/PhysRevD.79.044015.89. Bertolami, O.; Frazao, P.; Paramos, J. Accelerated expansion from a nonminimal gravitationalcoupling to matter. Phys. Rev. D , , doi:10.1103/PhysRevD.81.104046.90. Thakur, S.; Sen, A.A.; Seshadri, T.R. Non-minimally coupled f ( R ) cosmology. Phys. Lett. B , , 309–314.91. Bisabr, Y. Modified gravity with a nonminimal gravitational coupling to matter. Phys. Rev. D , , doi:10.1103/PhysRevD.86.044025.92. Bertolami, O.; Frazao, P.; Paramos, J. Reheating via a generalized non-minimal coupling ofcurvature to matter. Phys. Rev. D , , doi:10.1103/PhysRevD.83.044010.93. Bertolami, O.; Paramos, J. Modified Friedmann Equation from Nonminimally Coupled Theoriesof Gravity. Phys. Rev. D , , doi:10.1103/PhysRevD.89.044012.94. Bertolami, O.; Frazao, P.; Paramos, J. Cosmological perturbations in theories with non-minimalcoupling between curvature and matter. J. Cosmol. Astropart. Phys. , , doi:10.1088/1475-7516/2013/05/029. alaxies ,
95. Thakur, S.; Sen, A.A. Can structure formation distinguish Λ CDM from nonminimal f ( R ) gravity? Phys. Rev. D , , doi:10.1103/PhysRevD.88.044043.96. Bertolami, O.; Paramos, J. Mimicking dark matter through a non-minimal gravitational couplingwith matter. J. Cosmol. Astropart. Phys. , , doi:10.1088/1475-7516/2010/03/009.97. Bertolami, O.; Frazao, P.; Paramos, J. Mimicking dark matter in galaxy clustersthrough a nonminimal gravitational coupling with matter. Phys. Rev. D , ,doi:10.1103/PhysRevD.86.044034.98. Bertolami, O.; Paramos, J.; Harko, T.; Lobo, F.S.N. Non-minimal curvature-matter couplings inmodified gravity. ArXiv E-Prints , , arXiv:0811.2876.99. Bertolami, O.; Paramos, J. Minimal extension of General Relativity: Alternative gravity modelwith non-minimal coupling between matter and curvature. Int. J. Geom. Methods Mod. Phys. , , doi:10.1142/S0219887814600032.100. Harko, T.; Lobo, F.S.N. Geodesic deviation, Raychaudhuri equation, and tidal forces inmodified gravity with an arbitrary curvature-matter coupling. Phys. Rev. D , ,doi:10.1103/PhysRevD.86.124034.101. Lobo, F.S.N.; Harko, T. Extended f ( R, L m ) theories of gravity. ArXiv E-Prints , ,arXiv:1211.0426.102. Wang, J.; Liao, K. Energy conditions in f ( R, L m ) gravity. Class. Quantum Gravity , ,doi:10.1088/0264-9381/29/21/215016.103. Huang, R.-N. The Wheeler-DeWitt equation of f ( R, L m ) gravity in minisuperspace. ArXivE-Prints , , arXiv:1304.5309.104. Tian, D.W.; Booth, I. Lessons from f ( R, R c , R m , L m ) gravity: Smooth Gauss-Bonnetlimit, energy-momentum conservation and nonminimal coupling. ArXiv E-Prints , ,arXiv:1404.7823.105. Milgrom, M. A modification of the Newtonian dynamics as a possible alternative to the hiddenmass hypothesis. Astrophys. J. , , 365–370.106. Bekenstein, J.D. Relativistic gravitation theory for the modified Newtonian dynamics paradigm. Phys. Rev. D , , doi.org/10.1103/PhysRevD.70.083509.107. Hawking, S.; Ellis, G.F.R. The Large Scale Structure of Space-Time ; Cambridge University Press:Cambridge, UK, 1973.108. Mashhoon, B. Tidal Gravitational Radiation.
Astrophys. J. , , 83–86.109. Mashhoon, B. On tidal phenomena in a strong gravitational field. Astrophys. J. , ,705–716.110. Mashhoon, B.; Theiss, D.S. Relativistic Tidal Forces and the Possibility of Measuring Them. Phys. Rev. Lett. , , doi:10.1103/PhysRevLett.49.1542.111. Ohanian, H.C. Gravitation and Spacetime ; Norton: New York, NY, USA, 1976.112. Gasperini, M.; Piazza, F.; Veneziano, G. Quintessence as a runaway dilaton.
Phys. Rev. D , , doi:10.1103/PhysRevD.65.023508.113. Damour, T.; Vilenkin, A. String theory and inflation. Phys. Rev. D , ,doi:10.1103/PhysRevD.53.2981. alaxies , Phys. Rev. D , , doi:10.1103/PhysRevD.66.064008.115. Damour, T.; Polyakov, A.M. The string dilation and a least coupling principle. Nucl. Phys. B , , 532–558.116. Damour, T.; Polyakov, A.M. String theory and gravity. Gen. Relativ. Gravit. , ,1171–1176.117. Damour, T.; Piazza, F.; Veneziano, G. Violations of the equivalence principle in a dilaton-runawayscenario. Phys. Rev. D , , doi:10.1103/PhysRevD.66.046007.118. Damour, T.; Donoghue, J.F. Equivalence principle violations and couplings of a light dilaton. Phys. Rev. D , , doi:10.1103/PhysRevD.82.084033.119. Damour, T. Theoretical aspects of the equivalence principle. Class. Quantum Gravity , ,doi:10.1088/0264-9381/29/18/184001.120. Damour, T.; Gibbons, G.W.; Gundlach, C.C. Dark matter, time-varying G , and a dilaton field. Phys. Rev. Lett. , , doi:10.1103/PhysRevLett.64.123.121. Casas, J.A.; García-Bellido, J.; Quiros, M. Scalar-tensor theories of gravity with Phi-dependentmasses. Class. Quantum Gravity , , doi:10.1088/0264-9381/9/5/018.122. Wetterich, C. An asymptotically vanishing time-dependent cosmological constant. Astron.Astrophys. , , 321–328.123. Amendola, L. Coupled quintessence. Phys. Rev. D , , doi:10.1103/PhysRevD.62.043511.124. Das, S.; Banerjee, N. Brans-Dicke scalar field as a chameleon. Phys. Rev. D , ,doi:10.1103/PhysRevD.78.043512.125. Minazzoli, O. The γ parameter in Brans-Dicke-like (light-)Scalar-Tensor theory with a universalscalar/matter coupling and a new decoupling scenario. ArXiv E-Prints , , arXiv:1208.2372.126. Aviles, A.; Cervantes-Cota, J.L. Dark degeneracy and interacting cosmic components. Phys. Rev. D , , doi.org/10.1103/PhysRevD.84.083515.127. Overduin, J.M.; Wesson, P.S. Kaluza-Klein gravity. Phys. Rep. , , 303–378.Fujii, Y.; Maeda, K.I. The Scalar-Tensor Theory of Gravitation ; Cambridge University Press:Cambridge, UK, 2003.128. Minazzoli, O.; Harko, T. New derivation of the Lagrangian of a perfect fluid with a barotropicequation of state.
Phys. Rev. D , , doi:10.1103/PhysRevD.86.087502.129. Harko, T.; Lobo, F.S.N. Generalized dark gravity. Int. J. Mod. Phys. D , , doi:10.1142/S0218271812420199.130. Poplawski, N.J. A Lagrangian description of interacting dark energy. ArXiv E-Prints , ,arXiv:gr-qc/0608031.131. Houndjo, M.J.S. Reconstruction of f ( R, T ) gravity describing matter dominated and acceleratedphases. Int. J. Mod. Phys. D , , doi:10.1142/S0218271812500034.132. Jamil, M.; Momeni, D.; Raza, M.; Myrzakulov, R. Reconstruction of some cosmological modelsin f ( R, T ) cosmology. Eur. Phys. J. C , , doi:10.1140/epjc/s10052-012-1999-9.133. Houndjo, M.J.S.; Piattella, O.F. Reconstructing f ( R, T ) gravity from holographic dark energy. Int. J. Mod. Phys. D , , doi:10.1142/S0218271812500241. alaxies , f ( R, T ) gravity and the effect of conformal anomaly. Can. J. Phys. , , 548–553.135. Sharif, M.; Zubair, M. Thermodynamics in f ( R, T ) Theory of Gravity.
J. Cosmol. Astropart.Phys. , , doi:10.1088/1475-7516/2012/03/028.136. Sharif, M.; Rani, S.; Myrzakulov, R. Analysis of f ( R, T ) gravity models through energyconditions. Eur. Phys. J. Plus , , doi:10.1140/epjp/i2013-13123-0.137. Sharif, M.; Zubair, M. Energy Conditions Constraints and Stability of Power Law Solutions in f ( R, T ) Gravity.
J. Phys. Soc. Jpn. , , doi:10.7566/JPSJ.82.014002.138. Alvarenga, F.G.; Houndjo, M.J.S.; Monwanou, A.V.; Orou, J.B.C. Testing some f ( R, T ) gravitymodels from energy conditions. J. Mod. Phys. , , 130–139.139. Santos, A.F. Gödel solution in f ( R, T ) gravity. Mod. Phys. Lett. A , , doi:10.1142/S0217732313501411.140. Priyanka, S.R. Some Kaluza-Klein cosmological models in f ( R, T ) gravity theory. Astrophys. Space Sci. , , 389–397.141. Naidu, R.L.; Reddy, D.R.K.; Ramprasad, T.; Ramana, K.V. Bianchi type-V bulk viscous stringcosmological model in f ( R, T ) gravity. Astrophys. Space Sci. , , 247–252.142. Sharif, M.; Zubair, M. Study of Bianchi I anisotropic model in f ( R, T ) gravity. Astrophys. Space Sci. , , 457–465.143. Chakraborty, S. An alternative f ( R, T ) gravity theory and the dark energy problem. Gen. Relativ. Gravit. , , 2039–2052.144. Alvarenga, F.G.; de la Cruz-Dombriz, A.; Houndjo, M.J.S.; Rodrigues, M.E.; Sáez-Gómez, D.Dynamics of scalar perturbations in f ( R, T ) gravity. Phys. Rev. D , , doi.org/10.1103/PhysRevD.87.103526.145. Sharif, M.; Zubair, M. Reconstruction and stability of f ( R, T ) gravity with Ricci and modifiedRicci dark energy. Astrophys. Space Sci. , , 529–537.146. Rubin, V.C.; Ford, W.K.; Thonnard, N. Rotational properties of 21 SC galaxies with a largerange of luminosities and radii, from NGC 4605 ( R = 4 kpc) to UGC 2885 ( R = 122 kpc). Astrophys. J. , , 471–487.147. Persic, M.; Salucci, P.; Stel, F. The universal rotation curve of spiral galaxies—I. The dark matterconnection. Mon. Not. R. Astron. Soc. , , 27–47.148. Borriello, A.; Salucci, P.; Borriello, A.; Salucci, P. The dark matter distribution in disc galaxies. Mon. Not. R. Astron. Soc. , , 285–292.149. Binney, J.; Tremaine, S. Galactic Dynamics ; Princeton, N.J., Ed.; Princeton University Press:Princeton, NJ, USA, 2008.150. Overduin, J.M.; Wesson, P.S. Dark matter and background light.
Phys. Rep. , , 267–406.151. Boehmer, C.G.; Harko, T.; Lobo, F.S.N. Dark matter as a geometric effect in f ( R ) gravity. Astropart. Phys. , , 386–392.152. Boehmer, C.G.; Harko, T.; Lobo, F.S.N. Generalized virial theorem in f ( R ) gravity. J. Cosmol.Astropart. Phys. , , doi:10.1088/1475-7516/2008/03/024. alaxies , J. Cosmol. Astropart. Phys. , ,doi:10.1088/1475-7516/2013/07/024.154. Capozziello, S.; Harko, T.; Koivisto, T.S.; Lobo, F.S.N.; Olmo, G.J. Galactic rotation curves inhybrid metric-Palatini gravity. Astropart. Phys. , , 65–75.155. Harko, T.; Koivisto, T.S.; Lobo, F.S.N.; Olmo, G.J. Metric-Palatini gravity unifying localconstraints and late-time cosmic acceleration. Phys. Rev. D , , doi:10.1103/PhysRevD.85.084016.156. Capozziello, S.; Harko, T.; Koivisto, T.S.; Lobo, F.S.N.; Olmo, G.J. Cosmology of hybridmetric-Palatini f ( X ) -gravity. J. Cosmol. Astropart. Phys. , , doi:10.1088/1475-7516/2013/04/011.157. Capozziello, S.; Harko, T.; Lobo, F.S.N.; Olmo, G.J. Hybrid modified gravity unifying local tests,galactic dynamics and late-time cosmic acceleration. Int. J. Mod. Phys. D , , doi:10.1142/S0218271813420066.158. Capozziello, S.; Harko, T.; Lobo, F.S.N.; Olmo, G.J.; Vignolo, S. The Cauchy problem in hybridmetric-Palatini f ( X ) -gravity. Int. J. Geom. Methods Mod. Phys. , , doi:10.1142/S021988781450042X.159. Nucamendi, U.; Salgado, M.; Sudarsky, D. Alternative approach to the galactic dark matterproblem. Phys. Rev. D , , doi:10.1103/PhysRevD.63.125016.160. Lake, K. Galactic Potentials. Phys. Rev. Lett. , , doi:10.1103/PhysRevLett.92.051101.161. Salucci, P.; Persic, M. Dark Halos around Galaxies, in Dark and Visible Matter in Galaxies ;Persic, M., Salucci, P., Eds.; ASP Conference Series; Astronomical Society of the Pacific:San Francisco, CA, USA, 1997; Volume 117, p. 1.162. Pun, C.S.J.; Kovacs, Z.; Harko, T. Thin accretion disks in f ( R ) modified gravity models. Phys. Rev. D , , doi:10.1103/PhysRevD.78.024043.163. Mamon, G.A.; Lokas, E.L. Dark matter in elliptical galaxies—II. Estimating the mass within thevirial radius. Mon. Not. R. Astron. Soc. , , 705–722.c (cid:13)(cid:13)