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

Theoretical Study on the Dynamics and Kinetics of the Reaction of CH2OH with OH

The stochastic one-dimensional chemical master equation (CME) simulation method was used to investigate the dynamics of the reaction of CH₂OH with OH. A multiwell multichannel potential energy surface (PES) was constructed at the CCSD(T)/Aug-cc-pVTZ//CBS-QB3 and QCISD(T)/Aug-cc-pVTZ//CBS-QB3 levels of theory. The constructed PES consisted of three chemically activated intermediates and two van der Waals complexes. The fractional population analysis unraveled the role of the energized intermediates and van der Waals complexes in the early stages of this complex reaction. The CME calculations provided the phenomenological rate constants through analysis of the eigenvalues and eigenvectors of collision matrices while Leonard-Jones potential was used to model the collisions. The CME results indicated that CH₂O and H₂O were the major products, in accordance with the literature. Also, the findings declared the temperature and pressure independence of the reaction over a wide range of temperature (250 to 2400 K) and pressure (0.1 to 7 atm). Furthermore, the efficiency of tunneling on the hydrogen transfer isomerization reaction of trans-HCOH to CH₂O was confirmed over the temperature range of 250 to 3000 K. The rate constants for different reaction channels are reported.

Double Photodetachment of F-·H2O: Experimental and Theoretical Studies of [F·H2O]. J Phys Chem Lett 2018;9:6808-6813. [PMID: 30433784 DOI: 10.1021/acs.jpclett.8b02562]

Double photodetachment of the cluster F^-·H₂O in a strong laser field is explored in a combined experimental-theoretical study. Products are observed experimentally by coincidence photofragment imaging following double ionization by intense laser pulses. Theoretically, equation of motion coupled cluster calculations (EOM-CC), suitable for modeling strong correlation effects in the electronic wave function, shed light on the Franck-Condon region, and ab initio molecular dynamics simulations also performed using EOM-CC methods reveal the fragmentation dynamics in time on the lowest-lying singlet and triplet states of [F·H₂O]⁺. The simulations show the formation of H₂O⁺ + F, which is the predominant experimentally observed product channel. Suggestions are proposed for the formation mechanisms of the minor products, for example, the very interesting H₂F⁺, which involves significant geometrical rearrangement. Analysis of the results suggests interesting future directions for the exploration of photodetachment of anionic clusters in an intense laser field.

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.

Quasiclassical trajectory calculations of CD3CHO dissociation to CD2H + DCO on a global potential energy surface

We present a quasiclassical trajectory study of the photodissociation of CD3CHO based on a global ab initio-based potential energy surface. Calculations are performed at the total energy corresponding to the photolysis wavelength of 280[Formula: see text]nm. In addition to the major radical and molecular products, CD[Formula: see text] and CD3H [Formula: see text] CO, respectively, this paper focuses on the unusual radical channel CD2H [Formula: see text] DCO, which requires a D/H exchange process before the conventional C–C bond cleavage. Five D/H exchange mechanisms are reported, which are related to the isomerizations from acetaldehyde to vinyl alcohol and back, to oxirane and back, and to the intermediate (CD–CHD–OD) and back. These D/H exchange mechanisms are in good agreement with the experimental findings [Heazlewood BR, Maccarone AT, Andrews DU, Osborn DL, Harding LB, Klippenstein SJ, Jordan MJT, Kable SH, Nat Chem 3:443, 2011].

Invited Review Article: Photofragment imaging

Photodissociation studies in molecular beams that employ position-sensitive particle detection to map product recoil velocities emerged thirty years ago and continue to evolve with new laser and detector technologies. These powerful methods allow application of tunable laser detection of single product quantum states, simultaneous measurement of velocity and angular momentum polarization, measurement of joint product state distributions for the detected and undetected products, coincident detection of multiple product channels, and application to radicals and ions as well as closed-shell molecules. These studies have permitted deep investigation of photochemical dynamics for a broad range of systems, revealed new reaction mechanisms, and addressed problems of practical importance in atmospheric, combustion, and interstellar chemistry. This review presents an historical overview, a detailed technical account of the range of methods employed, and selected experimental highlights illustrating the capabilities of the method.

