Accurate global fit of the H4potential energy surface

New Full-Dimensional Reactive Potential Energy Surface for the H4 System

As the most abundant molecule in the universe, collisions involving H2 have important implications in astrochemistry. Collisions between hydrogen molecules also represent a prototype for assessing various dynamic methods for understanding fundamental few-body processes. In this work, we develop a new and highly accurate full-dimensional potential energy surface (PES) covering all reactive channels of the H2 + H2 system, which extends our previously reported H2 + H2 nonreactive PES [J. Chem. Theory Comput., 2021, 17, 6747] by adding 39,538 additional ab initio points calculated at the MRCI/AV5Z level in the reactive channels. The global PES is represented with high fidelity (RMSE = 0.6 meV for a total of 79,000 points) by a permutation invariant polynomial neural network (PIP-NN) and is suitable for studying collision-induced dissociation, single-exchange, as well as four-center exchange reactions. Preliminary quasi-classical trajectory studies on the new PIP-NN PES reveal strong vibrational enhancement of all reaction channels.

Recent advances in quantum theory on ro-vibrationally inelastic scattering

Molecular collisions are of fundamental importance in understanding intermolecular interaction and dynamics. Its importance is accentuated in cold and ultra-cold collisions because of the dominant quantum mechanical nature of the scattering. We review recent advances in the time-independent approach to quantum mechanical characterization of non-reactive scattering in tetratomic systems, which is ideally suited for large collisional de Broglie wavelengths characteristic in cold and ultracold conditions. We discuss quantum scattering algorithms between two diatoms and between a triatom and an atom and their implementation, as well as various approximate schemes. They not only enable the characterization of collision dynamics in realistic systems but also serve as benchmarks for developing more approximate methods.

Neural network potential energy surface for the low temperature ring polymer molecular dynamics of the H2CO + OH reaction

A new potential energy surface (PES) and dynamical study of the reactive process of H₂CO + OH toward the formation of HCO + H₂O and HCOOH + H are presented. In this work, a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, which prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many-body term has been used, fitted with a novel artificial neural network scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory (QCT) and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated with excitations of the normal modes of the ring polymer.

Zero- and high-pressure mechanisms in the complex forming reactions of OH with methanol and formaldehyde at low temperatures

A recent Ring Polymer Molecular Dynamics study of the reactions of OH with methanol and formaldehyde, at zero pressure and below 100 K, has shown the formation of long lived complexes, with long lifetimes, longer than 100 ns for the lower temperatures studied, 20-100 K (del Mazo-Sevillano et al., 2019). These long lifetimes support the existence of multi collision events with the He buffer-gas atoms under experimental conditions, as suggested by several transition state theory studies of these reactions. In this work we study these secondary collisions, as a dynamical approach to study pressure effects on these reactions. For this purpose, the potential energy surfaces of He with H2CO, OH, H2O and HCO are calculated at highly accurate ab initio level. The stability of some of the complexes is studied using Path Integral Molecular dynamics techniques, determining that OH-H2CO complexes can be formed up to 100 K or higher temperatures, while the weaker He-H2CO complexes dissociate at approximately 50 K. The predicted IR intensity spectra shows new features which could help the identification of the OH-H2CO complex. Finally, the He-H2CO + OH and OH-H2CO + He collisions are studied using quassi-classical trajectories, finding that the cross section to produce HCO + H2O products increases with decreasing collision energy, and that it is ten times higher in the He-H2CO + OH case.

Non-adiabatic Quantum Dynamics of the Dissociative Charge Transfer He++H2 → He+H+H. Front Chem 2019;7:249. [PMID: 31041310 PMCID: PMC6477054 DOI: 10.3389/fchem.2019.00249]

We present the non-adiabatic, conical-intersection quantum dynamics of the title collision where reactants and products are in the ground electronic states. Initial-state-resolved reaction probabilities, total integral cross sections, and rate constants of two H₂ vibrational states, v₀ = 0 and 1, in the ground rotational state (j₀ = 0) are obtained at collision energies E_coll ≤ 3 eV. We employ the lowest two excited diabatic electronic states of HeH2+ and their electronic coupling, a coupled-channel time-dependent real wavepacket method, and a flux analysis. Both probabilities and cross sections present a few groups of resonances at low E_coll, whose amplitudes decrease with the energy, due to an ion-induced dipole interaction in the entrance channel. At higher E_coll, reaction probabilities and cross sections increase monotonically up to 3 eV, remaining however quite small. When H₂ is in the v₀ = 1 state, the reactivity increases by ~2 orders of magnitude at the lowest energies and by ~1 order at the highest ones. Initial-state resolved rate constants at room temperature are equal to 1.74 × 10⁻¹⁴ and to 1.98 × 10⁻¹² cm³s⁻¹ at v₀ = 0 and 1, respectively. Test calculations for H₂ at j₀ = 1 show that the probabilities can be enhanced by a factor of ~1/3, that is ortho-H₂ seems ~4 times more reactive than para-H₂.

