Scale-Dependent Growth of Large-Scale Structure
from Infrared Modifications of Gravity

Vitantonio Castronuovo

Independent Researcher

ORCID iD 0009-0003-0066-0738

Published in Annalen der Physik

Graphical abstract of the protected infrared kernel
Graphical Abstract. Visual representation of the theoretical framework, showing the link between microscopic ultraviolet completions and macroscopic cosmological observables via the infrared single-pole kernel.

Start Here: The Idea in Plain Language

This work asks a simple question: if gravity is almost exactly Einstein’s theory, could it still behave slightly differently on the largest cosmic scales?

The framework does not change how light bends, and it does not change the speed of gravitational waves. Those sectors are protected because observations already constrain them very strongly. Instead, the theory focuses on the growth of cosmic structure: how galaxies, clusters, and large-scale matter patterns assemble over time.

The central idea is that late-time deviations satisfying the protected-branch assumptions are highly restricted. In the leading effective description used in the paper, the possible growth modification is compressed into two quantities: an amplitude, which tells how large the effect is, and a physical scale, which tells where the effect turns on.

Study Summary

The study proposes a theoretical model for controlled discrepancies in the growth of cosmic structure through infrared modifications of gravity. It introduces a two-parameter kernel, defined by a coupling amplitude and a mass scale, that can modify galaxy formation and large-scale clustering while preserving the observational constraints on gravitational-wave speed and GR-like lensing.

Within the stated assumptions, representative physical theories, from extra-dimensional scenarios to screening mechanisms and hidden-sector mediators, can reduce to the same effective mathematical description over the observable linear window. Benchmark forecasts for missions such as DESI and Euclid outline a strategy for testing these deviations from General Relativity.

The broader program is falsifiable rather than confirmatory: cosmological data, spectroscopic probes, void statistics, varying-constant searches, collider limits, and fifth-force constraints can progressively distinguish or exclude the relevant ultraviolet classes. The framework does not identify the origin of dark matter; it provides a controlled route for testing whether late-time structure growth carries evidence of infrared gravitational physics.

The Three Key Formula Ideas

The paper uses equations, but the logic can be read without following every mathematical step. The main symbols have a direct physical meaning.

1. \(\mu(k,a)\): how gravity affects matter clustering If \(\mu = 1\), structure grows as in General Relativity. If \(\mu\) differs from 1, matter clustering changes in a scale-dependent way.
2. \(\beta_0\): the strength of the effect \(\beta_0 = 0\) means no deviation from Einstein gravity in the growth sector. A nonzero value would mean that cosmic structure feels an additional effective response.
3. \(m_\ast\): the scale of the effect The mass scale \(m_\ast\) sets the characteristic cosmic distance where the effective modification becomes relevant. The framework therefore tests a specific turnover scale, not an arbitrary deviation.

Two further conditions keep the proposal tightly constrained: \(\Sigma = 1\) means that light bending remains GR-like, and \(c_T = 1\) means that gravitational waves still travel at the speed of light.

The protected single-pole kernel
Conceptual illustration. The protected single-pole kernel summarizes the two-parameter growth response under the conditions \(\Sigma=1\) and \(c_T=1\). It is used here as a visual guide to the theoretical structure, not as an additional numerical result.

What the Theory Is Not Claiming

This framework does not claim that General Relativity has already been observed to fail. It does not introduce an arbitrary new force, and it does not identify a specific dark-matter particle.

Its value is its falsifiability. Future surveys can either detect the predicted scale-dependent pattern or rule out this protected single-pole branch within the assumptions and analysis window stated in the paper.

What Publication in Annalen der Physik Means

The paper has been published in Annalen der Physik. This means that the theoretical framework and its stated assumptions have passed peer review. It does not mean that the effect has been observed; rather, it means that the proposed testable structure is considered suitable for publication as a theoretical physics contribution.

Cosmological Context: Why Test Gravity?

The ΛCDM model remains the reference framework of modern cosmology. It successfully describes the expansion history of the Universe, the cosmic microwave background, and the formation of large-scale structure. However, several tensions between early-universe inferences and late-universe observations continue to motivate sharper tests of gravity.

Two examples are the Hubble tension, related to the present expansion rate, and the S8 tension, related to the degree of matter clustering. These tensions do not by themselves prove that General Relativity is wrong, but they identify the regimes where new theories must be tested most carefully.

The work presented here belongs to that program: it does not propose an arbitrary modification of gravity, but a restricted and falsifiable class of deviations in the growth sector. Light bending and gravitational-wave propagation remain protected by existing observations; the possible deviation is searched for in the growth of matter on large scales.

Surveys such as DESI, Euclid, Rubin/LSST, and future CMB experiments will help determine whether these tensions reflect observational systematics, dark-sector physics, or a controlled modification of gravity on cosmological scales.

Related Direction: f(Q) Gravity A related line of research studies symmetric teleparallel gravity, or f(Q) gravity, where gravity is described through non-metricity rather than curvature. This direction is distinct from the main paper presented here, but it belongs to the same broader program: building modified-gravity models that are mathematically controlled, observationally falsifiable, and compatible with lensing and gravitational-wave constraints.

