f(R,L_m,T) Gravity: Unified Matter–Geometry Framework

• f(R,L_m,T) gravity is a unified theory where the gravitational Lagrangian depends on the Ricci scalar, matter Lagrangian, and the trace of the energy-momentum tensor.
• The theory unifies f(R,T) and f(R,L_m) models by introducing nonminimal couplings that modify the field equations and break energy-momentum conservation.
• Its applications span FLRW cosmology, compact objects, wormholes, and gravastars, offering new insights into cosmic acceleration and stellar structure.

Searching arXiv for core papers on $f(R,L_m,T)$ gravity and recent applications. $f(R,L_m,T)$ gravity is a generalized matter–geometry coupling theory in which the gravitational Lagrangian depends simultaneously on the Ricci scalar $R$, the matter Lagrangian density $L_m$, and the trace $T=g^{\mu\nu}T_{\mu\nu}$ of the matter energy-momentum tensor. In the notation of some applications, $L\equiv L_m$, so $f(R,L,T)$ and $f(R,L_m,T)$ denote the same framework. The theory was formulated as a unification of $f(R,T)$ and $f(R,L_m)$ gravity, and it has since been used as a common language for curvature–matter couplings in weak-field dynamics, FLRW cosmology, compact stars, wormholes, and gravastars (Haghani et al., 2021).

1. Definition, action, and theory space

The standard action used in the unified formulation is

$f(R,L_m,T)$0

with $f(R,L_m,T)$1 defined by variation of the matter action,

$f(R,L_m,T)$2

and, when $f(R,L_m,T)$3 depends on the metric but not on its derivatives,

$f(R,L_m,T)$4

This construction contains $f(R,L_m,T)$5, $f(R,L_m,T)$6, and $f(R,L_m,T)$7 as limiting cases, with GR recovered for $f(R,L_m,T)$8 (Haghani et al., 2021, Harko et al., 2010).

A structurally important point is that $f(R,L_m,T)$9 and $R$0 are not interchangeable variables. In $R$1 gravity, the trace can be written as

$R$2

so $R$3 depends not only on $R$4 but also on how $R$5 varies with the metric. This distinction remains central in $R$6: $R$7 is the matter Lagrangian variable entering the action, whereas $R$8 is a derived scalar built from the stress tensor (Lobo et al., 2012).

The literature uses more than one normalization convention. For example, the wormhole study in the notation $R$9 adopts

$L_m$0

in units $L_m$1, while still displaying $L_m$2 explicitly in the field equations. That paper states explicitly that $L_m$3 is used instead of the more standard $L_m$4 only for simplicity (Moraes et al., 2024).

2. Variational structure and field equations

With

$L_m$5

metric variation gives the field equations

$L_m$6

where

$L_m$7

Because $L_m$8 varies nontrivially with the metric, the formalism also introduces

$L_m$9

These tensors encode the genuinely new matter-response sector of the theory (Haghani et al., 2021).

The trace equation is

$T=g^{\mu\nu}T_{\mu\nu}$0

For weakly perturbed backgrounds in the class $T=g^{\mu\nu}T_{\mu\nu}$1, the linearized trace dynamics yields a Dolgov–Kawasaki criterion identical in form to the $T=g^{\mu\nu}T_{\mu\nu}$2 case,

$T=g^{\mu\nu}T_{\mu\nu}$3

so the sign of the second derivative with respect to $T=g^{\mu\nu}T_{\mu\nu}$4 remains the leading stability discriminator in that sector (Haghani et al., 2021).

Many phenomenological applications adopt algebraic ansätze with $T=g^{\mu\nu}T_{\mu\nu}$5. In that case the higher-derivative operator $T=g^{\mu\nu}T_{\mu\nu}$6 vanishes identically, and the modification relative to GR is carried entirely by explicit matter couplings. This is true for the commonly studied models $T=g^{\mu\nu}T_{\mu\nu}$7, $T=g^{\mu\nu}T_{\mu\nu}$8, and $T=g^{\mu\nu}T_{\mu\nu}$9 (Mota et al., 2024, Moraes et al., 2024).

3. Matter sector, non-conservation, and motion

A defining feature of $L\equiv L_m$0 gravity is that $L\equiv L_m$1 is generally not covariantly conserved. The unified formulation gives

$L\equiv L_m$2

with

$L\equiv L_m$3

Conservation is restored only when the gravitational Lagrangian has no explicit $L\equiv L_m$4- or $L\equiv L_m$5-dependence, namely $L\equiv L_m$6 (Haghani et al., 2021).

For perfect fluids,

$L\equiv L_m$7

