Chevallier-Polarski-Linder (CPL) Model

• CPL is a two-parameter phenomenological model of dark-energy evolution defined by w₀ and wₐ, capturing the present value and first-order slope around today.
• It provides a closed-form expression for dark-energy density, facilitating comparisons with BAO, supernova, and CMB observations.
• CPL serves as a baseline in cosmology to project more elaborate models and assess tensions like H₀ and neutrino-mass bounds, despite its high-redshift extrapolation limitations.

The Chevallier–Polarski–Linder (CPL) parametrization is a two-parameter phenomenological description of a time-varying dark-energy equation of state, conventionally written as

$w(a)=w_0+w_a(1-a),$

or equivalently,

$w(z)=w_0+w_a\frac{z}{1+z}.$

In this form, $w_0$ is the present-day value of the dark-energy equation of state and $w_a$ controls its first-order evolution away from the present epoch. Across the recent literature, CPL functions as a standard low-dimensional language for late-time expansion history, a baseline dynamical-dark-energy model, and a common projection space in which more elaborate scenarios are compared with geometric data such as BAO, supernovae, and CMB distance information (Barua et al., 22 Aug 2025).

1. Formal definition and kinematical meaning

CPL is most commonly introduced as a first-order expansion around today. In scale factor,

$w(a)=w_0+w_a(1-a),$

while in redshift,

$w(z)=w_0+w_a\frac{z}{1+z}.$

The two forms are equivalent through $a=(1+z)^{-1}$. The present epoch corresponds to $a=1$ and $z=0$, so $w(1)=w(0)=w_0$, while the formal high-redshift limit is

$w(z)=w_0+w_a\frac{z}{1+z}.$0

The $w(z)=w_0+w_a\frac{z}{1+z}.$1CDM limit is the special point

$w(z)=w_0+w_a\frac{z}{1+z}.$2

This is the standard interpretation used in observational analyses, critiques of low-order truncations, and model-comparison studies (Nesseris et al., 28 Mar 2025).

Several papers emphasize why this form remains pervasive. It has a 2-dimensional phase space, exhibits linear behaviour in redshift, and has a straightforward physical interpretation. In that sense, CPL is not usually presented as a microphysical theory; it is a compact phenomenological basis for late-time dark-energy inference (Barua et al., 22 Aug 2025). A recurring interpretation in the literature is

$w(z)=w_0+w_a\frac{z}{1+z}.$3

so $w(z)=w_0+w_a\frac{z}{1+z}.$4 encodes the present-day equation of state and $w(z)=w_0+w_a\frac{z}{1+z}.$5 the local slope around today (Giarè et al., 2024).

A frequent misconception is to treat CPL as if it were model-independent in a strong sense. The literature instead describes it as a truncated Taylor-like ansatz around $w(z)=w_0+w_a\frac{z}{1+z}.$6. One explicit higher-order extension is

$w(z)=w_0+w_a\frac{z}{1+z}.$7

for which CPL is the restriction $w(z)=w_0+w_a\frac{z}{1+z}.$8. This framing matters because CPL’s simplicity comes from setting higher derivatives to zero rather than measuring them directly (Nesseris et al., 28 Mar 2025).

2. Background evolution and the dark-energy density

In flat FLRW analyses, CPL enters the background through the continuity equation for dark energy. A standard expression used repeatedly is

$w(z)=w_0+w_a\frac{z}{1+z}.$9

For the CPL form this integrates to

$w_0$0

or, in redshift,

$w_0$1

This closed form is one of the practical reasons the parametrization remains standard in late-time cosmology (Nesseris et al., 28 Mar 2025).

The corresponding background expansion is commonly written as

$w_0$2

for a flat universe, or equivalently through the integral form

$w_0$3

This structure underlies distance-redshift calculations, BAO observables, luminosity distances, and compressed CMB constraints (Giarè et al., 2024).

A major interpretive point in recent work is that cosmological observables depend directly on the dark-energy density rather than directly on $w_0$4. One paper makes this criticism explicit, noting that the mapping from $w_0$5 to $w_0$6 is nonlinear and that the expansion history is only weakly sensitive to $w_0$7. The sensitivity kernel quoted there is

$w_0$8

This is used to argue that broad degeneracies in the $w_0$9 plane are structural rather than merely accidental properties of current data (Montefalcone et al., 26 Mar 2026).

3. Observational role and representative constraints

CPL is often the baseline evolving-dark-energy model in DESI-era analyses. In one broad comparison of two-parameter dark-energy parametrizations, the reported 68% constraints for CPL were

$w_a$0

for Planck+DESI+PantheonPlus, and

$w_a$1

