Search

Abstract We present a formally exact and simulation-free approach for the normalization of X-ray Thomson scattering (XRTS) spectra based on the f-sum rule of the imaginary-time correlation function (ITCF). Our method works for any degree of collectivity, over a broad range of temperatures, and is applicable even in nonequilibrium situations. In addition to giving us model-free access to electronic correlations, this new approach opens up the intriguing possibility to extract a plethora of physical properties from the ITCF based on XRTS experiments.

We demonstrate a significantly simplified experimental approach for investigating liquid metallic hydrogen, which is crucial to understand the internal structure and evolution of giant planets. Plastic samples were shockcompressed and then probed by short pulses of X-rays generated by free electron lasers. By comparison with ab initio simulations, we provide indirect evidence for the creation of elemental hydrogen in shock-compressed plastics at ∼150GPa and ∼5,000K and thus in a regime where hydrogen is predicted to be metallic. Being the most common form of condensed matter in our solar system, and ostensibly the simplest of all elements, hydrogen is the model case for many theoretical studies and we provide a new possibility to benchmark models for conditions with extreme pressures and temperatures. Moreover, this approach will also allow to probe the chemical behavior of metallic hydrogen in mixture with other elements, which, besides its importance for planetary physics, may open up promising pathways for the synthesis of new materials.

Charged particle interaction and energy dissipation in plasma is fundamentally interesting. Here the authors study proton stopping in laser-produced plasma for the moderate to strong coupling with electrons.

Abstract We present structure and equation of state (EOS) measurements of biaxially orientated polyethylene terephthalate (PET, $$({\hbox {C}}_{10} {\hbox {H}}_8 {\hbox {O}}_4)_n$$ ( C 10 H 8 O 4 ) n , also called mylar) shock-compressed to ( $$155 \pm 20$$ 155 ± 20 ) GPa and ( $$6000 \pm 1000$$ 6000 ± 1000 ) K using in situ X-ray diffraction, Doppler velocimetry, and optical pyrometry. Comparing to density functional theory molecular dynamics (DFT-MD) simulations, we find a highly correlated liquid at conditions differing from predictions by some equations of state tables, which underlines the influence of complex chemical interactions in this regime. EOS calculations from ab initio DFT-MD simulations and shock Hugoniot measurements of density, pressure and temperature confirm the discrepancy to these tables and present an experimentally benchmarked correction to the description of PET as an exemplary material to represent the mixture of light elements at planetary interior conditions.

Angularly resolved X-ray scattering measurements from fs-laser heated hydrogen have been used to determine the equilibration of electron and ion temperatures in the warm dense matter regime. The relaxation of rapidly heated cryogenic hydrogen is visualized using 5.5 keV X-ray pulses from the Linac Coherent Light (LCLS) source in a 1 Hz repetition rate pump-probe setting. We demonstrate that the electron-ion energy transfer is faster than quasi-classical Landau-Spitzer models that use ad hoc cutoffs in the Coulomb logarithm.

It is challenging to reliably probe the miscibility behavior of elements in extreme conditions. Here, the authors use X-ray Thomson scattering and compare to the X-ray diffraction method in order to determine mixing of different atomic species in warm dense matter conditions.

The energy loss of ions in plasma is a challenging issue in inertial confinement fusion and many theoretical models exist on ion-stopping power. Here, the authors use laser-generated plasma probed by accelerator-produced ions in experiments to discriminate various ion stopping models near the Bragg peak.

Studying the properties of dense plasmas is challenging due to strong interactions between electrons and ions, and numerical methods overcome this difficulty using a static thermostat. Here the authors predict a strong diffusive ion mode at low energy by including dissipative processes in the model.

Average atom (AA) models allow one to efficiently compute electronic and optical properties of materials over a wide range of conditions and are often employed to interpret experimental data. However, at high pressure, predictions from AA models have been shown to disagree with results from ab initio computer simulations. Here we reconcile these deviations by developing an innovative type of AA model, AvIon, that computes the electronic eigenstates with novel boundary conditions within the ion sphere. Bound and free states are derived consistently. We drop the common AA image that the free-particle spectrum starts at the potential threshold, which we found to be incompatible with ab initio calculations. We perform ab initio simulations of crystalline and liquid carbon and aluminum over a wide range of densities and show that the computed band structure is in very good agreement with predictions from AvIon.