Identifying the Multielectron Effect on Chemical Bond Rearrangement of CH3Cl Molecules in Strong Laser Fields

Strong field double ionization that triggers the chemical bond rearrangement of CH₃Cl is investigated by impulsive control of the alignment of molecules. The alignment and laser intensity dependent H₂⁺ and H₃⁺ yields in linearly polarized femtosecond laser have been measured, and the obtained data show that the maximum signal of H₂⁺ appears at the laser polarization parallel to the C-Cl axis of molecules and H₃⁺ species are more likely to eject at the laser polarization parallel to the C-Cl axis at low laser intensity while the H₃⁺ signal peaks at laser polarization perpendicular to the C-Cl axis at high laser intensity. The measurements indicate that electrons from HOMO - 1 and HOMO - 2 orbitals have been ionized for the generation of bond rearrangement at different laser intensity. Our results demonstrate the importance of multielectron effects and also provide an effective control method in the process of chemical bond rearrangement of the molecules in strong laser fields.

Low temperature reaction dynamics for CH3OH + OH collisions on a new full dimensional potential energy surface

Is the rise of the rate constant measured in laval expansion experiments of OH with organic molecules at low temperatures due to the reaction between the reactants or due to the formation of complexes with the buffer gas? This question has importance for understanding the evolution of prebiotic molecules observed in different astrophysical objects. Among these molecules methanol is one of the most widely observed, and its reaction with OH has been studied by several groups showing a fast increase in the rate constant under 100 K. Transition state theory doesn't reproduce this behavior and here dynamical calculations are performed on a new full dimensional potential energy surface developed for this purpose. The calculated classical reactive cross sections show an increase at low collision energies due to a complex forming mechanism. However, the calculated rate constant at temperatures below 100 K remains lower than the observed one. Quantum effects are likely responsible for the measured behavior at low temperatures.

Femtosecond photoexcitation dynamics inside a quantum solvent

The observation of chemical reactions on the time scale of the motion of electrons and nuclei has been made possible by lasers with ever shortened pulse lengths. Superfluid helium represents a special solvent that permits the synthesis of novel classes of molecules that have eluded dynamical studies so far. However, photoexcitation inside this quantum solvent triggers a pronounced response of the solvation shell, which is not well understood. Here, we present a mechanistic description of the solvent response to photoexcitation of indium (In) dopant atoms inside helium nanodroplets (He_N), obtained from femtosecond pump–probe spectroscopy and time-dependent density functional theory simulations. For the In–He_N system, part of the excited state electronic energy leads to expansion of the solvation shell within 600 fs, initiating a collective shell oscillation with a period of about 30 ps. These coupled electronic and nuclear dynamics will be superimposed on intrinsic photoinduced processes of molecular systems inside helium droplets.

Femtosecond laser spectroscopy has contributed to our understanding of structure and function of matter. Here, the authors explore the applicability of superfluid helium nanodroplets as a sample preparation method that allows investigation of previously inaccessible classes of tailor-made or fragile molecular systems.

Influence of mass and potential energy surface geometry on roaming in Chesnavich's CH4+ model

Chesnavich's model Hamiltonian for the reaction CH4+→ CH3+ + H is known to exhibit a range of interesting dynamical phenomena including roaming. The model system consists of two parts: a rigid, symmetric top representing the CH3+ ion and a free H atom. We study roaming in this model with focus on the evolution of geometrical features of the invariant manifolds in phase space that govern roaming under variations of the mass of the free atom m and a parameter a that couples radial and angular motion. In addition, we establish an upper bound on the prominence of roaming in Chesnavich's model. The bound highlights the intricacy of roaming as a type of dynamics on the verge between isomerisation and nonreactivity as it relies on generous access to the potential wells to allow reactions as well as a prominent area of high potential that aids sufficient transfer of energy between the degrees of freedom to prevent isomerisation.

Conical-intersection dynamics and ground-state chemistry probed by extreme-ultraviolet time-resolved photoelectron spectroscopy

