Kolmogorov-type Non-Thermal Fixed Points

• Kolmogorov-type non-thermal fixed points are far-from-equilibrium states characterized by universal power laws and self-similar scaling behavior.
• They bridge quantum turbulence and classical fluid dynamics, linking features like quantized vortices and acoustic excitations through directed cascades.
• Experimental and theoretical studies reveal cascade-driven energy and particle transport, establishing these fixed points as robust attractors in many-body dynamics.

Searching arXiv for relevant papers on Kolmogorov-type non-thermal fixed points. Kolmogorov-type non-thermal fixed points are far-from-equilibrium scaling states in which a many-body system exhibits universal, non-thermal power laws together with directed transport of conserved quantities across scales. In dilute Bose gases, the concept was made concrete by relating non-thermal fixed points to superfluid turbulence, quantized vortices, and acoustic excitations, with the single-particle momentum distribution $n(k)\sim k^{-\zeta}$ serving as a central diagnostic (Nowak et al., 2010). In this framework, the fixed point is “non-thermal” because the distribution is not of Rayleigh–Jeans or Bose–Einstein form, yet it is “fixed” in the sense of stationary or quasi-stationary scaling under the nonequilibrium evolution equations. The qualifier “Kolmogorov-type” refers to the structural analogy with turbulence theory: universal scaling laws emerge from conservation laws, dimensionality, and scale-to-scale transport, and in suitable observables the incompressible kinetic-energy spectrum exhibits the Kolmogorov law $E(k)\sim k^{-5/3}$ (Nowak et al., 2010).

1. Definition and scaling characterization

A non-thermal fixed point is a metastable far-from-equilibrium configuration of an isolated many-body system characterized by stationary scaling of correlation functions and momentum distributions in some infrared or ultraviolet regime, together with self-similarity of the evolution near the fixed point (Nowak et al., 2013). For a Bose gas, the central observable is the single-particle momentum distribution

$$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$

and its radial average $n(k,t)$ with $k=|\mathbf{k}|$ (Nowak et al., 2010). The fixed-point condition is expressed by the scaling law

$n(sk) = s^{-\zeta} n(k),$

equivalently $n(k)\sim k^{-\zeta}$ over a finite range of momenta (Nowak et al., 2010).

In a broader dynamical formulation, non-thermal fixed points are associated with self-similar spatio-temporal scaling,

$$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$

or, equivalently in related formulations,

$n(k,t)=t^\alpha f(t^\beta k),$

with universal exponents $\alpha$ and $E(k)\sim k^{-5/3}$0 constrained by conservation laws (Madeira et al., 2022). In Bose-gas experiments and simulations, this scaling is interpreted as an attractor solution in the space of nonequilibrium states, closely analogous to renormalization-group fixed points in equilibrium critical phenomena, but reached from strongly out-of-equilibrium initial conditions rather than by fine tuning control parameters (Madeira et al., 2022).

Kolmogorov-type behavior enters when the scaling is understood as the manifestation of a cascade. In this sense, a Kolmogorov-type non-thermal fixed point is a stationary or quasi-stationary scaling regime maintained by a constant directional flux of a conserved quantity, typically particle number or energy, across scales (Madeira et al., 2022). This definition encompasses both classical and quantum turbulence, while distinguishing non-thermal fixed points from thermal fixed points, where detailed balance eliminates net fluxes and the distribution reduces to a Gibbsian form (Nowak et al., 2013).

2. Field-theoretic origin and universality

The modern formulation of non-thermal fixed points was developed in nonequilibrium quantum field theory using the two-particle irreducible effective action, Kadanoff–Baym equations, and later functional renormalization group methods (0809.5208). In this language, non-thermal fixed points are stationary solutions of the full dynamical equations for statistical and spectral correlators that do not satisfy a fluctuation–dissipation relation and therefore lie outside thermal equilibrium (0809.5208). For an $E(k)\sim k^{-5/3}$1 scalar theory, the functional renormalization group reveals a hierarchy of fixed-point solutions that includes vacuum, thermal equilibrium, and nonequilibrium fixed points, with the nonequilibrium solutions admitting scaling exponents that explicitly depend on spatial dimension $E(k)\sim k^{-5/3}$2 (0809.5208).

For dilute Bose gases, the nonperturbative 2PI analysis yields a central infrared prediction: for constant radial quasiparticle flux, the occupation spectrum scales as

$E(k)\sim k^{-5/3}$3

where $E(k)\sim k^{-5/3}$4 is the spatial dimension (Nowak et al., 2010). In contrast, thermal Rayleigh–Jeans scaling for quadratic dispersion gives $E(k)\sim k^{-5/3}$5 in the ultraviolet (Nowak et al., 2010). The infrared exponent $E(k)\sim k^{-5/3}$6 is therefore the characteristic signature of a non-thermal rather than thermal fixed point. The same exponent appears in relativistic scalar field models, which supports universality across relativistic and nonrelativistic systems (Nowak et al., 2010).

