---
title: f(R,L_m) Gravity Insights
url: https://www.emergentmind.com/topics/f-r-l_m-gravity
type: topic
---

# f(R,L_m) Gravity Insights

$f(R,L_m)$ gravity is a metric modified-gravity framework in which the gravitational Lagrangian is taken to be an arbitrary function of the Ricci scalar $R$ and the matter Lagrangian density $L_m$, rather than a sum of a purely geometric Einstein–Hilbert term and a minimally coupled matter term. In its basic form, the theory is defined by
$$
S=\int d^4x\,\sqrt{-g}\,f(R,L_m),
$$
and reduces to standard General Relativity (GR) for $f(R,L_m)=R/2+L_m$ in units $8\pi G=c=1$ [1008.4193]. The explicit curvature–matter dependence generically modifies the field equations, changes the covariant balance laws, and can induce non-geodesic motion through an extra force orthogonal to the four-velocity. In the literature represented here, the framework has been applied to late-time cosmology, dynamical-systems analyses, perturbative growth of structure, bouncing and baryogenic early-universe scenarios, white-dwarf equilibrium, and several wormhole constructions [1008.4193].

## 1. General action, field equations, and formal structure

The general metric variation of the action yields the field equations
$$
f_R\,R_{\mu\nu}+\bigl(g_{\mu\nu}\Box-\nabla_\mu\nabla_\nu\bigr)f_R
-\tfrac12\bigl(f-f_{L_m}L_m\bigr)g_{\mu\nu}
=\tfrac12\,f_{L_m}\,T_{\mu\nu},
$$
with
$$
f_R\equiv \frac{\partial f}{\partial R},\qquad
f_{L_m}\equiv \frac{\partial f}{\partial L_m},
$$
and
$$
T_{\mu\nu}=-\frac{2}{\sqrt{-g}}\,
\frac{\delta(\sqrt{-g}\,L_m)}{\delta g^{\mu\nu}}.
$$
The corresponding trace equation is
$$
R\,f_R+3\Box f_R-2\bigl(f-f_{L_m}L_m\bigr)=\tfrac12\,f_{L_m}\,T,
$$
where $T=g^{\mu\nu}T_{\mu\nu}$ [1008.4193; 2308.06519].

This structure interpolates between several familiar limits. The linear choice $f(R,L_m)=R/2+L_m$ reproduces the Einstein equations exactly, and in that case the generalized non-conservation law collapses to the standard $\nabla^\mu T_{\mu\nu}=0$ [2306.09387]. By contrast, nonlinear matter sectors such as $f(R,L_m)=R/2+L_m^\alpha$, mixed couplings such as $f(R,L_m)=R/2+(1+\lambda R)L_m$, and cosmological forms such as $f(R,L_m)=\lambda R+\beta L_m^\alpha+\eta$ introduce effective source terms that cannot be reabsorbed into standard GR without redefining the matter sector [2307.02498; 2403.17037; 2406.08303].

A persistent technical point is that the theory depends not only on the functional choice of $f$, but also on the prescription for $L_m$. Across the studies summarized here, both $L_m=-p$ and $L_m=\rho$ are used, with materially different consequences for conservation laws and effective dynamics [2206.09278; 2406.08303].

## 2. Energy–momentum balance, extra force, and the role of $L_m$

A defining property of generic $f(R,L_m)$ models is the modified covariant balance equation
$$
\nabla^\mu T_{\mu\nu}=
\bigl(g_{\mu\nu}L_m-T_{\mu\nu}\bigr)\nabla^\mu\ln f_{L_m},
$$
or, equivalently in another common presentation,
$$
\nabla^\mu T_{\mu\nu}
=2\nabla^\mu\ln(f_{L_m})\,\frac{\partial L_m}{\partial g^{\mu\nu}}.
$$
Hence, unless $f_{L_m}$ is constant or a special matter prescription is chosen, matter and geometry exchange energy–momentum [1008.4193; 2308.06519].

For a perfect fluid with $L_m=L_m(\rho)$, the non-conservation law implies non-geodesic motion. The equation of motion may be written as
$$
u^\mu\nabla_\mu u^\nu=f^\nu,
$$
with an extra force $f^\nu$ orthogonal to the four-velocity, $f^\nu u_\nu=0$ [1008.4193]. In the weak-field, slow-motion limit, this becomes an Euler-type equation with an additional acceleration
$$
a_E^i=-\partial^i\Bigl[U(R,L_m)\,\frac{dL_m}{d\rho}\Bigr],
$$
beyond the Newtonian gravitational acceleration and the ordinary pressure-gradient term [1008.4193]. A related Newtonian-limit treatment gives an effective Poisson relation with $G_{\rm eff}=f_L/(2f_R)$ under the assumptions used there, together with an extra acceleration $a_E=-\nabla U$ [2106.10644].

