Ab Initio Potential Energy Surface Research Articles

An Ab Initio Neural Network Potential Energy Surface for the Dimer of Formic Acid and Further Quantum Tunneling Dynamics.

We construct a full-dimensional ab initio neural network potential energy surface (PES) for the isomerization system of the formic acid dimer (FAD). This is based upon ab initio calculations using the DLPNO-CCSD(T) approach with the aug-cc-pVTZ basis set, performed at over 14000 symmetry-unique geometries. An accurate fit to the obtained energies is generated using a general neural network fitting procedure combined with the fundamental invariant method, and the overall energy-weighted root-mean-square fitting error is about 6.4 cm-1. Using this PES, we present a multidimensional quantum dynamics study on tunneling splittings with an efficient theoretical scheme developed by our group. The ground-state tunneling splitting of FAD calculated with a four-mode coupled method is in good agreement with the most recent experimental measurements. The PES can be applied for further dynamics studies. The effectiveness of the present scheme for constructing a high-dimensional PES is demonstrated, and this scheme is expected to be feasible for larger molecular systems.

Ab initio potential energy surface and vibration-rotation energy levels of magnesium monohydroxide revisited

The accurate potential energy surface of magnesium monohydroxide, MgOH, in its ground electronic state X˜2Σ+ has been determined from ab initio calculations using the coupled-cluster approach in conjunction with the correlation-consistent basis sets up to septuple-zeta quality. The core-electron correlation, higher-order electron correlation, scalar relativistic, and adiabatic effects were taken into account. The equilibrium configuration of the MgOH molecule was confirmed to be linear, although with the bending potential energy function being very flat near the minimum. The vibration–rotation-spin energy levels of the MgOH, MgOD, 25MgOH, and 26MgOH isotopologues were predicted using a variational approach. The spectroscopic constants of these isotopologues were determined to high accuracy.

Rotational state-changes in C5N− by collisions with He and H2

ABSTRACT The anion C5N− is one of the largest linear (C,N)-bearing chains detected in the interstellar medium. Here we present and discuss the general features of new ab initio potential energy surfaces describing the interaction of this linear anion with He and H2. We employ a Legendre Polynomials expansion representation for the former and an artificial neural network fit for the latter. We then carry out quantum scattering calculations to yield rotationally inelastic cross-sections for collisions with He and H2, using relative translational energy values in the range of 0.1–300 cm−1. We then obtained the corresponding inelastic rate coefficients as a function of temperature covering the range from 1 to 100 K. The results for these two systems are compared with each other, as well as with the earlier results on the C3N− colliding with the same partners. We found that the final inelastic rate coefficients for this anion are all fairly large, those from collisions with H2 being the largest. The consequences of such findings on their non-equilibrium rotational populations in interstellar environments are discussed in our conclusions.

Tomography of Feshbach resonance states

Feshbach resonances are fundamental to interparticle interactions and become particularly important in cold collisions with atoms, ions, and molecules. In this work, we present the detection of Feshbach resonances in a benchmark system for strongly interacting and highly anisotropic collisions: molecular hydrogen ions colliding with noble gas atoms. The collisions are launched by cold Penning ionization, which exclusively populates Feshbach resonances that span both short- and long-range parts of the interaction potential. We resolved all final molecular channels in a tomographic manner using ion-electron coincidence detection. We demonstrate the nonstatistical nature of the final-state distribution. By performing quantum scattering calculations on ab initio potential energy surfaces, we show that the isolation of the Feshbach resonance pathways reveals their distinctive fingerprints in the collision outcome.

Signatures of Non-universal Quantum Dynamics of Ultracold Chemical Reactions of Polar Alkali Dimer Molecules with Alkali Metal Atoms: Li(2S) + NaLi(a3Σ+) → Na(2S) + Li2(a3Σu+).

Ultracold chemical reactions of weakly bound triplet-state alkali metal dimer molecules have recently attracted much experimental interest. We perform rigorous quantum scattering calculations with a new ab initio potential energy surface to explore the chemical reaction of spin-polarized NaLi(a3Σ+) and Li(2S) to form Li2(a3Σu+) and Na(2S). The reaction is exothermic and proceeds readily at ultralow temperatures. Significantly, we observe strong sensitivity of the total reaction rate to small variations of the three-body part of the Li2Na interaction at short range, which we attribute to a relatively small number of open Li2(a3Σu+) product channels populated in the reaction. This provides the first signature of highly non-universal dynamics seen in rigorous quantum reactive scattering calculations of an ultracold exothermic insertion reaction involving a polar alkali dimer molecule, opening up the possibility of probing microscopic interactions in atom+molecule collision complexes via ultracold reactive scattering experiments.

