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

Three-state surface hopping calculations of acetaldehyde photodissociation to CH3 + HCO on ab initio potential surfaces

We report Trajectory Surface Hopping (TSH) calculations of CH3CHO photodissociation involving three electronic states, S1, T1, and S0, with a focus on the radical products CH3 + HCO, which can be formed from both T0 and S0. We use previously reported potential energy surfaces and spin-orbit couplings for T1 and S0 and report a new potential and spin-orbit coupling for S1 here. Roughly 32 000 trajectories are performed at energies corresponding to seven photolysis wavelengths between 372 and 230 nm. Motivated by recent experiments, we examine the branching ratio of the T1 to S0 pathways as a function of photolysis energy. We also present the relative translational energy and CH3 vibrational energy distributions from these pathways at a photolysis energy of 100 kcal mol(-1), formed from both the T1 and S0 potentials. As with standard quasiclassical trajectory calculations, violation of zero-point energy for products also occurs in TSH calculations. This is shown to be a serious issue for this branching ratio and one of several methods considered to deal with this issue is shown to give satisfactory results.

Electronic Population Inversion in HCCO/DCCO Products from Hyperthermal Collisions of O((3)P) with HCCH/DCCD

The dynamics of hyperthermal O((3)P) reactions with acetylene have been investigated with the use of crossed molecular beams techniques, employing both mass spectrometric and optical detection of products. With collision energies of 40-150 kcal mol(-1), O((3)P) + HCCH/DCCD → HCCO/DCCO + H/D may follow multiple pathways to form the ketenyl radical (HCCO or DCCO) in ground doublet states or in electronically excited quartet and doublet states. Theoretical calculations support the assignment of the various reaction pathways. The fraction of electronic excitation is substantial. At the highest collision energy studied, ∼65% of the ketenyl radical products that survive are electronically excited, with the majority of the excited products in a quartet state. In this case, a population inversion exists between the electronically excited quartet and ground doublet states of the ketenyl product. Such significant electronic excitation in products is unusual in bimolecular reactions, especially when ground-state products are accessible by spin-allowed pathways.

Indirect dynamics in a highly exoergic substitution reaction

The highly exoergic nucleophilic substitution reaction F(-) + CH3I shows reaction dynamics strikingly different from that of substitution reactions of larger halogen anions. Over a wide range of collision energies, a large fraction of indirect scattering via a long-lived hydrogen-bonded complex is found both in crossed-beam imaging experiments and in direct chemical dynamics simulations. Our measured differential scattering cross sections show large-angle scattering and low product velocities for all collision energies, resulting from efficient transfer of the collision energy to internal energy of the CH3F reaction product. Both findings are in strong contrast to the previously studied substitution reaction of Cl(-) + CH3I [Science 2008, 319, 183-186] at all but the lowest collision energies, a discrepancy that was not captured in a subsequent study at only a low collision energy [J. Phys. Chem. Lett. 2010, 1, 2747-2752]. Our direct chemical dynamics simulations at the DFT/B97-1 level of theory show that the reaction is dominated by three atomic-level mechanisms, an indirect reaction proceeding via an F(-)-HCH2I hydrogen-bonded complex, a direct rebound, and a direct stripping reaction. The indirect mechanism is found to contribute about one-half of the overall substitution reaction rate at both low and high collision energies. This large fraction of indirect scattering at high collision energy is particularly surprising, because the barrier for the F(-)-HCH2I complex to form products is only 0.10 eV. Overall, experiment and simulation agree very favorably in both the scattering angle and the product internal energy distributions.

Quasiclassical trajectory study of CH3NO2 decomposition via roaming mediated isomerization using a global potential energy surface

We report a global potential energy surface (PES) for CH3NO2 based on fitting roughly 114,000 density functional theory (UB3LYP/6-311+g(d,p)) electronic energies. The PES is a precise, permutationally invariant fit to these energies. Properties of the PES are described, as are some preliminary quasiclassical trajectory calculations. An isomerization-roaming pathway to the CH3ONO isomer and then to the CH3O + NO products is found. Although the pathway occurs at larger distances than a related loose saddle-point on the PES, the pathway supports the supposition of such a pathway based on locating a loose first-order saddle point and associated IRC, reported previously by Zhu and Lin [Zhu, R. S. and Lin, M. C. Chem. Phys. Lett. 2009, 478, 11].

