Four-Velocity Normalization Resolves Matter Lagrangian Ambiguity in Modified Gravity

Theoretical physicists B. N. Jayawiguna and A. Sulaksono have resolved a mathematical ambiguity in modified gravity theories by introducing a consistent variational procedure that accounts for the normalization condition of fluid four-velocity. In general relativity, using either isotropic pressure or negative energy density as the matter Lagrangian density yields identical field equations and preserves thermodynamic consistency. However, in contemporary modified gravity theories featuring nonminimal couplings between spacetime geometry and matter, the resulting field equations depend directly on the second-order derivative of the matter Lagrangian density with respect to the metric tensor. Consequently, selecting pressure versus negative density historically produced divergent field equations and contradictory physical predictions across astrophysical and cosmological contexts. By demonstrating that fluid four-velocity and the metric tensor cannot be varied as independent degrees of freedom, the authors established a unified framework wherein field equations become strictly invariant under the choice of matter Lagrangian.

The theoretical necessity for modifying general relativity arises from persistent limitations encountered in both strong-gravity ultraviolet and large-distance infrared regimes. While general relativity remains exceptionally successful across solar system tests, it struggles to account for cosmic late-time acceleration without introducing an ad hoc cosmological constant that suffers from fine-tuning and coincidence dilemmas. To resolve this cosmic expansion puzzle, theorists frequently explore modified gravitational actions such as f(R) gravity models, which incorporate inverse-curvature corrections that dominate as cosmic curvature diminishes over time. Simultaneously, in the ultraviolet domain, general relativity predicts unphysical spacetime singularities inside black holes and fails to account for observations of super-Chandrasekhar white dwarfs possessing masses between 2.1 and 2.8 solar masses in overluminous Type Ia supernovae, exceeding the classical 1.4 solar mass Chandrasekhar limit. These observational challenges motivated theories that couple geometry directly to matter, but such couplings produce higher-order metric variations requiring precise mathematical treatment.

To establish the breadth of this ambiguity, the researchers systematically analyzed four prominent families of curvature-matter modified gravity: f(R, T) gravity, energy-momentum squared gravity denoted as f(R, tau), f(R, T, P) gravity, and the generalized coupling model f(R, TG, TGD). In f(R, T) gravity, the action combines the Ricci scalar with an arbitrary function of the energy-momentum trace, generating field equations governed by an auxiliary matter tensor Theta_mu_nu containing second metric derivatives of the matter Lagrangian. In energy-momentum squared gravity, the action incorporates a scalar contraction of squared energy-momentum tensors, producing an auxiliary tensor theta_mu_nu that similarly encapsulates second-order metric variations. In f(R, T, P) theories, the coupling introduces a scalar invariant contracting the Ricci tensor directly with the energy-momentum tensor, yielding a correction tensor Gamma_mu_nu. Finally, f(R, TG, TGD) gravity incorporates contractions between the Einstein tensor and the energy-momentum tensor, producing auxiliary tensors Theta_mu_nu and Xi_mu_nu.

The root of this divergence stems from how fluid variables were traditionally handled during variational derivation. Assuming a standard isotropic perfect fluid defined by energy density rho, pressure p, and four-velocity u_mu, the first derivative of pressure with respect to the inverse metric tensor produces a four-velocity dyad, while the first derivative of density contains both four-velocity and metric tensor components. Under the conventional assumption that the four-velocity and the metric tensor represent independent bases, the second derivative of pressure vanishes completely because pressure lacks explicit metric dependence. In stark contrast, the second derivative of energy density retains metric dependence, generating a non-zero tensor proportional to metric tensor products. This discrepancy meant that setting the matter Lagrangian to pressure versus negative density produced unequal auxiliary tensors. Previous literature handled this inconsistently, with some authors dropping the second derivative of density entirely while others retained divergent field equations that yielded conflicting physical predictions.

To overcome this longstanding inconsistency, Jayawiguna and Sulaksono reformulated the variational process by incorporating the four-velocity normalization condition directly into the derivatives. Utilizing the relation where the variation of the four-velocity dyad equals the fourfold product of four-velocities, the researchers computed the second metric derivatives without relying on an isentropic chain rule or introducing an adiabatic sound speed. Under this direct closure, the second derivative of pressure becomes non-zero and incorporates metric and four-velocity contractions, while the second derivative of density is expressed as the negative second derivative of pressure augmented by symmetric metric products. When these revised second-order derivatives are substituted into the field equations of f(R, T), f(R, tau), f(R, T, P), and f(R, TG, TGD) gravity, the resulting auxiliary tensors evaluate to identical forms regardless of whether one chooses pressure or negative density as the matter Lagrangian, restoring foundational invariance.

To assess the physical impact of their revised prescription in the ultraviolet strong-gravity regime, the authors examined relativistic stellar equilibrium in compact stars using the f(R, T) = R + 2 chi T model with chi equal to one in geometrized units where G and c equal one. Adopting a standard spherically symmetric metric governed by radial metric functions nu(r) and lambda(r), the revised formalism establishes that the auxiliary tensor Theta_mu_nu exactly cancels against the matter tensor, vanishing from the right-hand side of the modified field equations. This algebraic cancellation simplifies the gravitational field equations into an Einstein tensor coupled strictly to the energy-momentum tensor and a trace term. From this simplified structure, the authors derived explicit modified Tolman-Oppenheimer-Volkoff hydrostatic equilibrium equations for the enclosed stellar mass m’, the metric derivative nu’, and the pressure gradient p’, providing a clean description of stellar balance that naturally recovers general relativity when chi approaches zero.

