Observation of ultrafast nonequilibrium collective dynamics in warm dense hydrogen

Influence of the two-temperature effect on ionization potential depression in hot dense plasma

In hot dense plasma, the interaction between charged particles leads to the ionization potential depression (IPD), which further affects the physical properties of plasma, such as opacity and equation of state. The experiment of IPD of solid-density Al plasma has indicated that present theoretical models cannot give reasonable descriptions of the IPD in hot dense plasma. So, reasonable theoretical methods are needed to describe the effects of hot dense environments on IPD, and the process of generating hot dense plasmas through the interaction between ultrashort laser pulse and solid-density matter also needs to be carefully considered. In the manuscript, two kinds of temperatures for ions and electrons are considered when we compute the ionization potential depression in hot dense plasma. And, the influences of the hot dense environments are included in the electronic structure calculations by using the modified flexible atomic code (FAC), which has included the screening of free electrons and the correlation of ions by correlation functions obtained from the hyper-netted chain (HNC) approximation. A self-consistent-field method is used to calculate the electronic structures. Based on the calculations, the IPD is obtained through the two-step model. Considering the interaction of the femtosecond laser on the solid-density Al plasma of Ciricosta's experiment, we use the two-temperature model to calculate the IPD in nonlocal thermodynamic equilibrium, and the theoretical results are in good agreement with the experimental results. In addition, we also calculated the electron collision ionization cross section and compared it with the results from other models.

Transmissions of an x-ray free electron laser pulse through Al: Influence of nonequilibrium electron kinetics

A theoretical model for investigating the radiative transfer of an x-ray free electron laser (XFEL) pulse is developed based on a one-dimensional radiative transfer equation. The population dynamics of energy levels is obtained by rate equation approximation coupling with the Fokker-Planck equation, in which the electron energy distribution function (EEDF) is self-consistently determined. As an illustrative example, XFEL pulse propagation through a solid-density aluminum (Al) is investigated. The characteristics of the temporal evolution of the x-ray pulse shape, level population, and EEDF are demonstrated. The EEDF usually has two parts in XFEL-Al interactions: the near equilibrium part in the lower energy regions and the nonequilibrium part in the higher energy region. The deep gap between the two parts is quickly filled in the solid-density Al plasma. The pulse shape is distorted and the duration shortens as the x-ray pulse propagates through the Al sample. The x-ray transmission spectra were compared with experimental and other theoretical results, and good agreement was found. There are slight discrepancies between the transmission obtained by solving the Fokker-Planck equation and Maxwellian assumptions because nonequilibrium electrons in the higher energy region account for only a small fraction of the total electrons.

Electron-electron scattering limits thermal conductivity of metals under extremely high electron temperatures

We report on the thermal transport properties of noble metals (gold, silver and copper) under conditions of extremely high electron temperatures (that are on the order of the Fermi energy). We perform parameter-free density functional theory calculations of the electron temperature-dependent electron-phonon coupling, electronic heat capacities, and thermal conductivities to elucidate the strong role played by the excitation of the low lyingd-bands on the transport properties of the noble metals. Our calculations show that, although the three metals have similar electronic band structures, the changes in their electron-phonon coupling at elevated electron temperatures are drastically different; while electron-phonon coupling decreases in gold, it increases in copper and, it remains relatively unperturbed for silver with increasing electron temperatures of up to ∼60 000 K (or 5 eV). We attribute this to the varying contributions from acoustic and longitudinal phonon modes to the electron-phonon coupling in the three metals. Although their electron-phonon coupling changes with electron temperature, the thermal conductivity trends with electron temperature are similar for all three metals. For instance, the thermal conductivities for all three metals reach their maximum values (on par with the room-temperature values of some of the most thermally conductive semiconductors) at electron temperatures of ∼6000 K, and thereafter monotonically decrease due to the enhanced effect of electron-electron scattering for electronic states that are further away from the Fermi energy. As such, only accounting for electron-phonon coupling and neglecting electron-electron scattering can lead to large over-predictions of the thermal conductivities at extremely high electron temperatures. Our results shed light on the microscopic understanding of the electronic scattering mechanisms and thermal transport in noble metals under conditions of extremely high electron temperatures and, as such, are significant for a plethora of applications such as in plasmonic devices that routinely leverage hot electron transport.

Direct Observation of Enhanced Electron-Phonon Coupling in Copper Nanoparticles in the Warm-Dense Matter Regime

