Quasi-Stationary State (QSS)

• QSS is a long-lived out-of-equilibrium macroscopic state emerging after rapid violent relaxation and is a stationary Vlasov solution in the mean-field limit.
• It is observed in various classical models such as the HMF and self-gravitating systems, with lifetimes that scale with finite particle corrections and perturbative effects.
• In quantum settings, QSS generalizes to conditioned survival states linked to the spectral properties of quantum Markov semigroups, providing operational measurement insights.

A quasi-stationary state (QSS) is a long-lived state that is stationary only in a qualified sense. In the long-range-interaction literature, it denotes a non-Boltzmann macroscopic state reached after violent, collisionless relaxation and persisting for a time that grows with particle number; in the $N\to\infty$ description, such states are typically stable stationary solutions of the Vlasov equation, whereas for finite $N$ they slowly drift under collisional or graininess effects (Buyl, 2012, Filho et al., 2013). In the quantum-Markov-semigroup setting, the term generalizes the classical quasi-stationary distribution to open quantum dynamics conditioned on survival, with an operational interpretation via direct and indirect quantum measurements and a connection to spectral properties of the semigroup (Dhahri et al., 8 Aug 2025).

1. Definition and conceptual scope

In long-range interacting Hamiltonian systems, QSSs are long-lived out-of-equilibrium macroscopic states formed after an initial rapid collective or violent relaxation and before eventual Boltzmann-Gibbs equilibration. The common formulation across the HMF, self-gravitating, spin, and spherical models is that the collisionless $N\to\infty$ dynamics is governed by the Vlasov equation, and the relevant QSS is a stable stationary solution of that equation; finite-$N$ effects then induce a secular evolution on a much longer timescale (Patelli et al., 2011, Gupta et al., 2010, Gupta et al., 2013).

This long-range usage already contains an important distinction between exact and quasi stationarity. A Vlasov-stationary distribution is exactly stationary in the mean-field limit, but only quasi-stationary for finite $N$, because finite-size corrections eventually drive relaxation. This is why the literature repeatedly characterizes QSSs as long-lived remnants of stable stationary Vlasov states rather than as exact finite-$N$ equilibria (Filho et al., 2013, Buyl, 2012).

In the quantum-Markov-semigroup setting represented here, the supplied background describes a QSS as the quantum analogue of a quasi-stationary distribution for a process with absorption. If $\mathcal S_t$ is the survival semigroup, the survival probability is

$p_\rho(t)=\operatorname{Tr}[\mathcal S_t(\rho)],$

and the conditioned state is

$$\widetilde \rho_t=\frac{\mathcal S_t(\rho)}{\operatorname{Tr}[\mathcal S_t(\rho)]}.$$

A state $\rho_{\mathrm{QSS}}$ is quasi-stationary when conditioning on survival leaves it invariant, equivalently when

$N$0

for some decay rate $N$1. The supplied abstract further states that QSSs in this setting admit an operational interpretation via direct and indirect quantum measurements and are connected to spectral properties of the semigroup (Dhahri et al., 8 Aug 2025).

2. Collisionless mean-field origin

The basic mean-field mechanism is clearest in the Vlasov framework. For a broad class of one-dimensional long-range Hamiltonians,

$N$2

the empirical one-particle measure converges, in the large-$N$3 limit, to a smooth distribution $N$4 solving a Vlasov equation. In this description, a QSS is an initial distribution $N$5 satisfying the stationary unperturbed Vlasov equation and, crucially, linearly stable under Vlasov dynamics (Patelli et al., 2011).

The HMF model provides the canonical example. With interaction $N$6, the order parameter is the magnetization

$N$7

Homogeneous QSSs are distributions $N$8 uniform in $N$9, for which the unperturbed Liouvillian reduces to free streaming. In this regime, linear response around the QSS yields a Vlasov analogue of a Kubo formula, and the dielectric function $N\to\infty$0 determines whether the response grows, decays, or oscillates; for a stable QSS one must have $N\to\infty$1 (Patelli et al., 2011).

