## theory with non-minimal coupling

## Cosmologies in f ( R, L m ) between geometry and matter

Sergio Bravo Medina ∗

Departamento de F´ ısica, Pontificia Universidad Javeriana, Cra.7 No.40-62, Bogot´, a Colombia

Marek Nowakowski †

ICTP-South American Institute for Fundamental Research, Rua Dr. Bento Teobaldo Ferraz 271, 01140-070 S˜o Paulo, SP Brazil a

Ronaldo V. Lobato ‡

ICRANet, Piazza della Repubblica 10, Pescara, 65122, Italy

Davide Batic §

Department of Mathematics, Khalifa University of Science and Technology, Sas Al Nakhl Campus, P.O. Box 2533 Abu Dhabi, United Arab Emirates (Dated: August 2, 2024)

## Abstract

Among the recent extensions to standard General Relativity, f ( R, L m ) gravity has risen an interest given the possibility of coupling between geometry and matter. We examine the simplest model with non-minimal coupling in the context of cosmology. We pay special attention to the question of how far this model could reproduce the observational fact of our universe.

PACS numbers:

Keywords:

∗ Electronic address: sergiobravom@javeriana.edu.co

† Electronic address: marek.nowakowski@ictp-saifr.org

‡ Electronic address: contact@rvlobato.com

§ Electronic address: davide.batic@ku.ac.ae

## I. INTRODUCTION

Since the discovery of the accelerated expansion of the universe, cosmologists have developed the Λ CDM model, a standard cosmological model rooted in Einstein's general relativity [1]. This model incorporates a positive cosmological constant Λ [2], the FriedmannRobertson-Walker metric [8], and the concept of Dark Matter (DM) [9]. The otherwise robust model is not completely without problems. To start with, despite numerous efforts, a candidate for DM has not yet been found. The second problem is coined as Hubble tension, careful phrasing of the fact that different measurements of the Hubble constant yield different results [10]. In the future, as we seek to explain these issues, we may require a new theoretical model. This could involve modifying Einstein's gravity [11], replacing the cosmological constant with different models of Dark Energy (DE) [12], and exploring alternative candidates for DM [13]. The number of extensions of general relativity as well as the number of DE models is, of course, quite large [4, 6, 14-33].

Despite the remarkable success of the Λ Cold Dark Matter (ΛCDM) model in explaining a wide array of cosmological observations, from the cosmic microwave background (CMB) anisotropies to the large-scale structure of the universe, several significant challenges remain unresolved. The nature of dark energy, epitomized by the cosmological constant Λ, poses a profound theoretical conundrum, known as the 'cosmological constant problem' [2]. The astonishingly small value of Λ, as required to explain the observed acceleration of the universe, stands in stark contrast to theoretical predictions from quantum field theory, typically larger by many orders of magnitude [2].

Moreover, the ΛCDM model does not escape from the Hubble tension, a growing discrepancy between the values of the Hubble constant, H 0 , measured directly from local astronomical observations and those inferred from the CMB under the ΛCDM framework [3]. These inconsistencies point to potential physics beyond the standard model, either in the form of new particles or fields, or through modifications to the theory of general relativity itself.

In response to these challenges, a plethora of modified gravity theories have been proposed as alternatives to general relativity, aiming to provide a more comprehensive theoretical framework that can naturally incorporate the phenomena attributed to dark energy and dark matter. Theories such as f ( R ) gravity, where the Ricci scalar R in the Einstein-Hilbert

action is replaced by a function f ( R ), offer promising ways to explain cosmic acceleration without the need for a cosmological constant [4, 14].

Expanding further, f ( R,T ) theories introduce a coupling between matter, represented by the trace of the energy-momentum tensor T , and geometry, providing a framework to explore the effects of such couplings on the dynamics of the universe [6]. The f ( R, L m ) models we focus on in this paper extend this idea by incorporating the matter Lagrangian L m directly into the gravitational action, potentially offering new insights into the interaction between dark matter and dark energy, and their impact on the evolution of cosmic structures [7].

Our motivation for selecting f ( R, L m ) gravity stems from its ability to unite the geometrical and material sectors of the universe in a single, coherent theoretical framework. This approach not only allows for the exploration of the cosmic acceleration and dark matter problems from a new angle but also provides a platform for testing the limits of Einstein's general relativity on cosmological scales. By investigating the cosmological implications of f ( R, L m ) gravity, we aim to contribute to the ongoing discussion on viable alternatives to ΛCDM, exploring whether these theories can offer a more satisfactory explanation of observational phenomena without some of the fine-tuning issues that plague the standard model.

It is therefore a priori not clear which class of models is theoretically preferred over the others. We think that simplicity of the modification or, in other words, Occam's razor, could serve us here as guiding principle. Within a range of models, this possibility seems quite likely. Among the wide class of models based on a Lagrangian function of the type f ( R ), where R is the Ricci scalar [34], one could choose to focus on R 2 gravity [35] and the resulting cosmological models [36]. Recently, another class has been examined, based on f ( R, L m ) where L m [7, 37, 38] represents the matter Lagrangian, often taken to be the energy density ρ , pressure p , or the trace of the energy-momentum tensor T -with all these quantities being diffeomorphically invariant. In this class, an additional term in the EinsteinHilbert Lagrangian of the form σR L m would be considered a mild extension of the previous model. Moreover, it provides a straightforward coupling of the geometry encoded in R with matter represented astrophysically by L m . The astrophysical and other implications of such a theory have been examined in [39-42] with consequences for Neutron star and White Dwarf physics. Some cosmological aspects with specific choices of L m have been considered in [43-48]. Except for [48], the choice of the function f ( R,L m ) is different from ours. The

f ( R, L m ) theories have their own theoretical peculiarities that require careful examination. One notable aspect is that choosing L m = ± ρ, ± p, T does not lead, through standard metric variation to the energy-momentum tensor of a perfect fluid T µν [49-52]. On the other hand, a constrained variation based on the Lagrangian with constraints can result in the prefect fluid energy-momentum tensor. However, this comes at the cost of T µν no longer being conserved. We show that explicitly using the new Friedmann equation in f ( R, L m ).

In solving the new Friedmann equations, we try to make the model resemble our present universe. While it seems possible to achieve this by choosing adequate initial values, the evolution into the past and future can bring some surprises, as will be shown below. One of the simplest and model independent feature of our universe is the lower limit on its age as given for example by the redshift of the oldest galaxies [53]. A second constraint independent of the model comes from uranium decays [54]. If a cosmological model, as appealing as it may be, cannot reproduce these facts, it is then certainly not a good candidate to describe our universe. We will pay a special attention to these questions while examining the details of the model under discussion.

The paper is organized as follows. In Section II, we outline the basics of f ( R, L m ) theory with minimal coupling, denoted as σR L m . In this section, we also present the new Friedmann equations with a cosmological constant in a dimensionless form suitable for numerical integration. In Section III, we present the generalization of one of the Friedmann equations. In its standard form, it reads Ω Λ 0 , +Ω m, 0 = 1. However, in the context of f ( R, L m ) theory, it generalizes to Ω Λ 0 , + Ω m, 0 + Ω σ, 0 = 1. In the same section, we discuss the possible initial values, a necessary undertaking, since the new Friedmann equation contains the second derivative of the density. We believe it is necessary to provide a brief overview of the standard cosmological model in its analytical form (see Section IV). This is important because in Section V, we will present the numerical results and compare the new model with the standard one. In Section VI, we will draw our conclusions.

## II. THE f ( R, L m ) THEORY

As a generalization of the standard Einstein-Hilbert Lagrangian in General Relativity, i.e.

$$\mathcal { L } _ { G R } = \frac { 1 } { 2 \kappa } R + \mathcal { L } _ { m },$$

