Publications

Results 1–50 of 54
Skip to search filters

Stopping of Deuterium in Warm Dense Deuterium from Ehrenfest Time-Dependent Density Functional Theory

Contributions to Plasma Physics

Magyar, Rudolph J.; Shulenburger, Luke N.; Baczewski, Andrew D.

In these proceedings, we show that time-dependent density functional theory is capable of stopping calculations at the extreme conditions of temperature and pressure seen in warm dense matter. The accuracy of the stopping curves tends to be up to about 20% lower than empirical models that are in use. However, TDDFT calculations are free from fitting parameters and assumptions about the model form of the dielectric function. This work allows the simulation of ion stopping in many materials that are difficult to study experimentally. (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim).

More Details
TYPE Conference Poster YEAR 2016
ScopusOSTIDOI

Ethane-xenon mixtures under shock conditions

Physical Review B - Condensed Matter and Materials Physics

Magyar, Rudolph J.; Root, Seth R.; Cochrane, Kyle C.; Mattsson, Thomas M.; Flicker, Dawn G.

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. 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. 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.

More Details
TYPE Journal Article YEAR 2015
ScopusOSTIDOI

Validating density-functional theory simulations at high energy-density conditions with liquid krypton shock experiments to 850 GPa on Sandia's Z machine

Physical Review B - Condensed Matter and Materials Physics

Mattsson, Thomas M.; Root, Seth R.; Mattsson, Ann E.; Shulenburger, Luke N.; Magyar, Rudolph J.; Flicker, Dawn G.

We use Sandia's Z machine and magnetically accelerated flyer plates to shock compress liquid krypton to 850 GPa and compare with results from density-functional theory (DFT) based simulations using the AM05 functional. We also employ quantum Monte Carlo calculations to motivate the choice of AM05. We conclude that the DFT results are sensitive to the quality of the pseudopotential in terms of scattering properties at high energy/temperature. A new Kr projector augmented wave potential was constructed with improved scattering properties which resulted in excellent agreement with the experimental results to 850 GPa and temperatures above 10 eV (110 kK). Finally, we present comparisons of our data from the Z experiments and DFT calculations to current equation of state models of krypton to determine the best model for high energy-density applications.

More Details
TYPE Journal Article YEAR 2014
ScopusOSTIDOI

Numerical implementation of time-dependent density functional theory for extended systems in extreme environments

Baczewski, Andrew D.; Shulenburger, Luke N.; Desjarlais, Michael P.; Magyar, Rudolph J.

In recent years, DFT-MD has been shown to be a useful computational tool for exploring the properties of WDM. These calculations achieve excellent agreement with shock compression experiments, which probe the thermodynamic parameters of the Hugoniot state. New X-ray Thomson Scattering diagnostics promise to deliver independent measurements of electronic density and temperature, as well as structural information in shocked systems. However, they require the development of new levels of theory for computing the associated observables within a DFT framework. The experimentally observable x-ray scattering cross section is related to the electronic density-density response function, which is obtainable using TDDFT - a formally exact extension of conventional DFT that describes electron dynamics and excited states. In order to develop a capability for modeling XRTS data and, more generally, to establish a predictive capability for rst principles simulations of matter in extreme conditions, real-time TDDFT with Ehrenfest dynamics has been implemented in an existing PAW code for DFT-MD calculations. The purpose of this report is to record implementation details and benchmarks as the project advances from software development to delivering novel scienti c results. Results range from tests that establish the accuracy, e ciency, and scalability of our implementation, to calculations that are veri ed against accepted results in the literature. Aside from the primary XRTS goal, we identify other more general areas where this new capability will be useful, including stopping power calculations and electron-ion equilibration.

More Details
TYPE SAND Report YEAR 2014
OSTIDOI

LDRD project 151362 : low energy electron-photon transport

Kensek, Ronald P.; Hjalmarson, Harold P.; Magyar, Rudolph J.; Bondi, Robert J.

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.

More Details
TYPE SAND Report YEAR 2013
OSTIDOI

