Self-Consistent Relaxation Theory Overview

• Self-Consistent Relaxation Theory is a framework that simultaneously determines the relaxing quantity and its mediating processes using exact evolution equations and moment constraints.
• It applies to diverse fields—from classical many-body systems to quantum kinetics and gravitational dynamics—by employing memory functions, recurrence relations, and collision dressing for self-consistency.
• The approach enables parameter-free predictions of dynamic structure factors, transport coefficients, and spectral properties by closing hierarchical equations without imposing external relaxation laws.

“Self-Consistent Relaxation Theory” is used in several technically distinct literatures, but the phrase consistently denotes a class of formulations in which relaxation is determined together with the quantities that mediate it—correlation functions, response kernels, collision terms, collective fields, or dressed states—rather than imposed by an external ansatz. In classical many-body theory it is a memory-function and recurrence-relation framework for time-correlation functions; in fluids and plasmas it generates $S(k,\omega)$, transport coefficients, and collective-mode spectra from equilibrium structure and exact moments; in kinetic theory it yields collision operators dressed by dielectric response and constrained by conservation laws; in gravitational dynamics it appears in collisionless and resonant collective relaxation; and in condensed-matter and electronic-structure settings it refers to self-consistent treatments of built-in stress, orbital relaxation, and environment-dressed dissipation (Mokshin, 2019, Mokshin et al., 2022, Bar-Or et al., 2020, Bevzenko et al., 2014). This suggests a methodological family rather than a single field-independent formalism.

1. Common structural features

Across its uses, the term denotes frameworks that begin from exact evolution equations or constrained variational principles and then close the problem by imposing self-consistency. In the classical correlation-function literature, self-consistency is encoded in recurrence coefficients $\Delta_n$ fixed by exact moments and sum rules. In kinetic theories, it appears through dielectric dressing, Landau matching, or collisional invariants. In gravitational and MHD problems, it is tied to self-generated potentials or invariant-constrained equilibria. In electronic-structure variants, it refers to orbital optimization or curvature cancellation performed together with the relaxation process itself (Mokshin, 2019, Wu et al., 2024, Hamilton et al., 2018).

Domain	Central objects	Self-consistency device
Classical multiparticle systems	$C_{AB}(t)$, $\phi(t)$, $\Delta_n$	moment constraints, continued fractions
Liquids and plasmas	$S(k,\omega)$, $C_T(k,\omega)$, $\Delta_n(k)$	structural input, exact sum rules, high-level closures
Quantum kinetic theory	Wigner components, collision terms, $\epsilon(k,\omega)$	collisional invariants, Landau matching, dielectric response
Gravitational dynamics	$f(\mathbf{x},\mathbf{v},t)$, $\Delta_n$0, response matrix	Vlasov–Poisson or Balescu–Lenard self-gravity
Elastic and electronic structure	$\Delta_n$1, $\Delta_n$2, orbital occupations, dressed jump operators	cavity closure, orbital relaxation, dressed-basis rates

A recurrent feature is the rejection of purely phenomenological relaxation laws. The different formalisms instead enforce exact short-time expansions, spectral positivity, hydrodynamic limits, conservation laws, or stationary-state conditions. In that restricted sense, “self-consistent relaxation” always means that the relaxing object and the effective medium through which it relaxes are solved together.

2. Memory functions, recurrence relations, and exact constraints

Mokshin’s formulation for classical multiparticle systems places time-correlation functions at the center of the theory (Mokshin, 2019). For observables $\Delta_n$3 and $\Delta_n$4,

$\Delta_n$5

and the normalized autocorrelation function is

$\Delta_n$6

Its short-time expansion,

$\Delta_n$7

already shows the role of even frequency moments $\Delta_n$8.

The Zwanzig–Mori projection formalism yields an exact hierarchy. For a chosen slow variable $\Delta_n$9, the normalized correlator $C_{AB}(t)$0 obeys

$C_{AB}(t)$1

where $C_{AB}(t)$2 is the first memory function and $C_{AB}(t)$3. The associated Krylov/Mori chain,

$C_{AB}(t)$4

produces the continued-fraction representation

$C_{AB}(t)$5

The theory becomes predictive because the $C_{AB}(t)$6 are not adjustable. They are reconstructed from moment ratios, for example

$C_{AB}(t)$7