Shock synthesis of diamond and even harder carbon polymorphs from graphite is of great interest for science and technology. Here, the authors present unprecedented in situmeasurements of the structural changes, showing ultrafast formation of diamond and, at higher pressures, evidence for a pure lonsdaleite structure.

Using the Linac Coherent Light Source facility at the Stanford Linac Coherent Light Source National Accelerator Laboratory, we have observed x-ray scattering from iron compressed with laser-driven shocks to earth-core-like pressures above 400 GPa. The data show cases where melting is incomplete and we observe hexagonal-close-packed crystal structure at shock compressed densities up to 14.0 g cm^{−3} but no evidence of a double-hexagonal-close-packed crystal. The observation of a crystalline structure at these densities, where shock heating is expected to be in excess of the equilibrium melt temperature, may indicate superheating of the solid. These results are important for equation of state modeling at high strain rates relevant for impact scenarios and laser-driven shock-wave experiments.

We discuss the possibility of obtaining highly precise measurements of the ionization potential depression in dense plasmas with spectrally resolved x-ray scattering, while simultaneously determining the electron temperature and the free electron density. A proof-of-principle experiment at the Linac Coherent Light Source, probing isochorically heated carbon samples, demonstrates the capabilities of this method and motivates future experiments at x-ray free electron laser facilities. [ABSTRACT FROM AUTHOR]

We present results for the ionic structure in dense, moderately to strongly coupled plasmas using two models: the mean spherical approximation (MSA) and the hypernetted chain (HNC) approach. While the first method allows for an analytical solution, the latter has to be solved iteratively. Independent of the coupling strength, the results show only small differences when the ions are considered to form an unscreened one-component plasma (OCP) system. If the electrons are treated as a polarizable background, the different ways to incorporate the screening yield, however, large discrepancies between the models, particularly for more strongly coupled plasmas. [ABSTRACT FROM AUTHOR]

We present a formally exact and simulation-free approach for the normalization of X-ray Thomson scattering (XRTS) spectra based on the f-sum rule of the imaginary-time correlation function (ITCF). Our method works for any degree of collectivity, over a broad range of temperatures, and is applicable even in nonequilibrium situations. In addition to giving us model-free access to electronic correlations, this new approach opens up the intriguing possibility to extract a plethora of physical properties from the ITCF based on XRTS experiments., (© 2024. The Author(s).)

The combination of isochoric heating of solids by free-electron lasers (FELs) and in situ diagnostics by X-ray Thomson scattering (XRTS) allows for measurements of material properties at warm dense matter (WDM) conditions relevant for astrophysics, inertial confinement fusion, and materials science. In the case of metals, the FEL beam pumps energy directly into electrons with the lattice structure of ions being nearly unaffected. This leads to a unique transient state that gives rise to a set of interesting physical effects, which can serve as a reliable testing platform for WDM theories. In this work, we present extensive linear-response time-dependent density functional theory (TDDFT) results for the electronic dynamic structure factor of isochorically heated copper with a face-centered cubic lattice. At ambient conditions, the plasmon is heavily damped due to the presence of d-band excitations, and its position is independent of the wavenumber. In contrast, the plasmon feature starts to dominate the excitation spectrum and has a Bohm-Gross-type plasmon dispersion for temperatures T ≥ 4 eV, where the quasi-free electrons in the interstitial region are in the WDM regime. In addition, we analyze the thermal changes in the d-band excitations and outline the possibility to use future XRTS measurements of isochorically heated copper as a controlled testbed for WDM theories., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

Optical femtosecond pump-probe experiments allow to measure the dynamics of ultrafast heating of metals with high accuracy. However, the theoretical analysis of such experiments is often complicated because of the indirect connection of the measured signal and the desired temperature transients. Establishing such a connection requires an accurate model of the optical constants of a metal, depending on both the electron temperature T e and the lattice temperature T l . In this paper, we present first-principles simulations of the two-temperature scenario with T e  ≫ T l , showing the optical response of hot electrons to laser irradiation in gold and ruthenium. Comparing our simulations with the Kubo-Greenwood approach, we discuss the influence of electron-phonon and electron-electron scattering on the intraband contribution to optical constants. Applying the simulated optical constants to the analysis of ultrafast heating of ruthenium thin films we highlight the importance of the latter scattering channel to understand the measured heating dynamics.

