Publications

This is the seminar I will present at WCCM conference highlighting our latest research work on incorporating genetic programming to obtain data-driven strength models for complex materials.

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.

This report is the final documentation for the one-year LDRD project 226360: Simulated X-ray Diffraction and Machine Learning for Optimizing Dynamic Experiment Analysis. As Sandia has successfully developed in-house X-ray diffraction tools for study of atomic structure in experiments, it has become increasingly important to develop computational analysis methods to support these experiments. When dynamically compressed lattices and orientations are not known a priori, the identification requires a cumbersome and sometimes intractable search of possible final states. These final states can include phase transition, deformation and mixed/evolving states. Our work consists of three parts: (1) development of an XRD simulation tool and use of traditional data science methods to match XRD patterns to experiments; (2) development of ML-based models capable of decomposing and identifying the lattice and orientation components of multicomponent experimental diffraction patterns; and (3) conducting experiments which showcase these new analysis tools in the study of phase transition mechanisms. Our target material has been cadmium sulfide, which exhibits complex orientation-dependent phase transformation mechanisms. In our current one-year LDRD, we have begun the analysis of high-quality c-axis CdS diffraction data from DCS and Thor experiments, which had until recently eluded orientation identification.

Abstract not provided.

Multi-phase problems have so many more unknowns, we’d like to have a tool to constrain some open questions related to microstructure and twin & dislocation behavior. We want an atomistic scale perspective on aspects of strength. Some multi-scale questions accessible to atomistic study: What lattice-specific behavior influences dislocation production/mobility and/or twinning? Do the phase transformations wipe-out, modify or preserve grain size and orientation? Does plastic strain reset at phase transition? If so under what conditions? Tin is the material chosen for the effort because it is non-hazardous and has multiple accessible solid phases at relatively low pressures.

Abstract not provided.

Abstract not provided.

Pulsed-power generators can produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies using the magnetic loading technique. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength–strain rate dependence, structural phase transitions, and density of crystal defects, such as dislocations. Here, we present a cost-effective, compact, pulsed x-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically driven ramp compression of materials with a single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials’ community to investigate in situ dynamic phase transitions critical to equation of states. Finally, we present results using this new diagnostic to evaluate lattice compression in Zr and Al and to capture signatures of phase transitions in CdS.

Abstract: Here, we report the high pressure phase and morphology behavior of ordered anatase titanium dioxide (TiO2) nanocrystal arrays. One-dimensional TiO2 nanorods and nanorices were synthesized and self-assembled into ordered mesostructures. Their phase and morphological transitions at both atomic scale and mesoscale under pressure were studied using in situ synchrotron wide- and small-angle x-ray scattering (WAXS and SAXS) techniques. At the atomic scale, synchrotron WAXS reveals a pressure-induced irreversible amorphization up to 35 GPa in both samples but with different onset pressures. On the mesoscale, no clear phase transformations were observed up to 20 GPa by synchrotron SAXS. Intriguingly, sintering of TiO2 nanorods at mesoscale into nano-squares or nano-rectangles, as well as nanorices into nanowires, were observed for the first time by transmission electron microscopy. Such pressure-induced nanoparticle phase-amorphization and morphological changes provide valuable insights for design and engineering structurally stable nanomaterials. Impact statement: The high pressure behavior of nanocrystals (NCs) continues to be of interest, as previous studies have demonstrated that an externally applied pressure can serve as an efficient tool to induce structural phase transitions of NC assemblies at both the atomic scale and mesoscale without altering any chemistry by manipulating NC interatomic and interparticle distances. In addition, the high pressure generated deviatoric stress has been proven to be able to force adjacent NCs to connect and fuse into new crystalline nanostructures. Although the atomic structural evolution of TiO2 NCs under pressure has been widely investigated in the past decades, open questions remain regarding the mesoscale phase transition and morphology of TiO2 NC assemblies as a function of pressure. Therefore, in this work, systemic high pressure experiments on ordered arrays of TiO2 nanorods and nanorices were conducted by employing wide/small angle x-ray scattering techniques. The sintering of TiO2 assemblies at mesoscale into various nanostructures under pressure were revealed by transmission electron microscopy. Overall, this high pressure work fills the current gap in research on the mesoscale phase behavior of TiO2 assemblies. The observed morphology tunability attained by applying pressure opens new pathways for engineering nanomaterials and optimizing their collective properties through mechanical compression stresses. Graphical abstract: [Figure not available: see fulltext.].

