Ab Initio Potential Energy Surface Research Articles

Dynamics of the HCl + C2H5 Multichannel Reaction on a Full-Dimensional Ab Initio Potential Energy Surface.

We report a full-dimensional ab initio analytical potential energy surface (PES), which accurately describes the HCl + C2H5 multichannel reaction. The new PES is developed by iteratively adding selected configurations along HCl + C2H5 quasi-classical trajectories (QCTs), thereby improving our previous Cl(2P3/2) + C2H6 PES using the Robosurfer program package. QCT simulations for the H'Cl + C2H5 reaction reveal hydrogen-abstraction, chlorine-abstraction, and hydrogen-exchange channels leading to Cl + C2H5H', H' + C2H5Cl, and HCl + C2H4H', respectively. Hydrogen abstraction dominates in the collision energy (Ecoll) range of 1-80 kcal/mol and proceeds with indirect isotropic scattering at low Ecoll and forward-scattered direct stripping at high Ecoll. Chlorine abstraction opens around 40 kcal/mol collision energy and becomes competitive with hydrogen abstraction at Ecoll = 80 kcal/mol. A restricted opening of the cone of acceptance in the Cl-abstraction reaction is found to result in the preference for a backward-scattering direct-rebound mechanism at all energies studied. Initial attack-angle distributions show mainly side-on collision preference of C2H5 for both abstraction reactions, and in the case of the HCl reactant, H/Cl-side preference for the H/Cl abstraction. For hydrogen abstraction, the collision energy transfer into the product translational and internal energy is almost equally significant, whereas in the case of chlorine abstraction, most of the available energy goes into the internal degrees of freedom. Hydrogen exchange is a minor channel with nearly constant reactivity in the Ecoll range of 10-80 kcal/mol.

Theoretical Insights into the Vibrational Structure of Carbon Dioxide Rare-Gas Complexes.

Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar-CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.

Feature selection for high-dimensional neural network potentials with the adaptive group lasso

Neural network potentials are a powerful tool for atomistic simulations, allowing to accurately reproduce ab initio potential energy surfaces with computational performance approaching classical force fields. A central component of such potentials is the transformation of atomic positions into a set of atomic features in a most efficient and informative way. In this work, a feature selection method is introduced for high dimensional neural network potentials, based on the adaptive group lasso (AGL) approach. It is shown that the use of an embedded method, taking into account the interplay between features and their action in the estimator, is necessary to optimize the number of features. The method’s efficiency is tested on three different monoatomic systems, including Lennard–Jones as a simple test case, Aluminium as a system characterized by predominantly radial interactions, and Boron as representative of a system with strongly directional components in the interactions. The AGL is compared with unsupervised filter methods and found to perform consistently better in reducing the number of features needed to reproduce the reference simulation data at a similar level of accuracy as the starting feature set. In particular, our results show the importance of taking into account model predictions in feature selection for interatomic potentials.

Quantum dynamics of O(3P) + HBr → OH + Br reaction: integral cross sections and rate constants and their initial state dependence.

The integral cross sections (ICSs) and rate constants (RCs) for the title reaction are calculated by means of accurate quantum wave packet method employing the recent ab initio potential energy surface developed by Peterson and co-workers. Hundreds of partial wave contributions (up to J=150) are calculated explicitly taking Coriolis coupling into account. The ICSs are found to increase monotonically with energy and their energy dependences are smooth except in the energy range below the potential energy barrier, where several resonance peaks are observed. The temperature dependences of the RCs obey in general the Arrhenius law. These behaviors manifest the dominance of abstract mechanism. Initial rotational state (ji) excitations are found to greatly enhance the reactivity starting from ji=4, while exhibiting a complicated dependence for ji<4. From Boltzmann average of the initial state specified RCs (for ji=0-15) we obtained thermal RCs which are in reasonable agreement with available experimental results.

Accurate reference spectra of HD in an H2--He bath for planetary applications