Time-resolved photoelectron spectroscopy (TRPES) is a useful approach to elucidate the coupled electronic-nuclear quantum dynamics underlying chemical processes, but has remained limited by the use of low photon energies. Here, we demonstrate the general advantages of XUV-TRPES through an application to NO₂, one of the simplest species displaying the complexity of a non-adiabatic photochemical process. The high photon energy enables ionization from the entire geometrical configuration space, giving access to the true dynamics of the system. Specifically, the technique reveals dynamics through a conical intersection, large-amplitude motion and photodissociation in the electronic ground state. XUV-TRPES simultaneously projects the excited-state wave packet onto many final states, offering a multi-dimensional view of the coupled electronic and nuclear dynamics. Our interpretations are supported by ab initio wavepacket calculations on new global potential-energy surfaces. The presented results contribute to establish XUV-TRPES as a powerful technique providing a complete picture of ultrafast chemical dynamics from photoexcitation to the final products.

Selective hydrogenation of 1,3-butadiene catalyzed by a single Pd atom anchored on graphene: the importance of dynamics

The active-site structure, reaction mechanism, and product selectivity of the industrially important selective hydrogenation of 1,3-butadiene are investigated using first principles for an emerging single-atom Pd catalyst anchored on graphene. Density functional theory calculations suggest that the mono-π-adsorbed reactant undergoes sequential hydrogenation by Pd-activated H2. Importantly, the high selectivity towards 1-butene is attributed to the post-transition-state dynamics in the second hydrogenation step, which leads exclusively to the desorption of the product. This dynamical event prevails despite the existence of energetically preferred 1-butene adsorption on Pd, which would eventually lead to complete hydrogenation to butane and be thus inconsistent with experimental observations. This insight underscores the importance of dynamics in heterogeneous catalysis, which has so far been underappreciated.

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.

Sub-500 fs electronically nonadiabatic chemical dynamics of energetic molecules from the S1 excited state: Ab initio multiple spawning study

Energetic materials store a large amount of chemical energy. Different ignition processes, including laser ignition and shock or compression wave, initiate the energy release process by first promoting energetic molecules to the electronically excited states. This is why a full understanding of initial steps of the chemical dynamics of energetic molecules from the excited electronic states is highly desirable. In general, conical intersection (CI), which is the crossing point of multidimensional electronic potential energy surfaces, is well established as a controlling factor in the initial steps of chemical dynamics of energetic molecules following their electronic excitations. In this article, we have presented different aspects of the ultrafast unimolecular relaxation dynamics of energetic molecules through CIs. For this task, we have employed ab initio multiple spawning (AIMS) simulation using the complete active space self-consistent field (CASSCF) electronic wavefunction and frozen Gaussian-based nuclear wavefunction. The AIMS simulation results collectively reveal that the ultrafast relaxation step of the best energetic molecules (which are known to exhibit very good detonation properties) is completed in less than 500 fs. Many, however, exhibit sub-50 fs dynamics. For example, nitro-containing molecules (including C-NO₂, N-NO₂, and O-NO₂ active moieties) relax back to the ground state in approximately 40 fs through similar (S₁/S₀)_CI conical intersections. The N₃-based energetic molecule undergoes the N₂ elimination process in 40 fs through the (S₁/S₀)_CI conical intersection. Nitramine-Fe complexes exhibit sub-50 fs Fe-O and N-O bond dissociation through the respective (S₁/S₀)_CI conical intersection. On the other hand, tetrazine-N-oxides, which are known to exhibit better detonation properties than tetrazines, undergo internal conversion in a 400-fs time scale, while the relaxation time of tetrazine is very long (about 100 ns). Many other characteristics of sub-500 fs nonadiabatic decay of energetic molecules are discussed. In the end, many unresolved issues associated with the ultrafast nonadiabatic chemical dynamics of energetic molecules are presented.

What is special about how roaming chemical reactions traverse their potential surfaces? Differences in geodesic paths between roaming and non-roaming events

With the notable exception of some illustrative two-degree-of-freedom models whose surprising classical dynamics has been worked out in detail, theories of roaming have largely bypassed the issue of when and why the counterintuitive phenomenon of roaming occurs. We propose that a useful way to begin to address these issues is to look for the geodesic (most efficient) pathways through the potential surfaces of candidate systems. Although roaming manifests itself in an unusual behavior at asymptotic geometries, we found in the case of formaldehyde dissociation that it was the pathways traversing the parts of the potential surface corresponding to highly vibrationally excited reactants that were the most revealing. An examination of the geodesics for roaming pathways in this region finds that they are much less tightly defined than the geodesics in that same region that lead directly to dissociation (whether into closed-shell products or into radical products). Thus, the broader set of options available to the roaming channel gives it an entropic advantage over more conventional reaction channels. These observations suggest that what leads to roaming in other systems may be less the presence of a localized "roaming transition state," than the existence of an entire region of the potential surface conducive to multiple equivalent pathways.