This report documents details of the microstructure and mechanical properties of -tin (Sn), that is used in the Tri-lab (Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), Sandia National Laboratories (SNL)) collaboration project on Multi-phase Tin Strength. We report microstructural features detailing the crystallographic texture and grain morphology of as-received -tin from electron back scatter diffraction (EBSD). Temperature and strain rate dependent mechanical behavior was investigated by multiple compression tests at temperatures of 200K to 400K and strain rates of 0.0001 /s to 100 /s. Tri-lab tin showed significant temperature and strain rate dependent strength with no significant plastic anisotropy. A sample to sample material variation was observed from duplicate compression tests and texture measurements. Compression data was used to calibrate model parameters for temperature and rate dependent strength models, Johnson-Cook (JC), Zerilli-Armstrong (ZA) and Preston-Tonks-Wallace (PTW) strength models.

Abstract not provided.

Abstract not provided.

Pulsed-power generators using the magnetic loading technique are able to produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength-strain rate dependence, structural phase transitions, and density of crystal defects such as dislocations. Here, we present a cost effective, compact X-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically-driven ramp compression of materials with single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials community. The success in fielding this new XRD diagnostic dramatically improves our predictive capability and understanding of rate-dependent behavior at or near phase transition. As Sandia plans the next-generation pulse-power driver platform, a key element needed to deliver new state-of-the-art experiments will be having the necessary diagnostic tools to probe new regimes and phenomena. These diagnostics need to be as versatile, compact, and portable as they are powerful. The development of a platform-independent XRD diagnostic gives Sandia researchers a new window to study the microstructure and phase dynamics of materials under load. This project has paved the way for phase transition research in a variety of materials with mission interest.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.

Here, we describe recent efforts to improve our predictive modeling of rate-dependent behavior at, or near, a phase transition using molecular dynamics simulations. Cadmium sulfide (CdS) is a well-studied material that undergoes a solid-solid phase transition from wurtzite to rock salt structures between 3 and 9 GPa. Atomistic simulations are used to investigate the dominant transition mechanisms as a function of orientation, size and rate. We found that the final rock salt orientations were determined relative to the initial wurtzite orientation, and that these orientations were different for the two orientations and two pressure regimes studied. The CdS solid-solid phase transition is studied, for both a bulk single crystal and for polymer-encapsulated spherical nanoparticles of various sizes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Pressure-driven assembly of ligand-grafted gold nanoparticle superlattices is a promising approach for fabricating gold nanostructures, such as nanowires and nanosheets. Optimizing this fabrication method will require extending our understanding of superlattice mechanics to regimes of high pressures. We use molecular dynamics simulations to characterize the response of alkanethiol-grafted gold nanoparticle superlattices to applied hydrostatic pressures up to 15 GPa. At low applied pressures, intrinsic voids govern the mechanics of compaction. As applied pressures increase, the void collapse and ligand compression depend significantly on the ligand length. These microstructural observations correlate directly with trends in bulk modulus and elastic constants. For short ligands, core-core contact between gold nanoparticles is observed at high pressures, which augurs irreversible response and eventual sintering. This presintering behavior was unexpected under hydrostatic loading and is observed only for the shortest ligands.

Abstract not provided.

Abstract not provided.

Deformation mechanisms in bcc metals, especially in dynamic regimes, show unusual complexity, which complicates their use in high-reliability applications. Here, we employ novel, high-velocity cylinder impact experiments to explore plastic anisotropy in single crystal specimens under high-rate loading. The bcc tantalum single crystals exhibit unusually high deformation localization and strong plastic anisotropy when compared to polycrystalline samples. Several impact orientations - [100], [110], [111] and [149] -Are characterized over a range of impact velocities to examine orientation-dependent mechanical behavior versus strain rate. Moreover, the anisotropy and localized plastic strain seen in the recovered cylinders exhibit strong axial symmetries which differed according to lattice orientation. Two-, three-, and four-fold symmetries are observed. We propose a simple crystallographic argument, based on the Schmid law, to understand the observed symmetries. These tests are the first to explore the role of single-crystal orientation in Taylor impact tests and they clearly demonstrate the importance of crystallography in high strain rate and temperature deformation regimes. These results provide critical data to allow dramatically improved high-rate crystal plasticity models and will spur renewed interest in the role of crystallography to deformation in dynamics regimes.

