Intersystem crossing and dynamics in O(3P) + C2H4 multichannel reaction: experiment validates theory

Dynamics of the HCl + C2H5 Multichannel Reaction on a Full-Dimensional Ab Initio Potential Energy Surface

We report a full-dimensional ab initio analytical potential energy surface (PES), which accurately describes the HCl + C2H5 multichannel reaction. The new PES is developed by iteratively adding selected configurations along HCl + C2H5 quasi-classical trajectories (QCTs), thereby improving our previous Cl(2P3/2) + C2H6 PES using the Robosurfer program package. QCT simulations for the H'Cl + C2H5 reaction reveal hydrogen-abstraction, chlorine-abstraction, and hydrogen-exchange channels leading to Cl + C2H5H', H' + C2H5Cl, and HCl + C2H4H', respectively. Hydrogen abstraction dominates in the collision energy (Ecoll) range of 1-80 kcal/mol and proceeds with indirect isotropic scattering at low Ecoll and forward-scattered direct stripping at high Ecoll. Chlorine abstraction opens around 40 kcal/mol collision energy and becomes competitive with hydrogen abstraction at Ecoll = 80 kcal/mol. A restricted opening of the cone of acceptance in the Cl-abstraction reaction is found to result in the preference for a backward-scattering direct-rebound mechanism at all energies studied. Initial attack-angle distributions show mainly side-on collision preference of C2H5 for both abstraction reactions, and in the case of the HCl reactant, H/Cl-side preference for the H/Cl abstraction. For hydrogen abstraction, the collision energy transfer into the product translational and internal energy is almost equally significant, whereas in the case of chlorine abstraction, most of the available energy goes into the internal degrees of freedom. Hydrogen exchange is a minor channel with nearly constant reactivity in the Ecoll range of 10-80 kcal/mol.

First-principles mode-specific reaction dynamics

Controlling the outcome of chemical reactions by exciting specific vibrational and/or rotational modes of the reactants is one of the major goals of modern reaction dynamics studies. In the present Perspective, we focus on first-principles vibrational and rotational mode-specific dynamics computations on reactions of neutral and anionic systems beyond six atoms such as X + C2H6 [X = F, Cl, OH], HX + C2H5 [X = Br, I], OH- + CH3I, and F- + CH3CH2Cl. The dynamics simulations utilize high-level ab initio analytical potential energy surfaces and the quasi-classical trajectory method. Besides initial state specificity and the validity of the Polanyi rules, mode-specific vibrational-state assignment for polyatomic product species using normal-mode analysis and Gaussian binning is also discussed and compared with experiment.

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

Accurate fundamental invariant-neural network representation of ab initio potential energy surfaces

Highly accurate potential energy surfaces are critically important for chemical reaction dynamics. The large number of degrees of freedom and the intricate symmetry adaption pose a big challenge to accurately representing potential energy surfaces (PESs) for polyatomic reactions. Recently, our group has made substantial progress in this direction by developing the fundamental invariant-neural network (FI-NN) approach. Here, we review these advances, demonstrating that the FI-NN approach can represent highly accurate, global, full-dimensional PESs for reactive systems with even more than 10 atoms. These multi-channel reactions typically involve many intermediates, transition states, and products. The complexity and ruggedness of this potential energy landscape present even greater challenges for full-dimensional PES representation. These PESs exhibit a high level of complexity, molecular size, and accuracy of fit. Dynamics simulations based on these PESs have unveiled intriguing and novel reaction mechanisms, providing deep insights into the intricate dynamics involved in combustion, atmospheric, and organic chemistry.

Sulfur (3P) Reaction with Conjugated Dienes Gives Cyclization to Thiophenes under Single Collision Conditions