We present extensive new ab initio path integral Monte Carlo (PIMC) results for a variety of structural properties of warm dense hydrogen and beryllium. To deal with the fermion sign problem-an exponential computational bottleneck due to the antisymmetry of the electronic thermal density matrix-we employ the recently proposed [Y. Xiong and H. Xiong, J. Chem. Phys. 157, 094112 (2022); T. Dornheim et al., J. Chem. Phys. 159, 164113 (2023)] ξ-extrapolation method and find excellent agreement with the exact direct PIMC reference data where available. This opens up the intriguing possibility of studying a gamut of properties of light elements and potentially material mixtures over a substantial part of the warm dense matter regime, with direct relevance for astrophysics, material science, and inertial confinement fusion research., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

The accurate description of non-ideal quantum many-body systems is of prime importance for a host of applications within physics, quantum chemistry, materials science, and related disciplines. At finite temperatures, the gold standard is given by ab initio path integral Monte Carlo (PIMC) simulations, which do not require any empirical input but exhibit an exponential increase in the required computation time for Fermionic systems with an increase in system size N . Very recently, computing Fermionic properties without this bottleneck based on PIMC simulations of fictitious identical particles has been suggested. In our work, we use this technique to perform very large ( N ≤ 1000) PIMC simulations of the warm dense electron gas and demonstrate that it is capable of providing a highly accurate description of the investigated properties, i.e., the static structure factor, the static density response function, and the local field correction, over the entire range of length scales.

The properties of plasmas in the low-density limit are described by virial expansions. Analytical expressions are known from Green's function approaches only for the first three virial coefficients. Accurate path-integral Monte Carlo (PIMC) simulations have recently been performed for the uniform electron gas, allowing the virial expansions to be analyzed and interpolation formulas to be derived. The exact expression for the second virial coefficient is used to test the accuracy of the PIMC simulations and the range of validity of the interpolation formula of Groth et al. [Phys. Rev. Lett. 119, 135001 (2017)0031-900710.1103/PhysRevLett.119.135001], and we discuss the fourth virial coefficient, which is not exactly known yet. Combining PIMC simulations with benchmarks from exact virial expansion results would allow us to obtain more accurate representations of the equation of state for parameter ranges of conditions which are of interest, e.g., for helioseismology.

Hydrogen at extreme temperatures and pressures is of key relevance for cutting-edge technological applications, with inertial confinement fusion research being a prime example. In addition, it is ubiquitous throughout our universe and naturally occurs in a variety of astrophysical objects. In the present work, we present exact ab initio path integral Monte Carlo (PIMC) results for the electronic density of warm dense hydrogen along a line of constant degeneracy across a broad range of densities. Using the well-known concept of reduced density gradients, we develop a new framework to identify the breaking of bound states due to pressure ionization in bulk hydrogen. Moreover, we use our PIMC results as a reference to rigorously assess the accuracy of a variety of exchange-correlation (XC) functionals in density functional theory calculations for different density regions. Here, a key finding is the importance of thermal XC effects for the accurate description of density gradients in high-energy-density systems. Our exact PIMC test set is freely available online and can be used to guide the development of new methodologies for the simulation of warm dense matter and beyond.

The ab initio path integral Monte Carlo (PIMC) method is one of the most successful methods in statistical physics, quantum chemistry and related fields, but its application to quantum degenerate Fermi systems is severely hampered by an exponential computational bottleneck: the notorious fermion sign problem. Very recently, Xiong and Xiong [J. Chem. Phys. 157, 094112 (2022)] have suggested to partially circumvent the sign problem by carrying out simulations of fictitious systems guided by an interpolating continuous variable ξ ∈ [-1, 1], with the physical Fermi- and Bose-statistics corresponding to ξ = -1 and ξ = 1. It has been proposed that information about the fermionic limit might be obtained by calculations within the bosonic sector ξ > 0 combined with an extrapolation throughout the fermionic sector ξ < 0, essentially bypassing the sign problem. Here, we show how the inclusion of the artificial parameter ξ can be interpreted as an effective penalty on the formation of permutation cycles in the PIMC simulation. We demonstrate that the proposed extrapolation method breaks down for moderate to high quantum degeneracy. Instead, the method constitutes a valuable tool for the description of large Fermi-systems of weak quantum degeneracy. This is demonstrated for electrons in a 2D harmonic trap and for the uniform electron gas (UEG), where we find excellent agreement (∼0.5%) with exact configuration PIMC results in the high-density regime while attaining a speed-up exceeding 11 orders of magnitude. Finally, we extend the idea beyond the energy and analyze the radial density distribution (2D trap), as well as the static structure factor and imaginary-time density-density correlation function (UEG)., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

