Publications

Abstract not provided.

Ab intio molecular dynamics simulations show that the electrical conductivity of liquid SiO2 is semimetallic at the conditions of the deep molten mantle of early Earth and super-Earths, raising the possibility of silicate dynamos in these bodies. Whereas the electrical conductivity increases uniformly with increasing temperature, it depends nonmonotonically on compression. At very high pressure, the electrical conductivity decreases on compression, opposite to the behavior of many materials. We show that this behavior is caused by a novel compression mechanism: the development of broken charge ordering, and its influence on the electronic band gap.

In recent years, α-quartz has been used prolifically as an impedance matching standard in shock wave experiments in the multi-Mbar regime (1 Mbar = 100 GPa = 0.1 TPa). This is due to the fact that above ∼90-100 GPa along the principal Hugoniot α-quartz becomes reflective, and thus, shock velocities can be measured to high precision using velocity interferometry. The Hugoniot and release of α-quartz have been studied extensively, enabling the development of an analytical release model for use in impedance matching. However, this analytical release model has only been validated over a range of 300-1200 GPa (0.3-1.2 TPa). Here, we extend this analytical model to 200-3000 GPa (0.2-3 TPa) through additional α-quartz Hugoniot and release measurements, as well as first-principles molecular dynamics calculations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Materials in nuclear and conventional weapons can reach multi-megabar pressures and 1000s of degree temperatures on timescales ranging from microseconds to nanoseconds. Understanding the response of complex materials under these conditions is important for designing and assessing changes to nuclear weapons. In the next few decades, a major concern will be evaluating the behavior of aging materials and remanufactured components. The science to enable the program to underwrite decisions quickly and confidently on use, remanufacturing, and replacement of these materials will be critical to NNSA’s new Stockpile Responsiveness Program. Material response is also important for assessing the risks posed by adversaries or proliferants. Dynamic materials research, which refers to the use of high-speed experiments to produce extreme conditions in matter, is an important part of NNSA’s Stockpile Stewardship Program.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Shock compression experiments in the few hundred GPa (multi-Mbar) regime were performed on Lithium Deuteride single crystals. This study utilized the high velocity flyer plate capability of the Sandia Z Machine to perform impact experiments at flyer plate velocities in the range of 17-32 km/s. Measurements included pressure, density, and temperature between ∼190 and 570 GPa along the Principal Hugoniot - the locus of end states achievable through compression by large amplitude shock waves - as well as pressure and density of reshock states up to ∼920 GPa. The experimental measurements are compared with density functional theory calculations, tabular equation of state models, and legacy nuclear driven results that have been reanalyzed using modern equations of state for the shock wave standards used in the experiments.

1- and 2-D simulations of 1-cm radius, gas-puff liners of Ne, Ar, Kr, and Xe imploding onto a deuterium target are conducted using the discharge parameters for the Zebra (1 MA, 130 ns) driver using the resistive MHD code MACH2. This is an implementation of the Staged Z-pinch concept, in which the target is driven to high-energy-density first by shock compression launched by a diffused azimuthal magnetic field ( J×B force), and then by the adiabatic compression as the liner converges on axis. During the run-in phase, the initial shock heating preheats the deuterium plasma, with a subsequent stable, adiabatic compression heating the target to high energy density. Shock compression of the target coincides with the development of a J×B force at the target/liner interface. Stronger B-field transport and earlier shock compression increases with higher-Z liners, which results in an earlier shock arrival on axis. Delayed shock formation in lower-Z liners yields a relative increase in shock heating, however, the 2-D simulations show an increased target isolation from magneto-Rayleigh-Taylor instability penetration, suggesting that an optimal balance between these two effects is reached in an Ar or Kr liner, rather than with Xe.

Accurate theoretical calculations of the nonlinear elastic response of strong solids (e.g., diamond) constitute a fundamental and important scientific need for understanding the response of such materials and for exploring the potential synthesis and design of novel solids. However, without corresponding experimental data, it is difficult to select between predictions from different theoretical methods. Recently the complete set of third-order elastic constants (TOECs) for diamond was determined experimentally, and the validity of various theoretical approaches to calculate the same may now be assessed. We report on the use of density functional theory (DFT) methods to calculate the six third-order elastic constants of diamond. Two different approaches based on homogeneous deformations were used: (1) an energy-strain fitting approach using a prescribed set of deformations, and (2) a longitudinal stress-strain fitting approach using uniaxial compressive strains along the [100], [110], and [111] directions, together with calculated pressure derivatives of the second-order elastic constants. The latter approach provides a direct comparison to the experimental results. The TOECs calculated using the energy-strain approach differ significantly from the measured TOECs. In contrast, calculations using the longitudinal stress-uniaxial strain approach show good agreement with the measured TOECs and match the experimental values significantly better than the TOECs reported in previous theoretical studies. Our results on diamond have demonstrated that, with proper analysis procedures, first-principles calculations can indeed be used to accurately calculate the TOECs of strong solids.

Abstract not provided.

Abstract not provided.

Experiments on the Sandia Z pulsed-power accelerator have demonstrated the ability to produce warm dense matter (WDM) states with unprecedented uniformity, duration, and size, which are ideal for investigations of fundamental WDM properties. For the first time, space-resolved x-ray Thomson scattering (XRTS) spectra from shocked carbon foams were recorded on Z. The large (>20 MA) electrical current produced by Z was used to launch Al flyer plates up to 25 km/s. The impact of the flyer plate on a CH2 foam target produced a shocked state with an estimated pressure of 0.75 Mbar, density of 0.52 g/cm3, and temperature of 4.3 eV. Both unshocked and shocked portions of the foam target were probed with 6.2 keV x-rays produced by focusing the Z-Beamlet laser onto a nearby Mn foil. The data are composed of three spatially distinct spectra that were simultaneously captured with a single spectrometer with high spectral (4.8 eV) and spatial (190 μm) resolutions. Detailed spectral information from three target locations is provided simultaneously: the incident x-ray source, the scattered signal from unshocked foam, and the scattered signal from shocked foam.

We construct thermodynamic potentials for two superionic phases of water [with body-centered cubic (bcc) and face-centered cubic (fcc) oxygen lattice] using a combination of density functional theory (DFT) and molecular dynamics simulations (MD). For this purpose, a generic expression for the free energy of warm dense matter is developed and parametrized with equation of state data from the DFT-MD simulations. A second central aspect is the accurate determination of the entropy, which is done using an approximate two-phase method based on the frequency spectra of the nuclear motion. The boundary between the bcc superionic phase and the ices VII and X calculated with thermodynamic potentials from DFT-MD is consistent with that directly derived from the simulations. Differences in the physical properties of the bcc and fcc superionic phases and their impact on interior modeling of water-rich giant planets are discussed.

Abstract not provided.

Abstract not provided.