Structural Limitations on Constraining the Time Evolution of Dark Energy

Cosmological constraints on a time-varying dark energy equation of state are fundamentally limited by the integral structure through which the equation of state enters cosmological observables. We rigorously derive the linear response kernel that maps perturbations in the equation of state ω(z) to comoving distance fluctuations δD(z). By adopting a Fourier mode expansion δω(z) = \sin(kz), we obtain the exact analytic form of the distance response in terms of Sine and Cosine integrals. We show that this mapping involves a double integration over redshift, which acts as an intrinsic low-pass filter with a characteristic \sim k^{-2} scaling in redshift space. This structural limitation is visualized in a schematic diagram and confirmed by observational verification using the full covariance matrix of the Pantheon+ supernova dataset. Our analysis reveals a steep hierarchy in Fisher eigenvalues where the information content drops by an order of magnitude already at the second eigenmode. Consequently, distance-based probes effectively constrain only a single dominant mode of ω(z). This implies that the difficulty in constraining dynamical parameters such as w_a is not due to data precision, but is a necessary consequence of the observable’s integral nature, which renders it structurally blind to the instantaneous rate of change dω/da.

💡 Research Summary

The paper addresses a long‑standing puzzle in cosmology: why distance‑based probes (type Ia supernovae, BAO, CMB angular distances) struggle to constrain a time‑varying dark‑energy equation‑of‑state ω(z) beyond a single effective mode. The authors demonstrate that the difficulty is not a matter of data quality or the choice of parametrisation, but a fundamental structural limitation arising from the double‑integral relationship between ω(z) and observable comoving distances D(z).

Starting from a flat FLRW background, they write ω(z)=−1+δω(z) with |δω|≪1. The continuity equation gives a first‑order density perturbation δρ_DE/ρ_DE=3∫₀^{z} δω(z′)/(1+z′) dz′. Substituting this into the Friedmann equation yields the fractional Hubble perturbation δH/H=(3/2) Ω_DE(z)∫₀^{z} δω(z′)/(1+z′) dz′, i.e. a single integral over the equation‑of‑state history. The comoving distance is defined as D(z)=∫₀^{z} c/H(z′) dz′, so the distance perturbation becomes δD(z)=−∫₀^{z}