Bimolecular Reaction Dynamics Research Articles

Probing the dynamics and bottleneck of the key atmospheric SO2 oxidation reaction by the hydroxyl radical

SO2 (Sulfur dioxide) is the major precursor to the production of sulfuric acid (H2SO4), contributing to acid rain and atmospheric aerosols. Sulfuric acid formed from SO2 generates light-reflecting sulfate aerosol particles in the atmosphere. This property has prompted recent geoengineering proposals to inject sulfuric acid or its precursors into the Earth's atmosphere to increase the planetary albedo to counteract global warming. SO2 oxidation in the atmosphere by the hydroxyl radical HO to form HOSO2 is a key rate-limiting step in the mechanism for forming acid rain. However, the dynamics of the HO + SO2 → HOSO2 reaction and its slow rate in the atmosphere are poorly understood to date. Herein, we use photoelectron spectroscopy of cryogenically cooled HOSO2- anion to access the neutral HOSO2 radical near the transition state of the HO + SO2 reaction. Spectroscopic and dynamic calculations are conducted on the first ab initio-based full-dimensional potential energy surface to interpret the photoelectron spectra of HOSO2- and to probe the dynamics of the HO + SO2 reaction. In addition to the finding of a unique pre-reaction complex (HO⋯SO2) directly connected to the transition state, dynamic calculations reveal that the accessible phase space for the HO + SO2 → HOSO2 reaction is extremely narrow, forming a key reaction bottleneck and slowing the reaction rate in the atmosphere, despite the low reaction barrier. This study underlines the importance of understanding the full multidimensional potential energy surface to elucidate the dynamics of complex bimolecular reactions involving polyatomic reactants.

Theoretical Insights into the Dynamics of Gas-Phase Bimolecular Reactions with Submerged Barriers.

Much attention has been paid to the dynamics of both activated gas-phase bimolecular reactions, which feature monotonically increasing integral cross sections and Arrhenius kinetics, and their barrierless capture counterparts, which manifest monotonically decreasing integral cross sections and negative temperature dependence of the rate coefficients. In this Perspective, we focus on the dynamics of gas-phase bimolecular reactions with submerged barriers, which often involve radicals or ions and are prevalent in combustion, atmospheric chemistry, astrochemistry, and plasma chemistry. The temperature dependence of the rate coefficients for such reactions is often non-Arrhenius and complex, and the corresponding dynamics may also be quite different from those with significant barriers or those completely dominated by capture. Recent experimental and theoretical studies of such reactions, particularly at relatively low temperatures or collision energies, have revealed interesting dynamical behaviors, which are discussed here. The new knowledge enriches our understanding of the dynamics of these unusual reactions.

Dynamics and thermodynamics of bimolecular reactions at heterogeneous biological interfaces.

The behaviors of biomolecules at heterogeneous biological aqueous-nonaqueous interfaces markedly differ from those in bulk water. We found that various biomolecular reactions, including protein unfolding, protein-ligan binding, chlorophyll demetallation, redox reactions, and crystallization, are accelerated by a factor of 103 or higher at the interface of micron-sized water droplets (microdroplets) (PNAS, 2015; QRevPhys, 2015 and 2017). The FRET analysis revealed that molecules are highly localized at the nanoscale regime near the water surface, which is in part responsible for the observed accelerated dynamics (JPCB, 2020). Besides, the thermodynamics of biomolecular reactions are also altered near the water surface. We found that molecules are in highly ordered states near the aqueous-nonaqueous interface (Sci Adv, 2020, Sci Rep, 2020). This minimized entropy lowers entropic barriers so that thermodynamically unfavorable reactions such as phosphorylation and RNA formation can proceed without any enzyme or ATP energy source at room temperature (PNAS 2017, 2018). Moreover, the charge exchange between water molecules and biomolecules spontaneously occurs at the water interface, leading to the oxidation of water molecules and the reduction of biomolecules in the Krebs cycle (JACS, 2019; Nat Commun, 2018; PNAS, 2018, 2019; Chem Sci, 2019). Vibrational Stark spectroscopy coupled with stimulated Raman excited fluorescence microscopy showed the presence of a strong intrinsic interfacial electric field formed at the surface of the microdroplet (JPCL, 2020). This strong electric field, combined with the partial desolvation effect, can be responsible for the unusual behaviors of biomolecules at the aqueous-nonaqueous interface (Annu Rev Phys Chem, 2020). These findings suggest that the dynamics and thermodynamics of biomolecular reactions occurring in cells with various aqueous-nonaqueous interfaces and compartments may significantly differ from those in bulk water.