Warm dense matter (WDM) represents a highly excited state that lies at the intersection of solids, plasmas, and liquids and that cannot be described by equilibrium theories. The transient nature of this state when created in a laboratory, as well as the difficulties in probing the strongly coupled interactions between the electrons and the ions, make it challenging to develop a complete understanding of matter in this regime. In this work, by exciting isolated ∼8  nm copper nanoparticles with a femtosecond laser below the ablation threshold, we create uniformly excited WDM. Using photoelectron spectroscopy, we measure the instantaneous electron temperature and extract the electron-ion coupling of the nanoparticle as it undergoes a solid-to-WDM phase transition. By comparing with state-of-the-art theories, we confirm that the superheated nanoparticles lie at the boundary between hot solids and plasmas, with associated strong electron-ion coupling. This is evidenced both by a fast energy loss of electrons to ions, and a strong modulation of the electron temperature induced by strong acoustic breathing modes that change the nanoparticle volume. This work demonstrates a new route for experimental exploration of the exotic properties of WDM.

Dynamic and transient processes in warm dense matter

In this paper, we discuss some of the key challenges in the study of time-dependent processes and non-equilibrium behaviour in warm dense matter. We outline some of the basic physics concepts that have underpinned the definition of warm dense matter as a subject area in its own right and then cover, in a selective, non-comprehensive manner, some of the current challenges, pointing along the way to topics covered by the papers presented in this volume. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Application of Boltzmann kinetic equations to model X-ray-created warm dense matter and plasma

In this review, we describe the application of Boltzmann kinetic equations for modelling warm dense matter and plasma formed after irradiation of solid materials with intense femtosecond X-ray pulses. Classical Boltzmann kinetic equations are derived from the reduced N-particle Liouville equations. They include only single-particle densities of ions and free electrons present in the sample. The first version of the Boltzmann kinetic equation solver was completed in 2006. It could model non-equilibrium evolution of X-ray-irradiated finite-size atomic systems. In 2016, the code was adapted to study plasma created from X-ray-irradiated materials. Additional extension of the code was then also performed, enabling simulations in the hard X-ray irradiation regime. In order to avoid treatment of a very high number of active atomic configurations involved in the excitation and relaxation of X-ray-irradiated materials, an approach called 'predominant excitation and relaxation path' (PERP) was introduced. It limited the number of active atomic configurations by following the sample evolution only along most PERPs. The performance of the Boltzmann code is illustrated in the examples of X-ray-heated solid carbon and gold. Actual model limitations and further model developments are discussed. 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.

Plasma environmental effects in the atomic structure for simulating x-ray free-electron-laser-heated solid-density matter

High energy density (HED) matter exists extensively in the Universe, and it can be created with extreme conditions in laboratory facilities such as x-ray free-electron lasers (XFEL). In HED matter, the electronic structure of individual atomic ions is influenced by a dense plasma environment, and one of the most significant phenomena is the ionization potential depression (IPD). Incorporation of the IPD effects is of great importance in accurate modeling of dense plasmas. All theoretical treatments of IPD so far have been based on the assumption of local thermodynamic equilibrium, but its validity is questionable in ultrafast formation dynamics of dense plasmas, particularly when interacting with intense XFEL pulses. A treatment of transient IPD, based on an electronic-structure calculation of an atom in the presence of a plasma environment described by classical particles, has recently been proposed [Phys. Rev. E 103, 023203 (2021)2470-004510.1103/PhysRevE.103.023203], but its application to and impact on plasma dynamics simulations have not been investigated yet. In this work, we extend XMDYN, a hybrid quantum-classical approach combining Monte Carlo and molecular dynamics, by incorporating the proposed IPD treatment into plasma dynamics simulations. We demonstrate the importance of the IPD effects in theoretical modeling of aluminum dense plasmas by comparing two XMDYN simulations: one with electronic-structure calculations of isolated atoms (without IPD) and the other with those of atoms embedded in a plasma (with IPD). At equilibrium, the mean charge obtained in the plasma simulation with IPD is in good agreement with the full quantum-mechanical average-atom model. The present approach promises to be a reliable tool to simulate the creation and nonequilibrium evolution of dense plasmas induced by ultraintense and ultrashort XFEL pulses.

Density functional tight binding approach utilized to study X-ray-induced transitions in solid materials

Intense X-ray pulses from free-electron lasers can trigger ultrafast electronic, structural and magnetic transitions in solid materials, within a material volume which can be precisely shaped through adjustment of X-ray beam parameters. This opens unique prospects for material processing with X rays. However, any fundamental and applicational studies are in need of computational tools, able to predict material response to X-ray radiation. Here we present a dedicated computational approach developed to study X-ray induced transitions in a broad range of solid materials, including those of high chemical complexity. The latter becomes possible due to the implementation of the versatile density functional tight binding code DFTB+ to follow band structure evolution in irradiated materials. The outstanding performance of the implementation is demonstrated with a comparative study of XUV induced graphitization in diamond.

