Neural network potential energy surfaces and dipole moment surfaces for SO2(H2O) and SO2(H2O)2 complexes

Full-dimensional, ab initio-based many-body potential energy surfaces and dipole moment surfaces constructed using the neural network method for SO2(H2O)n (n = 1,2) complexes are reported. The database of the SO2 1-body PES, SO2(H2O) 2-body PES and SO2(H2O)2 3-body PES consists of 11 952, 79 882 and 84 159 ab initio energies, respectively. All 1-body energies were calculated at the CCSD(T)/CBS(AVTZ:AVQZ) level and all 2,3-body energies were calculated at the DSD-PBEP86/AVTZ level. The database of DMSs is the same as that of PESs and all dipole moments were calculated at the MP2/AVTZ level. Harmonic frequencies and dissociation energies of SO2(H2O) and SO2(H2O)2 were calculated on these PESs and compared with ab initio results to examine the fidelity of these PESs.

The MD17 datasets from the perspective of datasets for gas-phase “small” molecule potentials

There has been great progress in developing methods for machine-learned potential energy surfaces. There have also been important assessments of these methods by comparing so-called learning curves on datasets of electronic energies and forces, notably the MD17 database. The dataset for each molecule in this database generally consists of tens of thousands of energies and forces obtained from DFT direct dynamics at 500 K. We contrast the datasets from this database for three “small” molecules, ethanol, malonaldehyde, and glycine, with datasets we have generated with specific targets for the potential energy surfaces (PESs) in mind: a rigorous calculation of the zero-point energy and wavefunction, the tunneling splitting in malonaldehyde, and, in the case of glycine, a description of all eight low-lying conformers. We found that the MD17 datasets are too limited for these targets. We also examine recent datasets for several PESs that describe small-molecule but complex chemical reactions. Finally, we introduce a new database, “QM-22,” which contains datasets of molecules ranging from 4 to 15 atoms that extend to high energies and a large span of configurations.

Theoretical study of the O(3P) + SiH4 reaction: global potential energy surface, kinetics and dynamics study

In order to understand the gas-phase hydrogen abstraction reaction between O(³P) and silane we began by developing the first full-dimensional analytical potential energy surface, named PES-2022. It is basically a valence bond function augmented with molecular mechanic terms describing in an intuitive way stretching and bending nuclei motions, and it is fitted to high level ab initio calculations. The surface presents continuous and smooth behaviour, with analytical first energy derivatives, on which the hydrogen atoms in silane are permutationally symmetric. Based on PES-2022, a kinetics study was performed using the variational transition-state theory with multidimensional tunnelling corrections in the temperature range of 300-1000 K. We observed that experimental and theoretical results show widely spread results, both in absolute value and temperature dependence, possibly because they include the reactivity from both O(³P) and O(¹D) electronic states, which present different mechanisms and multiple channels. When the comparison is performed on the same footing, O(³P) + SiH₄ → HO + SiH₃, the present results agree with Ding and Marshall's experiments and with Zhang et al.'s theoretical rate constants. The kinetic isotope effects (KIEs) reproduced the only experimental value, improving previous theoretical results. Finally, a dynamics study was performed on PES-2022 using quasi-classical trajectory calculations under two different initial conditions, at fixed room temperature and at a fixed collision energy of 8.0 kcal mol^-1. In the first case, the available energy deposited as HO(v) vibration was 47%, with population inversion, P(v = 0)/P(v = 1) = 11/89%, reproducing the experimental evidence. In the second case, the experimental product translational distribution was reasonably simulated, while the angular product distribution presented opposite behaviour, backward versus forward. On analysing this discrepancy, we found that while in the present work the O(³P) + SiH₄ reaction was reported, in the experiment both O(³P) and O(¹D) electronic states are reported. So, the comparison was not performed on the same footing. In sum, agreement of the present results with experiments permits us to be reasonably optimistic about the quality and accuracy of the new PES, and at the same time to highlight the fact that theory/experiment comparisons must be performed on the same footing.

Mechanism and Kinetics of Methane Combustion. Part II: Potential Energy Surface for Hydrogen-Abstraction Reaction of CH4 + O(3P)

