Real-time kadanoff-baym approach to plasma oscillations in a correlated electron Gas

Momentum distribution of the uniform electron gas at finite temperature: Effects of spin polarization

We carry out extensive direct path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) at finite temperature for different values of the spin-polarization ξ. This allows us to unambiguously quantify the impact of spin effects on the momentum distribution function n(k) and related properties. We find that interesting physical effects like the interaction-induced increase in the occupation of the zero-momentum state n(0) substantially depend on ξ. Our results further advance the current understanding of the UEG as a fundamental model system, and are of practical relevance for the description of transport properties of warm dense matter in an external magnetic field. All PIMC results are freely available online and can be used as a benchmark for the development of methods and applications.

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green's function methods: Theory and application to GW and GF2

One-particle Green's function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green's functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as ⟨S²⟩ and ⟨N²⟩, can be challenging because they require a solution of the computationally expensive Bethe-Salpeter equation to find two-particle Green's functions. We demonstrate that the solution of the Bethe-Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann-Feynman theorem to self-consistent one-particle Green's function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of ⟨S²⟩ and ⟨N²⟩ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

Real-time GW-BSE investigations on spin-valley exciton dynamics in monolayer transition metal dichalcogenide

We develop an ab initio nonadiabatic molecular dynamics (NAMD) method based on GW plus real-time Bethe-Salpeter equation (GW + rtBSE-NAMD) for the spin-resolved exciton dynamics. From investigations on MoS₂, we provide a comprehensive picture of spin-valley exciton dynamics where the electron-phonon (e-ph) scattering, spin-orbit interaction (SOI), and electron-hole (e-h) interactions come into play collectively. In particular, we provide a direct evidence that e-h exchange interaction plays a dominant role in the fast valley depolarization within a few picoseconds, which is in excellent agreement with experiments. Moreover, there are bright-to-dark exciton transitions induced by e-ph scattering and SOI. Our study proves that e-h many-body effects are essential to understand the spin-valley exciton dynamics in transition metal dichalcogenides and the newly developed GW + rtBSE-NAMD method provides a powerful tool for exciton dynamics in extended systems with time, space, momentum, energy, and spin resolution.

Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte Carlo simulations

We systematically investigate finite-size effects in the dynamic structure factor S(q,ω) of the uniform electron gas obtained via the analytic continuation of ab initio path integral Monte Carlo data for the imaginary-time density-density correlation function F(q,τ). Using the recent scheme by Dornheim et al. [Phys. Rev. Lett. 121, 255001 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.255001], we find that the reconstructed spectra are not afflicted with any finite-size effects for as few as N=14 electrons both at warm dense matter (WDM) conditions and at the margins of the strongly correlated electron liquid regime. Our results further corroborate the high quality of our current description of the dynamic density response of correlated electrons, which is of high importance for many applications in WDM theory and beyond.

Nonlinear Electronic Density Response in Warm Dense Matter

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 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.

Ultrafast dynamics of strongly correlated fermions-nonequilibrium Green functions and selfenergy approximations

This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). We discuss applications of NEGF simulations to describe the femtosecond dynamics of various finite fermionic systems following an excitation out of equilibrium. This includes the expansion dynamics of ultracold atoms in optical lattices following a confinement quench and the excitation of strongly correlated electrons in a solid by the impact of a charged particle. NEGF, presently, are the only ab initio quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the second-order Born approximation (SOA). While they correctly capture the qualitative trends of the relaxation towards equilibrium, the reliability and accuracy of these NEGF simulations has remained open, for a long time. Here we report on recent tests of NEGF simulations for finite lattice systems against exact-diagonalization and density-matrix-renormalization-group benchmark data. The results confirm the high accuracy and predictive capability of NEGF simulations-provided selfenergies are used that go beyond the SOA and adequately include strong correlation and dynamical-screening effects. With an extended arsenal of selfenergies that can be used effectively, the NEGF approach has the potential of becoming a powerful simulation tool with broad areas of new applications including strongly correlated solids and ultracold atoms. The present review aims at making such applications possible. To this end we present a selfcontained introduction to the theory of NEGF and give an overview on recent numerical applications to compute the ultrafast relaxation dynamics of correlated fermions. In the second part we give a detailed introduction to selfenergies beyond the SOA. Important examples are the third-order approximation, the [Formula: see text] approximation, the T-matrix approximation and the fluctuating-exchange approximation. We give a comprehensive summary of the explicit selfenergy expressions for a variety of systems of practical relevance, starting from the most general expressions (general basis) and the Feynman diagrams, and including also the important cases of diagonal basis sets, the Hubbard model and the differences occuring for bosons and fermions. With these details, and information on the computational effort and scaling with the basis size and propagation duration, readers will be able to choose the proper basis set and straightforwardly implement and apply advanced selfenergy approximations to a broad class of systems.