Pressure-driven assembly of ligand-grafted gold nanoparticle superlattices is a promising approach for fabricating gold nanostructures, such as nanowires and nanosheets. However, optimizing this fabrication method requires an understanding of the mechanics of their complex hierarchical assemblies at high pressures. We use molecular dynamics simulations to characterize the response of alkanethiol-grafted gold nanoparticle superlattices to applied hydrostatic pressures up to 15 GPa, and demonstrate that the internal mechanics significantly depend on ligand length. At low pressures, intrinsic voids govern the mechanics of pressure-induced compaction, and the dynamics of collapse of these voids under pressure depend significantly on ligand length. These microstructural observations correlate well with the observed trends in bulk modulus and elastic constants. For the shortest ligands at high pressures, coating failure leads to gold core-core contact, an augur of irreversible response and eventual sintering. This behavior was unexpected under hydrostatic loading, and was only observed for the shortest ligands.

Abstract not provided.

The shock response of porous amorphous silica was investigated using classical molecular dynamics, over a range of porosity ranging from fully dense (2.21 g/cc) down to 0.14 g/cc. We observed an enhanced densification in the Hugoniot response at initial porosities above 50%, and the effect increased with increasing porosity. In the lowest initial densities, after an initial compression response, the systems expanded with increased pressure. These results show good agreement with experiments. We explored mechanisms leading to enhanced densification which appear to differ from mechanisms observed in similar studies in silicon.

The shock response of porous amorphous silica was investigated using classical molecular dynamics, over a range of porosity ranging from fully dense (2.21 g/cc) down to 0.14 g/cc. We observed an enhanced densification in the Hugoniot response at initial porosities above 50%, and the effect increased with increasing porosity. In the lowest initial densities, after an initial compression response, the systems expanded with increased pressure. These results show good agreement with experiments. We explored mechanisms leading to enhanced densification which appear to differ from mechanisms observed in similar studies in silicon.

Material strength and moduli can be determined from dynamic high-pressure ramp-release experiments using an indirect method of Lagrangian wave profile analysis of surface velocities. This method, termed self-consistent Lagrangian analysis (SCLA), has been difficult to calibrate and corroborate with other experimental methods. Using nonequilibrium molecular dynamics, we validate the SCLA technique by demonstrating that it accurately predicts the same bulk modulus, shear modulus, and strength as those calculated from the full stress tensor data, especially where strain rate induced relaxation effects and wave attenuation are small. We show here that introducing a hold in the loading profile at peak pressure gives improved accuracy in the shear moduli and relaxation-adjusted strength by reducing the effect of wave attenuation. When rate-dependent effects coupled with wave attenuation are large, we find that Lagrangian analysis overpredicts the maximum unload wavespeed, leading to increased error in the measured dynamic shear modulus. These simulations provide insight into the definition of dynamic strength, as well as a plausible explanation for experimental disagreement in reported dynamic strength values.

Abstract not provided.

Thermal degradation of polyethylene is studied under the extremely high rate temperature ramps expected in laser-driven and X-ray ablation experiments - from 1010 to 1014 K/s in isochoric, condensed phases. The molecular evolution and macroscopic state variables are extracted as a function of density from reactive molecular dynamics simulations using the ReaxFF potential. The enthalpy, dissociation onset temperature, bond evolution, and observed cross-linking are shown to be rate dependent. These results are used to parametrize a kinetic rate model for the decomposition and coalescence of hydrocarbons as a function of temperature, temperature ramp rate, and density. The results are contrasted to first-order random-scission macrokinetic models often assumed for pyrolysis of linear polyethylene under ambient conditions.

