Reaction Rates in Nitromethane under High Pressure from Density Functional Tight Binding Molecular Dynamics Simulations.

We use density functional tight binding (DFTB) molecular dynamics (MD) simulations to determine the reaction rates of nitromethane CH3NO2 (NM) under high pressure (P = 14-28 GPa), and temperature (T = 1450-1850 K). DFTB-MD simulations performed with the same initial conditions (P0, T0) reveal a stochastic behavior, both in terms of reaction times and chemical paths. By running series of MD simulations, we are able to obtain average reaction times with quantified errors and devise a simple two-step model for NM explosion: ignition/explosion. While our model bypasses the chemical complexity due to the numerous reaction paths and intermediates observed during reactions, the chemistry is accounted for via the accurate parameterization of the DFTB model, and our results suggest a single main reaction pathway for the pressure range considered here, dominated in the earlier stages by the formation of the aci-ion, CH2NOO-. By fitting our data to a Frank-Kamenetskii model, we extract prefactors and pressure-independent activation energies and volumes for the ignition and explosion stages. A two-step model is then built and compared to experimental observations. Single and two-step Arrhenius models are also provided for comparison with literature data. This work presents an efficient way of investigating the reactivity of high explosives by performing electronic structure-based MD simulations and provides reaction rates for simplified models that can be implemented into hydrocodes.

Probing ultrafast shock-induced chemistry in liquids using broad-band mid-infrared absorption spectroscopy.

We probe shock-induced chemistry in two organic liquids by measuring broadband, midinfrared absorption in the 800-1400 cm-1 frequency range. To test this new method and understand the signatures of chemical reactions in time resolved vibrational spectra, we compared liquid benzene shocked to unreactive conditions (shocked to a pressure of 18 GPa for a duration of 300 ps) to nitromethane under reactive conditions (25 GPa). We see clear signatures of shock-induced chemistry that are distinguishable from the pressure- and temperature-induced changes in vibrational mode shapes. While shocked benzene shows primarily a broadening and shifting of the vibrational modes, the nitromethane vibrational modes vanish once the shock wave enters the liquid and simultaneously, a spectrally broad feature appears that we interpret as the infrared spectrum of the complex mixture of product and intermediate species. To further interpret these measurements, we compare them to reactive quantum molecular dynamics simulations, which gives qualitatively consistent results. This work demonstrates a promising method for time resolving shock-induced chemistry, illustrating that chemical reactions produce distinct changes in the vibrational spectra.

Ultrafast Photodissociation Dynamics of Nitromethane.

Nitromethane (NM), a high explosive (HE) with low sensitivity, is known to undergo photolysis upon ultraviolet (UV) irradiation. The optical transparency, homogeneity, and extensive study of NM make it an ideal system for studying photodissociation mechanisms in conventional HE materials. The photochemical processes involved in the decomposition of NM could be applied to the future design of controllable photoactive HE materials. In this study, the photodecomposition of NM from the nπ* state excited at 266 nm is being investigated on the femtosecond time scale. UV femtosecond transient absorption (TA) spectroscopy and excited state femtosecond stimulated Raman spectroscopy (FSRS) are combined with nonadiabatic excited state molecular dynamics (NA-ESMD) simulations to provide a unified picture of NM photodecomposition. The FSRS spectrum of the photoproduct exhibits peaks in the NO2 region and slightly shifted C-N vibrational peaks pointing to methyl nitrite formation as the dominant photoproduct. A total photolysis quantum yield of 0.27 and an nπ* state lifetime of ∼20 fs were predicted from NA-ESMD simulations. Predicted time scales revealed that NO2 dissociation occurs in 81 ± 4 fs and methyl nitrite formation is much slower having a time scale of 452 ± 9 fs corresponding to the excited state absorption feature with a decay of 480 ± 17 fs observed in the TA spectrum. Although simulations predict C-N bond cleavage as the primary photochemical process, the relative time scales are consistent with isomerization occurring via NO2 dissociation and subsequent rebinding of the methyl radical and nitrogen dioxide.

Advances in explosives analysis--part II: photon and neutron methods.

The number and capability of explosives detection and analysis methods have increased dramatically since publication of the Analytical and Bioanalytical Chemistry special issue devoted to Explosives Analysis [Moore DS, Goodpaster JV, Anal Bioanal Chem 395:245-246, 2009]. Here we review and critically evaluate the latest (the past five years) important advances in explosives detection, with details of the improvements over previous methods, and suggest possible avenues towards further advances in, e.g., stand-off distance, detection limit, selectivity, and penetration through camouflage or packaging. The review consists of two parts. Part I discussed methods based on animals, chemicals (including colorimetry, molecularly imprinted polymers, electrochemistry, and immunochemistry), ions (both ion-mobility spectrometry and mass spectrometry), and mechanical devices. This part, Part II, will review methods based on photons, from very energetic photons including X-rays and gamma rays down to the terahertz range, and neutrons.

