Effects of molecular vacancies on the decomposition mechanisms and reaction dynamics of condensed-phase β-HMX at various temperatures were studied using ReaxFF molecular dynamics simulations. Results show that three primary initial decomposition mechanisms, namely, N-NO(2) bond dissociation, HONO elimination, and concerted ring fission, exist at both high and lower temperatures. The contribution of the three mechanisms to the initial decomposition of HMX is influenced by molecular vacancies, and the effects vary with temperature. At high temperature (2500 K), molecular vacancies remarkably promote N-N bond cleavage and concerted ring breaking but hinder HONO formation. N-N bond dissociation and HONO elimination are two primary competing reaction mechanisms, and the former is dominant in the initial decomposition. Concerted ring breaking of condensed-phase HMX is not favored at high temperature. At lower temperature (1500 K), the most preferential initial decomposition pathway is N-N bond dissociation followed by the formation of NO(3) (O migration), although all three mechanisms are promoted by molecular vacancies. The promotion effect on concerted ring breaking is considerable at lower temperature. Products resulting from concerted ring breaking appear in the defective system but not in the perfect crystal. The mechanism of HONO elimination is less important at lower temperature. We also estimated the reaction rate constant and activation barriers of initial decomposition with different vacancy concentrations. Molecular vacancies accelerate the decomposition of condensed-phase HMX by increasing the reaction rate constant and reducing activation barriers.

We report here the predictions on anisotropy of shock sensitivity and of chemical process initiation in single crystal b-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (b-HMX) using compressive shear reactive dynamics (CS-RD) model with ReaxFF-lg reactive force field. Analysis of resolved shear stress induced by uniaxial compression along three shock directions normal to (110), (011), and (010) planes leads to identify eight slip systems as candidates for shear deformation. For each of the eight slip systems, non-equilibrium reactive dynamics simulations were carried out to determine thermal, mechanical, and chemical responses to shear deformation. Shock direction normal to (010) plane exhibits large shear stress barriers arising from steric hindrance between molecules of adjacent layers leading to local dramatic energy and temperature increases under shear flow that in turn accelerate chemical bond breaking and initial product formation processes, promoting further molecular decomposition and eventually transition to detonation. This suggests that single crystal b-HMX is sensitive to shocks in direction normal to (010) plane. Shock directions normal to (110) and (011) planes reveal significantly less steric hindrance, leading to more modest energy and temperature increases followed by slower chemical reaction initiation. Thus, shock directions normal to (110) and (011) planes are less sensitive than shock direction normal to (010) plane, which agree with interpretations from currently available plate impact experiments on HMX. This validation of CS-RD and ReaxFF for characterizing sensitivity of single crystal energetic materials indicates that these methods can be applied to study sensitivity for more complex polymer bonded explosives and solid composite propellants having complex microstructures, corrugated interfaces, as well as defects. V

We report reactive molecular dynamics simulations using the ReaxFF reactive force field to examine shock-induced hot-spot formation followed by detonation initiation in realistic (2.7 million atoms) models of polymer bonded explosives (PBX) with nonplanar interfaces. We considered here two energetic materials (EMs) pentaerythritol tetranitrate (PETN), a common EM for PBX, and silicon pentaerythritol tetranitrate (Si-PETN), which is so extremely sensitive that it has not been possible to characterize its shock properties experimentally. In each case the EM was embedded in a hydroxyl-terminated polybutadiene (HTPB) based polymer binder matrix to form a model of PBX that has a periodic sawtooth nonplanar interface. For the cases in which the shock wave propagates from the EM to polymer (EM→poly), we observed that a hot spot arises from shear localization at the convex polymer asperity. For the case in which the shock direction is inverted (shock wave propagates from the polymer to the EM, EM←poly), we find that a hot spot is initiated at the concave polymer asperity and a second more significant hot spot forms at the convex polymer asperity. This second hot spot is enhanced due to converging shock wave interactions with the nonplanar interface. Under the same shock conditions, the first step in the Si-PETN decomposition is the Si−C−O−X rearrangement to Si−O−C−X through a five centered transition state on the Si that releases 45 kcal/mol of energy that leads to a continuous increase of temperature and pressure in the hot-spot region, until detonation. By contrast, the first step for PETN is NO 2 release, which is endothermic by 39 kcal/mol, with the consequence that the hot spot is attenuated by the polymer binder, reaching a steady temperature state involving NO 2 dissociation and HONO formation.

To gain an atomistic-level understanding of the thermal and chemical responses of condensed energetic materials under thermal shock, we developed a thermal shock reactive dynamics (TS-RD) computational protocol using molecular dynamics simulation coupled with ReaxFF force field. β-Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) was selected as a a target explosive due to its wide usage in the military and industry. The results show that a thermal shock initiated by a large temperature gradient between the "hot" region and the "cold" region results in thermal expansion of the particles and induces a thermal-mechanical wave propagating back and forth in the system with an averaged velocity of 3.32 km s(-1). Heat propagating along the direction of thermal shock leads to a temperature increment of the system and thus chemical reaction initiation. Applying a continuum reactive heat conduction model combined with the temperature distribution obtained from the RD simulation, a heat conduction coefficient is derived as 0.80 W m(-1) K(-1). The chemical reaction mechanisms during thermal shock were analyzed, showing that the reaction is triggered by N-NO2 bond breaking followed by HONO elimination and ring fission. The propagation rates of the reaction front and reaction center are obtained to be 0.069 and 0.038 km s(-1), based on the time and spatial distribution of NO2. The pressure effect on the thermal shock was also investigated by employing uniaxial compression before the thermal shock. We find that compression significantly accelerates thermal-mechanical wave propagation and heat conduction, resulting in higher temperature and more excited molecules and thus earlier initiation and faster propagation of chemical reactions.