⚛️ general relativity

Density reconstruction from biased tracers: Testing the equivalence principle through consistency relations

This paper proposes and validates a practical quadratic estimator method for testing the weak equivalence principle using biased tracers in large-scale structure surveys, demonstrating that it achieves competitive constraints on equivalence principle violations compared to direct bispectrum measurements by leveraging mildly nonlinear scales.

Original authors: Lawrence Dam, Omar Darwish

Published 2026-07-17

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Lawrence Dam, Omar Darwish

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

Imagine the universe as a giant, invisible trampoline made of space and time. According to our best understanding of physics, gravity is just the curve in that trampoline caused by heavy objects like stars and galaxies. A fundamental rule of this cosmic game is the Weak Equivalence Principle. Think of it as a cosmic "fairness rule": if you drop a bowling ball and a feather in a vacuum, they fall at the exact same speed, regardless of what they are made of or how heavy they are. They follow the same path because gravity treats them all equally. We've tested this rule in labs on Earth and with satellites, and it works perfectly. But the universe is huge, and on the scale of billions of light-years, we haven't really checked if this fairness rule still holds. Could there be a hidden "fifth force" that pushes different types of matter in slightly different directions? If so, the universe's structure would look subtly different than we expect.

To find out, scientists look at the Large Scale Structure of the universe—the vast web of galaxies stretching across the cosmos. They use a mathematical tool called a consistency relation, which is like a cosmic recipe. It predicts exactly how the clustering of galaxies should behave if gravity is fair. If the recipe is broken, it means the Equivalence Principle is violated. The paper you are about to read explores a clever new way to check this recipe without getting bogged down in the messy, complicated math usually required.

The Cosmic Tug-of-War

The authors of this paper, Lawrence Dam and Omar Darwish, are tackling a tricky problem: how do we test if gravity treats all galaxies the same way when they are being pulled by a massive, invisible gravitational wave stretching across the universe?

Imagine two distinct types of galaxies, let's call them "Red Giants" and "Blue Dwarfs." If the Equivalence Principle holds true, and a giant gravitational wave (a "long-wavelength mode") passes through them, both types of galaxies should get pushed by the exact same amount. They would drift together, maintaining their distance, just like two passengers in an elevator accelerating upward; they don't feel the acceleration relative to each other.

However, if the Equivalence Principle is violated, the "Red Giants" might get pushed harder than the "Blue Dwarfs." This would cause them to drift apart or closer together in a very specific way. The authors call this a "dipole" or an "anti-symmetric shift." It's like a cosmic tug-of-war where one team is secretly stronger than the other, causing the rope to stretch unevenly.

The Problem with the Old Way

Usually, to spot this subtle drift, astronomers try to measure something called the bispectrum. In simple terms, the bispectrum measures how three different points in space are connected. It's a powerful tool, but it's also incredibly messy and computationally expensive. It's like trying to hear a whisper in a hurricane by recording every single sound wave and analyzing them one by one. It takes a lot of time and computing power, and it's easy to get lost in the noise.

The New Trick: The Quadratic Estimator

Instead of listening to the whole hurricane, the authors propose a smarter, simpler method using something they call a Quadratic Estimator (QE).

Think of the QE as a specialized filter. Instead of trying to measure the complex three-way connection directly, the QE looks at how the pattern of galaxy clustering changes when a large-scale gravitational wave is present. It's like noticing that the ripples in a pond change shape when a big boat passes by, without needing to track every single water molecule.

The authors realized that if the Equivalence Principle is violated, this "ripple pattern" (specifically the cross-correlation between two different types of galaxies) will develop a unique, tell-tale signature: a dipole that scales in a very specific way (proportional to $1/K$ , where $K$ is the size of the gravitational wave).

What They Did and Found

To test their idea, the team didn't just do math on paper; they ran simulations using a massive supercomputer and data that mimics the DESI (Dark Energy Spectroscopic Instrument) survey, which is currently mapping millions of galaxies.

