Maxwell–Jüttner Distribution Overview

• Maxwell–Jüttner distribution is the unique equilibrium momentum distribution for a relativistic ideal gas, incorporating special relativity into classical thermodynamics.
• Derived by maximizing Boltzmann–Gibbs entropy under energy and mass-shell constraints, it employs modified Bessel functions for normalization and moment evaluation.
• Its applications span relativistic plasmas, high-energy astrophysics, and numerical simulations, providing accurate models for particle dynamics at extreme temperatures.

The Maxwell–Jüttner distribution is the unique equilibrium single-particle distribution for a classical ideal gas in the relativistic regime, generalizing the Maxwell–Boltzmann law to accommodate the kinematics and thermodynamics dictated by special relativity. Its domain spans relativistic kinetic theory, high-energy astrophysics, plasma physics, and statistical mechanics. Its properties, generalizations, and numerical implementation are foundational in accurate modeling of relativistic gases and plasmas in both laboratory and astrophysical contexts.

1. Definition, Norms, and Functional Form

The Maxwell–Jüttner (MJ) distribution describes the equilibrium momentum or velocity distribution for a dilute gas of classical, identical, non-interacting particles at temperature $T$. In the fiducial rest frame, with particle rest mass $m$ and speed of light $c$, denote the Lorentz factor as

$\gamma = \frac{1}{\sqrt{1 - (v/c)^2}}$

and define the dimensionless temperature

$\theta = \frac{k_B T}{mc^2}$

where $k_B$ is Boltzmann’s constant. The canonical one-dimensional MJ distribution in $\gamma$ reads

$$f(\gamma) = A \gamma \sqrt{\gamma^2 - 1} \exp\bigl(-\gamma / \theta\bigr), \quad \gamma \geq 1,$$

where the normalization constant is

$A = \frac{1}{\theta K_2(1/\theta)}$

and $K_2$ is the modified Bessel function of the second kind (Zenitani, 2024, Zaninetti, 2020). The $d$-dimensional velocity form is

$$f_\mathrm{MJ}(\mathbf v) = A(T) \gamma^{d+2} \exp\left[-\frac{\gamma m c^2}{k_B T}\right]$$

with normalization chosen via

$$A(T) = \left[ S_{d-1} \int_0^c v^{d-1} \gamma^{d+2} \exp\left(-\frac{\gamma m c^2}{k_B T}\right) dv \right]^{-1},$$

where $S_{d-1}$ is the angular measure (Mendoza et al., 2012, Curado et al., 2022).

2. Derivation and Theoretical Context

The MJ distribution is obtained by maximizing the Boltzmann–Gibbs entropy under mass-shell and energy constraints, or by factorizing the single-particle canonical Gibbs measure,

$$f(\mathbf p) \propto \exp\left(-\beta p^\mu U_\mu\right)$$

where $\beta = 1/(k_B T)$, $p^\mu = (\gamma m c, \mathbf p)$, and $U^\mu = (c, 0, 0, 0)$ in the rest frame. The distribution is properly Lorentz covariant when written in terms of the invariant momentum-space measure, but the standard MJ form, $$f(\mathbf v) \propto \gamma^{d+2} e^{-\gamma/\theta}$$, is only manifestly invariant under spatial rotations, not Lorentz boosts (Damião et al., 2022, Curado et al., 2022, Johnson et al., 21 Mar 2025).

The relativistic partition function is constructed as

$$Z_1(\beta) = V \int \frac{d^3p}{(2\pi)^3} \exp(-\beta E(p))$$

where $E(p) = \sqrt{p^2c^2 + m^2c^4}$. Normalization constants and all thermal moments can be expressed in terms of Bessel functions, e.g.,

$$Z_1(\beta) = V \frac{1}{2\pi^2} \int_{m}^{\infty} E \sqrt{E^2 - m^2} e^{-\beta E} dE$$

(Moradpour et al., 2021, Zaninetti, 2020).