Local field correction to ionization potential depression of ions in warm or hot dense matter

An analytical self-consistent approach was recently established to predict the ionization potential depression (IPD) in multicomponent dense plasmas, which is achieved by considering the self-energy of ions and electrons within the quantum statistical theory. In order to explicitly account for the exchange-correlation effect of electrons, we incorporate the effective static approximation of local field correction (LFC) within our IPD framework through the connection of dynamical structure factor. The effective static approximation poses an accurate description for the asymptotic large wave number behavior with the recently developed machine learning representation of static LFC induced from the path-integral Monte Carlo data. Our calculation shows that the introduction of static LFC through dynamical structure factor brings a nontrivial influence on IPD at warm/hot dense matter conditions. The correlation effect within static LFC could provide up to 20% correction to free-electron contribution of IPD in the strong coupling and degeneracy regime. Furthermore, a new screening factor is obtained from the density distribution of free electrons calculated within the average-atom model, with which excellent agreements are observed with other methods and experiments at warm/hot dense matter conditions.

First-Principles Determination of Electron-Ion Couplings in the Warm Dense Matter Regime

We present first-principles calculations of the rate of energy exchanges between electrons and ions in nonequilibrium warm dense plasmas, liquid metals, and hot solids, a fundamental property for which various models offer diverging predictions. To this end, a Kubo relation for the electron-ion coupling parameter is introduced, which includes self-consistently the quantum, thermal, nonlinear, and strong coupling effects that coexist in materials at the confluence of solids and plasmas. Most importantly, like other Kubo relations widely used for calculating electronic conductivities, the expression can be evaluated using quantum molecular dynamics simulations. Results are presented and compared to experimental and theoretical predictions for representative materials of various electronic complexity, including aluminum, copper, iron, and nickel.

A high resolution extreme ultraviolet spectrometer system optimized for harmonic spectroscopy and XUV beam analysis

We present a modular extreme ultraviolet (XUV) spectrometer system optimized for a broad spectral range of 12-41 nm (30-99 eV) with a high spectral resolution of λ/Δλ ≳ 784 ± 89. The spectrometer system has several operation modes for (1) XUV beam inspection, (2) angular spectral analysis, and (3) imaging spectroscopy. These options allow for a versatile use in high harmonic spectroscopy and XUV beam analysis. The high performance of the spectrometer is demonstrated using a novel cross-sectional imaging method called XUV coherence tomography.

A sensitive EUV Schwarzschild microscope for plasma studies with sub-micrometer resolution

We present an extreme ultraviolet (EUV) microscope using a Schwarzschild objective which is optimized for single-shot sub-micrometer imaging of laser-plasma targets. The microscope has been designed and constructed for imaging the scattering from an EUV-heated solid-density hydrogen jet. Imaging of a cryogenic hydrogen target was demonstrated using single pulses of the free-electron laser in Hamburg (FLASH) free-electron laser at a wavelength of 13.5 nm. In a single exposure, we observe a hydrogen jet with ice fragments with a spatial resolution in the sub-micrometer range. In situ EUV imaging is expected to enable novel experimental capabilities for warm dense matter studies of micrometer-sized samples in laser-plasma experiments.

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.

Short-wavelength free-electron laser sources and science: a review

This review is focused on free-electron lasers (FELs) in the hard to soft x-ray regime. The aim is to provide newcomers to the area with insights into: the basic physics of FELs, the qualities of the radiation they produce, the challenges of transmitting that radiation to end users and the diversity of current scientific applications. Initial consideration is given to FEL theory in order to provide the foundation for discussion of FEL output properties and the technical challenges of short-wavelength FELs. This is followed by an overview of existing x-ray FEL facilities, future facilities and FEL frontiers. To provide a context for information in the above sections, a detailed comparison of the photon pulse characteristics of FEL sources with those of other sources of high brightness x-rays is made. A brief summary of FEL beamline design and photon diagnostics then precedes an overview of FEL scientific applications. Recent highlights are covered in sections on structural biology, atomic and molecular physics, photochemistry, non-linear spectroscopy, shock physics, solid density plasmas. A short industrial perspective is also included to emphasise potential in this area.

Kinetic Boltzmann approach adapted for modeling highly ionized matter created by x-ray irradiation of a solid

