⚛️ general relativity

Observational constraints on a damped harmonic oscillator model of dark energy

This paper constrains a damped harmonic oscillator model of dark energy using CMB, BAO, and multiple Type Ia supernova datasets, revealing that the model mimics $\Lambda$ CDM at high redshifts while exhibiting compilation-dependent oscillatory or overdamped behaviors at low redshifts that yield varying $H_0$ values and require a reduced sound speed to resolve numerical stiffness.

Original authors: Saddam Hussain, Muhammad Ahmad, Ming Zhang, Ayodeji Ibitoye, Lang Cui

Published 2026-09-10

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Saddam Hussain, Muhammad Ahmad, Ming Zhang, Ayodeji Ibitoye, Lang Cui

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

The universe is expanding, and for decades, astronomers have watched that expansion speed up. To explain this acceleration, they propose the existence of "dark energy," a mysterious force that pushes galaxies apart. The simplest and most successful idea so far is that this force is a constant, unchanging property of space itself, like a fixed pressure that never varies. However, recent, incredibly precise measurements of the cosmos have hinted that this force might not be so static. Instead of being a steady push, it might be changing over time, perhaps even wobbling or oscillating like a spring. Understanding whether dark energy is a constant or a shifting force is one of the biggest challenges in modern physics, because the answer determines the ultimate fate of the universe and reveals the true nature of the invisible energy that dominates our reality.

A team of researchers recently put a specific, dynamic model of this wobbling energy to the test against the most detailed cosmic maps available. They focused on a theory where dark energy behaves like a damped harmonic oscillator—a system that swings back and forth around a stable point but gradually loses energy and settles down, much like a door closer that swings shut and then stops. In this model, the strength of the push from dark energy is not fixed; it oscillates around a central value before eventually calming down. The scientists wanted to see if the universe's history, as recorded in the light from the early cosmos and the brightness of distant exploding stars, matched this swinging behavior better than the standard, unchanging model.

To investigate, the researchers combined data from three major sources: the cosmic microwave background, which is the afterglow of the Big Bang; measurements of how galaxies are spaced apart, known as baryon acoustic oscillations; and observations of Type Ia supernovae, which serve as cosmic mile markers. They analyzed these datasets using three different catalogs of supernovae, each compiled by different teams using slightly different methods. By feeding the "damped harmonic oscillator" rules into powerful computer simulations, they tracked how the universe would evolve if dark energy were indeed swinging and settling, and then compared those predictions to the actual observations.

The results revealed a fascinating split in the data. When the team used two of the supernova catalogs, the model suggested that dark energy is currently in a state of "underdamping." In this scenario, the energy is still swinging back and forth with noticeable amplitude around its equilibrium point, creating a wobble in the expansion rate that is most visible in the relatively recent universe. However, when they used the third, most extensive catalog of supernovae, the data pointed toward an "overdamped" solution. Here, the energy does not swing at all; it simply drifts smoothly toward its final state without any oscillation. This difference highlights that the specific way we measure and compile the brightness of distant stars can change our interpretation of how dark energy behaves, suggesting that the apparent "wobble" might depend heavily on which dataset we trust most.

Despite these differences, the model showed that any significant deviation from the standard, constant dark energy theory is confined to the very recent past of the universe. For most of cosmic history, the swinging model looks almost identical to the standard model, only diverging noticeably in the last few billion years. The researchers also encountered a technical hurdle: in the regions of the model where the energy swings most rapidly, the computer equations became too difficult to solve, causing the simulations to crash. To get around this, they adjusted a numerical setting related to how the energy clumps together, allowing the simulations to run without changing the physical predictions for the observable universe. This adjustment was a computational fix rather than a change in the laws of physics, ensuring the results remained reliable.

Ultimately, the study found that while the oscillating model is mathematically interesting and fits the data reasonably well, the standard model of a constant dark energy remains the most statistically favored explanation. The data did not provide strong enough evidence to prove that dark energy is swinging, and in fact, the simpler, constant model was preferred by the statistical analysis. However, the study did uncover a crucial insight: the preference for a swinging versus a settling dark energy depends entirely on which catalog of supernovae is used. This suggests that the mystery of dark energy's nature is tightly linked to the details of how we observe the distant universe. The findings do not solve the puzzle of what dark energy is, but they provide a clear framework for testing future observations, showing that as we gather more precise data on the recent universe, we will be able to tell if the cosmic push is a steady hand or a trembling one.

Technical Summary: Observational constraints on a damped harmonic oscillator model of dark energy

