The Roaming Atom: Straying from the Reaction Path in Formaldehyde Decomposition

Molecular halogen elimination from halogen-containing compounds in the atmosphere

Atmospheric halogen chemistry has drawn much attention, because the halogen atom (X) playing a catalytic role may cause severe stratospheric ozone depletion. Atomic X elimination from X-containing hydrocarbons is recognized as the major primary dissociation process upon UV-light irradiation, whereas direct elimination of the X2 product has been seldom discussed or remained a controversial issue. This account is intended to review the detection of X2 primary products using cavity ring-down absorption spectroscopy in the photolysis at 248 nm of a variety of X-containing compounds, focusing on bromomethanes (CH2Br2, CF2Br2, CHBr2Cl, and CHBr3), dibromoethanes (1,1-C2H4Br2 and 1,2-C2H4Br2) and dibromoethylenes (1,1-C2H2Br2 and 1,2-C2H2Br2), diiodomethane (CH2I2), thionyl chloride (SOCl2), and sulfuryl chloride (SO2Cl2), along with a brief discussion on acyl bromides (BrCOCOBr and CH2BrCOBr). The optical spectra, quantum yields, and vibrational population distributions of the X2 fragments have been characterized, especially for Br2 and I2. With the aid of ab initio calculations of potential energies and rate constants, the detailed photodissociation mechanisms may be comprehended. Such studies are fundamentally important to gain insight into the dissociation dynamics and may also practically help to assess the halogen-related environmental variation.

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.

Identification of atomic-level mechanisms for gas-phase X- + CH3Y SN2 reactions by combined experiments and simulations

For the traditional model of gas-phase X(-) + CH3Y SN2 reactions, C3v ion-dipole pre- and postreaction complexes X(-)---CH3Y and XCH3---Y(-), separated by a central barrier, are formed. Statistical intramolecular dynamics are assumed for these complexes, so that their unimolecular rate constants are given by RRKM theory. Both previous simulations and experiments have shown that the dynamics of these complexes are not statistical and of interest is how these nonstatistical dynamics affect the SN2 rate constant. This work also found there was a transition from an indirect, nonstatistical, complex forming mechanism, to a direct mechanism, as either the vibrational and/or relative translational energy of the reactants was increased. The current Account reviews recent collaborative studies involving molecular beam ion-imaging experiments and direct (on-the-fly) dynamics simulations of the SN2 reactions for which Cl(-), F(-), and OH(-) react with CH3I. Also considered are reactions of the microsolvated anions OH(-)(H2O) and OH(-)(H2O)2 with CH3I. These studies have provided a detailed understanding of the atomistic mechanisms for these SN2 reactions. Overall, the atomistic dynamics for the Cl(-) + CH3I SN2 reaction follows those found in previous studies. The reaction is indirect, complex forming at low reactant collision energies, and then there is a transition to direct reaction between 0.2 and 0.4 eV. The direct reaction may occur by rebound mechanism, in which the ClCH3 product rebounds backward from the I(-) product or a stripping mechanism in which Cl(-) strips CH3 from the I atom and scatters in the forward direction. A similar indirect to direct mechanistic transition was observed in previous work for the Cl(-) + CH3Cl and Cl(-) + CH3Br SN2 reactions. At the high collision energy of 1.9 eV, a new indirect mechanism, called the roundabout, was discovered. For the F(-) + CH3I reaction, there is not a transition from indirect to direct reaction as Erel is increased. The indirect mechanism, with prereaction complex formation, is important at all the Erel investigated, contributing up ∼60% of the reaction. The remaining direct reaction occurs by the rebound and stripping mechanisms. Though the potential energy curve for the OH(-) + CH3I reaction is similar to that for F(-) + CH3I, the two reactions have different dynamics. They are akin, in that for both there is not a transition from an indirect to direct reaction. However, for F(-) + CH3I indirect reaction dominates at all Erel, but it is less important for OH(-) + CH3I and becomes negligible as Erel is increased. Stripping is a minor channel for F(-) + CH3I, but accounts for more than 60% of the OH(-) + CH3I reaction at high Erel. Adding one or two H2O molecules to OH(-) alters the reaction dynamics from that for unsolvated OH(-). Adding one H2O molecule enhances indirect reaction at low Erel, and changes the reaction mechanism from primarily stripping to rebound at high Erel. With two H2O molecules the dynamics is indirect and isotropic at all collision energies.

