Hot-Electron Limit (HEL) Closure

• Hot-Electron Limit (HEL) closure is a modeling principle that simplifies non-equilibrium electron kinetics by assuming rapid electron thermalization relative to slower subsystems like phonons or ions.
• It underpins various formulations—from two-temperature models in metals to electron-ion coupling and gyrokinetic moment closures—bridging microscopic kinetics to tractable macroscopic equations.
• HEL closure ensures analytical consistency via conservation laws, response functions, and sum rules, while exposing limitations in capturing non-thermal effects and nonlinear dynamics.

Hot-Electron Limit (HEL) closure denotes a class of reduced descriptions for non-equilibrium systems in which the electron subsystem is treated as internally equilibrated, or asymptotically simplified, while remaining out of equilibrium with phonons, ions, or other energy reservoirs. In ultrafast metals this usually means a transient regime with $T_e \gg T$ and a hot Fermi-Dirac electron distribution coupled to a colder Bose-Einstein phonon bath; in two-temperature electron-ion theories it means evolution equations for $T_e$ and $T_i$ closed by an electron-ion coupling factor; in gyrokinetics it means an asymptotic closure of the moment hierarchy in the limit $\tau=T_i/T_e\ll 1$; and in molecular plasmas it can take the form of a closed superelastic gain relation. Across these settings, HEL closure is the step that replaces full kinetic dynamics by a finite-dimensional transport, response, or moment system while retaining the relevant statistics, response functions, and limiting behaviors (Singh, 2024, Hoffmann et al., 18 Sep 2025, Parent, 28 Jan 2026).

1. Definition and constitutive assumptions

In the metallic two-temperature setting, the HEL is a transient state created after ultrafast energy input, typically by a femtosecond laser, in which electrons heat rapidly and their temperature becomes substantially greater than that of the lattice. The closure assumes that electron-electron interactions are much faster than electron-phonon scattering, so electrons quickly reach a quasi-equilibrium Fermi-Dirac distribution at temperature $T_e$, while phonons remain described by a Bose-Einstein distribution at temperature $T$ (Singh, 2024). In this form, the reduced dynamics are written as

$$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$

In the broader electron-ion formulation, the same logic appears as a two-temperature decomposition into electronic and ionic subsystems, with Maxwellian ion velocities at $T_i$ and quantum electrons at $T_e$. The temperature relaxation is then written as

$$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$

where $T_e$0 is the electron-ion coupling factor. The underlying assumptions stated for this closure are homogeneity, no external sources, classical ions, quantum electrons, adiabatic separation between electron and ion timescales, and internal equilibration of each subsystem so that the distribution functions are characterized solely by $T_e$1 and $T_e$2 (Daligault et al., 2019).

A recurrent point in the literature is that HEL closure is not a single formula but a closure principle. Depending on the problem, it may identify a heat-transfer coefficient $T_e$3, a coupling factor $T_e$4, a memory function $T_e$5, a nonlocal heat-flux closure, a finite set of gyrokinetic moments, or an analytic gain function for superelastic heating.

2. Two-temperature metals and memory-function closure

A direct microscopic realization of HEL closure in metals is the generalization of the Götze-Wölfle memory function formalism to non-equilibrium hot-electron systems, where the electron subsystem is at temperature $T_e$6 and the phonon bath at temperature $T_e$7. The key step is to treat the Fermi-Dirac distributions for electrons and the Bose-Einstein distributions for phonons independently, as in the two-temperature model. The imaginary part of the generalized memory function is written as

$T_e$8

with $T_e$9 the Bose distribution at lattice temperature $T_i$0, $T_i$1 the Bose distribution at electron temperature $T_i$2, $T_i$3 the Debye momentum cutoff, and $T_i$4. For $T_i$5, this reduces to the standard Götze-Wölfle result (Das et al., 2015).

