Reduced dimensionality quantum dynamics of Cl + CH4? HCl + CH3 on an ab initio potential

Isotope Effect and Heavy-Light-Heavy Reactivity Oscillation in the Cl + CHD3/CHT3 Reaction

Previous experiments and theories have shown the existence of heavy-light-heavy (HLH) reactivity oscillation in the Cl + CH4 reaction and anticipated that similar oscillations should exist in many HLH reactions involving polyatomic reagents. However, the total reaction probabilities for the Cl + CHD3 → HCl + CD3 reaction exhibit only a step-like feature, and the total reaction probabilities for Cl + CHT3 → HCl + CT3 do not show any structure at all. Here, we report seven-dimensional state-to-state quantum dynamics studies for this reaction on the FI-NN PES, and we demonstrate that HLH reactivity oscillations also exist in these two reactions, manifesting as peaks in the reaction probabilities for low product rotational states. These oscillations, however, are obscured in the total reaction probability because of the higher excitation of j ≥ 2 product rotational states. Furthermore, the isotope replacement of nonreactive hydrogen with deuterium and tritium significantly enhances reactivity at collision energies above 0.112 eV, indicating an inverse secondary isotope effect on the probabilities, which is proved to be also caused by HLH mass combination. We also demonstrate that the highly rotational excitation of CHD3 substantially enhances reactivity and the HLH oscillations, similar to HLH triatomic reactions. These observations are completely different from those in the H + CHD3 reaction, which is also a late-barrier reaction. Therefore, the HLH mass combination is very important, which affects not only the reactivity oscillation but also the amplitude and product rotational state distribution and makes the initial rotation excitation play a pivotal role in the reaction.

ManyHF-based full-dimensional potential energy surface development and quasi-classical dynamics for the Cl + CH3NH2 reaction

A full-dimensional spin-orbit (SO)-corrected potential energy surface (PES) is developed for the Cl + CH3NH2 multi-channel system. Using the new PES, a comprehensive reaction dynamics investigation is performed for the most reactive hydrogen-abstraction reactions forming HCl + CH2NH2/CH3NH. Hartree-Fock (HF) convergence problems in the reactant region are handled by the ManyHF method, which finds the lowest-energy HF solution considering several different initial guess orbitals. The PES development is carried out with the Robosurfer program package, which iteratively improves the surface. Energy points are computed at the ManyHF-UCCSD(T)-F12a/cc-pVDZ-F12 level of theory combined with basis set (ManyHF-RMP2-F12/cc-pVTZ-F12 - ManyHF-RMP2-F12/cc-pVDZ-F12) and SO (MRCI+Q/aug-cc-pwCVDZ) corrections. Quasi-classical trajectory simulations show that the CH3-side hydrogen abstraction occurs more frequently in contrast to the NH2-side reaction. In both cases, the integral cross sections decrease with increasing collision energy (Ecoll). A reaction mechanism shifting from indirect to direct stripping can be observed from the opacity functions, scattering angle, and translation energy distributions as Ecoll increases. Initial attack angle distributions reveal that chlorine prefers to abstract hydrogen from the approached functional group. The collision-energy dependence of the product energy distributions shows that the initial translational energy mainly transfers to product recoil. The HCl vibrational and rotational energy values are comparable and nearly independent of collision energy, while the CH2NH2 and CH3NH co-products' vibrational energy values are higher than the rotational energy values with more significant Ecoll dependence. The HCl(v = 0) rotational distributions are compared with experiment, setting the direction for future investigations.

First-Principles Theory for the Reaction of Chlorine with Methane

A full-dimensional quantum dynamics simulation of the Cl + CH₄ → HCl + CH₃ reaction based on first-principles theory is reported. Accurate thermal rate constants are calculated, and perfect agreement with experiment is obtained. Despite the heavy atoms present in both reactants, the passage of the reaction barrier is found to occur within only a few tens of femtoseconds. This surprisingly short time scale results from correlated motion of the transferring hydrogen atom and the hydrogen atoms in the methyl fragment which facilitates irreversible barrier passage without relevant participation of heavy atoms. Resonance effects resulting from the heavy-light-heavy characteristics of the reaction system, which were observed in reactive scattering studies, do not affect the thermal rate constant.

Reactivity oscillation in the heavy-light-heavy Cl + CH4 reaction