A common misconception is that $\nabla^\mu T_{\mu\nu}\neq0$ is unavoidable in every concrete realization. The white-dwarf study based on
$$
f(R,L_m)=\frac{R}{2}+L_m+\sigma R L_m
$$
instead chooses $L_m=-p$, specifically so that $\nabla^\mu T_{\mu\nu}=0$ remains valid on shell. In that setting, the modified dynamics enter through altered field equations and hydrostatic balance rather than through explicit matter non-conservation [2206.09278]. This suggests that the choice of matter Lagrangian is not a secondary convention but a central structural ingredient of the theory.

## 3. FLRW cosmology, autonomous systems, and late-time acceleration

In a spatially flat FLRW background,
$$
ds^2=-dt^2+a^2(t)\,d\vec x^2,\qquad H=\dot a/a,
$$
the modified Friedmann equations can be written as
$$
3H^2f_R+\tfrac12\bigl(f-f_RR-f_{L_m}L_m\bigr)+3H\dot f_R
=\tfrac12 f_{L_m}\rho,
$$
$$
\dot H f_R+3H^2f_R-\ddot f_R-3H\dot f_R
+\tfrac12\bigl(f_{L_m}L_m-f\bigr)
=\tfrac12 f_{L_m}p,
$$
with model-dependent matter evolution [2406.08303].

Several late-time cosmological realizations have been analyzed. For
$$
f(R,L_m)=\Lambda+\frac{\alpha}{2}R+\beta L_m^n,\qquad L_m=-p,
$$
a dynamical-systems reformulation yields three critical points: a matter-dominated decelerating point $A$ with $a\propto t^{2/3}$, $q=1/2$, and $w_{\rm eff}=0$; a de Sitter attractor $B$ with $q=-1$ and $w_{\rm eff}=-1$; and a scaling solution $C$ whose stability depends on parameter ranges. In that analysis, $B$ is stable for
$$
0<n<\frac{2+\omega}{1+\omega},
$$
and the phase-plane trajectories begin near $A$, pass close to $C$, and end in $B$ [2308.06519].

A broader nonlinear matter class,
$$
f(R,L_m)=\frac{R}{2}+c_1L_m+c_nL_m^n+c_0,
$$
has been split into two cosmological cases with $L_m=\rho_m$ and an uncoupled radiation sector. In Case A,
$$
f(R,L_m)=\frac{R}{2}+\beta\rho_m^n+\gamma,
$$
the standard radiation $\to$ matter $\to$ de Sitter sequence occurs only for $n\gtrsim4/5$, with acceleration essentially enforced by the vacuum term. In Case B,
$$
f(R,L_m)=\frac{R}{2}+\beta\rho_m+\gamma\rho_m^n,
$$
the interval $0<n<1/2$ produces a physical scaling de Sitter future attractor inside the bounded simplex, with $q=-1$ and $\omega_{\rm eff}=-1$ and without introducing $c_0$ [2601.10699].

Other background constructions emphasize phenomenological reconstruction rather than phase-space completeness. A transit-dark-energy model based on
$$
f(R,L_m)=\frac{R}{2}+\alpha L_m^n-\beta
$$
finds that as $t\to\infty$ (or $z\to0$) the cosmographic set approaches
$$
\{q,j,s,l,m\}\to\{-1,1,1,1,1\},
$$
with a transition redshift reported in the range $z_t\approx0.45\text{–}0.70$ depending on the dataset [2209.14269]. A constrained quintessence model with
$$
f(R,L_m)=\frac{1}{2}R+L_m^n+\zeta
$$
and a parametrized deceleration function yields late-time quintessence behavior, with joint best-fit values
$$
n=0.4999\pm0.0001,\quad
\alpha=1.3019\pm0.0001,\quad
\beta=2.7284\pm0.0001,\quad
t_0=13.1006\pm0.0001\ {\rm Gyr}
$$
for the four parameters of that construction [2212.12321]. A separate freezing-quintessence analysis of
$$
f(R,L_m)=\frac{R}{2}+L_m^\alpha
$$
with a sinh ansatz for $a(t)$ reports that the model approaches a de Sitter epoch identical to $\Lambda$CDM as $z\to-1$, and for the quoted best-fit values the adiabatic sound speed remains in $0<\nu_s^2<1$ for all $z\ge0$ [2412.10518].

## 4. Perturbations, growth observables, and observational status