The properties of hydrogen at warm dense matter (WDM) conditions are of high importance for the understanding of astrophysical objects and technological applications such as inertial confinement fusion. In this work, we present extensive ab initio path integral Monte Carlo results for the electronic properties in the Coulomb potential of a fixed ionic configuration. This gives us unique insights into the complex interplay between the electronic localization around the protons with their density response to an external harmonic perturbation. We find qualitative agreement between our simulation data and a heuristic model based on the assumption of a local uniform electron gas model, but important trends are not captured by this simplification. In addition to being interesting in their own right, we are convinced that our results will be of high value for future projects, such as the rigorous benchmarking of approximate theories for the simulation of WDM, most notably density functional theory.

Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [ Nat. Commun. 13 , 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density-density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Linear-response time-dependent density functional theory (LR-TDDFT) simulations of disordered extended systems require averaging over different snapshots of ion configurations to minimize finite size effects due to the snapshot-dependence of the electronic density response function and related properties. We present a consistent scheme for the computation of the macroscopic Kohn-Sham (KS) density response function connecting an average over snapshot values of charge density perturbations to the averaged values of KS potential variations. This allows us to formulate the LR-TDDFT within the adiabatic (static) approximation for the exchange-correlation (XC) kernel for disordered systems, where the static XC kernel is computed using the direct perturbation method [Moldabekov et al., J. Chem. Theory Comput. 19, 1286 (2023)]. The presented approach allows one to compute the macroscopic dynamic density response function as well as the dielectric function with a static XC kernel generated for any available XC functional. The application of the developed workflow is demonstrated for the example of warm dense hydrogen. The presented approach is applicable for various types of extended disordered systems, such as warm dense matter, liquid metals, and dense plasmas., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

The gravitational pressure in many astrophysical objects exceeds one gigabar (one billion atmospheres) 1-3 , creating extreme conditions where the distance between nuclei approaches the size of the K shell. This close proximity modifies these tightly bound states and, above a certain pressure, drives them into a delocalized state 4 . Both processes substantially affect the equation of state and radiation transport and, therefore, the structure and evolution of these objects. Still, our understanding of this transition is far from satisfactory and experimental data are sparse. Here we report on experiments that create and diagnose matter at pressures exceeding three gigabars at the National Ignition Facility 5 where 184 laser beams imploded a beryllium shell. Bright X-ray flashes enable precision radiography and X-ray Thomson scattering that reveal both the macroscopic conditions and the microscopic states. The data show clear signs of quantum-degenerate electrons in states reaching 30 times compression, and a temperature of around two million kelvins. At the most extreme conditions, we observe strongly reduced elastic scattering, which mainly originates from K-shell electrons. We attribute this reduction to the onset of delocalization of the remaining K-shell electron. With this interpretation, the ion charge inferred from the scattering data agrees well with ab initio simulations, but it is significantly higher than widely used analytical models predict 6 ., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

We study the linear energy response of the uniform electron gas to an external harmonic perturbation with a focus on resolving different contributions to the total energy. This has been achieved by carrying out highly accurate ab initio path integral Monte Carlo (PIMC) calculations for a variety of densities and temperatures. We report a number of physical insights into effects such as screening and the relative importance of kinetic and potential energies for different wave numbers. A particularly interesting finding is obtained from the observed non-monotonic behavior of the induced change in the interaction energy, which becomes negative for intermediate wave numbers. This effect is strongly dependent on the coupling strength and constitutes further direct evidence for the spatial alignment of electrons introduced in earlier works [T. Dornheim et al., Commun. Phys. 5, 304 (2022)]. The observed quadratic dependence on the perturbation amplitude in the limit of weak perturbations and the quartic dependence of perturbation amplitude corrections are consistent with linear and nonlinear versions of the density stiffness theorem. All PIMC simulation results are freely available online and can be used to benchmark new methods or as input for other calculations., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

We assess the accuracy of common hybrid exchange-correlation (XC) functionals (PBE0, PBE0-1/3, HSE06, HSE03, and B3LYP) within the Kohn-Sham density functional theory for the harmonically perturbed electron gas at parameters relevant for the challenging conditions of the warm dense matter. Generated by laser-induced compression and heating in the laboratory, the warm dense matter is a state of matter that also occurs in white dwarfs and planetary interiors. We consider both weak and strong degrees of density inhomogeneity induced by the external field at various wavenumbers. We perform an error analysis by comparing with the exact quantum Monte Carlo results. In the case of a weak perturbation, we report the static linear density response function and the static XC kernel at a metallic density for both the degenerate ground-state limit and for partial degeneracy at the electronic Fermi temperature. Overall, we observe an improvement in the density response when the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are used, compared with the previously reported results for the PBE, PBEsol, local-density approximation, and AM05 functionals; B3LYP, on the other hand, does not perform well for the considered system. Additionally, the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are more accurate for the density response properties than SCAN in the regime of partial degeneracy.

The electronic exchange─correlation (XC) kernel constitutes a fundamental input for the estimation of a gamut of properties such as the dielectric characteristics, the thermal and electrical conductivity, or the response to an external perturbation. In this work, we present a formally exact methodology for the computation of the system specific static XC kernel exclusively within the framework of density functional theory (DFT) and without employing functional derivatives─no external input apart from the usual XC-functional is required. We compare our new results with exact quantum Monte Carlo (QMC) data for the archetypical uniform electron gas model under both ambient and warm dense matter conditions. This gives us unprecedented insights into the performance of different XC functionals, and it has important implications for the development of new functionals that are designed for the application at extreme temperatures. In addition, we obtain new DFT results for the XC kernel of warm dense hydrogen as it occurs in fusion applications and astrophysical objects. The observed excellent agreement to the QMC reference data demonstrates that presented framework is capable to capture nontrivial effects such as XC-induced isotropy breaking in the density response of hydrogen at large wave numbers.

We present an analysis of the static exchange-correlation (XC) kernel computed from hybrid functionals with a single mixing coefficient such as PBE0 and PBE0-1/3. We break down the hybrid XC kernels into the exchange and correlation parts using the Hartree-Fock functional, the exchange-only PBE, and the correlation-only PBE. This decomposition is combined with exact data for the static XC kernel of the uniform electron gas and an Airy gas model within a subsystem functional approach. This gives us a tool for the non-empirical choice of the mixing coefficient under ambient and extreme conditions. Our analysis provides physical insights into the effect of the variation of the mixing coefficient in hybrid functionals, which is of immense practical value. The presented approach is general and can be used for other types of functionals like screened hybrids.

We combine ab initio path integral Monte Carlo (PIMC) simulations with fixed ion configurations from density functional theory molecular dynamics (DFT-MD) simulations to solve the electronic problem for hydrogen under warm dense matter conditions [Böhme et al., Phys. Rev. Lett. 129, 066402 (2022)0031-900710.1103/PhysRevLett.129.066402]. The problem of path collapse due to the Coulomb attraction is avoided by utilizing the pair approximation, which is compared against the simpler Kelbg pair potential. We find very favorable convergence behavior towards the former. Since we do not impose any nodal restrictions, our PIMC simulations are afflicted with the notorious fermion sign problem, which we analyze in detail. While computationally demanding, our results constitute an exact benchmark for other methods and approximations within DFT. Our setup gives us the unique capability to study important properties of warm dense hydrogen such as the electronic static density response and exchange-correlation kernel without any model assumptions, which will be very valuable for a variety of applications such as the interpretation of experiments and the development of new XC functionals.

The experimental investigation of matter under extreme densities and temperatures, as in astrophysical objects and nuclear fusion applications, constitutes one of the most active frontiers at the interface of material science, plasma physics, and engineering. The central obstacle is given by the rigorous interpretation of the experimental results, as even the diagnosis of basic parameters like the temperature T is rendered difficult at these extreme conditions. Here, we present a simple, approximation-free method to extract the temperature of arbitrarily complex materials in thermal equilibrium from X-ray Thomson scattering experiments, without the need for any simulations or an explicit deconvolution. Our paradigm can be readily implemented at modern facilities and corresponding experiments will have a profound impact on our understanding of warm dense matter and beyond, and open up a variety of appealing possibilities in the context of thermonuclear fusion, laboratory astrophysics, and related disciplines., (© 2022. The Author(s).)

Extreme conditions inside ice giants such as Uranus and Neptune can result in peculiar chemistry and structural transitions, e.g., the precipitation of diamonds or superionic water, as so far experimentally observed only for pure C─H and H 2 O systems, respectively. Here, we investigate a stoichiometric mixture of C and H 2 O by shock-compressing polyethylene terephthalate (PET) plastics and performing in situ x-ray probing. We observe diamond formation at pressures between 72 ± 7 and 125 ± 13 GPa at temperatures ranging from ~3500 to ~6000 K. Combining x-ray diffraction and small-angle x-ray scattering, we access the kinetics of this exotic reaction. The observed demixing of C and H 2 O suggests that diamond precipitation inside the ice giants is enhanced by oxygen, which can lead to isolated water and thus the formation of superionic structures relevant to the planets' magnetic fields. Moreover, our measurements indicate a way of producing nanodiamonds by simple laser-driven shock compression of cheap PET plastics.

The properties of hydrogen under extreme conditions are important for many applications, including inertial confinement fusion and astrophysical models. A key quantity is given by the electronic density response to an external perturbation, which is probed in x-ray Thomson scattering experiments-the state of the art diagnostics from which system parameters like the free electron density n&#95;{e}, the electronic temperature T&#95;{e}, and the charge state Z can be inferred. In this work, we present highly accurate path integral Monte Carlo results for the static electronic density response of hydrogen. We obtain the static exchange-correlation (XC) kernel K&#95;{XC}, which is of central relevance for many applications, such as time-dependent density functional theory. This gives us a first unbiased look into the electronic density response of hydrogen in the warm-dense matter regime, thereby opening up a gamut of avenues for future research.

The rigorous description of correlated quantum many-body systems constitutes one of the most challenging tasks in contemporary physics and related disciplines. In this context, a particularly useful tool is the concept of effective pair potentials that take into account the effects of the complex many-body medium consistently. In this work, we present extensive, highly accurate ab initio path integral Monte Carlo (PIMC) results for the effective interaction and the effective force between two electrons in the presence of the uniform electron gas. This gives us a direct insight into finite-size effects, thereby, opening up the possibility for novel domain decompositions and methodological advances. In addition, we present unassailable numerical proof for an effective attraction between two electrons under moderate coupling conditions, without the mediation of an underlying ionic structure. Finally, we compare our exact PIMC results to effective potentials from linear-response theory, and we demonstrate their usefulness for the description of the dynamic structure factor. All PIMC results are made freely available online and can be used as a thorough benchmark for new developments and approximations.

Ion stopping in warm dense matter is a process of fundamental importance for the understanding of the properties of dense plasmas, the realization and the interpretation of experiments involving ion-beam-heated warm dense matter samples, and for inertial confinement fusion research. The theoretical description of the ion stopping power in warm dense matter is difficult notably due to electron coupling and degeneracy, and measurements are still largely missing. In particular, the low-velocity stopping range, that features the largest modelling uncertainties, remains virtually unexplored. Here, we report proton energy-loss measurements in warm dense plasma at unprecedented low projectile velocities. Our energy-loss data, combined with a precise target characterization based on plasma-emission measurements using two independent spectroscopy diagnostics, demonstrate a significant deviation of the stopping power from classical models in this regime. In particular, we show that our results are in closest agreement with recent first-principles simulations based on time-dependent density functional theory., (© 2022. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

We explore a new formalism to study the nonlinear electronic density response based on Kohn-Sham density functional theory (KS-DFT) at partially and strongly quantum degenerate regimes. It is demonstrated that the KS-DFT calculations are able to accurately reproduce the available path integral Monte Carlo simulation results at temperatures relevant for warm dense matter research. The existing analytical results for the quadratic and cubic response functions are rigorously tested. It is demonstrated that the analytical results for the quadratic response function closely agree with the KS-DFT data. Furthermore, the performed analysis reveals that currently available analytical formulas for the cubic response function are not able to describe simulation results, neither qualitatively nor quantitatively, at small wavenumbers q < 2 q F , with q F being the Fermi wavenumber. The results show that KS-DFT can be used to describe warm dense matter that is strongly perturbed by an external field with remarkable accuracy. Furthermore, it is demonstrated that KS-DFT constitutes a valuable tool to guide the development of the nonlinear response theory of correlated quantum electrons from ambient to extreme conditions. This opens up new avenues to study nonlinear effects in a gamut of different contexts at conditions that cannot be accessed with previously used path integral Monte Carlo methods.