It has long been predicted that oscillatory behavior exists in reactivity as a function of collision energy for heavy-light-heavy (HLH) chemical reactions in which a light atom is transferred between two heavy atoms or groups of atoms, but direct observation of such a behavior in bimolecular reactions remains a challenge. Here we report a joint theoretical and crossed-molecular-beam study on the Cl + CH₄ → HCl + CH₃ reaction. A distinctive peak at a collision energy of 0.15 eV for the CH₃(v = 0) product was experimentally detected in the backward scattering direction. Detailed quantum-dynamics calculations on a highly accurate potential energy surface revealed that this feature originates from the reactivity oscillation in this HLH polyatomic reaction. We anticipate that such reactivity oscillations exist in many HLH reactions involving polyatomic reagents.

Benchmark ab initio and dynamical characterization of the stationary points of reactive atom + alkane and SN2 potential energy surfaces

We review composite ab initio and dynamical methods and their applications to characterize stationary points of atom/ion + molecule reactions.

A Machine Learning Approach for Prediction of Rate Constants

We report a machine learning approach to train and predict bimolecular thermal rate constants over a large temperature range. The approach uses Gaussian process (GP) regression to evaluate the difference between accurate quantum results and Eckart-corrected conventional transition state theory, mostly for collinear reactions. Training is done on a database of rate constants for 13 reaction/potential surface combinations, and testing is performed on a set of 39 reaction/potential surface combinations. Averaged over all test reactions, the GP method is within 80% of the accurate answer, whereas transition state theory (TST) is only within 330% and Eckart-corrected TST (ECK) is within 110%. In the tunneling region, GP is generally (with a few exceptions) more accurate and sometimes much more accurate. In the high-temperature recrossing region, GP is significantly more accurate than either TST or ECK, neither of which addresses the possibility of recrossing. The GP predictions for the 3D reactions O(³P) + H₂, OH + H₂, O(³P) + CH₄, and H + CH₄, for which accurate quantum results are available, provide further encouragement to the machine learning approach.

Molecular Dynamics Simulations on Relaxed Reduced-Dimensional Potential Energy Surfaces

Molecular dynamics (MD) simulations with full-dimensional potential energy surfaces (PESs) obtained from high-level ab initio calculations are frequently used to model reaction dynamics of small molecules (i.e., molecules with up to 10 atoms). Construction of full-dimensional PESs for larger molecules is, however, not feasible since the number of ab initio calculations required grows rapidly with the increase of dimension. Only a small number of coordinates are often essential for describing the reactivity of even very large systems, and reduced-dimensional PESs with these coordinates can be built for reaction dynamics studies. While analytical methods based on transition-state theory framework are well established for analyzing the reduced-dimensional PESs, MD simulation algorithms capable of generating trajectories on such surfaces are more rare. In this work, we present a new MD implementation that utilizes the relaxed reduced-dimensional PES for standard microcanonical (NVE) and canonical (NVT) MD simulations. The method is applied to the pyramidal inversion of a NH₃ molecule. The results from the MD simulations on a reduced, three-dimensional PES are validated against the ab initio MD simulations, as well as MD simulations on full-dimensional PES and experimental data.

Spiers Memorial Lecture : Introductory lecture: quantum dynamics of chemical reactions

This Spiers Memorial Lecture discusses quantum effects that can be calculated and observed in the chemical reactions of small molecules.

Kinetic isotope effects and how to describe them

We review several methods for computing kinetic isotope effects in chemical reactions including semiclassical and quantum instanton theory. These methods describe both the quantization of vibrational modes as well as tunneling and are applied to the ⋅H + H2 and ⋅H + CH4 reactions. The absolute rate constants computed with the semiclassical instanton method both using on-the-fly electronic structure calculations and fitted potential-energy surfaces are also compared directly with exact quantum dynamics results. The error inherent in the instanton approximation is found to be relatively small and similar in magnitude to that introduced by using fitted surfaces. The kinetic isotope effect computed by the quantum instanton is even more accurate, and although it is computationally more expensive, the efficiency can be improved by path-integral acceleration techniques. We also test a simple approach for designing potential-energy surfaces for the example of proton transfer in malonaldehyde. The tunneling splittings are computed, and although they are found to deviate from experimental results, the ratio of the splitting to that of an isotopically substituted form is in much better agreement. We discuss the strengths and limitations of the potential-energy surface and based on our findings suggest ways in which it can be improved.

Angular Scattering Dynamics of the CH4 + Cl → CH3 + HCl Reaction Using Nearside-Farside, Local Angular Momentum, and Resummation Theories

