Validated first-principles radiative opacity of warm dense hydrogen: extending the opacity table beyond the degenerate regime and benchmarking it experimentally

Warm dense hydrogen isotopes are the working fluid of inertial confinement fusion (ICF) and a benchmark case for warm dense matter (WDM) theory, yet their radiative opacity in the strongly coupled, electron-degenerate regime is still neither completely calculated nor experimentally validated. A first-principles opacity table (FPOT) for deuterium, built from quantum molecular dynamics combined with the Kubo-Greenwood linear-response formula, covers densities of 0.5-673.5 g/cm^3 only up to the Fermi temperature and exceeds the tabulated astrophysical opacity used in hydrocodes by factors of 3-100 at ICF-relevant photon energies; radiation-hydrodynamics simulations with these opacities shift predicted neutron yields by tens of percent to a factor of about 2.5. The open problem is to complete the first-principles opacity description of dense hydrogen - beyond the Fermi temperature, with quantified uncertainties including collisional corrections to the independent-particle response - and to validate it against a laboratory measurement, which does not yet exist for hydrogen at any warm dense condition.

Background

Radiative opacity quantifies how a plasma absorbs and scatters photons, and radiation transport in hydrodynamic simulations is governed by the Rosseland mean opacity, a weighted harmonic average of the spectral opacity over photon energy. In inertial confinement fusion (ICF), laser-driven implosion of a cryogenic deuterium-tritium (DT) capsule compresses the fuel through the warm dense matter (WDM) regime - densities from below solid density to several hundred g/cm^3 at temperatures of a few eV up to the electron Fermi temperature $T_F$ (the temperature scale above which electrons cease to be quantum-degenerate) - where strong Coulomb coupling and degeneracy invalidate classical plasma models. Because the DT opacity controls how much radiation from the hot coronal plasma is deposited in the fuel, ICF target-design predictions are sensitive to it.

The standard inputs were atomic-physics tables built for astrophysical plasmas, such as the Los Alamos astrophysical opacity table (AOT) of Huebner and coworkers (1977) used in ICF hydrocodes. Hu and collaborators instead derived a first-principles opacity table (FPOT) for deuterium (S. X. Hu, L. A. Collins, V. N. Goncharov, T. R. Boehly, R. Epstein, R. L. McCrory, and S. Skupsky, Phys. Rev. E 90, 033111 (2014)): quantum molecular dynamics (QMD) - molecular dynamics with electrons treated by Kohn-Sham density-functional theory - samples ionic configurations, the Kubo-Greenwood formula (linear-response theory evaluating the frequency-dependent conductivity $\sigma(\omega)$ from the Kohn-Sham orbitals) yields the optical response, and Kramers-Kronig relations convert it into the refractive index and frequency-grouped Rosseland mean opacities, with DT obtained by mass scaling (deuterium-tritium opacities follow from the deuterium calculations by scaling the ionic mass, on which the electronic structure depends only weakly). Covering $\rho = 0.5$-$673.5\,\mathrm{g/cm^3}$ and $T = 5000\,\mathrm{K}$ up to $T_F$ at each density, the FPOT exceeded the AOT by factors of 3-100 in the ICF-relevant photon-energy range under strongly coupled, degenerate conditions (illustrated at $5.388\,\mathrm{g/cm^3}$, $10.8\,\mathrm{eV}$), converging to the AOT only as classical-plasma conditions were approached (the $43\,\mathrm{eV}$ comparison). Implementing FP-based DT properties (equation of state, thermal conductivity, opacity) in radiation-hydrodynamics codes changed predicted ICF neutron yields by roughly 10% up to a factor of about 2.5, with the largest variation at the lowest fuel adiabats (the adiabat sets the fuel entropy during compression; lower-adiabat implosions reach more strongly degenerate conditions) (S. X. Hu et al., Phys. Plasmas 22, 056304 (2015)).