for Planck+DESI+DESY5, with both combinations favoring $w_a$2 and $w_a$3 (Giarè et al., 2024). A later combined DES Y3 $w_a$4pt + DESI DR2 BAO + DES Y5 SN + Planck analysis reported

$w_a$5

again placing CPL in the now-familiar DESI-preferred region (Tsedrik et al., 19 Dec 2025).

Representative CPL constraints reported in recent analyses are summarized below.

Data combination	Reported CPL constraint	Source
Planck+DESI+PantheonPlus	$w_a$6	(Giarè et al., 2024)
Planck+DESI+DESY5	$w_a$7	(Giarè et al., 2024)
DES+DESI+SN+CMB	$w_a$8	(Tsedrik et al., 19 Dec 2025)

These contours are often interpreted as favoring a present-day quintessence-like state with evolution toward the phantom regime at earlier times. In the Planck+DESI+PantheonPlus case, the quoted CPL crossing redshift was

$w_a$9

while for Planck+DESI+DESY5 it was

$w(a)=w_0+w_a(1-a),$0

This yields the specific narrative of a recent crossing of $w(a)=w_0+w_a(1-a),$1 near $w(a)=w_0+w_a(1-a),$2 in that phenomenological basis (Giarè et al., 2024).

CPL also matters in parameter-inference problems not primarily aimed at dark energy itself. In a study of neutrino-mass bounds, it served as the standard dynamical-dark-energy extension relative to $w(a)=w_0+w_a(1-a),$3CDM and a second-order EXP variant. The main conclusion there was that CPL generally yields slightly tighter neutrino-mass upper bounds than EXP, while both relax the bounds relative to $w(a)=w_0+w_a(1-a),$4CDM (Barua et al., 22 Aug 2025). In another late-time study using Bayesian Physics-Informed Neural Networks, allowing CPL dynamics shifted the posterior of $w(a)=w_0+w_a(1-a),$5 upward relative to $w(a)=w_0+w_a(1-a),$6CDM, and the CPL+$w(a)=w_0+w_a(1-a),$7 case reduced tension levels with Planck and SH0ES into the $w(a)=w_0+w_a(1-a),$8-$w(a)=w_0+w_a(1-a),$9 range for several dataset combinations (Yarahmadi et al., 1 Jan 2026).

4. Critiques, degeneracies, and known limitations

The sharpest recent critique is statistical rather than algebraic. One analysis argues that truncating the dark-energy equation of state at CPL order, instead of allowing higher-order terms and then marginalizing over them, injects information not present in the data. In that study, standard CPL could suggest several-$w(z)=w_0+w_a\frac{z}{1+z}.$0 evidence against the cosmological-constant point $w(z)=w_0+w_a\frac{z}{1+z}.$1, but once one more Taylor coefficient was included and marginalized over, the evidence dropped below $w(z)=w_0+w_a\frac{z}{1+z}.$2; with two higher-order coefficients marginalized over, it dropped below $w(z)=w_0+w_a\frac{z}{1+z}.$3 (Nesseris et al., 28 Mar 2025). A central statement there is that fixing $w(z)=w_0+w_a\frac{z}{1+z}.$4 is not statistically equivalent to allowing them to vary.

A second limitation concerns high-redshift extrapolation. In a DESI DR2 BAO-only assessment of flat $w(z)=w_0+w_a\frac{z}{1+z}.$5CDM, CPL fits to different BAO observable subsets pulled the parameters in different directions. The reported combined BAO-only constraint was

$w(z)=w_0+w_a\frac{z}{1+z}.$6

and the paper interpreted such behavior as “unphysical parameter compensation at high redshifts,” with elevated $w(z)=w_0+w_a\frac{z}{1+z}.$7 and very negative $w(z)=w_0+w_a\frac{z}{1+z}.$8 compensating for CPL’s limitations when high-redshift sensitivity matters (Lee, 23 Jun 2025).

A third line of criticism concerns physical mapping. Exact CPL reconstruction into minimally coupled quintessence or barotropic fluids does not uniformly span broad physical model space. One reconstruction study found that over most of the $w(z)=w_0+w_a\frac{z}{1+z}.$9 plane, CPL maps onto a fairly narrow form of quintessence potential $a=(1+z)^{-1}$0, while for barotropic models the functional dependence $a=(1+z)^{-1}$1, up to a multiplicative constant, depends only on

$a=(1+z)^{-1}$2

and not on $a=(1+z)^{-1}$3 and $a=(1+z)^{-1}$4 separately (Scherrer, 2015). This suggests that the nominal two-parameter CPL plane does not translate cleanly into two independent physical degrees of freedom for those model classes.