where κ = 8 πG N , a new Lagrangian has been proposed by making use of the invariants R and L m . It makes use of a general function, f ( R, L m ), and reads simply with an action [7, 37, 38]

$$S = \int d ^ { 4 } x \sqrt { - g } f ( R, \mathcal { L } _ { m } ).$$

From the formal variation of the action ( δ S g = 0) with respect to the metric, the following field equations are obtained [43, 55]

$$f _ { R } R _ { \mu \nu } + ( g _ { \mu \nu } \Box - \nabla _ { \mu } \nabla _ { \nu } ) f _ { R } - \frac { 1 } { 2 } ( f - f _ { L _ { m } } \mathcal { L } _ { m } ) g _ { \mu \nu } = \frac { 1 } { 2 } f _ { L _ { m } } T _ { \mu \nu },$$

where f R = ∂f ∂R , f L m = ∂f ∂L m and the energy-momentum tensor T µν is given by

$$\psi \varkappa _ { m } & & T _ { \mu \nu } = \frac { - 2 } { \sqrt { - g } } \frac { \delta ( \sqrt { - g } \mathcal { L } _ { m } ) } { \delta g ^ { \mu \nu } }. & &$$

An alternative form of these field equations may be obtained by taking the trace of (3) and solving for □ f R , namely

$$\Box f _ { R } = \frac { 1 } { 6 } f _ { L _ { m } } T + \frac { 2 } { 3 } ( f - f _ { L _ { m } } L _ { m } ) - \frac { 1 } { 3 } f _ { R } R.$$

When replacing it back into (3), we obtain the equation

$$f _ { R } G _ { \mu \nu } + \frac { 1 } { 6 } ( f + f _ { R } R - f _ { L _ { m } } L _ { m } ) g _ { \mu \nu } - \frac { 1 } { 2 } f _ { L _ { m } } \left ( T _ { \mu \nu } - \frac { 1 } { 3 } T g _ { \mu \nu } \right ) - \nabla _ { \mu } \nabla _ { \nu } f _ { R } = 0.$$

Here, G µν = R µν -1 2 g µν R represents the Einstein tensor. As mentioned in the Introduction, the particular form of f ( R, L m ) we wish to work with is the same as the one presented in [40]

$$f ( R, \mathcal { L } _ { m } ) = \frac { R } { 2 \kappa } + \mathcal { L } _ { m } + \sigma R \mathcal { L } _ { m },$$

where σ is the parameter which determines the coupling between matter and geometry.

With respect to the choice of the matter Lagrangian L m that reproduces the perfect fluid Energy-Momentum tensor there is a current debate [49, 51]. While some authors make the choice L m = ± ρ , others suggest L m = ± P [50] with particular values for the polytropic index. With the choice in [40], i.e. L m = -P , the field equations read

$$\left ( \frac { 1 } { 2 \kappa } - \sigma P \right ) R _ { \mu \nu } - \sigma \left [ g _ { \mu \nu } \Box - \nabla _ { \mu } \nabla _ { \nu } \right ] P - \frac { 1 } { 4 \kappa } R g _ { \mu \nu } = \frac { 1 } { 2 } \left ( 1 + \sigma R \right ) T _ { \mu \nu },$$

or in the alternative form

$$( 1 - 2 \kappa \sigma P ) G _ { \mu \nu } + \frac { 1 } { 3 } R g _ { \mu \nu } - \frac { \kappa \sigma P } { 3 } R g _ { \mu \nu } - \kappa ( 1 + \sigma R ) \left ( T _ { \mu \nu } - \frac { 1 } { 3 } T g _ { \mu \nu } \right ) + 2 \kappa \sigma \nabla _ { \mu } \nabla _ { \nu } P = 0, \ ( 9 )$$

which is consistent with the equations displayed above. The conservation of the energymomentum tensor is given in general by [55]

$$\nabla ^ { \mu } T _ { \mu \nu } = 2 \nabla ^ { \mu } \ln \left [ f _ { L _ { m } } ( R, \mathcal { L } _ { m } ) \right ] \frac { \partial \mathcal { L } _ { m } } { \partial g ^ { \mu \nu } }.$$

This suggests that choosing the matter Lagrangian to be proportional to ρ or P leads to the conservation of the energy-momentum tensor. However, this is only the case when relying directly on equation (4). Unfortunately, given the choices of the matter Lagrangian as mentioned above, it does not reproduce the perfect fluid energy-momentum tensor. To illustrate this point, if we start with L m = ± ρ , the energy-momentum tensor defined in (4) will not contain the pressure term. On the other hand, by using a constrained variation in the metric [50-52] or alternatively, a constrained Lagrangian formalism, we can explicitly show that, for example, choosing the energy density as the matter Lagrangian yields a perfect fluid energy-momentum tensor. However, its conservation no longer holds, as we introduce the metric through the constraints.

There is a certain confusion on this issue in the literature [48]. This is a point of divergence between our approach and others. Therefore, we explicitly demonstrate the non-conservation of T µν . The equations of motion, as presented in (3), (8), and (9), remain valid when using the constrained formalism. In this formalism, T µν is the energy-momentum tensor of perfect fluid, which, in general, is not conserved and is defined as

$$T _ { \mu \nu } = ( \rho + P ) u _ { \mu } u _ { \nu } - P g _ { \mu \nu }.$$

̸

As mentioned earlier, the model under discussion displays a non-conservation of the energymomentum tensor, specifically ∇ µ T µν = 0 (at least in general). This non-conservation is also a characteristic of other gravity models such as Rastall Gravity [56, 57], and has been advocated in modular gravity models [58, 59]. From the field equations (3) we can obtain the modified Friedmann equations in the flat ( k = 0) Friedmann-Lemaitre-Robertson-Walker metric

$$d s ^ { 2 } = g _ { \mu \nu } d x ^ { \mu } d x ^ { \nu } = d t ^ { 2 } - a ^ { 2 } ( t ) \left ( d r ^ { 2 } + r ^ { 2 } d \theta ^ { 2 } + r ^ { 2 } \sin ^ { 2 } \theta d \phi ^ { 2 } \right ).$$

By taking the 0-0 components of the new field equations, one obtains [43]

$$3 H ^ { 2 } f _ { R } + \frac { 1 } { 2 } ( f - f _ { R } R - f _ { L _ { m } } \mathcal { L } _ { m } ) + 3 H \dot { f } _ { R } = \frac { 1 } { 2 } f _ { L _ { m } } \rho,$$

where H is the Hubble parameter given by H = ˙ a/a . The space components (i.e. µ = , i ν = ) give j

$$\L ^ { \prime \prime } \circ \L ^ { \prime } \dot { H } f _ { R } + 3 H ^ { 2 } f _ { R } - \ddot { f } _ { R } - 3 H \dot { f } _ { R } + \frac { 1 } { 2 } ( f _ { L _ { m } } \mathcal { L } _ { m } - f ) = \frac { 1 } { 2 } f _ { L _ { m } } P,$$

where the ¨ f R term comes from the double derivatives in the field equations ( □ and ∇ ∇ µ ν ). For the model we have chosen, we have

$$f _ { R } = \frac { 1 } { 2 \kappa } + \sigma \mathcal { L } _ { m }, \quad f _ { L _ { m } } = 1 + \sigma R, \quad \dot { f } _ { R } = \sigma \dot { \rho }, \quad \ddot { f } _ { R } = \sigma \ddot { \rho }.$$

Together with the choice L m = ρ , the new Friedmann equations read

$$H ^ { 2 } _ { \underset { \vdots } { \dots } } \, = \, \kappa \frac { \rho } { 3 } + \sigma \kappa \left [ 4 \dot { H } \rho + 6 H ^ { 2 } \rho - 2 H \dot { \rho } \right ],$$

