⚛️ general relativity

Lagrangian perturbation theory for modified gravity

This paper presents a third-order Lagrangian perturbation theory formalism for modified gravity models with scale-dependent gravitational strengths, highlighting the necessity of frame-lagging terms to recover the $\Lambda$ CDM limit and applying the framework to compute 1-loop power spectra and correlation functions in $f(R)$ gravity.

Original authors: Alejandro Aviles, Jorge L. Cervantes-Cota

Published 2026-10-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Alejandro Aviles, Jorge L. Cervantes-Cota

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Gravity is the invisible architect of our universe, the force that pulls dust together to form stars and shepherds galaxies into vast clusters. For over a century, Albert Einstein's theory of General Relativity has been the gold standard for describing how this force works, passing every test scientists have thrown at it, from the motion of planets to the bending of light. Yet, when astronomers look at the universe on its grandest scales, they see something strange: the cosmos is not just expanding, but accelerating. To explain this, the standard model of cosmology includes a mysterious component called dark energy, often represented by a cosmological constant. However, this solution comes with its own puzzles, leading some researchers to wonder if gravity itself might behave differently over cosmic distances, perhaps requiring a modification to Einstein's laws rather than a new form of energy.

The challenge in testing these ideas is that gravity is notoriously difficult to study in the messy, crowded environment of the real universe. While we can simulate the smooth expansion of the cosmos, the formation of structures like galaxies involves complex, non-linear interactions where matter clumps together in ways that are hard to predict with simple math. To understand how these structures form, scientists rely on theories that describe how tiny fluctuations in the early universe grew into the massive web of matter we see today. One powerful approach is to track the journey of individual particles of dark matter as they move from their starting positions to their current locations. This method, known as Lagrangian perturbation theory, offers a clearer picture of how these particles drift and cluster, especially when trying to distinguish between standard gravity and alternative theories that propose a fifth force acting alongside gravity.

In a recent study, researchers Alejandro Aviles and Jorge L. Cervantes-Cota developed a new mathematical framework to track these particle movements within the context of modified gravity theories. Their work focuses on a specific class of theories where the strength of gravity changes depending on the scale, a feature designed to hide the effects of this extra force in our local solar system while allowing it to influence the universe at large. These theories often involve a scalar field, a kind of invisible energy field that mediates a fifth force, which must be "screened" or suppressed in dense regions to avoid contradicting local experiments. The authors aimed to create a tool that could calculate how matter moves and clusters under these conditions, pushing their calculations to a high level of precision to see how well these theories hold up against the standard model.

The team discovered that when they translated the equations governing this scalar field from a stationary viewpoint to a moving one that follows the particles, a new and previously overlooked effect emerged. They found that the transformation itself generates terms that act like a lag or a delay in the reference frame. If these terms are ignored, the theory fails to revert to the standard model of gravity on the largest scales, where the fifth force should naturally disappear. By including these frame-lagging terms, the researchers showed that their modified gravity model correctly smooths out to match the standard cosmological model at vast distances, ensuring consistency with what we observe in the local universe. This correction was essential to making the theory work as intended, acting as a necessary bridge between the complex, non-linear behavior of the scalar field and the familiar laws of gravity.

Using this refined formalism, the authors applied their method to a popular modified gravity model known as Hu-Sawicki f(R) gravity, which is designed to mimic the acceleration of the universe without dark energy. They calculated the statistical patterns of how matter is distributed, specifically looking at the power spectrum, which measures how much structure exists at different sizes, and the correlation function, which describes how likely galaxies are to be found near one another. They tested three variations of this model, labeled F4, F5, and F6, which represent different strengths of the deviation from standard gravity. The results showed that while these models produce distinct signatures in the way matter clusters, particularly at intermediate scales, they must carefully balance their non-linear effects to avoid breaking the rules of gravity in our local neighborhood.

A key finding of the study was that the behavior of these modified theories changes significantly depending on the scale of the observation. On small scales, where matter is dense, the screening mechanism effectively hides the fifth force, making the universe look like it is governed by standard gravity. However, on larger scales, the force becomes active, altering the growth of cosmic structures. The researchers demonstrated that their new approach could accurately capture these transitions, providing a more reliable way to predict the distribution of matter than previous methods. They also addressed potential mathematical problems that arise when calculating these statistics, showing that their method remains stable and does not produce infinite or nonsensical results, a common issue in complex gravitational calculations.

The study concludes that while modified gravity offers a compelling alternative to dark energy, any viable theory must include these specific frame-lagging corrections to be consistent with observations. Without them, the theory would predict a universe that looks different from what we see, even on the largest scales. By providing a robust mathematical tool to compute these effects up to the third order of precision, Aviles and Cervantes-Cota have given cosmologists a sharper instrument to test these ideas against the flood of data coming from modern sky surveys. Their work suggests that the path to understanding the accelerating universe may lie not in adding new ingredients to the cosmic recipe, but in refining our understanding of how gravity itself operates across the vast expanse of space and time.

Technical Summary: Lagrangian Perturbation Theory for Modified Gravity

Problem Statement
General Relativity (GR) remains the most successful theory of gravity, yet the discovery of the universe's accelerated expansion has spurred interest in Modified Gravity (MG) theories as alternatives to the cosmological constant ( $\Lambda$ CDM). A primary challenge in MG is constructing theories that satisfy local gravitational constraints while driving cosmic acceleration. This is typically achieved via screening mechanisms (e.g., chameleon, Vainshtein, kinetic) that suppress the "fifth force" mediated by an additional scalar degree of freedom in high-density regions, recovering GR locally.

