Two-Temperature Kinetic Equations

Updated 8 February 2026

Two-temperature kinetic equations are a theoretical framework distinguishing separate energy modes with distinct relaxation times in nonequilibrium systems.
They are derived from the Boltzmann equation using a multiscale Chapman–Enskog expansion, incorporating both resonant and inelastic collision operators.
These models underpin hydrodynamic closures in applications such as hypersonic flows, plasmas, and dense shockwave phenomena while ensuring energy conservation and an H-theorem.

Two-temperature kinetic equations constitute a theoretical framework for nonequilibrium systems in which distinct relaxation time scales exist for different energy modes—typically, translational/translational-rotational and internal (rotational, vibrational, or electronic) degrees of freedom. These equations underpin the derivation of hydrodynamic models with separate mode temperatures and explicit relaxation terms, offering a rigorous alternative to the classical single-temperature kinetic theory for polyatomic gases and gas mixtures. Their development connects collisional kinetic theory, multiscale Chapman–Enskog expansions, and applications ranging from hypersonic flow to plasmas and condensed-matter systems.

1. Fundamental Structure of Two-Temperature Kinetic Equations

The two-temperature framework originates from the Boltzmann equation, in which the single-particle distribution $f=f(t,\mathbf{x},\boldsymbol{\xi},I)$ —with $\mathbf{x}\in\mathbb{R}^3$ , $\boldsymbol{\xi}\in\mathbb{R}^3$ (particle velocity), and $I\geq 0$ (internal energy)—obeys

$\partial_t f + \boldsymbol{\xi}\cdot \nabla_{\mathbf{x}} f = Q_\theta(f, f),$

where $Q_\theta$ is a composite collision operator that interpolates between inelastic and resonant (elastic) collisions,

$Q_\theta(f, f) = \theta Q_s(f, f) + (1-\theta) Q_r(f, f).$

Here, $Q_s$ is the standard inelastic operator (with translational–internal exchanges), and $Q_r$ accounts for resonant collisions that conserve both kinetic and internal energies separately. The parameter $\theta$ modulates the strength of mode coupling, and $\delta$ denotes the number of internal degrees of freedom (Aoki et al., 14 Oct 2025).

The fundamental distinction is that resonant (or quasi-resonant) collisions are modeled not as rare, but as dominant, with weak inelasticity $\theta\ll 1$ serving as a small parameter.

Macroscopic fields (number density $n$ , bulk velocity $\mathbf{u}$ , translational energy $e_\mathrm{tr}$ , internal energy $e_\mathrm{int}$ , mode temperatures $T_\mathrm{tr}$ , $T_\mathrm{int}$ ) are defined as velocity or internal-energy moments of $f$ . For instance,

$T_\mathrm{tr} = \frac{2}{3k_B} e_\mathrm{tr}, \qquad T_\mathrm{int} = \frac{2}{\delta k_B} e_\mathrm{int},$

with total temperature $T = [3 T_\mathrm{tr} + \delta T_\mathrm{int}]/(3+\delta)$ .

The local quasi-equilibrium (in the regime of dominant resonant collisions) is a bi-Maxwellian,

$M_r = n\, I^{\delta/2-1} (2\pi k_B T_\mathrm{tr}/m)^{-3/2} (k_B T_\mathrm{int})^{-\delta/2} / \Gamma(\delta/2) \ \exp\left[ -\frac{m|\boldsymbol{\xi} - \mathbf{u}|^2}{2k_B T_\mathrm{tr}} - \frac{I}{k_B T_\mathrm{int}} \right].$

2. Chapman–Enskog Expansion and Multiscale Regimes

To systematically derive two-temperature hydrodynamics, the Chapman–Enskog method is applied with both Knudsen number $\varepsilon$ and the small inelasticity parameter $\theta\ll 1$ . Two low-coupling regimes emerge:

Very weak coupling: $\theta = \kappa \varepsilon^2$ (inelasticity enters at Navier–Stokes order).
Moderately weak coupling: $\theta = \bar{\kappa} \varepsilon$ (mode coupling enters already at Euler order).

The distribution is expanded as $f = f^{(0)} + \varepsilon f^{(1)} + \varepsilon^2 f^{(2)} + \cdots$ , with $f^{(0)} = M_r$ the two-temperature Maxwellian, and solubility conditions ensuring the orthogonality of corrections to collision invariants (Aoki et al., 14 Oct 2025).

At Euler order, the system closes on mass, momentum, and two separate energy balances (translational and internal), coupled by a relaxation source term $R$ . For $T_\mathrm{tr}\neq T_\mathrm{int}$ , $R$ describes Landau–Teller-type relaxation, scaling as either $O(\varepsilon^2)$ (very weak) or $O(\varepsilon)$ (moderately weak) depending on the regime.

3. Hydrodynamic Limits and Constitutive Laws

The Euler-level two-temperature model has the structure: $\begin{aligned} & \partial_t \rho + \nabla \cdot (\rho \mathbf{u}) = 0, \ & \partial_t (\rho \mathbf{u}) + \nabla \cdot (\rho \mathbf{u} \otimes \mathbf{u} + p) = 0, \ & \partial_t \mathcal{E}_\mathrm{tr} + \nabla \cdot (\cdots) = -R, \ & \partial_t \mathcal{E}_\mathrm{int} + \nabla \cdot (\cdots) = +R, \end{aligned}$ with $R$ determined by the coupling regime,

$R = \begin{cases} 0 + O(\varepsilon^2), & \theta = O(\varepsilon^2) \ -\kappa F(\rho, T_\mathrm{tr}, T_\mathrm{int}) (T_\mathrm{tr}-T_\mathrm{int}), & \theta = O(\varepsilon) \end{cases}$

for an explicit $F(\cdot)$ depending on collision kernel exponents $(\alpha, \beta)$ . At the Navier–Stokes level, corrections to the stress, translational and internal heat fluxes, and to the relaxation source (via a linearized collision operator) are constructed, featuring explicit positive-definite integral representations for transport and exchange coefficients (Aoki et al., 14 Oct 2025).

The entire system is closed, reduces to single-temperature hydrodynamics when $T_\mathrm{tr}=T_\mathrm{int}$ , and possesses an H-theorem structure at the kinetic level.

4. Comparison with Alternative Two-Temperature Kinetic Models

BGK/ES-BGK Approaches

In BGK-type models, two-temperature dynamics are encoded via relaxation to a composite Maxwellian with two differing temperatures, with further source terms characterizing the relaxation between $T_\mathrm{tr}$ and $T_\mathrm{int}$ or, in mixtures, between the species' temperatures. These models provide conservation of mass, momentum, energy, H-theorems, and positive decay to a global equilibrium with $T_\mathrm{tr}=T_\mathrm{int}$ (Pirner, 2018, Klingenberg et al., 2018).

Fokker–Planck Models

Fokker–Planck approximations also permit two-temperature closures, with drift-diffusion terms in both translational and internal variables, and collision frequencies (including cross-couplings) engineered for entropy dissipation and positivity of all temperatures (Pirner, 2024, Pirner, 2024).

Quasi-Resonant Boltzmann Kinetics

Recent works introduce quasi-resonant kernels that interpolate between strictly resonant and standard inelastic processes. This yields initial rapid relaxation towards a two-temperature Maxwellian (at $O(1)$ time), followed by slow, Landau–Teller-type equilibration of $T_\mathrm{tr}$ and $T_\mathrm{int}$ (at $O(\varepsilon^{-2})$ time). The model possesses a rigorous H-theorem and demonstrates by simulation that the system remains close to such a multi-temperature quasi-equilibrium throughout the relaxation process (Borsoni et al., 4 Jun 2025).

5. Physical Realizations and Applications

Two-temperature kinetic equations are critical in describing systems such as:

Polyatomic gases subject to weak collisional exchange between modes, especially in high-speed (hypersonic) flows involving thermal nonequilibrium in shock-layers (Aoki et al., 14 Oct 2025, Gao et al., 3 Jan 2026).
Plasmas with separation of electron and heavy particle temperatures, where mass disparity ensures slow electron–ion energy exchange and two-temperature closure is essential for both kinetic and fluid-level plasma modeling (Orlac'H et al., 2017).
Dense fluid shockwaves and MD-informed models, where strong spatial anisotropy leads to distinct longitudinal and transverse temperatures, necessitating a two-temperature coupled continuum system for accurate macroscopic shock-shape prediction (Hoover et al., 2010).
Bose–Einstein condensate mixtures, where inter-component and intra-component kinetic processes must accommodate distinct mode and component temperatures, and explicit collision integrals ensure the proper two-temperature relaxation and hydrodynamics (Edmonds et al., 2014).

6. H-Theorem, Positivity, and Structural Properties

All reputable two-temperature kinetic models guarantee mass, momentum, and energy conservation. The equilibrium state is a Maxwellian with coincident mode temperatures and velocities. The structural features include:

H-theorem: Monotonic decay of entropy, with equality only at global equilibrium with all $T_\mathrm{tr}=T_\mathrm{int}$ (and $u_1=u_2$ in mixtures) (Aoki et al., 14 Oct 2025, Pirner, 2018, Borsoni et al., 4 Jun 2025).
Positivity of all mode temperatures: Ensured by parameter choices in both BGK and Fokker–Planck models and, where appropriate, by the positivity of the collision operators (Klingenberg et al., 2018, Pirner, 2024).
Closure relations: Explicit constitutive relations at the hydrodynamic level, with all transport and relaxation coefficients derivable from the underlying kinetic model.

7. Multiscale Nature, Regime Classification, and Model Selection

The regime distinction—determined by the scaling of the inelastic coupling parameter $\theta$ relative to the Knudsen number—leads to qualitatively different behaviors:

Regime	Coupling Parameter	Source $R$ scaling	Physical Consequence
Very weak coupling (A)	$\theta = \kappa\varepsilon^2$	$R=O(\varepsilon^2)$	$T_\mathrm{tr}$ , $T_\mathrm{int}$ convected almost independently (Euler order decoupled)
Moderately weak coupling (B)	$\theta = \bar\kappa\varepsilon$	$R=O(1)$	Mode temperatures equilibrate at convective time scale

This classification informs the choice of model for simulation of rarefied nonequilibrium flows, hypersonic regimes, or the analysis of relaxation phenomena in gases and plasmas. The multiscale approach ensures that only relevant couplings contribute at any given order, and that macroscopic closures consistently reflect kinetic multitemperature effects (Aoki et al., 14 Oct 2025, Pirner, 2018).

References

"Two-temperature fluid models for a polyatomic gas based on kinetic theory for nearly resonant collisions" (Aoki et al., 14 Oct 2025)
"A kinetic model for polyatomic gas with quasi-resonant collisions leading to bi-temperature relaxation processes" (Borsoni et al., 4 Jun 2025)
"A BGK model for gas mixtures of polyatomic molecules allowing for slow and fast relaxation of the temperatures" (Pirner, 2018)
"A consistent kinetic model for a two-component mixture of polyatomic molecules" (Klingenberg et al., 2018)
"Kinetic theory of two-temperature polyatomic plasmas" (Orlac'H et al., 2017)
"Well-Posed Two-Temperature Constitutive Equations for Stable Dense Fluid Shockwaves using Molecular Dynamics …" (Hoover et al., 2010)
"A two-temperature gas-kinetic scheme for hypersonic nonequilibrium flow computations" (Gao et al., 3 Jan 2026)
"A consistent non-linear Fokker-Planck model for a gas mixture of polyatomic molecules" (Pirner, 2024)
"Kinetic Model of Trapped Finite Temperature Binary Condensates" (Edmonds et al., 2014)