The hydrogen deuteride (HD) molecule is an important deuterium tracer in astrophysical studies. The atmospheres of gas giants are dominated by molecular hydrogen, and the simultaneous observation of H$_2$ and HD lines provides reliable information on the D/H ratios on these planets. The reference spectroscopic parameters play a crucial role in such studies. Under the thermodynamic conditions encountered in these atmospheres, spectroscopic studies of HD require not only the knowledge of line intensities and positions but also accurate reference data on pressure-induced line shapes and shifts. Our aim is to provide accurate collision-induced line-shape parameters for HD lines that cover any thermodynamic conditions relevant to the atmospheres of giant planets, namely any relevant temperature, pressure, and perturbing gas composition (the H$_2$--He mixture). We performed quantum-scattering calculations on our new, highly accurate ab initio potential energy surface (PES) and we used scattering S matrices obtained in this way to determine the collision-induced line-shape parameters. We used cavity ring-down spectroscopy to validate our theoretical methodology. We report accurate collision-induced line-shape parameters for the pure rotational R(0), R(1), and R(2) lines, the most relevant HD lines for investigations of the atmospheres of the giant planets. Besides the basic Voigt-profile collisional parameters (i.e., the broadening and shift parameters), we also report their speed dependences and the complex Dicke parameter, which can influence the effective width and height of the HD lines up to almost a factor of 2 for giant planet conditions. The sub-percent-level accuracy reached in this work is a considerable improvement over previously available data. All the reported parameters (and their temperature dependences) are consistent with the HITRAN database format, hence allowing for the use of the HITRAN Application Programming Interface (HAPI) for generating the beyond-Voigt spectra of HD.

Global Potential Energy Surfaces by Compressed-State Multistate Pair-Density Functional Theory for Hyperthermal Collisions in the O2+O2 System.

Interactions between oxygen molecules play an important role in atmospheric chemistry and hypersonic flow chemistry in atmospheric entries. Recently, high-quality ab initio potential energy surface (PES) of the quintet O4 was reported by Paukku et al. [J. Chem. Phys. 147, 034301 (2017)]. 10543 configurations were sampled and calculated at the level of MS-CASPT2/maug-cc-pVTZ with scaled external correlation. The PES was fitted to a many-body (MB) form with the many-body part described by the permutationally invariant polynomial approach (MB-PIP). In this work, the PIP-Neural Network (PIP-NN) and MB-PIP-NN methods were used to refit the PES based on the same data by Paukku et al. Three PESs were compared. It was found that the performances differ significantly in the O+O3 region as well as in the long-range region. Therefore, additional 1300 points were sampled, and the efficient compressed-state multistate pair-density functional theory (CMS-PDFT) was used to calculate the electronic structure of these 1300 points and 10543 points by Paukku et al. Then, a completely new quintet PES was fitted using the MB-PIP-NN method. Based on this PES, the quasi-classical trajectory (QCT) approach was used to reveal all possible reaction channels for hyperthermal O2-O2 collisions.

State-to-state study of non-equilibrium recombination of oxygen and nitrogen molecules.

Rapidly cooled mixtures are of interest for several applications, including hypersonic flows due to the presence of strong cooling temperature gradients in regions such as hypersonic boundary layers and expanding nozzles. There have been very few studies of rapidly cooled mixtures using the high-fidelity rovibrational databases afforded by abinitio potential energy surfaces. This work makes use of existing rovibrational state-specific databases to study rapidly cooled mixtures. In particular, we seek to understand the importance of thermal non-equilibrium in recombining mixtures using both rovibrational and vibrational state-to-state methods for oxygen and nitrogen molecules. We find that although there is significant non-equilibrium during recombination, it is well captured by the vibrational state-specific approach. Finally, we compare the global recombination rate computed based on the state-specific recombination rate coefficients and the global recombination rate computed based on the time local dissociation rate coefficient, which is reversed using the principle of detailed balance. The local dissociation rate coefficient is computed by weighting the state-specific dissociation rate coefficients with the state-specific distribution of energy states. We find a large difference between these rates, highlighting a potential source of errors in hypersonic flow predictions.

A new ab initio potential energy surface and rovibrational spectra for the N2–N2O complex

We propose a new potential energy surface for the N2–N2O complex, which encompasses the Q3 normal mode for the ν3 antisymmetric stretching coordinate of N2O. The interaction potential energies were obtained at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations [CCSD(T)-F12a] level with aug-cc-pVTZ basis set plus mid-bond functions. Four-dimensional vibrationally averaged potentials with N2O in the ground and ν3 excited states were generated through integration over the Q3 coordinate. Both potentials exhibit two equivalent skew T-shaped global minima. The rovibrational energies were calculated by employing the radial discrete variable representation/angular finite basis representation approach. The determined fundamental band origin shift of the infrared spectra in the N2O ν3 region, transition frequencies and rotational constants are all consistent with the experimental values.