Beyond the background, $f(R,L_m)$ gravity modifies the evolution of matter perturbations through altered effective couplings and matter-dilution laws. In the model
$$
f(R,L_m)=\lambda R+\beta L_m^\alpha+\eta,
$$
the background was constrained using 57 OHD points, 1048 SNIa distance moduli, and their combination, while perturbations were analyzed with 14 growth-rate data points and 30 $f\sigma_8$ measurements. The reported OHD+SNIa best fit is
$$
\alpha=1.091^{+0.035}_{-0.042},\quad
\beta=1.237^{+0.056}_{-0.16},\quad
\lambda=0.630^{+0.031}_{-0.050},\quad
\Omega_m=0.287\pm0.031,\quad
H_0=71.72^{+0.26}_{-0.23},
$$
while the $f\sigma_8$ fit gives
$$
\Omega_m=0.284^{+0.035}_{-0.049},\quad
\sigma_8=0.799^{+0.045}_{-0.086},\quad
\alpha=0.766^{+0.026}_{-0.064},\quad
\beta=1.08^{+0.40}_{-0.16},\quad
\lambda=0.279^{+0.078}_{-0.11}.
$$
Using AIC and BIC, that study finds substantial support from OHD alone but weaker support once SNIa or combined datasets are included [2406.08303].

A more direct growth-rate study considered
$$
f(R,L_m)=\alpha R+L_m^\beta+\gamma,\qquad L_m=\rho,
$$
derived a modified Poisson law with
$$
G_{\rm eff}(z)=\frac{\beta(2\beta-1)\rho^{\beta-1}(z)}{8\pi\alpha},
$$
and used 23 $f\sigma_8(z)$ measurements over $0.02\le z\le2.6$. The quoted best-fit values are
$$
H_0=73.75\pm0.16\ {\rm km\,s^{-1}Mpc^{-1}},\qquad
\lambda=0.262\pm0.007,\qquad
w=-0.005\pm0.001,
$$
with equivalent parameters $\beta=1.00505$, $\alpha\simeq4.51\times10^5$, and $\gamma\simeq-1.14\times10^{-29}$. Relative to $\Lambda$CDM, the model produces small $(\sim1\%\text{–}5\%)$ deviations in $f\sigma_8$ at intermediate and high redshifts, a difference of $\sim0.05$ near $z\sim0.8$, and up to $\sim10\%$ suppression of growth at $z\gtrsim1.5$ [2507.04079].

At the purely background level, multiple fits cluster around $\Lambda$CDM-like kinematics but with model-dependent interpretations. The dynamical-systems study reports
$$
H_0=67.2^{+3.5}_{-4.0}\ {\rm km\,s^{-1}Mpc^{-1}},\quad
n=-1.24\pm0.15,\quad
\alpha=-0.059^{+0.021}_{-0.034},\quad
\beta=0.050^{+0.024}_{-0.030},\quad
\Omega_{m0}=0.25^{+0.13}_{-0.15},\quad
\omega=-1.278\pm0.034,
$$
together with $\chi^2_{\min}=16.49$, and interprets the resulting $w_{\rm eff}< -1$ as phantom dark energy [2308.06519]. By contrast, the 2026 autonomous-systems analysis finds Case A consistent with $\Lambda$CDM and Case B within the accelerating window, both at the background level, with $n=1.08\pm0.05$ and $n=0.05\pm0.10$, respectively [2601.10699]. A plausible implication is that background observables alone do not isolate a unique $f(R,L_m)$ phenomenology; perturbative probes are structurally more discriminating.

## 5. Compact stars, wormholes, and other strong-field realizations

The compact-object literature shows that curvature–matter couplings can be important in regimes where pressure gradients or effective exotic stresses are large. In the white-dwarf application of
$$
f(R,L_m)=\frac{R}{2}+L_m+\sigma R L_m,
$$
the static spherically symmetric system was recast into TOV-like equations using the on-shell choice $L_m=-p$ and the zero-temperature Hamada–Salpeter equation of state. The resulting equilibrium sequences predict maximum masses above the Chandrasekhar mass limit. The most important reported effect is a significant increase of the mass for stars with radius $<2000$ km, while GR is recovered for stars with radii larger than 3000 km independently of $\sigma$. For a helium core, the sequence goes from $M_{\max}\approx1.40\,M_\odot$ at $\sigma=0$ to $M_{\max}\gtrsim2.0\,M_\odot$ at $\sigma=0.5$, and the quoted white-dwarf bound is
$$
0<\sigma\lesssim1\ {\rm km}^2
$$
[2206.09278].

Wormhole studies have largely focused on nonlinear matter couplings of the form
$$
f(R,L_m)=\frac{R}{2}+L_m^\alpha.
$$
With non-commutative Gaussian and Lorentzian smearing, tideless Morris–Thorne wormholes are obtained that are traversable, horizon-free, and asymptotically flat. In both distributions, the reported energy-condition pattern is qualitatively the same: $\rho>0$, the radial NEC is violated, the tangential NEC holds, the SEC holds, and both radial and tangential DEC are violated. A stability analysis based on the TOV equation and the sound-speed criterion yields
$$
v_s^2=\frac{2}{3\alpha}-1,\qquad
\frac13\le \alpha <\frac23,
$$
as the dynamically stable range for the solutions discussed there [2307.02498].