Liberation of H2 from (o-C6H4Me)3P-H(+) + (-)H-B(p-C6F4H)3 ion-pair: A transition-state in the minimum energy path versus the transient species in Born-Oppenheimer molecular dynamics

Using Born-Oppenheimer molecular dynamics (BOMD) with density functional theory, transition-state (TS) calculations, and the quantitative energy decomposition analysis (EDA), we examined the mechanism of H₂-liberation from LB-H⁽⁺⁾ + ^(-)H-LA ion-pair, 1, in which the Lewis base (LB) is (o-C₆H₄Me)₃P and the Lewis acid (LA) is B(p-C₆F₄H)₃. BOMD simulations indicate that the path of H₂ liberation from the ion-pair 1 goes via the short-lived transient species, LB⋯H₂⋯LA, which are structurally reminiscent of the TS-structure in the minimum-energy-path describing the reversible reaction between H₂ and (o-C₆H₄Me)₃P/B(p-C₆F₄H)₃ frustrated Lewis pair (FLP). With electronic structure calculations performed on graphics processing units, our BOMD data-set covers more than 1 ns of evolution of the ion-pair 1 at temperature T ≈ 400 K. BOMD simulations produced H₂-recombination events with various durations of H₂ remaining fully recombined as a molecule within a LB/LA attractive "pocket"-from very short vibrational-time scale to time scales in the range of a few hundred femtoseconds. With the help of perturbational approach to trajectory-propagation over a saddle-area, we directly examined dynamics of H₂-liberation. Using EDA, we elucidated interactions between the cationic and anionic fragments in the ion-pair 1 and between the molecular fragments in the TS-structure. We have also considered a model that qualitatively takes into account the potential energy characteristics of H-H recombination and H₂-release plus inertia of molecular motion of the (o-C₆H₄Me)₃P/B(p-C₆F₄H)₃ FLP.

Finite slice analysis (FINA)-A general reconstruction method for velocity mapped and time-sliced ion imaging

Since the advent of ion imaging, one of the key issues in the field has been creating methods to reconstruct the initial 3D distribution of particles from its 2D projection. This has led to the development of a number of different numerical methods and fitting techniques to solve this fundamental issue in imaging. In recent years, slice-imaging methods have been developed that permit direct recording of the 3D distribution, i.e., a thin slice of the recoiling fragment distribution. However, in practice, most slice imaging experiments achieve a velocity slice width of around 10%-25% around the center of the distribution. This still carries significant out-of-plane elements that can blur the spectrum, lose fine resolution, and underestimate the contribution from slow recoiling products. To overcome these limitations, we developed a new numerical method to remove these out-of-plane elements from a sliced image. The finite sliced analysis method models the off-axis elements of the 3D particle distribution through the use of radial basis functions. Once applied, the method reconstructs the underlying central slice of the 3D particle distribution. The approach may be applied to arbitrarily sliced or unsliced data and has the further advantage that it neither requires nor enforces full cylindrical symmetry of the data. We demonstrate this reconstruction approach with a broad range of synthetic and experimental data that, at the same time, allows us to examine the impact of finite slicing on the recovered distributions in detail.

Semiclassical initial value theory of rotationally inelastic scattering: Some remarks on the phase index in the interaction picture

This paper deals with the treatment of quantum interferences in the semiclassical initial value theory of rotationally inelastic scattering in the interaction picture. Like many semiclassical methods, the previous approach involves a phase index related to sign changes of a Jacobian whose square root is involved in the calculations. It is shown that replacing the original phase index by a new one extends the range of applicability of the theory. The resulting predictions are in close agreement with exact quantum scattering results for a model of atom-rigid diatom collision involving strong interferences. The developments are performed within the framework of the planar rotor model, but are readily applicable to three-dimensional collisions.