We combine crossed-beam velocity map imaging with high-level ab initio/transition state theory modeling of the reaction of S(3P) with 1,3-butadiene and isoprene under single collision conditions. For the butadiene reaction, we detect both H and H2 loss from the initial adduct, and from reaction with isoprene, we see both H loss and methyl loss. Theoretical calculations confirm these arise following intersystem crossing to the singlet surface forming long-lived intermediates. For the butadiene reaction, these lose H2 to form thiophene as the dominant channel, H to form the detected 2H-thiophenyl radical, or ethene, giving thioketene. For isoprene, additional reaction products are suggested by theory, including the observed H and methyl loss radicals, but also methyl thiophene, thioformaldehyde, and thioketene. The results for S(3P) + 1,3-butadiene, showing direct cyclization to the aromatic product and yielding few bimolecular product channels, are in striking contrast to those for the analogous O(3P) reaction.

A combined crossed molecular beam and theorerical study of the O(3P,1D) + acrylonitrile (CH2CHCN) reactions and implications for combustion and extraterrestrial environments

Acrylonitrile (CH₂CHCN) is ubiquitous in space (molecular clouds, solar-type star forming regions, and circumstellar envelopes) and is also abundant in the upper atmosphere of Titan. The reaction O(³P) + CH₂CHCN can be of relevance in the chemistry of the interstellar medium because of the abundance of atomic oxygen. The oxidation of acrylonitrile is also important in combustion as the thermal decomposition of pyrrolic and pyridinic structures present in fuel-bound nitrogen generates many nitrogen-bearing compounds, including acrylonitrile. Despite its relevance, limited information exists on this reaction. We report a combined experimental and theoretical investigation of the reactions of acrylonitrile with both ground ³P and excited ¹D atomic oxygen. From product angular and time-of-flight distributions in crossed molecular beam experiments with mass spectrometric detection at a collision energy, E_c, of 31.4 kJ mol^-1, we have identified the primary reaction products and determined their branching fractions (BFs). Theoretical calculations of the relevant triplet and singlet potential energy surfaces (PESs) were performed to interpret the experimental results and elucidate the reaction mechanism. Adiabatic statistical calculations of product BFs for the decomposition of the main triplet and singlet intermediates have been carried out. Combining the experimental and theoretical results, we conclude that the O(³P) reaction leads to two main product channels: (i) CH₂CNH (ketenimine) + CO (dominant with a BF of 0.87 ± 0.05), formed via efficient intersystem crossing from the entrance triplet PES to the underlying singlet PES, and (ii) HCOCHCN + H (minor, with a BF of 0.13 ± 0.05), occurring adiabatically on the triplet PES. Our study suggests the inclusion of this reaction as a possible destruction pathway of CH₂CHCN and a possible formation route of CH₂CNH in the interstellar medium. The O(¹D) + CH₂CHCN reaction mainly leads to the formation of CH₂CNH + CO adiabatically on the singlet PES. This result can improve models related to the chemistry of interstellar ice and cometary comas, where O(¹D) reactions can play a role. Overall, our results are expected to be useful for improving the models of combustion and extraterrestrial environments.

A highly accurate full-dimensional ab initio potential surface for the rearrangement of methylhydroxycarbene (H3C-C-OH)

We report here a full-dimensional machine learning global potential surface (PES) for the rearrangement of methylhydroxycarbene (H₃C-C-OH, 1t). The PES is trained with the fundamental invariant neural network (FI-NN) method on 91 564 ab initio energies calculated at the UCCSD(T)-F12a/cc-pVTZ level of theory, covering three possible product channels. FI-NN PES has the correct symmetry properties with respect to permutation of four identical hydrogen atoms and is suitable for dynamics studies of the 1t rearrangement. The averaged root mean square error (RMSE) is 11.4 meV. Six important reaction pathways, as well as the energies and vibrational frequencies at the stationary geometries on these pathways are accurately preproduced by our FI-NN PES. To demonstrate the capacity of the PES, we calculated the rate coefficient of hydrogen migration in -CH₃ (path A) and hydrogen migration of -OH (path B) with instanton theory on this PES. Our calculations predicted the half-life of 1t to be 95 min, which is excellent in agreement with experimental observations.

