Thawing Quintessence Model

Updated 31 January 2026

Thawing quintessence is a dynamic dark energy model where a canonical scalar field remains frozen (w ≈ -1) at early times and begins to evolve as Hubble friction decreases.
It employs nearly flat potentials in the slow-roll regime and utilizes parametrizations such as CPL to differentiate its monotonic rise in w from freezing models.
Observations from CMB, BAO, and supernovae constrain its parameters tightly, while theoretical priors and microphysical frameworks ensure consistency with quantum field and Swampland conjectures.

A thawing quintessence model constitutes a dynamically evolving dark energy scenario in which a canonical scalar field, minimally coupled to gravity, remains effectively frozen at early cosmic times due to Hubble friction, with its equation of state $w$ stuck at $w\approx-1$ . Only as the Hubble parameter $H$ drops below the field's effective mass scale does the field begin to evolve, leading to a late-time increase of $w$ away from $-1$ . This signature differentiates thawing models from "freezing" models, where the field rolls during the matter era and $w(a)$ decays toward $-1$ . Thawing dynamics are most commonly analyzed with nearly flat potentials in the slow-roll limit, where the kinetic energy is initially negligible compared to the potential energy. The thawing quintessence paradigm is currently a central focus of both theoretical analysis and observational confrontation, leveraging datasets encompassing the cosmic microwave background (CMB), baryon acoustic oscillations (BAO), Type Ia supernovae, structure growth, the ISW–tSZ cross-correlation, and more (Durrive et al., 2018, Lima et al., 2015, Wolf et al., 2024, Shlivko, 23 Dec 2025, Ibitoye et al., 24 Jan 2026).

1. Theoretical Foundations and Model Structure

Thawing quintessence operates within the canonical scalar-field action: $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ The equation of motion in a flat FLRW background is

$\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$

The corresponding energy density and pressure are

$\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$

with equation-of-state parameter

$w\approx-1$ 0

Thawing models are defined by initial conditions $w\approx-1$ 1 at early times (Hubble suppression), so $w\approx-1$ 2, and only as $w\approx-1$ 3 decreases does the field "thaw", such that $w\approx-1$ 4 increases ( $w\approx-1$ 5) (Durrive et al., 2018, Lima et al., 2015, Pantazis et al., 2016, Tsujikawa, 2013).

A generic analytic form for the thawing equation of state is obtained by Taylor expanding the potential near the initial field value and working in the slow-roll limit ( $w\approx-1$ 6), yielding (Durrive et al., 2018, Tsujikawa, 2013, Chiba et al., 2012): $w\approx-1$ 7 with

$w\approx-1$ 8

$w\approx-1$ 9

where $H$ 0 is the present-day value, $H$ 1 is the current dark-energy fraction, and $H$ 2 parametrizes potential curvature at $H$ 3.

2. Phenomenological Parametrizations and Model Classes

In standard thawing models, the potential is sufficiently flat to ensure $H$ 4 at high redshift, with $H$ 5 remaining near $H$ 6 until the recent epoch. The archetype of thawing potentials is the pseudo–Nambu–Goldstone boson (PNGB, or "hilltop") form

$H$ 7

and quadratic ("chaotic") models: $H$ 8 Different parameterizations of $H$ 9 for thawing models have been developed for robust model-independent analysis. Two-parameter families like the Chevallier–Polarski–Linder (CPL) form: $w$ 0 are convex in $w$ 1 and well-suited for thawing dynamics ( $w$ 2) (Pantazis et al., 2016, García-García et al., 2019). Generalizations, such as the $w$ 3CPL family

$w$ 4

allow interpolation between thawing ( $w$ 5) and freezing ( $w$ 6) classes (Pantazis et al., 2016).

The evolution in the thawing regime is typically monotonic and convex: $w$ 7 increases from $w$ 8 at early times, rising toward less negative values today (but rarely above $w$ 9)(Chiba et al., 2012, Tsujikawa, 2013). Taylor expansions of the quintessence potential up to quadratic order around today, as well as Padé-type approximations for the deviation $-1$ 0, further enable Bayesian model assessment and theoretical prior construction (Shlivko, 23 Dec 2025, García-García et al., 2019).