Observation of resonances in the transition state region of the F + NH3 reaction using anion photoelectron spectroscopy.

The transition state of a chemical reaction is a dividing surface on the reaction potential energy surface (PES) between reactants and products and is thus of fundamental interest in understanding chemical reactivity. The transient nature of the transition state presents challenges to its experimental characterization. Transition-state spectroscopy experiments based on negative-ion photodetachment can provide a direct probe of this region of the PES, revealing the detailed vibrational structure associated with the transition state. Here we study the F + NH3 → HF + NH2 reaction using slow photoelectron velocity-map imaging spectroscopy of cryogenically cooled FNH3- anions. Reduced-dimensionality quantum dynamical simulations performed on a global PES show excellent agreement with the experimental results, enabling the assignment of spectral structure. Our combined experimental-theoretical study reveals a manifold of vibrational Feshbach resonances in the product well of the F + NH3 PES. At higher energies, the spectra identify features attributed to resonances localized across the transition state and into the reactant complex that may impact the bimolecular reaction dynamics.

Bimolecular reactions of S2+ with Ar, H2 and N2: reactivity and dynamics.

The reactivity, energetics and dynamics of bimolecular reactions between S2+ and three neutral species (Ar, H2 and N2) have been studied using a position-sensitive coincidence methodology at centre-of-mass collision energies below 6 eV. This is the first study of bimolecular reactions involving S2+, a species detected in planetary ionospheres, the interstellar medium, and in anthropogenic manufacturing processes. The reactant dication beam employed consists predominantly of S2+ in the ground 3P state, but some excited states are also present. Most of the observed reactions involve the ground state of S2+, but the dissociative electron transfer reactions appear to exclusively involve excited states of this atomic dication. We observe exclusively single electron-transfer between S2+ and Ar, a process which exhibits strong forward scatting typical of the Landau-Zener style dynamics observed for other dicationic electron transfer reactions. Following collisions between S2+ + H2, non-dissociative and dissociative single electron-transfer reactions were detected. The dynamics here show evidence for the formation of a long-lived collision complex, [SH2]2+, in the dissociative single electron-transfer channel. The formation of SH+ was not observed. In contrast, the collisions of S2+ + N2 result in the formation of SN+ + N+ in addition to the products of single electron-transfer reactions.

Probing the Exit Channel of the OH + CH3OH → H2O + CH3O Reaction by Photodetachment of CH3O-(H2O).

Transition state dynamics of bimolecular reactions can be probed by photodetachment of a precursor anion when the Franck-Condon region of the corresponding neutral potential energy surface is near a saddle point. In this study, photodetachment of anions at m/z = 49 enabled investigation of the exit channel of the OH + CH3OH → H2O + CH3O reaction using photoelectron-photofragment coincidence spectroscopy. High-level coupled-cluster calculations of the stationary points on the anion surface show that the methoxide-water cluster CH3O-(H2O) is the stable minimum on the anion surface. Photodetachment at a 3.20 eV photon energy leads to long-lived H2O(CH3O) complexes and H2O + CH3O products consistent with both direct dissociative photodetachment and resonance mediated processes on the neutral surface. The partitioning of total kinetic energy in the system indicates that water stretch and bend excitation is induced in dissociative photodetachment and evidence for long-lived complexes consistent with vibrational Feshbach resonances is reported.

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.