Related work has asked whether CPL’s strong early-time phantom behavior is observationally robust. One such study constructed alternative two-parameter families that match CPL near $a=(1+z)^{-1}$5 but soften its asymptotic phantomness, then found that current data cannot decisively distinguish those alternatives from CPL. That result supports a limited but important claim: CPL is robust as a low-redshift effective description, while its more dramatic high-redshift extrapolations are not yet uniquely demanded by late-time data (Artola et al., 5 Oct 2025).

5. Extensions, embeddings, and nonstandard uses

A large fraction of current work uses CPL less as a theory than as a benchmark language into which other dynamics are projected. One explicit extension is the “EXP” model,

$a=(1+z)^{-1}$6

described as a second-order correction to CPL that preserves the same two-parameter space $a=(1+z)^{-1}$7 (Barua et al., 22 Aug 2025). Closely related analyses considered a generalized exponential parent model whose first-order truncation is CPL, with higher-order truncations denoted Ext2 and Ext3. Those studies found that current data constrain these extensions similarly to CPL at the background level, while perturbative residuals reveal the nonlinear corrections more clearly (Pan et al., 2019).

CPL also appears inside interacting dark-energy models. In one analytic treatment, the intrinsic dark-energy equation of state was taken as

$a=(1+z)^{-1}$8

and embedded into the interaction models

$a=(1+z)^{-1}$9

For the $a=1$0 case, the dark-energy density became

$a=1$1

and the analysis concluded that this interacting CPL scenario modestly improves the fit by AIC relative to non-interacting CPL, while BIC continues to favor $a=1$2CDM (Neumann et al., 24 Apr 2026).

A different use appears in models where CPL describes an effective fluid rather than fundamental dark energy. In Buchert backreaction cosmology, the effective equation of state of the backreaction-induced fluid was assumed to take the form

$a=1$3

and the resulting averaged Hubble law was fitted to SN Ia and OHD data (Yao et al., 2024). In Extended Uncertainty Principle cosmology, CPL was used as an observational bridge between an exact Lambert-$a=1$4-based effective equation of state and the standard $a=1$5 plane. There the authors explicitly defined

$a=1$6

and used CPL only to “facilitate comparison with observational constraints” rather than to derive the dynamics (Roushan et al., 26 Nov 2025).

CPL has likewise been used as the shared background sector in interacting dark scattering, as the equation of state inside 4D Einstein–Gauss–Bonnet dark-energy phenomenology, and as the late-time evolving-dark-energy component in Hubble-tension studies with free neutrino mass (Tsedrik et al., 19 Dec 2025, Mukherjee et al., 7 Jul 2025, Yarahmadi et al., 1 Jan 2026). Across these applications, a consistent theme is that CPL supplies a common parameter plane for comparison, even when the underlying theory is structurally different.

6. Diagnostics, stability, and computational practice

Because many dark-energy parametrizations are nearly degenerate in $a=1$7 and $a=1$8, higher-order diagnostics have been used to discriminate CPL from alternatives. One Statefinder-hierarchy study concluded that $a=1$9 can differentiate CPL and JBP from $z=0$0CDM and Padé parameterizations, while

$z=0$1

has greater discriminatory power. In that work, ordinary background quantities such as $z=0$2 and $z=0$3 were found to be too degenerate for clean separation, whereas higher derivatives of the scale factor were more informative (Qi et al., 2015).

Thermodynamic analyses have also been applied to CPL. In a flat FRW universe with dark energy and cold dark matter treated as perfect fluids, the CPL dark-energy equation of state was written as

$z=0$4

leading to

$z=0$5

The stability criterion was then formulated in terms of the total equation-of-state parameter $z=0$6 and $z=0$7. A central conclusion was that phantom total behavior,

$z=0$8

is thermodynamically unstable, and that CPL stability in the examples considered favored negative $z=0$9 when $w(1)=w(0)=w_0$0 (Mamon et al., 2018).

On the computational side, CPL’s closed-form density still leaves many observables numerically defined. For luminosity distances in flat CPL cosmology,

$w(1)=w(0)=w_0$1

one numerical study compared Romberg integration with composite Gauss–Legendre quadrature and found the latter substantially faster for comparable accuracy over large redshift ranges (Yue et al., 2011). More recent inference work has moved in the opposite direction, using Bayesian Physics-Informed Neural Networks to reconstruct $w(1)=w(0)=w_0$2 while enforcing the CPL Friedmann equation as a soft constraint (Yarahmadi et al., 1 Jan 2026).

Taken together, these analyses support a narrow but durable conclusion. CPL remains central because it is simple, interpretable, and directly connected to late-time observables; yet its parameters are not themselves fundamental physical invariants, its high-redshift extrapolation is fragile, and its strongest empirical use is often as a common phenomenological interface across otherwise distinct cosmological models.