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

Dynamics of a gas-phase SNAr reaction: non-concerted mechanism despite the Meisenheimer complex being a transition state

The commonly accepted mechanism of the nucleophilic aromatic substitution (SNAr) reaction has been found to be governed by the nature of the Meisenheimer structure on the potential energy surface. A stable Meisenheimer intermediate favors a stepwise mechanism, while a Meisenheimer transition state favors a concerted mechanism. Here, we show by using a detailed potential energy map (using the DFT and DLPNO-CCSD(T)/CBS methods) and ab initio classical trajectory simulations that the F- + C6H5NO2 SNAr reaction involves a Meisenheimer transition state and follows a stepwise mechanism in contrast to the expected concerted pathway. The stepwise mechanism observed in the trajectory simulations takes place by the formation of various ion-dipole and σ-complexes. While the majority of the trajectories follow the multi-step mechanism and avoid the minimum energy path, a considerable fraction exhibit a roaming atom mechanism where the F atom hovers around the phenyl ring before the formation of the products.

Advances and New Challenges to Bimolecular Reaction Dynamics Theory

Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.

Degradation Mechanism of Benzo[a]pyrene Initiated by the OH Radical and 1O2: An Insight from Density Functional Theory Calculations

The degradation mechanism of benzo[a]pyrene (BaP) initiated by ^•OH and ¹O₂ in aqueous solution is investigated by density functional theory calculations. The main degradation products are BaP-1,6-quinone, BaP-3,6-quinone, BaP-4,6-quinone, and BaP-6,12-quinone. ^•OH and HO₂ are the main intermediate radical species. At a low initial concentration of ^•OH, ¹O₂ could be a primary driver for BaP degradation. The degradation mechanism includes six consecutive elementary reactions: (1) ¹O₂ initiation forming BaP-6-OO. (2) 1,3 H-shift (H atom shifts to the OO group) that is promoted by H₂O, forming BaP-6-OOH. (3) BaP-6-OOH decomposes into the ^•OH radical and BaP-6-O. (4) ^•OH addition to BaP-6-O forming BaP-6-O-1(3,4,12)-OH. (5) Extracting the H atom from the carbon with the OH group by ¹O₂. (6) Extracting the H atom from the OH group by HO₂. At a high initial concentration of ^•OH, the ^•OH-initiated and ¹O₂-initiated degradation reactions of BaP are both feasible. The degradation mechanism includes six consecutive elementary reactions: (1) ^•OH initiation forming BaP-6-OH or ¹O₂ initiation forming BaP-6-OO. (2) ¹O₂ addition to BaP-6-OH forming BaP-6-OH-12(1,3,4)-OO or ^•OH addition to BaP-6-OO forming BaP-6-OO-12(1,3,4)-OH. (3) Extracting the H atom from the carbon with the OH group by ¹O₂, forming HO₂. (4) 1,3 H-shift (H-shift from the carbon to the OO group), promoted by H₂O. (5) The loss of the OH radical. (6) Abstracting the H atom from the OH group by HO₂. In this paper, the formation of BaP-4,6-quinone via the BaP degradation is first reported. Water participates in the elementary reaction in which the H atom attached on the aromatic ring shifts to the OO group, serving as a bridge that stabilizes the transition state and transports the proton. A comprehensive investigation explains the degradation mechanism of BaP initiated by ^•OH and ¹O₂ in aqueous solution.

State-of-the-Art Quantum Chemistry Meets Variable Reaction Coordinate Transition State Theory to Solve the Puzzling Case of the H2S + Cl System

The atmospheric reaction of H₂S with Cl has been reinvestigated to check if, as previously suggested, only explicit dynamical computations can lead to an accurate evaluation of the reaction rate because of strong recrossing effects and the breakdown of the variational extension of transition state theory. For this reason, the corresponding potential energy surface has been thoroughly investigated, thus leading to an accurate characterization of all stationary points, whose energetics has been computed at the state of the art. To this end, coupled-cluster theory including up to quadruple excitations has been employed, together with the extrapolation to the complete basis set limit and also incorporating core-valence correlation, spin-orbit, and scalar relativistic effects as well as diagonal Born-Oppenheimer corrections. This highly accurate composite scheme has also been paralleled by less expensive yet promising computational approaches. Moving to kinetics, variational transition state theory and its variable reaction coordinate extension for barrierless steps have been exploited, thus obtaining a reaction rate constant (8.16 × 10^-11 cm³ molecule^-1 s^-1 at 300 K and 1 atm) in remarkable agreement with the experimental counterpart. Therefore, contrary to previous claims, there is no need to invoke any failure of the transition state theory, provided that sufficiently accurate quantum-chemical computations are performed. The investigation of the puzzling case of the H₂S + Cl system allowed us to present a robust approach for disclosing the thermochemistry and kinetics of reactions of atmospheric and astrophysical interest.