The differential cross section (DCS) for the CH4 + Cl → CH3 + HCl reaction is studied at six total energies where all of the species are in their ground states. The scattering (S) matrix elements have been calculated by the rotating line umbrella method for a dual-level ab initio analytic potential energy surface. We make the first application to this reaction of nearside-farside (NF) and local angular momentum (LAM) techniques, including resummation orders (r) of 0, 1, 2, and 3 for the partial-wave series representation of the full scattering amplitude. We find that resummation usually cleans the NF r = 0 DCSs of unphysical oscillations, especially at small angles. This cleaning effect is typically most pronounced when changing from no resummation (r = 0) to r = 1; further resummations from r = 1 to r = 2 and from r = 2 to r = 3 have smaller effects. The NF DCS analyses show that the reaction is N-dominated at sideward and large angles, whereas at small angles there are oscillations caused by NF interference. The NF LAM analysis provides consistent and complementary information, in particular for the total angular momenta that contribute to the reaction at different scattering angles. The NF analyses also provide justification for simpler N-dominant dynamical theories such as the semiclassical optical model, which provides an explanation for the distorted mirror image effect for the moduli of the S matrix elements and the DCSs, as well as the use of a hard-sphere DCS over limited angular ranges.

Rate constant calculations of the GeH4 + OH/OD → GeH3 + H2O/HOD reactions using an ab initio based full-dimensional potential energy surface

2D representation of the analytical potential energy surface. The saddle point and the complexes in the entry and exit channels are included.

The hydrogen abstraction reaction O(3P) + CH4: a new analytical potential energy surface based on fit to ab initio calculations

Based exclusively on high-level ab initio calculations, a new full-dimensional analytical potential energy surface (PES-2014) for the gas-phase reaction of hydrogen abstraction from methane by an oxygen atom is developed. The ab initio information employed in the fit includes properties (equilibrium geometries, relative energies, and vibrational frequencies) of the reactants, products, saddle point, points on the reaction path, and points on the reaction swath, taking especial caution respecting the location and characterization of the intermediate complexes in the entrance and exit channels. By comparing with the reference results we show that the resulting PES-2014 reproduces reasonably well the whole set of ab initio data used in the fitting, obtained at the CCSD(T) = FULL/aug-cc-pVQZ//CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical surface we perform an extensive dynamics study using quasi-classical trajectory calculations, comparing the results with recent experimental and theoretical data. The excitation function increases with energy (concave-up) reproducing experimental and theoretical information, although our values are somewhat larger. The OH rotovibrational distribution is cold in agreement with experiment. Finally, our results reproduce experimental backward scattering distribution, associated to a rebound mechanism. These results lend confidence to the accuracy of the new surface, which substantially improves the results obtained with our previous surface (PES-2000) for the same system.

Ab initio based potential energy surface and kinetics study of the OH + NH3 hydrogen abstraction reaction

A full-dimensional analytical potential energy surface (PES) for the OH + NH3 → H2O + NH2 gas-phase reaction was developed based exclusively on high-level ab initio calculations. This reaction presents a very complicated shape with wells along the reaction path. Using a wide spectrum of properties of the reactive system (equilibrium geometries, vibrational frequencies, and relative energies of the stationary points, topology of the reaction path, and points on the reaction swath) as reference, the resulting analytical PES reproduces reasonably well the input ab initio information obtained at the coupled-cluster single double triple (CCSD(T)) = FULL/aug-cc-pVTZ//CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical PES we perform an extensive kinetics study using variational transition-state theory with semiclassical transmission coefficients over a wide temperature range, 200-2000 K. The forward rate constants reproduce the experimental measurements, while the reverse ones are slightly underestimated. However, the detailed analysis of the experimental equilibrium constants (from which the reverse rate constants are obtained) permits us to conclude that the experimental reverse rate constants must be re-evaluated. Another severe test of the new surface is the analysis of the kinetic isotope effects (KIEs), which were not included in the fitting procedure. The KIEs reproduce the values obtained from ab initio calculations in the common temperature range, although unfortunately no experimental information is available for comparison.

Quantum effects in the abstraction reaction by H atoms of primary and secondary hydrogens in n-C4H10: a test of a new potential energy surface construction method