1. Introduction: Cosmic Tensions and the Standard Cosmological Model

The accelerated expansion of the late-time Universe, discovered through supernova observations nearly three decades ago, remains one of the central open problems in fundamental physics. Although General Relativity (GR) paired with a cosmological constant (\(\Lambda\)) provides an excellent fit to current data, it raises deep theoretical questions. The extreme smallness of the vacuum energy, combined with the presence of mild but persistent tensions in low-redshift observables—such as discrepancies between early-universe and late-universe measurements of the Hubble constant (\(H_0\)) and the clustering amplitude (\(S_8\))—suggests that our understanding of cosmic evolution is incomplete.

These persistent tensions have motivated a wide class of modified-gravity and dark-energy scenarios. Until recently, many such models remained viable, as scalar-tensor extensions, phenomenological parameterizations, and extra-dimensional constructions could reproduce late-time observations while evading local solar-system constraints. However, increasingly precise cosmological data have made the allowed theoretical window much narrower.

2. The Rules of the Game: The Protected Branch of Gravity

Modifying General Relativity is a highly constrained theoretical problem. Any viable theory must reproduce the stringent tests that Einstein’s theory has already passed. The watershed moment came with the detection of gravitational waves from neutron-star mergers and their electromagnetic counterparts (GW170817), which established that gravitational waves propagate at the speed of light to very high precision. This ruled out large classes of theories with modified tensor propagation, enforcing the condition \(c_T = 1\) within the observational accuracy.

Together with weak-lensing constraints on the relation between the two gravitational potentials, these results impose a strong restriction on viable late-time modifications of gravity. The protected branch considered in the paper imposes luminal tensor propagation and an operationally GR-like lensing sector (\(\Sigma = 1\)). Within that branch and within the quasi-static linear regime, departures from General Relativity are confined to the Newtonian potential governing the large-scale clustering of matter.

3. Theoretical Framework: Trace-Coupled Mediators

Within this protected branch, the paper studies a controlled class of trace-coupled scalar mediators. Late-time gravity is represented by an effective scalar degree of freedom (\(\phi\)) coupled universally to matter through the trace of the energy-momentum tensor (\(T \equiv T^\mu_\mu\)). A representative low-energy action has the canonical structure:

$$ S = \int d^4x \sqrt{-g} \left[ \frac{1}{2} Z(\phi)(\partial\phi)^2 - V(\phi) + \frac{\phi}{M_{\rm eff}} T \right] $$
What does this mean in simple terms?
This action is the effective starting point of the model class. It states that matter, through the trace \(T\), can source a scalar mediator \(\phi\). The scale \(M_{\rm eff}\) controls the coupling strength in this low-energy description.

By expanding this theory around a homogeneous cosmological background and working on sub-horizon scales in the quasi-static regime, the scalar perturbation satisfies the Klein-Gordon equation:

$$ \left(\nabla^2 + a^2 m_{\rm eff}^2\right)\varphi = \frac{a^2}{M_{\rm eff}} \delta T $$
What does this mean in simple terms?
This describes the linear response of the mediator in the quasi-static regime. The symbol \(\nabla^2\) measures spatial variation, while \(m_{\rm eff}\) sets a finite range. The result is a Yukawa-type response rather than an unrestricted long-range force.

4. The Modified Poisson Equation & The Single-Pole Kernel

On the scales probed by modern large-scale structure surveys, time derivatives are subdominant within the stated quasi-static regime. Solving the field equations then yields a modified Poisson equation for the Newtonian gravitational potential (\(\Psi\)):

$$ -k^2\Psi = 4\pi G a^2 \rho_m \delta_m \, \mu(k,a) $$
What does this mean in simple terms?
The original Poisson equation states that matter sources the Newtonian gravitational potential (\(\Psi\)). Here, \(k\) is the Fourier wave number, which labels the physical scale of a perturbation. The factor \(\mu\) describes a scale-dependent effective response, so the clustering strength can differ from its GR value in a controlled way.

Where \(\mu(k,a)\) represents the Single-Pole Infrared Kernel, defined as:

$$ \mu(k,a) = 1 + \beta_0(a) \frac{k^2}{k^2 + a^2 m_\ast^2} $$
What does this mean in simple terms?
This is the central object of the framework. If \(\beta_0\) is zero, the growth response reduces to its GR value. For positive \(\beta_0\), the effective clustering response is enhanced above the turnover scale. The mass scale \(m_\ast\) sets the physical distance at which the additional response becomes relevant. Within the declared assumptions, the modification is compressed into a controlled two-parameter form rather than an arbitrary function.
Conceptual explanation of the infrared growth kernel
Conceptual illustration. Original website diagram of the single-pole infrared response: growth remains GR-like in the low-response regime and may become scale dependent across a turnover while the protected branch maintains \(\Sigma=1\) and \(c_T=1\). This is an explanatory schematic, not a reproduced article figure or a numerical result.