$$\frac { \ddot { a } } { a } \, = \, - \frac { \check { \kappa } } { 6 } \, ( \rho + 3 P ) - \sigma \kappa [ \dot { H } ( 3 P + \rho ) + 6 H ^ { 2 } P + \ddot { \rho } + 2 H \dot { \rho } ].$$

In general relativity, the two Friedmann equations lead to the conservation law ˙ + 3 ρ H ρ ( + P ) = 0. It is known that the same conservation law follows from ∇ µ T µν = 0 [1]. It is straightforward to derive from (16) and (17) the following relation

$$\dot { \rho } + 3 H ( \rho + P ) = - 6 \sigma \left [ H \dot { H } ( 3 P + 1 1 \rho ) + \dot { \rho } ( 3 H ^ { 2 } + \dot { H } ) + 6 H ^ { 3 } ( P + \rho ) + 2 \ddot { H } \rho \right ], \quad ( 1 8 )$$

which clearly gives back the conservation law (and therefore, the conservation of the energymomentum tensor) in the case σ = 0. The nature of the non-conservation of the energymomentum tensor is explicitly proved here. It should have been addressed also in [45, 46, 48] where the authors consider cosmological models of f ( R,L m ) theories.

Commonly, equations (16) and (17) are known as (modified) Friedmann equations. In some papers [60] equation (17) is also called Raychaudhuri equation. Both equations are related by (18).

In rewriting the Friedmann equations in a dimensionless form, specific choices were made to align the theoretical framework with observational benchmarks and simplify the numerical analysis. By introducing χ = ρ ρ crit , we scale the density to the critical density, which is pivotal in determining the curvature of the universe. Similarly, scaling the Hubble parameter by its current value H 0 through h 2 = H 2 H 2 0 allows us to examine the evolution of the expansion rate relative to its present value. The introduction of ˜ = σ H σ 2 0 helps in assessing the impact of the coupling constant σ within the typical energy scales of the current universe. Lastly, by using η = H t 0 , time derivatives are normalized to the current expansion rate, making

them dimensionless and easier to handle numerically. These transformations are not only mathematically convenient but also provide a direct link between theoretical predictions and observable quantities. To this end, we introduce

$$\chi \equiv \frac { \rho } { \rho _ { \text{crit} } } = \frac { \kappa \rho } { 3 H _ { 0 } ^ { 2 } }, \quad h ^ { 2 } \equiv \frac { H ^ { 2 } } { H _ { 0 } ^ { 2 } }, \quad \tilde { \sigma } \equiv H _ { 0 } ^ { 2 } \sigma, \quad \eta = H _ { 0 } t, \quad \rho _ { \text{crit} } = 3 H _ { 0 } ^ { 2 } / \kappa, \quad ( 1 9 )$$

where ρ 0 is a reference value of ρ at the present epoch. Similar comment applies to the Hubble parameter H 0 . In addition, we write the time derivatives in terms of derivatives with respect to the dimensionless parameter η by means of

$$\mathring { A } \equiv \frac { d A } { d \eta } = \frac { d A } { d t } \frac { d t } { d \eta } = \dot { A } \frac { 1 } { H _ { 0 } }.$$

In the context of our model, the barotropic index γ characterizes the equation of state (EOS) of the cosmological fluid, namely P = ( γ -1) ρ . This parameter is crucial for defining the thermodynamic properties of the universe's contents, where γ = 1 corresponds to a universe dominated by non-relativistic matter (the choice we make here for the calculations), γ = 4 3 to one dominated by radiation, and / γ = 0 to the cosmological constant scenario with vacuum energy. Each value of γ distinctly affects the evolution dynamics of the universe, as reflected in the modified Friedmann equations. By means of the aforementioned (EOS), equations (16), (17) and (18) take the form

$$h ^ { 2 } \, = \, \chi + \tilde { \sigma } [ 1 8 h ^ { 2 } \chi + 1 2 \mathring { h } \chi - 6 h \mathring { \chi } ],$$

$$h ^ { 2 } + \mathring { h } \ = \ - \frac { 1 } { 2 } ( 3 \gamma - 2 ) \chi - \underset { \circ } { \tilde { \sigma } } [ 3 ( 3 \gamma - 2 ) \mathring { h } \chi + 1 8 ( \gamma - 1 ) h ^ { 2 } \chi + 3 \overset { \circ } { \chi } ^ { \circ } + 6 h \overset { \circ } { \chi } ], \quad \ ( 2 2 )$$

$$3 \mathring { \chi } + 9 \gamma h \chi \, = \, - 6 \tilde { \sigma } [ ( 9 \gamma + 2 4 ) h \mathring { h } \chi + 3 \mathring { \chi } ( 3 h ^ { 2 } + \mathring { h } ) + 1 8 \gamma h ^ { 3 } \chi + 6 \mathring { h } \chi ].$$

Since we want to compare the modified cosmological model with standard cosmology we should also include the cosmological constant Λ. This can be done without much effort by considering the action

$$S _ { f ( R, \mathcal { L } _ { m } ) + \Lambda } = \int d ^ { 4 } x \sqrt { - g } \left ( f _ { \Lambda } ( R, \mathcal { L } _ { m } ) - \frac { \Lambda } { \kappa } \right )$$

with

$$f _ { \Lambda } ( R, \mathcal { L } _ { m } ) = \frac { R } { 2 \kappa } + \mathcal { L } _ { m } + \sigma R \mathcal { L } _ { m }.$$

This gives the new Einstein equations

$$f _ { R } R _ { \mu \nu } + ( g _ { \mu \nu } \Box - \nabla _ { \mu } \nabla _ { \nu } ) f _ { R } - \frac { 1 } { 2 } ( f - f _ { L _ { m } } L _ { m } ) g _ { \mu \nu } + \frac { \Lambda } { 2 \kappa } g _ { \mu \nu } = \frac { 1 } { 2 } f _ { L _ { m } } T _ { \mu \nu }$$

as well as the generalized Friedmann equations

$$3 H ^ { 2 } f _ { R } + \frac { 1 } { 2 } ( f - f _ { L _ { m } } \mathcal { L } _ { m } - f _ { R } R ) + 3 H \dot { f } _ { R } - \frac { \Lambda } { 2 \kappa } \, = \, \frac { 1 } { 2 } f _ { L _ { m } } \rho,$$

$$\dot { H } f _ { R } + 3 H ^ { 2 } f _ { R } - \ddot { f } _ { R } - 3 H \dot { f } _ { R } + \frac { 1 } { 2 } ( f _ { L _ { m } } \mathcal { L } _ { m } - f ) + \frac { \ddot { \Lambda } } { 2 \kappa } \, = \, \frac { \dot { 1 } } { 2 } f _ { L _ { m } } P.$$

With the explicit choice of f Λ the modified Friedmann equations simplify to

$$H ^ { 2 } _ { \substack { \dots \\ \div } } \, = \, \kappa _ { \frac { \rho } { 3 } } + \sigma \kappa \left [ 4 \dot { H } \rho + 6 H ^ { 2 } \rho - 2 H \dot { \rho } \right ] + \frac { \Lambda } { 3 }, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,$$

$$\frac { \ddot { a } } { a } \ = \ - \frac { \lambda } { 6 } \left ( \rho + 3 P \right ) - \sigma \kappa [ \dot { H } ( 3 P + \rho ) + 6 H ^ { 2 } P + \ddot { \rho } + 2 H \dot { \rho } ] + \frac { \Lambda } { 3 }.$$

Their dimensionless form can be written by using the dimensionless parameter Ω Λ = ρ vac /ρ crit as follows