The resulting asymptotics distinguish sharply between regimes where both subsystem temperatures matter and regimes where the scattering rate is controlled only by the lattice.

Regime	$T_i$6 behavior
$T_i$7	$T_i$8
$T_i$9	$\tau=T_i/T_e\ll 1$0 with no $\tau=T_i/T_e\ll 1$1 dependence
$\tau=T_i/T_e\ll 1$2	const $\tau=T_i/T_e\ll 1$3 with no $\tau=T_i/T_e\ll 1$4 dependence
$\tau=T_i/T_e\ll 1$5	const $\tau=T_i/T_e\ll 1$6 with no $\tau=T_i/T_e\ll 1$7 dependence

For dc transport, the low-temperature result

$\tau=T_i/T_e\ll 1$8

recovers Bloch $\tau=T_i/T_e\ll 1$9 behavior at equilibrium and gives additive $T_e$0 and $T_e$1 contributions away from equilibrium. At high temperatures,

$T_e$2

with $T_e$3 and $T_e$4, so the dc scattering rate is independent of the electronic temperature $T_e$5 and dictated by the phonon temperature $T_e$6 only. For optical scattering at $T_e$7, the high-frequency rate is likewise independent of the electron temperature and depends only on the lattice temperature (Das et al., 2015).

These results are used in the source paper as closure information for hot-electron kinetics. The analytical expressions provide “the explicit form for the electron-phonon coupling kernel, valid away from equilibrium,” enabling rigorous and consistent closure of coupled two-temperature equations such as

$T_e$8

with $T_e$9 related to $T$0 (Das et al., 2015).

3. Electron-ion coupling factor as a general HEL closure

A more general HEL closure is provided by the electron-ion coupling factor theory spanning hot solid metals, liquid metals, warm dense matter, and plasmas. In this formulation, the rate of energy exchange is tied to an effective friction between ions and electrons. The basic relation is

$T$1

where the average friction coefficient is defined from a Kubo-type force autocorrelation. For a homogeneous system, the closure may be written in response-function form as

$T$2

This representation explicitly embeds the electron-ion interaction potential and the full quantum many-body electronic density response function (Daligault et al., 2019).

The theory is stated to self-consistently incorporate quantum mechanics and statistics for electrons, strong correlations and screening through $T$3, thermal and disorder effects, and particle correlations. It is also described as especially suitable for Averaged Atom Models, which offer an effective route to include non-ideal interaction effects at much reduced computational cost relative to first-principles quantum molecular dynamics simulations (Daligault et al., 2019).

A defining feature of this framework is that it unifies limiting cases that had previously been treated separately. In the weak-coupling hot plasma limit it reduces to the Spitzer result; in the limit of hot solids with lattice and electronic temperatures much greater than the Debye temperature, and with ionic motion harmonicized into phonons, it reduces to the Allen electron-phonon coupling formula

$T$4

For intermediate regimes such as liquid metals and warm dense matter, the full Kubo/response formulation is retained (Daligault et al., 2019).

The same work frames HEL closure as a regularization problem. Previous models such as Fermi’s golden rule can break down in the hot-electron limit and yield divergent results if screening and many-body correlations are not correctly included. The response-based closure remains finite and well-behaved for $T$5, $T$6, and $T$7 because it includes full screening, Fermi-Dirac statistics, and finite-temperature response functions (Daligault et al., 2019).

4. Conserving closures in ultrafast metals and dense hydrogen

The historical two-temperature closure in metals is now understood to be conditional. The review of ultrafast non-equilibrium electron relaxation states that the ideal HEL and TTM closure rests on $T$8, but pump-probe and photoemission experiments showed non-thermal electron distributions persisting for hundreds of femtoseconds, comparable to $T$9. In those cases a two-subsystem picture is insufficient and a three-subsystem description is needed, separating nonthermal electrons, thermalized electrons, and phonons. The same review also states that full Boltzmann calculations confirmed that the standard TTM/HEL closure is invalid unless the initial excitation is very large, whereas at high fluence the TTM closure is approximately valid (Singh, 2024).