Low temperature reaction dynamics for CH3OH + OH collisions on a new full dimensional potential energy surface

Is the rise of the rate constant measured in laval expansion experiments of OH with organic molecules at low temperatures due to the reaction between the reactants or due to the formation of complexes with the buffer gas? This question has importance for understanding the evolution of prebiotic molecules observed in different astrophysical objects. Among these molecules methanol is one of the most widely observed, and its reaction with OH has been studied by several groups showing a fast increase in the rate constant under 100 K. Transition state theory doesn't reproduce this behavior and here dynamical calculations are performed on a new full dimensional potential energy surface developed for this purpose. The calculated classical reactive cross sections show an increase at low collision energies due to a complex forming mechanism. However, the calculated rate constant at temperatures below 100 K remains lower than the observed one. Quantum effects are likely responsible for the measured behavior at low temperatures.

Full dimensional potential energy surface and low temperature dynamics of the H2CO + OH → HCO + H2O reaction

A new method is proposed to analytically represent the potential energy surface of reactions involving polyatomic molecules capable of accurately describing long-range interactions and saddle points, needed to describe low-temperature collisions. It is based on two terms, a reactive force field term and a many-body term. The reactive force field term accurately describes the fragments, long-range interactions among them and the saddle points for reactions. The many-body term increases the desired accuracy everywhere else. This method has been applied to the OH + H2CO → H2O + HCO reaction, giving a barrier of 27.4 meV. The simulated classical rate constants with this potential are in good agreement with recent experimental results [Ocaña et al., Astrophys. J., 2017, submitted], showing an important increase at temperatures below 100 K. The reaction mechanism is analyzed in detail here, and explains the observed behavior at low energy by the formation of long-lived collision complexes, with roaming trajectories, with a capture observed for very long impact parameters, >100 a.u., determined by the long-range dipole-dipole interaction.

Non-adiabatic couplings and dynamics in proton transfer reactions of Hn (+) systems: Application to H2+H2 (+)→H+H3 (+) collisions

Analytical derivatives and non-adiabatic coupling matrix elements are derived for Hn (+) systems (n = 3-5). The method uses a generalized Hellmann-Feynman theorem applied to a multi-state description based on diatomics-in-molecules (for H3 (+)) or triatomics-in-molecules (for H4 (+) and H5 (+)) formalisms, corrected with a permutationally invariant many-body term to get high accuracy. The analytical non-adiabatic coupling matrix elements are compared with ab initio calculations performed at multi-reference configuration interaction level. These magnitudes are used to calculate H2(v(')=0,j(')=0)+H2 (+)(v,j=0) collisions, to determine the effect of electronic transitions using a molecular dynamics method with electronic transitions. Cross sections for several initial vibrational states of H2 (+) are calculated and compared with the available experimental data, yielding an excellent agreement. The effect of vibrational excitation of H2 (+) reactant and its relation with non-adiabatic processes are discussed. Also, the behavior at low collisional energies, in the 1 meV-0.1 eV interval, of interest in astrophysical environments, is discussed in terms of the long range behaviour of the interaction potential which is properly described within the triatomics-in-molecules formalism.

A high-level ab initio study of the N2 + N2 reaction channel

A new six-dimensional (6D) global potential energy surface (PES) is proposed for the full range description of the interaction of the N2(1Σg+)+N2(1Σg+) system governing collisional processes, including N atom exchange. The related potential energy values were determined using high-level ab initio methods. The calculations were performed at a coupled-cluster with single and double and perturbative triple excitations level of theory in order to have a first full range picture of the PES. Subsequently, in order to accurately describe the stretching of the bonds of the two interacting N2 molecules by releasing the constraints of being considered as rigid rotors, for the same molecular geometries higher level of theory multi reference calculations were performed. Out of the calculated values a 6D 4-atoms global PES was produced for use in dynamical calculations. The ab initio calculations were made possible by the combined use of High Throughput Computing and High Performance Computing techniques within the frame of a computing grid empowered molecular simulator.

Full-dimensional time-dependent wave packet dynamics of H2 + D2 reaction

Collision induced dissociation (CID), four center reaction (4C), and single exchange reaction (SE) in H(2) (v(1) = high) + D(2) (v(2) = low) were studied by means of time-dependent wave packet approach within a full-dimensional model. Initial state-selected total reaction probabilities for the three competitive processes have been computed on two realistic global potential energy surfaces of Aguado-Suárez-Paniagua and Boothroyd-Martin-Keogh-Peterson (BMKP) with the total angular momentum J = 0. The role of both vibrationally excited and rotationally excited reagents was examined by varying the initial vibrational and rotational states. The vibrational excitation of the hot diatom gives an enhancement effect on the CID process, while the vibrational excitation of the cold diatom gives an inhibition effect. The rotational excitation of both reagents has a significant effect on the reaction process. The 4C and SE probabilities are at least one order of magnitude smaller than the CID probabilities over the energy range considered. Isotope substitution effects were also studied by substituting the collider D(2) by H(2) and HD on the BMKP potential energy surfaces. The CID process is most efficient for the H(2) + D(2) combination and least efficient for the H(2) + H(2) combination and is different for the 4C and SE processes.