$$h ^ { 2 } \, = \, \chi + \tilde { \sigma } [ 1 8 h ^ { 2 } \chi + 1 2 \mathring { h } \chi - 6 h \mathring { \chi } ] + \Omega _ { \Lambda },$$

$$h ^ { 2 } + \mathring { h } \, = \, - \frac { 1 } { 2 } ( 3 \gamma - 2 ) \chi - \tilde { \sigma } [ 3 ( 3 \gamma - 2 ) \mathring { h } \chi + 1 8 ( \gamma - 1 ) h ^ { 2 } \chi + 3 \mathring { \chi } + 6 h \mathring { \chi } ] + \Omega _ { \Lambda }. \quad ( 3 2 )$$

These equations are now suitable for numerical integration.

Normally, the inclusion of Λ in a modified gravity theory might not be well justified as the modification itself can account for Dark Energy. But it appears that this is a too global statement and the role of Λ can vary from model to model. In the section on numerical results we will demonstrate that the cosmological constant does not affect the characteristic feature of the model, i.e., its shortcoming to produce a reasonable lifetime of the universe.

## III. THE INITIAL VALUES AND THE ∑ i Ω = 1 i RELATION

The Friedman equations are now first order in h and second order in χ . Therefore, we need three initial values, which we provide at the time t 0 corresponding to the present epoch. We have

$$h _ { 0 } \equiv h ( t _ { 0 } ) = \frac { H _ { 0 } } { H _ { 0 } } = 1, \quad \chi _ { 0 } \equiv \chi ( t _ { 0 } ) = \frac { \kappa \rho _ { 0 } } { 3 H _ { 0 } ^ { 2 } }, \quad \mathring { h } _ { 0 } \equiv \mathring { h } ( t _ { 0 } ) = \frac { \dot { H } _ { 0 } } { H _ { 0 } ^ { 2 } }, \quad \mathring { \chi } _ { 0 } \equiv \mathring { \chi } ( t _ { 0 } ) = \frac { \kappa \dot { \rho } _ { 0 } } { 3 H _ { 0 } ^ { 3 } }.$$

We look for model independent measurements for the density of the Universe ( ρ 0 ), the Hubble parameter ( H 0 ) and its first and second derivatives ( ˙ H 0 , ¨ H 0 ). For ˙ H we use the deceleration parameter q , given by [62]

$$q \equiv - \frac { \ddot { a } a } { \dot { a } ^ { 2 } }, \quad \rightarrow \quad \frac { \dot { H } } { H ^ { 2 } } = - ( 1 + q )$$

and thus,

$$\dot { H } _ { 0 } = - H _ { 0 } ^ { 2 } ( 1 + q _ { 0 } ),$$

where q 0 = ( q t = t 0 ) has been obtained in a model independent way [61]. The jerk parameter will be useful for ¨ H 0 appearing in the second modified Friedmann equation. This is given by [62, 63]

$$j = \frac { \ddot { a } } { a H ^ { 3 } } \ \rightarrow j = \frac { \ddot { H } } { H ^ { 3 } } + 3 \frac { \dot { H } } { H ^ { 2 } } + 1$$

and therefore,

$$\ddot { H } _ { 0 } = ( j _ { 0 } - 1 ) H _ { 0 } ^ { 3 } - 3 \dot { H } _ { 0 } H _ { 0 } = ( j _ { 0 } + q _ { 0 } + 2 ) H _ { 0 } ^ { 3 }.$$

̸

To ensure that our exploration of the modified f ( R, L m ) gravity model remains grounded in observational reality, we adopt standard cosmological parameters that are widely accepted within the community. This choice allows us to rigorously test whether the modifications introduced by the theoretical framework can provide a viable alternative to the ΛCDM model, particularly in light of recent tensions and discrepancies such as those in measurements of the Hubble constant. By using these parameters, the model's predictions can be directly compared with those derived from both local and cosmological scales, providing a comprehensive evaluation of its empirical adequacy. For the numerical values, it is important to note that due to recent discrepancies in measurements of the Hubble constant, known as the Hubble Tension [10], there are several values for H 0 . For example we have H + 0 = 73,52 ± 1,62 km/s/Mpc and H -0 = 67,4 ± 0,5 km/s/Mpc, respectively. As for the current energy density of the Universe and the deceleration parameter, we use values found in [44, 45, 61]: χ 0 = ρ 0 ρ crit = 0 285 . ± 0 012, . q 0 = -0.545 ± 0.107, j 0 = 1 30 . ± 0 37 with the value for the . Hubble parameter taken as H 0 = 71 34 . ± 1 74 km s . -1 Mpc -1 . It is worth noting that these are specific parameter choices, but other sets of parameters are available in [61]. Possible values for these parameters are displayed in Table I. This particular choice of parameters that as of today we are as close as possible to the standard cosmology. For a standard reference value of ◦ χ 0 , we consider the standard cosmological model and use the continuity equation (23) with ˜ = 0, even in the presence of Λ = 0. σ This leads to

$$3 \mathring { \chi } + 9 \gamma h \chi = 0.$$

For the initial value we then obtain

$$\mathring { \chi } _ { 0 } = - 3 \gamma h _ { 0 } \chi _ { 0 } = - 3 \gamma \chi _ { 0 }.$$

TABLE I: Specific choices of the parameters according to [10, 44, 45, 61].

| Density Parameters   | Density Parameters   | Density Parameters   | Density Parameters   | Density Parameters   |
|----------------------|----------------------|----------------------|----------------------|----------------------|
| Data                 | SNIa                 | SNIa                 | Hubble               | Hubble               |
| Ω m, 0               | 0.285 ± 0.012        | 0.285 ± 0.012        | 0.239 ± 0.015        | 0.239 ± 0.015        |
| Ω Λ , 0              | 0.715                | 0.715                | 0.761                | 0.761                |
| ∆                    | 0.012                | 0.012                | 0.015                | 0.015                |
| Hubble data          | Hubble data          | Hubble data          | Hubble data          | Hubble data          |
| Model                | H exp                | GP                   | GA                   | ΛCDM                 |
| H 0                  | 73.88 ± 1.34         | 73.44 ± 1.40         | 71.34 ± 1.74         | 72.08 ± 1.06         |
| q 0                  | -1.070 ± 0.093       | -0.856 ± 0.111       | -0.545 ± 0.107       | -0.645 ± 0.023       |
| j 0                  | 3.00 ± 0.62          | 1.30 ± 0.37          | 0.52 ± 0.24          | 1.00                 |
| Pantheon data        | Pantheon data        | Pantheon data        | Pantheon data        | Pantheon data        |
| Model                | H exp                | GP                   | GA                   | ΛCDM                 |
| H 0                  | 71.13 ± 0.46         | 71.92 ± 0.38         | 71.81 ± 1.14         | 71.84 ± 0.22         |
| q 0                  | -0.616 ± 0.105       | -0.558 ± 0.040       | -0.466 ± 0.244       | -0.572 ± 0.018       |
| j 0                  | 1.56 ± 0.74          | 0.85 ± 0.12          | 0.55 ± 1.65          | 1.00                 |

Since h 0 = 1, in the case of dust γ = 1, we have

$$\mathring { \chi } _ { 0 } = - 3 \chi _ { 0 }.$$

In mathematics, when solving differential equations, one typically assumes that all parameters and functions entering the equation, as well as the initial values are known and given. However, in physics, the situation can become more complex because these parameters and initial values may not be precisely known. When considering a differential equation at the point where the initial values are applied, it leads to a relationship between these initial values and the parameters of the equation. The Friedmann equations provide a notable example of this phenomenon.

The first modified Friedmann equation can be expressed in terms of the density parameters Ω . i To this purpose, let us start by defining the Ω m parameter as