A recurrent conclusion in this literature is that the collisionless stage should not be confused with large inter-particle correlation buildup. One explicit interpretation is that the Vlasov equation arises in the Kac-scaled $N\to\infty$2 limit and that dynamically generated correlations remain of order at most $N\to\infty$3; finite-$N\to\infty$4 deviations are then attributed mainly to the secular accumulation of small collisional corrections rather than to macroscopically large correlations (Filho et al., 2013).

3. Representative classical realizations

Several distinct long-range models realize the same QSS logic while differing sharply in geometry, observables, and stability thresholds. In the anisotropic mean-field Heisenberg model of classical spins,

$N\to\infty$5

a nonmagnetic waterbag is an exact stationary Vlasov solution. Linearization yields

$N\to\infty$6

hence the threshold

$N\to\infty$7

For $N\to\infty$8, the state is Vlasov-unstable and the magnetization grows as $N\to\infty$9, implying $N$0. For $N$1, it is Vlasov-stable and the QSS lifetime scales numerically as $N$2 (Gupta et al., 2010).

For particles moving on a sphere with an attractive infinite-range Heisenberg-like interaction,

$N$3

the non-magnetized waterbag state is Vlasov-stationary and linearly stable for

$N$4

In this isotropic case the QSS lifetime again scales as $N$5, whereas the unstable regime yields $N$6. With global anisotropy $N$7, the thresholds shift to

$N$8

and the QSS window widens (Gupta et al., 2013).

The one-dimensional self-gravitating sheet model places the emphasis on violent relaxation and coarse-grained structure. For one-level waterbag initial conditions, Lynden-Bell predictions are reasonably good in the low-energy region, but at higher energies they are generally not even qualitatively correct. When they fail, the resulting QSS is characterized by a degenerate core-halo structure, which the paper associates with dynamical resonances. The order parameters

$N$9

are used to quantify non-separability and to distinguish QSSs from thermal equilibrium (Joyce et al., 2010).

Three-dimensional self-gravitating collapse adds another layer: the QSS is not universal but depends on the formation path, especially on the amount of energy exchanged during relaxation. Violent top-down collapse produces a flat central core and an outer density profile $N$0 with radially elongated orbits, whereas less violent bottom-up dynamics produces weaker anisotropy and density profiles well fitted by the Navarro-Frenk-White form. The paper’s practical test of collisionless quasi-equilibrium is the Jeans diagnostic

$N$1

rather than a claim of thermodynamic equilibrium (Labini et al., 2020).

4. Robustness, perturbations, and generalized QSSs

A major refinement of QSS theory concerns the conditions under which such states remain robust in the mean-field limit. For particle systems with attractive power-law interactions

$N$2

the decisive criterion is not integrability of the potential, but integrability of the force at large distances. The paper argues that robust QSSs exist only for

$N$3

because then the collisionality parameter $N$4 as $N$5. For $N$6, QSSs are not robust in the unregularized model; if they appear, their existence is strongly conditioned by short-range regularization (Marcos et al., 2017).

Dissipation does not automatically destroy QSS phenomenology. For long-range systems with viscous damping or quasi-elastic inelastic collisions, the relevant object is a scaling quasi-stationary state rather than a stationary state in the original variables. The central ansatz is

$N$7

with virial scaling

$N$8

In this regime the system remains quasi-stationary in rescaled variables, and the paper emphasizes that the velocity distributions never show any tendency to evolve towards a Maxwell-Boltzmann form (Joyce et al., 2013).

Weak local perturbations can also select rather than destroy QSSs. In the one-dimensional self-gravitating sheet model with internal local perturbations acting at collisions, very different initial QSSs are driven toward a unique non-equilibrium stationary state that is itself a QSS of the unperturbed long-range dynamics. If the perturbation is removed after this state is reached, the system remains there. This suggests a robustness of the QSS framework at the level of motion through a manifold of QSSs, even when the original unperturbed QSS is not preserved (Joyce et al., 2016).

5. Statistical, kinetic, and dynamical interpretations

The principal theoretical approaches differ in what they assume about mixing, stationarity, and resonant structure. Lynden-Bell theory treats violent relaxation as an entropy-maximization problem subject to Vlasov constraints. For one-level waterbag initial data in the sheet model, the predicted coarse-grained state has the Fermi-Dirac-like form