Quantum state-resolved molecular dipolar collisions over four decades of energy.

Collisions between cold polar molecules represent a fascinating research frontier but have proven hard to probe experimentally. We report measurements of inelastic cross sections for collisions between nitric oxide (NO) and deuterated ammonia (ND3) molecules at energies between 0.1 and 580 centimeter-1, with full quantum state resolution. At energies below the ~100-centimeter-1 well depth of the interaction potential, we observed backward glories originating from peculiar U-turn trajectories. At energies below 0.2 centimeter-1, we observed a breakdown of the Langevin capture model, which we interpreted in terms of a suppressed mutual polarization during the collision, effectively switching off the molecular dipole moments. Scattering calculations based on an ab initio NO-ND3 potential energy surface revealed the crucial role of near-degenerate rotational levels with opposite parity in low-energy dipolar collisions.

Low-Energy Collisions of Zeeman-Decelerated NH Radicals with He Atoms.

We report an experimental study of state-to-state inelastic scattering of NH (X 3Σ-, N = 0, j = 1) radicals with He atoms. Using a crossed molecular beam apparatus that combines a Zeeman decelerator and velocity map imaging, we study both integral and differential cross sections in the N = 0, j = 1 → N = 2, j = 3 inelastic channel. We developed various new REMPI schemes to state-selectively detect NH radicals, and tested their performance in terms of sensitivity and ion recoil velocity. We found a 1 + 2' + 1' REMPI scheme using the A 3Π ← X 3Σ- resonant transition, which yields acceptable recoil velocities and is more than an order of magnitude more sensitive than conventional one-color REMPI schemes to detect NH. We used this REMPI scheme to probe state-to-state integral and differential cross sections around the channel opening at 97.7 cm-1, as well as at higher energies where structure in the scattering images could be resolved. The experimental results are in excellent agreement with the predictions from quantum scattering calculations which are based on an ab initio NH-He potential energy surface.

Theoretical Calculations of the Thermal Rate Coefficients for the Interstellar NH3+ + H2 → NH4+ + H Reaction on a New Δ-Machine Learning Potential Energy Surface

The NH3+ + H2 → NH4+ + H ion–molecule hydrogen-exchange reaction is important for understanding the mechanisms of ammonia formation in interstellar medium and circumstellar environments. The thermal rate coefficients of this reaction are calculated using the ring polymer molecular dynamics approach, which can account for the nuclear quantum effects, on a newly developed ab initio potential energy surface (PES). The PES was constructed by fitting a large number of ab initio points at the MP2/cc-pVDZ and CCSD(T)/cc-pVTZ levels using a Δ-machine learning approach. The obtained rate coefficients exhibited a minimum at T ∼ 130 K, and this characteristic temperature-dependent behavior qualitatively agreed with the experimental measurements. This study confirms that quantum tunneling plays a significant role in the reaction dynamics along with the existence of the metastable NH3+···H2 pre-reactive complex in the entrance channel of the PES.

Collisional excitation of HCNH+ by He and H2: New potential energy surfaces and inelastic rate coefficients.

Protonated molecules have been increasingly detected in the interstellar medium (ISM), and usually astrochemical models fail at reproducing the abundances derived from observational spectra. Rigorous interpretation of the detected interstellar emission lines requires prior calculations of collisional rate coefficients with H2 and He, i.e., the most abundant species in the ISM. In this work, we focus on the excitation of HCNH+ induced by collision with H2 and He. Therefore, we first calculate ab initio potential energy surfaces (PESs) using the explicitly correlated and standard coupled cluster method with single, double, and non-iterative triple excitation in conjunction with the augmented-correlation consistent-polarized valence triple zeta basis set. Both the HCNH+-H2 and HCNH+-He potentials are characterized by deep global minima of 1426.60 and 271.72 cm-1, respectively, and large anisotropies. From these PESs, we derive state-to-state inelastic cross sections for the 16 low-lying rotational energy levels of HCNH+ using the quantum mechanical close-coupling approach. The differences between cross sections due to ortho- and para-H2 impacts turn out to be minor. Using a thermal average of these data, we retrieve downward rate coefficients for kinetic temperatures of up to 100K. As it could be anticipated, differences of up to two orders of magnitude exist between the rate coefficients induced by H2 and He collisions. We expect that our new collision data will help to improve the disagreement between abundances retrieved from observational spectra and astrochemical models.

