A global potential energy surface for the H2 + OH ↔ H2O + H reaction using neural networks

Recent Advances in Quantum Dynamics of Bimolecular Reactions

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

Methane dissociation on Ni(111): A fifteen-dimensional potential energy surface using neural network method

In the present work, we develop a highly accurate, fifteen-dimensional potential energy surface (PES) of CH4 interacting on a rigid flat Ni(111) surface with the methodology of neural network (NN) fit to a database consisted of about 194 208 ab initio density functional theory (DFT) energy points. Some careful tests of the accuracy of the fitting PES are given through the descriptions of the fitting quality, vibrational spectrum of CH4 in vacuum, transition state (TS) geometries as well as the activation barriers. Using a 25-60-60-1 NN structure, we obtain one of the best PESs with the least root mean square errors: 10.11 meV for the entrance region and 17.00 meV for the interaction and product regions. Our PES can reproduce the DFT results very well in particular for the important TS structures. Furthermore, we present the sticking probability S0 of ground state CH4 at the experimental surface temperature using some sudden approximations by Jackson's group. An in-depth explanation is given for the underestimated sticking probability.

Mode specific dynamics in the H2 + SH → H + H2S reaction

Full-dimensional quantum dynamics and quasi-classical trajectory studies indicate strong mode selectivity in the H₂ + SH reaction.

An accurate multi-channel multi-reference full-dimensional global potential energy surface for the lowest triplet state of H2O2

A full dimensional potential energy surface for the lowest triplet state of H₂O₂ was developed at the MRCI-F12 level.

A local mode picture for H atom reaction with vibrationally excited H2O: a full dimensional state-to-state quantum dynamics investigation

A local mode picture was postulated some time ago to explain the negligible OD product excitation from H reacting with D₂O in fundamental asymmetrical excitation: the D₂O asymmetric stretch state is best thought as a linear combination of local mode stretches with the H atom reacting preferentially with the excited OD bond, but has never been verified. Here we report the first full-dimensional state-to-state study for the title reaction with H₂O in the ground and the first symmetric and asymmetric stretching excited states. It is found that symmetric and asymmetric stretching excited states, which vibrate in very different ways in the normal mode picture, behave essentially identically in the reaction. More importantly, our calculations revealed that the reaction produces a small fraction of OH in the v = 1 state, with the population close to the relative reactivity between the ground and vibrationally excited states, and therefore it confirms the local mode picture for H₂O in symmetric or asymmetric stretch states when reacting with an H atom.

Dynamics of water dissociative chemisorption on Ni(111): effects of impact sites and incident angles

The dissociative chemisorption of water on rigid Ni(111) is investigated using a quasiclassical trajectory method on a nine-dimensional global potential energy surface based on a faithful permutation invariant fit of ∼25 000 density functional theory points. This full-dimensional model not only confirms the validity of our earlier reduced-dimensional model with 6 degrees of freedom, but also allows the examination of the influence of impact sites and incident angles. It is shown that the reactivity depends on the site of impact in a complex fashion controlled by the topography of the potential energy surface rather than the barrier height alone. In addition, the reaction is promoted by momenta both parallel and perpendicular to the surface, as predicted by the recently proposed sudden vector projection model.

CH4 dissociation on Ni(111): a quantum dynamics study of lattice thermal motion

Lattice thermal motion is of great importance because it has a significant effect on molecule activation on metal surfaces.

Extremely short-lived reaction resonances in Cl + HD (v = 1) → DCl + H due to chemical bond softening

The Cl + H2 reaction is an important benchmark system in the study of chemical reaction dynamics that has always appeared to proceed via a direct abstraction mechanism, with no clear signature of reaction resonances. Here we report a high-resolution crossed–molecular beam study on the Cl + HD (v = 1, j = 0) → DCl + H reaction (where v is the vibrational quantum number and j is the rotational quantum number). Very few forward scattered products were observed. However, two distinctive peaks at collision energies of 2.4 and 4.3 kilocalories per mole for the DCl (v′ = 1) product were detected in the backward scattering direction. Detailed quantum dynamics calculations on a highly accurate potential energy surface suggested that these features originate from two very short-lived dynamical resonances trapped in the peculiar H-DCl (v′ = 2) vibrational adiabatic potential wells that result from chemical bond softening. We anticipate that dynamical resonances trapped in such wells exist in many reactions involving vibrationally excited molecules.

Representing potential energy surfaces by high-dimensional neural network potentials

The development of interatomic potentials employing artificial neural networks has seen tremendous progress in recent years. While until recently the applicability of neural network potentials (NNPs) has been restricted to low-dimensional systems, this limitation has now been overcome and high-dimensional NNPs can be used in large-scale molecular dynamics simulations of thousands of atoms. NNPs are constructed by adjusting a set of parameters using data from electronic structure calculations, and in many cases energies and forces can be obtained with very high accuracy. Therefore, NNP-based simulation results are often very close to those gained by a direct application of first-principles methods. In this review, the basic methodology of high-dimensional NNPs will be presented with a special focus on the scope and the remaining limitations of this approach. The development of NNPs requires substantial computational effort as typically thousands of reference calculations are required. Still, if the problem to be studied involves very large systems or long simulation times this overhead is regained quickly. Further, the method is still limited to systems containing about three or four chemical elements due to the rapidly increasing complexity of the configuration space, although many atoms of each species can be present. Due to the ability of NNPs to describe even extremely complex atomic configurations with excellent accuracy irrespective of the nature of the atomic interactions, they represent a general and therefore widely applicable technique, e.g. for addressing problems in materials science, for investigating properties of interfaces, and for studying solvation processes.

A nine-dimensional global potential energy surface for NH4(X2A1) and kinetics studies on the H + NH3 ↔ H2 + NH2 reaction

A nine-dimensional global potential energy surface (PES) for the NH₄ system is developed from ∼10⁵ high-level ab initio points and the hydrogen abstraction kinetics on the PES agree with experiment.

SEVEN-DIMENSIONAL QUANTUM DYNAMICS STUDY OF THE H2 + NH2 → H + NH3 REACTION ON AN INTERPOLATED POTENTIAL ENERGY SURFACE

Initial-state-selected time-dependent wave packet dynamics studies have been performed for the H 2 + NH 2 → H + NH 3 reaction with a seven-dimensional model on a new interpolated ab initio potential energy surface (PES). The PES is constructed using modified Shepard interpolation Scheme and contains 1967 data points with ab initio calculations carried out on UCCSD(T)/aug-cc-pVTZ level. In the seven-dimensional model, NH 2 group keeps C2v symmetry and two NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated when (1) the two reactants are initially at their ground states; (2) NH 2 bending mode is excited, and (3) H 2 is on its first vibrational excited state. The integral cross sections are also reported for these initial states with centrifugal-sudden approximation. Thermal rate constants are calculated for the temperature range of 200–2000 K and compared with the previous calculated values and available experimental data. Good agreements between theory and experiments for the rate constants at intermediate temperature are achieved on this PES.