$N$9

and the corresponding entropy is

$N$0

This framework remains the broadest “ensemble view” of HMF QSSs and of the associated out-of-equilibrium phase diagram, but it presupposes sufficiently efficient mixing and stationary coarse-grained outcomes (Joyce et al., 2010, Buyl, 2012).

A distinct line of work uses exact stationary Vlasov states of the form $N$1, where

$N$2

for HMF-type models. This BGK-like perspective captures some dynamically selected stationary regimes, especially at sufficiently large initial magnetization, but it does not resolve oscillatory QSSs and is not fully predictive of the function $N$3 selected by the dynamics (Buyl, 2012).

Core-halo theory starts instead from the empirical observation that violent relaxation often produces a dense low-energy core plus a dilute high-energy halo. One widely used ansatz is

$N$4

with $N$5. A more refined formulation for the HMF model describes the QSS as a superposition of two independent Lynden-Bell distributions, one for the core and one for the halo. The review’s central claim is that this “double Lynden-Bell” structure captures low-energy core-halo QSSs better than a single Lynden-Bell distribution, while still showing that complete collisionless relaxation is usually not achieved (Konishi, 2016).

No single framework gives a full account of the numerically observed dynamics. A recurring source of disagreement is that QSSs may be oscillatory rather than stationary, so observable choice matters. In particular,

$N$6

and therefore $N$7, $N$8, and $N$9 need not yield the same apparent phase diagram when the QSS oscillates. This is one reason the review concludes that several theories exist but none yet provides a full account of the dynamics found in simulations (Buyl, 2012).

6. Spectral and quantum extensions, and terminological scope

A spectral interpretation of QSSs appears in both stochastic and quantum settings. In a conformally invariant one-dimensional stochastic model with a non-local perturbation, part of the finite-size scaling spectrum remains unchanged while some levels go exponentially to zero with system size. The resulting QSSs have relaxation times that grow exponentially with $\mathcal S_t$0, and the dominant QSS reproduces, in the finite-size scaling limit, the same spatial properties as the stationary state of the original conformal model (Alcaraz et al., 2011).

For quantum many-body systems with power-law couplings, the proposed mechanism is even sharper. The paper argues that when $\mathcal S_t$1, the spectrum remains discrete in the thermodynamic limit. As a consequence, the conventional dephasing logic based on continuous spectral components fails, Poincaré recurrence times remain finite, and long-lived metastable behavior can be traced to this pure-point spectral structure. In the long-range Ising and spherical models, this is used to interpret anomalous magnetization dynamics and the failure of observable equilibration as manifestations of QSS behavior (Defenu, 2020).

The quantum-Markov-semigroup literature represented here extends the term in a different direction. The abstract of “Quasi-stationary normal states for quantum Markov semigroups” states that QSSs generalize quasi-stationary distributions for classical Markov chains, have an operational interpretation through direct and indirect quantum measurements, and are connected to spectral properties of the semigroup (Dhahri et al., 8 Aug 2025). A computationally distinct use appears in the work on quantum state preparation, where a QSS is a pure many-body state with almost all spectral weight in a narrow energy window,

$\mathcal S_t$2

so that it behaves approximately as an eigenstate up to times $\mathcal S_t$3. The proposed algorithm prepares such states in polynomial time using quantum signal processing and quantum singular value transformations (Garratt et al., 2024).

The acronym also has distinct meanings outside the quasi-stationary-state literature proper. In power systems, QSS often denotes a quasi steady-state model or a quasi steady-state frequency, not a long-lived metastable state (Wang et al., 2014, Gutierrez-Florensa et al., 27 May 2025). In chemical kinetics and reaction-network theory, QSS denotes quasi-steady-state reduction, where one sets selected derivatives to zero and studies when the resulting reduced system agrees with a Tikhonov–Fenichel singular perturbation reduction (Feliu et al., 2019). This terminological spread suggests that “QSS” names a family of reduced or conditioned stationarity concepts rather than a single universally shared object.