The shock Hugoniot for full-density and porous CeO2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ=2.5 to 20g/cm3. The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO2. In addition, a revised multiphase SESAME-Type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

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.

Gold nanostructured materials exhibit important size- and shape-dependent properties that enable a wide variety of applications in photocatalysis, nanoelectronics and phototherapy. Here we show the use of superfast dynamic compression to synthesize extended gold nanostructures, such as nanorods, nanowires and nanosheets, with nanosecond coalescence times. Using a pulsed power generator, we ramp compress spherical gold nanoparticle arrays to pressures of tens of GPa, demonstrating pressure-driven assembly beyond the quasi-static regime of the diamond anvil cell. Our dynamic magnetic ramp compression approach produces smooth, shockless (that is, isentropic) one-dimensional loading with low-temperature states suitable for nanostructure synthesis. Transmission electron microscopy clearly establishes that various gold architectures are formed through compressive mesoscale coalescences of spherical gold nanoparticles, which is further confirmed by in-situ synchrotron X-ray studies and large-scale simulation. This nanofabrication approach applies magnetically driven uniaxial ramp compression to mimic established embossing and imprinting processes, but at ultra-short (nanosecond) timescales.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The aim of this study was to alter polymerization chemistry to improve network homogeneity in free-radical crosslinked systems. It was hypothesized that a reduction in heterogeneity of the network would lead to improved mechanical performance. Experiments and simulations were carried out to investigate the connection between polymerization chemistry, network structure and mechanical properties. Experiments were conducted on two different monomer systems - the first is a single monomer system, urethane dimethacrylate (UDMA), and the second is a two-monomer system consisting of bisphenol A glycidyl dimethacrylate (BisGMA) and triethylene glycol dimethacrylate (TEGDMA) in a ratio of 70/30 BisGMA/TEGDMA by weight. The methacrylate systems were crosslinked using traditional radical polymeriza- tion (TRP) with azobisisobutyronitrile (AIBN) or benzoyl peroxide (BPO) as an initiator; TRP systems were used as the control. The monomers were also cross-linked using activator regenerated by electron transfer atom transfer radical polymerization (ARGET ATRP) as a type of controlled radical polymerization (CRP). FTIR and DSC were used to monitor reac- tion kinetics of the systems. The networks were analyzed using NMR, DSC, X-ray diffraction (XRD), atomic force microscopy (AFM), and small angle X-ray scattering (SAXS). These techniques were employed in an attempt to quantify differences between the traditional and controlled radical polymerizations. While a quantitative methodology for characterizing net- work morphology was not established, SAXS and AFM have shown some promising initial results. Additionally, differences in mechanical behavior were observed between traditional and controlled radical polymerized thermosets in the BisGMA/TEGDMA system but not in the UDMA materials; this finding may be the result of network ductility variations between the two materials. Coarse-grained molecular dynamics simulations employing a novel model of the CRP reaction were carried out for the UDMA system, with parameters calibrated based on fully atomistic simulations of the UDMA monomer in the liquid state. Detailed metrics based on network graph theoretical approaches were implemented to quantify the bond network topology resulting from simulations. For a broad range of polymerization parameters, no discernible differences were seen between TRP and CRP UDMA simulations at equal conversions, although clear differences exist as a function of conversion. Both findings are consistent with experiments. Despite a number of shortcomings, these models have demonstrated the potential of molecular simulations for studying network topology in these systems.

We have conducted molecular dynamics (MD) simulations of quasi-isentropic ramp-wave compression to very high pressures over a range of strain rates from 1011 down to 108 1/s. Using scaling methods, we collapse wave profiles from various strain rates to a master profile curve, which shows deviations when material response is strain-rate dependent. Thus, we can show with precision where, and how, strain-rate dependence affects the ramp wave. We find that strain rate affects the stress-strain material response most dramatically at strains below 20%, and that above 30% strain the material response is largely independent of strain rate. We show good overall agreement with experimental stress-strain curves up to approximately 30% strain, above which simulated response is somewhat too stiff. We postulate that this could be due to our interatomic potential or to differences in grain structure and/or size between simulation and experiment. Strength is directly measured from per-atom stress tensor and shows significantly enhanced elastic response at the highest strain rates. This enhanced elastic response is less pronounced at higher pressures and at lower strain rates.