Later work sharpened the limitations rather than removing them. A 2024 community review of first-principles simulations of dense hydrogen (M. Bonitz, J. Vorberger, et al., Phys. Plasmas 31, 110501 (2024)) reproduces the FPOT comparison and identifies the open issues: Kohn-Sham QMD with Kubo-Greenwood cannot practically reach temperatures beyond $T_F$ because converging the required empty bands becomes prohibitive; the high-frequency tail of the conductivity relies on extrapolation, adding uncontrolled error to grouped Rosseland means; and the Kubo-Greenwood formula is an independent-particle approximation that omits electron-electron collision effects beyond the ground-state exchange-correlation functional. Density-functional average-atom models - in which a single ion embedded in its average plasma environment supplies the electronic structure - agree with QMD at low and mid photon energies and extend to high temperature, and multiple-scattering and excited-state opacity methods have been developed, but these approximations themselves require benchmarks. On the experimental side, no measurement of hydrogen opacity exists at any WDM condition. A platform proposed at the National Ignition Facility (J. Lütgert et al., Phys. Plasmas 29, 083301 (2022)) - x-ray absorption radiography of a strongly compressed hydrogen-tritium ice layer at about 150 g/cm^3 and 200 eV, conditions of red-dwarf (low-mass star) interiors, disentangling free-free absorption (inverse bremsstrahlung, in which a photon is absorbed while a free electron accelerates in the Coulomb field of an ion) from Thomson scattering (elastic scattering of photons by free electrons) using two probe photon energies - was accepted as a very first test of free-free opacity models (the plasma Gaunt factor, the quantum correction to classical free-free absorption), with supporting theory still being prepared (S. Schumacher et al., Contrib. Plasma Phys. 2025).

Problem Statement

Complete and validate the first-principles radiative opacity description of dense hydrogen isotopes in the warm dense matter regime. Concretely, for deuterium and deuterium-tritium mixtures over the density-temperature domain traversed by ICF implosions (densities from about 0.5 to several hundred g/cm^3, temperatures from thousands of kelvin through and beyond the electron Fermi temperature):

1. Extend first-principles opacity calculations (quantum molecular dynamics with Kubo-Greenwood linear response, or demonstrably equivalent first-principles methods) to temperatures above the Fermi temperature, where the current table stops, with quantified uncertainties from each source of error: exchange-correlation functional, finite cell size, k-point and empty-band convergence, high-frequency extrapolation of the conductivity, and collisional (electron-electron) corrections to the independent-particle response.

2. Provide a definite, defensible design for the first experimental benchmark of dense hydrogen opacity at a warm dense condition, separating the free-free absorption contribution from Thomson scattering, with in-situ determined density and temperature - as envisioned by the proposed National Ignition Facility red-dwarf-interior platform.

3. Use the validated data to determine, across this regime, where the tabulated opacities used in radiation-hydrodynamics codes (such as the Los Alamos astrophysical opacity table) are wrong and by how much, and propagate the corrected opacities into radiation-transport predictions for ICF target performance and for energy transport in dense-star interiors.

These three deliverables serve one overarching objective - a validated first-principles opacity database for warm dense hydrogen - and all three are required for completion. The problem is complete when a first-principles opacity table exists with stated uncertainties over the full regime, when the opacities tabulated in radiation-hydrodynamics codes have been compared against it with the corrected opacities propagated into radiation-transport predictions, and when a definite, defensible experimental benchmark design at a warm dense condition exists - the executed facility measurement itself remaining a longer-term aspiration that lies beyond the completion criterion.

References