Universal crossed beam imaging studies of polyatomic reaction dynamics.

The marriage between high level quantum calculations and experimental advances in laser technology, quantum state control, and detection techniques have opened the door to the study of molecular collision dynamics at a new level of detail. However, one current challenge lies in adapting these powerful strategies to address questions beyond the scope of the small ground state systems that have largely been the focus of reaction dynamics investigations to-date. For molecules with intermediate or large size (more than 6 atoms), lack of spectroscopic information and spectral congestion limit quantum state preparation, control and detection for experiment, and the large number of degrees of freedom of the system makes accurate quantum dynamics calculations prohibitively expensive. Nevertheless, studies of the chemical dynamics of such systems can reveal novel aspects of reactivity not anticipated based upon the behavior of smaller model systems. This Perspective will highlight applications of soft vacuum ultraviolet photoionization at 157 nm as a universal probe in combination with crossed beams and DC slice velocity map ion imaging to study bimolecular reaction dynamics of molecules of intermediate or large size, illuminated with support of high-level ab initio calculations. Here, we report on the chemical dynamics of atomic oxygen or chlorine reactions with organic compounds: propanol isomers, alkylamines (N(CH3)3 and NH(CH3)2), and isobutene ((CH3)2CCH2) studied using this approach. The polyatomic radical products from the hydrogen abstraction process have been detected by 157 nm photoionization and their slice ion images embody translational energy and angular information that directly reflect the underlying collision dynamics. Various reaction mechanisms (such as direct abstraction and addition-elimination) along with the involvement of roaming dynamics and novel intersystem crossing pathways are presented. These demonstrate the power of this technique to reveal fundamentally new aspects of reaction dynamics that arise in larger reaction systems.

Electronic spectroscopy of methyl vinyl ketone oxide: A four-carbon unsaturated Criegee intermediate from isoprene ozonolysis.

Ozonolysis of isoprene, one of the most abundant volatile organic compounds in the atmosphere, proceeds through methyl vinyl ketone oxide (MVK-oxide), methacrolein oxide, and formaldehyde oxide (CH2OO) Criegee intermediates. The present study focuses on MVK-oxide, a four-carbon unsaturated carbonyl oxide intermediate, using vacuum ultraviolet photoionization at 118 nm and UV-visible induced depletion of the m/z = 86 mass channel to characterize its first π* ← π electronic transition. The electronic spectrum is broad and unstructured with its peak at 388 nm (3.2 eV). The MVK-oxide spectrum is shifted to a significantly longer wavelength than CH2OO and alkyl-substituted Criegee intermediates studied previously due to extended conjugation across the vinyl and carbonyl oxide groups. Electronic excitation results in rapid dissociation at λ ≤ 430 nm to methyl vinyl ketone and O 1D products, the latter detected by 2 + 1 resonance enhanced multiphoton ionization using velocity map imaging. Complementary electronic structure calculations (CASPT2(12,10)/AVDZ) predict two π* ← π transitions with significant oscillator strength for each of the four conformers of MVK-oxide with vertical excitation energies (and corresponding wavelengths) in the 3.1-3.6 eV (350-400 nm) and 4.5-5.5 eV (220-280 nm) regions. The computed electronic absorption profile of MVK-oxide, based on a Wigner distribution of ground state configurations and summed over the four conformers, is predicted to peak at 397 nm. UV-visible spectroscopy on the first π* ← π transition is shown by a combination of experiment and theory to provide a sensitive method for detection of the MVK-oxide Criegee intermediate that will enable further studies of its photochemistry and unimolecular and bimolecular reaction dynamics.

Feshbach resonances in the exit channel of the F + CH3OH → HF + CH3O reaction observed using transition-state spectroscopy.

