First principles simulations of dense hydrogen

Michael Bonitz,Jan Vorberger,Mandy Bethkenhagen,Maximilian Böhme,David Ceperley,Alexey Filinov,Thomas Gawne,Frank Graziani,Gianluca Gregori,Paul Hamann,Stephanie Hansen,Markus Holzmann,S.X. Hu,Hanno Kählert,Valentin Karasiev,Uwe Kleinschmidt,Linda Kordts,Christopher Makait,Burkhard Militzer,Zhandos Moldabekov,Carlo Pierleoni,Martin Preising,Kushal Ramakrishna,Ronald Redmer,Sebastian Schwalbe,Pontus Svensson,Tobias Dornheim

2024-05-17

Abstract:Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data. rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange-correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principle PIMC and DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions.

Computational Physics,Plasma Physics

What problem does this paper attempt to address?

This paper focuses on the properties of hydrogen under high pressure, which is crucial for astrophysics (such as the interiors of planets and stars, brown dwarfs, and compact star atmospheres) as well as laboratory experiments (including inertial confinement fusion). Despite the existence of experimental data on the equation of state, conductivity, and Thomson scattering spectra, the measurement analysis under extreme pressure and temperature conditions often involves additional model assumptions, making it difficult to rigorously assess the accuracy of experimental data. Theoretical and simulation approaches have generated a large amount of data, but each method has its inherent limitations, such as the fermion sign problem in path integral Monte Carlo (PIMC) and the exchange-correlation function approximation in density functional theory (DFT).

The paper discusses recent breakthroughs in PIMC and DFT simulations and uses these results to benchmark different simulation methods. It updates the phase diagram of hydrogen under high pressure, expected phase transitions, and thermodynamic properties, including the equation of state and momentum distribution. Additionally, it discusses the dynamic properties of warm dense hydrogen, such as conductivity, dynamic structure factor, plasma dispersion, imaginary time structure, and density response function. The paper concludes by proposing a strategy to combine different simulations in order to achieve accurate theoretical predictions.

In summary, the problem that this paper attempts to address is how to improve the understanding of the properties of hydrogen under high pressure through theoretical and simulation methods, particularly in terms of phase transitions, thermodynamic properties, and dynamic characteristics, and how to combine different methods to overcome the limitations of existing approaches and provide more reliable and precise theoretical predictions.