Mesoscale Modeling of Explosives at Sandia National Laboratories: Past and Future Directions a GA Tech Graduate Seminar Series Lecture
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Physical Review B
In both continuum hydrodynamics simulations and also multimillion atom reactive molecular dynamics simulations of shockwave propagation in single crystal pentaerythritol tetranitrate (PETN) containing a cylindrical void, we observed the formation of an initial radially symmetric hot spot. By extending the simulation time to the nanosecond scale, however, we observed the transformation of the small symmetric hot spot into a longitudinally asymmetric hot region extending over a much larger volume. Performing reactive molecular dynamics shock simulations using the reactive force field (ReaxFF) as implemented in the LAMMPS molecular dynamics package, we showed that the longitudinally asymmetric hot region was formed by coalescence of the primary radially symmetric hot spot with a secondary triangular hot zone. We showed that the triangular hot zone coincided with a double-shocked region where the primary planar shockwave was overtaken by a secondary cylindrical shockwave. The secondary cylindrical shockwave originated in void collapse after the primary planar shockwave had passed over the void. A similar phenomenon was observed in continuum hydrodynamics shock simulations using the CTH hydrodynamics package. The formation and growth of extended asymmetric hot regions on nanosecond timescales has important implications for shock initiation thresholds in energetic materials.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
15th International Detonation Symposium
Abstract not provided.
Abstract not provided.
Density Functional Theory (DFT) has emerged as an indispensable tool in materials research, since it can accurately predict properties of a wide variety of materials at both equilibrium and extreme conditions. However, for organic molecular crystal explosives, successful application of DFT has largely failed due to the inability of current exchange-correlation functionals to correctly describe intermolecular van der Waals (vdWs) forces. Despite this, we have discovered that even with no treatment of vdWs bonding, the AM05 functional and DFT based molecular dynamics (MD) could be used to study the properties of molecular crystals under compression. We have used DFT-MD to predict the unreacted Hugoniots for PETN and HNS and validated the results by comparison with crystalline and porous experimental data. Since we are also interested in applying DFT methods to study the equilibrium volume properties of explosives, we studied the nature of the vdWs bonding in pursuit of creating a new DFT functional capable of accurately describing equilibrium bonding of molecular crystals. In this report we discuss our results for computing shock Hugoniots of molecular crystals and also what was learned about the nature of bonding in these materials.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Density Functional Theory (DFT) has over the last few years emerged as an indispensable tool for understanding the behavior of matter under extreme conditions. DFT based molecular dynamics simulations (MD) have for example confirmed experimental findings for shocked deuterium, enabled the first experimental evidence for a triple point in carbon above 850 GPa, and amended experimental data for constructing a global equation of state (EOS) for water, carrying implications for planetary physics. The ability to perform high-fidelity calculations is even more important for cases where experiments are impossible to perform, dangerous, and/or prohibitively expensive. For solid explosives, and other molecular crystals, similar success has been severely hampered by an inability of describing the materials at equilibrium. The binding mechanism of molecular crystals (van der Waals forces) is not well described within traditional DFT. Among widely used exchange-correlation functionals, neither LDA nor PBE balances the strong intra-molecular chemical bonding and the weak inter-molecular attraction, resulting in incorrect equilibrium density, negatively affecting the construction of EOS for undetonated high explosives. We are exploring a way of bypassing this problem by using the new Armiento-Mattsson 2005 (AM05) exchange-correlation functional. The AM05 functional is highly accurate for a wide range of solids, in particular in compression. In addition, AM05 does not include any van der Waals attraction, which can be advantageous compared to other functionals: Correcting for a fictitious van der Waals like attraction with unknown origin can be harder than correcting for a complete absence of all types of van der Waals attraction. We will show examples from other materials systems where van der Waals attraction plays a key role, where this scheme has worked well, and discuss preliminary results for molecular crystals and explosives.
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.
Abstract not provided.
Sandia journal manuscript; Not yet accepted for publication
In this work, we examine the formation energies of interstitials in semiconductors obtained with four different pure functionals. Explicitely we investigate three silicon self-interstitials. All functionals give the same trend among those interstitials; the lowest energy being for formation of the <110>-split, somewhat higher energy for the formation of the hexagonal interstitial, while highest energy among the three is obtained for the meta-stable tetragonal configuration. However, the value for the formation energy for a specific interstitial differs substantially in calculations using different functionals. It is shown that the main contribution to these differences is stemming from the functionals different surface intrinsic errors. We also discuss the puzzle that the values obtained with the surface intrisic error free AM05 functional (Armiento and Mattsson, Phys. Rev. B 72, 085108 (2006)) gives values substantially lower than Quantum Monte Carlo results
Abstract not provided.