Search

High pressure shear band formation is a critical phenomenon in energetic materials due to its influence on both mechanical strength and mechanochemical activation. While shear banding is know to occur in a variety of these materials, the governing dynamics of the mechanisms is not well defined for molecular crystals. We conduct molecular dynamics simulations of shock wave induced shear band formation in the energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) to assess shear band nucleation processes. We find, that at high pressures, the initial formation sites for shear bands 'over-nucleate' and rapidly lower deviatoric stresses prior to shear band formation and growth. This results in the suppression of plastic deformation. A local cluster analysis is used to quantify and contrast this mechanism with a more typical shear banding seen at lower pressures. These results demonstrate a mechanism that is reversible in nature and that supersedes shear band formation at increased pressures. We anticipate that these results will have a broad impact on the modeling and development of high strain rate application materials such as those for high explosives and hypersonic systems.

The rapid compaction of granular media results in localized heating that can induce chemical reactions, phase transformations, and melting. However, there are numerous mechanisms in play that can be dependent on a variety of microstructural features. Machine learning techniques such as neural networks offer a ubiquitous method to develop models for physical processes. Limiting what kinds of microstructural information is used as input and assessing normalized changes in network error, the relative importance of different mechanisms can be inferred. Here we utilize binned, initial density information as network inputs to predict local shock heating in a granular high explosive trained from large scale, molecular dynamics simulations. The spatial extend of the density field used in the network is altered to assess the importance and relevant length scales of the physical mechanisms in play, where different microstructural features result in different predictive capability.

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosive initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. We assess hotspots formed at shocked crystal-crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intra-molecular strain energy localization in the hotspot is assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5-20 nm of the crystal-crystal interface. Compressive work plays a considerably larger role in localizing temperature and intra-molecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal-crystal interfaces in coarse-grained models of polycrystalline high explosive materials.

Reactive force fields for molecular dynamics have enabled a wide range of studies in numerous material classes. These force fields are computationally inexpensive as compared to electronic structure calculations and allow for simulations of millions of atoms. However, the accuracy of traditional force fields is limited by their functional forms, preventing continual refinement and improvement. Therefore, we develop a neural network based reactive interatomic potential for the prediction of the mechanical, thermal, and chemical response of energetic materials at extreme conditions for energetic materials. The training set is expanded in an automatic iterative approach and consists of various CHNO materials and their reactions under ambient and under shock loading conditions. This new potential shows improved accuracy over the current state of the art force fields for a wide range of properties such as detonation performance, decomposition product formation, and vibrational spectra under ambient and shock loading conditions.

he manifestation of intra-molecular strains in covalent systems is widely known to accelerate chemical reactions and open alternative reaction paths. This process is moderately understood for isolated molecules and uni-molecular processes. However, in condensed matter processes such as phase transformations, material properties and structure may influence typical mechanochemical effects. Therefore, we utilize steered molecular dynamics to induce out of plane strains in graphite and compress the system under a constant strain rate to induce phase transformation. We show that the out of plane strain allows for phase transformations to initiate at lower amounts of compressive strain. Yet, in contrast to typical mechanochemical results, the sum of compressive and out of plane work needed to form diamond has a local minimum due to altered defect formation processes during phase transformation. Additionally, these altered processes slow the kinetics of the phase transformation, taking longer from initiation to total material transformation.

Mechanophores, stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing mechanophores. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics simulations of a PMMA-spiropyran co-polymer under dynamic mechanical loading and studied the activation of the mechanophore under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the mechanophore blocks experience significant many-body intra-molecular distortion that can significantly decrease the activation barrier as compared to when deformation rates are slow relative to molecular relaxation timescales. We find that the reactivity of mechanophores under material compression states is governed by many-body effects of intra-molecular torsions, whereas under tension the reactions are governed by tensile stresses.