Rotational Excitation of NCCN by p-H2(jc = 0) at Low Temperatures

NCCN is a potential candidate for detection in interstellar space, but to model its emission spectra, accurate collisional rates are required. Here, we present a first truncated ab initio four-dimensional potential energy surface of NCCN–H2, generated using explicitly correlated coupled-cluster singles and doubles with the perturbative triple excitations [CCSD(T)-F12b] method and aug-cc-pVTZ basis set. The surface allows NCCN to rotate in the x–z plane, and H2 is restricted to three orientations [(a) parallel to the z axis, (b) perpendicular to the z axis but still in the x–z plane, and (c) perpendicular to the z axis and also perpendicular to the x–z plane] and averaged to obtain a reduced two-dimensional (2D) surface. To augment data points, the surface is fitted into a neural networks model. Multipole expansion coefficients for this reduced 2D surface of NCCN–H2 have been calculated and utilized in determining collisional cross sections using the close-coupling method for collision energies until 1500 cm–1. Finally, the rate coefficients for reduced NCCN–H2 have been computed. These rate coefficients closely approximate the collision of NCCN with p-H2(jc = 0), and they are greater by 1.8–2.8 times than the rates of NCCN–He previously investigated.

Predicting Structural Properties of Pure Silica Zeolites Using Deep Neural Network Potentials.

Machine learning potentials (MLPs) capable of accurately describing complex ab initio potential energy surfaces (PESs) have revolutionized the field of multiscale atomistic modeling. In this work, using an extensive density functional theory (DFT) data set (denoted as Si-ZEO22) consisting of 219 unique zeolite topologies (350,000 unique DFT calculations) found in the International Zeolite Association (IZA) database, we have trained a DeePMD-kit MLP to model the dynamics of silica frameworks. The performance of our model is evaluated by calculating various properties that probe the accuracy of the energy and force predictions. This MLP demonstrates impressive agreement with DFT for predicting zeolite structural properties, energy-volume trends, and phonon density of states. Furthermore, our model achieves reasonable predictions for stress-strain relationships without including DFT stress data during training. These results highlight the ability of MLPs to capture the flexibility of zeolite frameworks and motivate further MLP development for nanoporous materials with near-ab initio accuracy.

Regularized by Physics: Graph Neural Network Parametrized Potentials for the Description of Intermolecular Interactions.

Simulations of molecular systems using electronic structure methods are still not feasible for many systems of biological importance. As a result, empirical methods such as force fields (FF) have become an established tool for the simulation of large and complex molecular systems. The parametrization of FF is, however, time-consuming and has traditionally been based on experimental data. Recent years have therefore seen increasing efforts to automatize FF parametrization or to replace FF with machine-learning (ML) based potentials. Here, we propose an alternative strategy to parametrize FF, which makes use of ML and gradient-descent based optimization while retaining a functional form founded in physics. Using a predefined functional form is shown to enable interpretability, robustness, and efficient simulations of large systems over long time scales. To demonstrate the strength of the proposed method, a fixed-charge and a polarizable model are trained on ab initio potential-energy surfaces. Given only information about the constituting elements, the molecular topology, and reference potential energies, the models successfully learn to assign atom types and corresponding FF parameters from scratch. The resulting models and parameters are validated on a wide range of experimentally and computationally derived properties of systems including dimers, pure liquids, and molecular crystals.

Elastic and glancing-angle rate coefficients for heating of ultracold Li and Rb atoms by collisions with room-temperature noble gases, H2, and N2.