We report on the kinetic Boltzmann approach adapted for simulations of highly ionized matter created from a solid by its x-ray irradiation. X rays can excite inner-shell electrons, which leads to the creation of deeply lying core holes. Their relaxation, especially in heavier elements, can take complicated paths, leading to a large number of active configurations. Their number can be so large that solving the set of respective evolution equations becomes computationally inefficient and another modeling approach should be used instead. To circumvent this complexity, the commonly used continuum models employ a superconfiguration scheme. Here, we propose an alternative approach which still uses "true" atomic configurations but limits their number by restricting the sample relaxation to the predominant relaxation paths. We test its reliability, performing respective calculations for a bulk material consisting of light atoms and comparing the results with a full calculation including all relaxation paths. Prospective application for heavy elements is discussed.

Theory of Thomson scattering in inhomogeneous media

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

Free-electron X-ray laser measurements of collisional-damped plasmons in isochorically heated warm dense matter

We present the first highly resolved measurements of the plasmon spectrum in an ultrafast heated solid. Multi-keV x-ray photons from the Linac Coherent Light Source have been focused to one micrometer diameter focal spots producing solid density aluminum plasmas with a known electron density of n_{e}=1.8×10^{23}  cm^{-3}. Detailed balance is observed through the intensity ratio of up- and down-shifted plasmons in x-ray forward scattering spectra measuring the electron temperature. The plasmon damping is treated by electron-ion collision models beyond the Born approximation to determine the electrical conductivity of warm dense aluminum.

Ab initio approach to model x-ray diffraction in warm dense matter

It is demonstrated how the static electron-electron structure factor in warm dense matter can be obtained from density functional theory in combination with quantum Monte Carlo data. In contrast to theories assuming well-separated bound and free states, this ab initio approach yields also valid results for systems close to the Mott transition (pressure ionization), where bound states are strongly modified and merge with the continuum. The approach is applied to x-ray Thomson scattering and compared to predictions of the Chihara formula whereby we use the ion-ion and electron-ion structure from the same simulations. The results show significant deviations of the screening cloud from the often applied Debye-like form.

Equilibration dynamics and conductivity of warm dense hydrogen

We investigate subpicosecond dynamics of warm dense hydrogen at the XUV free-electron laser facility (FLASH) at DESY (Hamburg). Ultrafast impulsive electron heating is initiated by a ≤ 300-fs short x-ray burst of 92-eV photon energy. A second pulse probes the sample via x-ray scattering at jitter-free variable time delay. We show that the initial molecular structure dissociates within (0.9 ± 0.2) ps, allowing us to infer the energy transfer rate between electrons and ions. We evaluate Saha and Thomas-Fermi ionization models in radiation hydrodynamics simulations, predicting plasma parameters that are subsequently used to calculate the static structure factor. A conductivity model for partially ionized plasma is validated by two-temperature density-functional theory coupled to molecular dynamic simulations and agrees with the experimental data. Our results provide important insights and the needed experimental data on transport properties of dense plasmas.

Resolving ultrafast heating of dense cryogenic hydrogen

We report on the dynamics of ultrafast heating in cryogenic hydrogen initiated by a ≲300  fs, 92 eV free electron laser x-ray burst. The rise of the x-ray scattering amplitude from a second x-ray pulse probes the transition from dense cryogenic molecular hydrogen to a nearly uncorrelated plasmalike structure, indicating an electron-ion equilibration time of ∼0.9  ps. The rise time agrees with radiation hydrodynamics simulations based on a conductivity model for partially ionized plasma that is validated by two-temperature density-functional theory.

Reduced coupled-mode approach to electron-ion energy relaxation

We present a reduced model for the energy transfer via coupled collective modes in two-temperature plasmas based on quantum statistical theory. The model is compared with exact numerical evaluations of the coupled-mode (CM) energy transfer rate and with alternative reduced approaches over a range of conditions in the warm dense matter (WDM) and inertial confinement fusion (ICF) regimes. Our approach shows excellent agreement with an exact treatment of the CM rate and supports the importance of the coupled-mode effect for the temperature and energy relaxation in WDM and ICF plasmas. We find that electronic damping of collective ion density fluctuations is crucial for correctly describing the mode spectrum and, thus, the energy exchange. The reduced CM approach is studied over a wide parameter space, enabling us to establish its limits of applicability.

Nonequilibrium electron dynamics in materials driven by high-intensity x-ray pulses

