Starrett and Saumon (High Energy Density Physics 10, 35, 2014) introduced a computationally cheap, self-consistent method for the ionic structure of warm dense matter: a finite-temperature density-functional average-atom model coupled to classical-fluid integral equations (later also to classical molecular dynamics with pseudoatom pair potentials). The method class carries stated limitations: spherical symmetry and pair interactions that exclude chemical bonding and molecule formation, the cost of the quantum version, and restriction to single-species plasmas. Follow-up work by the same group and others resolved the multi-species restriction (2014-2015) and largely mitigated the cost problem, and the method has been validated against ab initio simulations and x-ray scattering experiments for simple-metal-like plasmas. The surviving open problem, reaffirmed as current state of the art by the 2026 Roadmap for warm dense matter physics (Plasma Physics and Controlled Fusion 68, 073501, arXiv 2025), is to determine the ionic structure of warm dense matter within this inexpensive self-consistent model class in regimes where chemical bonding and molecular formation significantly shape that structure.

Background

Warm dense matter (WDM) is the regime of temperatures of order 1-100 eV and densities around solid density where ions, free electrons, and bound electrons are all moderately to strongly coupled; it occurs in inertial confinement fusion (ICF) capsules, giant-planet interiors, and laser-shocked solids. A standard observable of this regime is the ion-ion pair distribution function gII(r)g_{II}(r) (the probability, relative to an ideal gas, of finding two ions a distance rr apart) and its Fourier transform, the static structure factor, which x-ray diffraction and x-ray Thomson scattering (XRTS) measure.

At one end of the modeling spectrum sit ab initio simulations - Kohn-Sham density-functional-theory molecular dynamics (DFT-MD) and path-integral Monte Carlo (PIMC) - which treat nuclei and electrons in a simulation box and capture bonding, clusters, and all correlations, but at a computational cost that grows steeply with temperature and system size. At the other end sit average-atom (AA) models: the nucleus of charge ZZ is placed at the center of a sphere (the ion sphere, whose radius is set by the mean ion density) and surrounded by ZZ electrons treated in finite-temperature (Mermin) density-functional theory, with the surrounding plasma represented as a spherically averaged background. The AA model yields the average electronic structure and mean ionization of an ion in the plasma at negligible cost, but by construction knows nothing about the spatial arrangement of the surrounding ions.

Starrett and Saumon closed this gap self-consistently (Phys. Rev. E 85, 026403 (2012); 87, 013104 (2013)): the screening electron cloud around the central ion, computed by the AA model, supplies an effective ion-ion pair potential, which feeds classical-fluid integral equations - the Ornstein-Zernike (OZ) equations closed in the hypernetted-chain (HNC) approximation, the workhorse closure of classical liquid-state theory - whose output, the pair distribution functions and correlation functions of the surrounding two-component (ion plus free-electron) fluid, is fed back into the AA model until the whole coupled nonlinear system converges. The 2014 High Energy Density Physics paper examined here presents this machinery as a simple method for the ionic structure of WDM, validates gII(r)g_{II}(r) against DFT-MD, PIMC, and Thomas-Fermi molecular dynamics for H, C, Al, and Fe, and later variants replaced the integral equations by classical molecular dynamics driven by the same pseudoatom pair potentials (pseudoatom molecular dynamics, PAMD), which also yields equations of state (EOS) and ion transport coefficients. For simple-metal-like plasmas the method is two to three orders of magnitude cheaper than DFT-MD, which is what makes it attractive for wide-range EOS tables, XRTS interpretation, and opacity and transport calculations.

The method's stated limitation traces to its structural assumption of spherical symmetry and pairwise interactions: the pseudoatom potential between two ions is a function of their separation only, so directional (covalent) bonds, transient bonding, and molecule formation - precisely the physics of dense molecular hydrogen, hydrocarbons at a few eV, water and silicate mixtures at planetary conditions, and compressed silicon near its liquid-liquid transitions - cannot be represented. Because the pseudoatom framework replaces the surrounding plasma by a spherically averaged field, the general question of how to build a comparably cheap, self-consistent ionic-structure theory when the ions themselves form non-spherical, short-lived bonded units remains open in the field.

Problem Statement