Recently, von Horsten et al. suggested an efficient method to construct two dimensional potential energy surfaces (PESs) for use in quantum scattering simulations, in which they utilised the minimum energy path (MEP) and assumed harmonic behaviour of the PES near the MEP. In the same paper, the authors applied this method to various H-abstraction reactions from C1-C3 alkane molecules. In this work we demonstrate an alternative PES construction method, and apply it to the more challenging H-abstraction from n-C(4)H(10) reactions. The geometry optimizations and frequency calculations are done at the MP2/cc-pVTZ level of theory, while the energies are calculated with the CCSD(T) method with the same basis set. The calculations give adiabatic energy barrier heights of 45.9 kJ mol(-1) and 34.4 kJ mol(-1) for the primary and secondary hydrogens in n-C(4)H(10). When compared to purely classical transition state theory, quantum scattering calculations show that quantum tunnelling and zero-point effects have large contributions at low temperatures, typically below 500 K. The branching ratio study suggests that the abstraction of secondary hydrogen in n-C(4)H(10) dominates the overall reaction rate at low temperatures. The rate constants for the two abstraction channels become more comparable as the temperature increases.

Constructing Potential Energy Surfaces for Polyatomic Systems: Recent Progress and New Problems

Different methods of constructing potential energy surfaces in polyatomic systems are reviewed, with the emphasis put on fitting, interpolation, and analytical (defined by functional forms) approaches, based on quantum chemistry electronic structure calculations. The different approaches are reviewed first, followed by a comparison using the benchmark H + CH4 and the H + NH3 gas-phase hydrogen abstraction reactions. Different kinetics and dynamics properties are analyzed for these reactions and compared with the available experimental data, which permits one to estimate the advantages and disadvantages of each method. Finally, we analyze different problems with increasing difficulty in the potential energy construction: spin-orbit coupling, molecular size, and more complicated reactions with several maxima and minima, which test the soundness and general applicability of each method. We conclude that, although the field of small systems, typically atom-diatom, is mature, there still remains much work to be done in the field of polyatomic systems.

The hydrogen abstraction reaction H + CH4. I. New analytical potential energy surface based on fitting to ab initio calculations

A new analytical potential energy surface is presented for the reaction of hydrogen abstraction from methane by a hydrogen atom. It is based on an analytical expression proposed by Jordan and Gilbert [J. Chem. Phys. 102, 5669 (1995)], and its fittable parameters were obtained by a multibeginning optimization procedure to reproduce high-level ab initio electronic structure calculations obtained at the CCSD(T)/cc-pVTZ level. The ab initio information employed in the fit includes properties (equilibrium geometries, relative energies, and vibrational frequencies) of the reactants, products, saddle point, points on the reaction path, and points on the reaction swath. No experimental information is used. By comparison with the reference results we show that the resulting surface reproduces well not only the ab initio data used in the fitting but also other thermochemical and kinetic results computed at the same ab initio level, such as equilibrium constants, rate constants, and kinetic isotope effects, which were not used in the fit. In this way we show that the new potential energy surface is correctly fitted and almost as accurate as the CCSD(T)/cc-pVTZ method in describing the kinetics of the reaction. We analyze the limitations of the functional form and the fitting method employed, and suggest some solutions to their drawbacks. In a forthcoming communication, we test the quality of the new surface by comparing its results with experimental values.

Chemical dynamics of vibrationally excited molecules: Controlling reactions in gases and on surfaces

Experimental studies of the chemical reaction dynamics of vibrationally excited molecules reveal the ability of different vibrations to control the course of a reaction. This Perspective describes those studies for the prototypical reaction of vibrationally excited methane and its isotopologues in gases and on surfaces and looks to the prospects of similar studies in liquids. The influences of vibrational excitation on the C-H bond cleavage in a single collision reaction with Cl and in dissociative adsorption on a Ni surface bear some striking similarities. Both reactions are bond-selective processes in which the initial preparation of a molecular eigenstate containing a large component of C-H stretching results in preferential cleavage of that bond. It is possible to cleave either the C-H bond or C-D bond in the reaction of Cl with CH3D, CH2D2, or CHD3 and, similarly, to use initial excitation of the C-H stretch to promote dissociation of CHD3 to CD3 and H on a Ni surface. Different vibrational modes, such as the symmetric and antisymmetric stretches in CH3D or CH4, lead to very different reactivities, and molecules with the symmetric stretching vibration excited can be as much as 10 times more reactive than ones with the antisymmetric stretch excited. The origin of this behavior lies in the change in the vibrational motion induced by the interaction with the atomic reaction partner or the surface.