Search

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.

The theory of Nonequilibrium Green functions (NEGF) has seen a rapid development over the recent three decades. Applications include diverse correlated many-body systems in and out of equilibrium. Very good agreement with experiments and available exact theoretical results could be demonstrated if the proper selfenergy approximations were used. However, full two-time NEGF simulations are computationally costly, as they suffer from a cubic scaling of the computation time with the simulation duration. Recently we have introduced the G1-G2 scheme that exactly reformulates the Kadanoff-Baym ansatz with Hartree-Fock propagators (HF-GKBA) into time-local equations, allowing for a dramatic reduction of the scaling to time-linear scaling [Schluenzen et al., Phys. Rev. Lett. \textbf{124}, 076601 (2020)]. Remarkably, this scaling is achieved quickly, and also for high-level selfenergies, including the nonequilibrium $GW$ and $T$-matrix approximations [Joost et al., Phys. Rev. B \textbf{101}, 245101 (2020)]. Even the dynamically screened ladder approximation is now feasible [Joost et al., Phys. Rev. B \textbf{105}, 165155 (2022)], and also applications to electron-boson systems were demonstrated. Here we present an overview on recent results that were achieved with the G1--G2 scheme. We discuss problems and open questions and present further ideas how to overcome the current limitations of the scheme.We illustrate the G1--G2 scheme by presenting applying it to the excitation dynamics of Hubbard clusters, to optical excitation of graphene, and to charge transfer during stopping of ions by correlated materials.