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.

We report million-atom reactive molecular dynamic simulations of shock initiation of β-cyclotetramethylene tetranitramine (β-HMX) single crystals containing nanometer-scale spherical voids. Shock induced void collapse and subsequent hot spot formation as well as chemical reaction initiation are observed which depend on the void size and impact strength. For an impact velocity of 1 km s(-1) and a void radius of 4 nm, the void collapse process includes three stages; the dominant mechanism is the convergence of upstream molecules toward the centerline and the downstream surface of the void forming flowing molecules. Hot spot formation also undergoes three stages, and the principal mechanism is kinetic energy transforming to thermal energy due to the collision of flowing molecules on the downstream surface. The high temperature of the hot spot initiates a local chemical reaction, and the breakage of the N-NO2 bond plays the key role in the initial reaction mechanism. The impact strength and void size have noticeable effects on the shock dynamical process, resulting in a variation of the predominant mechanisms leading to void collapse and hot spot formation. Larger voids or stronger shocks result in more intense hot spots and, thus, more violent chemical reactions, promoting more reaction channels and generating more reaction products in a shorter duration. The reaction products are mainly concentrated in the developed hot spot, indicating that the chemical reactivity of the hmx crystal is greatly enhanced by void collapse. The detailed information derived from this study can aid a thorough understanding of the role of void collapse in hot spot formation and the chemical reaction initiation of explosives.

Selective excitation of C-H stretching vibrational modes, detection of intramolecular vibrational energy redistribution (IVR), and vibrational modes coupling in the electronic ground state of benzene are performed by using femtosecond time- and frequency-resolved coherent anti-Stokes Raman scattering (CARS) spectroscopy. Both of the parent modes in the Raman-active bands are coherently excited by an ultrafast stimulated Raman pump, giving initial excitations of 3056 cm (A) and 3074 cm (E) and subsequent IVR from the parent modes to daughter modes of 1181 and 992 cm, and the coherent vibrational coupling of the relevant modes is tracked. The directionality and selectivity of IVR and coherent coupling among all of the relevant vibrational modes are discussed in the view of molecular symmetry.

We employ DFT/B3LYP method to investigate linear open-chain clusters (n = 2-8) of the cis-triaziridine molecule that is a candidate molecule for high energy density materials (HEDM). Our calculations indicate that the pervasive phenomena of cooperative effects are observed in the clusters of n = 3-8, which are reflected in changes in lengths of N...H hydrogen bonds, stretching frequencies, and intensities of N-H bonds, dipole moments, and charge transfers as cluster size increases. The n(N) --> sigma*(N-H) interactions, i.e., the charge transfers from lone pairs (n(N)) of the N atoms into antibonds (sigma*) of the N-H bonds acting as H-donors, can be used to explain the observed cooperative phenomena. The approaches based upon natural bond orbital (NBO) method and theory of atoms in molecule (AIM) to evaluating N...H strengths are found to be equivalent. In the process of N...H bonding, cooperative nature of n(N) --> sigma*(N-H) interactions promotes formation of stronger N...H bonds as reflected in increases in the capacities of cis-triaziridine clusters to concentrate electrons at the bond critical points of N...H bonds. The calculated nonadditive energies also show that the cooperative effects due to n(N) --> sigma*(N-H) interactions indeed provide additional stabilities for the clusters.