skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction

Abstract

An accurate full-dimensional global potential energy surface (PES) is developed for the title reaction. While the long-range interactions in the reactant asymptote are represented by an analytical expression, the interaction region of the PES is fit to more than 81 000 of ab initio points at the UCCSD(T)-F12b/AVTZ level using the permutation invariant polynomial neural network approach. Fully symmetric with respect to permutation of all four hydrogen atoms, the PES provides a faithful representation of the ab initio points, with a root mean square error of 1.8 meV or 15 cm{sup −1}. The reaction path for this exoergic reaction features an attractive and barrierless entrance channel, a submerged saddle point, a shallow H{sub 4}O{sup +} well, and a barrierless exit channel. The rate coefficients for the title reaction and kinetic isotope effect have been determined on this PES using quasi-classical trajectories, and they are in good agreement with available experimental data. It is further shown that the H{sub 2}O{sup +} rotational enhancement of reactivity observed experimentally can be traced to the submerged saddle point. Using our recently proposed Sudden Vector Projection model, we demonstrate that a rotational degree of freedom of the H{sub 2}O{sup +} reactant is strongly coupled with themore » reaction coordinate at this saddle point, thus unraveling the origin of the pronounced mode specificity in this reaction.« less

Authors:
;  [1]
  1. Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131 (United States)
Publication Date:
OSTI Identifier:
22420101
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Physics; Journal Volume: 140; Journal Issue: 22; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; 36 MATERIALS SCIENCE; DEGREES OF FREEDOM; HYDROGEN; INTERACTION RANGE; INTERACTIONS; ISOTOPE EFFECTS; MEV RANGE; OXONIUM IONS; POLYNOMIALS; POTENTIAL ENERGY; SURFACES; WATER

Citation Formats