A complementary low-energy effective field theory for multicomponent $E(k)\sim k^{-5/3}$7 Bose gases integrates out density fluctuations and describes interacting Goldstone modes of total and relative phases. In the large-$E(k)\sim k^{-5/3}$8 limit, this theory predicts the same self-similar scaling exponents associated with the non-thermal fixed point and identifies a Gaussian fixed point dominated by relative phases, while anticipating a non-Gaussian fixed point for evolution with a distinctly lower power of time (Mikheev et al., 2018). This suggests that Kolmogorov-type non-thermal fixed points are not tied to a single microscopic derivation, but recur across kinetic, effective-theory, and renormalization-group descriptions.

3. Cascades, spectra, and Kolmogorov-type laws

The canonical turbulent interpretation starts from continuity equations in momentum space for radial particle and energy densities. In weak wave turbulence, constant particle and energy fluxes define Kolmogorov–Zakharov spectra; for a $E(k)\sim k^{-5/3}$9-dimensional Bose gas with quadratic dispersion, the weak-turbulence ultraviolet exponents are

$$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$0

while the strong-turbulence infrared exponents become

$$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$1

with dynamical exponent $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$2 (Nowak et al., 2013). For the nonrelativistic Bose gas, $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$3, so the infrared particle-cascade exponents are $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$4 in $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$5 and $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$6 in $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$7 (Nowak et al., 2013).

The concrete Bose-gas simulations that established the link to superfluid turbulence observed precisely these values: in $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$8, an infrared law $$n(\mathbf{k}, t) = \langle \psi^\dagger(\mathbf{k}, t)\psi(\mathbf{k}, t)\rangle ,$$9 and in $n(k,t)$0, $n(k,t)$1, together with a thermalized high-$n(k,t)$2 tail $n(k,t)$3 (Nowak et al., 2010). The ultraviolet sector can additionally display a weak-wave-turbulence scaling with $n(k,t)$4 in suitable parameter regimes (Nowak et al., 2010).

The Kolmogorov analogy becomes explicit upon decomposing the kinetic energy into incompressible, compressible, and quantum-pressure parts. With $n(k,t)$5 and superfluid velocity $n(k,t)$6, one defines generalized velocities

$n(k,t)$7

and decomposes $n(k,t)$8 into incompressible and compressible components, $n(k,t)$9 and $k=|\mathbf{k}|$0 (Nowak et al., 2010). The associated radial energy spectra are related to occupation-like spectra through

$k=|\mathbf{k}|$1

In the turbulent regime, the incompressible component develops a kinetic-energy spectrum consistent with

$k=|\mathbf{k}|$2

which is the direct quantum analogue of Kolmogorov’s inertial-range law (Nowak et al., 2010).

This structural relation between fixed-point scaling and cascade transport is also borne out experimentally. In a homogeneous three-dimensional Bose gas, the ultraviolet scaling exponents $k=|\mathbf{k}|$3 and $k=|\mathbf{k}|$4 satisfy $k=|\mathbf{k}|$5, in agreement with the energy-conservation constraint $k=|\mathbf{k}|$6 for $k=|\mathbf{k}|$7 and $k=|\mathbf{k}|$8, while the measured $k=|\mathbf{k}|$9 is consistent with the weak-wave-turbulence prediction $n(sk) = s^{-\zeta} n(k),$0 (Madeira et al., 2022). This identifies a direct energy cascade toward high momenta, coexisting with an infrared particle cascade.

4. Vortices, defect statistics, and defect-controlled scaling

A defining advance in the Bose-gas literature was the identification of the infrared non-thermal fixed point with a gas of quantized vortices. In the Madelung representation, vortices correspond to phase singularities around points in two dimensions or lines and loops in three dimensions where the density vanishes (Nowak et al., 2010). Their velocity profile scales as $n(sk) = s^{-\zeta} n(k),$1, and after regularization by the factor $n(sk) = s^{-\zeta} n(k),$2 the resulting generalized velocity $n(sk) = s^{-\zeta} n(k),$3 remains finite (Nowak et al., 2010).

For a finite density of independent vortices and antivortices in $n(sk) = s^{-\zeta} n(k),$4, the superposition of long-range vortex velocity fields yields

$n(sk) = s^{-\zeta} n(k),$5

while in $n(sk) = s^{-\zeta} n(k),$6 independent vortex lines and rings yield

$n(sk) = s^{-\zeta} n(k),$7

