A strong diffusive ion mode in dense ionized matter predicted by Langevin dynamics

Variance-reduced random batch Langevin dynamics

The random batch method is advantageous in accelerating force calculations in particle simulations, but it poses a challenge of removing the artificial heating effect in application to the Langevin dynamics. We develop an approach to solve this issue by estimating the force variance, resulting in a variance-reduced random batch Langevin dynamics. Theoretical analysis shows the high-order local truncation error of the time step in the numerical discretization scheme, consistent with the fluctuation-dissipation theorem. The numerical results indicate that the method can achieve a significant variance reduction since a smaller batch size provides accurate approximation, demonstrating the attractive feature of the variance-reduced random batch method for Langevin dynamics.

Hyperballistic transport in dense systems of charged particles under ac electric fields

The Langevin equation is ubiquitously employed to numerically simulate plasmas, colloids, and electrolytes. However, the usual assumption of white noise becomes untenable when the system is subject to an external ac electric field. This is because the charged particles in the system, which provide the thermal bath for the particle transport, become themselves responsive to the ac field and the thermal noise is field dependent and non-Markovian. We theoretically study the particle diffusivity in a Langevin transport model for a tagged charged particle immersed in a dense system of charged particles (plus also, possibly, other neutral particles) that act as the thermal bath, under an external ac electric field. This is done by properly accounting for the effects of the ac field on the thermal bath statistics. We analytically derive the time-dependent generalized diffusivity D(t) for different initial conditions. The generalized diffusivity exhibits damped oscillatory-like behavior with initial very large peaks, where the generalized diffusion coefficient is enhanced by orders of magnitude with respect to the infinite-time steady-state value. The latter coincides with the Stokes-Einstein diffusivity in the absence of an external field. For initial conditions where the external field is already on at t=0 and the system is thermalized under dc conditions for t≤0, the short-time behavior is hyperballistic, MSD∼t^{4} (where MSD is the mean-squared displacement), leading to giant enhancement of the particle transport. Finally, the theory elucidates the role of medium polarization on the local Lorentz field, and allows for estimates of the effective electric charge due to polarization by the surrounding charges.

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

Disentangling the effects of non-adiabatic interactions upon ion self-diffusion within warm dense hydrogen

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. We implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Temperature relaxation in strongly-coupled binary ionic mixtures

New facilities such as the National Ignition Facility and the Linac Coherent Light Source have pushed the frontiers of high energy-density matter. These facilities offer unprecedented opportunities for exploring extreme states of matter, ranging from cryogenic solid-state systems to hot, dense plasmas, with applications to inertial-confinement fusion and astrophysics. However, significant gaps in our understanding of material properties in these rapidly evolving systems still persist. In particular, non-equilibrium transport properties of strongly-coupled Coulomb systems remain an open question. Here, we study ion-ion temperature relaxation in a binary mixture, exploiting a recently-developed dual-species ultracold neutral plasma. We compare measured relaxation rates with atomistic simulations and a range of popular theories. Our work validates the assumptions and capabilities of the simulations and invalidates theoretical models in this regime. This work illustrates an approach for precision determinations of detailed material properties in Coulomb mixtures across a wide range of conditions.

Influence of dissipation and effective interaction on the dense plasma dynamic structure factor

The ionic dynamic structure factor is examined to assess the relative roles of dissipation and the effective ionic interaction. Two disparate physically based models of dissipation, which can differ numerically by orders of magnitude, are used in molecular dynamics. We find a negligible impact on the amplitudes of the dynamic structure factors for physically realistic parameter values. We then examine the effective ionic interaction by varying its strength, the size of the atomic core (through a pseudopotential), and the screening model. We find that "diffusive" peaks in the dynamic structure factor are very sensitive to the form of the ionic interaction, and this sensitivity arises primarily from atomic physics through the pseudopotential. This suggests that it would be useful to employ the measured zero-frequency dynamic structure factor S_{ii}(k,0) as a constraint on the effective interaction, which in turn can be used to compute physical properties.