For dense two-temperature hydrogen, the closure question is formulated in terms of energy relaxation rates. The work examines the Fermi Golden Rule in its $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$0-sum form,

$$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$1

and contrasts it with coupled-mode formulations in which electron and ion collective excitations enter through a coupled denominator:

$$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$2

The paper advocates the $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$3-sum version of the Fermi Golden Rule formula as the most convenient method for calculating the rate of cooling of hot electrons where energy is transferred to cold ions (Dharma-wardana, 2017).

A central issue is the interaction potential. The same source states that coupled-mode calculations using the simple Coulomb potential $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$4 within RPA will greatly over-estimate the coupled-mode contribution, whereas a weak pseudopotential $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$5 would probably bring the estimated coupled-mode contribution into agreement with simulations. It further argues that reduced models often neglect important analytical constraints such as the $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$6-sum rule, higher-moment sum rules, and Kramers-Krönig relations. In that account, satisfaction of sum rules and Kramers-Krönig relations is the litmus test for acceptability of a physical closure (Dharma-wardana, 2017).

Taken together, these works distinguish between a merely reduced closure and a conserving closure. A conserving HEL closure is one that remains compatible with spectral sum rules, correct static structure, and properly screened interaction physics.

5. HEL closure in nonlocal heat transport and gyrokinetic hierarchies

In nonlocal electron thermal conduction, HEL closure appears as the assumption that the hot-electron population adjusts instantly to the macroscopic plasma profiles. The comparison of nonlocal transport models against Vlasov-Fokker-Planck calculations tests three such closures: the eigenvector integral closure (EIC), the non-Fourier Landau-fluid (NFLF) model, and the Schurtz-Nicolaï-Busquet multigroup diffusion model (SNB). In the small-amplitude temperature-perturbation problem, EIC and NFLF predict the damping rate within $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$7 at moderate collisionalities. In the large-temperature-difference problem, however, EIC and NFLF overestimate the peak heat flow by as much as $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$8 and do not predict preheat. SNB agrees better with VFP results for that problem if care is taken with the definition of the mean free path. In a hohlraum-relevant case with inhomogeneous ionisation, SNB overestimates the heat flow in the helium gas-fill by a factor of $$\gamma T_e \frac{dT_e}{dt} = -\alpha (T_e - T), \qquad C_{ph} \frac{dT}{dt} = +\alpha (T_e - T).$$9 despite predicting the peak heat flux to within $T_i$0 (Brodrick et al., 2017).

In gyrokinetics, HEL closure has a different but closely related meaning: an asymptotic truncation of the gyromoment hierarchy in the small temperature-ratio limit $T_i$1. By expanding gyroaveraging kernels and retaining only essential $T_i$2 terms, the hierarchy is closed on four moments: density $T_i$3, parallel velocity $T_i$4, parallel temperature $T_i$5, and perpendicular temperature $T_i$6. The resulting quasi-neutrality relation is

$T_i$7

In Z-pinch geometry the closed gyromoment system is analytically equivalent to the model of Ivanov et al.; numerical benchmarks reproduce established linear growth rates, nonlinear heat transport, and low-collisionality dynamics (Hoffmann et al., 18 Sep 2025).

The source also states the limits of this truncation. In tokamak-relevant $T_i$8-$T_i$9 geometry, transport levels and temporal dynamics are qualitatively preserved even at $T_e$0, but the absence of higher-order kinetic moments prevents accurate Dimits shift prediction and transport suppression. The missing physics is associated with higher-order gyromoments, particularly parallel and perpendicular heat fluxes and the energy-weighted pressure tensor (Hoffmann et al., 18 Sep 2025).

6. Semi-analytical and analytic application-specific closures