Methane combustion plays an important role in various fields such as combustion chemistry and atmospheric chemistry of the stratosphere. Highly accurate study of its initial reaction remains a key challenge. Here, through extensive studies with a state-of-the-art ab initio and neural network method, we present a potential energy surface of the O(³P) + CH₄ → OH + CH₃ reaction on the ground state 1³A and the first excited state 2³A. In this work, the energies of 10 167 points covering all important regions are obtained with state-averaged complete active space self-consistent field calculations and then fitted using the Levenberg-Marquardt algorithm with a root-mean-square error of 0.391 and 0.442 kcal/mol for the 1³A and 2³A states, respectively. This study explores the characteristics of the radical van der Waals (VdW) complex and reveals a detailed mechanism of the methane combustion initial reaction. Within the scope of this mechanism, this surface gives a fairly accurate description of the regions around the saddle point, conical intersection, and vdW wells in the entrance for efficient computational simulations. As a theoretical study on a prototypical polyatomic reaction, it is hopeful that this work will modify our understanding of the primary process in hydrocarbon combustion.

When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111)

The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.

A global ab initio potential energy surface and dynamics of the proton-transfer reaction: OH− + D2 → HOD + D−

The reaction mechanisms of OH⁻ + D₂ → HOD + D⁻ were first revealed by theory, based on an accurate full-dimensional PES.

When Classical Trajectories Get to Quantum Accuracy: The Scattering of H2 on Pd(111)

When elementary reactive processes occur at such low energies that only a few states of reactants and/or products are available, quantum effects strongly manifest and the standard description of the dynamics within the classical framework fails. We show here, for H₂ scattering on Pd(111), that by pseudoquantizing in the spirit of Bohr the relevant final actions of the system, along with adequately treating the diffraction-mediated trapping of the incoming wave, classical simulations achieve an unprecedented agreement with state-of-the-art quantum dynamics calculations.

Theoretical Study of Barrierless Chemical Reactions Involving Nearly Elastic Rebound: The Case of S(1D) + X2, X = H, D

For some values of the total angular momentum consistent with reaction, the title processes involve nonreactive trajectories proceeding through a single rebound mechanism during which the internal motion of the reagent diatom is nearly unperturbed. When such paths are in a significant amount, the classical reaction probability is found to be markedly lower than the quantum mechanical one. This finding was recently attributed to an unusual quantum effect called diffraction-mediated trapping, and a semiclassical correction was proposed in order to take into account this effect in the classical trajectory method. In the present work, we apply the resulting approach to the calculation of opacity functions as well as total and state-resolved integral cross sections (ICSs) and compare the values obtained with exact quantum ones, most of which are new. As the title reactions proceed through a deep insertion well, mean potential statistical calculations are also presented. Seven values of the collision energy, ranging from 30 to 1127 K, are considered. Two remarkable facts stand out: (i) The corrected classical treatment strongly improves the accuracy of the opacity function as compared to the usual classical treatment. When the entrance transition state is tight, however, those trajectories crossing it with a bending vibrational energy below the zero point energy must be discarded. (ii) The quantum opacity function, particularly its cutoff, is finely reproduced by the statistical approach. Consequently, the total ICS is also very well described by the two previous approximate methods. These, however, do not predict state-resolved ICSs with the same accuracy, proving thereby that (i) one or several genuine quantum effects involved in the dynamics are missed by the corrected classical treatment and (ii) the dynamics are not fully statistical.

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.

An accurate full-dimensional potential energy surface and quasiclassical trajectory dynamics of the H + H2O2 two-channel reaction

Quasiclassical trajectory calculations reveal interesting dynamics features based on an accurate FI-NN PES for the H + H₂O₂ two-channel reaction.

Thermal Rate Constants for the O(3P) + CH4 → OH + CH3 Reaction: The Effects of Quantum Tunneling and Potential Energy Barrier Shape

The rate constants and kinetic isotope effects for the O(³P) + CH₄ reaction have been investigated with the quantum instanton method in full dimensionality. The calculated rate constants are in good agreement with the experimental values above 400 K, below which the measured values are scattered. Compared to other theoretical approaches, the quantum instanton method predicts the largest quantum tunneling effect, so it gives the largest rate constants at low temperatures. The calculated kinetic isotope effects are always much larger than 1 and increase with decreasing temperature, due to the zero-point energy and quantum tunneling. Our calculations on different potential energy surfaces demonstrate that the potential energy barrier shape dominates the magnitude of quantum tunneling and has a great effect on the kinetic isotope effect.

Theoretical and experimental investigations of rate coefficients of O(1D) + CH4 at low temperature

The O(¹D) + CH₄ rates from RPMD and VUV-LIF are close to each other in the temperature range of the atmosphere.