Search

The kinetic energy dependence of the title reaction is examined using guided ion beam tandem mass spectrometry. Because this reaction is spin-forbidden, crossings between octet and sextet hypersurfaces presumably must occur. Furthermore, Sm+ must transition from a 4f66s1 configuration in the reactant to 4f55d2 in order to have the orbital occupancy required to form the triple bond in SmO+ (6Δ). Despite being strongly exothermic (∼4 eV), the reaction proceeds with low efficiency (18% ± 4%) via a barrierless process at low energies. Below ∼0.3 eV, the cross section follows a kinetic energy dependence that roughly parallels that of the collision rate for ion–dipole reactions. At higher collision energies, the reaction cross section increases until it follows the trajectory cross section closely from 3 to 5 eV, indicating that another pathway opens on the reaction hypersurface. Modeling this increase yields a threshold energy for this new pathway at 0.54 ± 0.05 eV. Theoretical potential energy surfaces that do not include spin–orbit interactions for the reaction show that there is a barrier of height 1.19 eV (MP2) or 0.49 eV [CCSD(T)] to insertion of Sm+ into the N2–O bond and that there are several places where octet and sextet surfaces can intersect and interact. By considering the distribution of spin–orbit states generated in the ion source, the internal energy of the N2O reactant, and the influence of coupling between electronic, orbital, and rotational angular momentum, the low-efficiency, exothermic behavior as well as the increase in efficiency at higher energies can plausibly be explained. [ABSTRACT FROM AUTHOR]

