$\kappa(R,T)$ gravity

In this note we explore a modified theory of gravitation that is not based on the least action principle, but on a natural generalization of the original Einstein's field equations. This approach leads to the non-covariant conservation of the stress-energy tensor, a feature shared with other Lagrangian theories of gravity such as the $f(R,T)$ case. We consider the cosmological implications of a pair of particular models within this theory, and we show that they have some interesting properties. In particular, for some of the studied models we find that the density is bounded from above, and cannot exceed a maximum value that depends on certain physical constants. In the last part of the work we compare the theory to the $f(R,T)$ case and show that they lead to different predictions for the motion of test particles.


I. INTRODUCTION
The least action principle has become one of the most powerful tools to build a physical theory and also their possible generalizations. Among the numerous advantages of the Lagrangian formalism we list the direct implementation of symmetries and the derivation of general conservation laws. Nevertheless, there is no reason to believe that ordinary symmetries and/or standard conservation laws will always hold in a final theory of Nature. In this sense, it is interesting to recall that Einstein did not originally followed a variational principle in the derivation of General Relativity (GR) [1][2][3]. Instead, he arrived to the correct field equations following a very different approach, one that succeeded with the eventual addition of a trace term directly in the field equations. Indeed, it is well known that the equivalence principle and general covariance were the foundational concepts of the theory, and the variational principle, i.e, the Einstein-Hilbert Action (EHA), was discovered and incorporated to the theory when the correct field equations had already been derived. [4] Similarly, the other classical field theory, namely, Maxwell's Eletrodynamics (ME), was only completed after the addition of a source term (the Maxwell displacement current), and it wasn't either originally conceived from any variational principle. [5] In light of these historical facts, there is no reason to reject the search for an alternative approach,one different from the Lagrangian formalism. In this sense, a more general theory could be formulated following an alternative but consistent line of reasoning, and the variational principle (or something alike) could be incorporated in the last stages of the completion of the theory in order to strengthen and to reinforce the formalism. The fact that both the GR and ME field equations were found without resorting to a variational principle sends us the message that maybe a different approach deserves to being taken * gines.landau@gmail.com into account. Regarding the case of GR and their extensions, there exists a high degree of arbitrariness in the choice of the specific generalized gravity Lagrangian. One of the most natural strategies is to replace the curvature scalar R, by a function f (R) in the action [6][7][8][9][10][11][12]. This approach leads to field equations that give rise to a rich phenomenology (both in the metric and metric-affine or Palatini formalisms), although some of the beauty and simplicity of the original theory are lost in the process. Regardless the specific choice among the innumerable possible generalized Lagrangians, the resulting field equations are much more complicated than those of GR. In view of the above, given the arbitrariness in the choice of the possible gravity Lagrangian, we do not begin from a Lagrangian formalism to look for a modification of GR. Rather, we follow Maxwell's and Einstein's original approaches of adding new possible source terms directly in the field equations, this means that we focus on a similar strategy that proved successful in the completion of the classical Electromagnetic theory, and in the first correct derivation of GR. In particular, the theory considered here is based on a quite natural extension of GR, where the modification of the field equations involves the addition of terms that only include the curvature scalar and the trace of the stress-energy tensor. Therefore, in vacuum the field equations boil down to those of GR, but in presence of matter there can be significant departures, in particular the stress-energy tensor will not be covariantly conserved in the general case. There are several examples in the literature of gravitational theories of this kind. An example is Rastall's gravitational theory [13,14], which is also non-conservative since the divergence of T µν does not vanish in general. Another more recent example is given by the so-called f (R, T ) modified theories of gravity [15,[19][20][21][22]. We should mention that this note presents an study of a modified theory of gravity at a preliminary level, it is a first step in a different direction than that adopted by the usual modified gravity theories, and further investigations are required to discuss in more detail some important aspects that are addressed in this work. It is becoming increasingly clear that quantum field theories without a traditional Lagrangian description are important and even populate much of the QFT landscape. They also offer new opportunities in the search of new type of 4-manifold invariants [23,24]. This fact represents another good motivation to explore here an example of a Non-Lagrangian modified gravitational theory.