Trapped ultracold alkali-metal atoms can be used to measure pressure in the ultra-high-vacuum and XHV pressure regimes, those with p < 10-6 Pa. This application for ultracold atoms relies on precise knowledge of collision rate coefficients of alkali-metal atoms with residual room-temperature atoms and molecules in the ambient vacuum or with deliberately introduced gasses. Here, we determine combined elastic and inelastic rate coefficients as well as glancing-angle rate coefficients for ultracold 7Li and 87Rb with room-temperature noble gas atoms as well as H2 and 14N2 molecules. Glancing collisions are those processes where only little momentum is transferred to the alkali-metal atom and this atom is not ejected from its trap. Rate coefficients are found by performing quantum close-coupling scattering calculations using ab initio ground-state electronic Born-Oppenheimer potential energy surfaces. The potentials for Li and Rb with noble gas atoms and also for Rb(2S)-H2(XΣg +) and Rb(2S)-N2(X1Σg +) systems are based on the non-relativistic spin-restricted coupled-cluster method with single, double, and noniterative triple excitations [RCCSD(T)]. For Li(2S)-N2(X1Σg +), the potential is computed at the explicitly correlated spin-restricted RCCSD(T)-F12 level. For Rb, Kr, and Xe atoms, scalar relativistic corrections to the core electrons have been included, while second-order spin-orbit corrections from the valence electrons have been estimated. Data for Li-H2 and Li-He were taken from the existing literature. We estimate standard uncertainties of the rate coefficients by comparing rate coefficients calculated using potentials found with electronic basis sets of increasing size, including estimates of relativistic spin-orbit corrections and the uncertainty of the van der Waals coefficients. The relative uncertainties of rate coefficients are 1%-2% with the exception of 7Li or 87Rb colliding with 20Ne. Those have relative uncertainties of 9% and 8%, respectively. We also show that a commonly used semiclassical approximation for the total elastic rate coefficient agrees with the quantum calculations to 10% with the exception of 7Li and 87Rb collisions with H2, where the semiclassical value underestimates the quantum value by 20%.

Machine learning of double-valued nonadiabatic coupling vectors around conical intersections.

In recent years, machine learning has had an enormous success in fitting ab initio potential-energy surfaces to enable efficient simulations of molecules in their ground electronic state. In order to extend this approach to excited-state dynamics, one must not only learn the potentials but also nonadiabatic coupling vectors (NACs). There is a particular difficulty in learning NACs in systems that exhibit conical intersections, as due to the geometric-phase effect, the NACs may be double-valued and are, thus, not suitable as training data for standard machine-learning techniques. In this work, we introduce a set of auxiliary single-valued functions from which the NACs can be reconstructed, thus enabling a reliable machine-learning approach.

Vibrational mode-specific quasi-classical trajectory studies for the two-channel HI + C2H5 reaction.

We report a detailed dynamics study on the mode-specificity of the HI + C2H5 two-channel reaction (H-abstraction and I-abstraction), through performing quasi-classical trajectory computations on a recently developed high-level ab initio full-dimensional spin-orbit-corrected potential energy surface, by exciting four different vibrational modes of reactants at five collision energies. The effect of the normal-mode excitations on the reactivity, the mechanism, and the post-reaction energy flow is investigated. Both reaction pathways are intensely promoted when the HI-stretching mode is excited while the excitations imposed on C2H5 somewhat surprisingly inhibit the dominant H-abstraction reaction pathway. The enhancement effect of the excitation in the HI vibrational mode is found to be much more effective than increasing the translational energy, similar to the HBr + C2H5 reaction. Not like the Br-abstraction pathway, however, the I-abstraction reaction pathway could be comparable to the dominant H-abstraction reaction pathway. The dominance of the direct stripping mechanism is indicated in H-abstraction while the direct rebound mechanism is observed in I-abstraction. The H-abstraction is much pickier about the initial attack angle distributions for HI than I-abstraction is, which leads to a decrease in reactivity in the H-abstraction reaction pathway. The dominance of side-on CH3CH2 attack in I-abstraction is more obvious than in the case of H-abstraction. In the case of the H-abstraction reaction pathway, the major part of the initial translational energy ends up in translational recoil, while for I-abstraction most energy excites the product C2H5I.

Potential energy surfaces for singlet and triplet states of the LiH2+ system and quasi-classical trajectory cross sections for H + LiH+ and H+ + LiH.

A new set of six accurate ab initio potential energy surfaces (PESs) is presented for the first three singlet and triplet states of LiH2+ (1,21A', 11A'', 1,23A', and 13A'' states, where four of them are investigated for the first time), which have allowed new detailed studies gaining a global view on this interesting system. These states are relevant for the study of the most important reactions of lithium chemistry in the early universe. More than 45 000 energy points were calculated using the multi-reference configuration interaction level of theory using explicitly correlated methods (ic-MRCI-F12), and the results obtained for each individual electronic state were fitted to an analytical function. Using quasiclassical trajectories and considering the initial diatomic fragment in the ground rovibrational state, we have determined the integral cross sections for the H + LiH+(X2Σ+, C2Π) and H+ + LiH(X1Σ+, B1Π) reactions. In these calculations all available reaction channels were considered: the chemically most important H or H+ transfer/abstraction as well as atom exchange and collision induced dissociation for up to 1.0 eV of collision energy.