Gas-Phase Synthesis of 3-Vinylcyclopropene via the Crossed Beam Reaction of the Methylidyne Radical (CH; X2 Π) with 1,3-Butadiene (CH2 CHCHCH2 ; X1 Ag )

The crossed molecular beam reactions of the methylidyne radical (CH; X² Π) with 1,3-butadiene (CH₂ CHCHCH₂ ; X¹ A_g ) along with their (partially) deuterated counterparts were performed at collision energies of 20.8 kJ mol^-1 under single collision conditions. Combining our laboratory data with ab initio calculations, we reveal that the methylidyne radical may add barrierlessly to the terminal carbon atom and/or carbon-carbon double bond of 1,3-butadiene, leading to doublet C₅ H₇ intermediates with life times longer than the rotation periods. These collision complexes undergo non-statistical unimolecular decomposition through hydrogen atom emission yielding the cyclic cis- and trans-3-vinyl-cyclopropene products with reaction exoergicities of 119±42 kJ mol^-1 . Since this reaction is barrierless, exoergic, and all transition states are located below the energy of the separated reactants, these cyclic C₅ H₆ products are predicted to be accessed even in low-temperature environments, such as in hydrocarbon-rich atmospheres of planets and cold molecular clouds such as TMC-1.

Three-dimensional covariance-map imaging of molecular structure and dynamics on the ultrafast timescale

Ultrafast laser pump-probe methods allow chemical reactions to be followed in real time, and have provided unprecedented insight into fundamental aspects of chemical reactivity. While evolution of the electronic structure of the system under study is evident from changes in the observed spectral signatures, information on rearrangement of the nuclear framework is generally obtained indirectly. Disentangling contributions to the signal arising from competing photochemical pathways can also be challenging. Here we introduce the new technique of three-dimensional covariance-map Coulomb explosion imaging, which has the potential to provide complete three-dimensional information on molecular structure and dynamics as they evolve in real time during a gas-phase chemical reaction. We present first proof-of-concept data from recent measurements on CF₃I. Our approach allows the contributions from competing fragmentation pathways to be isolated and characterised unambiguously, and is a promising route to enabling the recording of 'molecular movies' for a wide variety of gas-phase chemical processes.

Identifying reaction pathways in phase space via asymptotic trajectories

In this paper, we revisit the concepts of the reactivity map and the reactivity bands as an alternative to the use of perturbation theory for the determination of the phase space geometry of chemical reactions. We introduce a reformulated metric, called the asymptotic trajectory indicator, and an efficient algorithm to obtain reactivity boundaries. We demonstrate that this method has sufficient accuracy to reproduce phase space structures such as turnstiles for a 1D model of the isomerization of ketene in an external field. The asymptotic trajectory indicator can be applied to higher dimensional systems coupled to Langevin baths as we demonstrate for a 3D model of the isomerization of ketene.

Time-resolving the ultrafast H2 roaming chemistry and H3+ formation using extreme-ultraviolet pulses

The time scales and formation mechanisms of tri-hydrogen cation products in organic molecule ionization processes are poorly understood, despite their cardinal role in the chemistry of the interstellar medium and in other chemical systems. Using an ultrafast extreme-ultraviolet pump and time-resolved near-IR probe, combined with high-level ab initio molecular dynamics calculations, here we report unambiguously that H₃⁺ formation in double-ionization of methanol occurs on a sub 100 fs time scale, settling previous conflicting findings of strong-field Coulomb explosion experiments. Our combined experimental-computational studies suggest that ultrafast competition, between proton-transfer and long-range electron-transfer processes, determines whether the roaming neutral H₂ dynamics on the dication result in [Formula: see text] or [Formula: see text] fragments respectively.

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.

Nonstatistical Reaction Dynamics