1. S. X. Hu, L. A. Collins, V. N. Goncharov, T. R. Boehly, R. Epstein, R. L. McCrory, S. Skupsky, 'First-principles opacity table of warm dense deuterium for inertial-confinement-fusion applications', Physical Review E 90, 033111 (2014). DOI: 10.1103/PhysRevE.90.033111, https://doi.org/10.1103/PhysRevE.90.033111
2. S. X. Hu, V. N. Goncharov, T. R. Boehly, R. L. McCrory, S. Skupsky, L. A. Collins, J. D. Kress, B. Militzer, 'Impact of first-principles properties of deuterium-tritium on inertial confinement fusion target designs', Physics of Plasmas 22, 056304 (2015). DOI: 10.1063/1.4917477, https://doi.org/10.1063/1.4917477
3. Michael Bonitz, Jan Vorberger, Mandy Bethkenhagen, Maximilian P. Böhme, David M. Ceperley, Alexey Filinov, Thomas Gawne, Frank Graziani, Gianluca Gregori, Paul Hamann, Stephanie B. Hansen, Markus Holzmann, S. X. Hu, Hanno Kählert, Valentin V. Karasiev, Uwe Kleinschmidt, Linda Kordels, Christopher Makait, Burkhard Militzer, Zhandos A. Moldabekov, Carlo Pierleoni, Martin Preising, Kushal Ramakrishna, Ronald Redmer, Sebastian Schwalbe, Pontus Svensson, Tobias Dornheim, 'Toward first principles-based simulations of dense hydrogen' (arXiv title: 'First principles simulations of dense hydrogen'), Physics of Plasmas 31, 110501 (2024). DOI: 10.1063/5.0219405, arXiv:2405.10627, https://doi.org/10.1063/5.0219405
4. J. Lütgert, M. Bethkenhagen, B. Bachmann, L. Divol, D. O. Gericke, S. H. Glenzer, G. N. Hall, N. Izumi, S. F. Khan, O. L. Landen, S. A. MacLaren, L. Masse, R. Redmer, M. Schörner, M. O. Schölmerich, S. Schumacher, N. R. Shaffer, C. E. Starrett, P. A. Sterne, C. Trosseille, T. Döppner, D. Kraus, 'Platform for probing radiation transport properties of hydrogen at conditions found in the deep interiors of red dwarfs', Physics of Plasmas 29, 083301 (2022). DOI: 10.1063/5.0094579, arXiv:2301.08610, https://doi.org/10.1063/5.0094579
5. Samuel Schumacher, Julian Lütgert, Dirk O. Gericke, Mandy Bethkenhagen, Tilo Döppner, Otto L. Landen, Nathaniel R. Shaffer, Charles E. Starrett, Dominik Kraus, 'Modelling Thomson Scattering in a Hydrogen Plasma at Stellar Interior Conditions Using the Hypernetted-Chain Approach', Contributions to Plasma Physics 65, e70002 (2025). DOI: 10.1002/ctpp.70002, https://doi.org/10.1002/ctpp.70002
6. S. X. Hu, L. A. Collins, V. N. Goncharov, J. D. Kress, T. R. Boehly, R. Epstein, R. L. McCrory, S. Skupsky, 'First-principles studies on the equation of state, thermal conductivity, and opacity of deuterium-tritium (DT) and polystyrene (CH) for inertial confinement fusion applications', Journal of Physics: Conference Series 717, 012064 (2016). DOI: 10.1088/1742-6596/717/1/012064, https://doi.org/10.1088/1742-6596/717/1/012064
7. Nathaniel R. Shaffer, Charles E. Starrett, 'Dense plasma opacity via the multiple-scattering method', Physical Review E 105, 015203 (2022). DOI: 10.1103/PhysRevE.105.015203, https://doi.org/10.1103/PhysRevE.105.015203
8. C. E. Starrett, C. J. Fontes, H. B. Tran Tan, J. M. Kasper, J. R. White, 'Dense plasma opacity from excited states method', Physical Review E 110, 045208 (2024). DOI: 10.1103/PhysRevE.110.045208, https://doi.org/10.1103/PhysRevE.110.045208
9. Mikael Tacu, Jean-Christophe Pain, Matthias Pautard, Christophe Blancard, 'Electrical conductivities and low frequency opacities in the warm dense matter regime', arXiv:2504.02368 (2025), https://arxiv.org/abs/2504.02368
10. Shen Zhang, Hongwei Wang, Wei Kang, Ping Zhang, X. T. He, 'Extended application of Kohn-Sham first-principles molecular dynamics method with plane wave approximation at high energy - From cold materials to hot dense plasmas', Physics of Plasmas 23, 042707 (2016). DOI: 10.1063/1.4947212, https://doi.org/10.1063/1.4947212