Decomposition of excited electronic state s-tetrazine and its energetic derivatives

A Photophysical Deactivation Channel of Laser-Excited TATB Based on Semiclassical Dynamics Simulation and TD-DFT Calculation

A deactivation channel for laser-excited 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) was studied by semiclassical dynamics. Results indicate that the excited state resulting from an electronic transition from the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular mrbital (LUMO) is deactivated via pyramidalization of the activated N atom in a nitro group, with a lifetime of 2.4 ps. An approximately 0.5-electron transfer from the aromatic ring to the activated nitro group led to a significant increase of the C–NO₂ bond length, which suggests that C–NO₂ bond breaking could be a trigger for an explosive reaction. The time-dependent density functional theory (TD-DFT) method was used to calculate the energies of the ground and S₁ excited states for each configuration in the simulated trajectory. The S₁←S₀ energy gap at the instance of non-adiabatic decay was found to be 0.096 eV, suggesting that the decay geometry is close to the conical intersection.

Sub-500 fs electronically nonadiabatic chemical dynamics of energetic molecules from the S1 excited state: Ab initio multiple spawning study

Energetic materials store a large amount of chemical energy. Different ignition processes, including laser ignition and shock or compression wave, initiate the energy release process by first promoting energetic molecules to the electronically excited states. This is why a full understanding of initial steps of the chemical dynamics of energetic molecules from the excited electronic states is highly desirable. In general, conical intersection (CI), which is the crossing point of multidimensional electronic potential energy surfaces, is well established as a controlling factor in the initial steps of chemical dynamics of energetic molecules following their electronic excitations. In this article, we have presented different aspects of the ultrafast unimolecular relaxation dynamics of energetic molecules through CIs. For this task, we have employed ab initio multiple spawning (AIMS) simulation using the complete active space self-consistent field (CASSCF) electronic wavefunction and frozen Gaussian-based nuclear wavefunction. The AIMS simulation results collectively reveal that the ultrafast relaxation step of the best energetic molecules (which are known to exhibit very good detonation properties) is completed in less than 500 fs. Many, however, exhibit sub-50 fs dynamics. For example, nitro-containing molecules (including C-NO₂, N-NO₂, and O-NO₂ active moieties) relax back to the ground state in approximately 40 fs through similar (S₁/S₀)_CI conical intersections. The N₃-based energetic molecule undergoes the N₂ elimination process in 40 fs through the (S₁/S₀)_CI conical intersection. Nitramine-Fe complexes exhibit sub-50 fs Fe-O and N-O bond dissociation through the respective (S₁/S₀)_CI conical intersection. On the other hand, tetrazine-N-oxides, which are known to exhibit better detonation properties than tetrazines, undergo internal conversion in a 400-fs time scale, while the relaxation time of tetrazine is very long (about 100 ns). Many other characteristics of sub-500 fs nonadiabatic decay of energetic molecules are discussed. In the end, many unresolved issues associated with the ultrafast nonadiabatic chemical dynamics of energetic molecules are presented.

Vacuum ultraviolet photoionization of carbohydrates and nucleotides