Spin-Orbit Couplings of Open-Shell Systems with Restricted Active Space Configuration Interaction

In this work we perform electronic structure calculations to unravel the origin of spin-orbit couplings (SOCs) in open-shell molecules. For that, we select systems displaying di or polyradical character, e.g., trimethylene, and analyze the changes in the magnitude of SOC constants along molecular distortions of ethylene and in the presence of intermolecular interactions between open and closed-shell moieties in the O₂-C₂H₄ system. Calculations were performed by using nonrelativistic wave functions obtained with the restricted active space configuration interaction (RASCI) method, in conjunction with a recent implementation for the calculation of SOC based on the spin-orbit mean field approximation. Our results demonstrate the suitability of RASCI in the calculation of SOCs of open-shell systems, while providing a deep understanding of the relationship between couplings and the nature of the electronic states. Moreover, we introduce a new definition of the SOC constant for the study of molecular aggregates.

Reactions O(3P, 1D) + HCCCN(X1Σ+) (Cyanoacetylene): Crossed-Beam and Theoretical Studies and Implications for the Chemistry of Extraterrestrial Environments

Cyanoacetylene (HCCCN), the first member of the cyanopolyyne family (HC_nN, where n = 3, 5, 7, ...), is of particular interest in astrochemistry being ubiquitous in space (molecular clouds, solar-type protostars, protoplanetary disks, circumstellar envelopes, and external galaxies) and also relatively abundant. It is also abundant in the upper atmosphere of Titan and comets. Since oxygen is the third most abundant element in space, after hydrogen and helium, the reaction O + HCCCN can be of relevance in the chemistry of extraterrestrial environments. Despite that, scarce information exists not only on the reactions of oxygen atoms with cyanoacetylene but with nitriles in general. Here, we report on a combined experimental and theoretical investigation of the reactions of cyanoacetylene with both ground ³P and excited ¹D atomic oxygen and provide detailed information on the primary reaction products, their branching fractions (BFs), and the overall reaction mechanisms. More specifically, the reactions of O(³P, ¹D) with HCCCN(X¹Σ⁺) have been investigated under single-collision conditions by the crossed molecular beams scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy, E_c, of 31.1 kJ/mol. From product angular and time-of-flight distributions, we have identified the primary reaction products and determined their branching fractions (BFs). Theoretical calculations of the relevant triplet and singlet potential energy surfaces (PESs) were performed to assist the interpretation of the experimental results and clarify the reaction mechanism. Adiabatic statistical calculations of product BFs for the decomposition of the main triplet and singlet intermediates have also been carried out. Merging together the experimental and theoretical results, we conclude that the O(³P) reaction is characterized by a minor adiabatic channel leading to OCCCN (cyanoketyl) + H (experimental BF = 0.10 ± 0.05), while the dominant channel (BF = 0.90 ± 0.05) occurs via intersystem crossing to the underlying singlet PES and leads to formation of ¹HCCN (cyanomethylene) + CO. The O(¹D) reaction is characterized by the same two channels, with the relative CO/H yield being slightly larger. Considering the recorded reactive signal and the calculated entrance barrier, we estimate that the rate coefficient for reaction O(³P) + HC₃N at 300 K is in the 10^-12 cm³ molec^-1 s^-1 range. Our results are expected to be useful to improve astrochemical and photochemical models. In addition, they are also relevant in combustion chemistry, because the thermal decomposition of pyrrolic and pyridinic structures present in fuel-bound nitrogen generates many nitrogen-bearing compounds, including cyanoacetylene.

Intersystem crossing in the entrance channel of the reaction of O(3P) with pyridine