5. Gravitational Slip and Observational Degeneracies

Maintaining \(\Sigma = 1\) (GR-like lensing) while simultaneously allowing for a modified growth rate (\(\mu \neq 1\)) requires a deeply correlated adjustment of the two scalar gravitational potentials, \(\Phi\) and \(\Psi\). Operationally, the Weyl potential sourcing light deflection remains identical to its General Relativity value: \( -k^2(\Phi + \Psi) = 8\pi G a^2 \rho_m \delta_m \).

This implies an effective gravitational slip between the space-curvature potential (\(\Phi\)) and the time-curvature potential (\(\Psi\)). We can parameterize this slip as \(\eta \equiv \Phi / \Psi\), which on the protected branch and at leading resolved order evaluates to:

$$ \eta = \frac{2-\mu}{\mu} \simeq 1 - \beta_0 \frac{k^2}{k^2 + a^2 m_\ast^2} + \mathcal{O}(\beta_0^2) $$
What does this mean in simple terms?
In GR without anisotropic stress, the two scalar potentials are equal. In the protected branch, the potentials may respond differently while their sum, which controls lensing, remains fixed. The slip parameter \(\eta\) is the bookkeeping relation that keeps growth modifications compatible with a GR-like Weyl potential.

6. Ultraviolet Completions: Where Does the Kernel Come From?

A central result of this theoretical program is that several distinct ultraviolet (UV) completions can reduce, at leading resolved order, to the same two-parameter infrared kernel on cosmological scales. We classify these UV origins into three representative benchmark classes:

  • Class I: Randall-Sundrum and Discrete Spectra. Warped extra-dimensional models decompose gravity into a massless graviton and a tower of massive Kaluza-Klein (KK) excitations. Stabilization can introduce a scalar radion. When the response is dominated by the lightest mass eigenstate, the effective kernel is discrete (\(m_\ast = m_1\)) and may be accompanied by configuration-space oscillatory features.
  • Class II: Chameleon and Scalar-Tensor Theories. Environmental screening mechanisms can render the scalar mass density dependent, generating a continuum spectral density. If the scalar also couples to Standard Model sectors, variations of fundamental constants may provide secondary tests.
  • Class III: Hidden-Sector Mediators. A light scalar may reside in a hidden sector and mix weakly with the gravitational trace. Such realizations can reproduce the same leading growth kernel while avoiding the secondary signatures associated with visible-sector couplings.
Property Class I (Extra-Dimensions) Class II (Chameleons) Class III (Hidden Sector)
Spectrum Discrete Continuum Continuum
Typical Amplitude \(10^{-8}\) – \(10^{-4}\) \(10^{-3}\) – \(10^{-1}\) \(10^{-2}\) – \(10^{-1}\)
Varying Constants No Yes No
Void Oscillations Yes No No

7. Looking to the Future: The Era of Precision Cosmology

The strength of this theoretical framework lies in its operational falsifiability. It provides a structured blueprint for future observations. The predicted scale-dependent growth gives a specific signature that can be searched for by the next generation of cosmological observatories.

Fisher-matrix benchmarks indicate that upcoming Stage-IV missions—such as the Dark Energy Spectroscopic Instrument (DESI) and the Euclid space telescope—may have enough precision to test this scale-dependent signal, provided that survey systematics and parameter degeneracies are controlled.

What does this mean in simple terms?
A Fisher matrix is a forecasting tool used to estimate parameter sensitivity before data are analyzed. It does not guarantee discovery; it indicates whether a survey such as Euclid could distinguish the predicted signal from noise and degeneracies under specified assumptions.
Conceptual path from infrared kernel to observational testing
Conceptual illustration. Original website diagram of the testing logic: an infrared response template induces a scale-dependent growth signature that can be confronted with large-scale-structure observables. It does not display measured data or reproduce a forecast figure from the article.

8. Supporting Information: Spectral Justification and UV Origins

To make the mathematical assumptions explicit and to show why the effective two-parameter formulation is not merely an ad-hoc phenomenological fit, we provide the formal spectral justification of the infrared kernel. The linear response of any weakly coupled sector entering the Newtonian potential is governed by a retarded propagator \(G_{\rm ret}\) that admits a Källén-Lehmann spectral representation:

$$ G_{\rm ret}(k,a) = \int_0^\infty dm^2\, \frac{\rho(m^2,a)}{k^2 + a^2 m^2} $$

where, in a controlled unitary spectral setting, the spectral density satisfies \(\rho(m^2,a) \ge 0\). The corresponding modification of the Poisson equation can then be written as an infinite-dimensional superposition of Yukawa terms:

$$ \mu(k,a) - 1 = \int_0^\infty dm^2\, \rho(m^2,a)\, \frac{k^2}{k^2+a^2 m^2} $$
What does this mean in simple terms?
This is a standard idea in quantum field theory: a linear response can often be decomposed into a spectrum of elementary contributions. In controlled unitary settings the spectral weight is non-negative, which supports a stable attractive response. Here the formula explains why a dominant infrared scale can reduce to the two-parameter kernel.