with analogous higher-order expressions. This construction guarantees the exact short-time expansion, positivity of spectra, and the relevant sum rules.

Several closures are then possible. A finite variable set gives undamped oscillatory solutions. Equal-time-scale closures such as $C_{AB}(t)$8 yield analytic Bessel-type memory functions. Arithmetic progression $C_{AB}(t)$9 gives Gaussian relaxation,

$\phi(t)$0

while $\phi(t)$1 for all $\phi(t)$2 produces damped oscillatory Bessel relaxation. The same framework also generalizes mode-coupling approximations and includes the Yulmetyev–Shurygin correlation approximations as special cases. For supercooled liquids, a scaling closure of the form

$\phi(t)$3

reproduces two-step relaxation, a nonergodicity plateau, and MCT-like $\phi(t)$4- and $\phi(t)$5-regime asymptotics.

3. Collective dynamics in liquids and plasmas

In strongly coupled Yukawa and Coulomb one-component plasmas, self-consistent relaxation theory is an explicitly parameter-free route from equilibrium structure to the dynamic structure factor, mode dispersion, and attenuation (Mokshin et al., 2022, Fairushin et al., 2023, Mokshin et al., 2021). The basic variables are the density correlator $\phi(t)$6, the dynamic structure factor $\phi(t)$7, and a hierarchy of frequency relaxation parameters $\phi(t)$8 defined by exact moment ratios.

For the Yukawa OCP, the first two relaxation parameters are

$\phi(t)$9

and

$\Delta_n$0

where $\Delta_n$1 is an explicit $\Delta_n$2-dependent integral. In the intermediate screening regime explored in the paper,

$\Delta_n$3

with $\Delta_n$4. This closes the hierarchy and yields a closed-form $\Delta_n$5 with a bicubic denominator, Brillouin peaks, sound speed, and sound attenuation, all without adjustable parameters. The theory is validated for $\Delta_n$6 and $\Delta_n$7, and reproduces the roton minimum near $\Delta_n$8 for $\Delta_n$9, $S(k,\omega)$0.

For the Coulomb OCP, the same architecture is adapted to the long-range case. The exact two-particle-level inputs are

$S(k,\omega)$1

and

$S(k,\omega)$2

with $S(k,\omega)$3. For fluid states $S(k,\omega)$4, the higher-order correlations are encoded empirically by

$S(k,\omega)$5

where $S(k,\omega)$6 and $S(k,\omega)$7. This reproduces plasma side peaks near $S(k,\omega)$8, their damping, the high-$S(k,\omega)$9 central peak, and the positive-to-negative dispersion crossover, with a low-$C_T(k,\omega)$0 estimate $C_T(k,\omega)$1.

A closely related transverse-current formulation describes equilibrium liquids. There the central object is

$C_T(k,\omega)$2

with a continued-fraction hierarchy in $C_T(k,\omega)$3. The theory recovers the hydrodynamic Lorentzian at small $C_T(k,\omega)$4, the Gaussian free-particle spectrum at large $C_T(k,\omega)$5, and the onset of propagating shear excitations when

$C_T(k,\omega)$6

The generalized kinematic viscosity is

$C_T(k,\omega)$7

and in liquid lithium near melting the theory yields $C_T(k,\omega)$8 and $C_T(k,\omega)$9. Taken together, these plasma and liquid applications show that the memory-function variant of self-consistent relaxation theory is a structure-based spectral theory spanning hydrodynamic, microscopic, and free-particle regimes.

4. Quantum-kinetic and dielectric formulations

A second major usage of the term is kinetic rather than correlational: relaxation is described by a self-consistent collision operator dressed by collective response (Bar-Or et al., 2020, Crowley, 2015, Wu et al., 2024). In fuzzy dark matter halos, the starting point is the Schrödinger–Poisson system and the Wigner distribution. Interference on the de Broglie scale produces persistent density granules, and the resulting stochastic gravitational potential leads to a Lenard–Balescu-type kinetic equation with bosonic enhancement: $\Delta_n(k)$0 In the full BL form the kernel contains $\Delta_n(k)$1, where $\Delta_n(k)$2 is the dielectric function of the halo. In the Landau limit the same dynamics reduces to a Fokker–Planck equation with effective mass $\Delta_n(k)$3, and the relaxation time scales as $\Delta_n(k)$4. The stationary solutions are Bose–Einstein distributions in velocity space, and for sufficiently low $\Delta_n(k)$5 the thermal component saturates at $\Delta_n(k)$6 and the remaining mass condenses into a soliton core.

