Shock-induced melting of (100)-oriented nitromethane: Energy partitioning and vibrational mode heating

High energy vibrational excitations of nitromethane in liquid water

The pathways and timescales of vibrational energy flow in nitromethane are investigated in both gas and condensed phases using classical molecular mechanics, with a particular focus on relaxation in liquid water. We monitor the flow of excess energy deposited in vibrational modes of nitromethane into the surrounding solvent. A marked energy flux anisotropy is found when nitromethane is immersed in liquid water, with a preferential flow to those water molecules in contact to the nitro group. The factors that permit such anisotropic energy relaxation are discussed, along with the potential implications on the molecule's non-equilibrium dynamics. In addition, the energy flux analysis allows us to identify the solvent motions responsible for the uptake of solute energy, confirming the crucial role of water librations. Finally, we also show that no anisotropic vibrational energy relaxation occurs when nitromethane is surrounded by argon gas.

Vibrational energy redistribution in crystalline nitromethane simulated by ab initio molecular dynamics

Ab initio molecular dynamics simulations (AIMD) are systematically performed to study the Vibrational Energy Redistribution (VER) in solid nitromethane (NM) by combining normal mode decomposition and short-time Fourier transform technique. After the selective excitations of all fourteen intramolecular vibrational modes above 400 cm⁻¹, four three-dimensional (3D) excitation and detected vibrational spectra are obtained. The evolution of the kinetic energy proportion of all vibrations are also given and discussed quantitatively. These results show that, as the daughter modes, NO₂ symmetric stretches, CH₃ stretches and bends are usually excited quickly and relatively conspicuously compared with the other vibrations. Interestingly, we found that, although the stretching vibration of the CN bond which is a bridge between the methyl and nitro group can not respond immediately to the selective excitations, it always accumulates the vibrational energy slowly and steadily. Then, the underlying mechanisms are discussed based on the response of vibrational modes in both the time and frequency domain. As a result, we found that anharmonic transfers following symmetry rules which involve the couplings assisted by the overtones and rotations, as well as the transfers among the adjacent modes, play important roles in the VER of solid NM.

Ab initio molecular dynamics simulations (AIMD) are systematically performed to study the Vibrational Energy Redistribution (VER) in solid nitromethane (NM) by combining normal mode decomposition and short-time Fourier transform technique.

Ab initio molecular dynamics simulation of vibrational energy redistribution of selective excitation of C-H stretching vibrations for solid nitromethane

Vibrational energy redistribution (VER) of energetic materials plays an important role in transferring the injected energy to the hot spots, but it is extremely challenging to understand the mechanism of VER from experimental or theoretical studies. Here, we combined nonequilibrium molecular dynamics with density functional theory to study the processes of VER for solid nitromethane after the selective excitation of the C-H stretching vibration. The VER processes are traced by monitoring the normal-mode kinetic energies of both excited and unexcited vibrations. To explore the underlying VER mechanism, we also analyzed the spectral energy density for the normal mode, obtained from the squared modulus of the short-time Fourier transition of their normal mode momentum. The results showed that the simulated VER progress was reproduced well compared with the previous 3D IR-Raman experiments of liquid nitromethane. Interestingly, the symmetric dependence of the coupling mechanism between the normal modes has been found.

Normal mode analysis on the relaxation of an excited nitromethane molecule in argon bath

In our previous work [Rivera-Rivera et al., J. Chem. Phys. 142, 014303 (2015)], classical molecular dynamics simulations followed the relaxation, in a 300 K Ar bath at a pressure of 10-400 atm, of nitromethane (CH₃NO₂) instantaneously excited by statistically distributing 50 kcal/mol among all its internal degrees of freedom. Both rotational and vibrational energies decayed with nonexponential curves. The present work explores mode-specific mechanisms at work in the decay process. With the separation of rotation and vibration developed by Rhee and Kim [J. Chem. Phys. 107, 1394 (1997)], one can show that the vibrational kinetic energy decomposes only into vibrational normal modes, while the rotational and Coriolis energies decompose into both vibrational and rotational normal modes. The saved CH₃NO₂ positions and momenta were converted into mode-specific energies whose decay was monitored over 1000 ps. The results identify vibrational and rotational modes that promote/resist energy lost and drive nonexponential behavior.