Carbohydrates (2-deoxyribose, ribose, and xylose) and nucleotides (adenosine-, cytidine-, guanosine-, and uridine-5(')-monophosphate) are generated in the gas phase, and ionized with vacuum ultraviolet photons (VUV, 118.2 nm). The observed time of flight mass spectra of the carbohydrate fragmentation are similar to those observed [J.-W. Shin, F. Dong, M. Grisham, J. J. Rocca, and E. R. Bernstein, Chem. Phys. Lett. 506, 161 (2011)] for 46.9 nm photon ionization, but with more intensity in higher mass fragment ions. The tendency of carbohydrate ions to fragment extensively following ionization seemingly suggests that nucleic acids might undergo radiation damage as a result of carbohydrate, rather than nucleobase fragmentation. VUV photoionization of nucleotides (monophosphate-carbohydrate-nucleobase), however, shows that the carbohydrate-nucleobase bond is the primary fragmentation site for these species. Density functional theory (DFT) calculations indicate that the removed carbohydrate electrons by the 118.2 nm photons are associated with endocyclic C-C and C-O ring centered orbitals: loss of electron density in the ring bonds of the nascent ion can thus account for the observed fragmentation patterns following carbohydrate ionization. DFT calculations also indicate that electrons removed from nucleotides under these same conditions are associated with orbitals involved with the nucleobase-saccharide linkage electron density. The calculations give a general mechanism and explanation of the experimental results.

Identification of the free radical produced in the photolysis of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB)

1,3,5-Triamino-2,4,6-trinitrobenzene (TATB) is a typical insensitive high explosive (IHE) that possesses excellent heat, impact, and shock stability. However, it is sensitive to light irradiation, which can produce a long-lived free radical. In this study, (time-dependent) density functional theory is employed to study the features of the ground state (S0), the first singlet excited state (S1), and the first triplet excited state (T1). Results indicate the progress of photolysis, which involves an intersystem crossing from S1 to T1 followed by NO2-ONO isomerization. The long-lived radical produced in TATB photolysis is identified as the phenoxyl radical R-NO by investigating the formation feasibility and the stability. The experimental observation of the changes in the ultraviolet absorption spectra of TATB supports the identification.

Experimental and theoretical studies of the decomposition of new imidazole based energetic materials: model systems

Decomposition of three imidazole based model energetic systems (2-nitroimidazole, 4-nitroimidazole, and 1-methyl-5-nitroimidazole) is investigated both experimentally and theoretically. The initial decomposition mechanism for these three nitroimidazoles is explored with nanosecond energy resolved spectroscopy, and quantum chemical theory at the complete active space self-consistent field (CASSCF) level. The NO molecule is observed as an initial decomposition product from these three nitroimidazoles subsequent to UV excitation. A unique, excitation wavelength independent dissociation channel is observed for these three nitroimidazoles that generates the NO product with a rotationally cold (∼50 K) and a vibrationally mildly hot (∼800 K) distribution. Potential energy surface calculations at the CASSCF∕6-31G(d) level of theory illustrate that conical intersections play an important and essential role in the decomposition mechanism. Electronically excited S(2) nitroimidazole molecules relax to the S(1) state through the (S(2)∕S(1))(CI) conical intersection, and undergo a nitro-nitrite isomerization to generate the NO product from the S(1) potential energy surface. Nevertheless, NO(2) elimination and nitro-nitrite isomerization are expected to be competitive reaction mechanisms for the decomposition of these molecules on the ground state potential energy surface from the Franck-Condon equilibrium geometry through thermal dissociation.

On the origin of regio- and stereoselectivity in singlet oxygen addition to enecarbamates

The reactions of excited state singlet molecular oxygen ((1)Δ(g),(1)O(2)) continue to witness interesting new developments. In the most recent manifestation, (1)O(2) is tamed to react with enecarbamates in a stereoselective manner, which is remarkable, in view of its inherently high reactivity (Acc. Chem. Res. 2008, 41, 387). Herein, we employed the CAS-MP2(8,7)/6-31G* as well as the CAS-MP2(10,8)/6-31G* computations to unravel the origin of (i) diastereoselectivities in dioxetane or hydroperoxide formation and (ii) regioselectivity leading to a [2 + 2] cycloadduct or an ene product when (1)O(2) reacts with an oxazolidinone tethered 2-phenyl-1-propenyl system. The computed Gibbs free energy profiles for E- and Z-isomers when (1)O(2) approaches through the hindered and nonhindered diastereotopic faces (by virtue of chiral oxazolidinone) of the enecarbamates exhibit distinct differences. In the case of E-isomer, the relative energies of the transition structures responsible for hydroperoxide (ene product) are lower than that for dioxetane formation. On the other hand, the ene pathway is predicted to involve higher barriers as compared to the corresponding dioxetane pathway for Z-isomer. The energy difference between the rate-determining diastereomeric transition structures involved in the most favored ene reaction for E-enecarbamate suggests high diastereoselectivity. In contrast, the corresponding energy difference for Z-enecarbamate in the ene pathway is found to be diminishingly close, implying low diastereoselectivity. However, the dioxetane formation from Z-enecarbamate is predicted to exhibit high diastereoselectivity. The application of activation strain model as well as the differences in stereoelectronic effects in the stereocontrolling transition structures is found to be effective toward rationalizing the origin of selectivities reported herein. These predictions are found to be in excellent agreement with the experimental observations.