Two quantum effects can enable reactions to take place at energies below the barrier separating reactants from products: tunnelling and intersystem crossing between coupled potential energy surfaces. Here we show that intersystem crossing in the region between the pre-reactive complex and the reaction barrier can control the rate of bimolecular reactions for weakly coupled potential energy surfaces, even in the absence of heavy atoms. For O(3P) plus pyridine, a reaction relevant to combustion, astrochemistry and biochemistry, crossed-beam experiments indicate that the dominant products are pyrrole and CO, obtained through a spin-forbidden ring-contraction mechanism. The experimental findings are interpreted-by high-level quantum-chemical calculations and statistical non-adiabatic computations of branching fractions-in terms of an efficient intersystem crossing occurring before the high entrance barrier for O-atom addition to the N-atom lone pair. At low to moderate temperatures, the computed reaction rates prove to be dominated by intersystem crossing.

Full-dimensional automated potential energy surface development and dynamics for the OH + C2H6 reaction

We develop a full-dimensional analytical potential energy surface (PES) for the OH + C2H6 reaction using the Robosurfer program system, which automatically (1) selects geometries from quasi-classical trajectories, (2) performs ab initio computations using a CCSD(T)-F12/triple-zeta-quality composite method, (3) fits the energies utilizing the permutationally-invariant monomial symmetrization approach, and iteratively improves the PES via steps (1)−(3). Quasi-classical trajectory simulations on the new PES reveal that hydrogen abstraction leading to H2O + C2H5 dominates in the collision energy range of 10−50 kcal/mol. The abstraction cross sections increase and the dominant mechanism shifts from rebound (small impact parameters and backward scattering) to stripping (larger impact parameters and forward scattering) with increasing collision energy as opacity functions and scattering angle distributions indicate. The abstraction reaction clearly favors side-on OH attack over O-side and the least-preferred H-side approach, whereas C2H6 behaves like a spherical object with only slight C−C-perpendicular side-on preference. Collision energy efficiently flows into the relative translation of the products, whereas product internal energy distributions show only little collision energy dependence. H2O/C2H5 vibrational distributions slightly/significantly violate zero-point energy and are nearly independent of collision energy, whereas the rotational distributions clearly blue-shift as collision energy increases.

The MD17 datasets from the perspective of datasets for gas-phase “small” molecule potentials

There has been great progress in developing methods for machine-learned potential energy surfaces. There have also been important assessments of these methods by comparing so-called learning curves on datasets of electronic energies and forces, notably the MD17 database. The dataset for each molecule in this database generally consists of tens of thousands of energies and forces obtained from DFT direct dynamics at 500 K. We contrast the datasets from this database for three “small” molecules, ethanol, malonaldehyde, and glycine, with datasets we have generated with specific targets for the potential energy surfaces (PESs) in mind: a rigorous calculation of the zero-point energy and wavefunction, the tunneling splitting in malonaldehyde, and, in the case of glycine, a description of all eight low-lying conformers. We found that the MD17 datasets are too limited for these targets. We also examine recent datasets for several PESs that describe small-molecule but complex chemical reactions. Finally, we introduce a new database, “QM-22,” which contains datasets of molecules ranging from 4 to 15 atoms that extend to high energies and a large span of configurations.

Exact Factorization Adventures: A Promising Approach for Non-Bound States

Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.

Dissection of the Multichannel Reaction O(3P) + C2H2: Differential Cross-Sections and Product Energy Distributions

The O(³P) + C₂H₂ reaction plays an important role in hydrocarbon combustion. It has two primary competing channels: H + HCCO (ketenyl) and CO + CH₂ (triplet methylene). To further understand the microscopic dynamic mechanism of this reaction, we report here a detailed quasi-classical trajectory study of the O(³P) + C₂H₂ reaction on the recently developed full-dimensional potential energy surface (PES). The entrance barrier TS1 is the rate-limiting barrier in the reaction. The translation of reactants can greatly promote reactivity, due to strong coupling with the reaction coordinate at TS1. The O(³P) + C₂H₂ reaction progress through a complex-forming mechanism, in which the intermediate HCCHO lives at least through the duration of a rotational period. The energy redistribution takes place during the creation of the long-lived high vibrationally (and rotationally) excited HCCHO in the reaction. The product energy partitioning of the two channels and CO vibrational distributions agree with experimental data, and the vibrational state distributions of all modes of products present a Boltzmann-like distribution.