Exploring Multiple Potential Energy Surfaces: Photochemistry of Small Carbonyl Compounds

In theoretical studies of chemical reactions involving multiple potential energy surfaces (PESs) such as photochemical reactions, seams of intersection among the PESs often complicate the analysis. In this paper, we review our recipe for exploring multiple PESs by using an automated reaction path search method which has previously been applied to single PESs. Although any such methods for single PESs can be employed in the recipe, the global reaction route mapping (GRRM) method was employed in this study. By combining GRRM with the proposed recipe, all critical regions, that is, transition states, conical intersections, intersection seams, and local minima, associated with multiple PESs, can be explored automatically. As illustrative examples, applications to photochemistry of formaldehyde and acetone are described. In these examples as well as in recent applications to other systems, the present approach led to discovery of many unexpected nonadiabatic pathways, by which some complicated experimental data have been explained very clearly.

Bromenium-catalysed tandem ring opening/cyclisation of vinylcyclopropanes and vinylcyclobutanes: a metal-free [3+2+1]/[4+2+1] cascade for the synthesis of chiral amidines and computational investigation

We present a detailed study of a [3+2+1] cascade cyclisation of vinylcyclopropanes (VCP) catalysed by a bromenium species (Br(δ+)-X(δ-)) generated in situ, which results in the synthesis of chiral bicyclic amidines in a tandem one-pot operation. The formation of amidines involves the ring-opening of VCPs with Br-X, followed by a Ritter-type reaction with chloramine-T and a tandem cyclisation. The reaction has been further extended to vinylcyclobutane systems and involves a [4+2+1] cascade cyclisation with the same reagents. The versatility of the methodology has been demonstrated by careful choice of VCPs and VCBs to yield bicyclo[4.3.0]-, -[4.3.1]- and -[4.4.0]amidines in enantiomerically pure form. On the basis of the experimental observations and DFT calculations, a reasonable mechanism has been put forth to account for the formation of the products and the observed stereoselectivity. We propose the existence of a π-stabilised homoallylic carbocation at the cyclopropane carbon as the reason for high stereoselectivity. DFT studies at B3LYP/6-311+G** and M06-2X/6-31+G* levels of theory in gas-phase calculations suggest the ring-opening of VCP is initiated at the π-complex stage (between the double bond and Br-X). This can be clearly perceived from the solution-phase (acetonitrile) calculations using the polarisable continuum model (PCM) solvation model, from which the extent of the ring opening of VCP was found to be noticeably high. Studies also show that the formation of zero-bridge bicyclic amidines is favoured over other bridged bicyclic amidines. The energetics of competing reaction pathways is compared to explain the product selectivity.

Automated Exploration of Photolytic Channels of HCOOH: Conformational Memory via Excited-State Roaming

To elucidate the photodissociation mechanism of HCOOH, we systematically explored reaction pathways starting from the first excited singlet state (S1) by using automated reaction path search methods. All critical points, that is, minima, transition states, minimum energy conical intersections, and minima on seam of crossing, for the S0, T1, and S1 potential energy surfaces (PESs) obtained in the present search were optimized at the CASPT2 level. The structure list obtained by the search explained all experimentally reported photolytic channels. A new mechanism for the previously suggested but unexplained conformational memory in the 193 nm photolysis is proposed, which involves two steps: partial dissociation and succeeding roaming of one of H atoms on the S1 PES, followed by intramolecular recombination on the S0 PES after radiationless transition through a conical intersection at a partially dissociated geometry. This is partially similar to the excited-state roaming recently discovered for the NO3 radical.

Bi-fidelity fitting and optimization