Quasi-Classical Trajectory Dynamics Study of the Cl(2P) + C2H6 → HCl(v,j) + C2H5 Reaction. Comparison with Experiment

To understand and simulate the dynamics behavior of the title reaction, QCT calculations were performed on a recently developed global analytical potential energy surface, PES-2017. These calculations combine the classical description of the dynamics with pseudoquantization in the reactants and products to perform a theoretical/experimental comparison on the same footing. Thus, in the products a series of constraints are included to analyze the HCl(v = 0,j) product, which is experimentally detected. At collision energies of 5.5 and 6.7 kcal mol^-1 the largest fraction of available energy is deposited as translation, 67%, while the ethyl radical shows significant internal energy, 27%, and so it does not act as a spectator of the reaction, thus reproducing recent experimental evidence. The HCl(v=0, j) rotational distribution is cold, peaking at j = 2, only one unit hotter than experiment, which represents an error of 0.12 kcal mol^-1. At a collision energy of 5.5 kcal mol^-1 product translational distribution is slightly hotter than experiment, but at 6.7 kcal mol^-1 agreement with recent experiments is practically quantitative, suggesting that the first experiments should be revised. In addition, we observe that the HCl(v=0, j) scattering distribution shifts from isotropic at low values of j to backward at high values of j, which is in agreement with experimental data. Finally, no evidence was found for the "chattering" mechanism suggested to explain the low translational energy of the HCl product in the backward scattering region. In sum, agreement with experiments of a series of sensible dynamic properties permits us to be optimistic on the quality and accuracy of the theoretical tools used in the present work, QCT and PES-2017.

The phase space geometry underlying roaming reaction dynamics

Recent studies have found an unusual way of dissociation in formaldehyde. It can be characterized by a hydrogen atom that separates from the molecule, but instead of dissociating immediately it roams around the molecule for a considerable amount of time and extracts another hydrogen atom from the molecule prior to dissociation. This phenomenon has been coined roaming and has since been reported in the dissociation of a number of other molecules. In this paper we investigate roaming in Chesnavich's CH 4 + model. During dissociation the free hydrogen must pass through three phase space bottleneck for the classical motion, that can be shown to exist due to unstable periodic orbits. None of these orbits is associated with saddle points of the potential energy surface and hence related to transition states in the usual sense. We explain how the intricate phase space geometry influences the shape and intersections of invariant manifolds that form separatrices, and establish the impact of these phase space structures on residence times and rotation numbers. Ultimately we use this knowledge to attribute the roaming phenomenon to particular heteroclinic intersections.

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.

Ab initio studies on the photodissociation dynamics of the 1,1-difluoroethyl radical

Born-Oppenheimer molecular dynamics trajectory calculations at the HCTH147/6-31G** level of theory simulate the dissociation dynamics of photolytically excited 1,1-difluoroethyl radicals. EOMCCSD/AUG-cc-pVDZ calculations show that an excitation energy of 94.82 kcal/mol is necessary to initiate photodissociation reactions. In contrast to photodissociation dynamics of ethyl radicals where a large discrepancy between actual dissociation rates and rates that are predicted by statistical rate theories, we find reaction rates of 5.1 × 10¹¹ s^-1 for the dissociation of an H atom, which is in perfect accord with what is predicted by Rice-Ramsperger-Kassel-Marcus (RRKM) calculations and there is no indication of any nonstatistical effects. However, our trajectory calculations show a much larger fraction of C-C bond breakage reaction of 56% occurring than that expected by RRKM (only 16%).

Full dimensional potential energy surface and low temperature dynamics of the H2CO + OH → HCO + H2O reaction

A new method is proposed to analytically represent the potential energy surface of reactions involving polyatomic molecules capable of accurately describing long-range interactions and saddle points, needed to describe low-temperature collisions. It is based on two terms, a reactive force field term and a many-body term. The reactive force field term accurately describes the fragments, long-range interactions among them and the saddle points for reactions. The many-body term increases the desired accuracy everywhere else. This method has been applied to the OH + H2CO → H2O + HCO reaction, giving a barrier of 27.4 meV. The simulated classical rate constants with this potential are in good agreement with recent experimental results [Ocaña et al., Astrophys. J., 2017, submitted], showing an important increase at temperatures below 100 K. The reaction mechanism is analyzed in detail here, and explains the observed behavior at low energy by the formation of long-lived collision complexes, with roaming trajectories, with a capture observed for very long impact parameters, >100 a.u., determined by the long-range dipole-dipole interaction.