In multicomponent Coulomb plasmas, the response is constructed within RPA supplemented by static local field corrections. The fundamental response equation is

$\Delta_n(k)$7

and the effective interaction is written as $\Delta_n(k)$8. The general temperature-relaxation formula derived in this framework is explicitly stated to be identical to the Daligault–Dimonte result for the two-component case, while the present approach differs in how the LFCs are obtained: they are reconstructed from static pair correlations, including quantal electrons and possible bound states. This gives a self-consistent route from $\Delta_n(k)$9, $\epsilon(k,\omega)$0, and $\epsilon(k,\omega)$1 to equilibration coefficients and low-frequency transport.

In relativistic quantum kinetic theory, the self-consistent element shifts from the dielectric kernel to the collision sector itself. Using the gauge-invariant Wigner function, collision terms are constrained so that they preserve charge and energy-momentum conservation. At the level of collisional invariants this is expressed as

$\epsilon(k,\omega)$2

with analogous axial constraints. The paper then constructs a relaxation-time approximation in terms of projected nonconserved deviations $\epsilon(k,\omega)$3 and $\epsilon(k,\omega)$4, with auxiliary vector functions $\epsilon(k,\omega)$5 and $\epsilon(k,\omega)$6, so that the RTA automatically satisfies the zeroth- and first-order Wigner constraints. Here self-consistency means that the collision term is not guessed independently of the Wigner hierarchy; it is derived so that lower-order on-shell and conservation conditions remain intact.

5. Gravitational and magnetohydrodynamic realizations

In self-gravitating dynamics, the phrase refers both to collisionless mixing and to collectively dressed secular transport (Labini, 2012, Hamilton et al., 2018, Servidio et al., 2014). For the collapse of an isolated uniform spherical cloud, the exact continuum dynamics is Vlasov–Poisson,

$\epsilon(k,\omega)$7

and the initial virial ratio $\epsilon(k,\omega)$8 controls the relaxation regime. Simulations identify a critical value $\epsilon(k,\omega)$9. For $f(\mathbf{x},\mathbf{v},t)$0, the evolution is a mild relaxation: the cloud roughly preserves its size, no particles are ejected, $f(\mathbf{x},\mathbf{v},t)$1, and the quasi-stationary state is well described by Lynden-Bell theory in a box, with a density profile

$f(\mathbf{x},\mathbf{v},t)$2

For $f(\mathbf{x},\mathbf{v},t)$3, the evolution is violent: the system contracts to $f(\mathbf{x},\mathbf{v},t)$4, ejects mass and energy, and develops a halo

$f(\mathbf{x},\mathbf{v},t)$5

with $f(\mathbf{x},\mathbf{v},t)$6 at $f(\mathbf{x},\mathbf{v},t)$7. In this regime LB-in-a-box fails because the bound component is effectively open.

For globular clusters, the relevant self-consistent relaxation theory is the Balescu–Lenard equation in angle–action space. It replaces the classical Coulomb-logarithm picture with a resonant flux

$f(\mathbf{x},\mathbf{v},t)$8

whose kernel is dressed by the response matrix $f(\mathbf{x},\mathbf{v},t)$9. The paper argues heuristically that relaxation is not predominantly two-particle scattering and is enhanced by self-gravity; low-order collective modes, especially $\Delta_n$00 and $\Delta_n$01, can dominate the secular flux. A complete theory is suggested to require a decomposition into a finite sum over small wavenumbers plus a large-wavenumber integral analogous to classical theory.

In incompressible MHD turbulence, self-consistent relaxation refers to equilibria that are simultaneously stationary solutions of the MHD equations and extrema of constrained energy minimization. Minimization of magnetic energy at fixed helicity gives the Taylor state,

$\Delta_n$02

while combined constraints on $\Delta_n$03, $\Delta_n$04, and $\Delta_n$05 yield the Beltrami family $\Delta_n$06. Cluster multi-spacecraft observations show that such relaxation does not occur only as a late-time global state: local patches in the solar wind display strong PDF peaks in $\Delta_n$07 near $\Delta_n$08, and corresponding Alfvénic alignments in $\Delta_n$09. This suggests that local suppression of nonlinear terms can occur within the cascade itself.