Complete active space second order perturbation theory (CASPT2) study of N(²D) + H₂O reaction paths on D₁ and D₀ potential energy surfaces: direct and roaming pathways

We report reaction paths starting from N((2)D) + H2O for doublet spin states, D0 and D1. The potential energy surfaces are explored in an automated fashion using the global reaction route mapping strategy. The critical points and reaction paths have been fully optimized at the complete active space second order perturbation theory level taking all valence electrons in the active space. In addition to direct dissociation pathways that would be dominant, three roaming processes, two roaming dissociation, and one roaming isomerization: (1) H2ON → H-O(H)N → H-HON → NO((2)Π) + H2, (2) cis-HNOH → HNO-H → H-HNO → NO + H2, (3) H2NO → H-HNO → HNO-H → trans-HNOH, are confirmed on the D0 surface.

Photochemistry. Evidence for direct molecular oxygen production in CO₂ photodissociation

Photodissociation of carbon dioxide (CO2) has long been assumed to proceed exclusively to carbon monoxide (CO) and oxygen atom (O) primary products. However, recent theoretical calculations suggested that an exit channel to produce C + O2 should also be energetically accessible. Here we report the direct experimental evidence for the C + O2 channel in CO2 photodissociation near the energetic threshold of the C((3)P) + O2(X(3)Σ(g)(-)) channel with a yield of 5 ± 2% using vacuum ultraviolet laser pump-probe spectroscopy and velocity-map imaging detection of the C((3)PJ) product between 101.5 and 107.2 nanometers. Our results may have implications for nonbiological oxygen production in CO2-heavy atmospheres.

Dynamics of chlorine atom reactions with hydrocarbons: insights from imaging the radical product in crossed beams

We present a comprehensive overview of our ongoing studies applying dc slice imaging in crossed molecular beams to probe the dynamics of chlorine atom reactions with polyatomic hydrocarbons. Our approach consists in measuring the full velocity-flux contour maps of the radical products using vacuum ultraviolet "soft" photoionization at 157 nm. Our overall goal is to extend the range of chemical dynamics investigations from simple triatomic or tetraatomic molecules to systematic investigations of a sequence of isomers or a homologous series of reactants of intermediate size. These experimental investigations are augmented by high-level ab initio calculations which, taken together, reveal trends in product energy and angular momentum partitioning and offer deep insight into the reaction mechanisms as a function of structure, bonding patterns, and kinematics. We explore these issues in alkanes, for which only direct reactive encounters are found, and in unsaturated hydrocarbons, for which an addition-elimination mechanism competes with direct abstraction. The results for alkene addition-elimination in particular suggest a new view of these reactions: The only pathway to HCl elimination is accessed by means of roaming excursions of the Cl atom from the strongly bound adduct.

Uncovering the role of intra- and intermolecular motion in frustrated Lewis acid/base chemistry: ab initio molecular dynamics study of CO2 binding by phosphorus/boron frustrated Lewis pair [tBu3P/B(C6F5)3]

The role of the intra- and intermolecular motion, i.e., molecular vibrations and the relative motion of reactants, remains largely unexplored in the frustrated Lewis acid/base chemistry. Here, we address the issue with the ab initio molecular dynamics (AIMD) study of CO2 binding by a Lewis acid (LA) and a Lewis base (LB), i.e., tBu3P + CO2 + B(C6F5)3 → tBu3P-C(O)O-B(C6F5)3 ([1]). Reasonably large ensemble of AIMD trajectories propagated at 300 K from structures in the saddle region as well as trajectories propagated directly from the reactants region revealed an effect arising from significant recrossing of the saddle area. The effect is that transient complexes composed of weakly interacting reactants nearly cease to progress along the segment of the minimum energy pathway (MEP) at the saddle region for a (subpicosecond) period of time during which the dominant factor is the light-to-heavy type of relative motion of the vibrating reactants, i.e., the "bouncing"-like movement of CO2 with respect to much heavier phosphine and borane as main contributor to the mode that is perpendicular to the MEP-direction. In terms of how P···C and B···O distances change with time, the roaming-like patterns of typical AIMD trajectories, reactive and nonreactive alike, extend far beyond the saddle region. In addition to the dynamical portrayal of [1], we provide the energy-landscape perspective that takes into account the hierarchy of time scales. The verifiable implication of the effect found here is that the isotopically substituted (heavier) LB/LA "pair" should be less reactive that the "normal" and thus lighter counterpart.