Applying these hydrostatic equilibrium equations to realistic equations of state revealed fundamentally different physical responses across distinct classes of compact stars. The authors modeled quark stars using the standard MIT bag model describing up, down, and strange quarks with a bag constant B of 60 MeV/fm^3, while modeling neutron stars using the realistic, density-dependent MPA1 equation of state. Because the MIT bag model possesses an approximately linear equation of state relating energy density and pressure, quark stars exhibit minimal sensitivity to the chi coupling corrections, producing mass-radius trajectories that align remarkably closely with standard general relativity while slightly shifting the maximum mass. In contrast, the MPA1 equation of state for neutron stars depends strongly on central density, causing the modified gravity coupling to produce substantial deviations in both stellar radius and maximum mass at high central densities, exceeding general relativity predictions. Both stellar profiles demonstrated robust dynamical stability up to their respective maximum masses.

Turning to the infrared cosmological domain, the researchers analyzed cosmic evolution within energy-momentum squared gravity across a spatially flat Friedmann-Lemaître-Robertson-Walker metric with scale factor a(t) and zero cosmological constant. In this revised framework, the auxiliary tensor theta_mu_nu evaluates identically to zero, modifying cosmological dynamics strictly through effective energy density and pressure without introducing extraneous source terms. Assuming a barotropic equation of state p = w rho, the authors derived the unified trio of cosmological equations across matter-dominated dust where w = 0, radiation where w = 1/3, and vacuum energy where w = -1. A crucial physical consequence of this derivation is that the quadratic energy density correction governed by coupling constant alpha undergoes an essential sign reversal compared to the traditional prescription. Under the revised formulation, the modified Friedmann equation incorporates a negative quadratic correction proportional to alpha rho squared, ensuring that deviations from general relativity become prominent at ultra-high densities.

This structural sign reversal produces profound differences in early-universe cosmology, most notably during the radiation-dominated era where energy densities reach extreme magnitudes. In the revised radiation field equations, the quadratic correction appears with a negative sign, establishing a rigorous upper bound on radiation energy density equal to 3/(8 alpha) for positive values of alpha, attained when the Hubble rate squared reaches pi/(2 alpha). This finding directly contrasts with previous energy-momentum squared gravity literature relying on the traditional pressure Lagrangian, which required a negative coupling parameter to obtain an unphysical bound of -1/(8 alpha). Evaluating cosmic expansion with a coupling parameter of 10^-37 cm^3/erg, the new prescription reveals that radiation density attained a finite maximum value in the early universe between 10^-5 and 10^-4 seconds rather than diverging into an initial singularity. Furthermore, analytical expressions for cosmic time demonstrate that the cosmological trajectories remain remarkably close to standard general relativity throughout subsequent cosmic expansion.

By demonstrating that four-velocity normalization restores matter Lagrangian invariance across diverse modified gravity theories, this study establishes a consistent mathematical framework for exploring nonminimal matter-geometry couplings. Nevertheless, the authors emphasize critical theoretical limitations and observational distinctions that define necessary boundaries for interpretation. Because previous observational constraints from Big Bang nucleosynthesis were derived exclusively under the traditional assumption of a pressure Lagrangian, those historical parameter limits do not directly apply to the revised cosmological equations and must be reassessed under the new sign structure. In addition, the direct sum closure employed in this work results in second metric derivatives that lack manifest index exchange symmetry across pairs of metric indices, presenting a technical contrast with earlier thermodynamic variational approaches that incorporated adiabatic sound speed. Overall, this universal variational formulation clarifies longstanding ambiguities in the matter sector, offering a reliable theoretical foundation for modeling extreme matter in neutron star interiors and formulating dynamical equations for infrared cosmology.

Subject of Research: Resolution of matter Lagrangian ambiguity in modified gravity theories through four-velocity normalization

Article Title: Ambiguity in matter sector for modified gravity involving \(\delta ^2 \mathcal {L}_{m}/\delta g^{\mu \nu }\delta g^{\alpha \beta }\) and its implications to astrophysics and cosmology

Article References: Jayawiguna, B. N., & Sulaksono, A. (2026). Ambiguity in matter sector for modified gravity involving $$\delta ^2 \mathcal {L}_{m}/\delta g^{\mu \nu }\delta g^{\alpha \beta }$$ and its implications to astrophysics and cosmology. The European Physical Journal C, 86(11), Article 1254. https://doi.org/10.1140/epjc/s10052-026-16443-z

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16443-z

Keywords: Modified gravity, General relativity, Matter Lagrangian density, Four-velocity normalization, Energy-momentum squared gravity, Neutron stars, Quark stars, MIT bag model, MPA1 equation of state, Cosmological expansion, Friedmann equations, Radiation-dominated universe, Ricci scalar, Tolman-Oppenheimer-Volkoff equations, Cosmic acceleration, Scalar-tensor theories, FLRW metric, Barotropic equation of state