3. Properties, Limits, and Moments

3.1 Qualitative Behavior

The MJ distribution exhibits a single pronounced peak for low temperatures ($\theta \ll 1$), with decay at high energies following a heavy but exponential tail,

$f(\gamma) \sim \gamma^{3/2} \exp(-\gamma/\theta)$

for $\gamma \gg \theta$ (Zenitani, 2024). At ultrarelativistic temperatures ($\theta \gg 1$), the distribution becomes broad, with the typical velocity approaching $c$ (Damião et al., 2022, Zaninetti, 2020).

3.2 Connection to Maxwell–Boltzmann Law

In the nonrelativistic limit ($\theta \ll 1$), expanding $\gamma$ and the exponential reproduces the Maxwell–Boltzmann result to leading order. For $\theta \to 0$, contributions above $v \sim 0$ are strongly suppressed, and all thermal moments take classical values.

3.3 Moments

Moments are efficiently evaluated with respect to the $\gamma$ representation. The mean energy is

$$\langle E \rangle = m c^2 \left[ \frac{K_1(1/\theta)}{K_2(1/\theta)} + 3\theta \right ]$$

while the variance follows from higher-order integrals involving $K_0, K_1, K_2$. In the ultrarelativistic limit, $\langle E \rangle \approx 3k_B T$ per particle (Zaninetti, 2020, Johnson et al., 21 Mar 2025).

4. Dimensionality, Geometry, and Lorentz-Invariance

Special relativity imposes a non-Euclidean (Lobachevsky) structure on velocity and momentum spaces. The canonical MJ prefactor, $\gamma^{d+2}$, arises from imposing the flat measure $d^d v$, whereas the invariant measure is $\gamma^{d+1} d^d v$ (Curado et al., 2022). The extra power is a relic of the entropy maximization procedure performed on non-invariant velocity-space; as a result, the standard MJ form is not Lorentz-invariant. This issue becomes manifest when considering rapidity space: at high temperature, the MJ rapidity distribution acquires a curvature change and becomes bimodal, contrasting with Lorentz-invariant frameworks and numerical dynamics data favoring unimodal rapidity PDFs at all temperatures.

Simulations indicate that for a truly Lorentz-invariant one-dimensional velocity PDF (LID), one must distribute rapidities rather than velocities, with the invariant measure and additivity of rapidity ensuring form invariance under boosts (Curado et al., 2022).

5. Bimodality Transition and Critical Temperature

Under increasing temperature, the MJ distribution (especially in velocity space) exhibits a morphological transition from single-peaked to double-peaked (“bimodal”) form. The critical temperature for this transition, in $d$ spatial dimensions, is

$k_B T_c = \frac{m c^2}{d + 2}$

For $d = 1,2,3$, this gives $T_c = m c^2/3 k_B$, $m c^2/4 k_B$, and $m c^2/5 k_B$, respectively (Mendoza et al., 2012). The transition is second-order (continuous) with respect to temperature at zero net flow, with order parameter proportional to $(T - T_c)^{1/2}$. Introduction of finite mean velocity drives the transition first-order, with discontinuous jump in the order parameter (Mendoza et al., 2012).

A key implication is that, for sufficiently hot relativistic gases, the most probable particle states lie near $|v| \to c$, with populations clustering near the velocity boundary in the form of beams, rings, or hollow shells depending on spatial dimension.

6. Generalizations: Anisotropy, Spin, Rotation, and Quantum Corrections

6.1 Anisotropic Distributions

The MJ law admits covariant generalization to anisotropic temperature tensors. Introducing a second-rank inverse-temperature tensor $\Theta^{\mu\nu}$, yields

$$f(x,p)\,\sqrt{-g} = C \exp\bigl[ -\beta_\perp \sqrt{1 + (p_\mu \Theta^{\mu\nu} p_\nu)/(m^2 c^2)} \bigr] \sqrt{-g}$$