Problem Statement
Current cosmological observations, particularly from the Dark Energy Spectroscopic Instrument (DESI) combined with Cosmic Microwave Background (CMB) data, suggest a preference for dynamical dark energy (DE) over the cosmological constant ( $\Lambda$ ) at the $2.6\sigma$ to $3.9\sigma$ level, depending on the Type Ia supernova compilation used. However, the standard phenomenological model for dynamical DE, the CPL ( $w_0 w_a$ CDM) parametrization, is limited by its inability to capture non-monotonic behaviors such as oscillations. Furthermore, previous studies of oscillatory DE models often relied on ad hoc trigonometric functions lacking a common dynamical principle. Additionally, while background-level constraints exist, a full treatment of linear perturbations for rapidly varying equations of state (EOS) crossing the phantom divide ( $w = -1$ ) remains challenging due to numerical instabilities in Boltzmann solvers.

Methodology
The authors constrain a Damped Harmonic Oscillator (DHO) model for the dark energy equation of state, $w_{de}$ , governed by a second-order differential equation with respect to the number of $e$ -folds ( $N$ ):
$\frac{d^2 w_{de}}{dN^2} = -f^2(w_{de} - w_m) + b \frac{dw_{de}}{dN}$
where $f$ is the oscillation frequency, $b$ is the damping rate, and $w_m = -1$ is the equilibrium value. The model is implemented directly into the Boltzmann solver CLASS, integrating the oscillator equation alongside background and linear perturbation equations.

The analysis utilizes:

Full Planck 2018 likelihoods (temperature, polarization, and lensing).
DESI BAO measurements (Data Release II).
Three distinct Type Ia supernova compilations: Pantheon+, DES-Dovekie (a recalibrated DES-SN5YR sample), and Union3.

A critical methodological contribution addresses a numerical pathology: in regions of the $(f, b)$ parameter space characterized by rapid variations of $w_{de}$ at low redshifts, the linear perturbation equations become numerically stiff, preventing integration with the canonical rest-frame sound speed ( $c_s^2 = 1$ ). The authors demonstrate that reducing the sound speed to $c_s^2 \sim 10^{-5}$ alleviates this stiffness without significantly altering cosmological observables (at the $10^{-3}$ level), treating this reduction as a numerical prescription rather than a physical modification.

Key Results

Parameter Space and Dynamics: The posterior distributions for $(f, b)$ $(f, b)$ vary significantly across supernova compilations.
- Pantheon+ favors an overdamped solution (large $b$ ), suppressing oscillations.
- DES-Dovekie and Union3 favor underdamped solutions, exhibiting pronounced oscillatory behavior at low redshifts ( $z \lesssim 0.6$ ).
Equation of State: For the underdamped cases (DES and Union3), the EOS reaches strongly phantom values ( $w_0 < -3.0$ ). For Pantheon+, $w_0$ is constrained to be less negative ( $w_0 \approx -0.5$ ).
Hubble Constant ( $H_0$ ): The model yields different $H_0$ $H_{0}$ values depending on the dataset:
- Pantheon+: $H_0 = 67.53^{+1.22}_{-1.18}$ km/s/Mpc.
- DES-Dovekie: $H_0 = 69.08^{+1.23}_{-1.16}$ km/s/Mpc.
- Union3: $H_0 = 70.67^{+1.89}_{-1.87}$ km/s/Mpc.
  The authors note that the higher $H_0$ values in DES and Union3 arise from a degeneracy with matter density ( $\Omega_{m0}$ ) while keeping the sound horizon ( $r_d$ ) fixed, rather than a modification of pre-recombination physics. Consequently, these results do not constitute a significant alleviation of the Hubble tension with SH0ES.
Model Comparison:
- Chi-Squared: The minimum $\chi^2$ values for the DHO model are intermediate between $\Lambda$ CDM and CPL but generally higher (less favorable) than CPL for all datasets.
- Bayesian Evidence: The logarithmic Bayes factors consistently favor the CPL model over the DHO model (moderate to strong evidence) and $\Lambda$ CDM over the DHO model (very strong evidence for Pantheon+).
Numerical Stability: The reduction of $c_s^2$ to $10^{-5}$ successfully allows the integration of the stiff perturbation equations for the Pantheon+ dataset without introducing significant dark energy clustering or altering the CMB temperature power spectrum ( $C_{\ell}^{TT}$ ) at the observable level.

Significance and Claims
The paper claims that the DHO model provides a physically motivated, dynamical framework for oscillatory dark energy, superior to ad hoc trigonometric parametrizations. Its primary significance lies in demonstrating that late-time cosmological inferences are highly sensitive to the specific supernova compilation used. The model's ability to exhibit distinct dynamical regimes (overdamped vs. underdamped) depending on the dataset highlights potential systematic differences or progenitor-age biases among the supernova samples.

While the DHO model is not statistically preferred over the CPL model or $\Lambda$ CDM by current data, the authors argue it serves as a valuable phenomenological tool. It allows researchers to probe the geometric information encoded in different supernova datasets and test the robustness of the apparent preference for dynamical dark energy. The study concludes that future, larger, and better-controlled supernova samples (e.g., from the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope) are essential to determine whether the inferred oscillatory features reflect genuine late-time dynamics or residual dataset differences.

More like this