Atomistic and continuum scale modeling efforts have shown that shock induced collapse of porosity can occur via a wide range of mechanisms dependent on pore morphology, shockwave pressure, and material properties. The mechanisms that occur under weaker shocks tend to be more efficient at localizing thermal energy, but do not result in high, absolute temperatures or spatially large localizations compared to mechanisms found under strong shock conditions. However, the energetic material TATB undergoes a wide range of collapse mechanisms that are not typical of similar materials, leaving the collapse mechanisms and the resultant energy localization from collapse, i.e., hotspots, relatively uncharacterized. Therefore, we present pore collapse simulations of cylindrical pores in TATB for a wide range of pore sizes and shock strengths that trigger viscoplastic collapses that occur almost entirely perpendicular to the shock direction for weak shocks and hydrodynamic-like collapses for strong shocks that do not break the strong hydrogen bonds of the TATB basal planes. The resulting hotspot temperature fields from these mechanisms follow trends that differ considerably from other energetic materials; hence we compare them under normalized temperature values to assess the relative efficiency of each mechanism to localize energy. The local intra-molecular strain energy of the hotspots is also assessed to better understand the physical mechanisms behind the phenomena that leads to a latent potential energy.

Regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shockwave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material TATB, we map molecular temperature and intramolecular strain energy prior to reaction to decomposition kinetics. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper can inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

The relaxation of polymers around and below their glass transition temperature is governed by a range of correlated unit processes with a wide range of timescales. The fast deformation rates of shock loading can negate a significant fraction of these processes resulting in the dynamical glass transition in rubbers. In this letter we report the inverse, a transient melting of glassy polymer under shock loading. The large deviatoric stresses near the shock front induce fast transitions in backbone dihedral angles and a stress relaxation characteristic of polymer melts. This is followed by the slower relaxation expected for glasses.

The shock to detonation transition in heterogeneous plastically bonded explosives is dominated by energy localization into hotspots that arise from the interaction of the shockwave with microstructural features and defects. The complex polycrystalline structure of these materials leads to a network of hotspot that can coalesce into deflagration and detonation waves. Significant progress has been made on the formation and potency of hotspots using atomistic simulations, but most of the work has focused on ideal and isolated defects. Hence, developed a method, denoted PBXGen, to build realistic PBX microstructures for all-atom simulations. PBXGen is generally applicable, and we demonstrate it with two systems: an RDX-polystyrene PBX with a 3D microstructure and a TATB-polystyrene with columnar grains. The resulting structure exhibit key features of PBXs, albeit at smaller scales, and are validated against experimental mechanical and shock properties.

In energetic materials, the localization of energy into "hotspots" when a shock wave interacts with the material's microstructure is known to dictate the initiation of chemical reactions and detonation. Recent results have shown that, following the shock-induced collapse of pores with circular cross-sections, more energy is localized as internal potential energy (PE) than can be inferred from the kinetic energy (KE) distribution. This leads to a complex thermo-mechanical state that is typically overlooked. The mechanisms associated with pore collapse and hotspot formation and the resulting energy localization are known to be highly dependent on material properties, especially its ability to deform plastically and alleviate strain energy, as well as the size and shape of the pore. Therefore, we use molecular dynamics simulations to characterize shock-induced pore collapse and the subsequent formation of hotspots in TATB, a highly anisotropic molecular crystal, for various defect shapes, shock strengths and crystallographic orientations. We find that the the localization of energy as PE is consistently higher and its extent larger than as localized as KE. A detailed analysis of the MD trajectories reveal the underlying molecular process that govern the effect of orientation and pore shape on the resulting hotspots. We find that the regions of highest PE for a given KE relate to not the impact front of the collapse, but the areas of maximum plastic deformation, while KE is maximized at the point of impact. A comparison with previous results in HMX reveal less energy localization in TATB which could be a contributing factor to its insensitivity.

Shockwave interactions with material microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists through thermal diffusion. The origin of the PE hotspot is traced to large intra-molecular strains, storing energy in modes readily available for chemical decomposition., Comment: Supplemental Material appended to bottom of manuscript document, 19 total pages, 8 total figures, manuscript is 12 pages 4 figures