Indirect dynamics in SN2@N: insight into the influence of central atoms

Central atoms have a significant influence on the reaction kinetics and dynamics of nucleophilic substitution (S_N2). Herein, atomistic dynamics of a prototype S_N2@N reaction F^- + NH₂Cl is uncovered employing direct dynamics simulations that show strikingly distinct features from those determined for a S_N2@C congener F^- + CH₃Cl. Indirect scattering is found to prevail, which proceeds predominantly through a hydrogen-bonded F^--HNHCl complex in the reactant entrance channel. This unexpected finding of a pronounced contribution of indirect reaction dynamics, even at a high collision energy, is in strong contrast to a general evolution from indirect to direct dynamics with enhanced energy that characterizes S_N2@C. This result suggests that the relative importance of different atomic-level mechanisms may depend essentially on the interaction potential of reactive encounters and the coupling between inter- and intramolecular modes of the pre-reaction complex. For F^- + NH₂Cl the proton transfer pathway is less competitive than S_N2. A remarkable finding is that the more favorable energetics for NH₂Cl proton transfer, as compared to that for CH₃Cl, does not manifest itself in the reaction dynamics. The present work sheds light on the underlying reaction dynamics of S_N2@N, which remain largely unclear compared to well-studied S_N2@C.

Photodissociation dynamics of acetone studied by time-resolved ion imaging and photofragment excitation spectroscopy

The photodissociation dynamics of acetone has been investigated using velocity-map ion imaging and photofragment excitation (PHOFEX) spectroscopy across a range of wavelengths spanning the first absorption band (236-308 nm). The radical products of the Norrish Type I dissociation, methyl and acetyl, as well as the molecular product ketene have been detected by single-photon VUV ionization at 118 nm. Ketene appears to be formed with non-negligible yield at all wavelengths, with a maximum value of Φ ≈ 0.3 at 280 nm. The modest translational energy release is inconsistent with dissociation over high barriers on the S₀ surface, and ketene formation is tentatively assigned to a roaming pathway involving frustrated dissociation to the radical products. Fast-moving radical products are detected at λ ≤ 305 nm with total translational energy distributions that extend to the energetic limit, consistent with dissociation occurring near-exclusively on the T₁ surface following intersystem crossing. At energies below the T₁ barrier a statistical component indicative of S₀ dissociation is observed, although dissociation via the S₁/S₀ conical intersection is absent at shorter wavelengths, in contrast to acetaldehyde. The methyl radical yield is enhanced over that of acetyl in PHOFEX spectra at λ ≤ 260 nm due to the onset of secondary dissociation of internally excited acetyl radicals. Time-resolved ion imaging experiments using picosecond duration pulses at 266 nm find an appearance time constant of τ = 1490 ± 140 ps for CH₃ radicals formed on T₁. The associated rate is representative of S₁ → T₁ intersystem crossing. At 284 nm, CH₃ is formed on T₁ with two distinct timescales: a fast <10 ns component is accompanied by a slower component with τ = 42 ± 7 ns. A two-step mechanism involving fast internal conversion, followed by slower intersystem crossing (S₁ → S₀ → T₁) is proposed to explain the slow component.

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.

Dynamic exit-channel pathways of the microsolvated HOO-(H2O) + CH3Cl SN2 reaction: Reaction mechanisms at the atomic level from direct chemical dynamics simulations

Microsolvated bimolecular nucleophilic substitution (S_N2) reaction of monohydrated hydrogen peroxide anion [HOO^-(H₂O)] with methyl chloride (CH₃Cl) has been investigated with direct chemical dynamics simulations at the M06-2X/6-31+G(d,p) level of theory. Dynamic exit-channel pathways and corresponding reaction mechanisms at the atomic level are revealed in detail. Accordingly, a product distribution of 0.85:0.15 is obtained for Cl^-:Cl^-(H₂O), which is consistent with a previous experiment [D. L. Thomsen et al. J. Am. Chem. Soc. 135, 15508 (2013)]. Compared with the HOO^- + CH₃Cl S_N2 reaction, indirect dynamic reaction mechanisms are enhanced by microsolvation for the HOO^-(H₂O) + CH₃Cl S_N2 reaction. On the basis of our simulations, further crossed molecular beam imaging experiments are highly suggested for the S_N2 reactions of HOO^- + CH₃Cl and HOO^-(H₂O) + CH₃Cl.

