Skip to main content

Three Model Families

This page defines the three flat-universe models fitted to Pantheon+SH0ES + DESI DR2 BAO: ΛCDM (2 free parameters), CPL (4) and Ξosc (6). All three share H₀ and Ωₘ and are sampled with the same top-hat priors on those parameters.


Expansion histories​

In a flat universe, E(z) = H(z)/H₀. Matter is the only other component, so only the dark-energy density term changes between models.

ModelwE²(z)
ΛCDM−1-1Ωm(1+z)3+(1−Ωm)\Omega_m (1+z)^3 + (1-\Omega_m)
CPLw0+wa z1+zw_0 + w_a\,\dfrac{z}{1+z}Ωm(1+z)3+(1−Ωm)(1+z)3(1+w0+wa)exp⁡ ⁣[−3waz1+z]\Omega_m (1+z)^3 + (1-\Omega_m)(1+z)^{3(1+w_0+w_a)}\exp\!\left[-\dfrac{3 w_a z}{1+z}\right]
Ξoscw0,mean+εcos⁡ ⁣(2πtT+φ0)w_{0,\mathrm{mean}} + \varepsilon\cos\!\left(\dfrac{2\pi t}{T} + \varphi_0\right)Ωm(1+z)3+(1−Ωm) ρDE(z)ρDE(0)\Omega_m (1+z)^3 + (1-\Omega_m)\,\dfrac{\rho_{\mathrm{DE}}(z)}{\rho_{\mathrm{DE}}(0)} with ρDE(z)ρDE(0)=exp⁡ ⁣[3∫0z1+w1+z′ dz′]\dfrac{\rho_{\mathrm{DE}}(z)}{\rho_{\mathrm{DE}}(0)} = \exp\!\left[3\int_0^z \dfrac{1+w}{1+z'}\,dz'\right]

How Ξosc is implemented​

The oscillation is a function of cosmic time t, but the likelihood is evaluated in redshift. The code (XiOscillator) therefore builds, once, a redshift grid and the cosmic age t(z) on that grid from a fixed reference cosmology — flat ΛCDM with H₀ = 70 and Ωₘ = 0.3 (astropy) — then evaluates w on it, integrates the density ratio with a trapezoidal rule and interpolates. The time axis is not recomputed from the cosmology being sampled; this is limitation 2 on Known Limitations.

Two consequences worth knowing

w(z=0) is not w₀,mean. The present-day value is w₀,mean + ε cos(2π t₀/T + φ₀), which depends on the phase. The w0_mean parameter is the average level around which w oscillates.

Ξosc does not nest ΛCDM. Setting ε → 0 gives a constant w = w₀,mean (a wCDM model), which equals ΛCDM only if w₀,mean = −1 as well. Neither is Ξosc nested in CPL. AIC and BIC remain usable, but a likelihood-ratio reading of "ε = 0" is not a test of ΛCDM.


Free parameters and priors​

All priors are uniform inside the bounds below and zero outside (open intervals; log_prior returns −∞ outside). They are read from the sampling script — see the Prior Specification page for the discussion of why they matter.

ParameterFirst used inPrior (uniform, open interval)Unit
H₀ΛCDM50 < x < 100km/s/Mpc
ΩₘΛCDM0.1 < x < 0.5—
w₀CPL−2 < x < 0—
wₐCPL−3 < x < 2—
w₀,meanΞosc−1.5 < x < −0.5—
εΞosc0 < x < 0.3—
TΞosc1 < x < 10Gyr
φ₀Ξosc−π < x < πrad

Why these three​

PropertyΛCDMCPLΞosc
Dark-energy behaviourconstantsmooth, monotonic in zperiodic in cosmic time
Rolebaselinestandard dynamical comparatorexploratory probe
Extra parameters vs ΛCDM—w₀, wₐw₀,mean, ε, T, φ₀
Weakly constrained parameters—wₐT, φ₀ (and ε near its bounds)