A common feature in computations of chemical and physical properties is the investigation of phenomena at different levels of computational accuracy. Less accurate computations are used to provide a relatively quick understanding of the behavior of a system and allow a researcher to focus on regions of initial conditions and parameter space where interesting phenomena are likely to occur. These inexpensive calculations are often discarded when more accurate calculations are performed. This paper demonstrates how computations at different levels of accuracy can be simultaneously incorporated to study chemical and physical phenomena with less overall computational effort than the most expensive level of computation. A smaller set of computationally expensive calculations is needed because the set of expensive calculations is correlated with the larger set of less expensive calculations. We present two applications. First, we demonstrate how potential energy surfaces can be fit by simultaneously using results from two different levels of accuracy in electronic structure calculations. In the second application, we study the optical response of metallic nanostructures. The optical response is generated with calculations at two different grid resolutions, and we demonstrate how using these two levels of computation in a correlated fashion can more efficiently optimize the response.

Dynamics, transition states, and timing of bond formation in Diels-Alder reactions

The time-resolved mechanisms for eight Diels-Alder reactions have been studied by quasiclassical trajectories at 298 K, with energies and derivatives computed by UB3LYP/6-31G(d). Three of these reactions were also simulated at high temperature to compare with experimental results. The reaction trajectories require 50-150 fs on average to transverse the region near the saddle point where bonding changes occur. Even with symmetrical reactants, the trajectories invariably involve unequal bond formation in the transition state. Nevertheless, the time gap between formation of the two new bonds is shorter than a C ─ C vibrational period. At 298 K, most Diels-Alder reactions are concerted and stereospecific, but at high temperatures (approximately 1,000 K) a small fraction of trajectories lead to diradicals. The simulations illustrate and affirm the bottleneck property of the transition state and the close connection between dynamics and the conventional analysis based on saddle point structure.

Separability of tight and roaming pathways to molecular decomposition

Recent studies have questioned the separability of the tight and roaming mechanisms to molecular decomposition. We explore this issue for a variety of reactions including MgH(2) → Mg + H(2), NCN → CNN, H(2)CO → H(2) + CO, CH(3)CHO → CH(4) + CO, and HNNOH → N(2) + H(2)O. Our analysis focuses on the role of second-order saddle points in defining global dividing surfaces that encompass both tight and roaming first-order saddle points. The second-order saddle points define an energetic criterion for separability of the two mechanisms. Furthermore, plots of the differential contribution to the reactive flux along paths connecting the first- and second-order saddle points provide a dynamic criterion for separability. The minimum in the differential reactive flux in the neighborhood of the second-order saddle point plays the role of a mechanism divider, with the presence of a strong minimum indicating that the roaming and tight mechanisms are dynamically distinct. We show that the mechanism divider is often, but not always, associated with a second-order saddle point. For the formaldehyde and acetaldehyde reactions, we find that the minimum energy geometry on a conical intersection is associated with the mechanism divider for the tight and roaming processes. For HNNOH, we again find that the roaming and tight processes are dynamically separable but we find no intrinsic feature of the potential energy surface associated with the mechanism divider. Overall, our calculations suggest that roaming and tight mechanisms are generally separable over broad ranges of energy covering most kinetically relevant regimes.

Observation of NH X(3)Σ(-) as a Primary Product of Methylamine Photodissociation: Evidence of Roaming-Mediated Intersystem Crossing?

3+1 Resonance-enhanced multiphoton ionization and photofragment excitation spectroscopy have been used to identify NH X(3)Σ(-) as a primary product of methylamine photodissociation after state-specific excitation to the S1 state. On the basis of standard thermochemical data, NH X(3)Σ(-) can be formed only in conjunction with closed-shell CH4 coproducts, indicating that dissociation must occur on the T1 surface. It is proposed that the mechanism for the formation of triplet NH and CH4 involves intramolecular abstraction between frustrated radical products and is an example of roaming-mediated intersystem crossing.

Imaging the molecular channel in acetaldehyde photodissociation: roaming and transition state mechanisms

