Non-linear evolution in $f(R)$ gravity: perturbative modelling of the Chameleon mechanism

1. Problem Statement

The standard $\Lambda$CDM cosmological model faces theoretical and observational challenges, including the nature of dark energy and tensions in the Hubble constant ($H_0$) and matter fluctuation amplitude ($\sigma_8$). Modified gravity theories, specifically $f(R)$ gravity, offer an alternative by replacing the Ricci scalar $R$ in the Einstein-Hilbert action with a function $f(R)$.

However, viable $f(R)$ models must include a screening mechanism (the Chameleon mechanism) to recover General Relativity (GR) in high-density regions (like the Solar System) while allowing deviations on cosmological scales. This mechanism introduces a scalar degree of freedom where the Compton wavelength ($\lambda_C$) of the scalar field depends on the local matter density.

The Core Challenge:
Modelling the non-linear growth of structure in $f(R)$ gravity is computationally expensive. In standard GR, the gravitational potentials $\Phi$ and $\Psi$ are equal, allowing for linear solutions. In $f(R)$, $\Phi \neq \Psi$, and the difference $\chi = \Phi - \Psi$ satisfies a non-linear differential equation due to the density-dependent Compton scale. Solving this requires complex numerical solvers (e.g., multi-grid methods or N-body simulations with adaptive mesh refinement), which are slow and computationally intensive. Previous perturbative approaches often ignored the Chameleon mechanism, treating the scalar field equation as linear.

2. Methodology

The authors propose a hybrid Eulerian-Lagrangian scheme incorporating a perturbative solution for the non-linear scalar field equation.

• Model Framework: They utilize the Hu-Sawicki $f(R)$ model, a specific functional form designed to recover $\Lambda$CDM expansion history at early epochs while satisfying Solar System constraints.
• Perturbative Approach for $\chi$:
  • The authors define $\chi = \Phi - \Psi$. The governing equation for $\chi$ is non-linear because the Compton scale $\bar{x}_C$ depends on the Ricci scalar $R$, which in turn depends on the local density.
  • They assume an ansatz $\chi = \sum \epsilon^n \chi^{(n)}$, where $\epsilon$ is a book-keeping parameter related to the density contrast.
  • They expand the non-linear term $f_{RR}(R)$ around the background value $R_b$.
  • This transforms the non-linear problem into a hierarchy of linear Poisson-like equations for each order of $\chi^{(n)}$.
  • The solution is computed up to the fourth order in the perturbation expansion. This allows them to capture the non-linear feedback of the Chameleon mechanism without solving the full non-linear ODE numerically at every step.
• Evolution Scheme:
  • Step I (Eulerian): Solve the spatial equations for potentials $\Phi_+$ and $\chi$ assuming the density is constant over a time step.
  • Step II (Lagrangian): Evolve the density and velocity fields by tracking spherical shells (shells of constant initial comoving coordinate $q$) using the Euler and continuity equations.
  • The algorithm iterates between these steps, updating the density, velocity, and the local Compton scale at each time step.
• Initial Conditions: The study focuses on smooth, compensated top-hat perturbations of varying length scales ($x_{top}$) and smoothing parameters.

3. Key Contributions

1. Novel Perturbative Solver: The paper introduces a computationally efficient perturbative method to solve the non-linear Chameleon equation. Unlike previous works that relied on expensive relaxation methods or full N-body simulations, this approach uses analytic solutions to linear operators iteratively.
2. Spatio-Temporal Mapping of the Compton Scale: The authors explicitly track the evolution of the Compton scale ($\bar{x}_C$) as a function of both space and time, demonstrating how the screening mechanism dynamically alters the force range within a perturbation.
3. Identification of a New Feature: They report a density enhancement near the origin (center) of the perturbation, a feature not previously highlighted in literature focusing on spherical collapse.
4. Non-Unique DVDR: They demonstrate that in the presence of Chameleon screening, the Density-Velocity Divergence Relation (DVDR) is no longer a one-to-one mapping, even for a single spherical profile.

4. Key Results

The authors define three regimes of evolution based on the parameter $Q = \bar{x}_C / x_{top}$ (ratio of Compton scale to perturbation size):

• Strong-Field Regime ($Q \gg 1$):
  • Scenario: Perturbation size is much smaller than the Compton scale (e.g., $x_{top} = 0.002 h^{-1}$ Mpc).
  • Result: The Chameleon mechanism is ineffective. The scalar field equation reduces to a Poisson-like equation, and the effective gravitational constant is enhanced by a factor of $4/3$ ($G_{eff} = 4/3 G$). Growth is faster than GR, and the profile shape is retained.
• Weak-Field Regime ($Q \ll 1$):
  • Scenario: Perturbation size is much larger than the Compton scale (e.g., $x_{top} = 200 h^{-1}$ Mpc).
  • Result: The modification is suppressed. The dynamics closely mimic GR. Effects are only visible where density gradients are extremely sharp (smaller than $\lambda_C$), which is not the case for smooth profiles.
• Intermediate Regime ($Q \sim 1$):
  • Scenario: Perturbation size is comparable to the Compton scale (e.g., $x_{top} = 2 h^{-1}$ Mpc).
  • Result: This is the most astrophysically relevant regime.
    • Edge Enhancement: As observed in previous studies, a density peak forms near the outer edge of the top-hat due to mass shells piling up as they encounter the screened inner region.
    • Central Enhancement (New Finding): The authors find a secondary density bump near the origin.
      • Mechanism: As shells fall inward, the density near the origin decreases slightly relative to the edge. This lower density increases the local Compton scale ($\bar{x}_C$), reducing the screening effect and increasing the fifth-force locally. This accelerates inner shells, causing a second accumulation of mass near the center.
    • Growth Suppression: Unlike the strong-field case, the overall growth of the perturbation is suppressed compared to GR due to the screening effect dominating the inner regions.

Phase-Space Evolution (DVDR):

• In GR and $f(R)$ without screening, the relation between density contrast ($\Delta$) and velocity divergence ($\delta_v$) is a unique, monotonic curve (the Zeldovich curve).
• In the Chameleon-screened intermediate regime, the DVDR becomes multi-valued. A single density value can correspond to multiple velocity divergence values depending on the radial location within the perturbation. This "hysteresis-like" behavior arises from the non-linear feedback loop where density dictates the Compton scale, which in turn dictates the force and subsequent density evolution.

5. Significance

• Computational Efficiency: The perturbative scheme offers a faster alternative to full N-body simulations for studying non-linear structure formation in $f(R)$ gravity, potentially enabling rapid exploration of parameter spaces.
• Observational Signatures: The identification of the central density bump and the multi-valued DVDR provides distinct observational signatures. Future surveys (e.g., Euclid, DESI, Rubin Observatory) measuring the matter distribution and velocity fields could potentially distinguish Chameleon-screened $f(R)$ models from standard $\Lambda$CDM or unscreened modified gravity.
• Theoretical Insight: The work clarifies the spatio-temporal interplay between density and the Compton scale, showing that the screening mechanism is not static but evolves dynamically with the collapse, leading to complex internal structures in halos that are not captured by linear theory or simplified collapse models.

In summary, the paper successfully models the non-linear evolution of matter in $f(R)$ gravity with Chameleon screening using an efficient perturbative approach, revealing new structural features in collapsing halos and breaking the one-to-one correspondence between density and velocity divergence.