Abstract not provided.

Nanoparticles (NPs) grafted with organic layers form hybrids able to retain their unique properties through integration into the mesoscopic scale. The organic layer structure and response often determine the functionality of the hybrids on the mesoscopic length scale. Using molecular dynamics (MD) simulations, we probe the conformation of luminescent rigid polymers, dialkyl poly(p-phenylene ethynylene)s (PPE), end-grafted onto a silica nanoparticle in different solvents as the molecular weights and polymer coverages are varied. We find that, in contrast to NP-grafted flexible polymers, the chains are fully extended independent of the solvent. In toluene and decane, which are good solvents, the grafted PPEs chains assume a similar conformation to that observed in dilute solutions. In water, which is a poor solvent for the PPEs, the polymer chains form one large cluster but remain extended. The radial distribution of the chains around the core of the nanoparticle is homogeneous in good solvents, whereas in poor solvents clusters are formed independent of molecular weights and coverages. The clustering is distinctively different from the response of grafted flexible and semiflexible polymers.

Abstract not provided.

Abstract not provided.

We studied dynamic ramp waves using molecular dynamics to understand the effects of strain rate and nanograin structure on material response and strength.

We have conducted extremely long molecular dynamics simulations of glasses to microsecond times, which close the gap between experimental and atomistic simulation time scales by two to three orders of magnitude. Static, thermal, and structural properties of silica glass are reported for glass cooling rates down to 5×109 K/s and viscoelastic response in silica melts and glasses are studied over nine decades of time. We present results from relaxation of hydrostatic compressive stress in silica and show that time-temperature superposition holds in these systems for temperatures from 3500 to 1000 K.

Abstract not provided.

The high mechanical stiffness of single-nanoparticle-thick membranes is believed to result from the local structure of ligand coatings that mediate interactions between nanoparticles. These ligand structures are not directly observable experimentally. We use molecular dynamics simulations to observe variations in ligand structure and simultaneously measure variations in membrane mechanical properties. We have shown previously that ligand end group has a large impact on ligand structure and membrane mechanical properties. Here we introduce and apply quantitative molecular structure measures to these membranes and extend analysis to multiple nanoparticle core sizes and ligand lengths. Simulations of nanoparticle membranes with a nanoparticle core diameter of 4 or 6 nm, a ligand length of 11 or 17 methylenes, and either carboxyl (COOH) or methyl (CH3) ligand end groups are presented. In carboxyl-terminated ligand systems, structure and interactions are dominated by an end-to-end orientation of ligands. In methyl-terminated ligand systems large ordered ligand structures form, but nanoparticle interactions are dominated by disordered, partially interdigitated ligands. Core size and ligand length also affect both ligand arrangement within the membrane and the membrane's macroscopic mechanical response, but are secondary to the role of the ligand end group. Moreover, the particular end group (COOH or CH3) alters the nature of how ligand length, in turn, affects the membrane properties. The effect of core size does not depend on the ligand end group, with larger cores always leading to stiffer membranes. Asymmetry in the stress and ligand density is observed in membranes during preparation at a water-vapor interface, with the stress asymmetry persisting in all membranes after drying.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Under shock compression, most porous materials exhibit lower densities for a given pressure than that of a full-dense sample of the same material. However, some porous materials exhibit an anomalous, or enhanced, densification under shock compression. We demonstrate a molecular mechanism that drives this behavior. We also present evidence from atomistic simulation that silicon belongs to this anomalous class of materials. Atomistic simulations indicate that local shear strain in the neighborhood of collapsing pores nucleates a local solid-solid phase transformation even when bulk pressures are below the thermodynamic phase transformation pressure. This metastable, local, and partial, solid-solid phase transformation, which accounts for the enhanced densification in silicon, is driven by the local stress state near the void, not equilibrium thermodynamics. This mechanism may also explain the phenomenon in other covalently bonded materials.

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.