Achieving the Scaling Limit for Nonequilibrium Green Functions Simulations

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 the nonequilibrium Green functions (NEGF) method. 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), for second-order Born (SOA) selfenergies, which has substantially extended the scope of NEGF simulations. Here we demonstrate that GKBA-NEGF simulations can be performed with order T^{1} scaling, both for SOA and GW selfenergies, and point out the remarkable capabilities of this approach.

Ab Initio Lifetime and Concomitant Double-Excitation Character of Plasmons at Metallic Densities

The accurate calculation of excited state properties of interacting electrons in the condensed phase is an immense challenge in computational physics. Here, we use state-of-the-art equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) to calculate the dynamic structure factor, which can be experimentally measured by inelastic x-ray and electron scattering. Our calculations are performed on the uniform electron gas at densities corresponding to Wigner-Seitz radii of r_{s}=5, 4, and 3 corresponding to the valence electron densities of common metals. We compare our results to those obtained using the random-phase approximation (RPA), which is known to provide a reasonable description of the collective plasmon excitation and which resums only a small subset of the polarizability diagrams included in EOM-CCSD. We find that EOM-CCSD, instead of providing a perturbative improvement on the RPA plasmon, predicts a many-state plasmon resonance, where each contributing state has a double-excitation character of 80% or more. This finding amounts to an ab initio treatment of the plasmon linewidth, which is in good quantitative agreement with previous diagrammatic calculations, and highlights the strongly correlated nature of lifetime effects in condensed-phase electronic structure theory.

Ab initio Path Integral Monte Carlo Results for the Dynamic Structure Factor of Correlated Electrons: From the Electron Liquid to Warm Dense Matter

The accurate description of electrons at extreme density and temperature is of paramount importance for, e.g., the understanding of astrophysical objects and inertial confinement fusion. In this context, the dynamic structure factor S(q,ω) constitutes a key quantity as it is directly measured in x-ray Thomson scattering experiments and governs transport properties like the dynamic conductivity. In this work, we present the first ab initio results for S(q,ω) by carrying out extensive path integral Monte Carlo simulations and developing a new method for the required analytic continuation, which is based on the stochastic sampling of the dynamic local field correction G(q,ω). In addition, we find that the so-called static approximation constitutes a promising opportunity to obtain high-quality data for S(q,ω) over substantial parts of the warm dense matter regime.

CHEERS: a tool for correlated hole-electron evolution from real-time simulations

We put forward a practical nonequilibrium Green's function (NEGF) scheme to perform real-time evolutions of many-body interacting systems driven out of equilibrium by external fields. CHEERS is a computational tool to solve the NEGF equation of motion in the so called generalized Kadanoff-Baym ansatz and it can be used for model systems as well as first-principles Hamiltonians. Dynamical correlation (or memory) effects are added to the Hartree-Fock dynamics through a many-body self-energy. Applications to time-dependent quantum transport, time-resolved photoabsorption and other ultrafast phenomena are discussed.

Charge Separation in Donor-C60 Complexes with Real-Time Green Functions: The Importance of Nonlocal Correlations

We use the nonequilibrium Green function (NEGF) method to perform real-time simulations of the ultrafast electron dynamics of photoexcited donor-C₆₀ complexes modeled by a Pariser-Parr-Pople Hamiltonian. The NEGF results are compared to mean-field Hartree-Fock (HF) calculations to disentangle the role of correlations. Initial benchmarking against numerically highly accurate time-dependent density matrix renormalization group calculations verifies the accuracy of NEGF. We then find that charge-transfer (CT) excitons partially decay into charge separated (CS) states if dynamical nonlocal correlation corrections are included. This CS process occurs in ∼10 fs after photoexcitation. In contrast, the probability of exciton recombination is almost 100% in HF simulations. These results are largely unaffected by nuclear vibrations; the latter become however essential whenever level misalignment hinders the CT process. The robust nature of our findings indicates that ultrafast CS driven by correlation-induced decoherence may occur in many organic nanoscale systems, but it will only be correctly predicted by theoretical treatments that include time-nonlocal correlations.

Quantum theory for the dynamic structure factor in correlated two-component systems in nonequilibrium: Application to x-ray scattering

We present a quantum theory for the dynamic structure factors in nonequilibrium, correlated, two-component systems such as plasmas or warm dense matter. The polarization function, which is needed as the input for the calculation of the structure factors, is calculated in nonequilibrium based on a perturbation expansion in the interaction strength. To make our theory applicable for x-ray scattering, a generalized Chihara decomposition for the total electron structure factor in nonequilibrium is derived. Examples are given and the influence of correlations and exchange on the structure and the x-ray-scattering spectrum are discussed for a model nonequilibrium distribution, as often encountered during laser heating of materials, as well as for two-temperature systems.

