Dynamical, spectroscopic and computational imaging of bond breaking in photodissociation: roaming and role of conical intersections

Direct tracking of H2 roaming reaction in real time

Roaming is an unconventional type of molecular reaction where fragments, instead of immediately dissociating, remain weakly bound due to long-range Coulombic interactions. Due to its prevalence and ability to form new molecular compounds, roaming is fundamental to photochemical reactions in small molecules. However, the neutral character of the roaming fragment and its indeterminate trajectory make it difficult to identify experimentally. Here, we introduce an approach to image roaming, utilizing intense, femtosecond IR radiation combined with Coulomb explosion imaging to directly reconstruct the momentum vector of the neutral roaming H2, a precursor toH 3 + formation, in acetonitrile, CH3CN. This technique provides a kinematically complete picture of the underlying molecular dynamics and yields an unambiguous experimental signature of roaming. We corroborate these findings with quantum chemistry calculations, resolving this unique dissociative process.

Vacuum UV photodesorption of organics in the interstellar medium: an experimental study of formic acid HCOOH and methyl formate HCOOCH3-containing ices

Being a potential process that could explain gas phase abundances of so-called complex organic molecules (COMs) in the cold interstellar medium (ISM), the UV photon-induced desorption from organics-containing molecular ices has been experimentally studied. In this work, we focused on the observation of the photodesorbed products and the measurement of the associated photodesorption yields from pure and mixed molecular ices, each containing organic molecules whose detection has been achieved in the gas phase of the cold ISM, namely formic acid HCOOH and methyl formate HCOOCH3. Each molecule, in pure ice or in ice mixed with CO or water, was irradiated at 15 K with monochromatic vacuum UV photons in the 7-14 eV range using synchrotron radiation from the SOLEIL synchrotron facility, DESIRS beamline. Photodesorption yields of the intact molecules and of the photoproducts were derived as a function of the incident photon energy. Experiments have revealed that the desorbing species match the photodissociation pattern of each isolated molecule, with little influence of the kind of ice (pure or mixed in CO or H2O-rich environment). For both species, the photodesorption of the intact organics is found to be negligible in our experimental conditions, resulting in yields typically below 10-5 ejected molecules per incident photon. The results obtained on HCOOH and HCOOCH3-containing ices are similar to what has already been found for methanol-containing ices, but contrast with the case of another complex molecule, CH3CN, photodesorption of which has been recently studied. Such experimental results may be linked to the observation of COMs in protoplanetary disks, in which CH3CN is commonly observed whereas HCOOH or methanol are detected only in some sources, HCOOCH3 not being detected at all.

Roaming pathways and survival probability in real-time collisional dynamics of cold and controlled bialkali molecules

Perfectly controlled molecules are at the forefront of the quest to explore chemical reactivity at ultra low temperatures. Here, we investigate for the first time the formation of the long-lived intermediates in the time-dependent scattering of cold bialkali \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{23}\hbox {Na}^{87}$$\end{document}23Na87Rb molecules with and without the presence of infrared trapping light. During the nearly 50 nanoseconds mean collision time of the intermediate complex, we observe unconventional roaming when for a few tens of picoseconds either NaRb or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Na}_2$$\end{document}Na2 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Rb}_2$$\end{document}Rb2 molecules with large relative separation are formed before returning to the four-atom complex. We also determine the likelihood of molecular loss when the trapping laser is present during the collision. We find that at a wavelength of 1064 nm the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Na}_2\hbox {Rb}_2$$\end{document}Na2Rb2 complex is quickly destroyed and thus that the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{23}\hbox {Na}^{87}$$\end{document}23Na87Rb molecules are rapidly lost.

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

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

Roaming Reactions and Dynamics in the van der Waals Region

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

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

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

Full Dimensional Potential Energy Function and Calculation of State-Specific Properties of the CO+N2 Inelastic Processes Within an Open Molecular Science Cloud Perspective

A full dimensional Potential Energy Surface (PES) of the CO + N₂ system has been generated by extending an approach already reported in the literature and applied to N₂-N₂ (Cappelletti et al., 2008), CO₂-CO₂ (Bartolomei et al., 2012), and CO₂-N₂ (Lombardi et al., 2016b) systems. The generation procedure leverages at the same time experimental measurements and high-level ab initio electronic structure calculations. The procedure adopts an analytic formulation of the PES accounting for the dependence of the electrostatic and non-electrostatic components of the intermolecular interaction on the deformation of the monomers. In particular, the CO and N₂ molecular multipole moments and electronic polarizabilities, the basic physical properties controlling the behavior at intermediate and long-range distances of the interaction components, were made to depend on relevant internal coordinates. The formulated PES exhibits substantial advantages when used for structural and dynamical calculations. This makes it also well suited for reuse in Open Molecular Science Cloud services.