Advances in explosives analysis--part I: animal, chemical, ion, and mechanical methods.

The number and capability of explosives detection and analysis methods have increased substantially since the publication of the Analytical and Bioanalytical Chemistry special issue devoted to Explosives Analysis (Moore and Goodpaster, Anal Bioanal Chem 395(2):245-246, 2009). Here we review and critically evaluate the latest (the past five years) important advances in explosives detection, with details of the improvements over previous methods, and suggest possible avenues towards further advances in, e.g., stand-off distance, detection limit, selectivity, and penetration through camouflage or packaging. The review consists of two parts. This part, Part I, reviews methods based on animals, chemicals (including colorimetry, molecularly imprinted polymers, electrochemistry, and immunochemistry), ions (both ion-mobility spectrometry and mass spectrometry), and mechanical devices. Part II will review methods based on photons, from very energetic photons including X-rays and gamma rays down to the terahertz range, and neutrons.

Photoactive high explosives: linear and nonlinear photochemistry of petrin tetrazine chloride.

Pentaerythritol tetranitrate (PETN), a high explosive, initiates with traditional shock and thermal mechanisms. In this study, the tetrazine-substituted derivative of PETN, pentaerythritol trinitrate chlorotetrazine (PetrinTzCl), is being investigated for a photochemical initiation mechanism that could allow control over the chemistry contributing to decomposition leading to initiation. PetrinTzCl exhibits a photochemical quantum yield (QYPC) at 532 nm not evident with PETN. Using static spectroscopic methods, we observe energy absorption on the tetrazine (Tz) ring that results in photodissociation yielding N2, Cl-CN, and Petrin-CN as the major photoproducts. The QYPC was enhanced with increasing irradiation intensity. Experiment and theoretical calculations imply this excitation mechanism follows sequential photon absorption. Dynamic simulations demonstrate that the relaxation mechanism leading to the observed photochemistry in PetrinTzCl is due to vibrational excitation during internal conversion. PetrinTzCl's single photon stability and intensity dependence suggest this material could be stable in ambient lighting, yet possible to initiate with short-pulsed lasers.

Quantitative tradeoffs between spatial, temporal, and thermometric resolution of nonresonant Raman thermometry for dynamic experiments.

The ratio of Stokes to anti-Stokes nonresonant spontaneous Raman can provide an in situ thermometer that is noncontact, independent of any material specific parameters or calibrations, can be multiplexed spatially with line imaging, and can be time resolved for dynamic measurements. However, spontaneous Raman cross sections are very small, and thermometric measurements are often limited by the amount of laser energy that can be applied without damaging the sample or changing its temperature appreciably. In this paper, we quantitatively detail the tradeoff space between spatial, temporal, and thermometric accuracy measurable with spontaneous Raman. Theoretical estimates are pinned to experimental measurements to form realistic expectations of the resolution tradeoffs appropriate to various experiments. We consider the effects of signal to noise, collection efficiency, laser heating, pulsed laser ablation, and blackbody emission as limiting factors, provide formulae to help choose optimal conditions and provide estimates relevant to planning experiments along with concrete examples for single-shot measurements.

Ultrafast chemical reactions in shocked nitromethane probed with dynamic ellipsometry and transient absorption spectroscopy.

Initiation of the shock driven chemical reactions and detonation of nitromethane (NM) can be sensitized by the addition of a weak base; however, the chemical mechanism by which sensitization occurs remains unclear. We investigated the shock driven chemical reaction in NM and in NM sensitized with diethylenetriamine (DETA), using a sustained 300 ps shock driven by a chirped Ti:sapphire laser. We measured the solutions' visible transient absorption spectra and measured interface particle and shock velocities of the nitromethane solutions using ultrafast dynamic ellipsometry. We found there to be a volume-increasing reaction that takes place around interface particle velocity up = 2.4 km/s and up = 2.2 km/s for neat NM and NM with 5% DETA, respectively. The rate at which transient absorption increases is similar in all mixtures, but with decreasing induction times for solutions with increasing DETA concentrations. This result supports the hypothesis that the chemical reaction mechanisms for shocked NM and NM with DETA are the same. Data from shocked NM are compared to literature experimental and theoretical data.

