A statistical theory of electronic degrees of freedom in wave packet molecular dynamics

We derive statistical distributions for the degrees of freedom in wave packet molecular dynamics models. Specifically, a theory is developed for the width distributions of Gaussian wavepackets in both isotropic and anisotropic formulations. The resulting distribution functions show good agreement with molecular dynamics data under warm dense matter conditions, providing practical guidance for constraining the confining potential, an empirical parameter in the model. We also discuss how these distributions influence the resulting effective Coulomb interactions.

💡 Research Summary

The paper presents a rigorous statistical framework for describing the electronic degrees of freedom in Wave‑Packet Molecular Dynamics (WP‑MD), a method that treats electrons as Gaussian wave packets characterized by a center, momentum, and a width (σ). While WP‑MD has become popular for simulating warm dense matter (WDM) and related high‑energy‑density systems, the width parameter has traditionally been set empirically through a confining potential V_conf(σ)=½kσ², with the stiffness constant k adjusted to match limited benchmark data. This empirical approach limits predictive power, especially under extreme temperature and density conditions where experimental validation is scarce.

The authors start by placing the ensemble of electronic wave packets in a microcanonical (or canonical) statistical‑mechanical setting. Treating σ as an independent generalized coordinate, they derive its probability density from the Boltzmann factor exp(‑βH), where the Hamiltonian includes the kinetic contribution of the electrons and the harmonic confining term. For the isotropic case, where all spatial directions share the same width, the Jacobian of the transformation to σ‑space contributes a factor σ^{d‑1} (d = 3 for three‑dimensional systems). The resulting distribution is

P(σ) ∝ σ^{d‑1} exp(‑βkσ²/2).

This is essentially a chi‑square distribution with a temperature‑dependent scale set by the stiffness k. The mean width ⟨σ⟩ = √