Nonstatistical dynamics is important for many chemical reactions. The Rice-Ramsperger-Kassel-Marcus (RRKM) theory of unimolecular kinetics assumes a reactant molecule maintains a statistical microcanonical ensemble of vibrational states during its dissociation so that its unimolecular dynamics are time independent. Such dynamics results when the reactant's atomic motion is chaotic or irregular. Intrinsic non-RRKM dynamics occurs when part of the reactant's phase space consists of quasiperiodic/regular motion and a bottleneck exists, so that the unimolecular rate constant is time dependent. Nonrandom excitation of a molecule may result in short-time apparent non-RRKM dynamics. For rotational activation, the 2J + 1 K levels for a particular J may be highly mixed, making K an active degree of freedom, or K may be a good quantum number and an adiabatic degree of freedom. Nonstatistical dynamics is often important for bimolecular reactions and their intermediates and for product-energy partitioning of bimolecular and unimolecular reactions. Post–transition state dynamics is often highly complex and nonstatistical.

Ultrafast Proton Transfer Dynamics on the Repulsive Potential of the Ethanol Dication: Roaming-Mediated Isomerization versus Coulomb Explosion

If a molecular dication is produced on a repulsive potential energy surface (PES), it normally dissociates. Before that, however, ultrafast nuclear dynamics can change the PES and significantly influence the fragmentation pathway. Here, we investigate the electron-impact-induced double ionization and subsequent fragmentation processes of the ethanol molecule using multiparticle coincident momentum spectroscopy and ab initio dynamical simulations. For the electronic ground state of the ethanol dication, we observe several fragmentation channels that cannot be reached by direct Coulomb explosion (CE) but require preceding isomerization. Our simulations show that ultrafast hydrogen or proton transfer (PT) can stabilize the repulsive PES of the dication before the direct CE and form intermediate H₂ or H₂O. These neutrals stay in the vicinity of the precursor, and roaming mechanisms lead to isomerization and finally PT resulting in emission of H₃⁺ or H₃O⁺. The present findings can help to understand the complex fragmentation dynamics of molecular cations.

Long-Range N-N Bonding by Rydberg Electrons

By using high-level ab initio methods, we examine the nature of bonding between Rydberg electrons hosted by two four-coordinate nitrogen centers embedded in a hydrocarbon scaffold. The electronic structure of these species resembles that of diradicals, yet the diffuse nature of the orbitals hosting the unpaired electrons results in unusual features. The unpaired Rydberg electrons exhibit long-range bonding interactions, leading to stabilization of the singlet state (relative to the triplet) and a reduced number of effectively unpaired electrons. However, thermochemical gains due to through-space bonding are offset by strong Coulomb repulsion between positively charged nitrogen cores. The kinetic stability of these Rydberg diradicals may be controlled by a judicious choice of the molecular scaffold, suggesting possible strategies for their experimental characterization.

Reaction Profiles and Kinetics for Radical-Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H₂+ Be, H + NH₂ → H₂ + NH, CH₃ + C₂H₅ → CH₄ + C₂H₄, H + C₂H₅ → H₂ + C₂H₄, and H + HCO → H₂ + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH₂, and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol^-1 and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH₃ + C₂H₅, H + C₂H₅, and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q, 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

Accurate Neural Network Representation of the Ab Initio Determined Spin-Orbit Interaction in the Diabatic Representation Including the Effects of Conical Intersections

A method for fitting ab initio determined spin-orbit coupling interactions, in the Breit-Pauli approximation, based on quasidiabatic representations using neural network fits is reported. The algorithm generalizes our recently reported neural network approach for representing the dipole interaction. The S₀, S₁, and T₁ states of formaldehyde are used as an example. First, the two singlet states S₀ and S₁ are diabatized with a modified Boys Localization diabatization method. Second, the spin-orbit coupling between singlet and triplet states is transformed to the diabatic representation. This removes the discontinuities in the adiabatic representation. The diabatized spin-orbit couplings are then fit with smooth neural network functions. The analytic representation of spin-orbit coupling interactions in a diabatic basis by neural networks will make accurate full-dimensional quantum dynamical treatment of both internal conversion and intersystem crossing possible, which will help us to gain better understanding of both processes.

Real-time tracking of the entangled pathways in the multichannel photodissociation of acetaldehyde