A New Ab Initio Potential Energy Surface and Rovibrational Spectra for the CO-N2O Complex.

We constructed a new ab initio potential energy surface (PES) for CO-N2O which includes the intramolecular Q3 normal coordinate for the N2O ν3 antisymmetric stretching vibration. The intermolecular potential was evaluated employing the supermolecular method at the [CCSD(T)]-F12a level, with the aug-cc-pVTZ basis set plus bond functions. By integral over the intramolecular Q3 coordinate, we obtained the vibrationally averaged PESs for the CO-N2O system in the ground and ν3 excited states of N2O. Each PES features one nearly T-shaped global minimum and one skewed T-shaped local minimum. Based on these obtained PESs of CO-N2O, the radial discrete variable representation/angle finite base representation method and the Lanczos algorithm were applied for the calculations of bound states and rovibrational energy levels. The calculated ν3 vibrational band origin shift of the N2O monomer in CO-N2O is 2.7570 cm-1, matching well with the observed value of 2.9048 cm-1. The computed microwave and infrared transition frequencies, as well as the rotational parameters, are consistent with the experimental observations.

Non-IRC Mechanism of Bimolecular Reactions with Submerged Barriers: A Case Study of Si+ + H2O Reaction.

Chemical reactions with submerged barriers may feature interesting dynamic behaviors that are distinct from those with substantial barriers or those entirely dominated by capture. The Si+ + H2O reaction is a prototypical example, involving even two submerged saddle points along the reaction path: one for the direct dissociation of H (H-dissociation SP) and another for H migration from the O-side to the Si-side (H-migration SP). We investigated the intricacies of this process by employing quasi-classical trajectory calculations on an accurate, full-dimensional ab initio potential energy surface. Through careful trajectory analysis, an interesting nonintrinsic reaction coordinate mechanism was found to play an important role in producing SiOH+ and H. This new pathway is featured as that the H atoms do not form HSiOH+ complexes along the minimum-energy path but directly dissociate into the products after passing through the H-migration SP. Furthermore, based on artificially scaled potential energy surfaces (PESs), the impact of barrier height on the reaction is also explored. This work provides new insights into the dynamics of the Si+ + H2O reaction and enriches our understanding of reactions with submerged barriers.

Translational eigenstates of He@C60 from four-dimensional abinitio potential energy surfaces interpolated using Gaussian process regression.

We investigate the endofullerene system 3He@C60 with a four-dimensional potential energy surface (PES) to include the three He translational degrees of freedom and C60 cage radius. We compare second order Møller-Plesset perturbation theory (MP2), spin component scaled-MP2, scaled opposite spin-MP2, random phase approximation (RPA)@Perdew, Burke, and Ernzerhof (PBE), and corrected Hartree-Fock-RPA to calibrate and gain confidence in the choice of electronic structure method. Due to the high cost of these calculations, the PES is interpolated using Gaussian Process Regression (GPR), owing to its effectiveness with sparse training data. The PES is split into a two-dimensional radial surface, to which corrections are applied to achieve an overall four-dimensional surface. The nuclear Hamiltonian is diagonalized to generate the in-cage translational/vibrational eigenstates. The degeneracy of the three-dimensional harmonic oscillator energies with principal quantum number n is lifted due to the anharmonicity in the radial potential. The (2l + 1)-fold degeneracy of the angular momentum states is also weakly lifted, due to the angular dependence in the potential. We calculate the fundamental frequency to range between 96 and 110 cm-1 depending on the electronic structure method used. Error bars of the eigenstate energies were calculated from the GPR and are on the order of ∼±1.5 cm-1. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two abinitio methods MP2 and RPA@PBE show the best agreement. While MP2 has better agreement than RPA@PBE, due to its higher computational efficiency and comparable performance, we recommend RPA as an alternative electronic structure method of choice to MP2 for these systems.

Hyperfine and Zeeman interactions in ultracold collisions of molecular hydrogen with atomic lithium.

