Ultrafast Photodissociation Dynamics of Nitromethane

Metal-Determined Explosive Characteristics of M(NO3)2(1-AT)x of Thermal and Laser Ignition (M2+ = Cr2+, Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+; x = 2 or 3)

Optical and thermal ignition are two common pathways to initiate the explosion of primary explosives, where laser ignition is a more reliable and safer initiation method. Caused by the current-applied laser igniter with the wavelength of 1064 or 915 nm, the energetic complexes with strong absorption in the near-infrared (NIR) region are possibly applied as laser-ignited explosives. Recently, [Cu(NO3)2(1-AT)3] complex has been synthesized with excellent NIR absorption properties, where 1-amino-5H-tetrazole (1-AT) has been proved to be a promising laser-ignited energetic ligand. To confirm the structure-thermal/optical explosive characteristics, based on the structure of synthesized [Cu(NO3)2(1-AT)3], the commonly used transition-metal cations (M2+ = Cr2+, Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+) have been selected to construct the series of complexes of M(NO3)2(1-AT)x (x = 2 or 3) theoretically. Car-Parrinello molecular dynamics (CPMD) method has been applied to unveil the role of center metals in the initiation and growth pathways. Time-dependent density functional theory (TD-DFT) method is used to explore their charge-transfer (CT) characteristics. The optical characteristic of the metal complex is mainly determined by the behaviors of the 3d electrons of center metals in excitation, where the activity of β-d electrons is an important factor to affect the NIR characteristic of complexes.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

On-the-Fly Nonadiabatic Dynamics Simulations of Single-Walled Carbon Nanotubes with Covalent Defects

Single-walled carbon nanotubes (SWCNTs) with covalent surface defects have been explored recently due to their promise for use in single-photon telecommunication emission and in spintronic applications. The all-atom dynamic evolution of electrostatically bound excitons (the primary electronic excitations) in these systems has only been loosely explored from a theoretical perspective due to the size limitations of these large systems (>500 atoms). In this work, we present computational modeling of nonradiative relaxation in a variety of SWCNT chiralities with single-defect functionalizations. Our excited-state dynamics modeling uses a trajectory surface hopping algorithm accounting for excitonic effects with a configuration interaction approach. We find a strong chirality and defect-composition dependence on the population relaxation (varying over 50-500 fs) between the primary nanotube band gap excitation E₁₁ and the defect-associated, single-photon-emitting E₁₁* state. These simulations give direct insight into the relaxation between the band-edge states and the localized excitonic state, in competition with dynamic trapping/detrapping processes observed in experiment. Engineering fast population decay into the quasi-two-level subsystem with weak coupling to higher-energy states increases the effectiveness and controllability of these quantum light emitters.

Triggering the mechanism of the initial reaction of energetic materials under pressure based on Raman intensity analysis

The Raman intensity and other stoichiometric calculations of nitromethane (NM) and 2-nitrimino-5-nitro-hexahydro-1,3,5-triazine (NNHT) have been made by using first-principles density functional theory. We propose a method to judge the initial reaction mechanism of NM and NNHT under pressure based on the Raman intensity. Both the resulting NM and NNHT undergo hydrogen transfer and conventional trigger bond cleavage. And the results obtained from the Raman peak intensities infer a reaction path that is not inferior to the traditional C-NO₂ and N-NO₂ bond cleavage, thus verifying our results.

Electronic State Spectroscopy of Nitromethane and Nitroethane

High-resolution photoabsorption cross-sections in the 3.7-10.8 eV energy range are reinvestigated for nitromethane (CH₃NO₂), while for nitroethane (C₂H₅NO₂), they are reported for the first time. New absorption features are observed for both molecules which have been assigned to vibronic excitations of valence, Rydberg, and mixed valence-Rydberg characters. In comparison with nitromethane, nitroethane shows mainly broad absorption bands with diffuse structures, which can be interpreted as a result of the side-chain effect contributing to an increased number of internal degrees of freedom. New theoretical quantum chemical calculations performed at the time-dependent density functional theory (TD-DFT) level were used to qualitatively help interpret the recorded photoabsorption spectra. From the photoabsorption cross-sections, photolysis lifetimes in the terrestrial atmosphere have been obtained for both compounds. Relevant internal conversion from Rydberg to valence character is noted for both molecules, while the nuclear dynamics of CH₃NO₂ and C₂H₅NO₂ along the C-N reaction coordinate have been evaluated through potential energy curves at the TD-DFT level of theory, showing that the pre-dissociative character is more prevalent in nitromethane than in nitroethane.

Modeling Photolytic Decomposition of Energetically Functionalized Dodecanes