An Investigation into the Approximations Used in Wave Packet Molecular Dynamics for the Study of Warm Dense Matter

Wave packet molecular dynamics (WPMD) has recently received a lot of attention as a computationally fast tool with which to study dynamical processes in warm dense matter beyond the Born–Oppenheimer approximation. These techniques, typically, employ many approximations to achieve computational efficiency while implementing semi-empirical scaling parameters to retain accuracy. We investigated three of the main approximations ubiquitous to WPMD: a restricted basis set, approximations to exchange, and the lack of correlation. We examined each of these approximations in regard to atomic and molecular hydrogen in addition to a dense hydrogen plasma. We found that the biggest improvement to WPMD comes from combining a two-Gaussian basis with a semi-empirical correction based on the valence-bond wave function. A single parameter scales this correction to match experimental pressures of dense hydrogen. Ultimately, we found that semi-empirical scaling parameters are necessary to correct for the main approximations in WPMD. However, reducing the scaling parameters for more ab-initio terms gives more accurate results and displays the underlying physics more readily.

High-resolution inelastic x-ray scattering at the high energy density scientific instrument at the European X-Ray Free-Electron Laser

We introduce a setup to measure high-resolution inelastic x-ray scattering at the High Energy Density scientific instrument at the European X-Ray Free-Electron Laser (XFEL). The setup uses the Si (533) reflection in a channel-cut monochromator and three spherical diced analyzer crystals in near-backscattering geometry to reach a high spectral resolution. An energy resolution of 44 meV is demonstrated for the experimental setup, close to the theoretically achievable minimum resolution. The analyzer crystals and detector are mounted on a curved-rail system, allowing quick and reliable changes in scattering angle without breaking vacuum. The entire setup is designed for operation at 10 Hz, the same repetition rate as the high-power lasers available at the instrument and the fundamental repetition rate of the European XFEL. Among other measurements, it is envisioned that this setup will allow studies of the dynamics of highly transient laser generated states of matter.

Probing the Electronic Structure of Warm Dense Nickel via Resonant Inelastic X-Ray Scattering

The development of bright free-electron lasers (FEL) has revolutionized our ability to create and study matter in the high-energy-density (HED) regime. Current diagnostic techniques have been successful in yielding information on fundamental thermodynamic plasma properties, but provide only limited or indirect information on the detailed quantum structure of these systems, and on how it is affected by ionization dynamics. Here we show how the valence electronic structure of solid-density nickel, heated to temperatures of around 10 of eV on femtosecond timescales, can be probed by single-shot resonant inelastic x-ray scattering (RIXS) at the Linac Coherent Light Source FEL. The RIXS spectrum provides a wealth of information on the HED system that goes well beyond what can be extracted from x-ray absorption or emission spectroscopy alone, and is particularly well suited to time-resolved studies of electronic-structure dynamics.

An approach for the measurement of the bulk temperature of single crystal diamond using an X-ray free electron laser

We present a method to determine the bulk temperature of a single crystal diamond sample at an X-Ray free electron laser using inelastic X-ray scattering. The experiment was performed at the high energy density instrument at the European XFEL GmbH, Germany. The technique, based on inelastic X-ray scattering and the principle of detailed balance, was demonstrated to give accurate temperature measurements, within [Formula: see text] for both room temperature diamond and heated diamond to 500 K. Here, the temperature was increased in a controlled way using a resistive heater to test theoretical predictions of the scaling of the signal with temperature. The method was tested by validating the energy of the phonon modes with previous measurements made at room temperature using inelastic X-ray scattering and neutron scattering techniques. This technique could be used to determine the bulk temperature in transient systems with a temporal resolution of 50 fs and for which accurate measurements of thermodynamic properties are vital to build accurate equation of state and transport models.