II. DEFINITIONS AND FIELD EQUATIONS
Our framework is based on the following field equations Where R µν is the Ricci tensor, g µν is the space-time metric, Λ is a cosmological constant, T µν the stress-energy tensor of the matter sources, and κ(R, T ) corresponds to the Einstein gravitational constant that we are promoting to the status of a function of the traces T ≡ g µν T µν , and R ≡ g µν R µν . The possible dependence of the gravitational constant κ on scalars means that we explore the possibility of a running gravitational constant, i.e. we generalize the original Einstein's gravitational constant, but not at the level of an action functional. A varying gravitational constant in the action leads to a Brans-Dicke type theory [25][26][27], with quite different field equations from (1). The field equations (1) imply the non-covariant conservation of T µν . Indeed, since the left hand side of these equations is divergence-free, we have Then, the non-conservation of the T µν can be expressed as In what follows, some cosmological implications (homogeneous and isotropic universe for a perfect fluid) of two particular cases are analyzed. The first model considered arises by setting, κ(T ) = 8πG − λT , and corresponds to a matter-matter coupling. The second model that will be studied is characterized by a gravitational "constant" that varies as κ (R) = 8πG + αR, which will provide a coupling between matter and curvature terms. We assume that the coupling constants λ, α are sufficiently small to be consistent with a small violation of the covariant conservation of the stress-energy tensor. Obviously, in the limit λ, α → 0, Einstein's GR is recovered. If we consider an homogeneous and isotropic universe filled by a perfect fluid as the matter source, the stressenergy tensor will be T µν = (p + ρ)u µ u ν − pg µν (4) where p, ρ, and u µ are the pressure, the density and the macroscopic speed of the medium, respectively. On the other hand, the standard FLRW metric for modeling the assumed properties of such a universe leads to the line element With these ingredients, the two independent Modified Friedmann Equations (MFE) for a general κ(R, T ) model are Two independent models of the form κ(T ) and κ(R) are analyzed in the next subsections.

B. Matter-matter coupling
Here we analyze the cosmological implications of the model, κ(T ) = k − λT , where k ≡ 8πG (c = 1) and λ is a constant with the appropriate units. The reasons to choose a negative sign will be understood later on. Therefore, with such a choice for κ(R, T ) the field equations 1 acquire the form The MFE that arise by solving the field equations are (10) Where, H =ȧ/a. In order to better understand the physical meaning of the MFE, and for the sake of simplicity it is convenient to make use of a equation of state of the type, p = wρ. Then, for the flat case (K −1 = 0), we can write It is therefore clear that the new additional term introduced vanishes for w = 1/3 (radiation-dominated universe). For w = −1/3 it vanishes for the second equation as well. For a radiation-dominated universe, the traceless of the stress-energy tensor implies that the MFE collapse to the standard solution of GR with a cosmological constant. On the other hand, the contribution of the new term is proportional to ρ 2 , for −1/3 < w < 1/3 this quadratic term will be positive in the acceleration equation, and this means that at sufficiently high densities it could contribute to the cosmic speed-up. However, to explain the late-time acceleration (low densities) we unavoidably need the inclusion of a cosmological constant.
The non-conservation of the stress-energy tensor implies a modification of the relativistic fluid equation characterized byρ The correction with respect to the GR case is represented by the presence of a certain function F (ρ) which is explicitly given by Where we have denoted a certain constant ρ m with units of density as In Fig.1, We plot the behavior of the correction factor F as function of ρ/ρ m . Notice that when ρ << ρ m , then, F (ρ) 1 and we essentially recover the same fluid equation of GR. However, for ultra high-densities, namely, ρ >> ρ m , F (ρ) 2. For such a regime, the relativistic fluid equation becomeṡ This implies that, ρ(t) ∼ a(t) −6(1+w) which suggests a much more rapid decrease of the density with the scale factor. However, the density cannot be arbitrarily large, and in general it will be of order ρ m . This is due to the requirement H 2 ≥ 0 which imposes an upper bound for the density. Indeed, suppose that at sufficiently high densities we can neglect the contribution of the cosmological constant term compared to the other two in the first MFE. Then, reality of H 2 requires which implies For w = 1/3. Notice that for ρ = ρ max , H ∼ √ Λ and therefore the Hubble parameter behaves such as de Sitter vacuum solution. Finally, the choice of the negative sign in the function κ(T ) = 8πG−λT is due to very good reasons that turn out to be evident now. If we had chosen the positive sign, we would have that at high densities H 2 ∼ ρ 2 , which is worst in terms of divergences than the GR case. [29] Field equations similar to (1), or (11)(12) in cosmology, could arise from Ricci and specially from generalized Ricci-Gauss-Bonnet holographic dark energy models, see for example Refs. [16][17][18] 1. Dust solution Let us solve the equations for the simplest cosmological model, which is obtained by assuming a dust universe, where p = 0. The MFE for such a universe become Using thatä/a =Ḣ + H 2 and combining both equations we obtainḢ These results mean that H(t) for dust evolves as in the GR case, and does not depend on the value of λ. For a matter-dominated universe where we can neglect the contribution of the cosmological constant, we have H(t) = 2/3t, a(t) = t 2/3 . Then, the dust solution shows that the λ parameter for this particular model has no effect in the expansion rate of the universe, which turns out to be identical to the GR case, and the main difference is the existence of a bound for the density, given by ρ max = 8πG/3λ. On the other hand, the dust solution for the f (R, T ) theory was studied in [15], where for a function of the form f (T ) = λT the Hubble parameter was found to be dependent on λ as Therefore, H(t) f > H(t) κ , i.e, the Hubble parameter in the f (R, T ) theory turns out to be a bit larger than the value of the Hubble parameter in κ(R, T ) for a similar functional dependence on T . Nevertheless, the exact discrepancy with respect to the GR case is not clear since the authors do not constrain the value of λ.
On the other hand, to see how the density evolves with time for a dust universe in the κ(R, T ) theory, we come back to the fluid equation (13) which for dust (w = 0) reduce to the expressioṅ Where we have employed the identity ρ m = 3ρ max . To solve this equation, first notice that Where the minimum value 1/2 of this function is reached for ρ = ρ max and the maximum value is approached for ρ << ρ max . Then, an acceptable solution can be provided if we set, N where N is a number be- Assuming that the density obeys a relation of the type ρ = Ct α , where C is a constant, and substituting this in Eq. (23), we obtain the condition α + 2N = 0, which implies At low densities N 1, and for such a regime the density evolves with time as ρ Ct −2 , like the standard dust solution of GR.