$$\Omega _ { m } = \frac { \rho } { \rho _ { \text{crit} } } = \frac { \kappa \rho } { 3 H ^ { 2 } }.$$

This, in turn, implies that χ = Ω m and ρ = Ω m ρ crit = Ω 3 H 2 κ . Thus, our first Friedmann equation can be reformulates as follows

$$H ^ { 2 } = \Omega _ { m } H ^ { 2 } + \Omega _ { m } H ^ { 2 } ( 1 8 \sigma H ^ { 2 } ) + \Omega _ { m } H ^ { 2 } ( 1 2 \sigma \dot { H } ) - \sigma \kappa H \dot { \rho } + \frac { \Lambda } { 3 }.$$

Recognizing this equation as

$$\Omega _ { m } + \Omega _ { \Lambda } + \tilde { \sigma } \left [ 1 8 \Omega _ { m } H ^ { 2 } + 1 2 \Omega _ { m } \dot { H } - \frac { \kappa \dot { \rho } } { H } \right ] = 1,$$

we can express it in a more compact form. Finally, defining the term proportional to σ as Ω , we arrive at σ

$$\Omega _ { m, 0 } + \Omega _ { \Lambda, 0 } + \Omega _ { \sigma, 0 } = 1,$$

where we have taken the values at t = t 0 . Making use of (35) and (40) the above equation reduces to

$$\Omega _ { m, 0 } + \Omega _ { \Lambda, 0 } + \tilde { \sigma } \Omega _ { m, 0 } ( 1 5 - 1 2 q _ { 0 } ) = 1.$$

Equation (45) establishes a relationship between initial values and parameters within the modified Friedmann equation. Moreover, it introduces the deceleration parameter q 0 , which is a measured quantity, though not strictly necessary for solving the first order modified Friedmann equation in H . Consequently, there are generally two approaches to solve the new Friedmann equations. The first one involves specifying h 0 = 1, χ 0 , ◦ χ 0 , and Ω Λ while varying ˜. σ Equation (45) can then be used to determine q 0 . The second approach relies on the measured value of q 0 and fixes the other parameters as previously described. In this case, equation (45) provides the parameter ˜. In particular, we note that in the case without σ coupling between R and L m , the relative densities today read

$$\Omega _ { m, 0 } + \Omega _ { \Lambda, 0 } = 1.$$

Numerically, however, there is an uncertainty in the obtained values, reflected in the relation Ω m, 0 +Ω Λ 0 , = 1 ± ∆, where ∆ represents this uncertainty. We can exploit it in our model by setting, as shown in equation (44),

$$\tilde { \sigma } = \pm \frac { \Delta } { \Omega _ { m, 0 } ( 1 5 - 1 2 q _ { 0 } ) }.$$

## IV. THE STANDARD COSMOLOGY

Before embarking on a discussion of the numerical solutions it is convenient to briefly outline the analytical picture of the standard cosmology. This makes sense since we will be comparing the two scenarios.

The dimensionless solution for h Λ+ for the standard cosmological model can be obtained by taking ˜ = 0 in equations (21) and (22) with σ γ = 1, namely

$$h ^ { 2 } \, = \, \chi + \Omega _ { \Lambda },$$

$$1$$

$$h ^ { 2 } + \mathring { h } \, = \, - \frac { 1 } { 2 } \chi + \Omega _ { \Lambda }.$$

Combining the above equations, we arrive at a Riccati equation of the form

$$\mathring { h } = - \frac { 3 } { 2 } h ^ { 2 } + \frac { 3 } { 2 } \Omega _ { \Lambda }.$$

We can solve the differential equation above by setting

$$\int _ { h ^ { \prime } = 1 } ^ { h ^ { \prime } = h } \frac { d h ^ { \prime } } { \Omega _ { \Lambda } - ( h ^ { \prime } ) ^ { 2 } } = \frac { 3 } { 2 } \int _ { \eta ^ { \prime } = \eta _ { 0 } } ^ { \eta ^ { \prime } = \eta } d \eta ^ { \prime },$$

where η 0 = H t 0 0 is the time at which the Hubble parameter is H 0 , which, in turn, is related to the present values of the parameters used. To solve the integral above, we redefine x = h ′ √ Ω Λ and write √

$$\frac { 1 } { \sqrt { \Omega _ { \Lambda } } } \int _ { x = 1 / \sqrt { \Omega _ { \Lambda } } } ^ { x = h / \sqrt { \Omega _ { \Lambda } } } \frac { d x } { 1 - x ^ { 2 } } = \frac { 3 } { 2 } ( \eta - \eta _ { 0 } ).$$

One possible solution for the integral on the left-hand side is

$$\int \frac { d x } { 1 - x ^ { 2 } } = \frac { 1 } { 2 } \ln \left ( \frac { 1 + x } { 1 - x } \right ).$$

However, such a solution is valid for -1 &lt; x &lt; 1. This would imply that h &lt; √ Ω Λ and thus, χ = h 2 -Ω Λ &lt; 0, which would suggest an unphysical negative density. Another solution is

$$\int \frac { d x } { 1 - x ^ { 2 } } = \frac { 1 } { 2 } \ln \left ( \frac { x + 1 } { x - 1 } \right ).$$

It is valid whenever | x &gt; | 1. Moreover, it implies h 2 &gt; Ω Λ and, consequently χ = h 2 -Ω Λ &gt; 0, ensuring a positive and physically meaningful density. If we also use the fact that

$$\coth ^ { - 1 } x = \frac { 1 } { 2 } \ln \left ( \frac { x + 1 } { x - 1 } \right )$$

and solve for the integral (52), we end up with the result

$$\frac { 1 } { \sqrt { \Omega _ { \Lambda } } } \left [ \coth ^ { - 1 } \left ( \frac { h } { \sqrt { \Omega _ { \Lambda } } } \right ) - \coth ^ { - 1 } \left ( \frac { 1 } { \sqrt { \Omega _ { \Lambda } } } \right ) \right ] = \frac { 3 } { 2 } ( \eta - \eta _ { 0 } ).$$

Solving for h in the γ = 1 case leads to

$$h _ { \Lambda + } ( \eta ) = \sqrt { \Omega _ { \Lambda } } \coth \left [ \frac { 3 } { 2 } \sqrt { \Omega _ { \Lambda } } ( \eta - \eta _ { 0 } ) + \coth ^ { - 1 } \left ( \frac { 1 } { \sqrt { \Omega _ { \Lambda } } } \right ) \right ].$$

For the general γ case, the equation above reads

$$h _ { \Lambda + } ( \eta ) = \sqrt { \Omega _ { \Lambda } } \coth \left [ \frac { 3 \gamma } { 2 } \sqrt { \Omega _ { \Lambda } } ( \eta - \eta _ { 0 } ) + \coth ^ { - 1 } \left ( \frac { 1 } { \sqrt { \Omega _ { \Lambda } } } \right ) \right ].$$

In the standard notation of the scale factor a Λ ( ) and the Hubble parameter t H Λ ( ), we get t

$$a _ { \Lambda + } ( t ) \, = \, \left ( \frac { 1 - \Omega _ { \Lambda } } { \Omega _ { \Lambda } } \right ) ^ { \frac { 1 } { 3 \gamma } } \sinh ^ { \frac { 2 } { 3 \gamma } } ( \mu t + \nu ),$$

$$H _ { \Lambda + } ( t ) \, = \, \sqrt { \frac { \Lambda } { 3 } } \coth ( \mu t + \nu ),$$

where

$$\mu = \frac { \gamma \sqrt { 3 \Lambda } } { 2 }, \quad \nu = \frac { 1 } { 2 } \ln \left ( \frac { 1 + \sqrt { \Omega _ { \Lambda } } } { 1 - \sqrt { \Omega _ { \Lambda } } } \right ) - \mu t _ { 0 }$$