The photolytic stability of explosives and energetic functional groups is of importance for those who regularly handle or are exposed to explosives in typical environmental conditions. This study models the photolytic degradation of dodecane substituted with various energetic functional groups: azide, nitro, nitrate ester, and nitramine. For the studied molecules, it was found that excitons localize on the energetic functional group, no matter where they were initially formed, and thus, the predominant degradation pathway involves the degradation of the energetic functional group. The relative trends for both 4 and 8 eV excitation energies followed with what is expected from the relative stability of the energetic functional groups to thermal and sub-shock degradation. The one notable exception was the azide functional group; more work should be done to further understand the photolytic effects on the azide functional group.

Ultrafast Dynamics of Nitro-Nitrite Rearrangement and Dissociation in Nitromethane Cation

We report new insights into the ultrafast rearrangement and dissociation dynamics of nitromethane cation (NM⁺) using pump-probe measurements, electronic structure calculations, and ab initio molecular dynamics simulations. The "roaming" nitro-nitrite rearrangement (NNR) pathway involving large-amplitude atomic motion, which has been previously described for neutral nitromethane, is demonstrated for NM⁺. Excess energy resulting from initial population of the electronically excited D₂ state of NM⁺ upon strong-field ionization provides the necessary energy to initiate NNR and subsequent dissociation into NO⁺. Both pump-probe measurements and molecular dynamics simulations are consistent with the completion of NNR within 500 fs of ionization with dissociation into NO⁺ and OCH₃ occurring ∼30 fs later. Pump-probe measurements indicate that NO⁺ formation is in competition with the direct dissociation of NM⁺ to CH₃⁺ and NO₂. Electronic structure calculations indicate that a strong D₀ → D₁ transition can be excited at 650 nm when the C-N bond is stretched from its equilibrium value (1.48 Å) to 1.88 Å. On the other hand, relaxation of the NM⁺ cation after ionization into D₀ occurs in less than 50 fs and results in observation of intact NM⁺. Direct dissociation of the equilibrium NM⁺ to produce NO₂⁺ and CH₃ can be induced with 650 nm excitation via a weakly allowed D₀ → D₂ transition.

New perspectives on the laser initiation for metal tetrazine complexes: a theoretical study

In order to understand the relationship between laser initiation and charge transfer of metal tetrazine complexes (MTCs), several sets of MTCs with different metals and ligands were designed and their charge transfer (CT) characters were examined using a time-dependent density functional theory method (TD-DFT) in combination with UV-vis spectra, hole-electron distribution, interfragment charge transition, and transition density matrix analyses. Results show that Fe(II), Mn(II), and Cu(II) are suitable divalent transition metal cores in constructing the optical initiation tetrazine complexes. By replacing the divalent metal cores with a monovalent center, new sets of complexes are proved to possess metal-to-ligand charge transfer (MLCT) character and stronger absorption intensity in the near-infrared (NIR) region, which implies that monovalent MTCs are more in favor of low-energy laser initiation than divalent MTCs. Reasonable tuning of the structure of pyrazole substituent can expect to enhance the explosive performance while preserving the optical characteristics, which is an important design principle. This work thoroughly depicts the photoactive states for MTCs and gives a train of thought to explore new desirable laser initiation explosives.

Femtosecond photodecarbonylation of photo-ODIBO studied by stimulated Raman spectroscopy and density functional theory

Photo-oxa-dibenzocyclooctyne (Photo-ODIBO) undergoes photodecarbonylation under UV excitation to its bright S₂ state, forming a highly reactive cyclooctyne, ODIBO. Following 321 nm excitation with sub-50 fs actinic pulses, the excited state evolution and cyclopropenone bond cleavage with CO release were characterized using femtosecond stimulated Raman spectroscopy and time-dependent density functional theory Raman calculations. Analysis of the photo-ODIBO S₂ CO Raman band revealed multi-exponential intensity, peak splitting and frequency-shift dynamics. This suggests a stepwise cleavage of the two C-C bonds in the cyclopropenone structure that is completed within <300 fs after excitation. Evidence of intramolecular vibrational relaxation on the S₂ state, concurrent with photodecarbonylation, with dynamics matching previous electronic transient absorption spectroscopy, was also observed. This confirms an excited state, as opposed to ground state, photodecarbonylation mechanism resulting in a vibronically excited photoproduct, ODIBO.

Dissociation of Singly and Multiply Charged Nitromethane Cations: Femtosecond Laser Mass Spectrometry and Theoretical Modeling

Dissociation pathways of singly- and multiply charged gas-phase nitromethane cations were investigated with strong-field laser photoionization mass spectrometry and density functional theory computations. There are multiple isomers of the singly charged nitromethane radical cation, several of which can be accessed by rearrangement of the parent CH₃-NO₂ structure with low energy barriers. While direct cleavage of the C-N bond from the parent nitromethane cation produces NO₂⁺ and CH₃⁺, rearrangement prior to dissociation accounts for fragmentation products including NO⁺, CH₂OH⁺, and CH₂NO⁺. Extensive Coulomb explosion in fragment ions observed at high laser intensity indicates that rapid dissociation of multiply charged nitromethane cations produces additional species such as CH₂⁺, H⁺, and NO₂²⁺.  On the basis of analysis of Coulomb explosion in the mass spectral signals and pathway calculations, sufficiently intense laser fields can remove four or more electrons from nitromethane.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Auto-ignition behaviors of nitromethane in diluted oxygen in a rapid compression machine: Critical conditions for ignition, ignition delay times measurements, and kinetic modeling interpretation

