Demonstration of X-ray Thomson scattering as diagnostics for miscibility in warm dense matter

Modeling of warm dense hydrogen via explicit real-time electron dynamics: Dynamic structure factors

We present two methods for computing the dynamic structure factor for warm dense hydrogen without invoking either the Born-Oppenheimer approximation or the Chihara decomposition, by employing a wave-packet description that resolves the electron dynamics during ion evolution. First, a semiclassical method is discussed, which is corrected based on known quantum constraints, and second, a direct computation of the density response function within the molecular dynamics. The wave-packet models are compared to PIMC and DFT-MD for the static and low-frequency behavior. For the high-frequency behavior the models recover the expected behavior in the limits of small and large momentum transfers and show the characteristic flattening of the plasmon dispersion for intermediate momentum transfers due to interactions, in agreement with commonly used models for x-ray Thomson scattering. By modeling the electrons and ions on an equal footing, both the ion and free electron part of the spectrum can now be treated within a single framework where we simultaneously resolve the ion-acoustic and plasmon mode, with a self-consistent description of collisions and screening.

Physisorption Behaviors of Organochlorine Pesticides on the InP3 Monolayer from Theoretical Insight

Dichlorodiphenyltrichloroethane (DDT), hexachlorocyclohexane (BHC), aldrin, and chlordimeform are ubiquitous organochlorine pesticide (OCP) residues in the environment, which pose a great threat to human health and ecosystems due to their high toxicity and easy accumulation. Based on the density functional theory (DFT) calculations, a two-dimensional InP3 monolayer was selected as a sensing material to study the sensitivity detection and adsorption behaviors toward BHC, aldrin, chlordimeform, and DDT. The calculation results show that four pesticide molecules are adsorbed on the InP3 surface by physical interaction. The identified response values (69.1, -43.1%) for DDT and chlordimeform reveal the potential of the InP3 monolayer as a sensing material for the detection of these two pesticides, accompanied by the achievement of cyclic utilization by heating to 498 K. The most satisfactory result is the adsorption of BHC, owing to the admirable sensing response (62.7%) and short recovery time (1.8 s) at room temperature, which makes InP3 a promising pesticide sensor for BHC. However, the InP3 surface is unsuitable for aldrin sensing due to poor response (-1.9%). Our work gives theoretical insight into the good sensitivity and recycling of the InP3 monolayer as a new pesticide sensor to detect DDT, BHC, and chlordimeform, which further broadens the application prospect of the InP3 nanosheet into the sensitive detection of organochlorine pesticides in the ecological environment.

Development of a new quantum trajectory molecular dynamics framework

An extension to the wave packet description of quantum plasmas is presented, where the wave packet can be elongated in arbitrary directions. A generalized Ewald summation is constructed for the wave packet models accounting for long-range Coulomb interactions and fermionic effects are approximated by purpose-built Pauli potentials, self-consistent with the wave packets used. We demonstrate its numerical implementation with good parallel support and close to linear scaling in particle number, used for comparisons with the more common wave packet employing isotropic states. Ground state and thermal properties are compared between the models with differences occurring primarily in the electronic subsystem. Especially, the electrical conductivity of dense hydrogen is investigated where a 15% increase in DC conductivity can be seen in our wave packet model compared with other models. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

X-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula

We study ab initio approaches for calculating x-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula that expresses the inelastic contribution in terms of the dielectric function. We study the electronic dynamic structure factor computed from the Mermin dielectric function using an ab initio electron-ion collision frequency in comparison to computations using a linear-response time-dependent density functional theory (LR-TDDFT) framework for hydrogen and beryllium and investigate the dispersion of free-free and bound-free contributions to the scattering signal. A separate treatment of these contributions, where only the free-free part follows the Mermin dispersion, shows good agreement with LR-TDDFT results for ambient-density beryllium, but breaks down for highly compressed matter where the bound states become pressure ionized. LR-TDDFT is used to reanalyze x-ray Thomson scattering experiments on beryllium demonstrating strong deviations from the plasma conditions inferred with traditional analytic models at small scattering angles.