8.1 Quasi-Adiabatic Evolution

In the protected branch, the scalar mediator is assumed to be sufficiently heavy compared to the Hubble expansion rate (\(m_{\rm eff} \gg H(a)\)). Its evolution is quasi-adiabatic and can be treated within a WKB-type approximation. The linearized mode equation takes the form:

$$ \ddot{\varphi}+3H\dot{\varphi} + \left(\frac{k^2}{a^2}+m_{\rm eff}^2(a)\right)\varphi = \text{source} $$

8.2 Randall-Sundrum: From 5D Action to the 4D Poisson Kernel

When explicitly deriving this from a higher-dimensional framework, we consider the original Randall-Sundrum geometry with a single warped extra dimension. The five-dimensional gravitational action and metric are described by:

$$ ds^2 = e^{-2ky}\left(\eta_{\mu\nu}+h_{\mu\nu}\right)dx^\mu dx^\nu - dy^2 $$
What does this mean in simple terms?
This equation describes a universe with a hidden, microscopic fifth dimension (\(y\)). The term \(e^{-2ky}\) acts like a warp factor. Stabilizing the extra dimension can introduce a radion-like scalar mode in the four-dimensional effective theory, producing a response of the same kernel form.

Linearized metric perturbations admit a Kaluza-Klein (KK) decomposition. By projecting the trace of the energy-momentum tensor and ensuring canonical normalization, we obtain the precise amplitude for the RS model:

$$ \beta_{0}^{\rm RS} = \frac{\tilde{\lambda}_1^{\,2}\, f_{\rm KK}}{16} \sim 3 \times 10^{-11} $$

8.3 Parametric Screening and Chameleon Theories

Alternatively, for continuum realizations like Chameleon screening, the local effective mass of the field is dynamically generated by an interacting potential. A representative chameleon-like effective potential takes the form:

$$ V_{\rm eff}(\phi;\rho) = \tfrac{1}{2} m_\ast^2 \phi^2 + \alpha \phi \rho + \tfrac{\kappa}{3}\phi^3 $$
What does this mean in simple terms?
This equation describes a chameleon-like scalar. Because it couples to the matter density (\(\rho\)), its effective mass can become large in dense environments, suppressing local deviations from GR. In low-density cosmic regions it can remain lighter and contribute to a long-range effective response.

In high-density regions (\(\rho \gg \rho_{\rm cosmo}\)), the effective mass is enhanced, \(m_{\rm eff}^2(\rho) \simeq m_\ast^2 + 2\kappa \bar{\phi}(\rho)\). Requiring Solar-System screening at the order-of-magnitude level gives \(\kappa \alpha \gtrsim 10^{-60} \, {\rm eV\;cm^3/g}\). For representative chameleon-like potentials this bound can be satisfied without extreme parameter choices, while a full PPN analysis remains a separate requirement.

9. Implications for the Scientific Community

If future large-scale-structure data were to favor the infrared single-pole kernel over \(\Lambda\)CDM, it would provide evidence for a protected-branch modification of late-time growth on cosmological scales. Such a result would require independent cross-checks of systematics, lensing consistency, and parameter degeneracies before any claim of new gravitational physics.

Detecting this scale-dependent signature would open a macroscopic window onto possible dark-sector or infrared gravitational physics. The microscopic origin would remain class-conditional and would have to be distinguished through secondary observables such as void correlations, varying-constant searches, fifth-force constraints, or collider limits.

About the Author & Research

Vitantonio Castronuovo

I am an independent researcher born in Turin, Italy, and currently based in Madrid, Spain. My work focuses on theoretical cosmology, modified gravity, and the possibility that late-time large-scale structure may carry controlled signatures of infrared physics beyond the standard cosmological model.

My paper Scale-Dependent Growth of Large-Scale Structure from Infrared Modifications of Gravity has been published in Annalen der Physik* This website presents the core idea in accessible language while preserving the distinction between a peer-reviewed theoretical framework, its observational tests, and any future empirical confirmation.

*Established in 1790, Annalen der Physik is one of the oldest and most prestigious physics journals in the world. It is historically renowned for publishing landmark scientific papers, including Albert Einstein's 1905 Annus Mirabilis papers (special relativity, photoelectric effect, Brownian motion, and mass-energy equivalence $E=mc^2$), as well as foundational works by Max Planck, Wilhelm Röntgen, and Heinrich Hertz.

Historic Contributions in Annalen der Physik

  • 1. Albert Einstein & the Annus Mirabilis (1905)
    Einstein published four historic papers that revolutionized classical physics: on the photoelectric effect, Brownian motion, special relativity, and mass-energy equivalence.
  • 2. Einstein's General Relativity (1916)
    The foundation of general relativity, describing gravity as spacetime curvature.
  • 3. Max Planck & the Birth of Quantum Physics (1900)
    The quantum theory of radiation, introducing Planck's constant.
  • 4. Wilhelm Röntgen & X-Rays (1896)
    The discovery of X-rays.
  • 5. Heinrich Hertz & Electromagnetic Waves (1887)
    The experimental proof of electromagnetic waves.
  • 6. Erwin Schrödinger & the Wave Equation (1926)
    The wave mechanics formulation.