Stationary solutions and exponential expansion
By setting H = H 0 = const. into Eq. (11) and rearranging terms we find an algebraic quadratic equation for the density given by The solutions of this equation are Where ρ max is the maximum density deduced in the previous subsection. Therefore, for H 0 = Λ/3 we have ρ = 0 which is the standard vacuum de Sitter solution, and ρ = ρ max . Moreover, the difference, 3H 2 0 − Λ is bounded from above. Indeed, we can rewrite the last equation in the form Therefore, Another interesting feature of this model is the prediction of a specific value of the density (depending on λ) for the exponential expansion governed by the cosmological constant, which will be exactly equal to ρ max . Indeed, the acceleration equation (12) reduces toä/a = Λ/3 (which implies a(t) ∼ exp( Λ/3t)) for ρ inf = 8πG/λ(1 − 3w) = ρ max . Then, the inflation takes place when the density reaches the maximum value.

C. Matter-curvature coupling
Here we study cosmological solutions for the theory κ(R) = 8πG + αR, where R is the curvature scalar and α a constant with units of the inverse of the density. The field equations with such a choice for κ(R, T ) are With a bit of algebra, we find that the expansion rate equation for this model is given by Where H ≡ȧ/a, and f (ρ) is a function of the density given by The constant ρ 0 is a certain density denoted by, ρ 0 = 1/α. Since we expect α to be very small, ρ 0 will presumably take a very large value. For ρ << ρ 0 , f (ρ) ≈ 0 and we recover the result of GR for the expansion rate equation. At ρ = ρ 0 , the function f (ρ) is always regular for all the physically acceptable values of w, and H 2 is also regular as required. Indeed, we have that f (ρ 0 ) = 1/(2 − 3w) which is singular only if w = 2/3. The expansion rate parameter at ρ = ρ 0 becomes It is important to note that this model also implies that in general there should exist bounds for the density, namely, the density cannot take an infinite value (with the exception of some particular case). Indeed, for ρ >> ρ 0 we have that f (ρ) ≈ 1/(1 − 3w). Therefore the first contribution due to the density in Eq.(31) vanishes identically, and the Hubble parameter becomes for such a regime which is always negative for w < 1/3. (The special case w = 1/3 will be analyzed apart). Therefore, the limit ρ >> ρ 0 lacks physical sense, and we conclude that the density cannot be much larger than ρ 0 = 1/α when w < 1/3. On the other hand, the case w = 1/3 corresponds to a radiation-dominated universe, the function f for such a value of the parameter w becomes, f (ρ) = ρ/ρ 0 . Taking this into account and setting w = 1/3 into Eq.(31), we have In this case, there are no bounds for the growth of the density because the assumed smallness of Λ/ρ 0 compared to 2πG assures that the quadratic density term will not acquire a negative sign.
Regarding the generalized Friedmann acceleration equation for this theory, we havë Where the scale factor H 2 is provided by Eq. (31). The acceleration equation is difficult to interpret given the quantity of terms involved. For ρ << ρ 0 we obviously recover the GR result, namely, the cosmological constant rules the late-time cosmic speed-up. It is worth noting that the term that accounts for the accelerated expansion vanishes when ρ = ρ 0 /(3w + 1). It turns out that for such a specific value of the density the effects of the accelerated expansion are null. Let us conclude with some brief comments regarding the theory κ(R) = 8πG−αR. By solving again the equations for the FLRW metric and the perfect fluid, we find a similar expression for the modified expansion rate equation which is explicitly given by (37) Where the function g(ρ) is denoted by The difference with respect to the theory κ(T ) = κ + αR discussed before, lies in the fact that the latter is better in terms of singularities. In particular, the main difference is due to the behavior of the functions f (ρ) and g(ρ) regarding regular properties. In contrast to f (ρ), that was regular everywhere and in particular at ρ = ρ 0 /(1 − 3w), the function g(ρ) is singular at such points and therefore H 2 for the model κ(R) = κ − αR will also diverge for that family of values of the density. Talking about the divergences of f (ρ) and g(ρ) is equivalent to talk about the divergences of the curvature scalar R for these models (see the appendix for details).