Controlling cross pumping between C-N and C-H vibration in nitromethane by selective fluorescence-enhanced stimulated Raman scattering

To investigate the vibrational features of nitromethane (NM), which is a kind of energy material and a well known low-sensitivity and high explosive, experiments are performed to obtain the stimulated Raman scattering (SRS) of NM by employing a 532 nm pulsed pump laser. The Raman signal involves two stimulated emissions at 918 and 2,963 cm^-1, attributed to the C-N and C-H stretching vibrations, respectively. To overcome the complexity of cross pump in the pure NM, one stimulated Raman radiation is chosen as a pump source to excite the other Raman mode. Two fluorescence dyes were added to selectively enhance each Raman cross section. By internally seeding the Raman gain medium with fluorescent photons, a significant modification in the stimulated Raman scattering spectrum has been observed. The enhanced Stokes emission at 918 cm^-1 was able to induce the 2,963 cm^-1 vibration mode when the all-trans-β-carotene was internal seeding in the NM, while the Raman radiation at 2,963 cm^-1 was enhanced to excite the C-N mode with the addition of m-Cresol purple. The output energy of both 918 and 2,963 cm^-1 under different input energy was also measured to illustrate this result.

Intermolecular Energy Transfer Dynamics at a Hot-Spot Interface in RDX Crystals

The phonon mediated vibrational up-pumping mechanisms assume an intact lattice and climbing of a vibrational ladder using strongly correlated multiphonon dynamics under equilibrium or near-equilibrium conditions. Important dynamic processes far from-equilibrium in regions of large temperature gradient after the onset of decomposition reactions in energetic solids are relatively unknown. In this work, we present a classical molecular dynamics (MD) simulation-based study of such processes using a nonreactive and a reactive potential to study a fully reacted and unreacted zone in RDX (1,3,5-trinitro-1,3,5-triazocyclohexane) crystal under nonequilibrium conditions. The energy transfer rate is evaluated as a function of temperature difference between the reacted and unreacted regions, and for different widths and cross-sectional area of unreacted RDX layers. Vibrational up-pumping processes probed using velocity autocorrelation functions indicate that the mechanisms at high-temperature interfaces are quite different from the standard phonon-based models proposed in current literature. In particular, the up-pumping of high-frequency vibrations are seen in the presence of small molecule collisions at the hot-spot interface with strong contributions from bending modes. It also explains some major difference in the order of decomposition of C-N and N-N bonds as seen in recent literature on initiation chemistry.

Prediction of the Chapman–Jouguet chemical equilibrium state in a detonation wave from first principles based reactive molecular dynamics

This Rx2CJ first principle based protocol for predicting the CJ state provides the matching point between atomistic reaction dynamic simulations and the macroscopic properties of detonation, and can be used as a measure of performance for in silico synthesis and characterization of new materials.

Understanding the desensitizing mechanism of olefin in explosives: shear slide of mixed HMX-olefin systems

We simulated the shear slide behavior of typical mixed HMX-olefin systems and the effect of thickness of olefin layers (4-22 Å) on the behavior at a molecular level by considering two cases: bulk shear and interfacial shear. The results show that: (1) the addition of olefin into HMX can reduce greatly the shear sliding barriers relative to the pure HMX in the two cases, suggesting that the desensitizing mechanism of olefin is controlled dominantly by its good lubricating property; (2) the change of interaction energy in both systoles of shear slide is strongly dominated by van der Waals interaction; and (3) the thickness of olefin layers in the mixed explosives can influence its desensitizing efficiency. That is, the excessive thinness of olefin layers in the mixed explosive systems, for example, several angstroms, can lead to very high sliding barriers.