In shock ignition, HEL closure is used in a diagnostic and parameterization sense. Semi-analytical approaches infer hot-electron temperature and conversion efficiency from measured $T_e$1 emission. The simplest model assumes a monoenergetic source and straight-line propagation,

$T_e$2

with the fitted range mapped to energy using ESTAR. For thick Cu targets, this gives $T_e$3 mg/cm$T_e$4, corresponding to $T_e$5 keV. The Harrach-Kidder model replaces this with a planar, isotropic Maxwellian source and an attenuation law

$T_e$6

For $T_e$7m, the fit gives $T_e$8 mg/cm$T_e$9 and $$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$0 keV; the inferred conversion efficiency for the highest-yielding targets is $$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$1 with upper limit $$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$2 (Afshari et al., 2018).

The same study states that hot electrons below $$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$3 keV predominantly deposit energy in the ablator and enhance shock formation, whereas higher-energy electrons can preheat the core and degrade compression. Under the reported third-harmonic irradiation conditions, the majority of hot electrons are below $$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$4 keV. In that sense, the semi-analytical closure translates experimental diagnostics into source terms for target design and integrated hydrodynamic calculations (Afshari et al., 2018).

A different application-specific HEL closure appears in non-equilibrium molecular plasmas with Treanor-Gordiets vibrational distributions. There the problem is not electron-ion or electron-phonon relaxation, but superelastic electron heating when the vibrational temperature exceeds the gas temperature. The unified analytic closure derives a thermodynamically consistent anharmonic gain function,

$$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$5

which enters the total superelastic heating rate

$$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$6

The formulation is based on detailed balance and a second-order Dunham expansion, identifies the kinetic crossover between vibrational-vibrational up-pumping and vibrational-translational relaxation, predicts the Treanor minimum, and recovers the accuracy of full state-to-state benchmarks at a fraction of the computational cost (Parent, 28 Jan 2026).

7. Validity limits, failure modes, and recurring disputes

The literature does not treat HEL closure as universally reliable. In ultrafast metals, its principal failure mode is incomplete electron thermalization: non-thermal electron distributions can persist on timescales comparable to electron-phonon relaxation, invalidating the assumption that the electron subsystem can always be represented by a single $$\frac{dT_i}{dt} = +\frac{g}{c_i}[T_e - T_i],\qquad \frac{dT_e}{dt} = -\frac{g}{c_e}[T_e - T_i],$$7 (Singh, 2024).

In electron-ion relaxation, a second failure mode is analytical rather than phenomenological. The general coupling-factor theory states that previous models such as Fermi’s golden rule can yield divergent results in the HEL if screening and many-body correlations are not included correctly. The dense-hydrogen analysis sharpens this objection by arguing that Coulomb-only coupled-mode calculations within RPA overestimate the coupled-mode contribution and that reduced models often fail basic sum-rule and Kramers-Krönig tests (Daligault et al., 2019, Dharma-wardana, 2017).

In transport closures, the failure mode is loss of nonlocal or nonlinear structure. The nonlocal heat-transport study shows that linearized stationary closures may reproduce damping rates for small perturbations while missing preheat and overpredicting peak heat flow for large temperature differences. The gyrokinetic study shows that a four-moment HEL closure can reproduce substantial Z-pinch physics yet fail to capture the Dimits shift and transport suppression in tokamak geometry because higher-order kinetic moments are absent (Brodrick et al., 2017, Hoffmann et al., 18 Sep 2025).

A common misconception is that HEL closure is synonymous with a single two-temperature model. The body of work summarized here suggests instead that HEL closure is a family of reductions whose validity depends on the subsystem decomposition, the interaction model, the analytical constraints imposed on the response functions, and the scale separation between fast electron equilibration and slower energy exchange. Under those conditions it can provide a rigorous bridge from microscopic kinetics to tractable macroscopic evolution equations; outside them, it can misrepresent the dominant relaxation channel or omit essential higher-order dynamics.