The transition state governs how chemical bonds form and cleave during a chemical reaction and its direct characterization is a long-standing challenge in physical chemistry. Transition state spectroscopy experiments based on negative-ion photodetachment provide a direct probe of the vibrational structure and metastable resonances that are characteristic of the reactive surface. Dynamical resonances are extremely sensitive to the topography of the reactive surface and provide an exceptional point of comparison with theory. Here we study the seven-atom F + CH3OH → HF + CH3O reaction using slow photoelectron velocity-map imaging spectroscopy of cryocooled CH3OHF- anions. These measurements reveal spectral features associated with a manifold of vibrational Feshbach resonances and bound states supported by the post-transition state potential well. Quantum dynamical calculations yield excellent agreement with the experimental results, allow the assignment of spectral structure and demonstrate that the key dynamics of complex bimolecular reactions can be captured with a relatively simple theoretical framework.

Taking the plunge: chemical reaction dynamics in liquids.

The dynamics of chemical reactions in liquid solutions are now amenable to direct study using ultrafast laser spectroscopy techniques and advances in computer simulation methods. The surrounding solvent affects the chemical reaction dynamics in numerous ways, which include: (i) formation of complexes between reactants and solvent molecules; (ii) modifications to transition state energies and structures relative to the reactants and products; (iii) coupling between the motions of the reacting molecules and the solvent modes, and exchange of energy; (iv) solvent caging of reactants and products; and (v) structural changes to the solvation shells in response to the changing chemical identity of the solutes, on timescales which may be slower than the reactive events. This article reviews progress in the study of bimolecular chemical reaction dynamics in solution, concentrating on reactions which occur on ground electronic states. It illustrates this progress with reference to recent experimental and computational studies, and considers how the various ways in which a solvent affects the chemical reaction dynamics can be unravelled. Implications are considered for research in fields such as mechanistic synthetic chemistry.

Perspective: Ultracold molecules and the dawn of cold controlled chemistry.

Ultracold molecules offer unprecedented opportunities for the controlled interrogation of molecular events, including chemical reactivity in the ultimate quantum regime. The proliferation of methods to create, cool, and confine them has allowed the investigation of a diverse array of molecular systems and chemical reactions at temperatures where only a single partial wave contributes. Here we present a brief account of recent progress on the experimental and theoretical fronts on cold and ultracold molecules and the opportunities and challenges they provide for a fundamental understanding of bimolecular chemical reaction dynamics.

Testing the Palma-Clary Reduced Dimensionality Model Using Classical Mechanics on the CH4 + H → CH3 + H2 Reaction.

Application of exact quantum scattering methods in theoretical reaction dynamics of bimolecular reactions is limited by the complexity of the equations of nuclear motion to be solved. Simplification is often achieved by reducing the number of degrees of freedom to be explicitly handled by freezing the less important spectator modes. The reaction cross sections obtained in reduced-dimensionality (RD) quantum scattering methods can be used in the calculation of rate coefficients, but their physical meaning is limited. The accurate test of the performance of a reduced-dimensionality method would be a comparison of the RD cross sections with those obtained in accurate full-dimensional (FD) calculations, which is not feasible because of the lack of complete full-dimensional results. However, classical mechanics allows one to perform reaction dynamics calculations using both the RD and the FD model. In this paper, an RD versus FD comparison is made for the 8-dimensional Palma-Clary model on the example of four isotopologs of the CH4 + H → CH3 + H2 reaction, which has 12 internal dimensions. In the Palma-Clary model, the only restriction is that the methyl group is confined to maintain C3v symmetry. Both RD and FD opacity and excitation functions as well as differential cross sections were calculated using the quasiclassical trajectory method. The initial reactant separation has been handled according to our one-period averaging method [ Nagy et al. J. Chem. Phys. 2016, 144, 014104 ]. The RD and FD excitation functions were found to be close to each other for some isotopologs, but in general, the RD reactivity parameters are lower than the FD reactivity parameters beyond statistical error, and for one of the isotopologs, the deviation is significant. This indicates that the goodness of RD cross sections cannot be taken for granted.

Recent Advances in Quantum Dynamics of Bimolecular Reactions.