the non-conservation law decomposes into modified energy and momentum balance equations, and test-particle motion becomes non-geodesic. The extra force can be written as

$L\equiv L_m$8

with $L\equiv L_m$9, so $f(R,L,T)$0 (Haghani et al., 2021).

The choice of matter Lagrangian is therefore physically consequential. In GR, the perfect-fluid prescriptions $f(R,L,T)$1 and $f(R,L,T)$2 are often interchangeable at the level of the stress tensor, but in nonminimally coupled theories they lead to inequivalent field equations, non-conservation laws, and stellar structure equations. This is emphasized both in compact-star applications and in the foundational $f(R,L,T)$3 literature on Lagrangian ambiguity (Mota et al., 2024, Otoniel et al., 24 Jul 2025, Carvalho et al., 2020).

Two criticisms sharpen this issue. First, the $f(R,L,T)$4 analysis of the Lagrangian ambiguity argues that once field equations depend explicitly on $f(R,L,T)$5, the choice of fluid Lagrangian becomes a physical input rather than a harmless convention. Second, the critique of separable $f(R,L,T)$6 concludes that $f(R,L,T)$7 should be absorbed into the matter Lagrangian and therefore has no independent gravitational significance. This suggests that in $f(R,L,T)$8, terms independent of curvature should be scrutinized for redefinition redundancy, whereas curvature–matter cross-couplings are the more robustly novel sector (Carvalho et al., 2020, Fisher et al., 2019).

4. Weak-field limit and cosmological dynamics

In the Newtonian, weak-field, slow-motion limit, the particle equation of motion becomes

$f(R,L,T)$9

with $f(R,L_m,T)$0. For the perturbative model $f(R,L_m,T)$1, $f(R,L_m,T)$2, and a linear barotropic equation of state $f(R,L_m,T)$3 with $f(R,L_m,T)$4, one finds

$f(R,L_m,T)$5

The generalized Poisson equation takes the form

$f(R,L_m,T)$6

with $f(R,L_m,T)$7 and $f(R,L_m,T)$8 determined by the background expansion of $f(R,L_m,T)$9 (Haghani et al., 2021).

FLRW cosmology has been studied both in the general unified theory and in concrete algebraic models. In the original unified treatment, the spatially flat metric

$f(R,T)$0

leads to modified Friedmann equations and a generalized acceleration condition

$f(R,T)$1

For the multiplicative model

$f(R,T)$2

the theory admits de Sitter solutions in the presence of ordinary matter only, both in radiation and dust eras, while an additive algebraic model

$f(R,T)$3

also supports exact FLRW solutions for radiation and dust (Haghani et al., 2021).

A later dynamical-systems study considered

$f(R,T)$4

and rewrote the flat-FLRW equations in terms of

$f(R,T)$5

with constraint $f(R,T)$6. The autonomous system has five critical points $f(R,T)$7, and the unique fully stable hyperbolic attractor is

$f(R,T)$8

corresponding to a de Sitter-like late-time state with $f(R,T)$9. The same study reports a transition from deceleration to acceleration and quotes $f(R,L_m)$0 and $f(R,L_m)$1 for the illustrated trajectories, while noting that the analysis remains background-level and does not include perturbations or detailed parameter estimation (Kshirsagar et al., 18 Nov 2025).

5. Compact objects and nontrivial spacetimes

Compact stars provide the most developed astrophysical applications. For isotropic neutron stars and quark stars in the model

$f(R,L_m)$2

the modified TOV equations depend strongly on whether one adopts $f(R,L_m)$3 or $f(R,L_m)$4. With realistic IU-FSU+BPS and MIT bag equations of state, negative $f(R,L_m)$5 generally favors heavier stars, while positive $f(R,L_m)$6 suppresses the maximum mass. The study emphasizes that the two Lagrangian prescriptions are appreciably different in the mass–radius diagrams and that, for the chosen EoSs, sufficiently negative $f(R,L_m)$7 can accommodate pulsars above $f(R,L_m)$8 even when the GR limit cannot (Mota et al., 2024).

The anisotropic quark-star analysis fixes $f(R,L_m)$9, uses the MIT bag model

$f(R,L_m,T)$00

and a quasi-local anisotropy relation

$f(R,L_m,T)$01

For $f(R,L_m,T)$02, $f(R,L_m,T)$03, negative $f(R,L_m,T)$04 increases the maximum mass and radius, whereas positive $f(R,L_m,T)$05 decreases them. In the tabulated benchmarks with $f(R,L_m,T)$06, the GR maximum mass $f(R,L_m,T)$07 rises to $f(R,L_m,T)$08 for $f(R,L_m,T)$09. Static stability, $f(R,L_m,T)$10, and $f(R,L_m,T)$11 are satisfied throughout the reported interiors (Tangphati et al., 2024).