D. Static spherically symmetric perfect fluid
Spherically symmetric scenarios in GR and their extensions are very important because a wide range of phenomena such as black holes, neutron/quark stars and gravitational collapse can be theoretically described by means of spherical symmetry. The aim of this subsection is to present the main equations for the model κ(T ) = 8πG − λT , and to show how the a unknown metric components can be expressed in terms of the modified sources. Nonetheless, the explicit solution of the equations is an open problem that we leave for future works. The line element for a static spherically symmetric spacetime takes the form By solving the field equations of the theory (8) for Λ = 0, we obtain the system Where we already denoted the constant ρ m as ρ m ≡ 8πG/λ. The non-conservation of the stress-energy tensor implies another first order differential equation given by Where p ≡ dp/dr. The presence of a non-null right hand side represents the departure with respect to the GR case. On the other hand, the integration of the first equation gives Where Represents the "effective mass" content of the distribution within the sphere of radius r. The remaining unknown metric component, the function A(r), can also be expressed as an integral over the matter content and their density and pressure. Indeed, combining the first and second equations, we obtain

This equation can be immediately integrated to give
(46) With C an arbitrary integration constant. Therefore, the unknown metric components have been expressed in terms of the modified sources, as required. It is important to recall that the internal Schwarzschild-type solution should match with the external one at r = R, being R the radius of the compact object. Since the external Schwarzschild solution satisfies A(r) = B(r), and B(r) = 1 − 2GM/r, we see that the obtained internal metric components have the appropriate structure to match with the external solution, and this fixes the value of the constant C to be

E. Comparison among κ(R, T ) and f (R, T ) theories
It would be interesting to compare the theory κ(R, T ) to a gravitational Lagrangian theory, for example, we can compare the special case κ(R, T ) = κ(T ) to a subclass of the variational theories f (R, T ), which are also non-conservative theories (in the sense that ∇ ν T µν = 0). One of the most natural choices is a model of the type f (R, T ) = R + f 2 (T ). For a perfect fluid, this model leads to the field equations [15] Where T ef f µν for this particular theory f (R, T ) = R + f 2 (T ) is given by Then, the comparison to the field equations leads to the conclusion that the field equations of both theories only match if we assume in the field equations of the κ(T ) theory 1 the possibility of a variable cosmological 1 Note that if we consider the variable cosmological term in the action we will have a different theory term Λ(T ) that depends on the matter sources, namely It has been pointed out that recent cosmological data favor a variable cosmological constant [28].