Permutation-blocking path-integral Monte Carlo approach to the static density response of the warm dense electron gas

The static density response of the uniform electron gas is of fundamental importance for numerous applications. Here we employ the recently developed ab initio permutation blocking path integral Monte Carlo (PB-PIMC) technique [T. Dornheim et al., New J. Phys. 17, 073017 (2015)10.1088/1367-2630/17/7/073017] 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 us 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 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.

Vertex Corrections for Positive-Definite Spectral Functions of Simple Metals

We present a systematic study of vertex corrections in a homogeneous electron gas at metallic densities. The vertex diagrams are built using a recently proposed positive-definite diagrammatic expansion for the spectral function. The vertex function not only provides corrections to the well known plasmon and particle-hole scatterings, but also gives rise to new physical processes such as the generation of two plasmon excitations or the decay of the one-particle state into a two-particle-one-hole state. By an efficient Monte Carlo momentum integration we are able to show that the additional scattering channels are responsible for a reduction of the bandwidth, the appearance of a secondary plasmon satellite below the Fermi level, and a substantial redistribution of spectral weights. The feasibility of the approach for first-principles band-structure calculations is also discussed.

Theory of Thomson scattering in inhomogeneous media

Thomson scattering of laser light is one of the most fundamental diagnostics of plasma density, temperature and magnetic fields. It relies on the assumption that the properties in the probed volume are homogeneous and constant during the probing time. On the other hand, laboratory plasmas are seldom uniform and homogeneous on the temporal and spatial dimensions over which data is collected. This is particularly true for laser-produced high-energy-density matter, which often exhibits steep gradients in temperature, density and pressure, on a scale determined by the laser focus. Here, we discuss the modification of the cross section for Thomson scattering in fully-ionized media exhibiting steep spatial inhomogeneities and/or fast temporal fluctuations. We show that the predicted Thomson scattering spectra are greatly altered compared to the uniform case, and may lead to violations of detailed balance. Therefore, careful interpretation of the spectra is necessary for spatially or temporally inhomogeneous systems.

Ion potential in warm dense matter: wake effects due to streaming degenerate electrons

The effective dynamically screened potential of a classical ion in a stationary flowing quantum plasma at finite temperature is investigated. This is a key quantity for thermodynamics and transport of dense plasmas in the warm-dense-matter regime. This potential has been studied before within hydrodynamic approaches or based on the zero temperature Lindhard dielectric function. Here we extend the kinetic analysis by including the effects of finite temperature and of collisions based on the Mermin dielectric function. The resulting ion potential exhibits an oscillatory structure with attractive minima (wakes) and, thus, strongly deviates from the static Yukawa potential of equilibrium plasmas. This potential is analyzed in detail for high-density plasmas with values of the Brueckner parameter in the range 0.1≤r(s)≤1 for a broad range of plasma temperature and electron streaming velocity. It is shown that wake effects become weaker with increasing temperature of the electrons. Finally, we obtain the minimal electron streaming velocity for which attraction between ions occurs. This velocity turns out to be less than the electron Fermi velocity. Our results allow for reliable predictions of the strength of wake effects in nonequilibrium quantum plasmas with fast streaming electrons showing that these effects are crucial for transport under warm-dense-matter conditions, in particular for laser-matter interaction, electron-ion temperature equilibration, and stopping power.

Quasiparticle parametrization of mean fields, Galilei invariance, and universal conserving response functions

The general possible form of mean-field parametrization in a running frame in terms of current, energy, and density functionals is examined under the restrictions of Galilean invariance. It is found that only two density-dependent parameters remain which are usually condensed in a position-dependent effective mass and the self-energy formed by current and mass. The position-dependent mass induces a position-dependent local current, which is identified for different nonlinear frames. In a second step the response to an external perturbation and relaxation towards a local equilibrium is investigated. The response function is found to be universal in the sense that the actual parametrization of the local equilibrium does not matter and is eliminated from the theory due to the conservation laws. The explicit form of the response with respect to density, momentum, and energy is derived. The compressibility sum rule as well as the sum rule by first- and third-order frequency moments are proved analytically to be fulfilled simultaneously. The results are presented for Bose or Fermi systems in one, two, and three dimensions.

Time-dependent screening of a point charge at a metal surface

The space-time evolution of the dynamical screening charge density caused by a suddenly created point charge at the Cu(111) surface is investigated in the linear response approximation. Considering a thin slab as a model for the Cu(111) surface, we investigate the confinement effects on dynamical screening as well. The results have been obtained on the basis of self-consistent evaluation of the energy-momentum-dependent response function, taking into account the realistic surface band structure of Cu(111). At the initial stage, we observe fast long-range charge density oscillations due to excitation of the surface plasmon modes. Then we observe the propagation of the shock wave of the electron-hole excitations along the slab with velocity determined by the Fermi velocity of bulk Cu. At longer times, we have identified the propagation along the two slab surfaces of a much slower (with velocity ∼ 0.3 au, close to the Fermi velocity of the Cu(111) surface state) charge disturbance due to acoustic surface plasmon. The role of the energy band gap in the direction perpendicular to the surface in establishing the screening is also addressed.