Shock Hugoniot equations of state for binary ideal (toluene/fluorobenzene) and nonideal (ethanol/water) liquid mixtures.

Laser shock Hugoniot data were obtained using ultrafast dynamic ellipsometry (UDE) for both nonideal (ethanol/water solutions with mole percent χ(ethanol) = 0%, 3.4%, 5.4%, 7.5%, 9.7%, 11%, 18%, 33%, 56%, 100%) and ideal liquid mixtures (toluene/fluorobenzene solutions with mole percent χ(toluene) = 0%, 26.0%, 49.1%, 74.9%, 100%). The shock and particle velocities obtained from the UDE data were compared to the universal liquid Hugoniot (ULH) and to literature shock (plate impact) data where available. It was found that the water UDE data fit to a ULH-form equation suggests an intercept of 1.32 km/s, lower than the literature ambient sound speed in water of 1.495 km/s (Mijakovic et al. J. Mol. Liq. 2011, 164, 66-73). Similarly, the ethanol UDE data fit to a ULH-form equation suggests an intercept of 1.45 km/s, which lies above the literature ambient sound speed in ethanol of 1.14 km/s. Both the literature plate impact and UDE Hugoniot data lie below the ULH for water. Likewise, the literature plate impact and UDE Hugoniot data lie above the ULH for ethanol. The UDE Hugoniot data for the mixtures of water and ethanol cross the predictions of the ULH near the same concentration where the sound speed reaches a maximum. In contrast, the UDE data from the ideal liquids and their mixtures are well behaved and agree with ULH predictions across the concentration range. The deviations of the nonideal ethanol/water data from the ULH suggest that complex hydrogen bonding networks in ethanol/water mixtures alter the compressibility of the mixture.

Quantum chemistry studies of electronically excited nitrobenzene, TNA, and TNT.

The electronic excitation energies and excited-state potential energy surfaces of nitrobenzene, 2,4,6-trinitroaniline (TNA), and 2,4,6-trinitrotoluene (TNT) are calculated using time-dependent density functional theory and multiconfigurational ab initio methods. We describe the geometrical and energetic character of excited-state minima, reaction coordinates, and nonadiabatic regions in these systems. In addition, the potential energy surfaces for the lowest two singlet (S(0) and S(1)) and lowest two triplet (T(1) and T(2)) electronic states are investigated, with particular emphasis on the S(1) relaxation pathway and the nonadiabatic region leading to radiationless decay of S(1) population. In nitrobenzene, relaxation on S(1) occurs by out-of-plane rotation and pyramidalization of the nitro group. Radiationless decay can take place through a nonadiabatic region, which, at the TD-DFT level, is characterized by near-degeneracy of three electronic states, namely, S(1), S(0), and T(2). Moreover, spin-orbit coupling constants for the S(0)/T(2) and S(1)/T(2) electronic state pairs were calculated to be as high as 60 cm(-1) in this region. Our results suggest that the S(1) population should quench primarily to the T(2) state. This finding is in support of recent experimental results and sheds light on the photochemistry of heavier nitroarenes. In TNT and TNA, the dominant pathway for relaxation on S(1) is through geometric distortions, similar to that found for nitrobenzene, of a single ortho-substituted NO(2). The two singlet and lowest two triplet electronic states are qualitatively similar to those of nitrobenzene along a minimal S(1) energy pathway.

Anharmonic vibrational properties of explosives from temperature-dependent Raman.

Raman spectra from 50 to 3500 cm(-1) and 4-296 K are analyzed for molecular crystal powders of the explosives pentaerythritol tetranitrate (PETN), beta-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX), and 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and the inert naphthalene. Temperature-dependent Raman spectroscopy is utilized for its sensitivity to anharmonic couplings between thermally populated phonons and higher frequency vibrations relevant to shock up-pumping. The data are analyzed with anharmonic perturbation theory, which is shown to have significant fundamental limitations in application to real data. Fitting to perturbation theory revealed no significant differences in averaged anharmonicities among the three explosives, all of which exhibited larger averaged anharmonicities than naphthalene by a factor of 3. Calculations estimating the multiphonon densities of states also failed to correlate clearly with shock sensitivity. However, striking differences in temperature-dependent lifetimes were obvious: PETN has long lived phonons and vibrons, HMX has long lived phonons but short lived vibrons, while TATB has short lived phonons and vibrons at low temperature. Naphthalene, widely used as a model system, has significantly different anharmonicities and density of states from any of the explosives. The data presented suggest the further hypothesis that hindered vibrational energy transfer in the molecular crystals is a significant factor in shock sensitivity.