White-dwarf structure has been examined in the same algebraic model $f(R,L_m,T)$12 for both $f(R,L_m,T)$13 and $f(R,L_m,T)$14, using a realistic zero-temperature EoS with a relativistic degenerate electron gas and ionic bcc lattice corrections. The coupling alters the mass–radius relation strongly at high central densities, permits stable super-Chandrasekhar configurations, and can remove the usual turning point $f(R,L_m,T)$15 in some branches. The sign of $f(R,L_m,T)$16 required to increase the maximum mass depends on the Lagrangian choice: for $f(R,L_m,T)$17, positive $f(R,L_m,T)$18 increases the mass, whereas for $f(R,L_m,T)$19, negative $f(R,L_m,T)$20 does so (Otoniel et al., 24 Jul 2025).

Nontrivial geometries have also been constructed. In the linear model

$f(R,L_m,T)$21

static Morris–Thorne wormholes with $f(R,L_m,T)$22 and barotropic radial equation of state $f(R,L_m,T)$23 admit the shape function

$f(R,L_m,T)$24

The flare-out condition still requires effective NEC violation, but the physical matter can satisfy $f(R,L_m,T)$25 and $f(R,L_m,T)$26 when $f(R,L_m,T)$27 and $f(R,L_m,T)$28. The central claim is therefore that ordinary matter can thread the wormhole while the effective exoticity is shifted to the modified-gravity sector (Moraes et al., 2024).

A gravastar model in $f(R,L_m,T)$29 with $f(R,L_m,T)$30 uses the Mazur–Mottola three-region decomposition $f(R,L_m,T)$31, $f(R,L_m,T)$32, and $f(R,L_m,T)$33. The interior remains regular and singularity-free, the shell density is expressed through a Lambert $f(R,L_m,T)$34-function,

$f(R,L_m,T)$35

and the exterior reduces to Schwarzschild because the matter-coupling modifications vanish in vacuum (Sinha et al., 2024).

6. Interpretation, reductions, and open problems

The existing $f(R,L_m,T)$36 literature is dominated by simple algebraic realizations in which $f(R,L_m,T)$37. This removes the higher-derivative sector characteristic of generic $f(R,L_m,T)$38 models and isolates the phenomenology of direct matter couplings. A plausible implication is that current applications map the consequences of nonminimal matter dependence more thoroughly than they map the full dynamical space of arbitrary $f(R,L_m,T)$39 functions (Moraes et al., 2024, Mota et al., 2024).

Three methodological issues recur. The first is the nontrivial status of the matter Lagrangian. Many analyses fix $f(R,L_m,T)$40, some compare it with $f(R,L_m,T)$41, and the resulting predictions are often qualitatively different. The second is that non-conservation is typically established, but its thermodynamic, microphysical, and observational consequences are only partially explored. The third is that many results are existence or background studies rather than complete viability analyses: wormhole work does not address dynamical stability or formation, cosmological dynamical-systems studies do not include perturbations or structure formation, and compact-star studies usually neglect rotation, magnetic fields, or finite temperature (Kshirsagar et al., 18 Nov 2025, Moraes et al., 2024, Tangphati et al., 2024).

A broader conceptual caution comes from adjacent $f(R,L_m,T)$42 and $f(R,L_m,T)$43 work. Explicit matter-sector dependence can make the theory sensitive to field redefinitions and to the arbitrary split between “gravity” and “matter.” This does not invalidate $f(R,L_m,T)$44 gravity, but it means that pure $f(R,L_m,T)$45- or $f(R,L_m,T)$46-dependence cannot automatically be treated as independent gravitational physics. By contrast, nonseparable curvature–matter couplings remain the least ambiguous source of genuinely new dynamics (Fisher et al., 2019, Carvalho et al., 2020).

In that sense, $f(R,L_m,T)$47 gravity is best characterized as a unified framework for simultaneous $f(R,L_m,T)$48-, $f(R,L_m,T)$49-, and $f(R,L_m,T)$50-dependence, with generically modified field equations, non-conserved $f(R,L_m,T)$51, non-geodesic motion, a modified Newtonian limit, and a strong dependence on the matter representation. Its phenomenology already spans FLRW acceleration, super-Chandrasekhar white dwarfs, heavy neutron and quark stars, gravastars, and wormholes supported by ordinary physical matter at the level of $f(R,L_m,T)$52. The main unresolved issue is no longer whether such effects can occur, but how much of them survive once the matter-sector ambiguities are fixed in a theoretically and observationally controlled way (Haghani et al., 2021).