Nonlinear magnetoplasmons in strongly coupled Yukawa plasmas

The existence of plasma oscillations at multiples of the magnetoplasmon frequency in a strongly coupled two-dimensional magnetized Yukawa plasma is reported, based on extensive molecular dynamics simulations. These modes are the analogues of Bernstein modes which are renormalized by strong interparticle correlations. Their properties are theoretically explained by a dielectric function incorporating the combined effect of a magnetic field, strong correlations and finite temperature.

Successes and failures of Kadanoff-Baym dynamics in Hubbard nanoclusters

We study the nonequilibrium dynamics of small, strongly correlated clusters, described by a Hubbard Hamiltonian, by propagating in time the Kadanoff-Baym equations within the Hartree-Fock, second Born, GW, and T-matrix approximations. We compare the results to exact numerical solutions. We find that the time-dependent T matrix is overall superior to the other approximations, and is in good agreement with the exact results in the low-density regime. In the long time limit, the many-body approximations attain an unphysical steady state which we attribute to the implicit inclusion of infinite-order diagrams in a few-body system.

Two-particle irreducible effective action approach to nonlinear current-conserving approximations in driven systems

Using closed-time-path two-particle irreducible coarse-grained effective action (CTP 2PI CGEA) techniques, we study the response of an open interacting electronic system to time-dependent external electromagnetic fields. We show that the CTP 2PI CGEA is invariant under a simultaneous gauge transformation of the external field and the full Schwinger-Keldysh propagator, and that this property holds even when the loop expansion of the CTP 2PI CGEA is truncated at arbitrary order. The effective action approach provides a systematic way of calculating the propagator and response functions of the system, via the Schwinger-Dyson equation and the Bethe-Salpeter equations, respectively. We show that, due to the invariance of the CTP 2PI CGEA under external gauge transformations, the response functions calculated from it satisfy the Ward-Takahashi hierarchy, thus warranting the conservation of the electronic current beyond the expectation value level. We also clarify the connection between nonlinear response theory and the WT hierarchy, and discuss an example of an ad hoc approximation that violates it. These findings may be useful in the study of current fluctuations in correlated electronic pumping devices.

Solving the Kadanoff-Baym equations for inhomogeneous systems: application to atoms and molecules

We implement time propagation of the nonequilibrium Green function for atoms and molecules by solving the Kadanoff-Baym equations within a conserving self-energy approximation. We here demonstrate the usefulness of time propagation for calculating spectral functions and for describing the correlated electron dynamics in a nonperturbative electric field. We also demonstrate the use of time propagation as a method for calculating charge-neutral excitation energies, equivalent to highly advanced solutions of the Bethe-Salpeter equation.

How many-particle interactions develop after ultrafast excitation of an electron-hole plasma

Electrostatic coupling between particles is important in many microscopic phenomena found in nature. The interaction between two isolated point charges is described by the bare Coulomb potential, but in many-body systems this interaction is modified as a result of the collective response of the screening cloud surrounding each charge carrier. One such system involves ultrafast interactions between quasi-free electrons in semiconductors-which are central to high-speed and future quantum electronic devices. The femtosecond kinetics of nonequilibrium Coulomb systems has been calculated using static and dynamical screening models that assume the instantaneous formation of interparticle correlations. However, some quantum kinetic theories suggest that a regime of unscreened bare Coulomb collisions might exist on ultrashort timescales. Here we monitor directly the temporal evolution of the charge-charge interactions after ultrafast excitation of an electron-hole plasma in GaAs. We show that the onset of collective behaviour such as Coulomb screening and plasmon scattering exhibits a distinct time delay of the order of the inverse plasma frequency, that is, several 10(-14) seconds.

Harmonics generation in electron-ion collisions in a short laser pulse

Anomalously high generation efficiency of coherent higher field harmonics in collisions between oppositely charged particles in the field of femtosecond lasers is predicted. This is based on rigorous numerical solutions of a quantum kinetic equation for dense laser plasmas that overcomes limitations of previous investigations.

Dynamical properties and plasmon dispersion of a weakly degenerate correlated one-component plasma

Classical Molecular Dynamics simulations for a one-component plasma are presented. Quantum effects are included in the form of the Kelbg potential. Results for the dynamical structure factor are compared with the Vlasov and random phase approximation theories. The influence of the coupling parameter Gamma, degeneracy parameter rho Lambda(3), and the form of the pair interaction on the optical plasmon dispersion is investigated. An improved analytical approximation for the dispersion of Langmuir waves is presented.