Photodissociation dynamics of nitromethane and methyl nitrite by infrared multiphoton dissociation imaging with quasiclassical trajectory calculations: Signatures of the roaming pathway

Matrix-formation dynamics dictate methyl nitrite conformer abundance

Methyl nitrite has two stable conformational isomers resulting from rotation about the primary C-O-N-O dihedral angle: cis-CH3ONO and trans-CH3ONO, with cis being more stable by ∼5 kJ/mol. The barrier to rotational interconversion (∼45 kJ/mol) is too large for isomerization to occur under ambient conditions. This paper presents evidence of a change in conformer abundance when dilute CH3ONO is deposited onto a cold substrate; the relative population of the freshly deposited cis conformer is seen to increase compared to its gas-phase abundance, measured by in situ infrared spectroscopy. We observe abundance changes depending on the identity of the bath gas (N2, Ar, and Xe) and deposition angle. The observations indicate that the surface properties of the growing matrix influence conformer abundance-contrary to the widely held assumption that conformer abundance in matrices reflects gas-phase abundance. We posit that differences in the angle-dependent host-gas deposition dynamics affect the growing surfaces, causing changes in conformer abundances. Quantum chemistry calculations of the binding energies between CH3ONO and a single bath-gas component reveal that significant energetic stabilization is not observed in 1:1 complexes of N2:CH3ONO, Ar:CH3ONO, or Xe:CH3ONO. From our results, we conclude that the growing surface plays a significant role in trapping cis-CH3ONO more effectively than trans-CH3ONO, likely because cis-CH3ONO is more compact. Taken together, the observations highlight the necessity for careful characterization of conformers in matrix-isolated systems, emphasizing a need for further study into the deposition dynamics and surface structure of chemically inert matrices.

Unusual In-plane Aromaticity Facilitates Intramolecular Hydrogen Transfer in Long-Bonded cis-Isonitrosyl Methoxide

The hydrogen-atom transfer from methoxy radical to nitric oxide, leading to the formation of formaldehyde and nitroxyl, represents a secondary reaction of photodissociation of methyl nitrite, which is used as rocket fuel. In this study, we explored the potential energy profile of the hydrogen-atom transfer using the electronic structure calculations at the DLPNO-CCSD(T)/aug-cc-pVTZ level of theory for two isomeric forms (cis and trans) of the pre-reaction complex. The cis-oriented pre-reaction complex has a weak elongated O─O bond, which gets further elongated in the hydrogen transfer transition state. This O─O bond stabilizes the pre-reaction complex by 32.9 kJ/mol. The O─O-induced stabilization is even greater for the transition state (48.2 kJ/mol), which was unexpected because of the larger O─O distance in the transition state structure. To address this paradox, we performed the electronic structure analysis of the reaction participants using the valence bond (VB) theory, natural resonance theory, topological analysis of the electron density and its derivatives, and analysis of the electron localization function distribution. This combined analysis led to the conclusion that the cis-transition state for hydrogen transfer, instead of being directly stabilized by the O─O interaction, gained substantial stabilization from the in-plane five-center six-electron aromaticity.

Insights into non-thermal plasma chemistry of acetone diluted in N2/O2 mixtures: a real-time MS experiment

Understanding non-thermal plasma reactivity is a complicated task as many reactions take place due to a large energy spectrum. In this work, we used a well-defined photo-triggered non-filamentous discharge to study acetone decomposition in N₂/O₂ gas mixtures. The plasma reactor is associated to a compact chemical ionization FTICR mass spectrometer (BTrap) in order to identify and quantify in real-time acetone and by-products in the plasma. Presence of oxygen (1 to 5%) decreased notably acetone degradation. A tremendous change is observed in the by-products distribution concomitantly to a global decrease of their total concentration. While main products observed in oxygen-free gas mix are nitrile compounds, in oxygenated media they are replaced by formaldehyde, methanol and ketene. Methanol is maximum for 1% of O₂ whereas formaldehyde and ketene concentration reach their maximum value at the highest oxygen concentration tested (5%). A number of nitrate, nitrite and isocyanate organic compounds (C1 and C2) are observed as well with HNO₂, HNO₃ and HNCO.

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.

Filming ultrafast roaming-mediated isomerization of bismuth triiodide in solution

Roaming reaction, defined as a reaction yielding products via reorientational motion in the long-range region (3 - 8 Å) of the potential, is a relatively recently proposed reaction pathway and is now regarded as a universal mechanism that can explain the unimolecular dissociation and isomerization of various molecules. The structural movements of the partially dissociated fragments originating from the frustrated bond fission at the onset of roaming, however, have been explored mostly via theoretical simulations and rarely observed experimentally. Here, we report an investigation of the structural dynamics during a roaming-mediated isomerization reaction of bismuth triiodide (BiI₃) in acetonitrile solution using femtosecond time-resolved x-ray liquidography. Structural analysis of the data visualizes the atomic movements during the roaming-mediated isomerization process including the opening of the Bi-I_b-I_c angle and the closing of I_a-Bi-I_b-I_c dihedral angle, each by ~40°, as well as the shortening of the I_b···I_c distance, following the frustrated bond fission.

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

We explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

Photodissociation Dynamics of CO-Forming Channels on the Ground-State Surface of Methyl Formate at 248 nm: Direct Dynamics Study and Assessment of Generalized Multicenter Impulsive Models

The photodissociation dynamics of methyl formate in the electronic ground state S₀, initiated by a 248 nm-wavelength laser, is studied by direct dynamics simulations. We analyze five channels, where four of them have as products CH₃OH + CO, one leading to the formation of three fragments, H₂CO + H₂ + CO, and a channel characterized by a roaming transition state. The analysis of energy distribution among the degrees of freedom of the product and the comparison with experimental results previously published by other groups provide the ingredients to distinguish the examined dissociation pathways. The interpretation of the results proves that the characterization of dissociation mechanisms must rely on a dynamics approach involving multiple electronic states, including considerations on the features of the S₁/S₀ conical intersection. Here, we also assess the generalized multicenter impulsive model, GMCIM, that has been designed for dissociation processes with exit barriers, and the energy distribution in the products is predicted on the basis of information from the saddle points and the intrinsic reaction coordinates. Main features, advantages, limits, and future perspectives of the method are reported and discussed.

Quasiclassical simulations based on cluster models reveal vibration-facilitated roaming in the isomerization of CO adsorbed on NaCl

The desire to better understand the quantum nature of isomerization led to recent experimental observations of the vibrationally induced isomerization of OC-NaCl(100) to CO-NaCl(100). To investigate the mechanism of this isomerization, we performed dynamics calculations using finite (CO-NaCl)_n cluster models. We constructed new potential energy surfaces for CO-NaCl and CO-CO interactions using high-level ab initio data and report key properties of the bare CO-NaCl potential energy surface, which show much in common with the experiment. We investigated the isomerization dynamics using several cluster models and, in all cases, isomerization was seen for highly excited CO vibrational states, in agreement with experiments. A detailed examination of the reaction trajectories indicates that isomerization occurs when the distance between CO and NaCl is larger than the distance at the conventional isomerization saddle point, which is a strong indicator of 'roaming'.

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.

Roaming Reactions and Dynamics in the van der Waals Region

Roaming reactions were first clearly identified in photodissociation of formaldehyde 15 years ago, and roaming dynamics are now recognized as a universal aspect of chemical reactivity. These reactions typically involve frustrated near-dissociation of a quasibound system to radical fragments, followed by reorientation at long range and intramolecular abstraction. The consequences can be unexpected formation of molecular products, depletion of the radical pool in chemical systems, and formation of products with unusual internal state distributions. In this review, I examine some current aspects of roaming reactions with an emphasis on experimental results, focusing on possible quantum effects in roaming and roaming dynamics in bimolecular systems. These considerations lead to a more inclusive definition of roaming reactions as those for which key dynamics take place at long range.

Roaming Dynamics and Conformational Memory in Photolysis of Formic Acid at 193 nm Using Time-resolved Fourier-transform Infrared Emission Spectroscopy

In photodissociation of trans-formic acid (HCOOH) at 193 nm, we have observed two molecular channels of CO + H₂O and CO₂ + H₂ by using 1 μs-resolved Fourier-transform infrared emission spectroscopy. With the aid of spectral simulation, the CO spectra are rotationally resolved for each vibrational state (v = 1-8). Each of the resulting vibrational and rotational population distributions is characteristic of two Boltzmann profiles with different temperatures, originating from either transition state pathway or OH-roaming to form the same CO + H₂O products. The H₂O roaming co-product is also spectrally simulated to understand the interplay with the CO product in the internal energy partitioning. Accordingly, this work has evaluated the internal energy disposal for the CO and H₂O roaming products; especially the vibrational-state dependence of the roaming signature is reported for the first time. Further, given a 1 μs resolution, the temporal dependence of the CO/CO₂ product ratio at v ≥ 1 rises from 3 to 10 of study, thereby characterizing the effect of conformational memory and well reconciling with the disputed results reported previously between absorption and emission methods.

Imaging the infrared multiphoton excitation and dissociation of propargyl chloride