A charged Casimir generalization in the same gravity sector adds a GUP correction and an electromagnetic contribution through
$$
f(R,L_m)=\frac{R}{2\kappa}+L_m^\alpha.
$$
That analysis reports that electric charges, GUP effects, and higher model parameter values increase the throat length, and that the NEC remains violated despite a positive contribution from the electromagnetic energy density. The deflection angle was studied via the Gauss–Bonnet theorem, with explicit charge- and GUP-dependent corrections in the large-impact-parameter expansion [2408.03969].

Wormholes sourced by galactic dark-matter profiles have also been constructed in both
$$
f(R,L_m)=\frac{R}{2}+L_m^\alpha
$$
and
$$
f(R,L_m)=\frac{R}{2}+(1+\lambda R)L_m.
$$
Using URC, NFW model-I, and NFW model-II halos, the reported solutions satisfy the standard wormhole requirements with appropriate parameter choices, while violating the null energy constraints in the radial sector. In that sense, dark-matter halos are found to support wormholes within the effective stress structure of these models [2403.17037].

## 6. Early-universe sectors, extensions, and unresolved issues

The framework has also been used in early-universe model building. In gravitational baryogenesis with
$$
f(R,L_m)=\frac{R}{2}+L_m^\alpha+\zeta,
$$
the effective CP-violating interaction proportional to $\partial_\mu f(R,L_m)$ generates a non-zero baryon-to-entropy ratio during radiation dominance. For the numerical choice
$$
g_{*s}=106,\quad g_B=1,\quad T_D=2\times10^{12}\ {\rm GeV},\quad M_*=2\times10^{16}\ {\rm GeV},
$$
the paper reports
$$
\frac{n_B}{s}\simeq7.29\times10^{-11}
$$
for $\alpha=0.79$ and $\zeta=2$, and states that acceptable values lie roughly in
$$
0.75\lesssim\alpha\lesssim0.85\qquad (\zeta=2)
$$
[2304.02482].

A bouncing cosmology based on
$$
f(R,L_m)=\frac{R}{2}+L_m^\gamma
$$
finds a non-singular bounce for the ansatz $a(t)=(\alpha+\beta t^2)^{1/n}$, with $H=0$ at the bounce point and an equation-of-state parameter that crosses the phantom divide line $\omega=-1$ when $\dot H=0$, i.e. at $t^2=\alpha/\beta$. In that treatment, a successful bounce requires NEC violation near $t=0$, the SEC is likewise violated around the bounce, and for $\gamma>1$ the model can satisfy $0\le c_s^2\le1$ in the immediate vicinity of the bounce, although late-time ghost-like behavior is reported as a possible instability [2403.15772].

The direct generalization to $f(R,L_m,T)$ gravity retains $f(R,L_m)$ as a special case and extends the explicit coupling to the trace $T$ of the energy–momentum tensor. In that larger framework, the Newtonian limit, generalized Poisson equation, and Dolgov–Kawasaki instability were re-examined, while the $f(R,L_m)$ sector continues to exhibit non-geodesic motion and matter–geometry exchange in generic circumstances [2106.10644].

Several open issues recur across the literature. One is the physical status of the matter Lagrangian: $L_m=-p$ and $L_m=\rho$ are both employed, and the choice can determine whether $\nabla^\mu T_{\mu\nu}$ vanishes or not in a given model [2206.09278; 2406.08303]. Another is stability: the original formulation explicitly notes the need to reanalyze ghost modes, Dolgov–Kawasaki instability, and singularity formation in the generalized setting [1008.4193]. A further unresolved direction is perturbation-level phenomenology. The 2026 phase-space analysis identifies tensor no-ghost and luminal propagation as minimal consistency conditions, but also argues that growth, lensing, and scalar-sector stability must be studied explicitly through observables such as $f\sigma_8$ and weak lensing before the viability of nonlinear matter branches can be decisively assessed [2601.10699].

In aggregate, $f(R,L_m)$ gravity is best understood not as a single model but as a formal umbrella for curvature–matter couplings whose phenomenology depends crucially on the choice of $f(R,L_m)$, the prescription for $L_m$, and the physical regime under study. The supplied literature shows that these theories can reproduce GR in controlled limits, generate late-time acceleration, alter structure growth, support super-Chandrasekhar white dwarfs, and sustain wormhole geometries, while simultaneously raising unresolved questions about conservation laws, stability, and precision observational discrimination.

Source: https://www.emergentmind.com/topics/f-r-l_m-gravity