Building the Tool: They built their Quadratic Estimator to specifically hunt for that "anti-symmetric shift." They fed it simulated data from the AbacusSummit simulations, which are like high-fidelity video games of the universe.
The Test: They checked if their estimator could successfully reconstruct the large-scale gravitational waves from the galaxy data. They found that it worked well, matching the theoretical predictions almost perfectly.
The Forecast: Since we don't have the final DESI data yet, they used a statistical method (Fisher forecast) to predict how well this tool would work with real data. They asked: "If we use this method on the DESI survey, how small of a violation can we detect?"

The Results:
Their new method is surprisingly effective. They found that their Quadratic Estimator is competitive with the much more complicated direct bispectrum method.

They forecast that using this method, a survey like DESI could constrain the "amplitude of EP violation" (a number they call $\tilde{\epsilon}$ ) to an uncertainty of roughly $6 \times 10^{-3}$ (or about 0.006).
This means they could detect a violation if it exists at a level of about 0.6% or larger.

Why This Matters

The paper doesn't claim to have found a violation of the Equivalence Principle. In fact, they assume the principle holds for their baseline simulations. Instead, they have built a better, faster, and more practical tool to look for it.

They argue that while the old way (direct bispectrum) is like using a sledgehammer to crack a nut, their Quadratic Estimator is like using a laser scalpel. It extracts the same crucial information about the "dipole" signature but is much easier to apply to real-world surveys like DESI, Euclid, and future projects.

The Catch

The authors are honest about the limitations. While their tool is great at finding the signal of a violation, there is a "degeneracy" (a confusion) with how galaxies form. We don't perfectly know how different types of galaxies "bias" (cluster) relative to the underlying dark matter. This uncertainty makes it hard to pin down the exact cause of any signal we find. It's like hearing a strange noise in your car; you know something is wrong, but without knowing the engine's exact specs, you can't be 100% sure if it's a broken belt or a loose bolt.

The Bottom Line

This paper is a "how-to" guide for the next generation of cosmic detectives. It shows that we don't need to wait for supercomputers to crunch billions of numbers to test the fundamental laws of gravity. By using a clever mathematical shortcut (the Quadratic Estimator), we can use current and upcoming galaxy surveys to test if gravity is truly fair to all matter in the universe. If the Equivalence Principle is broken, this tool will be the first to spot the cosmic drift.

Technical Summary: Density Reconstruction from Biased Tracers: Testing the Equivalence Principle Through Consistency Relations

Problem Statement
The weak equivalence principle (EP)—the assertion that test bodies follow identical trajectories in a gravitational field regardless of their internal structure—is a cornerstone of general relativity, rigorously tested on Solar System scales. However, its validity on cosmological scales remains less certain, with few direct tests. In the context of Large-Scale Structure (LSS), a violation of the EP implies that different galaxy tracers (e.g., distinct populations A and B) would experience different accelerations in response to a uniform, long-wavelength gravitational field.