The problem is the surviving open core of the limitations stated by Starrett and Saumon (2014) for their self-consistent ionic-structure method for warm dense matter. The primary objective: extend the average-atom-type model class - a finite-temperature density-functional average atom coupled self-consistently to a classical-fluid description of the surrounding plasma, via Ornstein-Zernike-type integral equations (the classical liquid-state equations relating pair correlation functions to effective pair interactions) or pseudoatom pair potentials (effective ion-ion potentials generated by the average atom's screening cloud) - so that it determines the ionic structure (ion-ion pair distribution functions and static structure factors) of warm dense matter with quantitative, validated accuracy in temperature and density regimes where chemical bonding, transient bonding, or molecule formation significantly affect that structure, while retaining the orders-of-magnitude cost advantage over Kohn-Sham molecular dynamics that is the reason for the model class's existence. This requires confronting the structural restriction the source identified: a spherically symmetric, pairwise pseudoatom cannot represent directional bonds or molecular units, so a solution must supply what replaces it - for example a closure or effective-interaction theory for bonded, many-center, short-lived correlations, together with a defensible definition of the 'ion' or bonded unit whose structure is being computed. Answers are to be validated against ab initio simulation (Kohn-Sham molecular dynamics, path-integral Monte Carlo) and, where available, x-ray scattering and diffraction experiments, at conditions where bonding is significant. Dense partially molecular hydrogen, hydrocarbon or water mixtures at a few eV, and shock-compressed silicon in the 100-200 GPa regime are representative cases, not an exhaustive list; the objective is the general theory, and a result on any single case does not resolve it. The two remaining items of the source's limitation list - single-species scope and the cost of the quantum version - were substantially resolved by 2014-2025 follow-up work; they are not part of the open problem. The common objective is a reliable, inexpensive, self-consistent ionic-structure theory of warm dense matter valid across the bonding and non-bonding regimes alike.

References

  1. C. E. Starrett and D. Saumon, 'A simple method for determining the ionic structure of warm dense matter', High Energy Density Physics 10, 35-42 (2014). DOI: 10.1016/j.hedp.2013.12.001, https://doi.org/10.1016/j.hedp.2013.12.001
  2. C. E. Starrett and D. Saumon, 'Electronic and ionic structures of warm and hot dense matter', Physical Review E 87, 013104 (2013). DOI: 10.1103/PhysRevE.87.013104, https://doi.org/10.1103/PhysRevE.87.013104
  3. D. Saumon, C. E. Starrett, J. A. Anta, W. Daughton, and G. Chabrier, 'The structure of warm dense matter modeled with an average atom model with ion-ion correlations', in proceedings of the Computational Challenges in Warm Dense Matter Workshop (Springer), arXiv:1311.3986 (2013). arXiv:1311.3986, https://arxiv.org/abs/1311.3986
  4. C. E. Starrett, D. Saumon, J. Daligault, and S. Hamel, 'Integral equation model for warm and hot dense mixtures', Physical Review E 90, 033110 (2014). DOI: 10.1103/PhysRevE.90.033110, arXiv:1408.3548, https://arxiv.org/abs/1408.3548
  5. C. E. Starrett, J. Daligault, and D. Saumon, 'Pseudoatom molecular dynamics', Physical Review E 91, 013104 (2015). DOI: 10.1103/PhysRevE.91.013104, https://doi.org/10.1103/PhysRevE.91.013104
  6. C. E. Starrett and D. Saumon, 'Models of the elastic x-ray scattering feature for warm dense aluminum', Physical Review E 92, 033101 (2015). DOI: 10.1103/PhysRevE.92.033101, https://doi.org/10.1103/PhysRevE.92.033101
  7. C. E. Starrett and D. Saumon, 'Equation of state of dense plasmas with pseudoatom molecular dynamics', Physical Review E 93, 063206 (2016). DOI: 10.1103/PhysRevE.93.063206, https://doi.org/10.1103/PhysRevE.93.063206
  8. A. A. Ovechkin, P. A. Loboda, A. L. Falkov, and F. A. Sapozhnikov, 'Equation of state modeling with pseudoatom molecular dynamics', Physical Review E 103, 053206 (2021). DOI: 10.1103/PhysRevE.103.053206, https://doi.org/10.1103/PhysRevE.103.053206
  9. T. J. Callow, S. B. Hansen, E. Kraisler, and A. Cangi, 'First-principles derivation and properties of density-functional average-atom models', Physical Review Research 4, 023055 (2022). DOI: 10.1103/PhysRevResearch.4.023055, arXiv:2103.09928, https://arxiv.org/abs/2103.09928
  10. R. Piron, 'Atomic models of dense plasmas, applications and current challenges', Atoms 12, 26 (2024). DOI: 10.3390/atoms12040026, arXiv:2311.15418 (2023), https://arxiv.org/abs/2311.15418
  11. J. Vorberger, F. Graziani, D. Riley, et al., 'Roadmap for warm dense matter physics', Plasma Physics and Controlled Fusion 68, 073501 (2026). DOI: 10.1088/1361-6587/ae672c, arXiv:2505.02494 (2025), https://arxiv.org/abs/2505.02494
  12. M. W. C. Dharma-wardana, D. D. Klug, H. Poole, and G. Gregori, 'Ionic structure, liquid-liquid phase transitions, x-ray diffraction, and x-ray Thomson scattering in shock-compressed liquid silicon in the 100-200 GPa regime', Physical Review E 111, 015205 (2025). DOI: 10.1103/PhysRevE.111.015205, https://doi.org/10.1103/PhysRevE.111.015205