Hydrocarbon polymers, foams and nanocomposites are increasingly being subjected to extreme environments. Molecular scale modeling of these materials offers insight into failure mechanisms and complex response. Prior classical molecular dynamics (MD) simulations of the principal shock Hugoniot for two hydrocarbon polymers, polyethylene (PE) and poly (4-methyl-1-pentene) (PMP) have shown good agreement with density functional theory (DFT) calculations and experiments conducted at Sandia National Laboratories. We extended these results to include low-density polymer foams using nonequilibrium MD techniques and found good quantitative agreement with experiment. Here, we have measured the local temperature during void collapse to investigate the formation of hot spots and their relationship to polymer dissociation in foams. © 2013 Elsevier B.V. All rights reserved.

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.

The response of beryllium to dynamic loading has been extensively studied, both experimentally and theoretically, due to its importance in several technological areas. We use a MEAM empirical potential to examine the melt transition. MD simulations of equilibrated two-phase systems were used to calculate the HCP melting curve up to 300 GPa. This was found to agree well with previous ab initio calculations. The Hugoniostat method was used to examine dynamic compression along the two principal orientations of the HCP crystal. In both directions, the melting transition occurred at 230 GPa and 5000 K, consistent with the equilibrium melting curve. Direct NEMD simulations of uniaxial compression show a transition to an amorphous material at shocked states that lie below the equilibrium melt curve. © 2012 American Institute of Physics.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Hydrocarbon polymers, foams and nanocomposites are increasingly being subjected to extreme environments. Molecular scale modeling of these materials offers insight into failure mechanisms and complex response. Classical molecular dynamics (MD) simulations of the principal shock Hugoniot were conducted for two hydrocarbon polymers, polyethylene (PE) and poly(4-methyl-1-pentene) (PMP). We compare these results with recent density functional theory (DFT) calculations and experiments conducted at Sandia National Laboratories. Here, we extend these results to include low-density polymer foams using nonequilibrium MD techniques. We find good quantitative agreement with experiment. Further, we have measured local temperatures to investigate the formation of hot spots and polymer dissociation near foam voids.

Abstract not provided.

Molecular dynamics simulation (MD) is an invaluable tool for studying problems sensitive to atomscale physics such as structural transitions, discontinuous interfaces, non-equilibrium dynamics, and elastic-plastic deformation. In order to apply this method to modeling of ramp-compression experiments, several challenges must be overcome: accuracy of interatomic potentials, length- and time-scales, and extraction of continuum quantities. We have completed a 3 year LDRD project with the goal of developing molecular dynamics simulation capabilities for modeling the response of materials to ramp compression. The techniques we have developed fall in to three categories (i) molecular dynamics methods (ii) interatomic potentials (iii) calculation of continuum variables. Highlights include the development of an accurate interatomic potential describing shock-melting of Beryllium, a scaling technique for modeling slow ramp compression experiments using fast ramp MD simulations, and a technique for extracting plastic strain from MD simulations. All of these methods have been implemented in Sandia's LAMMPS MD code, ensuring their widespread availability to dynamic materials research at Sandia and elsewhere.

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 have conducted a molecular dynamics (MD) simulation study of water confined between methyl-terminated and carboxyl-terminated alkylsilane self-assembled monolayers (SAMs) on amorphous silica substrates. In doing so, we have investigated the dynamic and structural behavior of the water molecules when compressed to loads ranging from 20 to 950 MPa for two different amounts of water (27 and 58 water molecules/nm{sup 2}). Within the studied range of loads, we observe that no water molecules penetrate the hydrophobic region of the carboxyl-terminated SAMs. However, we observe that at loads larger than 150 MPa water molecules penetrate the methyl-terminated SAMs and form hydrogen-bonded chains that connect to the bulk water. The diffusion coefficient of the water molecules decreases as the water film becomes thinner and pressure increases. When compared to bulk diffusion coefficients of water molecules at the various loads, we found that the diffusion coefficients for the systems with 27 water molecules/nm{sup 2} are reduced by a factor of 20 at low loads and by a factor of 40 at high loads, while the diffusion coefficients for the systems with 58 water molecules/nm{sup 2} are reduced by a factor of 25 at all loads.

Abstract not provided.

Abstract not provided.