Standard consistency relations in LSS dictate that the equal-time bispectrum in the squeezed limit (where one mode $K$ is much smaller than the others $k$ ) should be free of a $1/K$ pole. If the EP holds, the relative positions of galaxies are invariant under displacements induced by long-wavelength modes (analogous to Einstein's elevator). A violation of the EP would manifest as a specific dipole modulation in the squeezed bispectrum, scaling as $1/K$ , which cannot be transformed away by a change of coordinates. While direct bispectrum measurements can detect this, they face significant practical challenges regarding computational cost, covariance estimation, and systematic mitigation.

Methodology
The authors propose an alternative approach using quadratic estimators (QEs) to probe EP violations, avoiding the direct estimation of the bispectrum. The methodology proceeds as follows:

Linear Response Framework: The paper utilizes linear response theory to describe how long-wavelength matter fluctuations ( $\delta_L$ ) modulate the statistics of short-wavelength tracer modes ( $\delta_A, \delta_B$ ). The response function $f_{AB}(k_1, k_2)$ is decomposed into symmetric and anti-symmetric components across three physical effects: Growth (G), Shift (S), and Tidal (T).
Identification of the EP-Violating Signal: The authors demonstrate that the EP protects the anti-symmetric shift component ( $S^-$ ) from acquiring a non-zero coefficient in standard bias expansions. Specifically, the coefficient $C^S_{[AB]}$ vanishes if the EP holds. If the EP is violated (parameterized by $\epsilon$ ), this coefficient becomes non-zero, leading to a response function with a characteristic $1/K$ dipole dependence in the squeezed limit:
$f^{(-)}_S(k, K-k) \sim \frac{K \cdot k}{K^2} P_L(k)$
Quadratic Estimator Construction: Based on this specific response, the authors construct a "displacement estimator" ( $\hat{b}^D_{XY}$ ) to reconstruct the long-wavelength mode $\delta_L(K)$ from the product of two distinct tracer fields. The estimator weights the product of short modes by the theoretical response function.
$\hat{b}^D_{XY}(K) = \int_k w^D_{XY}(k, K-k) \delta_X(k) \delta_Y(K-k)$
The estimator is designed to be unbiased for the long mode but is sensitive to the anti-symmetric shift.
Forecasting: The authors apply this formalism to a survey configuration similar to DESI (Dark Energy Spectroscopic Instrument). They perform Fisher matrix forecasts using cross-correlations between the reconstructed long modes and the original galaxy tracers (and the auto-spectrum of the reconstruction) to constrain the EP violation amplitude.

Key Contributions

Theoretical Decomposition: The paper rigorously decomposes the gravitational response of biased tracers into a basis of symmetric and anti-symmetric modes (G, S, T). It identifies the anti-symmetric shift ( $S^-$ ) as the unique "smoking gun" for EP violation, distinguishing it from other effects that are either symmetric or suppressed by the EP.
Practical Estimator: The authors develop a sub-optimal but computationally efficient quadratic estimator tailored to isolate the $1/K$ dipole signature. This offers a practical alternative to the exhaustive direct bispectrum analysis, which is difficult to implement for large datasets.
Validation: The estimator is validated against AbacusSummit N-body simulations, showing that the reconstructed power spectra match theoretical predictions, including the noise scaling expected from Gaussian variance.

Results

Forecasted Constraints: Applying the estimator to a DESI-like survey (volume $\sim 58 h^{-3} \text{Gpc}^3$ ), the authors forecast constraints on the combined EP violation parameter $\tilde{\epsilon} = \epsilon \times C^S_{[AB]}$ .
Competitiveness: By utilizing modes up to mildly nonlinear scales ( $k_{\text{max,rec}} \approx 0.15 h \text{Mpc}^{-1}$ ), the quadratic estimator achieves a constraint of $\sigma(\tilde{\epsilon}) \sim 6 \times 10^{-3}$ . This is competitive with forecasts from direct full-bispectrum analyses (e.g., Graham et al. 2025) and comparable to constraints derived from CMB+DESI BAO data regarding modified growth rates.
Bias Degeneracy: The authors highlight a significant limitation: the inability to disentangle the fundamental EP violation parameter $\epsilon$ from the unknown astrophysical bias parameters ( $b_{\epsilon, X}$ ) without additional modeling. Consequently, they constrain the degenerate product $\epsilon \times C^S_{[AB]}$ . They note that without knowledge of the relative bias parameters, direct constraints on $\epsilon$ alone are weaker than those from CMB-based growth rate analyses.

Significance and Claims
The paper claims that quadratic estimators provide a robust and practical tool for testing the Equivalence Principle on cosmological scales.

Robustness: Unlike growth-based tests, the $1/K$ consistency relation signature is robust against nonlinear gravitational evolution, baryonic physics, and redshift-space distortions, as it relies on the symmetry properties of the squeezed limit rather than detailed modeling of structure formation.
Feasibility: The authors argue that current and upcoming surveys (like DESI) can immediately benefit from incorporating this QE approach, as it extracts information from the squeezed bispectrum without the computational overhead of direct bispectrum estimation.
Distinctiveness: The method targets a signature (the anti-symmetric shift) that is distinct from primordial non-Gaussianity (which typically scales as $K$ or $1/K^2$ but lacks the specific dipolar structure associated with EP violation in the shift sector), although the authors acknowledge that a combined treatment of PNG and EP violation is necessary for a complete analysis.

The work concludes that while the degeneracy with unknown bias parameters limits the precision of $\epsilon$ alone, the detection of a non-zero anti-symmetric shift amplitude would constitute direct evidence of EP violation, offering a unique probe of new physics on large scales.