Other notable historical contributions:

  • Gustav Kirchhoff (1859): Formulated the law of thermal radiation, stating that a good absorber of heat is also a good emitter.
  • Friedrich Kohlrausch (1874): Established the laws of electrical conduction in electrolytes.

Multimedia & Resources

Explore the Theory

Explore the theory through ten core resources: an interactive simulator, a glossary, a physical analogy, an observational timeline, a derivation guide, hi_class integration notes, the validation roadmap, reproducibility tools, and direct article links.

🔗

Official Record via DOI

DOI assigned: 10.1002/andp.70230. Opens the Wiley record when available online.

Validation Protocol

Guidelines for integrating and testing the infrared kernel in Large-Scale Structure (LSS) analysis pipelines.

1. Theoretical Foundations and Protected Branch Constraints

Implementation in large-scale-structure analysis pipelines must operate inside the protected branch. This preserves consistency with GW170817 and weak-lensing constraints, while isolating departures from General Relativity in the linear growth sector.

1.1 Invariance Assumptions and Coupling

The framework assumes Lorentz invariance and universal metric coupling. The scalar mediator couples to matter via the trace of the energy-momentum tensor (\(T = T^\mu_\mu\)). This guarantees a local, unitary EFT free from ghost instabilities, restricting modifications to a leading-order IR single-pole response.

1.2 Observational Constraints and Gravitational Slip

The protected branch imposes \(c_T=1\) and \(\Sigma=1\). Tensor propagation remains luminal, and the Weyl potential remains GR-like. Modified growth is then carried by \(\mu(k,a)\), with the scalar potentials adjusted so that the lensing sector is not independently modified.

Operationally one may write \(\Psi=\mu\Psi_{\rm GR}\), \(\Phi=(2-\mu)\Psi_{\rm GR}\), and \(\eta\equiv\Phi/\Psi=(2-\mu)/\mu\). This keeps the growth modification tied to the protected lensing condition.

1.3 Physical Origin of the Mediator

The same infrared response can arise as the leading observable projection of distinct ultraviolet classes, including warped extra-dimensional modes, screened scalar sectors, and weakly coupled hidden mediators.

2. Single-Pole Kernel Architecture

The kernel is a leading infrared response designed to minimize parametric freedom while retaining a controlled physical scale.

2.1 Mathematical Form and Distinction from Minimal \(f(Q)\)

The modified Poisson response is \(\mu(k,a)=1+\beta_0(a)k^2/(k^2+a^2m_\ast^2)\). Here \(\beta_0\) controls the coupling amplitude and \(m_\ast\) sets the effective transition scale.

Unlike many minimal metric-affine \(f(Q)\) implementations, which are often nearly scale independent in the quasi-static regime, this kernel predicts a single turnover together with infrared recovery and ultraviolet saturation.

2.2 Quasi-Static and Adiabatic Evolution

The response is Yukawa-like: the modification is suppressed for \(k\ll a m_\ast\). The template is stable when \(m_{\rm eff}\gg H\), so the mediator evolves adiabatically and \(m_\ast\) can be treated as constant over the declared survey window, with residual evolution absorbed into \(\beta_0(a)\).

3. Signal Isolation and Degeneracy Breaking

The statistical task is to separate \(\beta_0\) and \(m_\ast\) from degeneracies with \(\sigma_8\), \(\Omega_m\), nuisance bias terms, and survey systematics.

3.1 Degeneracy Logic and Stage-IV Benchmarks

The primary discriminator is scale dependence. A change in \(\sigma_8\) mostly rescales the spectrum, while the kernel produces a localized turnover \(k_{\rm turn}(z)=a(z)m_\ast\). Representative Fisher-level benchmarks are \(\sigma(\beta_0)\simeq0.008\), \(\sigma(m_\ast)/m_\ast\simeq0.19\), and \(\rho(\beta_0,m_\ast)\simeq0.68\).

3.2 Multi-Probe Isolation

Redshift-space distortions map the scale dependence of \(f\sigma_8(k,z)\). CMB priors from Planck or Simons Observatory restrict background parameters, limiting the freedom in \(\Omega_m\) and \(h\).

4. Numerical Implementation: System Requirements

Pipeline implementation should be performed inside modified Einstein-Boltzmann solvers or validated forward models, with the protected-branch assumptions imposed explicitly.

4.1 hi_class Configuration

In an EFT/Horndeski implementation, the protected submanifold may be represented through the Bellini-Sawicki relation \(\alpha_B=-\alpha_M/2\), which enforces the lensing condition while allowing the growth response to reproduce the single-pole structure.

4.2 Fast Growth Template

For preliminary sensitivity checks, one may use a schematic template \(P_{\rm matter}(k,a)\simeq P_{\Lambda{\rm CDM}}(k,a)[\mu(k,a)]^2\). This is only a fast diagnostic; final validation must use the full linearized system and survey covariance.

