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.

Yuri L'vovich Klimontovich (28.09.1924--26.10.2002) was an outstanding theoretical physicist who made major contributions to kinetic theory. On the occasion of his 100th birthday we recall his main scientific achievements.

Correlated classical and quantum many-particle systems out of equilibrium are of high interest in many fields, including dense plasmas, correlated solids, and ultracold atoms. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. While for classical systems, in principle, exact simulations are possible via molecular dynamics, this is not the case for quantum systems. Alternatively, one can use many-particle approaches such as hydrodynamics, kinetic theory or nonequilibrium Green functions (NEGF). However, NEGF exhibit a very unfavorable cubic scaling of the CPU time with the number of time steps. An alternative is the G1--G2 scheme [N. Schl\"unzen et al., Phys. Rev. Lett. \textbf{124}, 076601 (2020)] which allows for NEGF simulations with time linear scaling, however, at the cost of large memory consumption. The reason is the need to store the two-particle correlation function. This problem can be overcome for a number of approximations by reformulating the kinetic equations in terms of fluctuations -- an approach that was developed, for classical systems, by Yu.L. Klimontovich [JETP \textbf{33}, 982 (1957)]. Here we present an overview of his ideas and extend them to quantum systems. In particular, we demonstrate that this quantum fluctuations approach can reproduce the nonequilibrium $GW$ approximation [E. Schroedter \textit{et al.}, Cond. Matt. Phys. \textbf{25}, 23401 (2022)] promising high accuracy at low computational cost which arises from an effective semiclassical stochastic sampling procedure. We also demonstrate how to extend the approach to the two-time exchange-correlation functions and the density response properties. [E. Schroedter \textit{et al.}, Phys. Rev. B \textbf{108}, 205109 (2023)].

Correlated quantum many-particle systems out of equilibrium are of high interest in many fields, including correlated solids, ultracold atoms or dense plasmas. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. We have recently presented a quantum fluctuations approach which is equivalent to the nonequilibrium $GW$ approximation [E. Schroedter \textit{et al.}, Cond. Matt. Phys. \textbf{25}, 23401 (2022)] that promises high accuracy at low computational cost. In a second publication [E. Schroedter \textit{et al.}, Phys. Rev. B \textbf{108}, 205109 (2023)], this approach was extended to the two-time exchange-correlation functions and the density response properties. Here, we analyze the properties of this approach in more detail. We demonstrate that the method is equivalent to the Bethe--Salpeter equation for the two-time exchange-correlation function when the generalized Kadanoff-Baym ansatz with Hartree-Fock propagators is applied.

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.

The quantum dynamics of correlated fermionic or bosonic many-body systems following external excitation can be successfully studied using nonequilibrium Green functions (NEGF) or reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle selfenergy or decoupling of the BBGKY-hierarchy, respectively. These approximations are based on Feynman's diagram approaches or on cluster expansions into single-particle and correlation operators. In a recent paper [E. Schroedter, J.-P. Joost, and M. Bonitz, Cond. Matt. Phys. \textbf{25}, 23401 (2022)] we have presented a different approach where, instead of equations of motion for the many-particle NEGF (or density operators), equations for the correlation functions of fluctuations are analyzed. In particular, we derived the stochastic GW and polarization approximations that are closely related to the nonequilibrium GW approximation. Here, we extend this approach to the computation of two-time observables depending on the specific ordering of the underlying operators. In particular, we apply this extension to the calculation of the density correlation function and dynamic structure factor of correlated Hubbard clusters in and out of equilbrium.

In a recent Letter [T. Dornheim \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 255001 (2018)], it was predicted on the basis of \textit{ab initio} quantum Monte Carlo simulations that, in a uniform electron gas, the peak $\omega_0$ of the dynamic structure factor $S(q,\omega)$ exhibits an unusual non-monotonic wave number dependence, where $d\omega_0/dq < 0$, at intermediate $q$, under strong coupling conditions. This effect was subsequently explained by the pair alignment of electrons %at an intermediate range of wave numbers [T. Dornheim \textit{et al.}, Comm. Phys. \textbf{5}, 304 (2022)]. Here we predict that this non-monotonic dispersion resembling the roton-type behavior known from superfluids should be observable in a dense, partially ionized hydrogen plasma. Based on a combination of path integral Monte Carlo simulations and linear response results for the density response function, we present the approximate range of densities, temperatures and wave numbers and make predictions for possible experimental observations.

We explore the non-equilibrium dynamics of a one-dimensional Fermi-Hubbard system as a sensitive testbed for the capabilities of the time-dependent two-particle reduced density matrix (TD2RDM) theory to accurately describe time-dependent correlated systems. We follow the time evolution of the out-of-equilibrium finite-size Fermi-Hubbard model initialized by a quench over extended periods of time. By comparison with exact calculations for small systems and with matrix product state (MPS) calculations for larger systems but limited to short times, we demonstrate that the TD2RDM theory can accurately account for the non-equilibrium dynamics in the regime from weak to moderately strong inter-particle correlations. We find that the quality of the approximate reconstruction of the three-particle cumulant (or correlation) required for the closure of the equations of motion for the reduced density matrix is key to the accuracy of the numerical TD2RDM results. We identify the size of the dynamically induced three-particle correlations and the amplitude of cross correlations between the two- and three-particle cumulants as critical parameters that control the accuracy of the TD2RDM theory when current state-of-the art reconstruction functionals are employed., Comment: 18 pages, 13 figures

Magnetized plasmas are well known to exhibit a rich spectrum of collective modes. Here, we focus on the density modes in dense or cold plasmas, where strong coupling effects alter the mode spectrum known from traditional weakly coupled plasmas. In particular, we study the dynamic structure factor (DSF) of the magnetized one-component plasma with molecular dynamics simulations. Extending our previous results [H.~K\"ahlert and M.~Bonitz, Phys. Rev. Research \textbf{2022}, 4, 013197], it is shown that Bernstein modes can be observed in the weakly magnetized regime, where they are found below the upper hybrid frequency, provided the coupling strength is sufficiently low. We investigate the DSF for a variety of different wave numbers and plasma parameters and show that even small magnetization can give rise to a strong zero-frequency mode perpendicular to the magnetic field and change the dispersion as well as the damping of the upper hybrid mode.

Matter at extreme temperatures and pressures -- commonly known as warm dense matter (WDM) in the literature -- is ubiquitous throughout our Universe and occurs in a number of astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion, and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are routinely probed in x-ray Thomson scattering (XRTS) experiments, and, in addition, are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear-response theory as well as nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques including ab initio path integral Monte Carlo (PIMC) simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrarily complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.

Real-time nonequilibrium Green functions (NEGF) have been very successful to simulate the dynamics of correlated many-particle systems far from equilibrium. However, NEGF simulations are computationally expensive since the effort scales cubically with the simulation duration. Recently we have introduced the G1--G2 scheme that allows for a dramatic reduction to time-linear scaling [Schl\"unzen, Phys. Rev. Lett. 124, 076601 (2020); Joost et al., Phys. Rev. B 101, 245101 (2020)]. Here we tackle another problem: the rapid growth of the computational effort with the system size. In many situations where the system of interest is coupled to a bath, to electric contacts or similar macroscopic systems for which a microscopic resolution of the electronic properties is not necessary, efficient simplifications are possible. This is achieved by the introduction of an embedding selfenergy -- a concept that has been successful in standard NEGF simulations. Here, we demonstrate how the embedding concept can be introduced into the G1--G2 scheme, allowing us to drastically accelerate NEGF embedding simulations. The approach is compatible with all advanced selfenergies that can be represented by the G1--G2 scheme [as described in Joost et al., Phys. Rev. B 105, 165155 (2022)] and retains the memory-less structure of the equations and their time linear scaling. As a numerical illustration we investigate the charge transfer between a Hubbard nanocluster and an additional site which is of relevance for the neutralization of ions in matter.

The quantum dynamics of fermionic or bosonic many-body systems following external excitation can be successfully studied using two-time nonequilibrium Green's functions (NEGF) or single-time reduced density matrix methods. Approximations are introduced via a proper choice of the many-particle self-energy or decoupling of the BBGKY hierarchy. These approximations are based on Feynman's diagram approaches or on cluster expansions into single-particle and correlation operators. Here, we develop a different approach where, instead of equations of motion for the many-particle NEGF (or density operators), single-time equations for the correlation functions of fluctuations are analyzed. We present a derivation of the first two equations of the alternative hierarchy of fluctuations and discuss possible decoupling approximations. In particular, we derive the polarization approximation (PA) which is shown to be equivalent to the single-time version [following by applying the generalized Kadanoff-Baym ansatz (GKBA)] of the nonequilibrium $GW$ approximation with exchange effects of NEGF theory, for weak coupling. The main advantage of the quantum fluctuations approach is that the standard ensemble average can be replaced by a semiclassical average over different initial realizations, as was demonstrated before by Lacroix and co-workers [see e.g. D. Lacroix et al., Phys. Rev. B, 2014, 90, 125112]. Here, we introduce the stochastic $GW$ (SGW) approximation and the stochastic polarization approximation (SPA) which are demonstrated to be equivalent to the single-time $GW$ approximation without and with exchange, respectively, in the weak coupling limit. Our numerical tests confirm that our approach has the same favorable linear scaling with the computation time as the recently developed G1-G2 scheme [Schluenzen et al., Phys. Rev. Lett., 2020, 124, 076601]., Comment: 36 pages, 14 figures

The study of matter under extreme densities and temperatures as they occur e.g. in astrophysical objects and nuclear fusion applications has emerged as one of the most active frontiers in physics, material science, and related disciplines. In this context, a key quantity is given by the dynamic structure factor $S(\mathbf{q},\omega)$, which is probed in scattering experiments -- the most widely used method of diagnostics at these extreme conditions. In addition to its crucial importance for the study of warm dense matter, the modelling of such dynamic properties of correlated quantum many-body systems constitutes one of the most fundamental theoretical challenges of our time. Here we report a hitherto unexplained \emph{roton feature} in $S(\mathbf{q},\omega)$ of the warm dense electron gas, and introduce a microscopic explanation in terms of a new \emph{electronic pair alignment} model. This new paradigm will be highly important for the understanding of warm dense matter, and has a direct impact on the interpretation of scattering experiments. Moreover, we expect our results to give unprecedented insights into the dynamics of a number of correlated quantum many-body systems such as ultracold helium, dipolar supersolids, and bilayer heterostructures.

Dynamical screening is a key property of charged many-particle systems. Its theoretical description is based on the $GW$ approximation that is extensively applied for ground-state and equilibrium situations but also for systems driven out of equilibrium. The main limitation of the $GW$ approximation is the neglect of strong electronic correlation effects that are important in many materials as well as in dense plasmas. Here we derive the dynamically screened ladder (DSL) approximation that selfconsistently includes, in addition to the $GW$ diagrams, also particle--particle and particle--hole $T$-matrix diagrams. The derivation is based on reduced-density-operator theory and the result is equivalent to the recently presented G1--G2 scheme [Schl\"unzen \textit{et al.}, Phys. Rev. Lett. \textbf{124}, 076601 (2020); Joost \textit{et al.}, Phys. Rev. B \textbf{101}, 245101 (2020)]. We perform extensive time-dependent DSL simulations for finite Hubbard clusters and present tests against exact results that confirm excellent accuracy as well as total energy conservation of the approximation. At strong coupling and for long simulation durations, instabilities are observed. These problems are solved by enforcing contraction consistency and applying a purification approach.

Plasmas in strong magnetic fields have been mainly studied in two distinct limiting cases--that of weak and strong nonideality with very different physical properties. While the former is well described by the familiar theory of Braginskii, the latter regime is closer to the behavior of a Coulomb liquid. Here we study in detail the transition between both regimes. We focus on the evolution of the dynamic structure factor of the magnetized one-component plasma from weak to strong coupling, which is studied with first-principle molecular dynamics simulations. The simulations show the vanishing of Bernstein modes and the emergence of higher harmonics of the upper hybrid mode across the magnetic field, a redistribution of spectral power between the two main collective modes under oblique angles, and a suppression of plasmon damping along the magnetic field. Comparison with results from various models, including the random phase approximation, a Mermin-type dielectric function, and the Quasi-Localized Charge Approximation show that none of the theories is capable of reproducing the crossover that occurs when the coupling parameter is on the order of unity. The findings are relevant to the scattering spectra, stopping power, and transport coefficients of correlated magnetized plasmas.

We present extensive new ab initio path integral Monte Carlo (PIMC) results for an electron gas at warm dense matter conditions that is subject to multiple harmonic perturbations. In addition to the previously investigated nonlinear effects at the original wave number [Dornheim \emph{et al.}, PRL \textbf{125}, 085001 (2020)] and the excitation of higher harmonics [Dornheim \emph{et al.}, PRR \textbf{3}, 033231 (2021)], the presence of multiple external potentials leads to mode-coupling effects, which constitute the dominant nonlinear effect and lead to a substantially more complicated density response compared to linear response theory. One possibility to estimate mode-coupling effects from a PIMC simulation of the unperturbed system is given in terms of generalized imaginary-time correlation functions that have been recently introduced by Dornheim \emph{et al.}~[JCP \textbf{155}, 054110 (2021)]. In addition, we extend our previous analytical theory of the nonlinear density response of the electron gas in terms of the static local field correction [Dornheim \emph{et al.}, PRL \textbf{125}, 235001 (2020)], which allows for a highly accurate description of the PIMC results with negligible computational cost.

In a recent Letter [Balzer \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 267602 (2018)] it was demonstrated that ions impacting a correlated graphene cluster can excite strongly nonequilibrium states. In particular, this can lead to an enhanced population of bound pairs of electrons with opposite spin -- doublons -- where the doublon number can be increased via multiple ion impacts. These predictions were made based on nonequilibrium Green functions (NEGF) simulations allowing for a time-dependent non-perturbative study of the energy loss of charged particles penetrating a strongly correlated system. Here we extend these simulations to larger clusters and longer simulation times, utilizing the recently developed G1--G2 scheme [Sch\"unzen \textit{et al.}, Phys. Rev. Lett. \textbf{124}, 076601 (2020)] which allows for a dramatic speedup of NEGF simulations. Furthermore, we investigate the dependence of the energy and doublon number on the time interval between ion impacts and on the impact point.

In a recent Letter, Dornheim et al. [PRL 125, 085001 (2020)] have investigated the nonlinear density response of the uniform electron gas in the warm dense matter regime. More specifically, they have studied the cubic response function at the first harmonic, which cannot be neglected in many situations of experimental relevance. In this work, we go one step further and study the full spectrum of excitations at the higher harmonics of the original perturbation based on extensive new ab initio path integral Monte Carlo (PIMC) simulations. We find that the dominant contribution to the density response beyond linear response theory is given by the quadratic response function at the second harmonic in the moderately nonlinear regime. Furthermore, we show that the nonlinear density response is highly sensitive to exchange-correlation effects, which makes it a potentially valuable new tool of diagnostics. To this end, we present a new theoretical description of the nonlinear electronic density response based on the recent effective static approximation to the local field correction [PRL 125, 235001 (2020)], which accurately reproduces our PIMC data with negligible computational cost.

We study the time-dependent neutralization of a slow highly charged ion that penetrates a hexagonal hollow-centred graphene nanoflake. To compute the ultrafast charge transfer dynamics, we apply an effective Hubbard nanocluster model and use the method of nonequilibrium Green functions (NEGF) in conjunction with an embedding self-energy scheme which allows one to follow the temporal changes of the number of electrons in the nanoflake. We perform extensive simulations of the charge transfer dynamics for a broad range of ion charge states and impact velocities. The results are used to put forward a simple semi-analytical model of the neutralization dynamics that is in very good agreement with transmission experiments, in which highly charged xenon ions pass through sheets of single-layer graphene.

In a classical plasma the momentum distribution, $n(k)$, decays exponentially, for large $k$, and the same is observed for an ideal Fermi gas. However, when quantum and correlation effects are relevant simultaneously, an algebraic decay, $n_\infty(k)\sim k^{-8}$ has been predicted. This is of relevance for cross sections and threshold processes in dense plasmas that depend on the number of energetic particles. Here we present extensive \textit{ab initio} results for the momentum distribution of the nonideal uniform electron gas at warm dense matter conditions. Our results are based on first principle fermionic path integral Monte Carlo (CPIMC) simulations and clearly confirm the $k^{-8}$ asymptotic. This asymptotic behavior is directly linked to short-range correlations which are analyzed via the on-top pair distribution function (on-top PDF), i.e. the PDF of electrons with opposite spin. We present extensive results for the density and temperature dependence of the on-top PDF and for the momentum distribution in the entire momentum range.

The screening of a test charge by partially degenerate non-ideal free electrons at conditions related to warm dense matter and dense plasmas is investigated using linear response theory and the local field correction based on ab inito Quantum Monte-Carlo simulations data. The analysis of the obtained results is performed by comparing to the random phase approximation and the Singwi-Tosi-Land-Sj\"olander approximation. The applicability of the long-wavelength approximation for the description of screening is investigated. The impact of electronic exchange-correlations effects on structural properties and the applicability of the screened potential from linear response theory for the simulation of the dynamics of ions are discussed.

Warm dense matter (WDM) is an exotic state on the border between condensed matter and dense plasmas. Important occurrences of WDM include dense astrophysical objects, matter in the core of our Earth, as well as matter produced in strong compression experiments. As of late, x-ray Thomson scattering has become an advanced tool to diagnose WDM. The interpretation of the data requires model input for the dynamic structure factor $S(q,\omega)$ and the plasmon dispersion $\omega(q)$. Recently the first \textit{ab initio} results for $S(q,\omega)$ of the homogeneous warm dense electron gas were obtained from path integral Monte Carlo simulations, [Dornheim \textit{et al.}, Phys. Rev. Lett. \textbf{121}, 255001 (2018)]. Here, we analyse the effects of correlations and finite temperature on the dynamic dielectric function and the plasmon dispersion. Our results for the plasmon dispersion and damping differ significantly from the random phase approximation and from earlier models of the correlated electron gas. Moreover, we show when commonly used weak damping approximations break down and how the method of complex zeros of the dielectric function can solve this problem for WDM conditions.

There is growing interest in warm dense matter (WDM) -- an exotic state on the border between condensed matter and plasmas. Due to the simultaneous importance of quantum and correlation effects WDM is complicated to treat theoretically. A key role has been played by \textit{ab initio} path integral Monte Carlo (PIMC) simulations, and recently extensive results for thermodynamic quantities have been obtained. The first extension of PIMC simulations to the dynamic structure factor of the uniform electron gas were reported by Dornheim \textit{et al.} [Phys. Rev. Lett. \textbf{121}, 255001 (2018)]. This was based on an accurate reconstruction of the dynamic local field correction. Here we extend this concept to other dynamical quantities of the warm dense electron gas including the dynamic susceptibility, the dielectric function and the conductivity.

Warm dense matter (WDM)---an extreme state with high temperatures and densities that occurs e.g. in astrophysical objects---constitutes one of the most active fields in plasma physics and materials science. These conditions can be realized in the lab by shock compression or laser excitation, and the most accurate experimental diagnostics is achieved with lasers and free electron lasers which is theoretically modeled using linear response theory. Here, we present first \textit{ab initio} path integral Monte Carlo results for the nonlinear density response of correlated electrons in WDM and show that for many situations of experimental relevance nonlinear effects cannot be neglected.

Finite graphene nanoribbon (GNR) heterostructures host intriguing topological in-gap states (Rizzo, D. J. et al.~\textit{Nature} \textbf{2018}, \textit{560}, 204]). These states may be localized either at the bulk edges, or at the ends of the structure. Here we show that correlation effects (not included in previous density functional simulations) play a key role in these systems: they result in increased magnetic moments at the ribbon edges accompanied by a significant energy renormalization of the topological end states -- even in the presence of a metallic substrate. Our computed results are in excellent agreement with the experiments. Furthermore, we discover a striking, novel mechanism that causes an energy splitting of the non-zero-energy topological end states for a weakly screened system. We predict that similar effects should be observable in other GNR heterostructures as well.

The dynamics of strongly correlated fermions following an external excitation reveals extremely rich collective quantum effects. Examples are fermionic atoms in optical lattices, electrons in correlated materials, and dense quantum plasmas. Presently, the only quantum-dynamics approach that rigorously describes these processes in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their $T^3$ scaling with the simulation duration $T$. Recently, $T^2$ scaling was achieved with the generalized Kadanoff--Baym ansatz (GKBA) which has substantially extended the scope of NEGF simulations. Here we present a novel approach to GKBA-NEGF simulations that is of order $T$, and demonstrate its remarkable capabilities.

The study of matter at extreme densities and temperatures as they occur in astrophysical objects and state-of-the art experiments with high-intensity lasers is of high current interest for many applications. While no overarching theory for this regime exists, accurate data for the density response of correlated electrons to an external perturbation are of paramount importance. In this context, the key quantity is given by the local field correction (LFC), which provides a wave-vector resolved description of exchange-correlation effects. In this work, we present extensive new path integral Monte Carlo (PIMC) results for the static LFC of the uniform electron gas, which are subsequently used to train a fully connected deep neural network. This allows us to present a continuous representation of the LFC with respect to wave-vector, density, and temperature covering the entire warm dense matter regime. Both the PIMC data and neural-net results are available online. Moreover, we expect the presented combination of ab initio calculations with machine-learning methods to be a promising strategy for many applications.

Being motivated by the surge of fermionic quantum Monte Carlo simulations at finite temperature, we present a detailed analysis of the permutation-cycle properties of path integral Monte Carlo (PIMC) simulations of degenerate electrons. Particular emphasis is put onto the uniform electron gas in the warm dense matter regime. We carry out PIMC simulations of up to $N=100$ electrons and investigate exchange-cycle frequencies, which are found not to follow any simple exponential law even in the case of ideal fermions due to the finite size of the simulation box. Moreover, we introduce a permutation-cycle correlation function, which allows us to analyse the joint probability to simultaneously find cycles of different lengths within a single configuration. Again, we find that finite-size effects predominate the observed behaviour. Finally, we briefly consider an inhomogeneous system, namely electrons in a $2D$ harmonic trap. We expect our results to be of interest for the further development of fermionic PIMC methods, in particular to alleviate the notorious fermion sign problem.

Leonid Keldysh -- one of the most influential theoretical physicists of the 20th century -- passed away in November 2016. Keldysh is best known for the diagrammatic formulation of real-time (nonequilibrium) Green functions theory and for the theory of strong field ionization of atoms. Both theories profoundly changed large areas of theoretical physics and stimulated important experiments. Both these discoveries emerged almost simultaneously -- like Einstein, also Keldysh had his \textit{annus mirabilis} -- the year 1964. But the list of his theoretical developments is much broader and is briefly reviewed here.

The self-diffusion phenomenon in a two-dimensional dusty plasma at extremely strong (effective) magnetic fields is studied experimentally and by means of molecular dynamics simulations. In the experiment the high magnetic field is introduced by rotating the particle cloud and observing the particle trajectories in a co-rotating frame, which allows reaching effective magnetic fields up to 3000 Tesla. The experimental results confirm the predictions of the simulations: (i) super-diffusive behavior is found at intermediate time-scales and (ii) the dependence of the self-diffusion coefficient on the magnetic field is well reproduced., Comment: accepted by Physical Review E

We review the uniform electron gas (UEG) at finite temperature and over a broad density range relevant for warm dense matter (WDM) applications. We provide an overview of different simulation techniques, focusing on recent developments in the dielectric formalism and quantum Monte Carlo (QMC). Our primary focus is on two novel QMC methods: Permutation blocking path integral MC (PB-PIMC) and configuration PIMC (CPIMC). In fact, a combination of PB-PIMC and CPIMC has allowed for a highly accurate description of the warm dense UEG over a broad density-temperature range. We are able to effectively avoid the notorious fermion sign problem, without invoking uncontrolled approximations such as the fixed node approximation. Furthermore, a new finite-size correction scheme is presented that makes it possible to treat the UEG in the thermodynamic limit without loss of accuracy. In addition, we in detail discuss the construction of a parametrization of the exchange-correlation free energy that provides a complete description of the UEG and is of crucial importance as input for the simulation of real WDM applications. Further, we test previous theories, including restricted PIMC, finite-temperature Green functions, the classical mapping by Perrot and Dharma-wardana, and various dielectric methods such as the random phase approximation, or the Singwi-Tosi-Land-Sj\"olander, Vashishta-Singwi and the recent Tanaka scheme for the local field correction. Thus, for the first time, thorough benchmarks of important approximation schemes regarding various quantities such as different energies, in particular the exchange-correlation free energy, and the static structure factor, are possible. Finally, we outline a way how to rigorously extend our QMC studies to the inhomogeneous electron gas and present first ab initio data for the static density response and for the static local field correction.

Strongly correlated systems of fermions have a number of exciting collective properties. Among them, the creation of a lattice that is occupied by doublons, i.e. two quantum particles with opposite spins, offers interesting electronic properties. In the past a variety of methods have been proposed to control doublon formation, both, spatially and temporally. Here, a novel mechanism is proposed and verified by exact diagonalization and nonequilibrium Green functions simulations---fermionic doublon creation by the impact of energetic ions. We report the formation of a nonequilibrium steady state with homogeneous doublon distribution. The effect should be observable in strongly correlated solids in contact with a high-pressure plasma and in fermionic atoms in optical lattices.

The uniform electron gas at finite temperature is of high current interest for warm dense matter research. The complicated interplay of quantum degeneracy and Coulomb coupling effects is fully contained in the pair distribution function or, equivalently, the static strucutre factor. By combining exact quantum Monte Carlo results for large wave vectors with the long-range behavior from the Singwi-Tosi-Land-Sj\"olander approximation, we are able to obtain highly accurate data for the static structure factor over the entire $k$-range. This allows us to gauge the accuracy of previous approximations and discuss their respective shortcomings. Further, our new data will serve as valuable input for the computation of other quantities.

Precise knowledge of the static density response function (SDRF) of the uniform electron gas (UEG) serves as key input for numerous applications, most importantly for density functional theory beyond generalized gradient approximations. Here we extend the configuration path integral Monte Carlo (CPIMC) formalism that was previously applied to the spatially uniform electron gas to the case of an inhomogeneous electron gas by adding a spatially periodic external potential. This procedure has recently been successfully used in permutation blocking path integral Monte Carlo simulations (PB-PIMC) of the warm dense electron gas [Dornheim \textit{et al.}, Phys. Rev. E in press, arXiv:1706.00315], but this method is restricted to low and moderate densities. Implementing this procedure into CPIMC allows us to obtain exact finite temperature results for the SDRF of the electron gas at \textit{high to moderate densities} closing the gap left open by the PB-PIMC data. In this paper we demonstrate how the CPIMC formalism can be efficiently extended to the spatially inhomogeneous electron gas and present the first data points. Finally, we discuss finite size errors involved in the quantum Monte Carlo results for the SDRF in detail and present a solution how to remove them that is based on a generalization of ground state techniques.

The static density response of the uniform electron gas is of fundamental importance for numerous applications. Here, we employ the recently developed \textit{ab initio} permutation blocking path integral Monte Carlo (PB-PIMC) technique [T.~Dornheim \textit{et al.}, \textit{New J.~Phys.}~\textbf{17}, 073017 (2015)] to carry out extensive simulations of the harmonically perturbed electron gas at warm dense matter conditions. In particular, we investigate in detail the validity of linear response theory and demonstrate that PB-PIMC allows to obtain highly accurate results for the static density response function and, thus, the static local field correction. A comparison with dielectric approximations to our new \textit{ab initio} data reveals the need for an exact treatment of correlations. Finally, we consider a superposition of multiple perturbations and discuss the implications for the calculation of the static response function.

In a recent Letter [T.~Dornheim \textit{et al.}, Phys. Rev. Lett. \textbf{117}, 156403 (2016)], we presented the first \textit{ab initio} quantum Monte-Carlo (QMC) results of the warm dense electron gas in the thermodynamic limit. However, a complete parametrization of the exchange-correlation free energy with respect to density, temperature, and spin polarization remained out of reach due to the absence of (i) accurate QMC results below $\theta=k_\text{B}T/E_\text{F}=0.5$ and (ii) of QMC results for spin polarizations different from the paramagnetic case. Here we overcome both remaining limitations. By closing the gap to the ground state and by performing extensive QMC simulations for different spin polarizations, we are able to obtain the first complete \textit{ab initio} exchange-correlation free energy functional; the accuracy achieved is an unprecedented $\sim 0.3\%$. This also allows us to quantify the accuracy and systematic errors of various previous approximate functionals.

This work examines the formation of wake fields caused by ions streaming around a charged dust particle, using three-dimensional particle-in-cell (PIC) simulations with charge-neutral collisions included. The influence of an external driving electric field, which leads to a non-Maxwellian distribution of ions, is investigated in detail. The wake features formed for non-Maxwellian ions exhibit significant deviations from those observed within the model of a shifted Maxwellian distribution. The dependence of the peak amplitude and position of the wake potential upon the degree of collisionality is analyzed for a wide range of streaming velocities (Mach numbers). In contrast to a shifted Maxwellian distribution of ions, the drift-driven non-Maxwellian distribution exhibits an increase of the wake amplitude of the first attractive peak with increase in collisionality for high streaming velocities. At very low Mach numbers, collision-induced amplification is observed for Maxwellian as well as non-Maxwellian distributions.

Time reversal symmetry is a fundamental property of many quantum mechanical systems. The relation between statistical physics and time reversal is subtle and not all statistical theories conserve this particular symmetry, most notably hydrodynamic equations and kinetic equations such as the Boltzmann equation. In this article it is shown analytically that quantum kinetic generalizations of the Boltzmann equation that are derived using the nonequilibrium Green functions formalism as well as all approximations that stem from $\Phi$-derivable selfenergies are time reversal invariant.

The uniform electron gas is a key model system in the description of matter, including dense plasmas and solid state systems. However, the simultaneous occurence of quantum, correlation, and thermal effects makes the theoretical description challenging. For these reasons, over the last half century many analytical approaches have been developed the accuracy of which has remained unclear. We have recently obtained the first \textit{ab initio} data for the exchange correlation free energy of the uniform electron gas [T. Dornheim \textit{et al.}, Phys.~Rev.~Lett.~\textbf{117}, 156403 (2016)] which now provides the opportunity to assess the quality of the mentioned approaches and parametrizations. Particular emphasis is put on the warm dense matter regime, where we find significant discrepancies between the different approaches.

Warm dense matter is one of the most active frontiers in plasma physics due to its relevance for dense astrophysical objects as well as for novel laboratory experiments in which matter is being strongly compressed e.g. by high-power lasers. Its description is theoretically very challenging as it contains correlated quantum electrons at finite temperature---a system that cannot be accurately modeled by standard analytical or ground state approaches. Recently several breakthroughs have been achieved in the field of fermionic quantum Monte Carlo simulations. First, it was shown that exact simulations of a finite model system ($30 \dots 100$ electrons) is possible that avoid any simplifying approximations such as fixed nodes [Schoof {\em et al.}, Phys. Rev. Lett. {\bf 115}, 130402 (2015)]. Second, a novel way to accurately extrapolate these results to the thermodynamic limit was reported by Dornheim {\em et al.} [Phys. Rev. Lett. {\bf 117}, 156403 (2016)]. As a result, now thermodynamic results for the warm dense electron gas are available that have an unprecedented accuracy on the order of $0.1\%$. Here we present an overview on these results and discuss limitations and future directions.

Comparative studies of capacitively coupled radio-frequency discharges in helium and argon at pressures between 10 and 80 Pa are presented applying two different fluid modelling approaches as well as two independently developed particle-in-cell/Monte Carlo collision (PIC/MCC) codes. The focus is on the analysis of the range of applicability of a recently proposed fluid model including an improved drift-diffusion approximation for the electron component as well as its comparison with fluid modelling results using the classical drift-diffusion approximation and benchmark results obtained by PIC/MCC simulations. Main features of this time- and space-dependent fluid model are given. It is found that the novel approach shows generally quite good agreement with the macroscopic properties derived by the kinetic simulations and is largely able to characterize qualitatively and quantitatively the discharge behaviour even at conditions when the classical fluid modelling approach fails. Furthermore, the excellent agreement between the two PIC/MCC simulation codes using the velocity Verlet method for the integration of the equations of motion verifies their accuracy and applicability.

We perform \emph{ab initio} quantum Monte Carlo (QMC) simulations of the warm dense uniform electron gas in the thermodynamic limit. By combining QMC data with linear response theory we are able to remove finite-size errors from the potential energy over the entire warm dense regime, overcoming the deficiencies of the existing finite-size corrections by Brown \emph{et al.}~[PRL \textbf{110}, 146405 (2013)]. Extensive new QMC results for up to $N=1000$ electrons enable us to compute the potential energy $V$ and the exchange-correlation free energy $F_{xc}$ of the macroscopic electron gas with an unprecedented accuracy of $|\Delta V|/|V|, |\Delta F_{xc}|/|F|_{xc} \sim 10^{-3}$. A comparison of our new data to the recent parametrization of $F_{xc}$ by Karasiev {\em et al.} [PRL {\bf 112}, 076403 (2014)] reveals significant deviations to the latter.

A combined nonequilibrium Green functions-Ehrenfest dynamics approach is developed that allows for a time-dependent study of the energy loss of a charged particle penetrating a strongly correlated system at zero and finite temperature. Numerical results are presented for finite inhomogeneous two-dimensional Fermi-Hubbard models, where the many-electron dynamics in the target are treated fully quantum mechanically and the motion of the projectile is treated classically. The simulations are based on the solution of the two-time Dyson (Keldysh-Kadanoff-Baym) equations using the second-order Born, third-order and T-matrix approximations of the self-energy. As application, we consider protons and helium nuclei with a kinetic energy between 1 and 500 keV/u passing through planar fragments of the two-dimensional honeycomb lattice and, in particular, examine the influence of electron-electron correlations on the energy exchange between projectile and electron system. We investigate the time dependence of the projectile's kinetic energy (stopping power), the electron density, the double occupancy and the photoemission spectrum. Finally, we show that, for a suitable choice of the Hubbard model parameters, the results for the stopping power are in fair agreement with ab-initio simulations for particle irradiation of single-layer graphene., Comment: 13 pages, 12 figures

The uniform electron gas (UEG) at finite temperature is of high current interest due to its key relevance for many applications including dense plasmas and laser excited solids. In particular, density functional theory heavily relies on accurate thermodynamic data for the UEG. Until recently, the only existing first-principle results had been obtained for $N=33$ electrons with restricted path integral Monte Carlo (RPIMC), for low to moderate density, $r_s = \overline{r}/a_B \gtrsim 1$. This data has been complemented by Configuration path integral Monte Carlo (CPIMC) simulations for $r_s \leq 1$ that substantially deviate from RPIMC towards smaller $r_s$ and low temperature. In this work, we present results from an independent third method---the recently developed permutation blocking path integral Monte Carlo (PB-PIMC) approach [T. Dornheim \textit{et al.}, NJP \textbf{17}, 073017 (2015)] which we extend to the UEG. Interestingly, PB-PIMC allows us to perform simulations over the entire density range down to half the Fermi temperature ($\theta=k_BT/E_F=0.5$) and, therefore, to compare our results to both aforementioned methods. While we find excellent agreement with CPIMC, where results are available, we observe deviations from RPIMC that are beyond the statistical errors and increase with density.

The effective potential $\Phi$ of a classical ion in a weakly correlated quantum plasma in thermodynamic equilibrium at finite temperature is well described by the RPA screened Coulomb potential. Additionally, collision effects can be included via a relaxation time ansatz (Mermin dielectric function). These potentials are used to study the quality of various statically screened potentials that were recently proposed by Shukla and Eliasson (SE) [Phys. Rev. Lett. {\bf 108}, 165007 (2012)], Akbari--Moghanjoughi (AM) [Phys. Plasmas {\bf 22}, 022103 (2015)] and Stanton and Murillo (SM) [Phys. Rev. E {\bf 91}, 033104 (2015)] starting from quantum hydrodynamic theory (QHD). Our analysis reveals that the SE potential is qualitatively different from the full potential, whereas the SM potential (at any temperature) and the AM potential (at zero temperature) are significantly more accurate. This confirms the correctness of the recently derived [Michta {\em et al.}, Contrib. Plasma Phys. {\bf 55}, (2015)] pre-factor $1/9$ in front of the Bohm term of QHD for fermions.

We develop a time-dependent theory to investigate electron dynamics and photoionization processes of diatomic molecules interacting with strong laser fields including electron-electron correlation effects. We combine the recently formulated time-dependent generalized-active-space configuration interaction theory [D. Hochstuhl and M. Bonitz, Phys. Rev. A 86, 053424 (2012); S. Bauch, et al., Phys. Rev. A 90, 062508 (2014)] with a prolate spheroidal basis set including localized orbitals and continuum states to describe the bound electrons and the outgoing photoelectron. As an example, we study the strong-field ionization of the two-center four-electron lithium hydride molecule in different intensity regimes. By using single-cycle pulses, two orientations of the asymmetric heteronuclear molecule are investigated: Li-H, with the electrical field pointing from H to Li, and the opposite case of H-Li. The preferred orientation for ionization is determined and we find a transition from H-Li, for low intensity, to Li-H, for high intensity. The influence of electron correlations is studied at different levels of approximation, and we find a significant change in the preferred orientation. For certain intensity regimes, even an interchange of the preferred configuration is observed, relative to the uncorrelated simulations. Further insight is provided by detailed comparisons of photoelectron angular distributions with and without correlation effects taken into account.