Analytic and numerical study of dark matter halos in the Horndeski theory

Technical Summary: Analytic and Numerical Study of Dark Matter Halos in the Horndeski Theory

Problem Statement
Modified gravity (MG) theories, such as the Horndeski theory, offer alternatives to the cosmological constant for explaining the late-time accelerated expansion of the Universe. However, to remain viable, these models must satisfy stringent local gravity tests (e.g., in the solar system) by suppressing "fifth forces" at small scales via screening mechanisms, specifically the Vainshtein mechanism. While the Vainshtein radius for stellar systems is small ($\sim 100$ pc), dark matter halos exist on intermediate scales where the transition between screened and unscreened regimes may occur. A critical open question is whether the efficiency of Vainshtein screening depends on the detailed density profile of the halo, particularly in the presence of "core-like" structures, and whether residual fifth-force signatures can survive the dynamical evolution of these halos to be detectable via observational mass proxies.

Methodology
The authors employ a dual approach combining analytical derivation and numerical simulation:

1. Analytical Framework:

   • The study utilizes the Horndeski theory with a spherically symmetric halo model.
   • The authors derive asymptotic solutions for the scalar field and gravitational potentials, defining the Vainshtein radius ($r_V$) where the screening condition $A(r_V)=1$ is met.
   • They analyze the scaling behavior of the fractional acceleration deviation ($\Delta_{acc}$) relative to the density profile slope ($\gamma$) and the core concentration parameter ($c_s$).
   • A key analytical focus is the behavior of the dimensionless source term $A(r)$ in the presence of finite-density cores, investigating whether $A(r)$ saturates to a value that prevents effective screening even near the halo center.

2. Numerical Simulations (N-body):

   • The authors perform N-body simulations using the Gadget-2 code for isolated, spherically symmetric halos.
   • Setup: The simulation assumes Newtonian two-body interactions for particle dynamics, while the fifth force is treated as an external central force derived from the enclosed mass. The scalar field equation is solved at each time step to determine the fifth-force acceleration.
   • Models: Simulations are run for both Navarro-Frenk-White (NFW) profiles (cuspy) and Burkert profiles (cored) with varying core concentration parameters ($c_s$).
   • Parameters: The study focuses on the regime where $\beta \approx 0$ (consistent with GW170817 constraints on gravitational wave speed) but explores the impact of the halo profile on screening efficiency.

3. Mass Proxies:
   To detect deviations from General Relativity (GR), the authors compute three distinct mass proxies from the simulation outputs and compare them against the "true" enclosed mass ($M_{true}$):

   • Hydrostatic Equilibrium Mass ($M_{HSE}$): Derived from the gas pressure gradient (assuming a polytropic relation).
   • Jeans Mass ($M_{Jeans}$): Derived from the radial velocity dispersion of dark matter particles.
   • Lensing Mass ($M_{WL}$): Derived from the weak gravitational lensing convergence (Weyl potential).

Key Contributions and Results

• Analytical Insight on Core Profiles: The authors demonstrate that the efficiency of Vainshtein screening is not solely determined by the total mass but is sensitive to the inner density profile. For cored profiles (where $\gamma \to 0$), the dimensionless source $A(r)$ saturates to a finite value ($A_{core}$) near the center. If $A_{core}$ is not sufficiently large ($A_{core} \sim \mathcal{O}(1)$), the screening mechanism fails to fully suppress the fifth force even within the nominal Vainshtein radius. This leads to a "residual" fifth force in the central regions of cored halos.
• Numerical Verification of Residual Effects: The N-body simulations confirm that halos with extended cores can maintain a dynamical (quasi-)equilibrium state where screening remains inefficient in the central regions. Consequently, the fifth force persists toward the halo center, contrary to the expectation that screening is always effective deep inside the Vainshtein radius.
• Bias in Mass Proxies:
  • Dynamical Proxies ($M_{HSE}$ and $M_{Jeans}$): In the unscreened outer regions, these proxies systematically overestimate the true mass by approximately 30% (consistent with $\Delta_{acc} \approx 1/3$). Crucially, in cored halos, this overestimation persists even into the inner regions where GR would predict effective screening.
  • Lensing Proxy ($M_{WL}$): The lensing mass remains in close agreement with the true mass ($M_{true}$) across all radii. This is because the Weyl potential ($\Phi_{Weyl} = (\Phi + \Psi)/2$) probed by lensing cancels out the leading-order scalar field contributions (proportional to $\xi x$) that drive the fifth force in the dynamical equations.
  • Conclusion on Bias: The study finds a systematic bias between dynamical mass proxies and lensing mass proxies in the presence of modified gravity. Specifically, $M_{HSE}$ and $M_{Jeans}$ deviate from $M_{WL}$ and $M_{true}$, providing a potential observational signature of the unscreened fifth force.

Significance and Claims
The paper claims to provide a comprehensive study linking the structural properties of dark matter halos (specifically core concentration) to the efficiency of Vainshtein screening in the Horndeski theory.

• Theoretical Implication: The authors assert that the assumption of a fixed density profile (e.g., NFW) is insufficient for testing MG; the dynamical evolution of the halo profile itself dictates the screening efficiency. They posit that "core-like" structures allow fifth-force signatures to survive the screening mechanism.
• Observational Utility: The results suggest that comparing multiple mass proxies (dynamical vs. lensing) in galaxy clusters can serve as a test for gravity on halo scales. A discrepancy where dynamical masses exceed lensing masses could indicate a violation of Newtonian gravity/GR due to a residual fifth force.
• Modesty and Limitations: The authors explicitly state that their implementation treats the fifth force as a central force and approximates two-body interactions as Newtonian, which is an approximation of the full Horndeski dynamics. They acknowledge that relaxing spherical symmetry (anisotropy) could further affect screening but argue their spherical analysis serves as a leading-order prescription. They do not claim to have calibrated these biases for real observational data, noting that error estimation and calibration are left for future work.

In summary, the paper argues that the interplay between halo density profiles and screening mechanisms creates a regime where fifth forces may be observable through systematic biases in mass measurements, offering a new avenue for testing the Horndeski theory with upcoming multi-wavelength observations of galaxy clusters.