The roaming mechanism, an unconventional reaction path, was discovered more than a decade ago in the studies of formaldehyde photodissociation, H₂CO → H₂ + CO. Since then, observations of roaming have been claimed in numerous photochemical processes. A closer examination of the presented data, however, revealed that evidence for roaming is not always unequivocal, and some of the conclusions could be misleading. We report here an in-depth, joint experimental and theoretical study of the title reaction. By tracking the time-evolution of the pair-correlated product state distributions, we decipher the competing, interwoven reaction pathways that lead to the radical (CH₃ + HCO) and molecular (CH₄ + CO) products. Possible roaming pathways are then elucidated and a more precise descriptor of the phenomenon is delineated.

Strong-field control of H3 + production from methanol dications: Selecting between local and extended formation mechanisms

Using the CD₃OH isotopologue of methanol, the ratio of D₂H⁺ to D₃ ⁺ formation is manipulated by changing the characteristics of the intense femtosecond laser pulse. Detection of D₂H⁺ indicates a formation process involving two hydrogen atoms from the methyl side of the molecule and a proton from the hydroxyl side, while detection of D₃ ⁺ indicates local formation involving only the methyl group. Both mechanisms are thought to involve a neutral D₂ moiety. An adaptive control strategy that employs image-based feedback to guide the learning algorithm results in an enhancement of the D₂H⁺/D₃ ⁺ ratio by a factor of approximately two. The optimized pulses have secondary structures 110-210 fs after the main pulse and result in photofragments that have different kinetic energy release distributions than those produced from near transform limited pulses. Systematic changes to the linear chirp and higher order dispersion terms of the laser pulse are compared to the results obtained with the optimized pulse shapes.

Strong inverse kinetic isotope effect observed in ammonia charge exchange reactions

Isotopic substitution has long been used to understand the detailed mechanisms of chemical reactions; normally the substitution of hydrogen by deuterium leads to a slower reaction. Here, we report our findings on the charge transfer collisions of cold [Formula: see text] ions and two isotopologues of ammonia, [Formula: see text] and [Formula: see text]. Deuterated ammonia is found to react more than three times faster than hydrogenated ammonia. Classical capture models are unable to account for this pronounced inverse kinetic isotope effect. Moreover, detailed ab initio calculations cannot identify any (energetically accessible) crossing points between the reactant and product potential energy surfaces, indicating that electron transfer is likely to be slow. The higher reactivity of [Formula: see text] is attributed to the greater density of states (and therefore lifetime) of the deuterated reaction complex compared to the hydrogenated system. Our observations could provide valuable insight into possible mechanisms contributing to deuterium fractionation in the interstellar medium.

Collision-induced and complex-mediated roaming dynamics in the H + C2H4 → H2 + C2H3 reaction

Roaming is a novel mechanism in reaction dynamics. It describes an unusual pathway, which can be quite different from the conventional minimum-energy path, leading to products. While roaming has been reported or suggested in a number of unimolecular reactions, it has been rarely reported for bimolecular reactions. Here, we report a high-level computational study of roaming dynamics in the important bimolecular combustion reaction H + C₂H₄ → H₂ + C₂H₃, using a new, high-level machine learning-based potential energy surface. In addition to the complex-mediated roaming mechanism, a non-complex forming roaming mechanism is found. It can be described as a direct inelastic collision where the departing H atom roams and then abstracts an H atom. We denoted this as “collision-induced” roaming. These two roaming mechanisms have different angular distributions; however, both produce highly internally excited C₂H₃. The roaming pathway leads to remarkably different dynamics as compared with the direct abstraction pathway. A clear signature of the roaming mechanism is highly internally excited C₂H₃, which could be observed experimentally.

Collision-induced and complex-mediated roaming mechanisms are revealed for an important bimolecular reaction in combustion.

Signature of a conical intersection in the dissociative photoionization of formaldehyde

Electron/ion coincidence experiments and ab initio calculations of the dissociative photoionization of formaldehyde reveal the presence of a conical intersection controlling the dynamics and favoring dissociation into the molecular channel, CO⁺ + H₂.

Universal crossed beam imaging studies of polyatomic reaction dynamics

Crossed-beam imaging studies of polyatomic reactions show surprising dynamics not anticipated by extrapolation from smaller model systems.

Long-range versus short-range effects in cold molecular ion-neutral collisions