The roaming dynamics in the photodissociation of acetaldehyde is studied through the first absorption band, in the wavelength interval ranging from 230 nm to 325 nm. Using a combination of the velocity-map imaging technique and rotational resonance enhanced multiphoton ionization (REMPI) spectroscopy of the CO fragment, the branching ratio between the canonical transition state and roaming dissociation mechanisms is obtained at each of the photolysis wavelengths studied. Upon one photon absorption, the molecule is excited to the first singlet excited S(1) state, which, depending on the excitation wavelength, either converts back to highly vibrationally excited ground S(0) state or undergoes intersystem crossing to the first excited triplet T(1) state, from where the molecule can dissociate over two main channels: the radical (CH(3) + HCO) and the molecular (CO + CH(4)) channels. Three dynamical regions are characterized: in the red edge of the absorption band, at excitation energies below the T(1) barrier, the ratio of the roaming dissociation channel increases, largely surpassing the transition state contribution. As the excitation wavelength is increased, the roaming propensity decreases reaching a minimum at wavelengths ∼308 nm. Towards the blue edge, at 230 nm, an upper limit of ∼50% has been estimated for the contribution of the roaming channel. The experimental results are interpreted in terms of the interaction between the different potential energy surfaces involved by means of ab initio stationary points and intrinsic reaction coordinate paths calculations.

Shock tube explorations of roaming radical mechanisms: the decompositions of isobutane and neopentane

The thermal decompositions of isobutane and neopentane have been studied using both shock tube experiments and ab initio transition state theory based master equation calculations. Dissociation rate constants for these molecules have been measured at high temperatures (1260-1566 K) behind reflected shock waves using high-sensitivity H-ARAS detection. The two major dissociation channels at high temperature are iso-C(4)H(10) → CH(3) + i-C(3)H(7) (1a) and neo-C(5)H(12) → CH(3) + t-C(4)H(9) (2a). Ultrahigh-sensitivity ARAS detection of H-atoms produced from the rapid decomposition of the product radicals, i-C(3)H(7) in (1a) and t-C(4)H(9) in (2a), through i-C(3)H(7) + M → H + C(3)H(6) + M (3a) and t-C(4)H(9) + M → H + i-C(4)H(8) + M (4a) allowed measurements of both the total decomposition rate constants, k(total), and the branching to radical products, which were observed to be equivalent in both systems, k(1a)/k(total) and k(2a)/k(total) = 0.79 ± 0.05. Theoretical analyses indicate that in isobutane, the non-H-atom fraction has two contributions, the dominant fraction being due to the roaming radical mechanism leading to molecular products through iso-C(4)H(10) → CH(4) + C(3)H(6) (1b) with k(1b)/k(total) = 0.16, and a minor fraction that involves the isomerization of i-C(3)H(7) to n-C(3)H(7) that then subsequently forms methyl radicals, i-C(3)H(7) + M → n-C(3)H(7) + M → CH(3) + C(2)H(4) + M (3b). In contrast to isobutane, in neopentane, the contribution to the non-H-atom fraction is exclusively through the roaming radical mechanism that leads to neo-C(5)H(12) → CH(4) + i-C(4)H(8) (2b) with k(2b)/k(total) = 0.21. These quantitative measurements of larger contributions from the roaming mechanism for larger molecules are in agreement with the qualitative theoretical arguments that suggest long-range dispersion interactions (which become increasingly important for larger molecules) may enhance roaming.

Developments in Laboratory Studies of Gas-Phase Reactions for Atmospheric Chemistry with Applications to Isoprene Oxidation and Carbonyl Chemistry

Laboratory studies of gas-phase chemical processes are a key tool in understanding the chemistry of our atmosphere and hence tackling issues such as climate change and air quality. Laboratory techniques have improved considerably with greater emphasis on product detection, allowing the measurement of site-specific rate coefficients. Radical chemistry lies at the heart of atmospheric chemistry. In this review we consider issues around radical generation and recycling from the oxidation of isoprene and from the chemical reactions and photolysis of carbonyl species. Isoprene is the most globally significant hydrocarbon, but uncertainties exist about its oxidation in unpolluted environments. Recent experiments and calculations that cast light on radical generation are reviewed. Carbonyl compounds are the dominant first-generation products from hydrocarbon oxidation. Chemical oxidation can recycle radicals, or photolysis can be a net radical source. Studies have demonstrated that high-resolution and temperature-dependent studies are important for some significant species.