In this review, we survey the latest advances in theoretical understanding of bimolecular reaction dynamics in the past decade. The remarkable recent progress in this field has been driven by more accurate and efficient ab initio electronic structure theory, effective potential-energy surface fitting techniques, and novel quantum scattering algorithms. Quantum mechanical characterization of bimolecular reactions continues to uncover interesting dynamical phenomena in atom-diatom reactions and beyond, reaching an unprecedented level of sophistication. In tandem with experimental explorations, these theoretical developments have greatly advanced our understanding of key issues in reaction dynamics, such as microscopic reaction mechanisms, mode specificity, product energy disposal, influence of reactive resonances, and nonadiabatic effects.

Oscillatory reaction cross sections caused by normal mode sampling in quasiclassical trajectory calculations.

The quasiclassical trajectory (QCT) method is an efficient and important tool for studying the dynamics of bimolecular reactions. In this method, the motion of the atoms is simulated classically, and the only quantum effect considered is that the initial vibrational states of reactant molecules are semiclassically quantized. A sensible expectation is that the initial ensemble of classical molecular states generated this way should be stationary, similarly to the quantum state it is supposed to represent. The most widely used method for sampling the vibrational phase space of polyatomic molecules is based on the normal mode approximation. In the present work, it is demonstrated that normal mode sampling provides a nonstationary ensemble even for a simple molecule like methane, because real potential energy surfaces are anharmonic in the reactant domain. The consequences were investigated for reaction CH4 + H → CH3 + H2 and its various isotopologs and were found to be dramatic. Reaction probabilities and cross sections obtained from QCT calculations oscillate periodically as a function of the initial distance of the colliding partners and the excitation functions are erratic. The reason is that in the nonstationary ensemble of initial states, the mean bond length of the breaking C-H bond oscillates in time with the frequency of the symmetric stretch mode. We propose a simple method, one-period averaging, in which reactivity parameters are calculated by averaging over an entire period of the mean C-H bond length oscillation, which removes the observed artifacts and provides the physically most reasonable reaction probabilities and cross sections when the initial conditions for QCT calculations are generated by normal mode sampling.

Dynamics of Complex-Forming Bimolecular Reactions: A Comparative Theoretical Study of the Reactions of H Atoms with O2((3)Σg(-)) and O2((1)Δg).

The atomic-level mechanism of the reaction of H atoms with triplet and singlet molecular oxygen, H((2)S) + O2((3)Σg(-)) → O((3)P) + OH((2)Πg) ( R1 ) and H((2)S) + O2((1)Δg) → O((3)P) + OH((2)Πg) ( R2 ) is analyzed in terms of the topology of the potential energy surfaces (PES) of the two reactions. Both PES exhibit a deep potential well corresponding to the ground and first excited electronic state of HO2. The ground-state reaction is endothermic with no barrier on either side of the well; the excited-state reaction is exothermic with a barrier in the entrance valley of the PES. The differences of the PES are manifested in properties such as the excitation functions, which show reaction R1 to be much slower and the effect of rotational excitation on reactivity, which speeds up reaction R1 and has little effect on R2 . Numerous common dynamics features arise from the presence of the deep potential well on the PES. Such are the significant role of isomerization (for example, 90% of reactive collisions in R2 involve at least one H atom transfer from one of the O atoms to the other in reaction R2 ), which is shown to give rise to a significant rotational excitation of the product OH radicals. Common is the significant sideways scattering of the products that originates from collisions in propeller-type arrangements induced by the presence of two bands of acceptance around the O2 molecule. The HO2 complex in both reactions proves to behave nonstatistically, with signatures of the dynamics in lifetime distributions, angular distributions, opacity functions, and product quantum-state distributions.

Vibrational Excitation of Both Products of the Reaction of CN Radicals with Acetone in Solution.