We calculated the evolution of the electron system in solid-density matter irradiated by high-intensity x-ray pulses between 2 and 8 keV using molecular dynamics. For pulses shorter than 40 fs, the kinetic energy distribution of the electrons is highly nonthermal during and right after the pulse, and a large fraction of the absorbed x-ray energy resides with the fast photoelectrons which equilibrate on the timescale of the pulse length. The average ionization and electron temperature of the bulk of the electrons are significantly lower than their equilibrium values.

Determination of the ion temperature in a stainless steel slab exposed to intense ultrashort laser pulses

We present an effective approach to determine the amount of energy absorbed by solid samples exposed to ultrashort laser pulses, thus, retrieving the maximum temperature attained by the ion lattice in the picosecond time scale. The method is based on the pyrometric detection of a slow temperature fluctuation on the rear side of a sample slab associated with absorption of the laser pulse on the front side. This approach, successfully corroborated by theoretical calculations, can provide a robust and practical diagnostic tool for characterization of laser-generated warm dense matter.

Testing quantum mechanics in non-Minkowski space-time with high power lasers and 4(th) generation light sources

A common misperception of quantum gravity is that it requires accessing energies up to the Planck scale of 10¹⁹ GeV, which is unattainable from any conceivable particle collider. Thanks to the development of ultra-high intensity optical lasers, very large accelerations can be now the reached at their focal spot, thus mimicking, by virtue of the equivalence principle, a non Minkowski space-time. Here we derive a semiclassical extension of quantum mechanics that applies to different metrics, but under the assumption of weak gravity. We use our results to show that Thomson scattering of photons by uniformly accelerated electrons predicts an observable effect depending upon acceleration and local metric. In the laboratory frame, a broadening of the Thomson scattered x ray light from a fourth generation light source can be used to detect the modification of the metric associated to electrons accelerated in the field of a high power optical laser.

Ultrafast transitions from solid to liquid and plasma states of graphite induced by x-ray free-electron laser pulses

We used photon pulses from an x-ray free-electron laser to study ultrafast x-ray-induced transitions of graphite from solid to liquid and plasma states. This was accomplished by isochoric heating of graphite samples and simultaneous probing via Bragg and diffuse scattering at high time resolution. We observe that disintegration of the crystal lattice and ion heating of up to 5 eV occur within tens of femtoseconds. The threshold fluence for Bragg-peak degradation is smaller and the ion-heating rate is faster than current x-ray-matter interaction models predict.

Analysis of Thomson scattering from nonequilibrium plasmas

We develop the theory for light scattering as a diagnostic method for plasmas in nonequilibrium states. We show how well-known nonequilibrium features, like beam acoustic modes, arise in the spectra. The analysis of an experiment with strongly driven electrons demonstrates the abilities of the new approach; we find qualitatively different scattering spectra for different times and excellent agreement with the experimental data after time integration. Finally, an analysis of data from dense beryllium suggests that an energetic electron component exists in this experiment as well.

Extent of validity of the hydrodynamic description of ions in dense plasmas

We show that the hydrodynamic description can be applied to modeling the ionic response in dense plasmas for a wide range of length scales that are experimentally accessible. Using numerical simulations for the Yukawa model, we find that the maximum wave number k(max) at which the hydrodynamic description applies is independent of the coupling strength, given by k(max)λ(s)≃0.43, where λ(s) is the ionic screening length. Our results show that the hydrodynamic description can be used for interpreting x-ray scattering data from fourth generation light sources and high power lasers. In addition, our investigation sheds new light on how the domain of validity of the hydrodynamic description depends on both the microscopic properties and the thermodynamic state of fluids in general.

Thomson scattering on inhomogeneous targets

The introduction of brilliant free-electron lasers enables new pump-probe experiments to characterize warm dense matter states. For instance, a short-pulse optical laser irradiates a liquid hydrogen jet that is subsequently probed with brilliant soft x-ray radiation. The strongly inhomogeneous plasma prepared by the optical laser is characterized with particle-in-cell simulations. The interaction of the soft x-ray probe radiation for different time delays between pump and probe with the inhomogeneous plasma is also taken into account via radiative hydrodynamic simulations. We calculate the respective scattering spectrum based on the Born-Mermin approximation for the dynamic structure factor considering the full density and temperature-dependent Thomson scattering cross section throughout the target. We can identify plasmon modes that are generated in different target regions and monitor their temporal evolution. Therefore, such pump-probe experiments are promising tools not only to measure the important plasma parameters density and temperature but also to gain valuable information about their time-dependent profile through the target. The method described here can be applied to various pump-probe scenarios by combining optical lasers and soft x ray, as well as x-ray sources.