Publications

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Significant progress has been made over the last few years in understanding properties of matter subject to strong shocks and other extreme conditions. High-accuracy multi-Mbar experiments and first-principles theoretical studies together provide detailed insights into the physics and chemistry of high energy-density matter. While comprehensive advances have been made for pure elements like deuterium, helium, and carbon, progress has been slower for equally important, albeit more challenging, materials like molecular crystals, polymers, and foams. Hydrocarbon based polymer foams are common materials and in particular they are used in designing shock- and inertial confinement fusion experiments. Depending on their initial density, foams shock to relatively higher pressure and temperature compared to shocked dense polymers/plastics. As foams and polymers are shocked, they exhibit both structural and chemical transitions. We will present experimental and theoretical results for shocked polymers in the Mbar regime. By shock impact of magnetically launched flyer plates on poly(4-methyl-1-pentene) foams, we create multi-Mbar pressures in a dense plasma mixture of hydrogen, carbon, at temperatures of several eV. Concurrently with executing experiments, we analyze the system by multi-scale simulations, from density functional theory to continuum magneto-hydrodynamics simulations. In particular, density functional theory (DFT) molecular dynamics (MD) and classical MD simulations of the principal shock Hugoniot will be presented in detail for two hydrocarbon polymers: polyethylene (PE) and poly(4-methyl-1-pentene) (PMP).

The difficulty of calculating the ambient properties of molecular crystals, such as the explosive PETN, has long hampered much needed computational investigations of these materials. One reason for the shortcomings is that the exchange-correlation functionals available for Density Functional Theory (DFT) based calculations do not correctly describe the weak intermolecular van der Waals' forces present in molecular crystals. However, this weak interaction also poses other challenges for the computational schemes used. We will discuss these issues in the context of calculations of lattice constants and structure of PETN with a number of different functionals, and also discuss if these limitations can be circumvented for studies at non-ambient conditions.

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 the United States Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000. © 2009 American Institute of Physics.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We present extensive simulations modeling the casting of multiblock polymer films by evaporation. The domain structure of the resulting film is strongly affected by varying the relative stiffness of the coblocks. The morphology changes from a bicontinuous lamellar phase when both blocks are flexible to a small-scale phase-separated phase with isolated domains as the stiffness of one of the blocks increases. As the relative stiffness of the blocks changes, the rate of evaporation, interfacial width, and morphology of the system changes. The findings can be used to tailor membrane morphology of interest to fuel-cell applications where the morphology is important for proton conduction.

Warm dense matter is the region in phase space of density and temperature where the thermal, Fermi, and Coulomb energies are approximately equal. The lack of a dominating scale and physical behavior makes it challenging to model the physics to high fidelity. For Sandia, a fundamental understanding of the region is of importance because of the needs of our experimental HEDP programs for high fidelity descriptive and predictive modeling. We show that multi-scale simulations of macroscopic physical phenomena now have predictive capability also for difficult but ubiquitous materials such as stainless steel, a transition metal alloy.

Abstract not provided.

Mechanisms for enhanced low-dose-rate sensitivity are described. In these mechanisms, bimolecular reactions dominate the kinetics at high dose rates thereby causing a sub-linear dependence on total dose, and this leads to a dose-rate dependence. These bimolecular mechanisms include electron-hole recombination, hydrogen recapture at hydrogen source sites, and hydrogen dimerization to form hydrogen molecules. The essence of each of these mechanisms is the dominance of the bimolecular reactions over the radiolysis reaction at high dose rates. However, at low dose rates, the radiolysis reaction dominates leading to a maximum effect of the radiation.

Density functional calculations show that the electric field effect on Si ad-dimer diffusion on Si(0 0 1) is largely a reflection of the position dependence of the ad-dimer's dipole moment. Surface diffusion barriers' dependence on perpendicular electric fields can be used to discriminate between diffusion mechanisms. Since the previously accepted mechanism for ad-dimer diffusion on Si(0 0 1) has the opposite field dependence to what is observed, it cannot be the one that dominates mass-transport. We identify an alternate process, with a similar barrier at zero electric field and field dependence in agreement with measurements. For rotation, calculations to date show linear field dependence, in contrast to experiments. © 2003 Elsevier Science B.V. All rights reserved.

Abstract not provided.

Abstract not provided.

We change the diffusion mechanism of adsorbed Ge-Si dimers on Si(001) using the electric field of a scanning tunneling microscope tip. By comparing the measured field dependence with first-principles calculations we conclude that, in negative field, i.e., when electrons are attracted towards the vacuum, the dimer diffuses as a unit, rotating as it translates, whereas, in positive field the dimer bond is substantially stretched at the transition state as it slides along the substrate. Furthermore, the active mechanism in positive fields facilitates intermixing of Ge in the Si lattice, whereas intermixing is suppressed in negative fields. © 2003 The American Physical Society.

Abstract not provided.