H2 roaming chemistry and the formation of H3+ from organic molecules in strong laser fields

Roaming mechanisms, involving the brief generation of a neutral atom or molecule that stays in the vicinity before reacting with the remaining atoms of the precursor, are providing valuable insights into previously unexplained chemical reactions. Here, the mechanistic details and femtosecond time-resolved dynamics of H₃⁺ formation from a series of alcohols with varying primary carbon chain lengths are obtained through a combination of strong-field laser excitation studies and ab initio molecular dynamics calculations. For small alcohols, four distinct pathways involving hydrogen migration and H₂ roaming prior to H₃⁺ formation are uncovered. Despite the increased number of hydrogens and possible combinations leading to H₃⁺ formation, the yield decreases as the carbon chain length increases. The fundamental mechanistic findings presented here explore the formation of H₃⁺, the most important ion in interstellar chemistry, through H₂ roaming occurring in ionic species.

H₂ roaming is associated with H₃⁺ formation when certain organic molecules are exposed to strong laser fields. Here, the mechanistic details and time-resolved dynamics of H₃⁺ formation from a series of alcohols were obtained and found that the product yield decreases as the carbon chain length increases.

Kinetics of the OH+HCl→H2 O+Cl reaction: Rate determining roles of stereodynamics and roaming and of quantum tunneling

The OH + HCl → H₂ O + Cl reaction is one of the most studied four-body systems, extensively investigated by both experimental and theoretical approaches. Here, as a continuation of our previous work on the OH + HBr and OH + HI reactions, which manifest an anti-Arrhenius behavior that was explained by stereodynamic and roaming effects, we extend the strategy to understand the transition to the sub-Arrhenius behavior occurring for the HCl case. As previously, we perform first-principles on-the-fly Born-Oppenheimer molecular dynamics calculations, thermalized at four temperatures (50, 200, 350, and 500 K), but this time we also apply a high-level transition-state-theory, modified to account for tunneling conditions. We find that the theoretical rate constants calculated with Bell tunneling corrections are in good agreement with extensive experimental data available for this reaction in the ample temperature range: (i) simulations show that the roles of molecular orientation in promoting this reaction and of roaming in finding the favorable path are minor than in the HBr and HI cases, and (ii) dominating is the effect of quantum mechanical penetration through the energy barrier along the reaction path on the potential energy surface. The discussion of these results provides clarification of the origin on different non-Arrhenius mechanisms observed along this series of reactions. © 2018 Wiley Periodicals, Inc.

A generalized unimolecular impulsive model for curved reaction path

This work aims to introduce a generalized impulsive model for unimolecular dissociation processes. This model allows us to take into account the curvature of the reaction path explicitly. It is a generalization of the previously developed multi-center impulsive model [P.-Y. Tsai and K.-C. Lin, J. Phys. Chem. A 119, 29 (2015)]. Several limitations of conventional impulsive models are eliminated by this study: (1) Unlike conventional impulsive models, in which a single molecular geometry is responsible for the impulse determination, the gradients on the whole dissociation path are taken into account. The model can treat dissociation pathways with large curvatures and loose saddle points. (2) The method can describe the vibrational excitation of polyatomic fragments due to the bond formation by multi-center impulse. (3) The available energy in conventional impulsive models is separated into uncoupled statistical and impulsive energy reservoirs, while the interplay between these reservoirs is allowed in the new model. (4) The quantum state correlation between fragments can be preserved in analysis. Dissociations of several molecular systems including the roaming pathways of formaldehyde, nitrate radical, acetaldehyde, and glyoxal are chosen as benchmarks. The predicted photofragment energy and vector distributions are consistent with the experimental results reported previously. In these examples, the capability of the new model to treat the curved dissociation path, loose saddle points, polyatomic fragments, and multiple-body dissociation is verified. As a cheaper computational tool with respect to ab initio on-the-fly direct dynamic simulations, this model can provide detailed information on the energy disposal, quantum state correlation, and stereodynamics in unimolecular dissociation processes.

Theories and simulations of roaming

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