Crossed-beam and theoretical studies of multichannel nonadiabatic reactions: branching fractions and role of intersystem crossing for O(3P)+1,3-butadiene

Atomic oxygen reactions can contribute significantly to the oxidation of unsaturated aliphatic and aromatic hydrocarbons. The reaction mechanism is started by electrophilic O atom addition to the unsaturated bond(s) to...

A full-dimensional ab initio potential energy and dipole moment surfaces for (NH3)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) for the ammonia dimer (NH₃)₂ are reported. The database of the PES consists of 27 736 ab initio energy points and all of these points were calculated at the UCCSD(T)-F12a/AVTZ level. The PES was fitted by using the fundamental invariant neural network (FI-NN) method that satisfies the permutational symmetry of identical atoms, and the root mean square fitting error for the PES is very small as low as 0.562 meV. The geometries for the (NH₃)₂ DMS are the same as those used for the PES and are calculated at the XYG3/AVTZ level. This PES can describe a variety of internal floppy motions, including all kinds of vibrational modes no matter intermolecular or intramolecular. The CCSD(T)-PES can dissociate correctly to two NH₃ monomers, with D_e = 1135.55 cm^-1 (13.58 kJ/mol) which agrees accurately with the 13.5 ± 0.3 kJ/mol predicted by previous work.

Vibrational mode-specific dynamics of the F(2P3/2) + C2H6 → HF + C2H5 reaction

We investigate the competing effect of vibrational and translational excitation and the validity of the Polanyi rules in the early- and negative-barrier F(²P_3/2) + C₂H₆ → HF + C₂H₅ reaction by performing quasi-classical dynamics simulations on a recently developed full-dimensional multi-reference analytical potential energy surface. The effect of five normal-mode excitations of ethane on the reactivity, the mechanism, and the post-reaction energy flow is followed through a wide range of collision energies. Promoting effects of vibrational excitations and interaction time, related to the slightly submerged barrier, are found to be suppressed by the early-barrier-induced translational enhancement, in contrast to the slightly late-barrier Cl + C₂H₆ reaction. The excess vibrational energy mostly converts into ethyl internal excitation while collision energy is transformed into product separation. The substantial reaction energy excites the HF vibration, which tends to show mode-specificity and translational energy dependence as well. With increasing collision energy, direct stripping becomes dominant over the direct rebound and indirect mechanisms, being basically independent of reactant excitation.

Crossed-Beam and Theoretical Studies of the O(3P, 1D) + Benzene Reactions: Primary Products, Branching Fractions, and Role of Intersystem Crossing