While $N$ -body simulations are the standard for studying structure formation in MG, they are computationally expensive and obscure the physical mechanisms of structure growth. Analytical methods, specifically Perturbation Theory (PT), offer a complementary approach. However, existing PT formalisms for MG have largely relied on Eulerian Standard Perturbation Theory (SPT). While SPT tracks the broadband power spectrum well, it struggles to reproduce the Baryon Acoustic Oscillation (BAO) features and the correlation function around the acoustic peak. Conversely, Lagrangian Perturbation Theory (LPT) excels at modeling BAO and correlation functions but has been less developed for MG, particularly regarding the non-linearities introduced by screening mechanisms and the transformation of field equations between coordinate systems.

Methodology
The authors develop a comprehensive LPT formalism capable of computing Lagrangian displacement fields up to third order for a wide range of MG theories, specifically those that can be formulated as scalar-tensor gravity. The core of the methodology involves:

Coordinate Transformation: The authors transform the equations of motion from Eulerian coordinates ( $x$ ) to Lagrangian coordinates ( $q$ ), where $x = q + \Psi(q,t)$ . A critical technical step involves the transformation of the Klein-Gordon equation for the scalar field. This transformation introduces a "frame-lagging" term ( $\nabla_x^2 \phi - \nabla_q^2 \phi$ ), which arises from the Jacobian of the coordinate transformation. The authors argue that neglecting this term prevents the theory from correctly reducing to $\Lambda$ CDM at large scales.
Screening Mechanisms: The formalism incorporates the non-linearities of the scalar field potential (the screening mechanism) via an expansion of the self-interaction term $\delta I(\phi)$ . This allows the fifth force to be suppressed in high-density regions.
Order-by-Order Expansion: The Lagrangian displacement $\Psi$ $Ψ$ and the scalar field $\phi$ $ϕ$ are expanded in powers of a perturbation parameter. The authors derive explicit kernels and growth functions for:
- Linear Order (1LPT): Deriving the scale-dependent growth function $D_+(k,t)$ .
- Second Order (2LPT): Computing the second-order displacement kernel, which includes contributions from the standard gravitational interaction, the scalar field self-interaction, and the newly identified frame-lagging term.
- Third Order (3LPT): Extending the expansion to third order to enable the calculation of 1-loop power spectra and correlation functions.
Application to $f(R)$ Gravity: The formalism is explicitly applied to the Hu-Sawicki $f(R)$ model. The authors derive the specific mass functions ( $M_1, M_2, M_3$ ) for this model and compute the resulting growth functions for different triangular configurations (equilateral, orthogonal, squeezed) in Fourier space.

Key Contributions and Results

Identification of Frame-Lagging Terms: The paper identifies that transforming the Klein-Gordon equation to Lagrangian coordinates generates non-negligible "frame-lagging" terms. The authors demonstrate that these terms are essential for the theory to recover GR (specifically the $\Lambda$ CDM limit) at large scales where the fifth force is screened. Without these terms, the theory fails to reduce to the standard model in the appropriate limit.
Resolution of Discrepancies: The authors compare their 2LPT results with previous works (specifically Refs. [31] and [32]). They find that their results align with Ref. [32] but differ from Ref. [31], attributing the discrepancy to the omission of frame-lagging terms and incorrect assumptions about the equality of certain growth function components in the previous literature.
Screening and Large Scale Behavior: The analysis shows that in the "squeezed" limit (large scales, $k \to 0$ ), the combination of growth functions including the frame-lagging term cancels out to zero, recovering the GR result for planar collapse. This confirms that the Zel'dovich approximation remains exact for 1D collapse even in MG, provided the frame-lagging terms are included.
Power Spectrum and Correlation Function: The formalism is used to compute the 1-loop power spectrum and correlation function for $f(R)$ gravity using different resummation schemes. The authors discuss the Infrared (IR) divergences inherent in these formalisms, noting that for theories preserving the equivalence principle, the matter and displacement statistics should be IR-safe.
Model Specifics: For the Hu-Sawicki model (specifically the F4, F5, and F6 variants), the authors provide explicit expressions for the growth functions. They show that while the linear growth function depends on the scale $k$ , the higher-order growth functions depend on the interaction between modes of different wavelengths.

Significance and Claims
The paper claims to provide a robust, general formalism for applying Lagrangian Perturbation Theory to modified gravity theories, particularly those with scale-dependent gravitational strengths. The primary significance lies in:

Correcting the Large-Scale Limit: By rigorously including the frame-lagging terms arising from the coordinate transformation, the formalism ensures that MG theories correctly reduce to $\Lambda$ CDM at large scales, a feature previously overlooked or mishandled in some MG-LPT implementations.
Improved Modeling of Large-Scale Structure: The LPT approach is highlighted as superior to SPT for modeling BAO features and the correlation function, which are critical for interpreting upcoming cosmological survey data.
Clarification of Existing Literature: The work clarifies discrepancies between recent 2LPT derivations in MG, offering a unified framework that can be integrated into semi-analytical codes like COLA to efficiently simulate large-scale structure in MG scenarios.

The authors maintain a modest stance regarding the scope, noting that while the formalism is general, they specialize in $f(R)$ gravity for illustration. They also note that while the quasi-static approximation is used, its validity depends on specific model parameters and wavenumbers, and they do not explore dark energy perturbations in this specific work, though the framework allows for their inclusion.

More like this