Photodissociation of CH3CHO at 248 nm: identification of the channels of roaming, triple fragmentation and the transition state

Quasi-classical trajectory (QCT) calculations are performed on the molecular products CO + CH₄via the tight transition state (TS) and global minimum configurations. With the aid of this theoretical evidence, we have re-examined the experimental results published previously to clarify the controversial issue of photodissociation dynamics of CH₃CHO at 248 nm. For the CO (v = 0 and 1) bimodal rotational distributions obtained previously [K.-C. Hung, P.-Y. Tsai, H.-K. Li, and K.-C. Lin, J. Chem. Phys., 2014, 140, 064313], the low-rotational (J) component is re-assigned to the contribution of triple fragmentation (H + CO + CH₃), whereas the high-J component is ascribed to the CH₃-roaming pathway. The H-roaming pathway is not found in the calculations. Further, the QCT results have confirmed that the CO vibrational population especially at higher states and the low-energy component of CH₄ vibrational bimodality obtained experimentally are mainly produced following the TS pathway, which has never been identified before. While taking into account both the theoretical and experimental results, the ratio of the molecular products (CO(v = 1) + CH₄) obtained by the triple fragmentation/roaming/TS processes is evaluated to be 0.23 : 1 : 0.29.

Stereodirectional Origin of anti-Arrhenius Kinetics for a Tetraatomic Hydrogen Exchange Reaction: Born-Oppenheimer Molecular Dynamics for OH + HBr

Among four-atom processes, the reaction OH + HBr → H2O + Br is one of the most studied experimentally: its kinetics has manifested an unusual anti-Arrhenius behavior, namely, a marked decrease of the rate constant as the temperature increases, which has intrigued theoreticians for a long time. Recently, salient features of the potential energy surface have been characterized and most kinetic aspects can be considered as satisfactorily reproduced by classical trajectory simulations. Motivation of the work reported in this paper is the investigation of the stereodirectional dynamics of this reaction as the prominent reason for the peculiar kinetics: we started in a previous Letter ( J. Phys. Chem. Lett. 2015 , 6 , 1553 - 1558 ) a first-principles Born-Oppenheimer "canonical" molecular dynamics approach. Trajectories are step-by-step generated on a potential energy surface quantum mechanically calculated on-the-fly and are thermostatically equilibrated to correspond to a specific temperature. Here, refinements of the method permitted a major increase of the number of trajectories and the consideration of four temperatures -50, +200, +350, and +500 K, for which the sampling of initial conditions allowed us to characterize the stereodynamical effect. The role is documented of the adjustment of the reactants' mutual orientation to encounter the entrance into the "cone of acceptance" for reactivity. The aperture angle of this cone is dictated by a range of directions of approach compatible with the formation of the specific HOH angle of the product water molecule; and consistently the adjustment is progressively less effective the higher the kinetic energy. Qualitatively, this emerging picture corroborates experiments on this reaction, involving collisions of aligned and oriented molecular beams, and covering a range of energies higher than the thermal ones. The extraction of thermal rate constants from this molecular dynamics approach is discussed and the systematic sampling of the canonical ensemble is indicated as needed for quantitative comparison with the kinetic experiments.

Regulation of nonadiabatic processes in the photolysis of some carbonyl compounds

An energy scheme involving S_o → S₁ excitation, followed by dissociation channels via diabatic coupling, internal conversion, transition state and roaming mechanisms.

Stereodynamical Origin of Anti-Arrhenius Kinetics: Negative Activation Energy and Roaming for a Four-Atom Reaction

The OH + HBr → H2O + Br reaction, prototypical of halogen-atom liberating processes relevant to mechanisms for atmospheric ozone destruction, attracted frequent attention of experimental chemical kinetics: the nature of the unusual reactivity drop from low to high temperatures eluded a variety of theoretical efforts, ranking this one among the most studied four-atom reactions. Here, inspired by oriented molecular-beams experiments, we develop a first-principles stereodynamical approach. Thermalized sets of trajectories, evolving on a multidimensional potential energy surface quantum mechanically generated on-the-fly, provide a map of most visited regions at each temperature. Visualizations of rearrangements of bonds along trajectories and of the role of specific angles of reactants' mutual approach elucidate the mechanistic change from the low kinetic energy regime (where incident reactants reorient to find the propitious alignment leading to reaction) to high temperature (where speed hinders adjustment of directionality and roaming delays reactivity).

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

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