We present a rigorous quantum scattering study of the effects of hyperfine and Zeeman interactions on cold Li-H2 collisions in the presence of an external magnetic field using a recent abinitio potential energy surface. We find that the low-field-seeking states of H2 predominantly undergo elastic collisions: the ratio of elastic-to-inelastic cross sections exceeds 100 for collision energies below 100 mK. Furthermore, we demonstrate that most inelastic collisions conserve the space-fixed projection of the nuclear spin. We show that the anisotropic hyperfine interaction between the nuclear spin of H2 and the electron spin of Li can have a significant effect on inelastic scattering in the ultracold regime, as it mediates two processes: the electron spin relaxation in lithium and the nuclear spin-electron spin exchange. Given the predominance of elastic collisions and the propensity of inelastic collisions to retain H2 in its low-field-seeking states, our results open up the possibility of sympathetic cooling of molecular hydrogen by atomic lithium, paving the way for future exploration of ultracold collisions and high-precision spectroscopy of H2 molecules.

Ab Initio Potential Energy Surface for NaCl-H2 with Correct Long-Range Behavior.

We report a full dimensional ab initio potential energy surface for NaCl-H2 based on precise fitting of a large data set of CCSD(T)/aug-cc-pVTZ energies. A major goal of this fit is to describe the very long-range interaction accurately. This is done in this instance via the dipole-quadrupole interaction. The NaCl dipole and the H2 quadrupole are available through previous works over a large range of internuclear distances. We use these to obtain exact effect charges on each atom. Diffusion Monte Carlo calculations are done for the ground vibrational state using the new potential.

Vibrational mode-specificity in the dynamics of the OH- + CH3I multi-channel reaction.

We report a comprehensive characterization of the vibrational mode-specific dynamics of the OH- + CH3I reaction. Quasi-classical trajectory simulations are performed at four different collision energies on our previously-developed full-dimensional high-level abinitio potential energy surface in order to examine the impact of four different normal-mode excitations in the reactants. Considering the 11 possible pathways of OH- + CH3I, pronounced mode-specificity is observed in reactivity: In general, the excitations of the OH- stretching and CH stretching exert the greatest influence on the channels. For the SN2 and proton-abstraction products, the reactant initial attack angle and the product scattering angle distributions do not show major mode-specific features, except for SN2 at higher collision energies, where forward scattering is promoted by the CI stretching and CH stretching excitations. The post-reaction energy flow is also examined for SN2 and proton abstraction, and it is unveiled that the excess vibrational excitation energies rather transfer into the product vibrational energy because the translational and rotational energy distributions of the products do not represent significant mode-specificity. Moreover, in the course of proton abstraction, the surplus vibrational energy in the OH- reactant mostly remains in the H2O product owing to the prevailing dominance of the direct stripping mechanism.

Combining Force Fields and Neural Networks for an Accurate Representation of Bonded Interactions.

We present a formalism of a neural network encoding bonded interactions in molecules. This intramolecular encoding is consistent with the models of intermolecular interactions previously designed by this group. Variants of the encoding fed into a corresponding neural network may be used to economically improve the representation of torsional degrees of freedom in any force field. We test the accuracy of the reproduction of the ab initio potential energy surface on a set of conformations of two dipeptides, methyl-capped ALA and ASP, in several scenarios. The encoding, either alone or in conjunction with an analytical potential, improves agreement with ab initio energies that are on par with those of other neural network-based potentials. Using the encoding and neural nets in tandem with an analytical model places the agreements firmly within "chemical accuracy" of ±0.5 kcal/mol.

Inelastic scattering of NO(A2Σ+) + CO2: rotation-rotation pair-correlated differential cross sections.