Modified geodesic equation of motion
The fact that the stress-energy tensor is not covariantly conserved has the effect of modifying the equations of motion of particles. In fact, as a direct consequence an extra force will arise in the geodesic equation. To see in detail how this happens, we can write rewrite Eq.(3) for the perfect fluid in the form Where ∇ ν g µν = 0. Let us introduce now an auxiliary metric h µλ defined by h µλ ≡ g µλ − u µ u λ . Then, multiplying the last equation by h µλ we have Where we have used the identities, h µλ T µν = −h ν λ p, and h µλ u µ = 0. Therefore, with the aid of the identity where Γ µ νλ is the Levi-Civita connection of g µν , the modified geodesic equation of motion (53) acquires the form Where f µ denotes a four-vector "force" given by Therefore, for dust (p ≈ 0), we recover the geodesic equation of GR. Moreover, when κ(R, T ) = const = 8πG, the standard result of GR for perfect fluids with pressure is recovered as well. Notice that the vector f µ is orthogonal to u µ , namely, f µ u µ = 0. Eq.(56) suggests to redefine the pressure and density as, ρ ef f = κ(R, T )ρ, and p ef f = κ(R, T )p. Doing this, the form of the geodesic equation is identical to the GR case with the new variables ρ ef f and p ef f playing the role of p and ρ. On the other hand, it is worth comparing this extra force with other external force that arises in the geodesic equation of motion of the Lagrangian theory f (R, T ). For this theory, the extra force is given by [15]: Where f T (R, T ) = ∂f (R, T )/∂T . In order to see explicitly the differences among the extra force in both theories, we can select the same dependence on the trace T to compare predictions. For example, setting κ(R, T ) = κ(T ) = 8πG − λT , and f (R, T ) = f 1 (R) − βT , where λ, β are constants and f 1 (R) is an arbitrary function of R. Then, a first important consequence is that for T = 0 (photons), both theories predict the same extra force, but when T = 0 such is the case of massive particles, the two forces are different. Therefore, a detailed investigation on the trajectory of massive particles in a gravitational field could help to find out which among these different theories represent the most viable generalization of Einstein's GR.

F. Generalized energy conditions
Generalized energy conditions (GEC) in Extended Theories of Gravity have been studied in detail in several works, see for example [30][31][32][33]. To study the role of the energy conditions in the κ(R, T ) theory, with the aim to investigate if they are violated or not, it is convenient to recast the field equations in the form: Where the effective energy momentum tensor T ef f µν is defined by The effective energy momentum tensor T ef f µν , in turn allows to define the effective pressure p ef f and density ρ ef f necessary to present the conditions required for realizing each type of the energy conditions. In fact, p ef f and ρ ef f will be very similar to those that arose in the previous subsection. Indeed, by assuming that the content of the universe behaves like a perfect fluid, and for a flat FRLW metric we have It is convenient to focus on the case κ(R, T ) = κ(T ) in order to simplify the analysis. This choice for κ(R, T ) allows one to define the effective pressure and density as ρ ef f = κ(T )ρ + Λ, p ef f = κ(T )p − Λ. By using these expressions for ρ ef f and p ef f , we get the null energy condition (NEC), the weak energy condition (WEC), the strong energy condition (SNC) and the dominant energy condition (DEC) [30] as: NEC: WEC: SEC: Then, for a general κ(T ) model, the GEC will acquire the explicit expressions: WEC: It is not difficult to prove that for the perfect fluid and the model studied in this work, namely κ(T ) = 8πG−λT , the GEC are satisfied, since the density is bounded. Indeed, the NEC will be NEC: Where we have neglected the contribution of the cosmological constant. For w > −1, the NEC is automatically fulfilled if ρ ≤ ρ max , where ρ max = 8πG/λ(1 − 3w) in agreement with Eq. (18). Regarding the WEK, it is also satisfied identically if ρ ≤ ρ max . As for the SEC, we have: For w > −1/3, ρ ≤ ρ max the SEC is realized as well. Finally, the DEC acquires the form (72) For w < 1 the DEC is also satisfied for ρ ≤ ρ max . Therefore, the presence of a bound for the density guarantees that all the energy conditions are satisfied in a consistent way for the model κ(T ) = 8πG − λT .