Reliable modeling of hydrocarbon oxidation relies critically on knowledge of the branching fractions (BFs) as a function of temperature (T) and pressure (p) for the products of the reaction of the hydrocarbon with atomic oxygen in its ground state, O(³P). During the past decade, we have performed in-depth investigations of the reactions of O(³P) with a variety of small unsaturated hydrocarbons using the crossed molecular beam (CMB) technique with universal mass spectrometric (MS) detection and time-of-flight (TOF) analysis, combined with synergistic theoretical calculations of the relevant potential energy surfaces (PESs) and statistical computations of product BFs, including intersystem crossing (ISC). This has allowed us to determine the primary products, their BFs, and extent of ISC to ultimately provide theoretical channel-specific rate constants as a function of T and p. In this work, we have extended this approach to the oxidation of one of the most important species involved in the combustion of aromatics: the benzene (C₆H₆) molecule. Despite extensive experimental and theoretical studies on the kinetics and dynamics of the O(³P) + C₆H₆ reaction, the relative importance of the C₆H₅O (phenoxy) + H open-shell products and of the spin-forbidden C₅H₆ (cyclopentadiene) + CO and phenol adduct closed-shell products are still open issues, which have hampered the development of reliable benzene combustion models. With the CMB technique, we have investigated the reaction dynamics of O(³P) + benzene at a collision energy (E_c) of 8.2 kcal/mol, focusing on the occurrence of the phenoxy + H and spin-forbidden C₅H₆ + CO and phenol channels in order to shed further light on the dynamics of this complex and important reaction, including the role of ISC. Concurrently, we have also investigated the reaction dynamics of O(¹D) + benzene at the same E_c. Synergistic high-level electronic structure calculations of the underlying triplet/singlet PESs, including nonadiabatic couplings, have been performed to complement and assist the interpretation of the experimental results. Statistical (RRKM)/master equation (ME) computations of the product distribution and BFs on these PESs, with inclusion of ISC, have been performed and compared to experiment. In light of the reasonable agreement between the CMB experiment, literature kinetic experimental results, and theoretical predictions for the O(³P) + benzene reaction, the so-validated computational methodology has been used to predict (i) the BF between the C₆H₅O + H and C₅H₆ + CO channels as a function of collision energy and temperature (at 0.1 and 1 bar), showing that their increase progressively favors radical (phenoxy + H)-forming over molecule (C₅H₆ + CO and phenol stabilization)-forming channels, and (ii) channel-specific rate constants as a function of T and p, which are expected to be useful for improved combustion models.

Electronic Structure and Excited States of the Collision Reaction O(3P) + C2H4: A Multiconfigurational Perspective

We present a study of the O(³P) + C₂H₄ scattering reaction, a process that takes place in the interstellar medium and is of relevance in atmospheric chemistry as well. A comprehensive investigation of the electronic properties of the system has been carried out based on multiconfigurational ab initio CASSCF/CASPT2 calculations, using a robust and consistent active space that can deliver accurate potential energy surfaces in the key regions visited by the system. The paper discloses detailed description of the primary reaction pathways and the relevant singlet and triplet excited states at the CASSCF and CASPT2 level, including an accurate description of the critical configurations, such as minima and transition states. The chosen active space and the CASSCF/CASPT2 computational protocol are assessed against coupled-cluster calculations to further check the stability and reliability of the entire multiconfigurational procedure.

Reaction Dynamics Study of the Molecular Hydrogen Loss Channel in the Elementary Reactions of Ground-State Silicon Atoms (Si(3P)) With 1- and 2-Methyl-1,3-Butadiene (C5H8)

The bimolecular gas-phase reactions involving ground-state atomic silicon (Si; ³P) and 1- and 2-methyl-1,3-butadiene were studied via crossed molecular beam experiments. Our data revealed indirect scattering dynamics through long-lived SiC₅H₈ collision complex(es) along with molecular hydrogen loss pathways, leading to facile formation of SiC₅H₆ isomer(s). We propose that the reactions of silicon with 1- and 2-methyl-1,3-butadiene possess reaction dynamics in an analogy to the silicon-1,3-butadiene system. This leads to cyclic methyl-substituted 2-methylene-1-silacyclobutene isomers via nonadiabatic reaction dynamics through intersystem crossing (ISC) from the triplet to the singlet surface in overall exoergic reactions through tight exit transition states and molecular hydrogen loss. Our study also suggests that the methyl group-although a spectator from the chemical viewpoint-can influence the disposal of the angular momentum into the rotational excitation of the final product.

Modeling Spin-Crossover Dynamics

In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.

First-Principles Reaction Dynamics beyond Six-Atom Systems

Moving beyond the six-atomic benchmark systems, we discuss the new age and future of first-principles reaction dynamics, which investigates complex, multichannel chemical reactions. We describe the methodology starting from the benchmark ab initio characterization of the stationary points, followed by full-dimensional potential energy surface (PES) developments and reaction dynamics computations. We highlight our composite ab initio approach providing benchmark stationary-point properties with subchemical accuracy, the Robosurfer program system enabling automatic PES development, and applications for the Cl + C2H6, F + C2H6, and OH- + CH3I post-six-atom reactions focusing on ab initio issues and their solutions as well as showing the excellent agreement between theory and experiment.

Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies

We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.

Full-dimensional potential energy surface for acetylacetone and tunneling splittings

New, full-dimensional potential energy surface for acetylacetone allows for description of H-tunneling dynamics and characterization of stationary points.

Theoretical Study of the Extent of Intersystem Crossing in the O(3P) + C6H6 Reaction with Experimental Validation

The extent of intersystem crossing in the O(³P) + C₆H₆ reaction, a prototypical system for spin-forbidden reactions in oxygenated aromatic molecules, is theoretically evaluated for the first time. Calculations are performed using nonadiabatic transition-state theory coupled with stochastic master equation simulations and Landau-Zener theory. It is found that the dominant intersystem crossing pathways connect the T2 and S0 potential energy surfaces through at least two distinct minimum-energy crossing points. The calculated channel-specific rate constants and intersystem crossing branching fractions differ from previous literature estimates and provide valuable kinetic data for the investigation of benzene and polycyclic aromatic hydrocarbons oxidation in interstellar, atmospheric, and combustion conditions. The theoretical results are supported by crossed molecular beam experiments with electron ionization mass-spectrometric detection and time-of-flight analysis at 8.2 kcal/mol collision energy. This system is a suitable benchmark for theoretical and experimental studies of intersystem crossing in aromatic species.

The reaction of O(3P) with alkynes: a dynamic and computational study focusing on formyl radical production

Production of formyl radical, HCO, from reactions of O(³P) with alkynes (acetylene, propyne, 1-butyne, and 1-pentyne) has been investigated using cavity ringdown laser absorption spectroscopy (CRDLAS) and computational methods. No HCO was detected from reaction with acetylene, while the amount of HCO increased for propyne and 1-butyne, dropping off somewhat for 1-pentyne. These results differ from trends previously observed for reactions of O(³P) with alkenes, which exhibit the largest HCO production for the smallest alkene and drop off as the alkene size increases. Computational studies employing density functional and coupled cluster methods have been employed to investigate the triplet and singlet state pathways for HCO production. Because intersystem crossing (ISC) has been shown to be important in these processes, the minimum energy crossing point (MECP) between the triplet and singlet surfaces has been studied. We find the MECP for propyne to possess C₁ symmetry and to lie lower in energy than previous studies have found. Natural Bond Orbital and Natural Resonance Theory analyses have been performed to investigate the changes in spin density and bond order along the reaction pathways for formation of HCO. Explanations are suggested for the trend in HCO formation observed for the alkynes. The trend in alkyne HCO yield also is compared and contrasted with the trend previously observed for the alkenes.

Internal Conversion and Intersystem Crossing with the Exact Factorization

We present a detailed derivation of the generalized coupled-trajectory mixed quantum-classical (G-CT-MQC) algorithm based on the exact-factorization equations. The ultimate goal is to propose an algorithm that can be employed for molecular dynamics simulations of nonradiative phenomena, as the spin-allowed internal conversions and the spin-forbidden intersystem crossings. Internal conversions are nonadiabatic processes driven by the kinetic coupling between electronic states, whereas intersystem crossings are mediated by the spin-orbit coupling. In this paper, we discuss computational issues related to the suitable representation for electronic dynamics and the different natures of kinetic and spin-orbit coupling. Numerical studies on model systems allow us to test the performance of the G-CT-MQC algorithm in different situations.

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.

Spin-Orbit Interactions in Ultrafast Molecular Processes

We investigate spin-orbit interactions in ultrafast molecular processes employing the exact factorization of the electron-nuclear wave function. We revisit the original derivation by including spin-orbit coupling, and show how the dynamics driven by the time-dependent potential energy surface alleviates inconsistencies arising from different electronic representations. We propose a novel trajectory-based scheme to simulate spin-forbidden non-radiative processes, and we show its performance in the treatment of excited-state dynamics where spin-orbit effects couple different spin multiplets.