The potential, limitations, and challenges of divide and conquer quantum electronic structure calculations on energetic materials

Tucker, Jon R.; Magyar, Rudolph J.

High explosives are an important class of energetic materials used in many weapons applications. Even with modern computers, the simulation of the dynamic chemical reactions and energy release is exceedingly challenging. While the scale of the detonation process may be macroscopic, the dynamic bond breaking responsible for the explosive release of energy is fundamentally quantum mechanical. Thus, any method that does not adequately describe bonding is destined to lack predictive capability on some level. Performing quantum mechanics calculations on systems with more than dozens of atoms is a gargantuan task, and severe approximation schemes must be employed in practical calculations. We have developed and tested a divide and conquer (DnC) scheme to obtain total energies, forces, and harmonic frequencies within semi-empirical quantum mechanics. The method is intended as an approximate but faster solution to the full problem and is possible due to the sparsity of the density matrix in many applications. The resulting total energy calculation scales linearly as the number of subsystems, and the method provides a path-forward to quantum mechanical simulations of millions of atoms.

More Details
TYPE SAND Report YEAR 2012
OSTIDOI

Physics of intense, high energy radiation effects

Hjalmarson, Harold P.; Magyar, Rudolph J.; Crozier, Paul C.; Hartman, Elmer F.

This document summarizes the work done in our three-year LDRD project titled 'Physics of Intense, High Energy Radiation Effects.' This LDRD is focused on electrical effects of ionizing radiation at high dose-rates. One major thrust throughout the project has been the radiation-induced conductivity (RIC) produced by the ionizing radiation. Another important consideration has been the electrical effect of dose-enhanced radiation. This transient effect can produce an electromagnetic pulse (EMP). The unifying theme of the project has been the dielectric function. This quantity contains much of the physics covered in this project. For example, the work on transient electrical effects in radiation-induced conductivity (RIC) has been a key focus for the work on the EMP effects. This physics in contained in the dielectric function, which can also be expressed as a conductivity. The transient defects created during a radiation event are also contained, in principle. The energy loss lead the hot electrons and holes is given by the stopping power of ionizing radiation. This information is given by the inverse dielectric function. Finally, the short time atomistic phenomena caused by ionizing radiation can also be considered to be contained within the dielectric function. During the LDRD, meetings about the work were held every week. These discussions involved theorists, experimentalists and engineers. These discussions branched out into the work done in other projects. For example, the work on EMP effects had influence on another project focused on such phenomena in gases. Furthermore, the physics of radiation detectors and radiation dosimeters was often discussed, and these discussions had impact on related projects. Some LDRD-related documents are now stored on a sharepoint site (https://sharepoint.sandia.gov/sites/LDRD-REMS/default.aspx). In the remainder of this document the work is described in catergories but there is much overlap between the atomistic calculations, the continuum calculations and the experiments.

More Details
TYPE SAND Report YEAR 2011
OSTIDOI

High fidelity equation of state for xenon : integrating experiments and first principles simulations in developing a wide-range equation of state model for a fifth-row element

Magyar, Rudolph J.; Root, Seth R.; Carpenter, John H.; Mattsson, Thomas M.

The noble gas xenon is a particularly interesting element. At standard pressure xenon is an fcc solid which melts at 161 K and then boils at 165 K, thus displaying a rather narrow liquid range on the phase diagram. On the other hand, under pressure the melting point is significantly higher: 3000 K at 30 GPa. Under shock compression, electronic excitations become important at 40 GPa. Finally, xenon forms stable molecules with fluorine (XeF{sub 2}) suggesting that the electronic structure is significantly more complex than expected for a noble gas. With these reasons in mind, we studied the xenon Hugoniot using DFT/QMD and validated the simulations with multi-Mbar shock compression experiments. The results show that existing equation of state models lack fidelity and so we developed a wide-range free-energy based equation of state using experimental data and results from first-principles simulations.

More Details
TYPE Conference YEAR 2010
OSTI
Results 1–50 of 54
Results 1–50 of 54