Ion energy-loss characteristics and friction in a free-electron gas at warm dense matter and nonideal dense plasma conditions

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-Sjölander (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.

Structure and dynamics of warm dense aluminum: a molecular dynamics study with density functional theory and deep potential

We perform a systematic study on the structure and dynamics of warm dense aluminum (Al) at temperatures ranging from 0.5 to 5.0 eV with molecular dynamics utilizing both density functional theory (DFT) and the deep potential (DP) method. On one hand, unlike the Thomas-Fermi kinetic energy density functional (KEDF), we find that the orbital-free DFT method with the Wang-Teter non-local KEDF yields properties of warm dense Al that agree well with the Kohn-Sham DFT method, enabling accurate orbital-free DFT simulations of warm dense Al at relatively low temperatures. On the other hand, the DP method constructs a deep neural network that has a high accuracy in reproducing short- and long-ranged properties of warm dense Al when compared to the DFT methods. The DP method is orders of magnitudes faster than DFT and is well-suited for simulating large systems and long trajectories to yield accurate properties of warm dense Al. Our results suggest that the combination of DFT methods and the DP model is a powerful tool for accurately and efficiently simulating warm dense matter.

Fast nonadiabatic dynamics of many-body quantum systems

Modeling many-body quantum systems with strong interactions is one of the core challenges of modern physics. A range of methods has been developed to approach this task, each with its own idiosyncrasies, approximations, and realm of applicability. However, there remain many problems that are intractable for existing methods. In particular, many approaches face a huge computational barrier when modeling large numbers of coupled electrons and ions at finite temperature. Here, we address this shortfall with a new approach to modeling many-body quantum systems. On the basis of the Bohmian trajectory formalism, our new method treats the full particle dynamics with a considerable increase in computational speed. As a result, we are able to perform large-scale simulations of coupled electron-ion systems without using the adiabatic Born-Oppenheimer approximation.

Ion friction at small values of the Coulomb logarithm

Transport properties of high-energy-density plasmas are influenced by the ion collision rate. Traditionally, this rate involves the Coulomb logarithm, lnΛ. Typical values of lnΛ are ≈10-20 in kinetic theories where transport properties are dominated by weak-scattering events caused by long-range forces. The validity of these theories breaks down for strongly coupled plasmas, when lnΛ is of order one. We present measurements and simulations of collision data in strongly coupled plasmas when lnΛ is small. Experiments are carried out in the first dual-species ultracold neutral plasma (UNP), using Ca^{+} and Yb^{+} ions. We find strong collisional coupling between the different ion species in the bulk of the plasma. We simulate the plasma using a two-species fluid code that includes Coulomb logarithms derived from either a screened Coulomb potential or a the potential of mean force. We find generally good agreement between the experimental measurements and the simulations. With some improvements, the mixed Ca^{+} and Yb^{+} dual-species UNP will be a promising platform for testing theoretical expressions for lnΛ and collision cross-sections from kinetic theories through measurements of energy relaxation, stopping power, two-stream instabilities, and the evolution of sculpted distribution functions in an idealized environment in which the initial temperatures, densities, and charge states are accurately known.

Dynamical structure factor of strongly coupled ions in a dense quantum plasma

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 (nonideality) 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ölander approximation, and it is compared to the MD simulation data obtained considering the electronic screening in the random phase approximation and using the Yukawa potential. 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 nonnegligible 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.

Setup for meV-resolution inelastic X-ray scattering measurements and X-ray diffraction at the Matter in Extreme Conditions endstation at the Linac Coherent Light Source

We describe a setup for performing inelastic X-ray scattering and X-ray diffraction measurements at the Matter in Extreme Conditions (MEC) endstation of the Linac Coherent Light Source. This technique is capable of performing high-, meV-resolution measurements of dynamic ion features in both crystalline and non-crystalline materials. A four-bounce silicon (533) monochromator was used in conjunction with three silicon (533) diced crystal analyzers to provide an energy resolution of ∼50 meV over a range of ∼500 meV in single shot measurements. In addition to the instrument resolution function, we demonstrate the measurement of longitudinal acoustic phonon modes in polycrystalline diamond. Furthermore, this setup may be combined with the high intensity laser drivers available at MEC to create warm dense matter and subsequently measure ion acoustic modes.

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.