Infrared multiphoton excitation is combined with UV excitation and state-resolved probes of Cl(²P_3/2), Cl*(²P_1/2), and HCl to study the photochemistry of propargyl chloride. The results show evidence both of infrared multiphoton dissociation on the ground electronic state and infrared multiphoton excitation followed by UV dissociation. The results are interpreted with the aid of a full characterization of the stationary points on the ground state using ab initio methods, as well as our recent experimental and theoretical characterization of the UV photochemistry of the molecule. The data suggest elimination of HCl on the ground electronic state produces linear propadienylidene as a coproduct over a roaming-like transition state that accesses the Cl-H-C abstraction geometry. This identification is supported by separate chirped-pulse microwave studies in a quasi-uniform flow also reported here.

Invited Review Article: Photofragment imaging

Photodissociation studies in molecular beams that employ position-sensitive particle detection to map product recoil velocities emerged thirty years ago and continue to evolve with new laser and detector technologies. These powerful methods allow application of tunable laser detection of single product quantum states, simultaneous measurement of velocity and angular momentum polarization, measurement of joint product state distributions for the detected and undetected products, coincident detection of multiple product channels, and application to radicals and ions as well as closed-shell molecules. These studies have permitted deep investigation of photochemical dynamics for a broad range of systems, revealed new reaction mechanisms, and addressed problems of practical importance in atmospheric, combustion, and interstellar chemistry. This review presents an historical overview, a detailed technical account of the range of methods employed, and selected experimental highlights illustrating the capabilities of the method.

Theories and simulations of roaming

The phenomenon of roaming in chemical reactions has now become both commonly observed in experiment and extensively supported by theory and simulations. Roaming occurs in highly-excited molecules when the trajectories of atomic motion often bypass the minimum energy pathway and produce reaction in unexpected ways from unlikely geometries. The prototypical example is the unimolecular dissociation of formaldehyde (H₂CO), in which the "normal" reaction proceeds through a tight transition state to yield H₂ + CO but for which a high fraction of dissociations take place via a "roaming" mechanism in which one H atom moves far from the HCO, almost to dissociation, and then returns to abstract the second H atom. We review below the theories and simulations that have recently been developed to address and understand this new reaction phenomenon.

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.

Perspective: chemical dynamics simulations of non-statistical reaction dynamics

Non-statistical chemical dynamics are exemplified by disagreements with the transition state (TS), RRKM and phase space theories of chemical kinetics and dynamics. The intrinsic reaction coordinate (IRC) is often used for the former two theories, and non-statistical dynamics arising from non-IRC dynamics are often important. In this perspective, non-statistical dynamics are discussed for chemical reactions, with results primarily obtained from chemical dynamics simulations and to a lesser extent from experiment. The non-statistical dynamical properties discussed are: post-TS dynamics, including potential energy surface bifurcations, product energy partitioning in unimolecular dissociation and avoiding exit-channel potential energy minima; non-RRKM unimolecular decomposition; non-IRC dynamics; direct mechanisms for bimolecular reactions with pre- and/or post-reaction potential energy minima; non-TS theory barrier recrossings; and roaming dynamics.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Electron Paramagnetic Resonance Spectroscopic Study on Nonequilibrium Reaction Pathways in the Photolysis of Solid Nitromethane (CH3NO2) and D3-Nitromethane (CD3NO2)

Thin films of nitromethane (CH3NO2) along with its isotopically labeled counterpart D3-nitromethane (CD3NO2) were photolyzed at discrete wavelength between 266 nm (4.7 eV) and 121 nm (10.2 eV) to explore the underlying mechanisms involved in the decomposition of model compounds of energetic materials in the condensed phase at 5 K. The chemical modifications of the ices were traced in situ via electron paramagnetic resonance, thus focusing on the detection of (hitherto elusive) reaction intermediates and products with unpaired electrons. These studies revealed the formation of two carbon-centered radicals [methyl (CH3), nitromethyl (CH2NO2)], one oxygen-centered radical [methoxy (CH3O)], two nitrogen-centered radicals [nitrogen monoxide (NO), nitrogen dioxide (NO2)], as well as atomic hydrogen (H). The decomposition products of these channels and the carbon-centered nitromethyl (CH2NO2) radical in particular represent crucial reaction intermediates leading via sequential molecular mass growth processes in the exposed nitromethane samples to complex organic molecules as predicted previously by dynamics calculations. The detection of the nitromethyl (CH2NO2) radical along with atomic hydrogen (H) demonstrated the existence of a high-energy decomposition pathway, which is closed under collisionless conditions in the gas phase.

Pruning the Hamiltonian Matrix in MULTIMODE: Test for C2H4 and Application to CH3NO2 Using a New Ab Initio Potential Energy Surface

