Ab initio simulations of the K-edge shift along the aluminum Hugoniot

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

We present a free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin-Kohn-Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo-Greenwood optic data with a single atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo-Greenwood calculations with AIMD snapshots. Kubo-Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomenon of redshift-to-blueshift in K-L (1s→2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s→2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

Time-Resolved XUV Opacity Measurements of Warm Dense Aluminum

The free-free opacity in plasmas is fundamental to our understanding of energy transport in stellar interiors and for inertial confinement fusion research. However, theoretical predictions in the challenging dense plasma regime are conflicting and there is a dearth of accurate experimental data to allow for direct model validation. Here we present time-resolved transmission measurements in solid-density Al heated by an XUV free-electron laser. We use a novel functional optimization approach to extract the temperature-dependent absorption coefficient directly from an oversampled pool of single-shot measurements, and find a pronounced enhancement of the opacity as the plasma is heated to temperatures of order of the Fermi energy. Plasma heating and opacity enhancement are observed on ultrafast timescales, within the duration of the femtosecond XUV pulse. We attribute further rises in the opacity on ps timescales to melt and the formation of warm dense matter.

Path integral Monte Carlo simulations of warm dense aluminum

We perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of 0.1-32.4gcm^{-3} and 10^{4}-10^{8} K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multishock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multishock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and 〈Z〉 to reveal the evolution of the plasma structure and ionization behavior.

Continuum Lowering and Fermi-Surface Rising in Strongly Coupled and Degenerate Plasmas

Continuum lowering is a well known and important physics concept that describes the ionization potential depression (IPD) in plasmas caused by thermal- or pressure-induced ionization of outer-shell electrons. The existing IPD models are often used to characterize plasma conditions and to gauge opacity calculations. Recent precision measurements have revealed deficits in our understanding of continuum lowering in dense hot plasmas. However, these investigations have so far been limited to IPD in strongly coupled but nondegenerate plasmas. Here, we report a first-principles study of the K-edge shifting in both strongly coupled and fully degenerate carbon plasmas, with quantum molecular dynamics calculations based on the all-electron density-functional theory. The resulting K-edge shifting versus plasma density, as a probe to the continuum lowering and the Fermi-surface rising, is found to be significantly different from predictions of existing IPD models. In contrast, a simple model of "single-atom-in-box," developed in this work, accurately predicts K-edge locations as ab initio calculations provide.

Chemical and Morphological Inhomogeneity of Aluminum Metal and Oxides from Soft X-ray Spectromicroscopy

Oxygen and aluminum K-edge X-ray absorption spectroscopy (XAS), imaging from a scanning transmission X-ray microscope (STXM), and first-principles calculations were used to probe the composition and morphology of bulk aluminum metal, α- and γ-Al₂O₃, and several types of aluminum nanoparticles. The imaging results agreed with earlier transmission electron microscopy studies that showed a 2 to 5 nm thick layer of Al₂O₃ on all the Al surfaces. Spectral interpretations were guided by examination of the calculated transition energies, which agreed well with the spectroscopic measurements. Features observed in the experimental O and Al K-edge XAS were used to determine the chemical structure and phase of the Al₂O₃ on the aluminum surfaces. For unprotected 18 and 100 nm Al nanoparticles, this analysis revealed an oxide layer that was similar to γ-Al₂O₃ and comprised of both tetrahedral and octahedral Al coordination sites. For oleic acid-protected Al nanoparticles, only tetrahedral Al oxide coordination sites were observed. The results were correlated to trends in the reactivity of the different materials, which suggests that the structures of different Al₂O₃ layers have an important role in the accessibility of the underlying Al metal toward further oxidation. Combined, the Al K-edge XAS and STXM results provided detailed chemical information that was not obtained from powder X-ray diffraction or imaging from a transmission electron microscope.

High-pressure structural changes in liquid silica

The structural properties of liquid silica at high pressure and moderate temperature conditions, also referred to as the warm dense matter regime, were investigated using time-resolved K-edge x-ray absorption spectroscopy and ab initio calculations. We used a nanosecond laser beam to compress uniformly a solid SiO_{2} target and a picosecond laser beam to generate a broadband x-ray source. We obtained x-ray absorption spectra at the Si K edge over a large pressure-temperature domain to probe the liquid phase up to 3.6 times the normal solid density. Using ab initio simulations, we are able to interpret the changes in the x-ray absorption near-edge structure with increasing densities as an increase in the coordination number of silicon by oxygen atoms from 4 to 9. This indicates that, up to significant temperatures, the liquid structure becomes akin to what is found in the solid SiO_{2} phases.

Equations of state, transport properties, and compositions of argon plasma: combination of self-consistent fluid variation theory and linear response theory

A consistent theoretical model that can be applied in a wide range of densities and temperatures is necessary for understanding the variation of a material's properties during compression and heating. Taking argon as an example, we show that the combination of self-consistent fluid variational theory and linear response theory is a promising route for studying warm dense matter. Following this route, the compositions, equations of state, and transport properties of argon plasma are calculated in a wide range of densities (0.001-20 g/cm(3)) and temperatures (5-100 kK). The obtained equations of state and electrical conductivities are found in good agreement with available experimental data. The plasma phase transition of argon is observed at temperatures below 30 kK and density about 2-6g/cm(3). The minimum density for the metallization of argon is found to be about 5.8 g/cm(3), occurring at 30-40 kK. The effects of many-particle correlations and dynamic screening on the electrical conductivity are also discussed through the effective potentials.

Metallization of warm dense SiO(2) studied by XANES spectroscopy

We investigate the evolution of the electronic structure of fused silica in a dense plasma regime using time-resolved x-ray absorption spectroscopy. We use a nanosecond (ns) laser beam to generate a strong uniform shock wave in the sample and a picosecond (ps) pulse to produce a broadband x-ray source near the Si K edge. By varying the delay between the two laser beams and the intensity of the ns beam, we explore a large thermodynamical domain with densities varying from 1 to 5  g/cm^{3} and temperatures up to 5 eV. In contrast to normal conditions where silica is a well-known insulator with a wide band gap of 8.9 eV, we find that shocked silica exhibits a pseudogap as a semimetal throughout this thermodynamical domain. This is in quantitative agreement with density functional theory predictions performed using the generalized gradient approximation.

Dynamic ionic clusters with flowing electron bubbles from warm to hot dense iron along the Hugoniot curve

Complex structures of warm and hot dense matter are essential to understanding the behavior of materials in high energy density processes and provide new features of matter constitutions. Here, around a new unified first-principles determined Hugoniot curve of iron from the normal condensed condition up to 1 Gbar, the novel structures characterized by the ionic clusters with electron bubbles are found using quantum Langevin molecular dynamics. Subsistence of complex clusters can persist in the time scale of 50 fs dynamically with quantum flowing bubbles, which are produced by the interplay of Fermi electron degeneracy, the ionic coupling, and the dynamical nature. With the inclusion of those complicated features in quantum Langevin molecular dynamics, the present equation of states could serve as a first-principles based database in a wide range of temperatures and densities.

X-ray diagnosis of the pressure induced Mott nonmetal-metal transition

The evolution of the K-edge x-ray absorption near-edge spectroscopy (XANES) spectrum is investigated for an aluminum plasma expanding from the solid density down to 0.5  g/cm{3}, with temperatures lying from 5 down to 2 eV. The dense plasma is generated by nanosecond laser-induced shock compression. These conditions correspond to the density-temperature region where a metal-nonmetal transition occurs as the density decreases. This transition is directly observed in XANES spectra measurements through the progressive formation of a preedge structure for densities around 1.6  g/cm{3}. Ab initio calculations based on density functional theory and a jellium model have been efficiently tested through direct comparison with the experimental measurements and show that this preedge corresponds to the relocalization of the 3p atomic orbital as the system evolves from a dense plasma toward a partially ionized atomic fluid.

Electronic structure investigation of highly compressed aluminum with K edge absorption spectroscopy

The electronic structure evolution of highly compressed aluminum has been investigated using time resolved K edge x-ray absorption spectroscopy. A long laser pulse (500 ps, I(L)≈8×10(13) W/cm(2)) was used to create a uniform shock. A second ps pulse (I(L)≈10(17)  W/cm(2)) generated an ultrashort broadband x-ray source near the Al K edge. The main target was designed to probe aluminum at reshocked conditions up to now unexplored (3 times the solid density and temperatures around 8 eV). The hydrodynamical conditions were obtained using rear side visible diagnostics. Data were compared to ab initio and dense plasma calculations, indicating potential improvements in either description. This comparison shows that x-ray-absorption near-edge structure measurements provide a unique capability to probe matter at these extreme conditions and severally constrains theoretical approaches currently used.

High-temperature high-pressure phases of lithium from electron force field (eFF) quantum electron dynamics simulations

We recently developed the electron force field (eFF) method for practical nonadiabatic electron dynamics simulations of materials under extreme conditions and showed that it gave an excellent description of the shock thermodynamics of hydrogen from molecules to atoms to plasma, as well as the electron dynamics of the Auger decay in diamondoids following core electron ionization. Here we apply eFF to the shock thermodynamics of lithium metal, where we find two distinct consecutive phase changes that manifest themselves as a kink in the shock Hugoniot, previously observed experimentally, but not explained. Analyzing the atomic distribution functions, we establish that the first phase transition corresponds to (i) an fcc-to-cI16 phase transition that was observed previously in diamond anvil cell experiments at low temperature and (ii) a second phase transition that corresponds to the formation of a new amorphous phase (amor) of lithium that is distinct from normal molten lithium. The amorphous phase has enhanced valence electron-nucleus interactions due to localization of electrons into interstitial locations, along with a random connectivity distribution function. This indicates that eFF can characterize and compute the relative stability of states of matter under extreme conditions (e.g., warm dense matter).

Pressure and electrical resistivity measurements on hot expanded nickel: comparisons with quantum molecular dynamics simulations and average atom approaches

We present experimental results on pressure and resistivity on expanded nickel at a density of 0.1 g/cm3 and temperature of a few eV. These data, corresponding to the warm dense matter regime, are used to benchmark different theoretical approaches. A comparison is presented between fully three-dimensional quantum molecular dynamics (QMD) methods, based on density functional theory, with average atom methods, that are essentially one dimensional. In this regime the evaluation of the thermodynamic properties as well as electrical properties is difficult due to the concurrence of density and thermal effects which directly drive the metal-nonmetal transition. Experimental pressures and resistivities are given in a tabular form with temperatures deduced from QMD simulations.

Unified first principles description from warm dense matter to ideal ionized gas plasma: electron-ion collisions induced friction

Electron-ion interactions are central to numerous phenomena in the warm dense matter (WDM) regime and at higher temperature. The electron-ion collisions induced friction at high temperature is introduced in the procedure of ab initio molecular dynamics using the Langevin equation based on density functional theory. In this framework, as a test for Fe and H up to 1000 eV, the equation of state and the transition of electronic structures of the materials with very wide density and temperature can be described, which covers a full range of WDM up to high energy density physics. A unified first principles description from condensed matter to ideal ionized gas plasma is constructed.

Electronic structure of an XUV photogenerated solid-density aluminum plasma

By use of high intensity XUV radiation from the FLASH free-electron laser at DESY, we have created highly excited exotic states of matter in solid-density aluminum samples. The XUV intensity is sufficiently high to excite an inner-shell electron from a large fraction of the atoms in the focal region. We show that soft-x-ray emission spectroscopy measurements reveal the electronic temperature and density of this highly excited system immediately after the excitation pulse, with detailed calculations of the electronic structure, based on finite-temperature density functional theory, in good agreement with the experimental results.

Picosecond short-range disordering in isochorically heated aluminum at solid density

Using ultrafast x-ray probing, we experimentally observed a progressive loss of ordering within solid-density aluminum as the temperature raises from 300 K to >10{4} K. The Al sample was isochorically heated by a short ( approximately ps), laser-accelerated proton beam and probed by a short broadband x-ray source around the Al K edge. The loss of short-range ordering is detected through the progressive smoothing of the time-resolved x-ray absorption near-edge spectroscopy (XANES) structure. The results are compared with two different theoretical models of warm dense matter and allow us to put an upper bound on the onset of ion lattice disorder within the heated solid-density medium of approximately 10 ps.