The investigation of cold interactions between ions and neutrals has recently emerged as a new scientific frontier at the interface of physics and chemistry. Here, we report a study of charge-transfer (CT) collisions of Rb atoms with N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}^{+}$$\end{document}2+ and O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}^{+}$$\end{document}2+ ions in the mK regime using a dynamic ion-neutral hybrid trapping experiment. We observe markedly different CT kinetics and dynamics for the different systems and reaction channels studied. While the kinetics in some channels are consistent with classical capture theory, others show distinct non-universal dynamics. The experimental results are interpreted with the help of classical-capture, quasiclassical-trajectory and quantum-scattering calculations using ab-initio potentials for the highly excited molecular states involved. The theoretical analysis reveals an intricate interplay between short- and long-range effects in the different reaction channels which ultimately determines the CT dynamics and rates. Our results illustrate salient mechanisms that determine the efficiency of cold molecular CT reactions.

Studies on reactions between cold molecular ions and neutral atoms provide insights into intermolecular interactions. Here the authors explore the kinetics and dynamics of charge-transfer collisions between the cold N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}^{+}$$\end{document}2+ and O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}^{+}$$\end{document}2+ ions and neutral Rb atoms and discuss the role of long- and short-range effects.

Self-supported hydrogenolysis of aromatic ethers to arenes

Arenes are widely used chemicals and essential components in liquid fuels, which are currently produced from fossil feedstocks. Here, we proposed the self-supported hydrogenolysis (SSH) of aromatic ethers to produce arenes using the hydrogen source within the reactants, and it was found that RuW alloy nanoparticles were very efficient catalyst for the reactions. This route is very attractive and distinguished from the reported studies on the cleavage of the C_Ar─O bonds. The unique feature of this methodology is that exogenous hydrogen or other reductant is not required, and hydrogenation of aromatic rings could be avoided completely. The selectivities to arenes could reach >99.9% at complete conversion of the ethers. Moreover, lignin could also be transformed into arenes efficiently over the RuW alloy catalyst. The mechanism studies showed that the neighboring Ru and W species in the RuW alloy nanoparticles worked synergistically to accelerate the SSH reaction.

Quasi-Classical Trajectory Study of NH(3∑-) + NH(3∑-) Reactive Collisions

A full-dimension quasi-classical trajectory study of collisions between two NH radicals is presented. Interatomic interactions are represented by a previously reported global six-dimensional potential energy surface for singlet electronic state of the N₂H₂ system. This study suggests that the formation of N₂ from the collision of two NH radicals may occur via a one-step (NH + NH → N₂ + H + H) or two-step (NH + NH → N₂H + H → N₂ + H + H) microscopic reaction mechanism. A fast vibrational energy redistribution is observed in the four-body complex in the latter mechanism. Excitation functions are presented and discussed. A variant of the vibrational energy quantum mechanical threshold method was used to correct the zero-point energy leakage in the classical calculations. The influence of reactant's rotational and vibrational energy on reactivity was investigated by state-specific calculations with one or both reactants excited. Reaction rate constants for the ground and some rotationally excited states are presented using an Arrhenius-Kooij-like functional form.

Observation of the Carbon Elimination Channel in Vacuum Ultraviolet Photodissociation of OCS

The textbook mechanism for OCS photodissociation mainly involves the CO + S or CS + O product channel via a single bond fission. However, a third dissociation channel concerning the cleavage of both C-S and C-O bonds yielding SO + C products, though thermodynamically allowed, has never been verified experimentally to date. By using a tunable vacuum ultraviolet laser light and time-sliced velocity map ion imaging technique, we have clearly observed the SO(X³Σ^-) + C(³P_J=0) products as the vacuum ultraviolet laser photon energy gradually exceeds its thermodynamic threshold. The corresponding SO(X³Σ^-) coproducts are highly vibrationally excited and show varying angular distributions from isotropic to anisotropic as the excitation photon energy increases. Theoretical analysis suggests that a fast nonadiabatic pathway plays a dominant role in the formation of the anisotropic SO products. That isotropic products arise as the excitation photon energies approach the thermodynamic threshold can be reasonably explained by the "roaming mechanism".

The Role of Ozone Vibrational Resonances in the Isotope Exchange Reaction 16O16O + 18O → 18O16O + 16O: The Time-Dependent Picture