Persistence of transition-state structure in chemical reactions driven by fields oscillating in time

Chemical reactions subjected to time-varying external forces cannot generally be described through a fixed bottleneck near the transition-state barrier or dividing surface. A naive dividing surface attached to the instantaneous, but moving, barrier top also fails to be recrossing-free. We construct a moving dividing surface in phase space over a transition-state trajectory. This surface is recrossing-free for both Hamiltonian and dissipative dynamics. This is confirmed even for strongly anharmonic barriers using simulation. The power of transition-state theory is thereby applicable to chemical reactions and other activated processes even when the bottlenecks are time dependent and move across space.

Mode specificity and product energy disposal in unimolecular reactions: insights from the sudden vector projection model

A simple model is proposed to predict mode specificity and product energy disposal in unimolecular dissociation reactions. This so-called Sudden Vector Projection (SVP) model quantifies the coupling of a reactant or product mode with the reaction coordinate at the transition state by projecting the corresponding normal mode vector onto the imaginary frequency mode at the saddle point. Due to the sudden assumption, SVP predictions for mode specificity are expected to be valid only when the reactant molecule has weak intermodal coupling. On the other hand, the sudden limit is generally satisfied for its predictions of product energy disposal in unimolecular reactions with a tight barrier. The SVP model is applied to several prototypical systems and the agreement with available experimental and theoretical results is satisfactory.

Biosynthetic consequences of multiple sequential post-transition-state bifurcations

Selectivity in chemical reactions that form complex molecular architectures from simpler precursors is usually rationalized by comparing competing transition-state structures that lead to different possible products. Herein we describe a system for which a single transition-state structure leads to the formation of many isomeric products via pathways that feature multiple sequential bifurcations. The reaction network described connects the pimar-15-en-8-yl cation to miltiradiene, a tricyclic diterpene natural product, and isomers via cyclizations and/or rearrangements. The results suggest that the selectivity of the reaction is controlled by (post-transition-state) dynamic effects, that is, how the carbocation structure changes in response to the distribution of energy in its vibrational modes. The inherent dynamical effects revealed herein (characterized through quasiclassical direct dynamics calculations using density functional theory) have implications not only for the general principles of selectivity prediction in systems with complex potential energy surfaces, but also for the mechanisms of terpene synthase enzymes and their evolution. These findings redefine the challenges faced by nature in controlling the biosynthesis of complex natural products.

Photodissociation of Propionaldehyde at 248 nm: Roaming Pathway as an Increasingly Important Role in Large Aliphatic Aldehydes

Time-resolved Fourier transform infrared emission spectroscopy is employed in the photolysis of propionaldehyde (CH3CH2CHO) at 248 nm to characterize the role of the roaming pathway. High-resolution spectra of CO are analyzed to yield a single Boltzmann rotational distribution for each vibrational level (ν = 1-4) with small rotational and large vibrational energy disposals. A roaming saddle point is found containing two far separated moieties of HCO and CH3CH2 with a weak interaction between them. Quasiclassical trajectory calculations on this configuration yield the CO energy flow behavior, consistent with the findings. The rate constant along the roaming pathway is evaluated to be larger by >1-2 orders of magnitude than those along tight transition state or three-body dissociation pathways. This work implies that the roaming mechanism plays an increasingly important role in aliphatic aldehydes as the molecular size becomes larger.

Two roaming pathways in the photolysis of CH3CHO between 328 and 308 nm

We attribute the two product-state distributions previously seen in CH₃CHO photodissociation to CH₃-roaming and H-roaming, unifying all previous experimental results.

Roaming Dissociation of Ethyl Radicals

Previous studies on the photodissociation of C2H5 reported rate constants for H-atom formation several orders of magnitude smaller than that predicted by Rice-Ramsperger-Kassel-Marcus (RRKM) theory. This Letter provides a potential explanation for this anomaly, based on direct trajectory calculations of C2H5 dissociation. The trajectories reveal the existence of a roaming dissociation channel that leads to the formation of C2H3 and H2. This channel is found to proceed over the ridge between the transition state of H-atom elimination and that of bimolecular H-abstraction. The formed C2H3 radical can subsequently dissociate to C2H2 and a H atom; this secondary dissociation is suggested to be a potential reason for the unexpectedly slow H-atom formation observed in the photodissociation experiments.