(Nowak et al., 2011). These exponents coincide with the strong-turbulence prediction $n(sk) = s^{-\zeta} n(k),$8 and thereby provide a direct physical interpretation of the infrared fixed point in terms of defect statistics (Nowak et al., 2011).

The same paper demonstrates that pair correlations alter the scaling. In two dimensions, a vortex–antivortex pair behaves as a dipole whose far-field velocity decays faster, leading to

$n(sk) = s^{-\zeta} n(k),$9

at sufficiently small $n(k)\sim k^{-\zeta}$0 (Nowak et al., 2011). In three dimensions, small rings or strongly correlated line segments produce analogous weakening of the infrared exponent, such as $n(k)\sim k^{-\zeta}$1 in appropriate regimes (Nowak et al., 2011). The long-time departure from the non-thermal fixed point is therefore associated with vortex–antivortex pairing and annihilation rather than with a mere loss of scaling (Nowak et al., 2011).

This defect-based interpretation was generalized in later work. In a two-component Bose gas, universal scaling exponents differ across the miscible–immiscible transition because the dominant defects change from vortices and skyrmions to domain walls and spin textures (Karl et al., 2013). In that setting, vortices generate incompressible spectra $n(k)\sim k^{-\zeta}$2, while domain walls yield spin-pressure spectra $n(k)\sim k^{-\zeta}$3, and the total single-particle spectrum reflects the superposition of these defect sectors (Karl et al., 2013). More generally, the geometry of the dominant defects determines the scaling law: vortices in $n(k)\sim k^{-\zeta}$4 dimensions produce $n(k)\sim k^{-\zeta}$5, while solitons or domains produce $n(k)\sim k^{-\zeta}$6 (Nowak et al., 2013).

5. Acoustic sectors, multibranch fixed points, and multiple attractors

Kolmogorov-type non-thermal fixed points are not exhausted by incompressible vortex turbulence. The compressible sector can itself exhibit a non-thermal power law. In the ultracold Bose-gas simulations, the compressible contribution follows

$n(k)\sim k^{-\zeta}$7

which is interpreted as acoustic wave turbulence driven by long-wavelength sound modes on top of the vortical background (Nowak et al., 2010). These acoustic excitations can survive even after vortices have annihilated, indicating a distinct but related universality sector of the non-thermal fixed-point dynamics (Nowak et al., 2010).

The possibility of multiple non-thermal fixed points for the same microscopic symmetry was demonstrated experimentally in a quasi-one-dimensional spinor Bose gas. Two distinct initial conditions in the same easy-plane ferromagnetic phase yield two different scaling regimes with constrained exponents

$n(k)\sim k^{-\zeta}$8

and

$n(k)\sim k^{-\zeta}$9

together with different spectral tail exponents and different dominant excitations, namely phase modes in one case and phase-amplitude defects in the other (Lannig et al., 2023). This establishes that symmetry alone does not uniquely determine far-from-equilibrium universality; basins of attraction and defect content matter as well (Lannig et al., 2023). A plausible implication is that “Kolmogorov-type” should be understood as a family of scaling attractors rather than a single fixed point.

A related, more abstract generalization appears in a 2026 multibranch shell model of turbulent cascades, where a continuum of stationary forward-cascade fixed points exists, including the Kolmogorov solution $$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$0 with $$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$1 (Tuteri et al., 8 Jan 2026). In that model, the Kolmogorov fixed point is one member of a manifold of non-thermal cascade states distinguished by multiplier statistics and intermittency exponents (Tuteri et al., 8 Jan 2026). While this shell-model analysis is not a Bose-gas result, it reinforces the broader interpretation of Kolmogorov-type non-thermal fixed points as stationary cascade solutions, often organized into a family or manifold rather than an isolated point (Tuteri et al., 8 Jan 2026).

6. Dynamics of approach, stability, and experimental realizations

The dynamical route to a non-thermal fixed point typically begins with a far-from-equilibrium state involving strong overpopulation of low-momentum modes or violent phase perturbations. In dilute Bose gases governed semiclassically by the Gross–Pitaevskii equation,

$$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$2

such initial conditions produce shock-wave-like fronts, defect nucleation, and the emergence of vortex tangles or soliton trains (Nowak et al., 2010). The system then enters a quasi-stationary regime with bimodal scaling, critical slowing down, and long-lived power laws before ultimately thermalizing (Nowak et al., 2013).

Experimental evidence for this picture is now extensive. A quasi-one-dimensional Bose gas of $$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$3Rb atoms showed infrared scaling exponents

$$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$4

with $$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$5 for $$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$6, consistent with particle-conserving infrared scaling (Madeira et al., 2022). A homogeneous three-dimensional Bose gas of $$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$7K exhibited simultaneous infrared and ultraviolet scaling,

$$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$8