We consider the time-dependent dynamics of the isotope exchange reaction in collisions between an oxygen molecule and an oxygen atom: ¹⁶O¹⁶O + ¹⁸O → ¹⁶O¹⁸O + ¹⁶O. A theoretical approach using the multiconfiguration time-dependent Hartree method was employed to model the time evolution of the reaction. Two potential surfaces available in the literature were used in the calculations, and the results obtained with the two surfaces are compared with each other as well as with results of a previous theoretical time-independent approach. A good agreement for the reaction probabilities with the previous theoretical results is found. Comparing the results obtained using two potential energy surfaces allows us to understand the role of the reef/shoulder-like feature in the minimum energy path of the reaction in the isotope exchange process. Also, it was found that the distribution of final products of the reaction is highly anisotropic, which agrees with experimental observations and, at the same time, suggests that the family of approximated statistical approaches, assuming a randomized distribution over final exit channels, is not applicable to this case.

Chemical Kinetics Approves the Occurrence of C (3Pj) Reaction with H2O

Although both atomic carbon and water are omnipresent in human life, there is a debate about the possibility of carbon reaction with water. Some low-temperature spectroscopic investigations have rejected the reaction, whereas some room-temperature experiments and theoretical studies have accepted the possibility of the reaction by reporting rate coefficients ranging from 10⁵ to 10⁹ L mol^-1 s^-1. This study provides new lines of evidence about the reaction through exploration of the reaction mechanism using the CCSD(T) method and solving the corresponding master equation by following two main approaches. According to the results, the rate coefficient of the reaction is significantly influenced by the tunneling and hindered rotation effects, in addition to the selected total angular momentum (J). Furthermore, the total rate coefficient of the reaction increases dramatically (from 10⁷ to 10¹¹ L mol^-1 s^-1) with the rise of temperature from 100 to 4000 K, while the total rate coefficient is insensitive to pressure (0.1-10 atm). Despite some differences between the results of the two approaches, the rate coefficients of both methods are consistent with the previously reported rate coefficients. Also, in agreement with the previous studies, the major products are ²HOC + ²H and ²HCO + ²H. In general, the findings approve the occurrence of the title reaction and indicate that the mentioned conflict is due to the sensitivity of the reaction to the investigated temperature and J level. The sensitivity does not permit low-temperature spectroscopic studies to detect any products and varies the measured and calculated rate coefficients.

Can reactions follow non-traditional second-order saddle pathways avoiding transition states?

We report here an ab initio (CASSCF/6-31+G*) trajectory simulation study on the mechanisms of the denitrogenation of 1-pyrazoline and its subsituted analogue that reveals reaction pathways via a high energy second-order saddle (SOS) region. This mechanism involves the molecule adopting a five-membered planar structure contrary to the traditional boat-like transition state. The SOS offers a trifurcation point where a pathway branches into three, different from the single pathway associated with a transitions state. We observe that the molecules following the SOS path exhibit distinctive dynamical features and form products with high translational energies and low rotational energies compared to those following the traditional pathways. In addition, the SOS pathway provides an alternative mechanism for the formation of stereo-selective products. Interestingly, although the reaction proceeds via a trimethylene diradical intermediate, the simulations show that the product cyclopropane is formed with a major single inversion of the configuration consistent with experimental observations. They also reveal mechanisms that do not follow the minimum energy paths and exhibit non-statistical dissociation dynamics.

Photochemical pathways in nucleobases measured with an X-ray FEL

The conversion of light energy into other molecular energetic degrees of freedom is often dominated by ultrafast, non-adiabatic processes. Femtosecond spectroscopy with optical pulses has helped in shaping our understanding of crucial processes in molecular energy-conversion. The advent of new, ultrashort and bright X-ray free electron laser sources opens the possibility to use X-ray-typical element and site sensitivity for ultrafast molecular research. We present two types of spectroscopy, ultrafast Auger and ultrafast X-ray absorption spectroscopy, and discuss their sensitivity to molecular processes. While Auger spectroscopy is able to monitor bond distance changes in the vicinity of an X-ray created core hole, near-edge absorption spectroscopy can deliver high-fidelity information on non-adiabatic transitions involving lone-pair orbitals. We demonstrate these features on the example of the UV-excited nucleobase thymine, investigated at the oxygen K-edge. We find a C-O bond elongation in the Auger data in addition to ππ*/ nπ* non-adiabatic transition in X-ray near-edge absorption. We compare the results from both methods and draw a conclusive scenario of non-adiabatic molecular relaxation after UV excitation. This article is part of the theme issue 'Measurement of ultrafast electronic and structural dynamics with X-rays'.

