Collisional stopping power of ions in warm dense matter

A model for the collisional stopping of ions on free electrons in warm dense matter is developed and explored. It is based on plasma kinetic theory, but with modifications to address the warm dense matter regime. Specifically, it uses the Boltzmann-Uehling-Uhlenbeck kinetic equation to incorporate effects of Fermi degeneracy of electrons. The cross section is computed from quantum scattering of electrons and ions occuring via the potential of mean force derived from an average atom model, which incorporates effects of strong Coulomb correlations. Predictions from this model show comparable accuracy to results from time-dependent density functional theory calculations for deuterium near solid density and a temperature of several electronvolts, at a fraction of the computational cost. Further, the model captures the transition of a plasma from the classical limit to the degenerate limit, including qualitative behaviors of solid state theory.

💡 Research Summary

The paper presents a novel theoretical framework for calculating the collisional stopping power of ions moving through warm dense matter (WDM), a regime where both strong Coulomb coupling and electron degeneracy are simultaneously important. Traditional approaches either belong to the classical plasma domain—weak coupling, non‑degenerate electrons—and use analytic formulas such as the Bethe‑Bloch or Coulomb‑log models, or they belong to the condensed‑matter domain—strong coupling, fully degenerate electrons—and rely on time‑dependent density functional theory (TDDFT). The former fails when Γ≈1–30 and Θ≈0.1–1, while the latter becomes computationally prohibitive because the number of required electronic states grows dramatically at the moderate temperatures typical of WDM.

To bridge this gap, the authors combine two key ingredients: (1) the Boltzmann‑Uehling‑Uhlenbeck (BUU) kinetic equation, which extends the classical Boltzmann collision operator by including Pauli‑blocking factors (1 + δ θ f) that enforce Fermi‑Dirac statistics for electrons, and (2) a quantum‑mechanical potential of mean force (V_MF) derived from an average‑atom two‑component plasma model together with the Ornstein‑Zernike equation and a quantum hyper‑netted‑chain closure. V_MF represents the effective interaction between two particles after averaging over the remaining N‑2 particles in thermal equilibrium; it naturally incorporates static screening, exchange‑correlation, and strong‑coupling effects that are absent from simple Debye or Thomas‑Fermi screened Coulomb potentials.

Using V_MF, the authors compute quantum scattering phase shifts and differential cross sections for electron‑ion collisions. These quantum cross sections are inserted into the BUU collision term, yielding an expression for the ion energy loss per unit path length, dE/dx. The resulting stopping power contains an “effective Coulomb logarithm” L(v, Θ, Γ) that smoothly interpolates between the classical logarithmic behavior at high projectile velocities (v ≫ v_F) and the v⁻² scaling characteristic of a degenerate Fermi liquid at low velocities (v ≪ v_F). In the limit Θ → ∞ and Γ → 0 the model reproduces the standard Bethe‑Bloch result; in the opposite limit Θ → 0 it reduces to the Landau‑Fokker‑Planck low‑velocity stopping coefficient, including the Barkas correction that accounts for the asymmetry between positive and negative projectiles in a degenerate electron gas.

The authors validate the model numerically for hydrogen and deuterium plasmas at a density of 1.67 g cm⁻³ (≈10²³ cm⁻³ electrons) and a temperature of 5 eV, corresponding to Γ ≈ 4.6 and Θ ≈ 0.14—well within the WDM regime. Stopping powers computed with the BUU‑MF approach are compared to results from state‑of‑the‑art TDDFT simulations. Across a broad range of projectile velocities, the two methods agree within 5–10 %, with the BUU‑MF model capturing the non‑linear low‑velocity behavior that TDDFT predicts. Importantly, the computational cost of the BUU‑MF calculation is reduced by three to four orders of magnitude relative to TDDFT, making it feasible to generate extensive stopping‑power tables covering the full (density, temperature, projectile energy) parameter space required for inertial confinement fusion (ICF) hydrodynamic simulations.

Beyond benchmarking, the paper explores asymptotic limits analytically. In the classical, non‑degenerate limit (Γ ≪ 1, Θ ≫ 1) the model reduces to the familiar Coulomb‑log expression, confirming consistency with established plasma theory. In the fully degenerate limit (Θ → 0) the stopping power follows the v⁻² scaling of a Fermi liquid and reproduces the Landau‑Fokker‑Planck collision frequency, including the reduction of the effective collision frequency (the Barkas effect) due to Pauli blocking. The intermediate regime, where both Γ and Θ are of order unity, exhibits a smooth transition between these two extremes, demonstrating that the model correctly captures the combined influence of strong coupling and quantum degeneracy.

The authors argue that the BUU‑MF framework offers a practical, physically transparent, and computationally efficient tool for ICF modeling, where stopping power is a critical input for predicting burn propagation, alpha‑particle heating, and beam‑plasma interactions. Because the same kinetic equation can be extended to calculate other transport coefficients (electrical and thermal conductivity, viscosity), the approach promises a unified description of electron‑ion transport in WDM.

In conclusion, the paper introduces a robust intermediate‑theory that fills the methodological gap between classical plasma stopping models and fully quantum TDDFT calculations. By incorporating Fermi‑Dirac statistics via the BUU equation and realistic ion‑electron interactions via the quantum potential of mean force, the model delivers accurate stopping powers across the entire WDM regime while maintaining a computational cost compatible with large‑scale hydrodynamic simulations. This work therefore represents a significant step forward for predictive modeling of high‑energy‑density experiments and inertial confinement fusion research.