5. Validation Roadmap and Falsification Hierarchy

The validation process follows a sequential logic. It first tests whether a nonzero growth response is present, then asks whether the response is compatible with a single physical scale, and finally uses secondary signatures to separate ultraviolet classes.

5.1 Likelihood Test for \(\beta_0\)

If the improvement over \(\Lambda\)CDM is below the predeclared threshold, the nonzero branch is not selected. A resolved \(\beta_0\neq0\) in growth, without an independent lensing anomaly, is the entry condition for the protected branch.

5.2 Spectral Test and Oscillatory Signature

Discrete spectra, such as Randall-Sundrum-like towers, can generate oscillatory residuals in configuration-space observables, schematically \(\Delta\xi(r)\propto\sin(m_\ast r)/(m_\ast r)\). A monotonic response instead points toward screened continua or hidden-sector mediators.

5.3 Matter Couplings and Fundamental Constants

Quasar spectroscopy and varying-constant searches can test visible-sector couplings. A correlation between \(\Delta\alpha/\alpha\) and local density would favor chameleon-like models, while null variation would favor purely gravitational or hidden-sector mediators.

6. Result Interpretation and UV-IR Synthesis

The final interpretation is class-conditional. The same infrared kernel can be shared by distinct ultraviolet mechanisms, so secondary observables are needed before assigning a microscopic origin.

UV Class Spectral Signature Fundamental Constants Collider Relevance
Randall-Sundrum Oscillations \(\sin(m_\ast r)/(m_\ast r)\) Constant Spin-2 or radion resonances
Chameleon Monotonic Variable \(\Delta\alpha/\alpha\) Typically silent because of screening
Hidden Sector Monotonic Constant Missing energy or contact interactions

6.2 Success Criteria and Quantum-Gravity Context

A successful protocol must discriminate these classes within the controlled linear or quasi-linear regime. In fundamental terms, \(m_\ast\) may be read as the first resolved infrared scale of a tower or mediator sector, while the framework does not require it to source the background acceleration.

The single-pole kernel is therefore the minimal stable theoretical baseline for testing GR on cosmological scales through large-scale-structure data.

FAQ

Direct answers to common theoretical questions, scope limits, and implementation caveats.

Why exactly \(\Sigma=1\)?

It is not a generic prediction of all scalar-tensor theories. It is the defining condition of the protected branch studied in the paper: light deflection remains GR-like while the possible deviation is confined to the growth sector. If data require a resolved lensing deviation, that result would move the signal outside this protected branch.

Why is \(c_T=1\) imposed?

The near-luminal propagation of gravitational waves is a hard observational constraint after GW170817 and its electromagnetic counterpart. The framework therefore does not use modified tensor speed as a degree of freedom; it asks what remains testable in late-time structure growth after this sector is protected.

Has the effect been observed?

No. The published result is a theoretical and phenomenological framework, together with benchmark-level sensitivity estimates and validation targets. Publication means the construction and assumptions passed peer review; it is not an observational detection of modified gravity.

Is the kernel an arbitrary fit?

No. Under the stated assumptions, the single-pole form is the minimal two-parameter effective response for one dominant resolved infrared scale. This is a controlled minimality statement, not a theorem that all possible ultraviolet theories must reduce to this kernel.

What does \(m_\ast\) represent?

It is the effective transition scale of the resolved growth response. Depending on the ultraviolet realization, it may correspond to a mediator mass, a screened scalar scale, or the first resolved infrared scale of a tower. The paper does not claim a unique microscopic origin for it.

Does this identify dark matter or dark energy?

No. The framework keeps the background close to \(\Lambda\)CDM and tests whether late-time growth can carry a protected infrared modification. It does not identify a dark-matter particle and does not require the same sector to source cosmic acceleration.

Is this already implemented in hi_class?

Not as a public one-line switch. The paper and supplement define a validation target: a consistent EFT/Horndeski or forward-model implementation should reproduce the target \(\mu(k,a)\) while preserving \(c_T=1\), \(\Sigma=1\), stability, and the declared analysis window.

How is this different from standard \(f(R)\) or \(f(Q)\) models?

The comparison is class-conditional. Many minimal \(f(Q)\) implementations are nearly scale independent in the quasi-static regime, while the protected single-pole branch predicts one resolved turnover. Standard \(f(R)\)-like scalar sectors can overlap only after their screening, stability, and lensing behavior are controlled.

Interactive Visual Sandbox

Use the two sliders to modify the amplitude and transition scale of the kernel. The chart compares the effective response with the GR reference line, making the spatial threshold visible.

UV Completion Class
0.40
0.08
0.00
GR baseline
Single-pole response

Small scales stay close to the GR line, while the deviation turns on around the transition set by \(m_\ast\).

Glossary of Core Concepts

Large-Scale Structure

Large-scale structure is the cosmic web: a vast network of galaxy filaments, clusters, and voids spread across the Universe.

Gravitational Lensing

Gravitational lensing is the magnifying-glass effect of curved spacetime, which bends the light coming from distant galaxies.