Roaming Dynamics in the Photodissociation of Formic Acid at 230 nm

Roaming dynamics is observed in the photodissociation of formic acid (HCOOH) at 230 nm by using the slice imaging method. In combination with rotational state selective (2 + 1) resonance-enhanced multiphoton ionization of the CO fragments, the speed distributions of the CO fragments exhibit a low recoil velocity at low rotational levels of J = 9 and 20, while the velocity distributions of CO at high rotational levels of J = 30 and 48 show a relatively large recoil velocity. The experimental results indicate that the roaming of OH radical should be related with the formation of CO + H₂O channel at the present photolysis energy. Unlike the roaming pathways occurring in H₂CO that can be described by loose flat potential, our CO speed distribution analysis suggests the presence of a "tight" flat potential in the roaming dynamics of HCOOH molecules.

Quantum Roaming in the Complex-Forming Mechanism of the Reactions of OH with Formaldehyde and Methanol at Low Temperature and Zero Pressure: A Ring Polymer Molecular Dynamics Approach

The quantum dynamics of the title reactions are studied using the ring polymer molecular dynamics (RPMD) method from 20 to 1200 K using recently proposed full dimensional potential energy surfaces which include long-range dipole-dipole interactions. A V-shaped dependence of the reaction rate constants is found with a minimum at 200-300 K, in rather good agreement with the current experimental data. For temperatures above 300 K the reaction proceeds following a direct H-abstraction mechanism. However, below 100 K the reaction proceeds via organic-molecule···OH collision complexes, with very long lifetimes, longer than 10^-7 s, associated with quantum roaming arising from the inclusion of quantum effects by the use of RPMD. The long lifetimes of these complexes are comparable to the time scale of the tunnelling to form reaction products. These complexes are formed at zero pressure because of quantum effects and not only at high pressure as suggested by transition state theory (TST) calculations for OH + methanol and other OH reactions. The zero-pressure rate constants reproduce quite well measured ones below 200 K, and this agreement opens the question of how important the pressure effects on the reaction rate constants are, as implied in TST-like formalisms. The zero-pressure mechanism is applicable only to very low gas density environments, such as the interstellar medium, which are not repeatable by experiments.

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.

Dynamics and quantum yields of H2 + CH2CO as a primary photolysis channel in CH3CHO

A new ketene + H₂ channel in CH₃CHO photolysis is not modelled by quasi-classical trajectories over the transition state.

Photoinduced C–H bond fission in prototypical organic molecules and radicals

We survey and assess current knowledge regarding the primary photochemistry of hydrocarbon molecules and radicals.

Substituent effects on H3 + formation via H2 roaming mechanisms from organic molecules under strong-field photodissociation

Roaming chemical reactions are often associated with neutral molecules. The recent findings of roaming processes in ionic species, in particular, ones that lead to the formation of H₃ ⁺ under strong-field laser excitation, are of considerable interest. Given that such gas-phase reactions are initiated by double ionization and subsequently facilitated through deprotonation, we investigate the strong-field photodissociation of ethanethiol, also known as ethyl mercaptan, and compare it to results from ethanol. Contrary to expectations, the H₃ ⁺ yield was found to be an order of magnitude lower for ethanethiol at certain laser field intensities, despite its lower ionization energy and higher acidity compared to ethanol. In-depth analysis of the femtosecond time-resolved experimental findings, supported by ab initio quantum mechanical calculations, provides key information regarding the roaming mechanisms related to H₃ ⁺ formation. Results of this study on the dynamics of dissociative half-collisions involving H₃ ⁺, a vital cation which acts as a Brønsted-Lowry acid protonating interstellar organic compounds, may also provide valuable information regarding the formation mechanisms and observed natural abundances of complex organic molecules in interstellar media and planetary atmospheres.

Post-transition state bifurcations induce dynamical detours in Pummerer-like reactions

A post-transition state bifurcation (PTSB) involved in a Pummerer-type rearrangement is characterized using density functional theory (DFT) calculations on potential energy stationary points and direct dynamics simulations. A sensitivity of the ratio of products produced via this PTSB to solvent dielectric constant is revealed and implications of such a dependence for selectivity control of organic reactions are discussed.