Time-Resolved Kinetic Chirped-Pulse Rotational Spectroscopy in a Room-Temperature Flow Reactor

Chirped-pulse Fourier transform millimeter-wave spectroscopy is a potentially powerful tool for studying chemical reaction dynamics and kinetics. Branching ratios of multiple reaction products and intermediates can be measured with unprecedented chemical specificity; molecular isomers, conformers, and vibrational states have distinct rotational spectra. Here we demonstrate chirped-pulse spectroscopy of vinyl cyanide photoproducts in a flow tube reactor at ambient temperature of 295 K and pressures of 1-10 μbar. This in situ and time-resolved experiment illustrates the utility of this novel approach to investigating chemical reaction dynamics and kinetics. Following 193 nm photodissociation of CH₂CHCN, we observe rotational relaxation of energized HCN, HNC, and HCCCN photoproducts with 10 μs time resolution and sample the vibrational population distribution of HCCCN. The experimental branching ratio HCN/HCCCN is compared with a model based on RRKM theory using high-level ab initio calculations, which were in turn validated by comparisons to Active Thermochemical Tables enthalpies.

Flash Pyrolysis of t-Butyl Hydroperoxide and Di-t-butyl Peroxide: Evidence of Roaming in the Decomposition of Organic Hydroperoxides

Thermal decomposition of t-butyl hydroperoxide and di-t-butyl peroxide was investigated using flash pyrolysis (in a short reaction time of <100 μs) and vacuum-ultraviolet (λ = 118.2 nm) single-photon ionization time-of-flight mass spectrometry (VUV-SPI-TOFMS) at temperatures up to 1120 K and quantum computational methods. Acetone and methyl radical were detected as the predominant products in the initial decomposition of di-t-butyl peroxide via O-O bond fission. In the initial dissociation of t-butyl hydroperoxide, acetone, methyl radical, isobutylene, and isobutylene oxide products were identified. The novel detection of the unimolecular formation of isobutylene oxide, as supported by the computational study, was found to proceed via a roaming hydroxyl radical facilitated by a hydrogen-bonded intermediate. This new pathway could provide a new class of reactions to consider in the modeling of the low temperature oxidation of alkanes.

Empirical Classification of Trajectory Data: An Opportunity for the Use of Machine Learning in Molecular Dynamics

Classical Hamiltonian trajectories are initiated at random points in phase space on a fixed energy shell of a model two degrees of freedom potential, consisting of two interacting minima in an otherwise flat energy plane of infinite extent. Below the energy of the plane, the dynamics are demonstrably chaotic. However, most of the work in this paper involves trajectories at a fixed energy that is 1% above that of the plane, in which regime the dynamics exhibit behavior characteristic of chaotic scattering. The trajectories are analyzed without reference to the potential, as if they had been generated in a typical direct molecular dynamics simulation. The questions addressed are whether one can recover useful information about the structures controlling the dynamics in phase space from the trajectory data alone, and whether, despite the at least partially chaotic nature of the dynamics, one can make statistically meaningful predictions of trajectory outcomes from initial conditions. It is found that key unstable periodic orbits, which can be identified on the analytical potential, appear by simple classification of the trajectories, and that the specific roles of these periodic orbits in controlling the dynamics are also readily discerned from the trajectory data alone. Two different approaches to predicting trajectory outcomes from initial conditions are evaluated, and it is shown that the more successful of them has ∼90% success. The results are compared with those from a simple neural network, which has higher predictive success (97%) but requires the information obtained from the "by-hand" analysis to achieve that level. Finally, the dynamics, which occur partly on the very flat region of the potential, show characteristics of the much-studied phenomenon called "roaming." On this potential, it is found that roaming trajectories are effectively "failed" periodic orbits and that angular momentum can be identified as a key controlling factor, despite the fact that it is not a strictly conserved quantity. It is also noteworthy that roaming on this potential occurs in the absence of a "roaming saddle," which has previously been hypothesized to be a necessary feature for roaming to occur.

