Cavitation and Hollow Foam-like Structure Formation in the Dense Core of Exploded Wire
Abstract not provided.
Abstract not provided.
Abstract not provided.
APS Shock Compression of Condensed Matter
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.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Location of the liquid-vapor critical point (c.p.) is one of the key features of equation of state models used in simulating high energy density physics and pulsed power experiments. For example, material behavior in the location of the vapor dome is critical in determining how and when coronal plasmas form in expanding wires. Transport properties, such as conductivity and opacity, can vary an order of magnitude depending on whether the state of the material is inside or outside of the vapor dome. Due to the difficulty in experimentally producing states near the vapor dome, for all but a few materials, such as Cesium and Mercury, the uncertainty in the location of the c.p. is of order 100%. These states of interest can be produced on Z through high-velocity shock and release experiments. For example, it is estimated that release adiabats from {approx}1000 GPa in aluminum would skirt the vapor dome allowing estimates of the c.p. to be made. This is within the reach of Z experiments (flyer plate velocity of {approx}30 km/s). Recent high-fidelity EOS models and hydrocode simulations suggest that the dynamic two-phase flow behavior observed in initial scoping experiments can be reproduced, providing a link between theory and experiment. Experimental identification of the c.p. in aluminum would represent the first measurement of its kind in a dynamic experiment. Furthermore, once the c.p. has been experimentally determined it should be possible to probe the electrical conductivity, opacity, reflectivity, etc. of the material near the vapor dome, using a variety of diagnostics. We propose a combined experimental and theoretical investigation with the initial emphasis on aluminum.
Abstract not provided.
Physical Review Letters
X-ray production by imploding wire-array Z pinches is studied using radiation magnetohydrodynamics simulation. It is found that the density distribution created by ablating wire material influences both x-ray power production, and how the peak power scales with applied current. For a given array there is an optimum ablation rate that maximizes the peak x-ray power, and produces the strongest scaling of peak power with peak current. This work is consistent with trends in wire-array Z pinch x-ray power scaling experiments on the Z accelerator. © 2009 The American Physical Society.
Physics of Plasmas
The implosion phase of a wire-array Z pinch is investigated using three-dimensional (3D) simulations, which model the mass ablation phase and its associated axial instability using a mass injection boundary condition. The physical mechanisms driving the trailing mass network are explored, and it is found that in 3D the current paths though the trailing mass can reduce bubble growth on the imploding plasma sheath, relative to the 2D (r,z) equivalent. Comparison between the simulations and a high quality set of experimental radiographs is presented. © 2008 American Institute of Physics.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Nature
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Proceedings APS Shock Compression of Condensed Matter 2007
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Proposed for publication in Physical Review Letters.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
A finite temperature version of 'exact-exchange' density functional theory (EXX) has been implemented in Sandia's Socorro code. The method uses the optimized effective potential (OEP) formalism and an efficient gradient-based iterative minimization of the energy. The derivation of the gradient is based on the density matrix, simplifying the extension to finite temperatures. A stand-alone all-electron exact-exchange capability has been developed for testing exact exchange and compatible correlation functionals on small systems. Calculations of eigenvalues for the helium atom, beryllium atom, and the hydrogen molecule are reported, showing excellent agreement with highly converged quantumMonte Carlo calculations. Several approaches to the generation of pseudopotentials for use in EXX calculations have been examined and are discussed. The difficult problem of finding a correlation functional compatible with EXX has been studied and some initial findings are reported.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Physical Review B - Condensed Matter and Materials Physics
We present an ab initio method for electronic structure calculations of materials at finite temperature (FT) based on the all-electron quasiparticle self-consistent GW (QPSCGW) approximation and Keldysh time-loop Green's function approach. We apply the method to Si, Ge, GaAs, InSb, and diamond and show that the band gaps of these materials universally decrease with increasing temperature in contrast with results of the local density approximation (LDA) of density functional theory where the band gaps universally increase. At temperatures of a few eV the difference between quasiparticle energies obtained in FT QPSCGW and FT LDA approaches is significantly reduced. This result suggests that existing simulations of very high-temperature materials based on the FT LDA are more justified than it might appear from the well-known LDA band gap errors at zero temperature. © 2006 The American Physical Society.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
We report the implementation of an iterative scheme for calculating the Optimized Effective Potential (OEP). Given an energy functional that depends explicitly on the Kohn-Sham wave functions, and therefore, implicitly on the local effective potential appearing in the Kohn-Sham equations, a gradient-based minimization is used to find the potential that minimizes the energy. Previous work has shown how to find the gradient of such an energy with respect to the effective potential in the zero-temperature limit. We discuss a density-matrix-based derivation of the gradient that generalizes the previous results to the finite temperature regime, and we describe important optimizations used in our implementation. We have applied our OEP approach to the Hartree-Fock energy expression to perform Exact Exchange (EXX) calculations. We report our EXX results for common semiconductors and ordered phases of hydrogen at zero and finite electronic temperatures. We also discuss issues involved in the implementation of forces within the OEP/EXX approach.
The optimized effective potential (OEP) method allows orbital-dependent functionals to be used in density functional theory (DFT), which, in particular, allows exact exchange formulations of the exchange energy to be used in DFT calculations. Because the exact exchange is inherently self-interaction correcting, the resulting OEP calculations have been found to yield superior band-gaps for condensed-phase systems. Here we apply these methods to the isolated atoms He and Be, and compare to high quality experiments and calculations to demonstrate that the orbital energies accurately reproduce the excited state spectrum for these species. These results suggest that coupling the exchange-only OEP calculations with proper (orbital-dependent or other) correlation functions might allow quantitative accuracy from DFT calculations.