for a positive cosmological constant Λ. For completeness, we also give the solution for Λ &lt; 0 which reads

$$a ( t ) _ { \Lambda - } \, = \, a _ { 0 } \, \underbrace { \left [ \frac { \cos \left ( \alpha - \beta t \right ) } { \cos \left ( \alpha - \beta t _ { 0 } \right ) } \right ] ^ { \frac { 2 } { 3 \gamma } },$$

$$H _ { \Lambda - } ( t ) \, = \, \sqrt { - \frac { \Lambda } { 3 } } \tan \left ( \alpha - \beta t \right )$$

with

$$\alpha = \beta t _ { 0 } + \tan ^ { - 1 } \left ( \sqrt { - \frac { 3 } { \Lambda } } H _ { 0 } \right ), \quad \beta = \frac { \gamma \sqrt { - 3 \Lambda } } { 2 }.$$

Without loss of generality we can set t 0 = 0. Then from (64) we can infer that 0 &lt; α &lt; π/ 2 and hence cos( α ) &gt; 0. To ensure that the scale factor is positive it suffices to choose t ∈ [ α -π/ 2 β , α + π/ 2 β ] corresponding to Big Bang and collapse.

We will explore later how this simple solutions evolves when we introduce a non-zero σ

.

Two important model-independent restrictions on a realistic cosmological model relate to the universe lifetime. In the standard cosmological model, we can calculate the lifetime

(neglecting the short radiation period) by setting the argument of the Hubble function to zero, which results in a singularity. Choosing η 0 = 0, we obtain

$$\lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad\lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \ \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \ \alpha \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \alpha \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \quad \lambda \ \lambda
 \quad \quad \quad$$

we get | η univ | ≃ 0 96 for Ω . ≃ 0 7. . This lifetime is not in conflict with model-independent estimates, such as those from the oldest galaxies [53] (which indicate the existence of galaxies some 400 million years after the Big Bang) and [54], which points to a lower limit of 12.5Gy.

## V. NUMERICAL RESULTS

̸

Assuming σ = 0, we have converted the system under discussion into a first order system which we solved by standard numerical methods. As an independent check, we used MAPLE18 routines and found that the results coincide.

When solving the modified Friedmann equations numerically, the question arises: What values should we take for ˜? σ The standard cosmology has two constants, the Newtonian one G N and the cosmological constant. Interestingly, the dimension of σ is the same as that of the Newtonian constant. If we take σ as nG N with n some number of order 1, ˜ would σ turn out extremely small. Therefore, it makes sense to initially try more moderate values for ˜. σ In the subsequent figures, the blue curve represents the standard cosmological model prediction.

In Figure 1, we have chosen the initial values as discussed in Section III, with Ω Λ = 0 7. . The chosen values of ˜ σ are relatively small to illustrate that the dimensionless Hubble function closely follows the one of the standard cosmological model around η close to zero. Several notable features can be observed in the plot, which are likely to persist for other values of ˜. σ First of all, the Hubble function reaches the singularity much earlier than in the standard model. This poses a challenge if we aim to adapt the model under discussion as a viable cosmological model of our universe. As mentioned earlier, one of the constraints on any cosmological model is its ability to be roughly compatible with a lower limit on the lifetime of the universe. This constraint arises from evidence such as the existence of galaxies just 400 million years after the Big Bang [53]. Therefore, any cosmological model must strive to be in agreement with the estimated lifetime of the universe as predicted by the standard cosmological model. Secondly, we notice another interesting feature. The Hubble function

FIG. 1: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η for a choice of the parameter ˜. σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-67b95cf3-f655-4a2b-88aa-d724b7040a7d.jpg)

becomes negative which implies a contraction of the universe. Therefore, the dimensionless density χ starts increasing again. This is known as the Big Crunch scenario which happens also in f ( R,L m ) models with the simple choice L m = ρ . The Big Crunch here is a result of the model and cannot be avoided.

FIG. 2: The same plot as in 1, but over a large timescale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-505515fb-b7be-4d6a-99c2-7fc2de79651c.jpg)

In Figure 2, we have plotted the same functions for two values of ˜ appearing also in σ figure 1, but over a larger scale of time. One can appreciate the steep rise of h at a relatively

modest negative values of η .

FIG. 3: The same plot as in 1, but for large values of ˜ . σ.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-714f4c61-9ab4-425a-ad17-57d8a455c22d.jpg)

In Figure 3, we have increased the values of ˜ but there is no visible effect on the time of σ the Big-Bang. However, for larger values of the parameter, the density becomes unphysically negative.

FIG. 4: The same plot as in 2 but form larger values of ˜. σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-e77963ea-319a-441e-9604-b106c2714e46.jpg)

This is confirmed in Figure 4 which shows the same functions, but over a larger timescale.

The differences between the models become then visible at large η .

Increasing the values of ˜ does not bring any new insight, but for the first time there σ

FIG. 5: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η for large values the parameter ˜. σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-a3e894ae-79a2-4893-a268-d818649ce855.jpg)

appears a cusp for ˜ = 5. σ MAPLE18 interpretes this as singularities. This is demonstrated in Figure 5.

FIG. 6: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η for small values of the parameter ˜. σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-b61c18e6-510a-4cb2-a930-13d78412f2c8.jpg)

In Figure 6, we return to relative small values of ˜ to show that the cusp behavior is σ quite common in the model. This is interesting from a purely theoretical point of view, but leaves doubts about the viability of the model over large time scales.

FIG. 7: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η for large negative values of the parameter ˜. σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-cf15830b-c474-4830-b1d7-2a35dafe7fee.jpg)

Finally, we turn our attention the the behaviour of h and χ for negative values of ˜ which σ is depicted in Figures 7, 8, and 9.

FIG. 8: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η for small negative values of the parameter ˜. σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-cb8264fa-12da-4514-9ece-1922b1286a56.jpg)

In Figure 7, we encounter a universe whose contraction slows down to turn into a expanding universe at η around zero. Again many cusps accompany this behavior which could well mean that universe undergoes a singularity there.

FIG. 9: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η for very small values of the parameter ˜. σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-e8e14331-d72d-46d5-abca-2a804832c88f.jpg)

For small and very small values of negative ˜ we have again negative densities. σ

We have also explored a second approach for selecting the parameter ˜. σ More precisely, we utilized (44) and considered the error bars associated with measured values in standard cosmology. Figure 10 provides an example where we calculated ˜ with high precision, staying σ within the range of uncertainty. This allowed us to examine how sensitive the model is to the choice of parameters. The lifetime of the universe comes out approximately | η univ | ≃ 0 4, . which is too small to account for the existence of the first galaxies in our universe. The Big Crunch appears to be unavoidable again.

In the third approach, we used equation (45) while keeping the standard value for Ω m, 0 , but varying Ω Λ and q 0 . We then solved the equation to find ˜. σ An example is illustrated in Figure 11. This variation led to a slight increase in the universe lifetime. Motivated by this observation, we explored scenarios with a negative cosmological constant. Two such examples are noteworthy. Setting q 0 = -0 5 . and Ω Λ = -50 resulted in ˜ σ ≃ 10 and a lifetime of | η univ | ≃ 0 51. . Alternatively, when choosing a deceleration parameter of q 0 = -1, we obtained ˜ σ ≃ 7 8 and a lifetime of 0 56. . . However, it appears challenging to significantly extend the model lifetime beyond these values.

In summary, this model represents a mild extension of Einstein gravity, particularly when viewed from the perspective of its Lagrangian. Surprisingly, the cosmology given by this