3. Observational Constraints and Data Analysis

Recent analyses confront thawing quintessence models with CMB (Planck, ACT), BAO (DESI DR2, SDSS, BOSS), SNIa (Pantheon+, JLA, DES-Dovekie, Union3), redshift-space distortions, and cross-correlation measurements (ISW–tSZ) (Ibitoye et al., 24 Jan 2026, Shlivko, 23 Dec 2025, Dinda et al., 21 Apr 2025, Felegary et al., 2024).

Key findings:

With full Planck 2015, JLA SNIa, and BAO data, 95% bounds under a thawing prior are $-1$ 1 and $-1$ 2; $-1$ 3 is not meaningfully constrained ( $-1$ 4 prior), indicating insensitivity to potential curvature (Durrive et al., 2018).
Growth of structure and $-1$ 5 measurements suppress the parameter space for significant deviations from $-1$ 6, but the improvement relative to decade-old studies is modest (Lima et al., 2015).
Bayesian model comparison with theoretically-motivated priors on thawing trajectories reveals that the preference over $-1$ 7CDM is typically modest and strongly contingent on the inclusion of SNe data; DESI+Planck+SNe can yield $-1$ 8, favoring thawing but not decisively (Shlivko, 23 Dec 2025).
ISW–tSZ cross-correlation measurements are consistent with a thawing exponential potential, constraining the slope $-1$ 9, and giving minimal $w(a)$ 0 among standard quintessence scenarios, but all quintessence models remain observationally close to $w(a)$ 1CDM at $w(a)$ 2 (Ibitoye et al., 24 Jan 2026).
Inclusion of data on cosmic curvature shows that physical thawing models do not exhibit phantom crossing (i.e., $w(a)$ 3), in contrast to the behavior allowed by flexible $w(a)$ 4 parameterizations; the preference for phantom crossing in data is a parametrization artifact (Dinda et al., 21 Apr 2025).

A characteristic feature is that most constraints force $w(a)$ 5 close to $w(a)$ 6—e.g., $w(a)$ 7 at 95% confidence (Tsujikawa, 2013, Chiba et al., 2012). Table summarizing key observational bounds (from (Durrive et al., 2018, Chiba et al., 2012, Tsujikawa, 2013)):

Parameter	95% C.L. bounds	Dataset
$w(a)$ 8	$w(a)$ 9	SNe+CMB+BAO
$-1$ 0	$-1$ 1	SNe+CMB+BAO
$-1$ 2	$-1$ 3 (prior)	SNe+CMB+BAO

4. Physical Implications, Degeneracies, and Classification

Thawing models are generically robust to radiative corrections if protected by shift symmetries, as in PNGB or supergravity-motivated constructions (Tsujikawa, 2013). Theoretically, freezing and thawing models can be mapped onto distinct convexity/concavity classes in the $-1$ 4 or $-1$ 5 parameter space (Pantazis et al., 2016, Hara et al., 2017).

Thawing models are more flexible under current observational constraints than most tracker/freezing classes, surviving as a broad swath in the allowed parameter space—particularly in the $-1$ 6 plane—whereas freezing solutions often require fine-tuned or bounce-like solutions to remain viable (Hara et al., 2017).

Dynamical signatures (such as a rapid change in $-1$ 7 or higher-order derivatives) remain weakly constrained due to limited leverage of current data at $-1$ 8, and disentangling thawing from freezing-class solutions at high confidence requires future high-precision probes targeting $-1$ 9 directly (Hara et al., 2017, Takeuchi et al., 2014, Chiba et al., 2012).

Degeneracies remain among $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 0, $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 1, potential curvature, and neutrino mass $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 2; current data only weakly break them (Durrive et al., 2018). Growth rate $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 3 and ISW–tSZ cross-correlation present promising avenues to further delineate thawing from a pure cosmological constant scenario (Ibitoye et al., 24 Jan 2026).