In this work, the weakest thermodynamic conditions for the auto-ignition of mixtures containing nitromethane were experimentally determined by using the rapid compression machine facility. Results show there is a narrow weak ignition region between ignition and non-ignition. The weak ignition region would disappear with the increase of the EOC (end of compression) pressure and nitromethane concentration. In addition, the ignition delay times for successful auto-ignition for different nitromethane concentrations and equivalence ratio mixtures were measured and compared. Results show that the dependence of nitromethane ignition on the equivalence ratio is weak. Subsequently, the measured ignition delay time data were employed to validate several kinetic models in literature and our previous model shows better agreement with experimental results, as well as other available literature data. Sensitivity analysis for the model reveals the importance of unimolecular decomposition and H-abstraction reactions for the ignition delay times in the temperature range studied herein. Finally, critical conditions for nitromethane ignition under extended conditions that are beyond the ability of the experimental facility were predicted.

Unraveling Photodimerization of Cyclohexasilane from Molecular Dynamics Studies

Photoinduced reactions of a pair of cyclohexasilane (CHS) monomers are explored by time-dependent excited-state molecular dynamics (TDESMD) calculations. In TDESMD trajectories, one observes vivid reaction events including dimerization and fragmentation. A general reaction pathway is identified as (i) ring-opening formation of a dimer, (ii) rearrangement induced by bond breaking, and (iii) decomposition through the elimination of small fragments. The identified pathway supports the chemistry proposed for the fabrication of silicon-based materials using CHS as a precursor. In addition, we find dimers have smaller HOMO-LUMO gaps and exhibit a red shift and line-width broadening in the computed photoluminescence spectra compared with a pair of CHS monomers.

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.

Unravelling the Electronic State of NO2 Product in Ultrafast Photodissociation of Nitromethane

The primary photochemical reaction of nitromethane (NM) after ππ* excitation is known to be C-N bond cleavage (CH₃NO₂ + hν → CH₃ + NO₂). On the other hand, NO₂ can be formed in both the ground and excited states, and identification of the electronic state of the NO₂ product has been a central subject in the experimental and theoretical studies. Here we present time-resolved photoelectron spectroscopy using vacuum-ultraviolet probe pulses to observe all transient electronic states of NM and the reaction products. The result indicates that ultrafast internal conversion occurs down to S₁ and S₀ within 24 fs, and the dissociation proceeds on the S₁ surface (τ_diss ≲ 50 fs), leading to comparable product yields of NO₂(A) and NO₂(X). The overall dissociation quantum yield within our observation time window (<2 ps) is estimated to be 0.29.

Photofragmentation of Tetranitromethane: Spin-Unrestricted Time-Dependent Excited-State Molecular Dynamics

In this study, the photofragmentation dynamics of tetranitromethane (TNM) is explored by a spin-unrestricted time-dependent excited-state molecular dynamics (u-TDESMD) algorithm based on Rabi oscillations and principles similar to trajectory surface hopping, with a midintensity field approximation. The leading order process is represented by the molecule undergoing cyclic excitations and de-excitations. During excitation cycles, the nuclear kinetic energy is accumulated to overcome the dissociation barriers in the reactant and a sequence of intermediates. The dissociation pathway includes the ejection of NO₂ groups followed by the formation of NO and CO. The simulated mass spectra at the ab initio level, based on the bond length in possible fragments, are extracted from simulation trajectories. The recently developed methodology has the potential to model and monitor photoreactions with open-shell intermediates and radicals.

Cooperative enhancement of the nonlinear optical response in conjugated energetic materials: A TD-DFT study

Conjugated energetic molecules (CEMs) are a class of explosives with high nitrogen content that posses both enhanced safety and energetic performance properties and are ideal for direct optical initiation. As isolated molecules, they absorb within the range of conventional lasers. Crystalline CEMs are used in practice, however, and their properties can differ due to intermolecular interaction. Herein, time-dependent density functional theory was used to investigate one-photon absorption (OPA) and two-photon absorption (TPA) of monomers and dimers obtained from experimentally determined crystal structures of CEMs. OPA scales linearly with the number of chromophore units, while TPA scales nonlinearly, where a more than 3-fold enhancement in peak intensity, per chromophore unit, is calculated. Cooperative enhancement depends on electronic delocalization spanning both chromophore units. An increase in sensitivity to nonlinear laser initiation makes these materials suitable for practical use. This is the first study predicting a cooperative enhancement of the nonlinear optical response in energetic materials composed of relatively small molecules. The proposed model quantum chemistry is validated by comparison to crystal structure geometries and the optical absorption of these materials dissolved in solution.