Rudolph J. Magyar

X-ray Thomson Scattering in Warm Dense Matter without the Chihara Decomposition.

X-ray Thomson scattering is an important experimental technique used to measure the temperature, ionization state, structure, and density of warm dense matter (WDM). The fundamental property probed in these experiments is the electronic dynamic structure factor. In most models, this is decomposed into three terms [J. Chihara, J. Phys. F 17, 295 (1987)] representing the response of tightly bound, loosely bound, and free electrons. Accompanying this decomposition is the classification of electrons as either bound or free, which is useful for gapped and cold systems but becomes increasingly questionable as temperatures and pressures increase into the WDM regime. In this work we provide unambiguous first principles calculations of the dynamic structure factor of warm dense beryllium, independent of the Chihara form, by treating bound and free states under a single formalism. The computational approach is real-time finite-temperature time-dependent density functional theory (TDDFT) being applied here for the first time to WDM. We compare results from TDDFT to Chihara-based calculations for experimentally relevant conditions in shock-compressed beryllium.

A two-temperature model of radiation damage in α-quartz

Two-temperature models are used to represent the physics of the interaction between atoms and electrons during thermal transients such as radiation damage, laser heating, and cascade simulations. We introduce a two-temperature model applied to an insulator, α-quartz, to model heat deposition in a SiO(2) lattice. Our model of the SiO(2) electronic subsystem is based on quantum simulations of the electronic response in a SiO(2) repeat cell. We observe how the parametrization of the electronic subsystem impacts the degree of permanent amorphization of the lattice, especially compared to a metallic electronic subsystem. The parametrization of the insulator electronic subsystem has a significant effect on the amount of residual defects in the crystal after 10 ps. While recognizing that more development in the application of two-temperature models to insulators is needed, we argue that the inclusion of a simple electronic subsystem substantially improves the realism of such radiation damage simulations.

Ethane-Xenon Mixtures under Shock Conditions.

Mixtures of light elements with heavy elements are important in inertial confinement fusion. We explore the physics of molecular scale mixing through a validation study of equation of state (EOS) properties. Density functional theory molecular dynamics (DFT-MD) at elevated temperature and pressure is used to obtain the thermodynamic state properties of pure xenon, ethane, and various compressed mixture compositions along their principal Hugoniots. In order to validate these simulations, we have performed shock compression experiments using the Sandia Z-Machine. A bond tracking analysis correlates the sharp rise in the Hugoniot curve with the completion of dissociation in ethane. Furthermore, the DFT-based simulation results compare well with the experimental data along the principal Hugoniots and are used to provide insight into the dissociation and temperature along the Hugoniots as a function of mixture composition. Interestingly, we find that the compression ratio for complete dissociation is similar for several compositions suggesting a limiting compression for C-C bonded systems.

Mixing of equations of state for xenon-deuterium using density functional theory

We report on a theoretical study of equation of state (EOS) properties of fluid and dense plasma mixtures of xenon and deuterium to explore and illustrate the basic physics of the mixing of a light element with a heavy element. Accurate EOS models are crucial to achieve high-fidelity hydrodynamics simulations of many high-energy-density phenomena, for example inertial confinement fusion and strong shock waves. While the EOS is often tabulated for separate species, the equation of state for arbitrary mixtures is generally not available, requiring properties of the mixture to be approximated by combining physical properties of the pure systems. Density functional theory (DFT) at elevated-temperature is used to assess the thermodynamics of the xenon-deuterium mixture at different mass ratios. The DFT simulations are unbiased as to elemental species and therefore provide comparable accuracy when describing total energies, pressures, and other physical properties of mixtures as they do for pure systems. The st...

Scaling relationships for nonadiabatic energy relaxation times in warm dense matter: toward understanding the equation of state

Understanding the dynamics of electron-ion energy transfer in warm dense (WD) matter is important to the measurement of equation of state (EOS) properties and for understanding the energy balance in dynamic simulations. In this work, we present a comprehensive investigation of nonadiabatic electron relaxation and thermal excitation dynamics in aluminum under high pressure and temperature. Using quantum-classical trajectory surface hopping approaches, we examine the role of nonadiabatic couplings and electronic decoherence in electron-nuclear energy transfer in WD aluminum. The computed timescales range from 400 fs to 4.0 ps and are consistent with existing experimental studies. We have derived general scaling relationships between macroscopic parameters of WD systems such as temperature or mass density and the timescales of energy redistribution between quantum and classical degrees of freedom. The scaling laws are supported by computational results. We show that electronic decoherence plays essential role and can change the functional dependencies qualitatively. The established scaling relationships can be of use in modelling of WD matter.

DENSITY FUNCTIONAL THEORY (DFT) SIMULATIONS OF SHOCKED LIQUID XENON

Xenon is not only a technologically important element used in laser technologies and jet propulsion, but it is also one of the most accessible materials in which to study the metal‐insulator transition with increasing pressure. Because of its closed shell electronic configuration, xenon is often assumed to be chemically inert, interacting almost entirely through the van der Waals interaction, and at liquid density, is typically modeled well using Leonard‐Jones potentials. However, such modeling has a limited range of validity as xenon is known to form compounds under normal conditions and likely exhibits considerably more chemistry at higher densities when hybridization of occupied orbitals becomes significant. We present DFT‐MD simulations of shocked liquid xenon with the goal of developing an improved equation of state. The calculated Hugoniot to 2 MPa compares well with available experimental shock data. Sandia is a mul‐tiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for...

LDRD project 151362 : low energy electron-photon transport.

At sufficiently high energies, the wavelengths of electrons and photons are short enough to only interact with one atom at time, leading to the popular %E2%80%9Cindependent-atom approximation%E2%80%9D. We attempted to incorporate atomic structure in the generation of cross sections (which embody the modeled physics) to improve transport at lower energies. We document our successes and failures. This was a three-year LDRD project. The core team consisted of a radiation-transport expert, a solid-state physicist, and two DFT experts.