Protected Branch

The protected branch is the narrow safe region where gravity can still be modified without contradicting present measurements, especially the speed of gravitational waves constrained by GW170817.

The \(S_8\) Tension

The \(S_8\) tension is the puzzle that the present Universe appears slightly less clustered than one would expect from the early-Universe signal seen in the cosmic microwave background.

Conceptual Analogy

The Infrared Single-Pole Kernel as a Filter

Think of the kernel as a filter. At short distances it lets standard gravity pass almost unchanged, which is why solar-system and local-galaxy tests remain close to Einstein gravity.

Only beyond a very large cosmic threshold does the response begin to shift. The theory is therefore not an arbitrary deviation everywhere; it is a controlled change that activates around one physical scale.

Observational Roadmap

Chronological Falsification Path

1

DESI

Measure scale-dependent growth through RSD and clustering over wide redshift ranges.

2

Euclid

Cross-check the growth signal with weak lensing and high-precision large-scale mapping.

3

LSST

Expand the statistical power with deep photometric clustering and multi-probe consistency tests.

4

Falsification

Either the same scale survives across probes, or the protected single-pole interpretation is ruled out.

How the Single-Pole Kernel Emerges

This formula is not introduced at random. It is the simplest effective response that survives once one requires a single physical scale, GR recovery on one side of the transition, and a controlled departure on the other.

Shared effective formula
$$ \mu(k,a) = 1 + \beta_0(a)\frac{k^2}{k^2 + a^2 m_\ast^2} $$
Minimal cosmological closure
$$ \mu(k,a) - 1 \;\propto\; \frac{k^2}{k^2 + a^2 m_\ast^2} \qquad \Longrightarrow \qquad \mu(k,a) = 1 + \beta_0(a)\frac{k^2}{k^2 + a^2 m_\ast^2} $$
1

Start from late-time growth

The cosmological side begins by asking how matter growth can deviate from GR at late times without changing the whole background model.