providing a clean experimental realization of bidirectional cascades and a Kolmogorov-type dual-fixed-point structure (Madeira et al., 2022). In turbulent harmonically trapped three-dimensional condensates, exponents extracted from the reconstructed three-dimensional distributions were

$$\bar n(k,t)=\left(\frac{t}{t_0}\right)^\alpha f\!\left[\left(\frac{t}{t_0}\right)^\beta k\right],$$9

indicating ultraviolet-dominated transport associated with quantum turbulence (Madeira et al., 2022).

More recently, a quasi-one-dimensional Bose gas of strongly interacting $n(k,t)=t^\alpha f(t^\beta k),$0Li$n(k,t)=t^\alpha f(t^\beta k),$1 Feshbach molecules, driven far from equilibrium by imprinting a white-noise phase profile, was found to relax with the same exponents as the earlier weakly interacting $n(k,t)=t^\alpha f(t^\beta k),$2Rb system, namely

$n(k,t)=t^\alpha f(t^\beta k),$3

in a common quasi-one-dimensional scaling window (Liang et al., 26 May 2025). This points to a single universal fixed point with a large basin of attraction governing the relaxation of quasi-one-dimensional bosonic systems, largely independent of microscopic details (Liang et al., 26 May 2025).

Stability analysis in other nonequilibrium quantum field theories supports the same attractor interpretation. In longitudinally expanding kinetic theory for a non-Abelian plasma, linearization about the BMSS non-thermal fixed point yields negative eigenvalues

$n(k,t)=t^\alpha f(t^\beta k),$4

demonstrating linear stability and algebraic relaxation toward the fixed point (Mikheev et al., 2022). Although this result concerns an expanding plasma rather than a Bose gas, it supplies a technical example of how non-thermal fixed points function as dynamical attractors with calculable relaxation rates (Mikheev et al., 2022).

7. Broader extensions, anomalous cases, and open problems

The concept of Kolmogorov-type non-thermal fixed points now spans nonrelativistic Bose gases, relativistic scalar field theories, non-Abelian plasmas, holographic superfluids, shell models, and even dilute Fermi gases (Berges et al., 2016). In a relativistic $n(k,t)=t^\alpha f(t^\beta k),$5 scalar theory solved beyond the weak-coupling classical-statistical limit, non-thermal fixed points persist up to couplings of order unity, with an infrared inverse particle cascade

$n(k,t)=t^\alpha f(t^\beta k),$6

and, in the symmetric case, a direct energy cascade

$n(k,t)=t^\alpha f(t^\beta k),$7

demonstrating that Kolmogorov-type scaling survives beyond weak coupling (Berges et al., 2016).

A holographic superfluid in $n(k,t)=t^\alpha f(t^\beta k),$8 dimensions supplies a strongly coupled realization. Starting from various vortex ensembles, the system exhibits intermediate-time Kolmogorov scaling whose emergence depends on initial conditions, followed by a universal late-time regime where the occupation spectrum scales as $n(k,t)=t^\alpha f(t^\beta k),$9 with $\alpha$0 and characteristic length scales grow as $\alpha$1 (Ewerz et al., 2014). The late-time universal regime is interpreted as a non-thermal fixed point of the dynamical evolution, and, in the gravity dual, as a stationary point of the classical bulk equations (Ewerz et al., 2014).

Recent work on decaying two-dimensional superfluid turbulence near an anomalous non-thermal fixed point sharpens the distinction between spatial and temporal universality. During a universal interval with inter-defect distance

$\alpha$2

moments of the superfluid velocity circulation obey

$\alpha$3

with intermittency corrections for higher $\alpha$4 consistent with values measured in fully developed classical turbulence (Rasch et al., 25 Sep 2025). The spatial structure is therefore Kraichnan–Kolmogorov-like, while the temporal coarsening exponent is anomalously slow and distinctly quantum (Rasch et al., 25 Sep 2025). This suggests that Kolmogorov-type non-thermal fixed points need not share classical coarsening exponents even when their spatial turbulence statistics do.

Several open issues recur across the literature. A complete classification of non-thermal universality classes remains unavailable (Madeira et al., 2022). The role of dimensionality, internal symmetry, and defect content remains only partially understood, especially in multicomponent and low-dimensional systems (Lannig et al., 2023). The crossover from the vicinity of a non-thermal fixed point to final thermalization is also an active problem, as is the precise relation between defect-mediated coarsening, wave turbulence, and hydrodynamic attractors (Nowak et al., 2013). Inference beyond the cited results suggests that the encyclopedia entry “Kolmogorov-type non-thermal fixed point” properly denotes a unifying concept rather than a single model-specific phenomenon: stationary or quasi-stationary cascade states, often defect-mediated, characterized by universal spatio-temporal scaling and organized by conservation laws rather than by thermal detailed balance.