Metallization of Shock-Compressed Liquid Ammonia

Release dynamics of nanodiamonds created by laser-driven shock-compression of polyethylene terephthalate

Laser-driven dynamic compression experiments of plastic materials have found surprisingly fast formation of nanodiamonds (ND) via X-ray probing. This mechanism is relevant for planetary models, but could also open efficient synthesis routes for tailored NDs. We investigate the release mechanics of compressed NDs by molecular dynamics simulation of the isotropic expansion of finite size diamond from different P-T states. Analysing the structural integrity along different release paths via molecular dynamic simulations, we found substantial disintegration rates upon shock release, increasing with the on-Hugnoiot shock temperature. We also find that recrystallization can occur after the expansion and hence during the release, depending on subsequent cooling mechanisms. Our study suggests higher ND recovery rates from off-Hugoniot states, e.g., via double-shocks, due to faster cooling. Laser-driven shock compression experiments of polyethylene terephthalate (PET) samples with in situ X-ray probing at the simulated conditions found diamond signal that persists up to 11 ns after breakout. In the diffraction pattern, we observed peak shifts, which we attribute to thermal expansion of the NDs and thus a total release of pressure, which indicates the stability of the released NDs.

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.

Imaging velocity interferometer system for any reflector (VISAR) diagnostics for high energy density sciences

Two variants of optical imaging velocimetry, specifically the one-dimensional streaked line-imaging and the two-dimensional time-resolved area-imaging versions of the Velocity Interferometer System for Any Reflector (VISAR), have become important diagnostics in high energy density sciences, including inertial confinement fusion and dynamic compression of condensed matter. Here, we give a brief review of the historical development of these techniques, then describe the current implementations at major high energy density (HED) facilities worldwide, including the OMEGA Laser Facility and the National Ignition Facility. We illustrate the versatility and power of these techniques by reviewing diverse applications of imaging VISARs for gas-gun and laser-driven dynamic compression experiments for materials science, shock physics, condensed matter physics, chemical physics, plasma physics, planetary science and astronomy, as well as a broad range of HED experiments and laser-driven inertial confinement fusion research.

High-P-T impedance measurements using a laser-heated diamond anvil cell

The electrical conductivity (EC) of minerals found on Earth and throughout the solar system is a fundamental transport property that is used to understand various dynamical phenomena in planetary interiors. High-pressure and high-temperature (P-T) EC measurements are also an important tool for observing phase transitions. Impedance measurements can accurately measure the EC of a nonmetallic sample. In previous measurements under static conditions using a laser-heated diamond-anvil cell (LHDAC), only direct current resistance is measured, but this method overestimates the bulk sample resistance. Moreover, the previous methodology could only be applied to nontransparent samples in an LHDAC using infrared lasers, limiting the range of measurable composition. To the best of our knowledge, no in situ high-P-T EC measurements of transparent materials have been reported using LHDAC techniques. We developed a novel impedance measurement technique under high-P-T conditions in an LHDAC that applies to transparent samples. As a validation, we measured the EC of Mg_0.9Fe_0.1SiO₃ bridgmanite up to 51 GPa and 2000 K and found that the results are consistent with those of previous studies. We also measured the EC values of sodium chloride to compare with those of previous studies, as well as those of cubic boron nitride and zirconia cement to quantify how well they insulate under high P-T conditions. This is the first report of the impedance and EC measurements of transparent minerals in an LHDAC, which allows the measurement of Fe-poor/-free materials, including the major constituents of the interiors of gas giants and icy planets, under extreme conditions.

Continuous Kubo-Greenwood formula: Theory and numerical implementation

In this paper, we present the so-called continuous Kubo-Greenwood formula intended for the numerical calculation of the dynamic Onsager coefficients and, in particular, the real part of dynamic electrical conductivity. In contrast to the usual Kubo-Greenwood formula, which contains the summation over a discrete set of transitions between electron energy levels, the continuous one is formulated as an integral over the whole energy range. This integral includes the continuous functions: the smoothed squares of matrix elements, D(ɛ,ɛ+ℏω), the densities of state, g(ɛ)g(ɛ+ℏω), and the difference of the Fermi weights, [f(ɛ)-f(ɛ+ℏω)]/(ℏω). The function D(ɛ,ɛ+ℏω) is obtained via the specially developed smoothing procedure. From the theoretical point of view, the continuous formula is an alternative to the usual one. Both can be used to calculate matter properties and produce close results. However, the continuous formula includes the smooth functions that can be plotted and examined. Thus, we can analyze the contributions of various parts of the electron spectrum to the obtained properties. The possibility of such analysis is the main advantage of the continuous formula. The continuous Kubo-Greenwood formula is implemented in the parallel code cubogram. Using the code we demonstrate the influence of technical parameters on the simulation results for liquid aluminum. We also analyze various methods of matrix elements computation and their effect on dynamic electrical conductivity.

Metallic Aluminum Suboxides with Ultrahigh Electrical Conductivity at High Pressure

Aluminum, as the most abundant metallic elemental content in the Earth's crust, usually exists in the form of alumina (Al₂O₃). However, the oxidation state of aluminum and the crystal structures of aluminum oxides in the pressure range of planetary interiors are not well established. Here, we predicted two aluminum suboxides (Al₂O, AlO) and two superoxides (Al₄O₇, AlO₃) with uncommon stoichiometries at high pressures using first-principle calculations and crystal structure prediction methods. We find that the P4/nmm Al₂O becomes stable above ~765 GPa and may survive in the deep mantles or cores of giant planets such as Neptune. Interestingly, the Al₂O and AlO are metallic and have electride features, in which some electrons are localized in the interstitials between atoms. We find that Al₂O has an electrical conductivity one order of magnitude higher than that of iron under the same pressure-temperature conditions, which may influence the total conductivity of giant planets. Our findings enrich the high-pressure phase diagram of aluminum oxides and improve our understanding of the interior structure of giant planets.

Measuring the structure and equation of state of polyethylene terephthalate at megabar pressures

We present structure and equation of state (EOS) measurements of biaxially orientated polyethylene terephthalate (PET, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({\hbox {C}}_{10} {\hbox {H}}_8 {\hbox {O}}_4)_n$$\end{document}(C10H8O4)n, also called mylar) shock-compressed to (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$155 \pm 20$$\end{document}155±20) GPa and (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6000 \pm 1000$$\end{document}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.