FIG. 10: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η for values of ˜ resulting from equation (44). σ

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-ff5169fa-3299-4ce7-9800-93c4adb2cfe8.jpg)

FIG. 11: The dimensionless Hubble function h and the dimensionless density χ versus the dimensionless time η using equation (45) with values of the parameters as indicated in the Figure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00580v1-figure-95a8aae0-1116-4946-b78c-fe0f83193487.jpg)

model for a spatially flat, homogenous universe is quite different from the standard case based on Einstein equations.

## VI. CONCLUSIONS

Proposing a new gravity theory which goes beyond Einstein leads inevitably to the examination of the new model in a cosmological context. In this context, theoretical advancements revealing new phenomena are equally important as empirical efforts to evaluate the model potential to supersede the standard cosmological framework. We have emphasized both aspects studying the extension f ( R, L m ) theory with an explicit coupling between geometry and matter in the form σRρ with ρ being the energy density. This appears as a natural choice used also in astrophysical context. We explored the model over a wide range of the coupling constant σ with initial data mimicking our present universe. Given that the model predicts a universe lifespan that is drastically shorter than the observed age of our universe (approximately 13.8 billion years), it fails to meet a fundamental criterion for a viable cosmological model. This discrepancy leads us to conclude that the model, in its current form, cannot adequately describe the observed universe. As such, the potential occurrence of a Big Crunch, while theoretically interesting, becomes a secondary concern because the model is already ruled out based on its inability to account for the observed age of the universe.

All in all, the findings highlight challenges in adapting the model to realistically represent the universe, particularly regarding the universe lifetime and the issue of negative densities in some scenarios.

Before comparing our work with existent literature, let us point out that once we derive the field equations from the Lagrangian principle, the Lagrangian should contain a matter part L m from which the perfect energy-momentum tensor can be derived. The latter seems only possible via a constrained variation with respect to the metric. The conservation of the energy-momentum tensor in theories of the type f ( R,L m ) is not guaranteed and requires special attention. In [48] the choice of the Lagrangian is f ( R,L m ) = R/ 2+ L m + σRL m with L m = -p . It is then difficult to study the dust scenario with p = 0 as it implies L m = 0. Furthermore, the energy-momentum tensor is not conserved in this choice in contrast to the authors' claim. In [46] the choice is f ( R,L m ) = R/ 2 + L α m with α being some number. It is not clear how we can derive from that a hydrodynamical energy-momentum tensor and, as explained above, whether we obtain a conservation law. The choice of the Lagrangian in [46] is also quite different from ours as it lacks the direct coupling between geometry (represented by R ) and matter (represented by L m ). This shows that such coupling has

special consequences. In [45] the Lagrangian reads f ( R,L m ) = R/ 2 + αL n m -β with α and β constants. Remarks similar to the case of [46] would apply also here.

- [1] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity , Wiley &amp; Sons, New York (1972).
- [2] S. Weinberg, The cosmological constant problem , Rev. Mod. Phys. 61 , 1 (1989); T. Padmanabhan, Cosmological Constant-the Weight of the Vacuum , Phys. Rept. 380 , 235 (2003).
- [3] L. Verde, T. Treu, and A.G. Riess, Tensions between the Early and the Late Universe , Nat. Astron. 3 , 891 (2019).
- [4] T. P. Sotiriou and V. Faraoni, f ( R ) theories of gravity , Rev. Mod. Phys. 82 , 451 (2010).
- [5] A. De Felice and S. Tsujikawa, f ( R ) theories , Living Rev. Rel. 13 : 3, (2010).
- [6] T. Harko, F. S.N. Lobo, S. Nojiri and S. D. Odintsov. f ( R,T ) gravity , Phys. Rev. D84 :024020 (2011).
- [7] T. Harko, and F.S.N. Lobo, f ( R,L m ) gravity Eur. Phys. J. C. 70 , 373 (2010).
- [8] H. P. Robertson, Kinematics and World-Structure , ApJ 82 , 284 (1935); A.G. Walker, On Milne's Theory of World-Structure , Proc. Lond. Math. Soc., 42 , 90 (1936).
- [9] G. Bertone, Particle Dark Matter: Observations, Models and Searches , Cambridge University Press, Cambridge (2013).
- [10] E. Di Valentino, O. Mena, S. Pan, L. Visinelli, We. Yang, A. Melchiorri, D. F. Mota, A. G. Riess, and J. Silk, In the Realm of Hubble Tension-a Review of Solutions , Class. Quantum Grav. 38 , 153001 (2021); W.D. Kenworthy, D. Scolnic, and A. Riess, The Local Perspective on the Hubble Tension: Local Structure Does Not Impact Measurement of the Hubble Constant , ApJ 875 , 145 (2019).
- [11] S. Capozziello, and V. Faraoni, Beyond Einstein Gravity: a Survey of Gravitational Theories for Cosmology and Astrophysics , Springer (2013).
- [12] V. Motta, M. A. Garc´ ıa-Aspeitia, A. Hern´ndez-Almada, J. Maga˜a, T. Verdugo, a n Taxonomy of Dark Energy Models , Universe 7 , 163 (2021).
- [13] N. C. M. Martens and M. King, Doing More with Less: Dark Matter and Modified Gravity in N. Mills Boyd, S. de Baerdemaeker, K. Heng and V. Matarese (eds.), Philosophy of Astrophysics: Stars, Simulations and the Struggle to Determine What is Out There , Springer (2023).

