Search

We present results of resonant ultrasound spectroscopy (RUS) measurements applied to granular high explosive materials at different bulk pressing densities and degree of thermal conditioning. The material chosen in this study is a ubiquitously used explosive material known as pentaerythritol tetranitrate (PETN), which is used commercially in civil and defense applications both as a binderized plastic bonded explosive material and an unbinderized neat material. However, changes in granular PETN bulk elastic properties due to thermal conditioning, which could have implications for better understanding environmental aging-related effects, have not been well studied even though it is believed that elasticity may play an important role in explosive material initiation mechanisms. Furthermore, monitoring elastic property changes in granular explosive pressings has not yet been demonstrated using RUS, which is an appealing non-destructive characterization tool that requires only dry point contact with the explosive material. To this end, we report the first study using RUS to quantify the elastic properties of binderized and neat PETN pressings as well as to quantify changes in elastic properties as a function of both thermal conditioning and bulk pressing density. Elastic stiffness coefficients, sometimes more commonly referred to as elastic constants, calculated from the RUS measurements on the different PETN-based materials show a significant increase for the post-conditioned samples compared to the pre-conditioned samples. This trend of increasing elastic properties with thermal conditioning was consistent for different density pressings, different thermal exposure conditions, and even different neat PETN pressings of differing average crystal sizes and/or specific surface areas. [ABSTRACT FROM AUTHOR]

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. This suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects. [ABSTRACT FROM AUTHOR]

We describe a machine learning approach to rapidly tune density functional tight binding models for the description of detonation chemistry in organic molecular materials. Resulting models enable simulations on the several 10s of ps scales characteristic to these processes, with "quantum-accuracy." We use this approach to investigate early shock chemistry in 3,4-bis(3-nitrofurazan-4-yl)furoxan, a hydrogen-free energetic material known to form onion-like nanocarbon particulates following detonation. We find that the ensuing chemistry is significantly characterized by the formation of large CxNyOz species, which are likely precursors to the experimentally observed carbon condensates. Beyond utility as a means of investigating detonation chemistry, the present approach can be used to generate quantum-based reference data for the development of full machine-learned interatomic potentials capable of simulation on even greater time and length scales, i.e., for applications where characteristic time scales exceed the reach of methods including Kohn–Sham density functional theory, which are commonly used for reference data generation. [ABSTRACT FROM AUTHOR]

HN3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN3. We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory. [ABSTRACT FROM AUTHOR]

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. Here, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation. [ABSTRACT FROM AUTHOR]

We describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm−3, affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials. [ABSTRACT FROM AUTHOR]

Siloxane systems consisting primarily of polydimethylsiloxane (PDMS) are versatile, multifaceted materials that play a key role in diverse applications. However, open questions exist regarding the correlation between their varied atomic-level properties and observed macroscale features. To this effect, we have created a systematic workflow to determine coarse-grained simulation models for crosslinked PDMS in order to further elucidate the effects of network changes on the system's rheological properties below the gel point. Our approach leverages a fine-grained united atom model for linear PDMS, which we extend to include crosslinking terms, and applies iterative Boltzmann inversion to obtain a coarse-grain "bead-spring-type" model. We then perform extensive molecular dynamics simulations to explore the effect of crosslinking on the rheology of silicone fluids, where we compute systematic increases in both density and shear viscosity that compare favorably to experiments that we conduct here. The kinematic viscosity of partially crosslinked fluids follows an empirical linear relationship that is surprisingly consistent with Rouse theory, which was originally derived for systems comprised of a uniform distribution of linear chains. The models developed here serve to enable quantitative bottom-up predictions for curing- and age-induced effects on macroscale rheological properties, allowing for accurate prediction of material properties based on fundamental chemical data. [ABSTRACT FROM AUTHOR]