Theoretical vibrational mode-specific dynamics studies for the HBr + C2H5 reaction.

A quasi-classical trajectory (QCT) study is performed for the HBr + C2H5 multi-channel reaction using a recently-developed high-level ab initio full-dimensional spin-orbit-corrected potential energy surface (PES) by exciting five different vibrational modes of reactants at five collision energies. The effect of the normal-mode excitations on the reactivity, the mechanism, and the post-reaction energy flow is followed. A significant decrease of the reactivity caused by the longer initial distances of the reactants for the vHBr = 1 reaction at low collision energy (Ecoll) is observed due to the intramolecular vibrational-energy redistribution and the classical nature of the QCT method. All of the three reaction pathways (H-abstraction, Br-abstraction, and H-exchange) are intensely promoted when the HBr-stretching mode is excited. No clear promotion is observed when excitation is imposed to C2H5 except that asymmetric CH-stretching helps the H-exchange process. The enhancement effect of the excitation in the HBr vibrational mode is found to be much more effective than increasing the translational energy, in contrast to the HBr + CH3 reaction. The forward scattering mechanism can be clearly promoted by the excitation of the HBr-stretching mode, or by the high collision energy, indicating the dominance of the direct stripping mechanism in these cases. At low collision energy with no excitation or excitation of any vibrational mode of C2H5, the forward scattering feature is less obvious. At Ecoll = 1 kcal mol-1, when HBr-stretching is excited, the product clearly gains more relative translational energy. However, it is interesting to see that when the excitation is in C2H5, the effect is the opposite, i.e., the product gains less relative translational energy compared to the ground-state reaction.

New ab initio potential energy surface of NaFH (1A') system and quantum dynamics studies for the Na + HF (v, j) → NaF + H reaction.

A global potential energy surface (PES) for the electronic ground state of the Na + HF reactive system is constructed by three-dimensional cubic spline interpolation of 37 000 ab initio points obtained using the multireference configuration interaction method including the Davidson's correction (MRCI + Q) with auc-cc-pV5Z basis set. The endoergicity, well depth and properties of the separated diatomic molecules are in good agreement with experimental estimations. Quantum dynamics calculations have been performed and compared with those of the previous MRCI PES as well as experimental values. The better agreement between theory and experiment indicates the accuracy of the new PES.

First-principles predictions for shear viscosity of air components at high temperature.

The direct molecular simulation (DMS) method is used to obtain shear viscosity data for non-reacting air and its components by simulating isothermal, plane Poiseuille subsonic flows. Shear viscosity is estimated at several temperatures, from 273 K to 10 000 K, by fitting the DMS velocity profiles using the analytic solution of the Navier-Stokes equations for this simple canonical flow. The ab initio potential energy surfaces (PESs) that describe the various atomic-level interactions are the only input in the simulations. Molecules involved in a collision within the flow can occupy any rovibrational state that is allowed by the effective diatomic potential. For molecular nitrogen, oxygen, and air at standard condition molar composition, the DMS shear viscosity predictions are in excellent agreement with the experimental data that are available up to about 2000 K. The results for pure molecular nitrogen and pure molecular oxygen also agree very well with previously published quasi-classical trajectory (QCT) calculations based on the same PESs. It is further shown that the ab initio shear viscosity data are generally lower than the corresponding values used in popular computational fluid dynamics codes, over a wide temperature range. Finally, Wilke's mixing rule is demonstrated to accurately predict the DMS air viscosity results from the pure molecular components data up to 4000 K.

Quantum dynamics calculations on isotope effects of hydrogen transfer isomerization in formic acid dimer

We present a quantum dynamics study on the isotope effects of hydrogen transfer isomerization in the formic acid dimer, and this is achieved by multidimensional dynamics calculations with an efficient quantum mechanical theoretical scheme developed by our group, on a full-dimensional neural network ab initio potential energy surface. The ground-state and fundamental tunneling splittings for four deuterium isotopologues of formic acid dimer are considered, and the calculated results are in very good general agreement with the available experimental measurements. Strong isotope effects are revealed, the mode-specific fundamental excitation effects on the tunneling rate are evidently influenced by the deuterium substitution of H atom with the substitution on the OH bond being more effective than on the CH bond. Our studies are helpful for acquiring a better understanding of isotope effects in the double-hydrogen transfer processes.