Assessing the accuracy of hybrid exchange-correlation functionals for the density response of warm dense electrons

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.

Non-empirical Mixing Coefficient for Hybrid XC Functionals from Analysis of the XC Kernel

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.

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/cm^{3} to 16 g/cm^{3} and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. We find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

Nonequilibrium band occupation and optical response of gold after ultrafast XUV excitation

Free electron lasers offer unique properties to study matter in states far from equilibrium as they combine short pulses with a large range of photon energies. In particular, the possibility to excite core states drives new relaxation pathways that, in turn, also change the properties of the optically and chemically active electrons. Here, we present a theoretical model for the dynamics of the nonequilibrium occupation of the different energy bands in solid gold driven by exciting deep core states. The resulting optical response is in excellent agreement with recent measurements and, combined with our model, provides a quantitative benchmark for the description of electron-phonon coupling in strongly driven gold. Focusing on sub-picosecond time scales, we find essential differences between the dynamics induced by XUV and visible light.

Towards a dynamic compression facility at the ESRF

Results of the 2018 commissioning and experimental campaigns of the new High Power Laser Facility on the Energy-dispersive X-ray Absorption Spectroscopy (ED-XAS) beamline ID24 at the ESRF are presented. The front-end of the future laser, delivering 15 J in 10 ns, was interfaced to the beamline. Laser-driven dynamic compression experiments were performed on iron oxides, iron alloys and bismuth probed by online time-resolved XAS.

Measurements of Non-Maxwellian Electron Distribution Functions and Their Effect on Laser Heating

Electron velocity distribution functions driven by inverse bremsstrahlung heating are measured to be non-Maxwellian using a novel angularly resolved Thomson-scattering instrument and the corresponding reduction of electrons at slow velocities results in a ∼40% measured reduction in inverse bremsstrahlung absorption. The distribution functions are measured to be super-Gaussian in the bulk (v/v_{th}<3) and Maxwellian in the tail (v/v_{th}>3) when the laser heating rate dominates over the electron-electron thermalization rate. Simulations with the particle code quartz show the shape of the tail is dictated by the uniformity of the laser heating.

Like dissolves like: A first-principles theory for predicting liquid miscibility and mixture dielectric constant

Liquid mixtures are ubiquitous. Miscibility and dielectric constant are fundamental properties that govern the applications of liquid mixtures. However, despite their importance, miscibility is usually predicted qualitatively based on the vaguely defined polarity of the liquids, and the dielectric constant of the mixture is modeled by introducing mixing rules. Here, we develop a first-principles theory for polar liquid mixtures using a statistical field approach, without resorting to mixing rules. With this theory, we obtain simple expressions for the mixture's dielectric constant and free energy of mixing. The dielectric constant predicted by this theory agrees well with measured data for simple binary mixtures. On the basis of the derived free energy of mixing, we can construct a miscibility map in the parameter space of the dielectric constant and molar volume for each liquid. The predicted miscibility shows remarkable agreement with known data, thus providing a quantitative basis for the empirical "like-dissolves-like" rule.

Effective Static Approximation: A Fast and Reliable Tool for Warm-Dense Matter Theory

We present an effective static approximation (ESA) to the local field correction (LFC) of the electron gas that enables highly accurate calculations of electronic properties like the dynamic structure factor S(q,ω), the static structure factor S(q), and the interaction energy v. The ESA combines the recent neural-net representation by T. Dornheim et al., [J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013] of the temperature-dependent LFC in the exact static limit with a consistent large wave-number limit obtained from quantum Monte Carlo data of the on-top pair distribution function g(0). It is suited for a straightforward integration into existing codes. We demonstrate the importance of the LFC for practical applications by reevaluating the results of the recent x-ray Thomson scattering experiment on aluminum by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.115001]. We find that an accurate incorporation of electronic correlations in terms of the ESA leads to a different prediction of the inelastic scattering spectrum than obtained from state-of-the-art models like the Mermin approach or linear-response time-dependent density functional theory. Furthermore, the ESA scheme is particularly relevant for the development of advanced exchange-correlation functionals in density functional theory.