Recoil Inversion in the Photodissociation of Carbonyl Sulfide near 234 nm

We report the observation of recoil inversion of the CO (v=0, J_{CO}=66) state in the UV dissociation of lab-frame oriented carbonyl sulfide (OCS). This state is ejected in the opposite direction with respect to all other (>30) states and in absence of any OCS rotation, thus resulting in spatial filtering of this particular high-J rovibrational state. This inversion is caused by resonances occurring in shallow local minima of the molecular potential, which bring the sulfur closer to the oxygen than the carbon atom, and is a striking example where such subtleties severely modify the photofragment trajectories. The resonant behavior is observed only in the photofragment trajectories and not in their population, showing that stereodynamic measurements from oriented molecules offer an indispensable probe for exploring energy landscapes.

Development of a linear-type double reflectron for focused imaging of photofragment ions from mass-selected complex ions

An ion imaging apparatus with a double linear reflectron mass spectrometer has been developed, in order to measure velocity and angular distributions of mass-analyzed fragment ions produced by photodissociation of mass-selected gas phase complex ions. The 1st and the 2nd linear reflectrons were placed facing each other and controlled by high-voltage pulses in order to perform the mass-separation of precursor ions in the 1st reflectron and to observe the focused image of the photofragment ions in the 2nd reflectron. For this purpose, metal meshes were attached on all electrodes in the 1st reflectron, whereas the mesh was attached only on the last electrode in the 2nd reflectron. The performance of this apparatus was evaluated using imaging measurement of Ca⁺ photofragment ions from photodissociation reaction of Ca⁺Ar complex ions at 355 nm photoexcitation. The focused ion images were obtained experimentally with the double linear reflectron at the voltages of the reflection electrodes close to the predictions by ion trajectory simulations. The velocity and angular distributions of the produced Ca⁺ ([Ar] 4p¹, ²P_3/2) ion were analyzed from the observed images. The binding energy D₀ of Ca⁺Ar in the ground state deduced in the present measurement was consistent with those determined theoretically and by spectroscopic measurements. The anisotropy parameter β of the transition was evaluated for the first time by this instrument.

A new (multi-reference configuration interaction) potential energy surface for H2CO and preliminary studies of roaming

We report a new global potential energy surface (PES) for H2CO, based on precise fitting of roughly 67 000 MRCI/cc-pVTZ energies. This PES describes the global minimum, the cis- and trans-HCOH isomers, and barriers relevant to isomerization, formation of the molecular (H2+CO) and radical (H+HCO) products, and the loose so-called roaming transition-state saddle point. The key features of the PES are reviewed and compared with a previous PES, denoted by PES04, based on five local fits that are 'stitched' together by switching functions (Zhang et al. 2004 J. Phys. Chem. A108, 8980-8986 (doi:10.1021/jp048339l)). Preliminary quasi-classical trajectory calculations are performed at the total energy of 36 233 cm-1 (103 kcal mol-1), relative to the H2CO global minimum, using the new PES, with a particular focus on roaming dynamics. When compared with the results from PES04, the new PES findings show similar rotational distributions, somewhat more roaming and substantially higher H2 vibrational excitation.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Perspective: chemical dynamics simulations of non-statistical reaction dynamics

Non-statistical chemical dynamics are exemplified by disagreements with the transition state (TS), RRKM and phase space theories of chemical kinetics and dynamics. The intrinsic reaction coordinate (IRC) is often used for the former two theories, and non-statistical dynamics arising from non-IRC dynamics are often important. In this perspective, non-statistical dynamics are discussed for chemical reactions, with results primarily obtained from chemical dynamics simulations and to a lesser extent from experiment. The non-statistical dynamical properties discussed are: post-TS dynamics, including potential energy surface bifurcations, product energy partitioning in unimolecular dissociation and avoiding exit-channel potential energy minima; non-RRKM unimolecular decomposition; non-IRC dynamics; direct mechanisms for bimolecular reactions with pre- and/or post-reaction potential energy minima; non-TS theory barrier recrossings; and roaming dynamics.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.