Interpretation of the vibrational relaxation of H2 in H2 within the semiclassical effective mass approach

The temperature dependence of the rate coefficients for vibrational relaxation of H2 in neat H2 is interpreted within the semiclassical effective mass approach. Across the temperature range of 80-3000 K, the experimental rate coefficients vary by five orders of magnitude and fall onto a strongly nonlinear Landau-Teller plot. This behavior is explained by the nonclassical nature of the energy release and by a substantial participation of rotation of the colliding partners in inducing the vibrational transition. A single fitting parameter, the optimal reduced mass, permits one to represent the temperature dependence of the rate coefficient within a factor of 2. This parameter is found to be close to that obtained from a simple model suggested by Sewell et al. [J. Chem. Phys. 99, 2567 (1993)].

Quasiclassical trajectory study of reactive and dissociative processes in H2+H2: Comparison with quantum-mechanical calculations

Quasiclassical trajectory calculations have been carried out for H(2)(v(1)=high)+H(2)(v(2)=low) collisions within a three degrees of freedom model where five different geometries of the colliding complex were considered. Within this approach, probabilities for different competitive processes are studied: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. The purpose is to compare in detail with equivalent quantum-mechanical wave packet calculations [Bartolomei et al., J. Chem. Phys 122, 064305 (2005)], especially the behavior of the probabilities near reaction thresholds. Quasiclassical calculations compare quite well with the quantum-mechanical ones for collision induced dissociation as well as for the four center reaction, although quantum effects become very important near thresholds, particularly for lower v(1)'s and for the four center process. Less quantitative agreement is found for reactive dissociation and three-body complex formation. It is found that most quantum effects are due to differences between quantum and classical vibrational distributions of H(2)(v(1)=high). Zero point energy violation has been found in the classical reactive-dissociative probabilities. Extension of these findings to full-dimensional treatments is examined.

A quasiclassical trajectory study of collisional energy transfer and dissociation in He + H2(v,j) using a new potential energy surface

Quasiclassical trajectories for He + H2 were carried out using the recent ab initio potential of Boothroyd, Martin, and Peterson (J. Chem. Phys. 2003, 119, 3187) and results for the 348 (v, j) states of H2 are compared to those of earlier calculations that used the potential of Wilson, Kapral, and Burns (Chem. Phys. Lett. 1974, 24, 4884). Examined are the cross sections for energy transfer and dissociation, the extent of threshold elevation, and the interconversion of vibrational and rotational energy. Implications for modeling the interstellar medium are discussed.

Simple three-parameter model potential for diatomic systems: from weakly and strongly bound molecules to metastable molecular ions

Based on a simplest molecular-orbital theory of H(2)(+), a three-parameter model potential function is proposed to describe ground-state diatomic systems with closed-shell and/or S-type valence-shell constituents over a significantly wide range of internuclear distances. More than 200 weakly and strongly bound diatomics have been studied, including neutral and singly charged diatomics (e.g., H(2), Li(2), LiH, Cd(2), Na(2)(+), and RbH(-)), long-range bound diatomics (e.g., NaAr, CdNe, He(2), CaHe, SrHe, and BaHe), metastable molecular dications (e.g., BeH(++), AlH(++), Mg(2)(++), and LiBa(++)), and molecular trications (e.g., YHe(+++) and ScHe(+++)).

Wave packet dynamics of H2(v1=8-14)+H2(v2=0-2): the role of the potential energy surface on different reactive and dissociative processes

A time-dependent wave packet method has been used to study different competing products of H(2)+H(2) collisions: four center reaction, collision induced dissociation, reactive dissociation, and three-body complex formation. A three-degree-of-freedom reduced dimensionality model has been used for five different geometries of the colliding complex (parallel H, crossed X, collinear L, and two T-shaped geometries T(I) and T(II)), with reactants in selected vibrational states with one diatom vibrationally "hot" and the other one vibrationally "cold." Product probabilities have been calculated using two potential energy surfaces [J. Chem. Phys. 101, 4004 (1994); J. Chem. Phys. 116, 666 (2002)] in order to compare their performance in the dynamics. The regions of the potential energy surfaces responsible of the threshold behavior of the probabilities have been identified. Overall, we have found that the most recent potential energy surface is less anisotropic, provides a smaller propensity for insertion-type processes, and gives lower energy thresholds.

Inelastic collisions of molecular hydrogen: A comparison of results from quantum and classical mechanics

Full-dimensional quantum and classical calculations have been carried out for inelastic (nonreactive) energy transfer in H2+H2 on the ab initio potential energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)]. State-to-state cross sections are determined and compared for transitions from H2(0,j(ab))+H2(1,j(cd)). While there is excellent agreement for transitions involving small Deltaj, for larger Deltaj and for vibrational relaxation, significant differences are observed which exhibit no systematic trends. Reasons for this disagreement are discussed.