5. Bayesian Model Selection, Priors, and Recent Results

Recent methodology leverages physically motivated priors in $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 4 or Padé-parameter space, often derived from ensembles of microphysical thawing potentials or attractor/Swampland UV constraints (Shlivko, 23 Dec 2025, García-García et al., 2019, Storm et al., 2020). Under such priors, Bayesian evidence for thawing relative to $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 5CDM becomes positive ( $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 6) only when strong supernova datasets are included; growth and BAO data alone do not favor dynamical dark energy (Shlivko, 23 Dec 2025).

Information criteria (AIC, BIC, DIC) generally agree in diagnostic power, but DIC aligns most consistently with Bayesian evidence (Shlivko, 23 Dec 2025). However, $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 7CDM is not decisively ruled out, and current data, even when analyzed with informed priors and alternative parametrizations, confine allowed deviations $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 8– $S = \int d^4x \sqrt{-g} \left[\frac{M_{\mathrm{pl}}^2}{2} R - \frac{1}{2} g^{\mu\nu} \partial_\mu \phi \partial_\nu \phi - V(\phi)\right] + S_\mathrm{matter}$ 9 (Wolf et al., 2024).

Theoretical priors restrict observable $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 0 combinations to a narrow locus—typically $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 1 for slow-roll thawers, $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 2 for fast/hilltop thawers—leaving most of the CPL plane disfavored and highlighting the importance of model-informed analysis (Wolf et al., 2024, García-García et al., 2019, Tsujikawa, 2013).

Recent DESI analyses indicate a statistically competitive fit of thawing trajectories to BAO, CMB, and SNe data, and show that physically consistent thawing models do not display phantom crossing; prior evidence for phantom crossing in $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 3 fits is attributed to parametrization artifacts (Dinda et al., 21 Apr 2025).

6. Microphysical Implications, Swampland Constraints, and Inflation

Thawing potentials can be directly mapped to particle physics frameworks such as supergravity and axion/PNGB-type constructions, where radiative stability of the tiny field mass is ensured by symmetry protection (Tsujikawa, 2013).

Swampland conjectures from string theory impose nontrivial constraints: the refined de Sitter conjecture requires either a sufficiently steep ( $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 4) or tachyonic ( $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 5) potential curvature. Allowed regions in the $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 6 parameter space for slow-roll thawing with current data are $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 7, $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 8 (Storm et al., 2020). For $\ddot\phi + 3H \dot\phi + V_{,\phi} = 0$ 9, initial-field tuning is required to avoid rapid field excursions.

The initial field value for observationally viable thawing scenarios is constrained to be trans-Planckian, $\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$ 0, which—assuming quantum-fluctuation seeding during inflation—implies a minimal inflationary duration $\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$ 1 e-folds (Gupta et al., 2014).

7. Thermodynamic and Cosmological Implications

Thermodynamic analysis reveals that thawing quintessence scenarios generically violate the Generalized Second Law (GSL) of thermodynamics in a finite future. The combined entropy of the apparent horizon and the fluid begins decreasing as the universe continues accelerating, with violation typically manifesting within a few $\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$ 2-folds in the future for representative thawing parameters (Duary et al., 2019). This highlights a tension between thawing phenomenology and fundamental thermodynamic consistency, unless modifications to gravity or field content intervene.

Cosmologically, thawing quintessence produces a late-time, monotonic rise in $\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$ 3, leads to slightly earlier acceleration than $\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$ 4CDM, and suppresses the linear growth rate $\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$ 5 at low redshift by up to 10–16% for the steepest allowed potentials (Felegary et al., 2024, Ibitoye et al., 24 Jan 2026). Most current bounds keep these deviations within a few percent, rendering thawing cosmologically degenerate with a cosmological constant at the background level for $\rho_\phi = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad p_\phi = \frac{1}{2}\dot{\phi}^2 - V(\phi)$ 6.