We report vibrational self-consistent field/virtual state configuration interaction energies of nitromethane using the code MULTIMODE and a new full-dimensional potential energy surface (PES). The PES is a precise, permutationally invariant linear least-squares fit to 17,049 electronic energies, using the CCSD(T)-F12b method with HaDZ basis (cc-pVDZ basis for H atoms, and aug-cc-pVDZ basis for C, O, N atoms). Nitromethane has 15 vibrational degrees of freedom, including one that is a nearly free internal methyl torsion, which is accurately described by the PES. This torsional mode makes vibrational calculations very challenging and here we present the results of calculations without it. Nevertheless, 14-mode calculations are still challenging and can lead to very large Hamiltonian matrices. To address this issue, we apply a pruning scheme, suggested previously by Handy and Carter, that reduces the size of matrix without sacrificing accuracy in the eigenvalues. The method is briefly described here in the context of partitioning theory. A new and more efficient implementation of it, coded in the latest version of MULTIMODE program, is described. The accuracy and efficiency are demonstrated for 12-mode C2H4 and then applied to CH3NO2. Agreement with available experimental values of the CH3NO2 14 fundamentals is very good. Diffusion Monte Carlo calculations in full dimensionality are done for the zero-point energy and wavefuction. These indicate that the torisonal motion is nearly a free-rotor in this state.

Characterization of molecular channel in photodissociation of SOCl2 at 248 nm: Cl2 probing by cavity ring-down absorption spectroscopy

A primary elimination channel of the chlorine molecule in the one-photon dissociation of SOCl2 at 248 nm was investigated using cavity ring-down absorption spectroscopy (CRDS). By means of spectral simulation, the ratio of the vibrational population in the v = 0, 1, and 2 levels was evaluated to be 1 : (0.10 ± 0.02) : (0.009 ± 0.005), corresponding to a Boltzmann vibrational temperature of 340 ± 30 K. The Cl2 molecular channel was obtained with a quantum yield of 0.4 ± 0.2 from the X(1)A' ground state of SOCl2via internal conversion. The dissociation mechanism differs from a prior study where a smaller yield of <3% was obtained, initiated from the 2(1)A' excited state. Temperature-dependence measurements of the Cl2 fragment turn out to support our mechanism. With the aid of ab initio potential energy calculations, two dissociation routes to the molecular products were found, including one synchronous dissociation pathway via a three-center transition state (TS) and the other sequential dissociation pathway via a roaming-mediated isomerization TS. The latter mechanism with a lower energy barrier dominates the dissociation reaction.

Insight into photofragment vector correlation by a multi-center impulsive model

A multi-center impulsive model has been recently developed to characterize the dynamic feature of fragment vector correlation in photodissociation of formaldehyde, H₂CO → CO + H₂, via both transition state and roaming pathways.

Lyman α photolysis of solid nitromethane (CH3NO2) and D3-nitromethane (CD3NO2) – untangling the reaction mechanisms involved in the decomposition of model energetic materials

Solid nitromethane and D3-nitromethane ices were exposed to Lyman α photons to investigate the mechanism involved in the decomposition of energetic materials in the condensed phase.

Dissociation Pathways of Benzylpyridinium "Thermometer" Ions Depend on the Activation Regime: An IRMPD Spectroscopy Study

The dissociation of benzylpyridinium "thermometer" ions is widely used to calibrate the internal energy of ions produced in mass spectrometry. The fragmentation mechanism is usually believed to yield a benzylium cation, although recent studies suggest the possibility of a rearrangement leading to the tropylium isomer, which would compromise the accuracy of energy calibrations. In this study, we used IRMPD spectroscopy to probe the dissociation pathways of the p-(tert-butyl)benzylpyridinium ion. Our results show that the formation of both benzylium and tropylium products is feasible depending on the activation regime and on the reaction time scale. Varying the trapping delays in the hexapole gives insight into a rearrangement mechanism occurring through consecutive reactions with an isomerization from benzylium to tropylium. Our work provides experimental validations for the established calibration procedure and highlights the adequacy of IRMPD spectroscopy to qualitatively resolve gas-phase rearrangement kinetics.

A Signature of Roaming Dynamics in the Thermal Decomposition of Ethyl Nitrite: Chirped-Pulse Rotational Spectroscopy and Kinetic Modeling

Chirped-pulse (CP) Fourier transform rotational spectroscopy is uniquely suited for near-universal quantitative detection and structural characterization of mixtures that contain multiple molecular and radical species. In this work, we employ CP spectroscopy to measure product branching and extract information about the reaction mechanism, guided by kinetic modeling. Pyrolysis of ethyl nitrite, CH3CH2ONO, is studied in a Chen type flash pyrolysis reactor at temperatures of 1000-1800 K. The branching between HNO, CH2O, and CH3CHO products is measured and compared to the kinetic models generated by the Reaction Mechanism Generator software. We find that roaming CH3CH2ONO → CH3CHO + HNO plays an important role in the thermal decomposition of ethyl nitrite, with its rate, at 1000 K, comparable to that of the radical elimination channel CH3CH2ONO → CH3CH2O + NO. HNO is a signature of roaming in this system.