Photodissociation of iso-propoxy (i-C3H7O) radical at 248 nm

Photodissociation of the i-C₃H₇O radical is investigated using fast beam photofragment translational spectroscopy.

A global ab initio potential energy surface and dynamics of the proton-transfer reaction: OH− + D2 → HOD + D−

The reaction mechanisms of OH⁻ + D₂ → HOD + D⁻ were first revealed by theory, based on an accurate full-dimensional PES.

Conical intersection-regulated intermediates in bimolecular reactions: Insights from C(1D) + HD dynamics

The importance of conical intersections (CIs) in electronically nonadiabatic processes is well known, but their influence on adiabatic dynamics has been underestimated. Here, through combined experimental and theoretical studies, we show that CIs induce a barrier and regulate conversion from a precursor metastable intermediate (CI-R) to a deep well one. This results in bond-selective activation, influencing the adiabatic dynamics markedly in the C(1D) + HD reaction. Theory is validated by experiment; quantum dynamics calculations on highly accurate ab initio potential energy surfaces yield rate coefficients and product branching ratios in excellent agreement with the experiment. Quasi-classical trajectory calculations reveal that the CI-R intermediate leads to unusual reaction mechanisms (designated as C─H activation complex conversion and cyclic complex), which are responsible for large branching ratios. We also reveal that CI-R intermediates exist in other reactive systems, and the dynamical effects uncovered here may have general significance.

First-principles dynamics of collisional intersystem crossing: resonance enhanced quenching of C(1D) by N2

Intersystem crossing is a common and important nonadiabatic process in molecular systems. Revealed with a first-principles investigation, the quenching of C(¹D) by N₂ is efficient due to multiple passages via long-lived collisional resonances, despite relatively small spin–orbit couplings.

Potential energy surface stationary points and dynamics of the F- + CH3I double inversion mechanism

Direct dynamics simulations were performed to study the S_N2 double inversion mechanism S_N2-DI, with retention of configuration, for the F^- + CH₃I reaction. Previous simulations identified a transition state (TS) structure, i.e. TS0, for the S_N2-DI mechanism, including a reaction path. However, intrinsic reaction coordinate (IRC) calculations from TS0 show it is a proton transfer (PT) TS connected to the F^-HCH₂I S_N2 pre-reaction complex and the FHCH₂I^- proton transfer post-reaction complex. Inclusion of TS0 in the S_N2-DI mechanism would thus involve non-IRC atomistic dynamics. Indeed, trajectories initiated at TS0, with random ensembles of energies as assumed by RRKM theory, preferentially form the S_N2-DI products and ∼70% follow the proposed S_N2-DI pathway from TS0 to the products. In addition, the Sudden Vector Projection (SVP) method was used to identify which CH₃I vibrational mode excitations promote access to TS0 and the S_N2-DI mechanism. Results of F^- + CH₃I simulations, with SVP specified mode excitations, are disappointing. With the CH₃ deformations of CH₃I excited, the S_N2 single inversion mechanism is the dominant pathway. If the CH stretch modes are also excited, proton transfer dominates the reaction. S_N2-DI occurs, but with a very small probability of ∼1%. The reasons behind these results are discussed.

An accurate full-dimensional potential energy surface and quasiclassical trajectory dynamics of the H + H2O2 two-channel reaction

Quasiclassical trajectory calculations reveal interesting dynamics features based on an accurate FI-NN PES for the H + H₂O₂ two-channel reaction.

The role of O(1D) in the oxidation mechanism of ethylene by iodosobenzene and other hypervalent molecules

The role of the first excited state of oxygen (¹D) is proven essential for the description of terminal iodine-oxygen chemical bonds. The description of the I-O bond as a dative one from iodine to O(¹D) provides a simple and accurate picture which explains the oxidation properties of iodosobenzene and similar in nature molecules.