A crossed beam velocity-map ion-imaging apparatus has been used to determine differential cross sections (DCSs) for the rotationally inelastic scattering of NO(A2Σ+, v = 0, j = 0.5) with CO2, as a function of both NO(A, v = 0, N') final state and the coincident final rotational energy of the CO2. The DCSs are dominated by forward-peaked scattering for all N', with significant rotational excitation of CO2, and a small backward scattered peak is also observed for all final N'. However, no rotational rainbow scattering is observed and there is no evidence for significant product rotational angular momentum polarization. New ab initio potential energy surface calculations at the PNO-CCSD(T)-F12b level of theory report strong attractive forces at long ranges with significant anisotropy relative to both NO and CO2. The absence of rotational rainbow scattering is consistent with removal of low-impact-parameter collisions via electronic quenching, in agreement with the literature quenching rates of NO(A) by CO2 and recent electronic structure calculations. We propose that high-impact-parameter collisions, that do not lead to quenching, experience strong anisotropic attractive forces that lead to significant rotational excitation in both NO and CO2, depolarizing product angular momentum while leading to forward and backward glory scattering.

Vibrational energy transfer in ammonia-helium collisions.

While the rotational energy transfer of ammonia by rare gas atoms and hydrogen molecules has been the focus of many studies, little is known about its vibrational relaxation, even though transitions involving the umbrella bending mode have been observed in many astrophysical environments. Here we explore the vibrational relaxation of the umbrella mode of ammonia induced by collisions with helium atoms by means of the close-coupling method on an ab initio potential energy surface. We compute cross sections up to kinetic energies of 1500 cm-1 and rate coefficients up to a temperature of 300 K for vibrational, rotational, and inversion transitions involving the lowest two vibrational states. We show that vibrational relaxation is much less efficient than rotation-inversion relaxation, although the rate coefficients for vibrational relaxation strongly increase with the temperature. We also observe important differences for vibrationally-elastic transitions within the lowest two vibrational states, i.e., for rotation-inversion transitions. These are a direct consequence of the difference in the tunnelling splitting of the lowest inversion levels.

Quantum rotational dynamics of l-C4(3Σ-g) by H2 at low temperatures employing a machine learning augmented potential energy surface.

A new four dimensional (4D) ab initio potential energy surface (PES) is generated for the collision of C4(3Σg-) with H2(1Σg), considering both molecules as rigid rotors. A supervised neural network model is created to augment the ab initio PES and to get the missing data points. Furthermore, space fixed expansion of the augmented PES is carried out using a least squares fit over two spherical harmonics terms, resulting in radial coefficients (λ1, λ2, and λ). The centre of symmetry in both C4 and H2 forces λ1 and λ2 to have even values, respectively. Moreover, the rotational states of C4 are only populated by odd levels due to its ground state triplet symmetry and the nuclear spin (I = 0) of 12C. The cross-sections and rate coefficients with para and ortho H2 partners are studied for various odd state transitions, where the rate coefficients of the ortho are 10-20% higher than those of the latter. The de-excitation rates obtained by the para H2 collisions are also compared to those of He and are found to be ∼1.7-2.8 times the He rates, across various order transitions. The simple scaling of He rates using a factor of 1.38 proves insufficient to describe para H2 rates. Therefore, these results show the importance of explicitly studying H2 as an important colliding partner, governing the kinetics of various rotational processes in the interstellar space.

Detailed quasiclassical dynamics of the F- + SiH3Cl multi-channel reaction.

We report a detailed quasiclassical trajectory study on the F- + SiH3Cl multi-channel reaction using a full-dimensional ab initio analytical potential energy surface. Reaction probabilities, cross sections, initial attack and scattering angle distributions as well as product relative translational, internal, vibrational, and rotational energy distributions are obtained in the collision energy range of 1-40 kcal mol-1 for the following channels: SiH3F + Cl-, SiH2Cl- + HF, SiH2F- + HCl, SiH2FCl + H-, SiH2 + FHCl-, and SiHFCl- + H2. All the channels are translationally cold indicating indirect mechanisms, except proton transfer (SiH2Cl- + HF), which shows mixed direct-indirect character. The angular distributions vary depending on collision energy and inversion/retention for SiH3F + Cl-. In the case of SiH2Cl- + HF front-side/back-side attack backward-forward/forward scattering preference is found at low/high collision energy. SiH2F- + HCl is formed with isotropic scattering and the preferred angle of attack is similar to the SiH3F + Cl- channel. SiH2FCl + H-/SiH2 + FHCl- favors back-side attack and isotropic/backward scattering, whereas SiHFCl- + H2 does not show significant angular preference.

Methane dimer rovibrational states and Raman transition moments.

Benchmark-quality rovibrational data are reported for the methane dimer from variational nuclear motion computations using an ab initio intermolecular potential energy surface reported by [M. P. Metz et al., Phys. Chem. Chem. Phys., 2019, 21, 13504-13525]. A simple polarizability model is used to compute Raman transition moments that may be relevant for future direct observation of the intermolecular dynamics. Non-negligible ΔK ≠ 0 transition moments arise in this symmetric top system due to strong rovibrational couplings.