Transient electronic and vibrational absorption spectroscopy unravel the mechanisms and dynamics of bimolecular reactions of CN radicals with acetone in deuterated chloroform solutions. The CN radicals are produced by ultrafast ultraviolet photolysis of dissolved ICN. Two reactive forms of CN radicals are distinguished by their electronic absorption bands: “free” (uncomplexed) CN radicals, and “solvated” CN radicals that are complexed with solvent molecules. The lifetimes of the free CN radicals are limited to a few picoseconds following their photolytic production because of geminate recombination to ICN and INC, complexation with CDCl3 molecules, and reaction with acetone. The acetone reaction occurs with a rate coefficient of (8.0 ± 0.5) × 1010 M–1 s–1 and transient vibrational spectra in the C=N and C=O stretching regions reveal that both the nascent HCN and 2-oxopropyl (CH3C(O)CH2) radical products are vibrationally excited. The rate coefficient for the reaction of solvated CN with acetone is 40 times slower than for free CN, with a rate coefficient of (2.0 ± 0.9) × 109 M–1 s–1 obtained from the rise in the HCN product v1(C=N stretch) IR absorption band. Evidence is also presented for CN complexes with acetone that are more strongly bound than the CN–CDCl3 complexes because of CN interactions with the carbonyl group. The rates of reactions of these more strongly associated radicals are slower still.

Dynamics of bimolecular reactions in solution.

Mechanisms of bimolecular chemical reactions in solution are amenable to study on picosecond timescales, both by transient absorption spectroscopy and by computer simulation. The dynamics of exothermic reactions of CN radicals and of Cl and F atoms with organic solutes in commonly used solvents are contrasted with the corresponding dynamics in the gas phase. Many characteristics of the gas-phase reaction dynamics persist in solution, such as efficient energy release to specific vibrational modes of the products. However, additional complexities associated with the presence of the solvent are open to investigation. These features of liquid-phase reactions include the role of solvent-solute complexes, solvent caging, coupling of the product motions to the solvent bath, thermalization of internally excited reaction products, incipient hydrogen bond formation, and involvement of charge-separated states that arise from proton transfer.

Formation of Stabilized Ketene Intermediates in the Reaction of O(3P) with Oligo(phenylene ethynylene) Thiolate Self-Assembled Monolayers on Au(111)

We have taken steps to develop a methodology for quantifying the kinetics and dynamics of bimolecular reactions through spectroscopic monitoring of reactants and products during exposure of well-ordered self-assembled monolayers (SAMs) to supersonic beams of atomic reagents. The use of a SAM stabilizes highly energetic intermediates formed from bimolecular reactions at the vacuum/film interface due to rapid thermal equilibration with the SAM matrix that are otherwise not readily observed under single-collision conditions in the gas phase. In this paper, we will discuss the elucidation of the mechanistic details for the fundamental reaction between O(3P) and phenyl-substituted alkyne bonds by monitoring chemical and structural changes in an oligo(phenylene ethynylene) SAM reacting with O(3P) under collision conditions having specified initial reaction orientation. Utilizing time-resolved reflection–absorption infrared spectroscopy (RAIRS) and scanning tunneling microscopy (STM) under ultrahigh vacuum condi...

Perspective: Bimolecular chemical reaction dynamics in liquids

Bimolecular reactions in the gas phase exhibit rich and varied dynamical behaviour, but whether a profound knowledge of the mechanisms of isolated reactive collisions can usefully inform our understanding of reactions in liquid solutions remains an open question. The fluctuating environment in a liquid may significantly alter the motions of the reacting particles and the flow of energy into the reaction products after a transition state has been crossed. Recent experimental and computational studies of exothermic reactions of CN radicals with organic molecules indicate that many features of the gas-phase dynamics are retained in solution. However, observed differences may also provide information on the ways in which a solvent modifies fundamental chemical mechanisms. This perspective examines progress in the use of time-resolved infra-red spectroscopy to study reaction dynamics in liquids, discusses how existing theories can guide the interpretation of experimental data, and suggests future challenges for this field of research.