6. Elastic, orbital, and excited-state variants

In glassy liquids and aperiodic solids, self-consistent relaxation theory is formulated as an Onsager-like cavity theory for built-in stress (Bevzenko et al., 2014). The effective elastic tensor satisfies

$\Delta_n$10

where $\Delta_n$11 is the static susceptibility associated with the response of the built-in stress to an imposed average strain. For spherical inclusions in an isotropic medium, the Eshelby tensor yields three fixed points in the renormalization flow: the uniform-liquid point, the infinitely compressible solid, and a nontrivial point at Poisson ratio

$\Delta_n$12

At $\Delta_n$13, $\Delta_n$14, so hydrostatic and deviatoric sectors renormalize identically. The theory predicts a discontinuous jump in the finite-frequency shear modulus at the crossover from collisional to activated transport, consistent with RFOT.

In excited-state electronic structure, the phrase refers to state-specific orbital relaxation followed by post-HF correlation (Ye et al., 2020). A target excited configuration is optimized directly at the SCF level—typically ROHF, HPHF, or SPHF—and then corrected by perturbation theory. For ROHF references, semicanonical orbitals are defined within the core, open-shell, and virtual subspaces, and the correlated excitation energy is obtained from

$\Delta_n$15

The paper finds that ROMP2 gives excitation energies with a mean unsigned error of about $\Delta_n$16 over 104 low-lying singlet vertical excitations at non-iterative $\Delta_n$17 cost. Here “relaxation” means explicit orbital adaptation to the excited-state occupation pattern, rather than linear response around the ground-state reference.

A related but distinct DFT+$\Delta_n$18 usage treats orbital self-energy and wavefunction relaxation as separate eigenvalue contributions (Huang, 2015): $\Delta_n$19 Using Janak’s theorem,

$\Delta_n$20

the self-consistent goal is to choose $\Delta_n$21 so that the curvature is cancelled and the exact piecewise-linear condition is approximated. The paper emphasizes that fully occupied shells admit a unique set of occupations for which the self-energy and relaxation residues offset, whereas partially occupied shells retain a nonzero residue. In that sense, self-consistent relaxation is a curvature-cancellation strategy performed together with the orbital optimization.

7. Nonequilibrium quantum dynamics and environment-dressed dissipation

In nonequilibrium disordered fermion systems, Weidinger, Gopalakrishnan, and Knap use a self-consistent Hartree–Fock approach to study far-from-equilibrium relaxation and the onset of many-body localization (Weidinger et al., 2018). The time-local HF evolution obeys

$\Delta_n$22

with

$\Delta_n$23

At weak disorder, the self-energy fluctuates strongly and can be interpreted as a self-consistent noise process; the local spectral function is broad, and in random systems the imbalance decays subdiffusively as $\Delta_n$24 with $\Delta_n$25. At strong disorder, the self-energy exhibits only a few coherent oscillations, the local spectral function resolves into sharp spikes, and memory of the initial state persists. In quasi-periodic potentials, the subdiffusive response ceases to exist because rare-region effects are absent.

In ultra- and deep-strong coupling, self-consistent relaxation theory means that dissipation must be computed in the dressed basis of a Hamiltonian that already includes counter-rotating and diamagnetic terms (Sergeev et al., 2023). For the bosonic Hopfield-type model, the system is diagonalized first, yielding dressed frequencies $\Delta_n$26, and only then are environment-induced decay rates evaluated: $\Delta_n$27 This avoids the failure of bare-mode Lindblad or RWA dissipators in the USC/DSC regime. For frequency-independent density of states of the environment, the relaxation rates decrease exponentially with increasing coupling,

$\Delta_n$28

whereas if the microscopic environment coupling grows with frequency, the rates can instead increase with $\Delta_n$29. In this usage, self-consistency lies in computing relaxation from the same dressed spectrum and operators that define the strongly coupled system.

Taken together, these literatures show that “Self-Consistent Relaxation Theory” is best understood as a recurrent strategy: exact or near-exact dynamical equations are retained as far as possible, and the closure is imposed on response functions, moments, collision kernels, invariant constraints, or dressed states rather than on relaxation laws themselves. The resulting theories differ sharply in ontology—correlators, Wigner functions, action-space fluxes, elastic susceptibilities, orbital occupations, or dressed modes—but they share the same organizing principle that relaxation and the medium of relaxation must be solved simultaneously.