Li, Anyang, and Guo, Hua, E-mail: hguo@unm.edu. A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction. United States: N. p., 2014. Web. doi:10.1063/1.4881943.
Li, Anyang, & Guo, Hua, E-mail: hguo@unm.edu. A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction. United States. doi:10.1063/1.4881943.
Li, Anyang, and Guo, Hua, E-mail: hguo@unm.edu. Sat . "A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction". United States. doi:10.1063/1.4881943.
@article{osti_22420101,
title = {A nine-dimensional ab initio global potential energy surface for the H{sub 2}O{sup +} + H{sub 2} → H{sub 3}O{sup +} + H reaction},
author = {Li, Anyang and Guo, Hua, E-mail: hguo@unm.edu},
abstractNote = {An accurate full-dimensional global potential energy surface (PES) is developed for the title reaction. While the long-range interactions in the reactant asymptote are represented by an analytical expression, the interaction region of the PES is fit to more than 81 000 of ab initio points at the UCCSD(T)-F12b/AVTZ level using the permutation invariant polynomial neural network approach. Fully symmetric with respect to permutation of all four hydrogen atoms, the PES provides a faithful representation of the ab initio points, with a root mean square error of 1.8 meV or 15 cm{sup −1}. The reaction path for this exoergic reaction features an attractive and barrierless entrance channel, a submerged saddle point, a shallow H{sub 4}O{sup +} well, and a barrierless exit channel. The rate coefficients for the title reaction and kinetic isotope effect have been determined on this PES using quasi-classical trajectories, and they are in good agreement with available experimental data. It is further shown that the H{sub 2}O{sup +} rotational enhancement of reactivity observed experimentally can be traced to the submerged saddle point. Using our recently proposed Sudden Vector Projection model, we demonstrate that a rotational degree of freedom of the H{sub 2}O{sup +} reactant is strongly coupled with the reaction coordinate at this saddle point, thus unraveling the origin of the pronounced mode specificity in this reaction.},
doi = {10.1063/1.4881943},
journal = {Journal of Chemical Physics},
number = 22,
volume = 140,
place = {United States},
year = {Sat Jun 14 00:00:00 EDT 2014},
month = {Sat Jun 14 00:00:00 EDT 2014}
}
  • This work performs a time-dependent wavepacket study of the H{sub 2} + C{sub 2}H → H + C{sub 2}H{sub 2} reaction on a new ab initio potential energy surface (PES). The PES is constructed using neural network method based on 68 478 geometries with energies calculated at UCCSD(T)-F12a/aug-cc-pVTZ level and covers H{sub 2} + C{sub 2}H↔H + C{sub 2}H{sub 2}, H + C{sub 2}H{sub 2} → HCCH{sub 2}, and HCCH{sub 2} radial isomerization reaction regions. The reaction dynamics of H{sub 2} + C{sub 2}H → H + C{sub 2}H{sub 2} are investigated using full-dimensional quantum dynamics method. The initial-state selected reactionmore » probabilities are calculated for reactants in eight vibrational states. The calculated results showed that the H{sub 2} vibrational excitation predominantly enhances the reactivity while the excitation of bending mode of C{sub 2}H slightly inhibits the reaction. The excitations of two stretching modes of C{sub 2}H molecule have negligible effect on the reactivity. The integral cross section is calculated with J-shift approximation and the mode selectivity in this reaction is discussed. The rate constants over 200-2000 K are calculated and agree well with the experimental measured values.« less
  • Ab initio MCSCF+CI calculations of the minimum energy path in the doublet surface for the N(/sup 4/S)+O/sub 2/(/sup 3/..sigma../sup -//sub g/) ..-->..NO(/sup 2/Pi/sub u/)+O(/sup 3/P) reaction were carried out using a 4-31 G basis set. The minimum path was calculated to pass through a bent N--O--O /sup 2/A' transition complex with a barrier of 21.7 kcal/mole. A high barrier was found for the doublet C/sub 2v/ approach of N(/sup 4/S) to O/sub 2/(/sup 3/..sigma../sup -//sub g/). Two distinct regions of orbital transformations, given by the MCSCF calculations were found along the minimum energy path.
  • We report the first global potential energy surface (PES) for the X{sup  2}A{sup ′} ground electronic state of the Si({sup 3}P) + OH(X{sup 2}Π) → SiO(X{sup 1}Σ{sub g}{sup +}) + H({sup 2}S) reaction. The PES is based on a large number of ab initio energies obtained from multireference configuration interaction calculations plus Davidson correction (MRCI+Q) using basis sets of quadruple zeta quality. Corrections were applied to the ab initio energies in the reactant channel allowing a proper description of long-range interactions between Si({sup 3}P) and OH(X{sup 2}Π). An analytical representation of the global PES has been developed by means ofmore » the reproducing kernel Hilbert space method. The reaction is found barrierless. Two minima, corresponding to the SiOH and HSiO isomers, and six saddle points, among which the isomerization transition state, have been characterized on the PES. The vibrational spectra of the SiOH/HSiO radicals have been computed from second-order perturbation theory and quantum dynamics methods. The structural, energetic, and spectroscopic properties of the two isomers are in good agreement with experimental data and previous high quality calculations.« less
  • New, full-dimensional potential energy surfaces (PESs), obtained using precise least-squares fitting of high-level electronic energy databases, are reported for intrinsic H{sub 2}(H{sub 2}O) two-body and H{sub 2}(H{sub 2}O){sub 2} three-body potentials. The database for H{sub 2}(H{sub 2}O) consists of approximately 44 000 energies at the coupled cluster singles and doubles plus perturbative triples (CCSD(T))-F12a/haQZ (aug-cc-pVQZ for O and cc-pVQZ for H) level of theory, while the database for the three-body interaction consists of more than 36 000 energies at the CCSD(T)-F12a/haTZ (aug-cc-pVTZ for O, cc-pVTZ for H) level of theory. Two precise potentials are based on the invariant-polynomial technique and are comparedmore » to computationally faster ones obtained via “purified” symmetrization. All fits use reduced permutational symmetry appropriate for these non-covalent interactions. These intrinsic potentials are employed together with existing ones for H{sub 2}, H{sub 2}O, and (H{sub 2}O){sub 2}, to obtain full PESs for H{sub 2}(H{sub 2}O) and H{sub 2}(H{sub 2}O){sub 2}. Properties of these full PESs are presented, including a diffusion Monte Carlo calculation of the zero-point energy and wavefunction, and dissociation energy of the H{sub 2}(H{sub 2}O) dimer. These PESs together with an existing one for water clusters are used in a many-body representation of the PES of hydrogen clathrate hydrates, illustrated for H{sub 2}@(H{sub 2}O){sub 20}. An analysis of this hydrate is presented, including the electronic dissociation energy to remove H{sub 2} from the calculated equilibrium structure.« less
  • A new global potential energy surface for the ground electronic state (1{sup 2}A Prime ) of the Ar+H{sub 2}{sup +}{yields}ArH{sup +}+H reaction has been constructed by multi-reference configuration interaction method with Davidson correction and a basis set of aug-cc-pVQZ. Using 6080 ab initio single-point energies of all the regions for the dynamics, a many-body expansion function form has been used to fit these points. The quantum reactive scattering dynamics calculations taking into account the Coriolis coupling (CC) were carried out on the new potential energy surface over a range of collision energies (0.03-1.0 eV). The reaction probabilities and integral crossmore » sections for the title reaction were calculated. The significance of including the CC quantum scattering calculation has been revealed by the comparison between the CC and the centrifugal sudden approximation calculation. The calculated cross section is in agreement with the experimental result at collision energy 1.0 eV.« less