III. SUMMARY AND CONCLUSIONS
The important degree of arbitrariness inherent in the choice of the gravity Lagrangian has lead to a large amount of different modified gravity proposals, many of which are so similar that it is difficult to distinguish one from the other. The Lagrangian formalism has undoubted advantages at the level of symmetries implementation and conservation-laws derivation, but possible theoretical alternatives to standard Lagrangian theories also deserve consideration. In this sense, the importance of Non-Lagrangian theories in other branches of theoretical physics such as quantum field theory is being acknowledged in the last years. Among their advantages, it seems increasingly clear that these theories offer new opportunities in the search of new types of invariants. In this work, and in absence of a foundational principle, we have explored an example of a Non-Lagrangian modified gravity theory inspired by Maxwell's approach to Electrodynamics, adding new possible source terms directly in the field equations, namely, we have investigated a gravitational analogue of the Maxwellian "displacement current" contribution. It should be noted that our approach does not mean that a variational formulation of the theory could not exist, but in this work we did not focus on that problem. In particular, we have analyzed some special cases that belong to the classification: κ(R, T ) = k + f (T ) and κ(R, T ) = k + f (R), which corresponds to matter-matter and matter-curvature couplings respectively. We carried out a preliminary study of some cosmological aspects of these models in a FLRW universe filled by a perfect fluid, and it was shown that the density in bounded from above in some of them. Furthermore, the formal similarities and differences among the theory κ(R, T ) and the Lagrangian theory f (R, T ) were also investigated. The field equations can match in some particular cases that imply a variable cosmological term that depends on the energy-matter content. However, both theories are essentially different at the level of the equations of motion for massive particles. Moreover, the generalized energy conditions were also investigated for the theory κ(T ) = 8πG − λT and we have shown that the existence of a maximum density ρ max guarantees that all the energy conditions are satisfied in a consistent way. In summary, we have presented in this work an example of a Non-Lagrangian modified gravity theory, which is a relatively unexplored research avenue in the field of modified gravity.

ACKNOWLEDGMENTS
It is a pleasure to thank Prof. S. Odintsov for pointing me out Ref. [12] and Prof. H. Shabani for bringing to my attention Refs. [20][21][22]. I thank Prof. G. J. Olmo and Prof. F. Rahaman for stimulating discussions and encouragement. I also thank the anonymous referees for comments that helped me to improve the quality of the manuscript.

IV. APPENDIX. DIVERGENCES AND ZEROS OF THE CURVATURE SCALAR
By contracting the field equations (1) with g µν it is easy to obtain a generic relation among the traces given by Then, the exact relation among R and T requires that we fix a particular κ(R, T ) model. Choosing a model of the type κ(T ) = 8πG − λT , we find that there exists an algebraic quadratic equation among R and T given by In the limit Λ → 0 we obtain R λT T − 8πG λ (75) Therefore R vanishes for T = 0 (the vacuum solution as in GR), and for T = 8πG/λ. In general, R vanishes for specific values of the pressure and the density that are solutions of the quadratic equation, λT 2 − 8πGT − 4Λ = 0. Given that ρ m ≡ 8πG/λ, T = ρ − 3p for the perfect fluid, and using a barotropic equation of state of the type p = wρ, we obtain such specific values of the density where R = 0 For w = 1/3. If ρ m >> 4 Λ/λ, we can approximate the solutions as Therefore, we see that the values of the density that vanish R are purely mathematical an not physical solutions. Recall that ρ is bounded and verifies, ρ ≤ ρ max , namely, ρ ≤ ρ m /(1 − 3w) according to Eq.(17) In the limit Λ/λ → 0, we obtain that ρ 2 = 0 (as in GR), and ρ 1 = ρ m /(1 − 3w) = ρ max , which is the extra solution with respect to the GR case. On the other hand, for the pair of models κ(R) = k ± αR, the exact relation among the traces is This relation becomes, with the assumed approximations, Where, ρ 0 ≡ 1/α. Then, R has a zero at a density given by ρ = Λ/2πG(3w − 1) which is negative for w < 1/3. Regarding the divergences, the curvature scalar is regular everywhere for the theory κ(R) = k + αR (with the exception of the special case w = 2/3, which contains a pole at ρ = ρ 0 ). On the contrary, the model κ(R) = k − αR yields to divergences for the family of values of the density given by ρ = ρ 0 /(1 − 3w)