- [14] A. De Felice and S. Tsujikawa, f ( R ) theories , Living Rev. Rel. 13 : 3, (2010).
- [15] G. J. Olmo Palatini Approach to Modified Gravity: f ( R ) Theories and Beyond , Int. J. Mod. Phys. D20 :413-462, (2011).
- [16] R. Myrzakulov, FRW cosmology in F R,T ( ) gravity , Eur. Phys. J. bf C 72, 2203 (2012).
- [17] J. Wu, G. Li, T. Harko and S. Liang, Palatini formulation of f ( R,T ) gravity theory, and its cosmological implications , Eur. Phys. J. C 78 , 430 (2018).
- [18] N. Katirci and M. Kavuk, f ( R,T µν T µν gravity and Cardassian-like expansion as one of its consequences , Eur. Phys. J. Plus 129: 163 (2014).
- [19] M. Roshan and F. Shojai, Energy-Momentum Squared Gravity , Phys. Rev. D 94 , 044002 (2016).
- [20] C. V. R. Board and J.D. Barrow, Cosmological Models in Energy-Momentum-Squared Gravity , Phys. Rev. D 96 , 123517 (2017).
- [21] Y. Cai, S. Capozziello, M. De Laurentis and E. N. Saridakis, f ( T ) teleparallel gravity and cosmology , Rept.Prog.Phys. 79 no.4, 106901 (2016).
- [22] R. Yang, New types of f ( T ) gravity , Eur. Phys. J. C 71 , 1797 (2011)
- [23] S. Capozziello, V. F. Cardone, H. Farajollahi and A. Ravanpak, Cosmography in f ( T ) -gravity , Phys.Rev. D84 :043527 (2011).
- [24] L. Heisenberg, Review on f ( Q ) Gravity , Phys.Rept. 1066 (2024) 1-78.
- [25] W. Khyllep, J. Dutta, E. N. Saridakis, K. Yesmakhanova, Cosmology in f ( Q ) gravity: A unified dynamical system analysis at background and perturbation levels , Phys. Rev. D 107 , 044022 (2023).
- [26] K. Koussour, K. El Bourakadi, S.H. Shekh, S.K.J. Pacif and M. Bennai, Late-time acceleration in f ( Q ) gravity: Analysis and constraints in an anisotropic background , Ann. Physics, 445 (2022), 169092.
- [27] J. B. Jimenez, L. Heisenberg, T. Koivisto, and S. Pekar, Cosmology in f ( Q ) geometry , Phys. Rev. D 101 , 103507 (2020).
- [28] Y. X., Guangjie Li, T. Harko and Shi-Dong Liang f ( Q,T ) gravity , Eur. Phys. J. C 79 , 708 (2019).
- [29] M. Shiravand, S. Fakhry, and M, Farhoud Cosmological Inflation in f ( Q,T ) Gravity , Physics of the Dark Universe 37 101106 (2022).
- [30] Y. Xu, T. Harko, S. Shahidi and S. Liang Weyl type f ( Q,T ) gravity, and its cosmological

implications Eur. Phys. J. C 80 , 449 (2020).

- [31] S. Bahamonde et al. Teleparallel Gravity: From Theory to Cosmology , Rep. Prog. Phys. 86 026901 (2023)
- [32] Yu. N. Obukhov and J. G. Pereira, Metric-affine approach to teleparallel gravity , Phys. Rev. D 67 , 044016 (2003)
- [33] L. Heisenberg, M. Hohmann and S. Kuhn Homogeneous and isotropic cosmology in general teleparallel gravity , Eur. Phys. J. C 83 , 315 (2023).
- [34] S. Nojiri, S.D. Odintsov, and V.K. Oikonomou, Modified Gravity Theories on a Nutshell: Inflation, Bounce and Late-time Evolution , Phys. Rept. 692 , 1 (2017).
- [35] F. Payadeh, and M. Fathi, R 2 Theory of Gravity , J. Phys. Conf. Ser. 442 , 012053 (2013).
- [36] E. V. Arbuzova, A. D. Dolgov, and L. Reverberi, Cosmological evolution in R 2 gravity , JCAP 02 , 049 (2012).
- [37] T. Harko, and F.S.N. Lobo, Extensions of f ( R ) Gravity: Curvature-Matter Couplings and Hybrid Metric-Palatini Theory , Cambridge Monographs on Mathematical Physics, Cambridge University Press, Cambridge (2018).
- [38] T. Harko, and F.S.N. Lobo, Generalized Curvature-Matter Couplings in Modified Gravity , Galaxies 2 (3), 410 (2014).
- [39] N. Montelongo Garcia and F. S. N. Lobo, 'Wormhole geometries supported by a nonminimal curvature-matter coupling', Phys. Rev. D 82 , 104018 (2010)
- [40] R.V. Lobato, A. Carvalho, N.G. Kelkar, and M. Nowakowski, Massive white dwarfs in f ( R,L m ) gravity , Eur. Phys. J. C. 82 , 540 (2022).
- [41] R. V. Lobato, G. A. Carvalho, N. G. Kelkar, and M. Nowakowski, In-Stabilities of massive white dwarfs in modified gravity , Acta Phys. Pol. 142 , 315 (2022).
- [42] R. V. Lobato, G. A. Carvalhob, and C. A. Bertulani, Neutron stars in f ( R,L m ) gravity with realistic equations of state: joint-constrains with GW170817, massive pulsars, and the PSR J0030+0451 mass-radius from NICER data , Eur. Phys. J. C 81 , 1013 (2021).
- [43] L.V. Jaybhe, R. Solanki, S. Mandal, and P.K. Sahoo, Cosmology in f ( R,L m ) gravity , Phys. Lett. B 831 , 137148 (2022).
- [44] A. Pradhan, D.C. Maurya, G. K. Goswani, and A. Beesham, Modeling transit Dark Energy in f ( R, L m ) -gravity , Int. J. Geom. Methods Mod. 20 , 2350105 (2023).
- [45] D.C. Mauyura, Accelerating scenarios of massive universe in f ( R, L m ) -gravity , New Astron.

100 , 101974 (2023).

- [46] L.V. Jaybhaye, R. Solanki, S. Mandal, and P.K. Sahoo, Constraining Viscous Dark Energy Equation of State in f ( R,L m ) Gravity , Universe 9 , 163 (2023).
- [47] N.S. Kavya, V. Venkatesha, S. Mandal, and P.K. Sahoo, Constraining Anisotropic Comological Model in f ( R,L m ) Gravity , Phys. Dark Universe 38 , 101126 (2022).
- [48] B.S. Gonclaves, P.H.R.S. Moraes, and B. Mishra, Cosmology from Non-Minimal GeometryMatter Coupling ,. Fortschr. Phys. 71 , 2200153 (2023).
- [49] T. Harko, The matter Lagrangian and the energy-momentum tensor in modified gravity with non-minimal coupling between matter and geometry , Phys. Rev. D81 , 044021, (2010).
- [50] S. Mendoza, and S. Silva, The matter Lagrangian of an ideal fluid , Int. J. Geom. Methods Mod. 18 , 2150059 (2020).
- [51] O. Minazolli, and T. Harko, New derivation of the Lagrangian of a perfect fluid with a barotropic equation of state , Phys. Rev. D 86 , 087502 (2012).
- [52] D. Brown, Action Functionals for Relativistic Perfect Fluids , Class. Quant. Grav. 10 , 1579 (1993).
- [53] B. E. Robertson et al., Identification and properties of intense star-forming galaxies at redshifts z &gt; 10, Nat. Astron. 7 , 6 (2023); E. Curtis-Lake et al., Spectroscopic confirmation of four metal-poor galaxies at z = 10.3-13.2 , Nat. Astron. 7 , 622 (2023).
- [54] C. Sneden, The age of the universe , Nature 409 , 643 (2001); R Cayrel et al., Measurement of stellar age from uranium decay , Nature 409 , 691 (2001).
- [55] T. Harko, Thermodynamic interpretation of the generalized gravity models with geometrymatter coupling , Phys. Rev. D 90 , 044067 (2014).
- [56] H. Moradpour, Y. Heydarzade, F. Darabi, and I. G. Salako, A Generalization to the Rastall Theory and Cosmic Eras , Eur. Phys. J. C 77 , 259 (2017).
- [57] N. Katirci, S. Kumar, R. C. Nunes, B. Ozturk, and S. Sharma, Rastall gravity extension of the standard Λ CDM model: theoretical features and observational constraints , Eur. Phys. J. C 80 , 1050 (2020).
- [58] T. Maudlin, E. Okon, and D. Sudarsky, On the Status of Conservation Laws in Physics: Implications for Semiclassical Gravity , Stud. Hist. Phil. Sci. B 69 , 67 (2020).
- [59] S. J. Landau, M. Benetti, A. Perez, and D. Sudarsky, Cosmological constraints on unimodular gravity models with diffusion , arXiv:2211.07424 [astro-ph.CO] (preprint).

- [60] Parampreet Singh and S. K. Soni, 'On the relationship between the modifications to the Raychaudhuri equation and the canonical Hamiltonian structures', Class. Quant. Grav 33 125001 (2016)
- [61] A. Mehrabi, and M. Rezaei, Cosmographic parameters in model-independent approaches , ApJ 923 , 274, (2021).
- [62] M. Visser, Jerk, snap and the cosmological equation of state , Class. Quant. Grav. 21 , 2603 (2004).
- [63] M. Visser, Cosmography: Cosmology without the Einstein equations , Gen. Rel. Grav. 37 , 1541, (2005).
Sergio Bravo Medina Marek Nowakowski Ronaldo V. Lobato Davide Batic
Cosmologies in $f(R,\mathcal{L}_m)$ theory with non-minimal coupling between geometry and matter
Among the recent extensions to standard General Relativity,
$f(R,\mathcal{L}_m)$ gravity has risen an interest given the possibility of
coupling between geometry and matter. We examine the simplest model with
non-minimal coupling in the context of cosmology. We pay special attention to
the question of how far this model could reproduce the observational fact of
our universe.