One of the most reactive intermediates for oxidative reactions is the oxyl radical, an electron-deficient oxygen atom. The discovery of a new vibration upon photoexcitation of the oxygen evolution catalysis detected the oxyl radical at the SrTiO3 surface. The vibration was assigned to a motion of the sub-surface oxygen underneath the titanium oxyl (Ti–O●−) created upon hole transfer to (or electron extraction from) a hydroxylated surface site. Evidence for such an interfacial mode is derived from its spectral shape, which exhibited a Fano resonance—a coupling of a sharp normal mode to continuum excitations. Here, this Fano resonance is utilized to derive precise formation kinetics of the oxyl radical and its associated potential energy surface (PES). From the Fano lineshape, the formation kinetics are obtained from the anti-resonance (the kinetics of the coupling factor), the resonance (the kinetics of the coupled continuum excitations), and the frequency integrated spectrum (the kinetics of the normal mode's cross-section). All three perspectives yield logistic function growth with a half-rise of 2.3 ± 0.3 ps and a time constant of 0.48 ± 0.09 ps. A non-equilibrium transient associated with photoexcitation is separated from the rise of the equilibrated PES. The logistic function characterizes the oxyl coverage at the very initial stages (t ∼ 0) to have an exponential growth rate that quickly decreases toward zero as a limiting coverage is reached. Such time-dependent reaction kinetics identify a dynamic activation barrier associated with the formation of a PES and quantify it for oxyl radical coverage. [ABSTRACT FROM AUTHOR]

A new multiple-time scale integration method is presented that propagates ab initio path integral molecular dynamics (PIMD). This method uses a large time step to generate an approximate geometrical configuration whose energy and gradient are evaluated at the level of an ab initio method, and then, a more precise integration scheme, e.g., the Bulirsch–Stoer method or velocity Verlet integration with a smaller time step, is used to integrate from the previous step using the computationally efficient interpolated potential energy surface constructed from two consecutive points. This method makes the integration of PIMD more efficient and accurate compared with the velocity Verlet integration. A Nos e ́ –Hoover chain thermostat combined with this new multiple-time scale method has good energy conservation even with a large time step, which is usually challenging in velocity Verlet integration for PIMD due to the very small chain mass when a large number of beads are used. The new method is used to calculate infrared spectra and free energy profiles to demonstrate its accuracy and capabilities. [ABSTRACT FROM AUTHOR]

A merged potential energy surface (PES) is introduced for CO + CO collisions by combining a recent full-dimensional ab initio PES [Chen et al. J. Chem. Phys. 153, 054310 (2020)] and analytical long-range multipolar interactions. This merged PES offers a double advantage: it retains the precision of the ab initio PES in describing the van der Waals well and repulsive short range while providing an accurate physical description of long-range interaction; it significantly reduces the computational time required for trajectory integration since the long-range portion of the ab initio PES (involving numerous neural network fitting parameters) is now replaced by the analytical model potential. Based on the present merged PES, mixed Quantum-Classical (MQC) calculations, which capture quantum effects related to vibrational motion, align with a range of experimental data, including transport properties, vibrational energy transfer between CO and its isotoplogues, as well as rate coefficients for V–V and V–T/R processes. Notably, the original ab initio PES yields V–T/R rate coefficients at low temperatures that are significantly higher than the experimental data due to the artificial contribution of its unphysical long-range potential. In addition to conducting extensive MQC calculations to obtain raw data for V–V and V–T/R rate coefficients, we employ Gaussian process regression to predict processes lacking computed MQC data, thereby completing the considered V–V and V–T/R datasets. These extensive rate coefficient datasets, particularly for V–T/R processes, are unprecedented and reveal the significant role played by V–T/R processes at high temperatures, emphasizing the necessity of incorporating both V–V and V–T/R processes in the applications. [ABSTRACT FROM AUTHOR]

This work is an extensive investigation of the F + CH3NH2 reaction dynamics using a newly-developed potential energy surface (PES). The full-dimensional spin–orbit (SO) corrected (MRCI+Q/aug-cc-pwCVDZ) PES is developed by the Robosurfer program package and the ManyHF method is used in order to fix the Hartree–Fock (HF) convergence issues in the entrance channel. On the surface, retrieved by the fitting of the iteratively extended set of the ManyHF-CCSD(T)-F12a/triple-zeta-quality and SO-corrected energy points, quasi-classical trajectory (QCT) simulations are run. By analyzing the opacity functions and integral cross sections (ICSs) of six reaction channels, the dynamics of the two most reactive hydrogen-abstraction reactions resulting in HF + CH2NH2/CH3NH products are selected for a thorough examination. Despite the statistically and thermodynamically expected results, the kinetically preferred amino hydrogen-abstraction is the dominant mechanism at low collision energies. The initial attack angle and scattering angle distributions are investigated as well. The post-reaction energy distributions show that the collision energy mostly converts into the translational energy of the products, while the reaction energy excites the vibration of the products. The computed vibrationally resolved rotational distributions and vibrational state distributions of the HF product are compared to experimental data, and the theory and experiment are found to be in good agreement. [ABSTRACT FROM AUTHOR]

In this paper, a global and full-dimensional potential energy surface at the 2A″ ground state for the Al + O2 → AlO + O reaction was constructed, for the first time, based on extensive electronic structure calculations using the doubly hybrid functional XYG3 and potential energy surface fittings by neural networks. Details of the reaction paths have been analyzed. It was found that both two intermediates, the cyclic-AlO2 and the linear-OAlO, were able to dissociate to the AlO + O products, and the isomerization process between these two intermediates was controlled by conical intersections between two 2A″ states. Ro-vibrational state resolved integral cross sections have also been calculated at collision energies from 1.0 to 10.0 kcal/mol. The results support the harpooning mechanism in this metal-oxidant-involved reaction. [ABSTRACT FROM AUTHOR]

An accurate potential energy surface (PES) for the lowest lying A ′ ′ 4 state of the CNO system is presented based on explicitly correlated multi-reference configuration interaction calculations with quadruple zeta basis set (MRCI-F12/cc-pVQZ-F12). The ab initio energies are fitted using the double many-body expansion method, thus incorporating long-range energy terms that can accurately describe the electrostatic and dispersion interactions with physically motivated decaying functions. Together with the previously fitted lowest A ′ 2 and A ′ ′ 2 states using the same theoretical framework, this constitutes a new set of PESs that are suitable to predict rate coefficients for all atom–diatom reactions of the CNO system. We use this set of PESs to calculate thermal rate coefficients for the C( P 3 ) + NO ( Π 2 ) reaction and compare the temperature dependence and product branching ratios with experimental results. The comparison between theory and experiment is shown to be improved over previous theoretical studies. We highlight the importance of the long-range interactions for low-temperature rate coefficients. [ABSTRACT FROM AUTHOR]

Symmetry, in particular permutational symmetry, of a potential energy surface (PES) is a useful property in quantum chemical calculations. It facilitates, in particular, state labelling and identification of degenerate states. In many practically important applications, however, these issues are unimportant. The imposition of exact symmetry and the perception that it is necessary create additional methodological requirements narrowing or complicating algorithmic choices that are thereby biased against methods and codes that by default do not incorporate symmetry, including most off-the-shelf machine learning methods that cannot be directly used if exact symmetry is demanded. By introducing symmetric and unsymmetric errors into the PES of H2CO in a controlled way and computing the vibrational spectrum with collocation using symmetric and nonsymmetric collocation point sets, we show that when the deviations from an ideal PES are random, imposition of exact symmetry does not bring any practical advantages. Moreover, a calculation ignoring symmetry may be more accurate. We also compare machine-learned PESs with and without symmetrization and demonstrate that there is no advantage of imposing exact symmetry for the accuracy of the vibrational spectrum. [ABSTRACT FROM AUTHOR]

In this work the H2O–HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O–HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O–HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O–HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former. [ABSTRACT FROM AUTHOR]

New six-dimensional ab initio potential energy surfaces (PESs) for the N2–CO2 complex, which involve the stretching vibration of N2 and the Q3 normal mode for the ν3 asymmetric stretching vibration of CO2, were constructed using the CCSD(T)-F12/AVTZ method with midpoint bond functions. Two vibrational averaged 4D interaction potentials were obtained by integrating over the two intramolecular coordinates. It was found that both PESs possess two equivalent T-shaped global minima as well as two in-plane and one out-of-plane saddle points. Based on these PESs, rovibrational bound states and energy levels were calculated applying the radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm. The splitting of the energy levels between oN2–CO2 and pN2–CO2 for the intermolecular vibrational ground state is determined to be only 0.000 09 cm−1 due to the higher barriers. The obtained band origin shift is about +0.471 74 cm−1 in the N2–CO2 infrared spectra with CO2 at the ν3 zone, which coincides with the experimental data of +0.483 74 cm−1. The frequencies of the in-plane geared-bending for N2–CO2 at the ν3 = 0 and 1 states of CO2 turn out to be 21.6152 and 21.4522 cm−1, the latter reproduces the available experimental 21.3793 cm−1 value with CO2 at the ν3 zone. The spectral parameters fitted from the rovibrational energy levels show that this dimer is a near prolate asymmetric rotor. The computed microwave transitions as well as the infrared fundamental and combination bands for the complex agree well with the observed data. [ABSTRACT FROM AUTHOR]

An accurate interatomic potential energy surface for ground state of SH2− is reported at the multireference configuration interaction level of theory using aug–cc–pVQZ basis set. The calculated ab initio energy points are utilized, after semi-empirically correcting the dynamical correlation by the many-body expansion-scaled external correlation method. To verify the accuracy of the potential energy surface, molecular reaction dynamics, including the relationship between reaction probability and collision energy, thermal rate constant, and integral scattering cross section are calculated by using the quantum time-dependent wave packet method. The final results show that our potential energy surface is accurate and the molecular reaction dynamics results are satisfactory, and the accurate potential energy surface and its dynamics results can also lay a foundation for the future experimental research of SH−-containing molecules. This study provides a reasonable ionic chemical dissociation channel for the dynamics, A global many-body expansion (MEB) PES is reported by fitting high-precision ab initio energies calculated at the multireference configuration interaction level with Davidson correction (MRCI(Q)) and the aug-cc-pVQZ basis set (AVQZ). By comparing with the integral cross section with experimental result, which have agreement with the previous experimental and theoretical values. Our theoretical results can provide effective support for future dynamic experiments. [ABSTRACT FROM AUTHOR]

A full-dimensional spin–orbit (SO)-corrected potential energy surface (PES) is developed for the Cl + CH3NH2 multi-channel system. Using the new PES, a comprehensive reaction dynamics investigation is performed for the most reactive hydrogen-abstraction reactions forming HCl + CH2NH2/CH3NH. Hartree–Fock (HF) convergence problems in the reactant region are handled by the ManyHF method, which finds the lowest-energy HF solution considering several different initial guess orbitals. The PES development is carried out with the Robosurfer program package, which iteratively improves the surface. Energy points are computed at the ManyHF-UCCSD(T)-F12a/cc-pVDZ-F12 level of theory combined with basis set (ManyHF-RMP2-F12/cc-pVTZ-F12 – ManyHF-RMP2-F12/cc-pVDZ-F12) and SO (MRCI+Q/aug-cc-pwCVDZ) corrections. Quasi-classical trajectory simulations show that the CH3-side hydrogen abstraction occurs more frequently in contrast to the NH2-side reaction. In both cases, the integral cross sections decrease with increasing collision energy (Ecoll). A reaction mechanism shifting from indirect to direct stripping can be observed from the opacity functions, scattering angle, and translation energy distributions as Ecoll increases. Initial attack angle distributions reveal that chlorine prefers to abstract hydrogen from the approached functional group. The collision-energy dependence of the product energy distributions shows that the initial translational energy mainly transfers to product recoil. The HCl vibrational and rotational energy values are comparable and nearly independent of collision energy, while the CH2NH2 and CH3NH co-products' vibrational energy values are higher than the rotational energy values with more significant Ecoll dependence. The HCl(v = 0) rotational distributions are compared with experiment, setting the direction for future investigations. [ABSTRACT FROM AUTHOR]

The remote sensing of abundance and properties of HCl—the main atmospheric reservoir of Cl atoms that directly participate in ozone depletion—is important for monitoring the partitioning of chlorine between "ozone-depleting" and "reservoir" species. Such remote studies require knowledge of the shapes of molecular resonances of HCl, which are perturbed by collisions with the molecules of the surrounding air. In this work, we report the first fully quantum calculations of collisional perturbations of the shape of a pure rotational line in H35Cl perturbed by an air-relevant molecule [as the first model system we choose the R(0) line in HCl perturbed by O2]. The calculations are performed on our new highly accurate HCl(X1Σ+)–O2( X 3 Σ g − ) potential energy surface. In addition to pressure broadening and shift, we also determine their speed dependencies and the complex Dicke parameter. This gives important input to the community discussion on the physical meaning of the complex Dicke parameter and its relevance for atmospheric spectra (previously, the complex Dicke parameter for such systems was mainly determined from phenomenological fits to experimental spectra and the physical meaning of its value in that context is questionable). We also calculate the temperature dependence of the line shape parameters and obtain agreement with the available experimental data. We estimate the total combined uncertainties of our calculations at 2% relative root-mean-square error in the simulated line shape at 296 K. This result constitutes an important step toward computational population of spectroscopic databases with accurate ab initio line shape parameters for molecular systems of terrestrial atmospheric importance. [ABSTRACT FROM AUTHOR]

The hydrogen abstraction reaction of the cyano radical with molecules of ethane presents some interesting points in the chemistry from ultra-cold to combustion environments especially with regard to HCN(v) product vibrational distribution. In order to understand its dynamics, a new analytical full-dimensional potential energy surface was developed, named PES-2023. It uses a combination of valence bond and mechanic molecular terms as the functional form, fitted to high-level ab initio calculations at the explicitly correlated CCSD(T)-F12/aug-cc-pVTZ level on a reduced and selected number of points describing the reactive process. The new surface showed a continuous and smooth behavior, describing reasonably the topology of the reaction: high exothermicity, low barrier, and presence of intermediate complexes in the entrance and exit channels. Using quasi-classical trajectory calculations (QCT) on the new PES-2023, a dynamics study was performed at room temperature with special emphasis on the HCN(v1,v2,v3) product stretching and bending vibrational excitations, and the results were compared with the experimental evidence, which presented discrepancies in the bending excitation. The available energy was mostly deposited as HCN(v) vibrational energy with the vibrational population inverted in the CH stretching mode and not inverted in the CN stretching and bending modes, thus simulating the experimental evidence. Other dynamics properties at room temperature were also analyzed; cold rotational energy distribution was found, associated with a linear and soft transition state, and backward scattering distribution was found, associated with a rebound mechanism. [ABSTRACT FROM AUTHOR]

The potential energy surface (PES) describing the interaction of the methanol molecule with molecular hydrogen has been calculated by the use of the explicitly correlated coupled cluster method, including single, double, and (perturbative) triple excitations [CCSD(T)-F12a] and a correlation-consistent aug-cc-pVTZ basis, with the assumption of fixed molecular geometries. The computed points were fit to a functional form appropriate for time-independent quantum scattering calculations of rotationally inelastic cross sections and rate coefficients. Stationary points on the PES were located, and the global minimum was found to have an energy equal to −254.7 cm−1 relative to the energy of the separated molecules. This PES was used in time-independent close coupling quantum scattering calculations to determine state-to-state cross sections and rate coefficients for rotational transitions within the A- and E-type nuclear spin torsional ground states. [ABSTRACT FROM AUTHOR]

The rotational excitation rates of NCCN species are studied for its collision with hydrogen (H2) in temperatures ranging from 1 to 100 K. Such collisions can occur in the interstellar medium with H2 in either para (p-) or ortho (o-) state, of which the p-H2 state can be approximated via its collision with He (using a scaling factor) or with a reduced rigid rotor-H2 surface (by averaging over various orientations of H2). In the current work, a four-dimensional (4D) ab initio potential energy surface (PES) is considered to study the collision dynamics of H2 in both p- and o-states and the results are compared with previous approximations. The 4D surface is constructed using the explicitly correlated coupled-cluster method CCSD(T)-F12b with the augmented triple zeta basis AVTZ and then fitted into an artificial neural networks (NN) model to augment the surface and account for missing data points. The radial coefficients are obtained from this NN fitted 4D PES via a least square fit over two spherical harmonics functions. The cross sections (σ) are computed using the close-coupling (CC) method (until 230 cm−1) for both p- and o-H2 collisions, and the rates are obtained by Boltzmann distribution over the translational energy of H2 until 100 K. The o-H2 rates are found to be higher by 25%–30% and 10%–20% compared to the p-H2 rates for Δj = 2 and higher order transitions, respectively. The coupled-state/centrifugal sudden approximated rates are also computed and found to have deviations as large as 40% when compared to CC rates, thus making quantitative descriptions unreliable. [ABSTRACT FROM AUTHOR]

In this work, the dynamics of a prototypical heavy–light–heavy abstract reaction, Cl(2P) + HCl → HCl + Cl(2P), is investigated both by constructing a new potential energy surface (PES) and by rate coefficient calculations. Both the permutation invariant polynomial neural network method and the embedded atom neural network (EANN) method, based on ab initio MRCI-F12+Q/AVTZ level points, are used for obtaining globally accurate full-dimensional ground state PES, with the corresponding total root mean square error being only 0.043 and 0.056 kcal/mol, respectively. In addition, this is also the first application of the EANN in a gas-phase bimolecular reaction. The saddle point of this reaction system is confirmed to be nonlinear. In comparison with both the energetics and rate coefficients obtained on both PESs, we find that the EANN is reliable in dynamic calculations. A full-dimensional approximate quantum mechanical method, ring-polymer molecular dynamics with a Cayley propagator, is employed to obtain the thermal rate coefficients and kinetic isotopic effects of the title reaction Cl(2P) + XCl→ XCl + Cl(2P) (H, D, Mu) on both new PESs, and the kinetic isotope effect (KIE) is also obtained. The rate coefficients reproduce the experimental results at high temperatures perfectly but with moderate accuracy at lower temperatures, but the KIE is with high accuracy. The similar kinetic behavior is supported by quantum dynamics using wave packet calculations as well. [ABSTRACT FROM AUTHOR]

A precise analytical potential energy surface (PES) of HCO(X2A′) is fitted from a great quantity of ab initio energy points computed with the multi-reference configuration interaction method and aug-cc-pV(Q/5)Z basis sets. The whole energy points extrapolated to the complete basis set limit are fitted by the many-body expansion formula. The calculated topographic characteristics are analyzed and compared with the existing work to prove the precision of the present HCO(X2A′) PES. By utilizing the time-dependent wave packet and quasi-classical trajectory methods, the reaction probabilities, integral cross sections, and rate constants are computed. The results are compared in detail with the former results carried out on the other PES. Moreover, the provided information on stereodynamics leads to an in-depth understanding of the role of collision energy in product distribution. [ABSTRACT FROM AUTHOR]

We develop a full-dimensional analytical potential energy surface (PES) for the F− + SiH3Cl reaction using Robosurfer for automatically sampling the configuration space, the robust [CCSD-F12b + BCCD(T) − BCCD]/aug-cc-pVTZ composite level of theory for computing the energy points, and the permutationally invariant polynomial method for fitting. Evolution of the fitting error and the percentage of the unphysical trajectories are monitored as a function of the iteration steps/number of energy points and polynomial order. Quasi-classical trajectory simulations on the new PES reveal rich dynamics resulting in high-probability SN2 (SiH3F + Cl−) and proton-transfer (SiH2Cl− + HF) products as well as several lower-probability channels, such as SiH2F− + HCl, SiH2FCl + H−, SiH2 + FHCl−, SiHFCl− + H2, SiHF + H2 + Cl−, and SiH2 + HF + Cl−. The Walden-inversion and front-side-attack-retention SN2 pathways are found to be competitive, producing nearly racemic products at high collision energies. The detailed atomic-level mechanisms of the various reaction pathways and channels as well as the accuracy of the analytical PES are analyzed along representative trajectories. [ABSTRACT FROM AUTHOR]

We have developed a full-dimensional analytical ab initio potential energy surface (PES) for the Cl− + CH3I reaction using the Robosurfer program system. The energy points have been computed using a robust composite method defined as CCSD-F12b + BCCD(T) – BCCD with the aug-cc-pVTZ(-PP) basis set and have been fitted by the permutationally invariant polynomial approach. Quasi-classical trajectory simulations on the new PES reveal that two product channels are open in the collision energy (Ecoll) range of 1–80 kcal/mol, i.e., SN2 leading to I− + CH3Cl and iodine abstraction (above ∼45 kcal/mol) resulting in ICl− + CH3. Scattering angle, initial attack angle, product translational energy, and product internal energy distributions show that the SN2 reaction is indirect at low Ecoll and becomes direct-rebound–back-side (CH3-side) attack-type, as Ecoll increases. Iodine abstraction mainly proceeds with direct stripping mechanism with side-on/back-side attack preference. Comparison with crossed-beam experiments and previous direct dynamics simulations shows quantitative or qualitative agreement and also highlights possible theoretical and/or experimental issues motivating further research. [ABSTRACT FROM AUTHOR]

Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl 2H 2 potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M 2H 2 (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist. [ABSTRACT FROM AUTHOR]

A simple and intuitive symmetry-based method for deriving the vibration–rotation kinetic energy operator (in polyspherical coordinates) for a general polyatomic molecule using the kinetic energy of a triatomic molecule was considered. Two vectors define the molecule-fixed coordinate system. For example, those vectors can determine positions of two atoms with respect to the third atom. By rotating the coordinate system, it is possible to obtain expressions that relate coefficients determining the kinetic energy in two coordinate systems. Those expressions let us obtain formulas for kinetic energy of general polyatomic molecules from kinetic energy of the triatomic molecules. Vibrational Energy levels of COF2 molecules were calculated from new ab initio potential energy surface. The comparison of the calculated and experimental vibrational energy levels is also reported. [ABSTRACT FROM AUTHOR]

Transition state (TS) on the potential energy surface (PES) plays a key role in determining the kinetics and thermodynamics of chemical reactions. Inspired by the fact that the dynamics of complex systems are always driven by rare but significant transition events, we herein propose a TS search method in accordance with the Q‐learning algorithm. Appropriate reward functions are set for a given PES to optimize the reaction pathway through continuous trial and error, and then the TS can be obtained from the optimized reaction pathway. The validity of this Q‐learning method with reasonable settings of Q‐value table including actions, states, learning rate, greedy rate, discount rate, and so on, is exemplified in 2 two‐dimensional potential functions. In the applications of the Q‐learning method to two chemical reactions, it is demonstrated that the Q‐learning method can predict consistent TS and reaction pathway with those by ab initio calculations. Notably, the PES must be well prepared before using the Q‐learning method, and a coarse‐to‐fine PES scanning scheme is thus introduced to save the computational time while maintaining the accuracy of the Q‐learning prediction. This work offers a simple and reliable Q‐learning method to search for all possible TS and reaction pathway of a chemical reaction, which may be a new option for effectively exploring the PES in an extensive search manner. [ABSTRACT FROM AUTHOR]

Using as functional form a combination of valence bond and mechanic molecular terms a new full-dimensional potential energy surface was developed for the title reaction, named PES-2022, which was fitted to high-level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples-F12 explicitly correlated level on a representative number of points describing the reactive system. This surface simultaneously describes the two reaction channels, hydrogen abstraction from the methyl group [(R1) path] and from the alcohol group [(R2) path] of methanol to form water. PES-2022 is a smooth and continuous surface, which reasonably describes the topology of this reactive system from reactants to products, including the intermediate complexes present in the system. Based on PES-2022 an exhaustive dynamics study was performed using quasi-classical trajectory calculations under two different initial conditions: at a fixed room temperature, for direct comparison with the experimental evidence and at different collision energies, to analyze possible mechanisms of reaction. In the first case, the available energy was mostly deposited as water vibrational energy, with the vibrational population inverted in the stretching modes and not inverted in the bending modes, reproducing the experimental evidence. In the second case, the analysis of different dynamics magnitudes (excitation functions, product energy partitioning, and product scattering distributions), allows us to suggest different mechanisms for both (R1) and (R2) paths: a direct mechanism for the (R2) path vs an indirect one, related with "nearly trapped" trajectories in the intermediate complexes, for the (R1) path. [ABSTRACT FROM AUTHOR]

To facilitate the data analysis of current and future high-resolution space telescope missions, we adopt 'Best Theory + Reliable High-resolution Experiment' (BTRHE) strategy to develop highly accurate infrared line lists for nitrous oxide (N2O). The 'Ames-1' potential energy surface (PES) is a CCSD(T)/aug-cc-pVQZ PES refined using selected HITRAN experimental data, with σrms = 0.02–0.03 cm-1 for five isotopologues. The 'Ames-1' dipole moment surface (DMS) is fitted from CCSD(T)/aug-cc-pV(T,Q,5)Z dipoles extrapolated to one electron basis set limit. Using the Ames-1 PES and DMS, Ames-296K line lists are computed in the full range of 0–15,000 cm-1 for 12 N2O isotopologues, with S296K ≥ 10−31 cm-1/molecule.cm-2. The reliability and consistency of Ames-296K intensity predictions (SAmes) are demonstrated through comparisons with HITRAN (SHITRAN), NOSL-296 (SNOSL), recent observed intensities (Sobs) and Effective Dipole Model (EDM) intensities (SEDM). Agreements and discrepancies are discussed, along with preliminary uncertainty estimate for SAmes. The SAmes provides a good constraint to prevent substantial errors in intensity predictions (e.g. for weak bands and minor isotopologues) and can be further improved. Ames-296K and NOSL-296 may complement each other to provide improved input for future database updates, combining the strengths of EH/EDM and BTRHE approaches. Data available at and https://doi.org/10.48667/9kmk-0334. [ABSTRACT FROM AUTHOR]

The fine structure excitation of the interstellar CCS radical induced by collisions with He is investigated. The first potential energy surface (PES) for the CCS–He van der Waals complex is presented. It was obtained from a highly correlated spin unrestricted coupled cluster approach with single double and perturbative triple excitations. The PES presents two shallow minima of 31.85 and 37.12 cm−1 for the linear (He facing S) and the nearly T-shaped geometries, respectively. The dissociation energy of the complex was calculated and found to be D0 = 14.183 cm−1. Inelastic scattering calculations were performed using the close-coupling approach. Cross-sections for transitions between the 61 first fine structure levels of CCS were obtained for energy up to 600 cm−1 and rate coefficients for the 5–50 K temperature range were derived. This set of collisional data can be used to model CCS emission spectra in dark molecular interstellar clouds and circumstellar envelopes and enable an accurate determination of CCS abundance in these astrophysical media. [ABSTRACT FROM AUTHOR]

We present ab initio calculations of the collisional broadening of the R(0) pure rotational line in CO (at 115 GHz) perturbed by O2. Our calculations are done in a fully quantum way by solving close-coupling quantum-scattering equations without any approximations. We also report a new, highly accurate CO–O2 potential energy surface on which we did the quantum-scattering calculations. The calculated collisional broadening agrees with the available experimental data in a wide temperature range. The calculated collisional shift is negligible compared to the broadening, which is also consistent with the experimental data. We combine this result with our previous calculations for the same line in CO perturbed by N2 [Jóźwiak et al., J. Chem. Phys. 154, 054314 (2021)]; the obtained air-perturbed broadening of the R(0) pure rotational line in CO and its temperature dependence perfectly agree with the HITRAN database. This result constitutes an important step toward developing a methodology for providing accurate ab initio reference data on spectroscopic collisional line-shape parameters for molecular systems relevant to the Earth's atmosphere and for populating spectroscopic line-by-line databases. [ABSTRACT FROM AUTHOR]

The status of machine learning in chemistry and related disciplines is becoming increasingly prominent, but its application in undergraduate chemistry experiments has yet to be seen. Existing computational chemistry experiments mostly focus on the application of quantum chemistry methods to study molecular properties and chemical reactions. In order to popularize machine learning as a powerful tool, this article will introduce a comprehensive computational chemistry experiment suitable for senior undergraduate students or graduate students. Through this experiment, students can gain a preliminary understanding of the basic principles and operational procedures of machine learning, while also develop their ability to apply quantum chemistry and machine learning to solve chemical problems. [ABSTRACT FROM AUTHOR]

The effect of para -hydrogen (j p = 0) collisions on the rotational de-excitation transitions of molecule NCCP is investigated in this study. The scattering information is obtained by spherically averaging a four-dimensional potential energy surface (4DPES) over various H2 molecule orientations. The calculations used the CCSD(T)-F12a method and aug-cc-pVTZ basis set to generate a 4DPES for the NCCP–H2 van der Waals system. Within the NCCP– para -H2 4DPES, a minimum energy point of 191.82 cm−1 is attained at a distance of 3.6 Å from the centre of mass of H2 and NCCP. To compute cross-sectional data for NCCP interacting with para -H2 (j p = 0), close coupling calculations are employed, encompassing total energies up to 600 cm−1. The resulting rate coefficients [ |$k_{j{\rightarrow }j^{\prime }}(T$|⁠)] are calculated across a temperature range spanning from 5 K to 200 K. In accordance with propensity, even Δ j = −2 transitions are highly preferred. Comparatively, the derived |$k_{j{\rightarrow }j^{\prime }}(T$|⁠) for NCCP–H2 are determined to be 1.5–4.5 times of NCCP–He. This observation implies that relying on a scaling factor of 1.38 to extrapolate rate coefficients for H2 collisions from those of NCCP–He is not a reliable approach. [ABSTRACT FROM AUTHOR]

The bimolecular reaction between OH− and CH3F is not just a prototypical SN2 process, but it has three other product channels. Here, we develop an accurate full-dimensional potential energy surface (PES) based on 191 193 points calculated at the level CCSD(T)-F12a/aug-cc-pVTZ. A detailed dynamics and mechanism analysis was carried out on this potential energy surface using the quasi-classical trajectory approach. It is verified that the trajectories do not follow the minimum energy path (MEP), but directly dissociate to F− and CH3OH. In addition, a new transition state for proton exchange and a new product complex CH2F−⋯H2O for proton abstraction were discovered. The trajectories avoid the transition state or this complex, instead dissociate to H2O and CH2F− directly through the ridge regions of the minimum energy path before the transition state. These non-MEP dynamics become more pronounced at high collision energies. Detailed dynamic simulations provide new insights into the atomic-level mechanisms of the title reaction, thanks to the new chemically accurate PES, with the aid of machine learning. [ABSTRACT FROM AUTHOR]

The effect of surface mode vibrations on the reactive scattering of D2, initialized in the ground rovibrational state (v = 0, j = 0), from a Cu(111) surface is investigated for different surface temperature situations. We adopt a time and temperature dependent effective Hamiltonian [Dutta et al., J. Chem. Phys. 154, 104103 (2021)] constructed by combining the linearly coupled many oscillator model [Sahoo et al., J. Chem. Phys. 136, 084306 (2012)] and the static corrugation model [M. Wijzenbroek and M. F. Somers, J. Chem. Phys. 137, 054703 (2012)] potential within the mean-field approach. Such an effective Hamiltonian is employed for six-dimensional quantum dynamical calculations to obtain temperature dependent reaction and state-to-state scattering probability profiles as a function of incidence energy of colliding D2 molecules. As reported in the experimental studies, the movements of surface atoms modify the dissociative scattering dynamics at higher surface temperature by exhibiting vibrational quantum and surface atoms' recoil effects in the low and high collision energy domains, respectively. Finally, we compare our present theoretical results with the experimental and other theoretical outcomes, as well as discuss the novelty of our findings. [ABSTRACT FROM AUTHOR]

We present a new ab initio five-dimensional potential energy surface for the CO–CO2 complex containing the Q3 normal mode for the ν3 asymmetric stretching vibration of the CO2 unit. The potential was calculated by the supermolecular approach at the explicitly correlated coupled cluster [CCSD (T)-F12a] level with aug-cc-pVTZ basis set plus midpoint bond functions. Two vibrationally averaged four-dimensional potentials for CO–CO2 with CO2 at the ground and ν3 excited states were generated by the integration of the five-dimensional potential over the Q3 intramolecular coordinate. Each potential displays a T-shaped global minimum with the C end in the CO unit pointing toward the C atom in the CO2 unit and a T-shaped local minimum but with the CO monomer rotated by 180°. The rovibrational bound states and energy levels for the CO–CO2 dimer were obtained employing the radial discrete variable representation/angular finite basis representation method in conjunction with the Lanczos algorithm. The vibrational ground and some lower excited states for CO–CO2 are localized around the global minimum because of the higher potential barriers. The band origin is blueshifted by 0.2089 cm−1 for CO–CO2 in the CO2ν3 range, which is consistent with the experimental result of 0.211 cm−1. The geared bending vibrational frequencies for CO–CO2 are 24.7101 and 24.5549 cm−1 at the ground and ν3 excited states of CO2, respectively. The predicted rovibrational frequencies, as well as spectral constants, coincide with the available observations, and these parameters show that the CO–CO2 complex is a nearly prolate asymmetric rotor. [ABSTRACT FROM AUTHOR]

We develop a full-dimensional analytical potential energy surface (PES) for the OH + C2H6 reaction using the Robosurfer program system, which automatically (1) selects geometries from quasi-classical trajectories, (2) performs ab initio computations using a coupled-cluster singles, doubles, and perturbative triples-F12/triple-zeta-quality composite method, (3) fits the energies utilizing the permutationally invariant monomial symmetrization approach, and (4) iteratively improves the PES via steps (1)–(3). Quasi-classical trajectory simulations on the new PES reveal that hydrogen abstraction leading to H2O + C2H5 dominates in the collision energy range of 10–50 kcal/mol. The abstraction cross sections increase and the dominant mechanism shifts from rebound (small impact parameters and backward scattering) to stripping (larger impact parameters and forward scattering) with increasing collision energy as opacity functions and scattering angle distributions indicate. The abstraction reaction clearly favors side-on OH attack over O-side and the least-preferred H-side approach, whereas C2H6 behaves like a spherical object with only slight C–C-perpendicular side-on preference. The collision energy efficiently flows into the relative translation of the products, whereas product internal energy distributions show only little collision energy dependence. H2O/C2H5 vibrational distributions slightly/significantly violate zero-point energy and are nearly independent of collision energy, whereas the rotational distributions clearly blue-shift as the collision energy increases. [ABSTRACT FROM AUTHOR]

A detailed exploration of the potential energy surface of quinoline cation (C9H7N·+) is carried out to extend the present understanding of its fragmentation mechanisms. Density functional theory calculations have been performed to explore new fragmentation schemes, giving special attention to previously unexplored pathways, such as isomerization and elimination of HNC. The isomerization mechanisms producing five- to seven-membered ring intermediates are described and are found to be a dominant channel both energetically and kinetically. Energetically competing pathways are established for the astrochemically important HNC-loss channel, which has hitherto never been considered in the context of the loss of a 27 amu fragment from the parent ions. Elimination of acetylene was also studied in great detail. Overall, the computational results are found to complement the experimental observations from the concurrently conducted PEPICO investigation. These could potentially open the doors for rich and interesting vacuum ultraviolet radiation-driven chemistry on planetary atmospheres, meteorites, and comets. [ABSTRACT FROM AUTHOR]

Due to the lack of specific collisional data, the abundance of NS+ in cold dense interstellar clouds was determined using collisional rate coefficients of CS as a substitute. To better understand the chemistry of sulfur in the interstellar medium, further abundance modeling using the actual NS+ collisional rate coefficients is needed. For this purpose, we have computed the first full 4D potential energy surface of the NS+–H2 van der Waals complex using the explicitly correlated coupled cluster approach with single, double, and non-iterative triple excitation in conjunction with the augmented-correlation consistent-polarized valence triple zeta basis set. The potential energy surface exhibits a global minimum of 848.24 cm−1 for a planar configuration of the complex. The long-range interaction energy, described using multipolar moments, is sensitive to the orientation of H2 up to radial distances of ∼50 a0. From this new interaction potential, we derived excitation cross sections, induced by collision with ortho- and para-H2, for the 15 low-lying rotational levels of NS+ using the quantum mechanical close-coupling approach. By thermally averaging these data, we determined downward rate coefficients for temperatures up to 50 K. By comparing them with the previous NS+–H2 data, we demonstrated that reduced dimensional approaches are not suited for this system. In addition, we found that the CS collisional data underestimate our results by up to an order of magnitude. The differences clearly indicate that the abundance of NS+, in cold dense clouds retrieved from observational spectra, must be reassessed using these new collisional rate coefficients. [ABSTRACT FROM AUTHOR]

Carbon disulfide dimer (CS2)2 is a model system that has been widely studied in the field of computational chemistry because of its relevance to a variety of chemical and biological processes. The (CS2)2 dimer is a relatively simple molecular system composed of two carbon disulfide (CS2) molecules interacting with each other through intermolecular forces. Despite its apparent simplicity, the (CS2)2 dimer exhibits a rich array of structural and dynamical properties that are of great interest to researchers. In this research, we present an ab initio study of the intermolecular interactions in the carbon disulfide dimer (CS2)2 using an improved Lennard–Jones (ILJ) potential with CCSD(T)/QZVPP calculations. The potential energy surface of (CS2)2 is calculated using high‐level quantum mechanical calculations based on the CCSD(T)/def2‐qzvpp method, which accurately accounts for electron correlation effects. The resulting potential energy surface is then fitted to an ILJ potential energy function, which includes both long‐range dipole–dipole interactions and short‐range repulsive interactions. The calculated potential energy surface reveals a rich variety of structural and dynamical properties of (CS2)2, including multiple minima and saddle points, which are sensitive to the relative orientation of the two CS2 molecules. It is essential to use extended basis sets to accurately incorporate the significant quadrupole moment of CS2, which we have calculated to be 2.44 a.u. The results of this study demonstrate the importance of using high‐level ab initio methods for the accurate calculation of potential energy surfaces in complex molecular systems such as (CS2)2. The use of an ILJ potential, which takes into account both dipole–dipole interactions and short‐range repulsive interactions, provides a more accurate and efficient approach for modeling intermolecular interactions in (CS2)2 and other similar systems. The results of this study will be useful for understanding the behavior of carbon disulfide dimers in different environments and for the development of new materials and chemical processes. [ABSTRACT FROM AUTHOR]

Based on an analytical full‐dimensional potential energy surface (PES), named PES‐2022, fitted to high‐level ab initio calculations previously developed by our group and specifically developed to describe this polyatomic reactive process, an exhaustive kinetics analysis was performed in the temperature range 50–2000 K, that is, interstellar, atmospheric and combustion conditions. Using the competitive canonical unified theory with multidimensional tunneling corrections of small curvature, CCUS/SCT, and low‐ and high‐pressure limit (LPL and HPL) models, in this wide temperature range we found that the overall rate constants increase with temperature at T > 300 K and T < 200 K, showing a V‐shaped temperature dependence, reproducing the experimental evidence when the HPL model was used. The increase of the rate constant with temperature at low temperatures was due to the strong contribution of the tunneling factor. The title reaction evolves by two paths, H2O + CH2OH (R1) and H2O + CH3O (R2), and the branching ratio analysis showed that the R2 path was dominant at T < 200 K while the R1 path dominated at T > 300 K, with a turnover temperature of ∼260 K, in agreement with previous theoretical estimations. Three kinetics isotope effects (KIEs), 13CH3OH, CH318OH, and CD3OH, were theoretically studied, reproducing the experimental evidence. The kinetics analysis in the present paper together with the dynamics study previously reported showed the capacity of the PES‐2022 to understand this important chemical process. [ABSTRACT FROM AUTHOR]

Intrinsic properties of metal-organic frameworks (MOFs), such as their ultra porosity and high surface area, deem them promising solutions for problems involving gas adsorption. Nevertheless, due to their combinatorial nature, a huge number of structures is feasible which renders cumbersome the selection of the best candidates with traditional techniques. Recently, machine learning approaches have emerged as efficient tools to deal with this challenge, by allowing researchers to rapidly screen large databases of MOFs via predictive models. The performance of the latter is tightly tied to the mathematical representation of a material, thus necessitating the use of informative descriptors. In this work, a generalized framework to predict gaseous adsorption properties is presented, using as one and only descriptor the capstone of chemical information: the potential energy surface (PES). In order to be machine understandable, the PES is voxelized and subsequently a 3D convolutional neural network (CNN) is exploited to process this 3D energy image. As a proof of concept, the proposed pipeline is applied on predicting CO 2 uptake in MOFs. The resulting model outperforms a conventional model built with geometric descriptors and requires two orders of magnitude less training data to reach a given level of performance. Moreover, the transferability of the approach to different host-guest systems is demonstrated, examining CH 4 uptake in COFs. The generic character of the proposed methodology, inherited from the PES, renders it applicable to fields other than reticular chemistry. [ABSTRACT FROM AUTHOR]

We report the derivation of rate coefficients for the rotational (de-)excitation of PO+ induced by collisions with H2. The calculations were performed on a 4D potential energy surface, obtained on top of highly accurate ab initio energy points. Preliminary tests pointed out the low influence of the coupling between j  = 0 and the higher rotational levels of H2 on the cross-sections values, thus allowing to neglect the rotational structure of H2. On this basis, state-to-state collisional rate coefficients were derived for temperatures ranging from 5 to 200 K. Radiative transfer calculations have been used to model the recent observation of PO+ in the G+0.693–0.027 molecular cloud, in order to evaluate the possible impact of non-LTE models on the determination of its physical conditions. The derived column density was found to be approximately ∼ 3.7 × 1011 cm−2, which is 60% (a factor of ∼ 1.7) smaller than the previously LTE-derived value. Extensive simulations show that PO+ low- j rotational lines exhibit maser behaviour at densities between 104 and 106 cm−3, thus highlighting the importance of a proper treatment of the molecular collisions to accurately model PO+ emissions in the interstellar medium. [ABSTRACT FROM AUTHOR]

In this work, we report rovibrational energy levels for four isotopologues of methylene (CH2, CHD, CD2, and 13CH2) in their ground triplet electronic state (X~3B1) from variational calculation up to ~10,000 cm−1 and using a new accurate ab initio potential energy surface (PES). Triplet methylene exhibits a large‐amplitude bending vibration and can reach a quasilinear configuration due to its low barrier (~2000 cm−1). To construct the ab initio PES, the Dunning's augmented correlation‐consistent core‐valence orbital basis sets were employed up to the sextuple‐ζ quality [aug‐cc‐pCVXZ, X = T, Q, 5, and 6] combined with the single‐ and double‐excitation unrestricted coupled cluster approach with a perturbative treatment of triple excitations [RHF‐UCCSD(T)]. We have shown that the accuracy of the ab initio energies is further improved by including the corrections due to the scalar relativistic effects, DBOC and high‐order electronic correlations. For the first time, all the available experimental rovibrational transitions were reproduced with errors less than 0.12 cm−1, without any empirical corrections. Unlike more "traditional" nonlinear triatomic molecules, we have shown that even the energies of the ground vibrational state (000) with rather small rotational quantum numbers are strongly affected by the very pronounced rovibrational resonance interactions. Accordingly, the polyad structure of the vibrational levels of CH2 and CD2 was analyzed and discussed. The comparison between the energy levels obtained from the effective Watson A‐reduced Hamiltonian, from the generating‐function approach and from a variational calculation was given. [ABSTRACT FROM AUTHOR]

A highly accurate, (HF)2 potential energy surface (PES) is constructed based on ab initio calculations performed at the coupled-cluster single double triple level of theory with an aug-cc-pVQZ-F12 basis set at about 152 000 points. A higher correlation correction is computed at coupled-cluster single double triple quadruple level for 2000 points and is considered alongside other more minor corrections due to relativity, core-valence correlation, and Born–Oppenheimer failure. The analytical surface constructed uses 500 constants to reproduce the ab initio points with a standard deviation of 0.3 cm−1. Vibration–rotation–inversion energy levels of the HF dimer are computed for this PES by variational solution of the nuclear-motion Schrödinger equation using the program WAVR4. Calculations over an extended range of rotationally excited states show very good agreement with the experimental data. In particular, the known empirical rotational constants B for the ground vibrational states are predicted to better than about 2 MHz. B constants for excited vibrational states are reproduced several times more accurately than by previous calculations. This level of accuracy is shown to extend to higher excited inter-molecular vibrational states v and higher excited rotational quantum numbers (J, Ka). [ABSTRACT FROM AUTHOR]

The CO2 molecule is of great interest for astrophysical studies since it can be found in a large variety of astrophysical media where it interacts with the dominant neutral species, such as He, H2, or H2O. The CO2–He collisional system was intensively studied over the last two decades. However, collisional data appear to be very sensitive to the potential energy surface (PES) quality. Thus, we provide, in this study, a new PES of the CO2–He van der Waals complex calculated with the coupled-cluster method and a complete basis set extrapolation in order to provide rotational rate coefficients that are as accurate as possible. The PES accuracy was tested through the calculations of bound state transition frequencies and pressure broadening coefficients that were compared to experimental data. An excellent agreement was globally found. Then, revised collisional data were provided for the 10–300 K temperature range. Rate coefficients were compared to previously computed ones and are found to be up to 50% greater than previously provided ones. These differences can induce non-negligible consequences for the modeling of CO2 abundance in astrophysical media. [ABSTRACT FROM AUTHOR]

With gates of a quantum computer designed to encode multi-dimensional vectors, projections of quantum computer states onto specific qubit states can produce kernels of reproducing kernel Hilbert spaces. We show that quantum kernels obtained with a fixed ansatz implementable on current quantum computers can be used for accurate regression models of global potential energy surfaces (PESs) for polyatomic molecules. To obtain accurate regression models, we apply Bayesian optimization to maximize marginal likelihood by varying the parameters of the quantum gates. This yields Gaussian process models with quantum kernels. We illustrate the effect of qubit entanglement in the quantum kernels and explore the generalization performance of quantum Gaussian processes by extrapolating global six-dimensional PESs in the energy domain. [ABSTRACT FROM AUTHOR]

The reaction HO + SO → H + SO2 (Rt) and its reverse (R‐t) play an important role in environment and the combustion of sulfur‐containing fuels. However, their kinetics is of high uncertainty as its reaction profile is complicated with multiple deep complexes and channels. In this work, the kinetics and mechanisms of Rt and R‐t are studied comprehensively based on a newly developed full‐dimensional accurate potential energy surface (PES) with the aid of machine learning. This highly accurate PES is interfaced with the software Gaussian. Then reliable information, including the energy, structures, and vibrational frequencies of the stationary points, as well as the minimum energy path and variational analysis can be efficiently determined. The variational transition state theory (VTST) and Rice−Ramsperger−Kassel−Marcus (RRKM) theory are employed to obtain the rate coefficients of each elementary reaction. The temperature‐ and pressure‐dependent rate coefficients of Rt are derived by the RRKM‐based master equation with hindered rotor and free rotor model considered. In addition, the effect of isotope substitution for the hydrogen is investigated on the reaction kinetics. Meanwhile, the quasi‐classical trajectory (QCT) calculation is performed on the PES‐2020 to obtain the temperature‐dependent reaction kinetics. [ABSTRACT FROM AUTHOR]

We present a 3D isotropic ab initio three-body (para-H2)3 interaction potential energy surface (PES). The electronic structure calculations are carried out at the correlated coupled-cluster theory level, with single, double, and perturbative triple excitations. The calculations use an augmented correlation-consistent triple zeta basis set and a supplementary midbond function. We construct the PES using the reproducing-kernel Hilbert space toolkit [O. T. Unke and M. Meuwly, J. Chem. Inf. Model. 57, 1923 (2017)] with phenomenological and empirical adjustments to account for short-range and long-range behaviors. The (para-H2)3 interaction energies deviate drastically from the Axilrod–Teller–Muto (ATM) potential at short intermolecular separations. We find that the configuration of three para-H2 molecules at the corners of an equilateral triangle is responsible for the majority of the (para-H2)3 interaction energy contribution in a hexagonal-close-packed lattice. In cases where two para-H2 molecules are close to one another while the third is far away, the (para-H2)3 interaction PES takes the form of a modified version of the ATM potential. We expect the combination of this PES together with a first-principles para-H2–para-H2 adiabatic hindered rotor potential to outperform a widely used effective pair potential for condensed many-body systems of para-H2. [ABSTRACT FROM AUTHOR]

The H-atom dissociation of formaldehyde on the lowest triplet state (T1) is studied by quasi-classical molecular dynamic simulations on the high-dimensional machine-learning potential energy surface (PES) model. An atomic-energy based deep-learning neural network (NN) is used to represent the PES function, and the weighted atom-centered symmetry functions are employed as inputs of the NN model to satisfy the translational, rotational, and permutational symmetries, and to capture the geometry features of each atom and its individual chemical environment. Several standard technical tricks are used in the construction of NN-PES, which includes the application of clustering algorithm in the formation of the training dataset, the examination of the reliability of the NN-PES model by different fitted NN models, and the detection of the out-of-confidence region by the confidence interval of the training dataset. The accuracy of the full-dimensional NN-PES model is examined by two benchmark calculations with respect to ab initio data. Both the NN and electronic-structure calculations give a similar H-atom dissociation reaction pathway on the T1 state in the intrinsic reaction coordinate analysis. The small-scaled trial dynamics simulations based on NN-PES and ab initio PES give highly consistent results. After confirming the accuracy of the NN-PES, a large number of trajectories are calculated in the quasi-classical dynamics, which allows us to get a better understanding of the T1-driven H-atom dissociation dynamics efficiently. Particularly, the dynamics simulations from different initial conditions can be easily simulated with a rather low computational cost. The influence of the mode-specific vibrational excitations on the H-atom dissociation dynamics driven by the T1 state is explored. The results show that the vibrational excitations on symmetric C–H stretching, asymmetric C–H stretching, and C=O stretching motions always enhance the H-atom dissociation probability obviously. [ABSTRACT FROM AUTHOR]

The automated development of a new ab initio full-dimensional potential energy surface (PES) is reported for the CH4·F− complex using the Robosurfer program package. The new potential provides a near-spectroscopic quality description over a broad configuration range including the methane-ion dissociation, as well as isolated methane vibrations. In particular, it improves upon the earlier [Czakó, Braams, Bowman (2008)] PES over intermediate methane-fluoride distances. Full-dimensional (12D) variational vibrational computations using the new PES and the GENIUSH-Smolyak algorithm show that tunnelling splittings larger than 0.1 cm−1 appear below the top of the interconversion barrier of the four equivalent minima of the complex. [ABSTRACT FROM AUTHOR]

We present a full-dimensional ab initio intermolecular potential energy surface (IPES) for the OC–HF van der Waals complex. 3167 ab initio points were computed at the frozen-core (FC) explicitly correlated coupled cluster [FC-CCSD(T)-F12b] level, with the augmented correlation-consistent polarized valence quadruple-zeta basis set plus bond functions. Basis set superposition error correction was also considered by the full counterpoise procedure. Gaussian process regression (GPR) was used to map out the potential energy surface, while a multipole expansion method was employed to smooth the ab initio noise of intermolecular potential in the long range. The global minimum of −1248.364 cm−1 was located at the linear configuration with the C atom pointing toward the H atom of the HF molecule. In addition, a local minimum of −602.026 cm−1 was found at another linear configuration with the O atom pointing toward the H atom of the HF molecule. The eigenstates were calculated on the vibrational averaged four-dimensional IPESs with the mixed radial discrete variable representation/angular finite basis representation method and Lanczos propagation algorithm. The dissociation energy D0 was calculated to be 701.827 cm−1, well reproducing the experimental value of 732 ± 2 cm−1. The dipole moment surfaces were also fitted by GPR from 3132 ab initio points calculated using the coupled cluster method [CCSD(T)] with AVTZ basis set plus bond functions. The frequencies and relative line intensities of rovibrational transitions in the HF (DF) and CO stretching bands were further calculated and compared well with the experimental results. These results indicate the high fidelity of the new IPES. [ABSTRACT FROM AUTHOR]

Developing global, high-dimensional potential energy surfaces (PESs) is a formidable task. Beside the challenges of PES fitting and fitting set generation, one also has to choose an electronic structure method capable of delivering accurate potential energy values for all geometries in the fitting set, even in regions far from equilibrium. Such regions are often plagued by Hartree–Fock (HF) convergence issues, and even if convergence is achieved, self-consistent field (SCF) procedures that are used to obtain HF solutions offer no guarantee that the solution found is the lowest-energy solution. We present a study of the reactant regions of CH3OH + OH·, C2H6 + F·, and CH3NH2 + Cl·, where the SCF procedure often converges to a higher-energy state or fails to converge, resulting in erratic post-HF energies and regions where no energy is obtained, both of which are major obstacles for PES development. We introduce a pragmatic method for automatically finding better HF solutions (dubbed ManyHF) and present evidence that it may extend the applicability of single-reference methods to some systems previously thought to require multireference methods. [ABSTRACT FROM AUTHOR]

Collisional data for the excitation of NH by H2 are key to accurately derive the NH abundance in astrophysical media. We present a new four-dimensional potential energy surface (PES) for the NH–H2 van der Waals complex. The ab initio calculations of the PES were carried out using the explicitly correlated partially spin-restricted coupled cluster method with single, double, and perturbative triple excitations [RCCSD(T)-F12a] with the augmented correlation-consistent polarized valence triple zeta basis set. The PES was represented by an angular expansion in terms of coupled spherical harmonics. The global minimum corresponds to the linear structure with a well depth De = 149.10 cm−1. The calculated dissociation energy D0 is found to be 30.55 and 22.11 cm−1 for ortho-H2 and para-H2 complexes, respectively. These results are in agreement with the experimental values. Then, we perform quantum close-coupling calculations of the fine structure resolved excitation cross sections of NH induced by collisions with ortho-H2 and para-H2 for collisional energies up to 500 cm−1. We find strong differences between collisions induced by ortho-H2 and para-H2. Propensity rules are discussed. The cross sections are larger for fine structure conserving transitions than for fine structure changing ones, as predicted by theory. These new results should help in interpreting NH interstellar spectra and better constrain the abundance of NH in interstellar molecular clouds. [ABSTRACT FROM AUTHOR]