$$ D'' + \left(2+\frac{H'}{H}\right)D' - \frac{3}{2}\Omega_m(a)\,\mu(k,a)\,D = 0 $$
2

Require one physical threshold

The modification is organized around a single transition scale. Below that threshold the response is suppressed; above it the extra growth channel can switch on.

$$ k_{\rm turn}(a) = a\,m_\ast $$
3

Demand controlled limits

The response must return smoothly to the GR baseline in one limit and saturate in a controlled way in the opposite limit, instead of growing without bound.

$$ \mu \to 1 \;\; (k \ll a m_\ast) $$ $$ \mu \to 1+\beta_0(a) \;\; (k \gg a m_\ast) $$
4

Compress to the minimal form

Once those conditions are imposed, the simplest stable two-parameter expression is the single-pole kernel with one amplitude \\(\\beta_0\\) and one scale \\(m_\\ast\\).

$$ \mu(k,a) = 1 + \beta_0(a)\frac{k^2}{k^2 + a^2 m_\ast^2} $$
Spectral reduction to one dominant pole
$$ \mu(k,a) - 1 = \int_0^\infty dm^2\,\rho(m^2,a)\,\frac{k^2}{k^2 + a^2 m^2} $$
$$ \rho(m^2,a)\;\rightarrow\;\beta_0(a)\,\delta\!\left(m^2-m_\ast^2\right) \qquad \Longrightarrow \qquad \mu(k,a) = 1 + \beta_0(a)\frac{k^2}{k^2 + a^2 m_\ast^2} $$
1

Start from a linear response

The quantum or spectral side begins with a propagator-like response, where gravity receives contributions from a spectrum of mediator masses.

$$ G_{\rm ret}(k,a) = \int_0^\infty dm^2\,\frac{\rho(m^2,a)}{k^2 + a^2 m^2} $$
2

Write the Yukawa superposition

Each mass contributes a Yukawa-type factor. The full modification is therefore a weighted superposition rather than a freely invented shape.

$$ \mu(k,a) - 1 = \int_0^\infty dm^2\,\rho(m^2,a)\,\frac{k^2}{k^2 + a^2 m^2} $$
3

Resolve one dominant infrared scale

If one infrared scale dominates the observable response, the full spectrum can be compressed into an effective single resolved pole.

$$ \rho(m^2,a) \approx \beta_0(a)\,\delta\!\left(m^2-m_\ast^2\right) $$
4

Recover the same effective kernel

That reduction again produces the same rational structure, with \\(\\beta_0\\) encoding the net coupling strength and \\(m_\\ast\\) the dominant mass scale.

$$ \mu(k,a) = 1 + \beta_0(a)\frac{k^2}{k^2 + a^2 m_\ast^2} $$
Common conclusion

Different microscopic stories can project onto the same infrared observable. The single-pole kernel is therefore presented here not as an arbitrary guess, but as the minimal effective form consistent with a controlled transition and a dominant resolved scale.

hi_class Integration Sketch

hi_class does not contain a built-in single_pole_kernel switch. A correct implementation must add a custom EFT/Horndeski branch, or an externally validated forward model, whose solved perturbations reproduce the target growth response while keeping \(c_T=1\) and \(\Sigma=1\). The snippets below are therefore an implementation contract, not drop-in public hi_class syntax.

1. Custom Model Parameters

After adding parser support in a private model branch, expose only the analysis-level amplitude and transition scale. Keep the background expansion fixed unless a separate background model is specified.

# Example only after adding these parser entries to a custom branch.
# These are not stock hi_class_public option names.

gravity_model = custom_eft_single_scale
single_scale_beta0 = 0.40
single_scale_mstar = 0.08        # h/Mpc
single_scale_background = lcdm   # fixed background for this test
single_scale_lensing = sigma_one # protected Weyl/lensing sector

# The model module must build alpha_i(a) functions, not assign mu by hand.
# Required protected-sector conditions:
#   alpha_T(a) = 0
#   Sigma_solver(k,a) = 1 within the declared analysis window

2. EFT/Horndeski Model Hook

Do not overwrite G_eff or the slip variables directly inside the generic perturbation solver. Implement a model hook that supplies stable EFT functions; then verify that the solved observables reproduce the target \(\mu(k,a)\), \(\Sigma=1\), and the inverse-slip relation.

// Target response used for validation, not as a direct solver overwrite.
static double mu_target(double k, double a, double beta0, double mstar) {
    double k2 = k * k;
    double am2 = a * a * mstar * mstar;
    return 1.0 + beta0 * k2 / (k2 + am2);
}

static double eta_target(double mu) {
    // Convention: eta = Phi / Psi and Sigma = mu * (1 + eta) / 2.
    return 2.0 / mu - 1.0;
}

// Correct integration point:
//   1. define alpha_T(a)=0;
//   2. choose stable alpha_M(a), alpha_B(a), alpha_K(a) for the custom branch;
//   3. impose the protected-submanifold relation used by the model;
//   4. evolve the full linear system;
//   5. compare solver outputs with the target functions below.

3. Validation Conditions

The implementation is acceptable only if the solver output, after the same scale cuts and windows used in the analysis, satisfies the protected-branch checks.

for each analysis point (k, a):
    mu_ref  = mu_target(k, a, beta0, mstar)
    eta_ref = eta_target(mu_ref)

    require_close(mu_solver(k,a),    mu_ref,  tolerance_mu)
    require_close(Sigma_solver(k,a), 1.0,     tolerance_sigma)
    require_close(eta_solver(k,a),   eta_ref, tolerance_eta)
    require(alpha_T(a) == 0)

# If these checks fail, the run is not an implementation of the protected
# single-scale branch; it is a different EFT/Horndeski model.

The single-pole expression defines the target growth response to be reproduced by a protected-branch implementation in hi_class; it is not a standalone replacement for the perturbation equations.

Technical Note

Protected-Branch hi_class Contract

This note defines what counts as a correct hi_class implementation of the protected single-scale branch. It is an implementation contract, not a claim that the public solver already provides this branch as a stock option.

1. Target response

The implementation must reproduce the resolved growth response \(\mu(k,a)=1+\beta_0\,k^2/(k^2+a^2m_\ast^2)\) within the declared analysis window.

The single-pole formula is therefore the validation target of the model output, not a direct replacement for the perturbation system.

2. Protected-sector conditions

A valid realization must keep \(c_T=1\) and \(\Sigma=1\), and use the same background assumptions adopted in the observational contract unless an enlarged background sector is introduced explicitly.

Under \(\Sigma=1\), the slip target is fixed by \(\eta=2/\mu-1\) in the convention used throughout the site.

3. Solver-level realization

Inside hi_class, the correct route is to define a custom EFT/Horndeski branch whose \(\alpha_i(a)\) functions evolve through the full linear system and then generate the target response as an output.

Directly overwriting \(\mu\), \(G_{\rm eff}\), or slip variables inside the generic perturbation solver is only schematic and does not by itself define a consistent protected-branch implementation.

4. Validation checklist

A run qualifies as a protected single-scale implementation only if it simultaneously reproduces the target \(\mu(k,a)\), preserves \(\Sigma=1\), enforces \(c_T=1\), and satisfies the inverse-slip relation after the same scale cuts and survey windows used in the analysis.

If one of these checks fails, the run may still define an EFT/Horndeski model, but it is not the protected branch tested in the paper.

Code & Reproducibility

For computational researchers: download a minimal Jupyter Notebook that plots the protected kernel across redshift.

import numpy as np
import matplotlib.pyplot as plt

def mu_kernel(k, a, beta0, m_star):
    return 1.0 + beta0 * (k**2 / (k**2 + (a * m_star)**2))

k = np.logspace(-3, 0, 100)
for z in [1.0, 0.5, 0.0]:
    a = 1.0 / (1.0 + z)
    plt.semilogx(k, mu_kernel(k, a, 0.03, 0.05), label=f'z={z}')
plt.legend()
plt.title('Protected single-pole kernel across redshift')
plt.show()
Download Jupyter Notebook (.ipynb)