where $\beta_\perp = m c^2/T_\perp$, and the distribution reduces to the isotropic MJ law for $T_\perp = T_\parallel$ (Treumann et al., 2015).

6.2 Spinning and Rotating Gases

The equilibrium distribution of a relativistic gas with spin and global rotation is a direct generalization: $$f(x,p,s) \propto \exp\Bigl\{ -\frac{mc^2 \gamma}{k_B T} + \dots + \frac{\hbar \Sigma}{k_B T} (\text{spin-rotation terms}) \Bigr\}$$ where $\Sigma$ is the spin magnitude, and corrections to the partition function, thermodynamic potentials, and spin polarization arise at $O(\Omega)$ (angular velocity) and $O(\Sigma)$. In the nonrotating, spinless limit, it reduces to MJ (Kaparulin et al., 2023).

6.3 Quantum Gravity: GUP/EUP Corrections

Generalized and extended uncertainty principles (GUP/EUP) can deform the MJ density of states. GUP modifies the fundamental phase-space cell, introducing an extra high-momentum suppression factor and reducing high-temperature densities and energies; EUP leaves the momentum measure unaffected, so the MJ distribution itself is unchanged (Moradpour et al., 2021).

7. Numerical Sampling and Simulation Implementation

Sampling from the MJ distribution in particle-based simulations is nontrivial due to the heavy tail and non-analytic inversion of its cumulative. The piecewise rejection sampling algorithm of Swisdak (2013), as modified for full self-containment, constructs an envelope from exponential and linear pieces tangential to the PDF at key $e$-folding points and the maximizer of $f(p)/p$. Approximations for boundary points and slopes, and a linear-slope left envelope, yield a root-finder-free, efficient scheme with acceptance rates up to 93%, representing a ~20% reduction in wasted draws compared to previous methods (Zenitani, 2024).

For grid-based kinetic codes (e.g., discontinuous Galerkin methods), direct projection of MJ onto a finite domain leads to truncation and aliasing errors. An iterative moment-correction algorithm adjusts the input moments and reprojects until density, energy, and bulk velocity match target values to machine precision. Nonlinear operations, Lorentz boost factors, and other quantities are handled via weak operations and numerical quadrature (Johnson et al., 21 Mar 2025).

8. Applications and Physical Relevance

The MJ distribution underpins the kinetic theory of relativistic plasmas, models for high-energy astrophysical environments (e.g., radio sources, quasar jets), and is essential for the construction of relativistic BGK models. Its moments and normalization constants enter directly into analytic expressions for synchrotron emissivity and spectral fitting of observed nonthermal emission (Zaninetti, 2020).

The distinction between MB and MJ becomes significant when the dimensionless temperature $\theta$ approaches $\gtrsim 0.01$. Errors in characteristic speed and density between MB and MJ rapidly escalate with increasing $T$, only the MJ respects the relativistic velocity bound and provides physically meaningful predictions in this domain (Damião et al., 2022).

• (Zenitani, 2024) Modifications to Swisdak (2013)'s rejection sampling algorithm for a Maxwell-Jüttner distribution in particle simulations
• (Kaparulin et al., 2023) Generalized Maxwell-Juttner distribution for rotating spinning particle gas
• (Mendoza et al., 2012) Single to Double Hump Transition in the Equilibrium Distribution Function of Relativistic Particles
• (Moradpour et al., 2021) A Note on Effects of Generalized and Extended Uncertainty Principles on Jüttner Gas
• (Damião et al., 2022) The Relativistic Maxwell-Jüttner Velocity Distribution Function
• (Johnson et al., 21 Mar 2025) Discontinuous Galerkin Representation of the Maxwell-Jüttner Distribution
• (Treumann et al., 2015) Anisotropic Jüttner (relativistic Boltzmann) distribution
• (Curado et al., 2022) Relativistic Gas: Invariant Lorentz Distribution for the velocities
• (Zaninetti, 2020) New probability distributions in astrophysics: IV. The relativistic Maxwell-Boltzmann distribution