Search

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.

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.

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.

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.

Solids facing a plasma are a common situation in many astrophysical systems and laboratory setups. Moreover, many plasma technology applications rely on the control of the plasma-surface interaction. However, presently often a fundamental understanding of them is missing, so most technological applications are being developed via trial and error. In the majority of plasma simulations surface processes are either neglected or treated via phenomenological parameters such as sticking coefficients, sputter rates or secondary electron emission coefficients. However, those parameters are known only in some cases and with very limited accuracy. Similarly, while surface physics simulations have often studied the impact of single ions or neutrals, so far, the influence of a plasma medium and correlations between successive impacts have not been taken into account. Such an approach cannot have predictive power. In this paper we discuss in some detail the physical processes a the plasma-solid interface which brings us to the necessity of coupled plasma-solid simulations. We briefly summarize relevant theoretical methods from solid state and surface physics that are suitable to contribute to such an approach and identify four methods. The first are mesoscopic simulations such as kinetic Monte Carlo (KMC) and molecular dynamics (MD) that are able to treat complex processes on large scales but neglect electronic effects. The second are quantum kinetic methods based on the quantum Boltzmann equation that give access to a more accurate treatment of surface processes using simplifying models for the solid. The third approach are ab initio simulations of surface process that are based on density functional theory (DFT) and time-dependent DFT. The fourths are nonequilibrium Green functions that able to treat correlation effects in the material and at the interface.

In a recent letter, Dornheim et al. [Phys. Rev. Lett. 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 paper, 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 tool of diagnostics. To this end, we present a theoretical description of the nonlinear electronic density response based on the recent effective static approximation to the local field correction [T. Dornheim et al., Phys. Rev. Lett. 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.

The energy gap of correlated Hubbard clusters is well studied for one-dimensional systems using analytical methods and density-matrix- renormalization-group (DMRG) simulations. Beyond 1D, however, exact results are available only for small systems by quantum Monte Carlo. For this reason and, due to the problems of DMRG in simulating 2D and 3D systems, alternative methods such as Green functions combined with many-body approximations (GFMBA), that do not have this restriction, are highly important. However, it has remained open whether the approximate character of GFMBA simulations prevents the computation of the Hubbard gap. Here we present new GFMBA results that demonstrate that GFMBA simulations are capable of producing reliable data for the gap which agrees well with the DMRG benchmarks in 1D. An interesting observation is that the accuracy of the gap can be significantly increased when the simulations give up certain symmetry restriction of the exact system, such as spin symmetry and spatial homogeneity. This is seen as manifestation and generalization of the “symmetry dilemma” introduced by Lowdin for Hartree–Fock wave function calculations.

We present a finite-temperature density functional theory investigation of the nonequilibrium transient electronic structure of warm dense Li, Al, Cu, and Au created by laser excitation. Photons excite electrons either from the inner shell orbitals or from the valence bands according to the photon energy, and give rise to isochoric heating of the sample. Localized states related to the 3d orbital are observed for Cu when the hole lies in the inner shell 3s orbital. The electrical conductivity for these materials at nonequilibrium states is calculated using the Kubo-Greenwood formula. The change of the electrical conductivity, compared to the equilibrium state, is different for the case of holes in inner shell orbitals or the valence band. This is attributed to the competition of two factors: the shift of the orbital energies due to reduced screening of core electrons, and the increase of chemical potential due to the excitation of electrons. The finite temperature effect of both the electrons and the ions on the electrical conductivity is discussed in detail. This work is helpful to better understand the physics of laser excitation experiments of warm dense matter., Comment: 111 pages, 9 figures

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.

Quantum Monte Carlo (QMC) belongs to the most accurate simulation techniques for quantum many-particle systems. However, for fermions, these simulations are hampered by the sign problem that prohibits simulations in the regime of strong degeneracy. The situation changed with the development of configuration path integral Monte Carlo (CPIMC) by Schoof et al. [Contrib. Plasma Phys. 51, 687 (2011)] that allowed for the first ab initio simulations for dense quantum plasmas [Schoof et al., Phys. Rev. Lett. 115, 130402 (2015)]. CPIMC also has a sign problem that occurs when the density is lowered, i.e., in a parameter range that is complementary to traditional QMC formulated in coordinate space. Thus, CPIMC simulations for the warm dense electron gas are limited to small values of the Brueckner parameter—the ratio of the interparticle distance to the Bohr radius—rs=r⎯⎯⎯/aB≲1 . In order to reach the regime of stronger coupling (lower density) with CPIMC, here we investigate additional restrictions on the Monte Carlo procedure. In particular, we introduce two different versions of “restricted CPIMC”—called RCPIMC and RCPIMC+—where certain sign changing Monte Carlo updates are being omitted. Interestingly, one of the methods (RCPIMC) has no sign problem at all, but it introduces a systematic error and is less accurate than RCPIMC+, which neglects only a smaller class of the Monte Carlo steps. Here, we report extensive simulations for the ferromagnetic uniform electron gas with which we investigate the properties and accuracy of RCPIMC and RCPIMC+. Furthermore, we establish the parameter range in the density–temperature plane where these simulations are both feasible and accurate. The conclusion is that RCPIMC and RCPIMC+ work best at temperatures in the range of Θ = kBT/EF ∼ 0.1…0.5, where EF is the Fermi energy, allowing to reach density parameters up to rs ∼ 3…5, thereby partially filling a gap left open by existing ab initio QMC methods.

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 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 was reported by Dornheim et al. [Phys. Rev. Lett. 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.

The time evolution in quantum many-body systems after external excitations is attracting high interest in many fields. The theoretical modeling of these processes is challenging, and the only rigorous quantum-dynamics approach that can treat correlated fermions 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), for the second-order Born (SOA) selfenergy, which has substantially extended the scope of NEGF simulations. In a recent Letter [Schlunzen \textit{et al.}, Phys. Rev. Lett. \textbf{124}, 076601 (2020)] we demonstrated that GKBA-NEGF simulations can be efficiently mapped onto coupled time-local equations for the single-particle and two-particle Green functions on the time diagonal, hence the method has been called G1--G2 scheme. This allows one to perform the same simulations with order $T^1$-scaling, both for SOA and $GW$ selfenergies giving rise to a dramatic speedup. Here we present more details on the G1--G2 scheme, including derivations of the basic equations including results for a general basis, for Hubbard systems and for jellium. Also, we demonstrate how to incorporate initial correlations into the G1--G2 scheme. Further, the derivations are extended to a broader class of selfenergies, including the $T$ matrix in the particle--particle and particle--hole channels, and the dynamically screened-ladder approximation. Finally, we demonstrate that, for all selfenergies, the CPU time scaling of the G1--G2 scheme with the basis dimension, $N_b$, can be improved compared to our first report: the overhead compared to the original GKBA, is not more than an additional factor $N_b$.

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.

We investigate the energy-loss characteristics of an ion in warm dense matter (WDM) and dense plasmas concentrating on the influence of electronic correlations. The basis for our analysis is a recently developed ab initio quantum Monte Carlo- (QMC) based machine learning representation of the static local field correction (LFC) [Dornheim et al., J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013], which provides an accurate description of the dynamical density response function of the electron gas at the considered parameters. We focus on the polarization-induced stopping power due to free electrons, the friction function, and the straggling rate. In addition, we compute the friction coefficient which constitutes a key quantity for the adequate Langevin dynamics simulation of ions. Considering typical experimental WDM parameters with partially degenerate electrons, we find that the friction coefficient is of the order of γ/ω_{pi}=0.01, where ω_{pi} is the ionic plasma frequency. This analysis is performed by comparing QMC-based data to results from the random-phase approximation (RPA), the Mermin dielectric function, and the Singwi-Tosi-Land-Sjolander (STLS) approximation. It is revealed that the widely used relaxation time approximation (Mermin dielectric function) has severe limitations regarding the description of the energy loss of ions in a correlated partially degenerate electrons gas. Moreover, by comparing QMC-based data with the results obtained using STLS, we find that the ion energy-loss properties are not sensitive to the inaccuracy of the static local field correction (LFC) at large wave numbers, k/k_{F}>2 (with k_{F} being the Fermi wave number), but that a correct description of the static LFC at k/k_{F}≲1.5 is important.

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.

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.

Motivated by experiments on the generation of streaming plasmas in high energy density facilities,industrial setups, and fundamental dusty plasma research, the plasma polarization around a test charge instreaming plasmas is considered in this work. The induced charge density distribution of the plasmaconstituents is discussed for the subsonic, sonic, and supersonic regime taking into account the non-Maxwellian distribution of the flowing ions. It is shown that the plasma polarization (the plasmawakefield) in the vicinity of the test charge shows different scaling in the subsonic and supersonicregimes, where Mach number is defined as the ratio of the ion streaming velocity and the ion soundspeed. In contrast to the wake potential, the density decays strongly monotonically in the plasma wakeand does not exhibit an oscillatory pattern or trailing maxima. Therefore, the picture of an ion focusingeffect creating a separated ion region downstream was not confirmed.

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 structure factor. By combining exact quantum Monte Carlo results for large wave vectors with the long-range behaviour from the Singwi-Tosi-Land-Sjolander 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.

From the static polarization function of electrons in the random phase approximation, the quantum Bohm potential for the quantum hydrodynamic description of electrons and the density gradient correction to the Thomas–Fermi free energy at a finite temperature for the two- and one-dimensional cases are derived. The behaviour of the Bohm potential and of the density gradient correction as a function of the degeneracy parameter is discussed. Based on recent developments in the fluid description of quantum plasmas, the Bohm potential for the high-frequency domain is presented.

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.

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

Leonid V. Keldysh–one of the most influential theoretical physicists of the 20th century–passed away in November 2016. L. V. 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 L. V. Keldysh had his annus mirabilis – the year1964. But the list of his theoretical developments is much broader and is briefly reviewed here.

The dynamical structure factor (DSF) of strongly coupled ions in dense plasmas with partially and strongly degenerate electrons is investigated. The main focus is on the impact of electronic correlations (non-ideality) on the ionic DSF. The latter is computed by carrying out molecular dynamics (MD) simulations with a screened ion-ion interaction potential. The electronic screening is taken into account by invoking the Singwi-Tosi-Land-Sj\"olander approximation, and compared to the MD simulation data obtained considering the electronic screening in the random phase approximation and using the Yukawa potential. %This allows us to gain insight into the impact of the electronic non-ideality on the ionic DSF. We find that electronic correlations lead to lower values of the ion-acoustic mode frequencies and to an extension of the applicability limit with respect to the wave-number of a hydrodynamic description. Moreover, we show that even in the limit of weak electronic coupling, electronic correlations have a non-negligible impact on the ionic longitudinal sound speed. Additionally, the applicability of the Yukawa potential with an adjustable screening parameter is discussed, which will be of interest, e.g., for the interpretation of experimental results for the ionic DSF of dense plasmas.

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 response of the uniform electron gas (UEG) to an external perturbation is of paramount importance for many applications. Recently, highly accurate results for the static density response function and the corresponding local field correction have been provided both for warm dense matter [J. Chem. Phys. 151, 194 104 (2019)] and strongly coupled electron liquid [Phys. Rev. B 101, 045 129 (2020)] conditions based on exact ab initio path integral Monte Carlo (PIMC) simulations. In the present work, we further complete our current description of the UEG by exploring the high energy density regime, which is relevant for, e.g. astrophysical applications and inertial confinement fusion experiments. To this end, we present extensive new PIMC results for the static density response in the range of 0.05 ≤ r s ≤ 0.5 and 0.85 ≤ θ ≤ 8. These data are subsequently used to benchmark the accuracy of the widely used random phase approximation and the dielectric theory by Singwi, Tosi, Land, and Sjölander (STLS). Moreover, we compare our results to configuration PIMC data where they are available and find perfect agreement with a relative accuracy of 0.001 − 0.01%. All PIMC data are available online.

The energy-autocorrelation function is calculated for a screened Coulomb system in equilibrium and the contributions of different energy transport channels to the total heat conductivity are explored for magnetized and unmagnetized systems. A special focus is on the time scales of the energy-autocorrelation function which contribute to the field-parallel enhancement of heat conduction in strongly coupled plasmas. The investigation is based on first-principle molecular dynamics simulations. (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

The density response function for streaming ions in homogeneous, parallel electric and magnetic fields is derived self-consistently from kinetic theory. Ion-neutral collisions are treated with the Bhatnagar-Gross-Krook collision operator assuming a constant ion-neutral collision frequency. The result accounts for the non-Maxwellian distribution function of the ions and is valid in the full range from weak to strong magnetization. It provides the basis for various linear response calculations in the context of magnetized complex plasmas, where streaming ions interact with highly charged dust particles under the influence of a strong external magnetic field. (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

The ion potential in the warm dense matter regime exhibits wake effects due to streaming degenerate electrons and has been discussed previously [Phys. Rev. E 91, 023102 (2015)]. Here, we extend the analysis with particular focus on anomalous wake effects that is (i) the collision-induced wake amplification, and (ii) the non-monotonic temperature dependence of the wake amplitude. (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)