A new potential surface and quasiclassical trajectory study of H+H2O→OH+H2

High vibrational excitation of the reagent transforms the late-barrier H + HOD reaction into an early-barrier reaction

Polanyi's rules predict that a late-barrier reaction yields vibrationally cold products; however, experimental studies showed that the H2 product from the late-barrier H + H2O(|04⟩-) and H + HOD(vOH = 4) reactions is vibrationally hot. Here, we report a potential-averaged five-dimensional state-to-state quantum dynamics study for the H + HOD(vOH = 0-4) → H2 + OD reactions on a highly accurate potential energy surface with the total angular momentum J = 0. It is found that with the HOD vibration excitation increasing from vOH = 1 to 4, the product H2 becomes increasingly vibrationally excited and manifests a typical characteristic of an early barrier reaction for vOH = 3 to 4. Analysis of the scattering wave functions revealed that vibrational excitation in the breaking OH bond moves the location of dynamical saddle point from product side to reactant side, transforming the reaction into an early barrier reaction. Interestingly, pronounced oscillatory structures in the total and product vibrational-state-resolved reaction probabilities were observed for the H + HOD(vOH = 3, 4) reactions, in particular at low collision energies, which originate from the Feshbach resonance states trapped in the bending/torsion excited vibrational adiabatic potential wells in the entrance region due to van der Waals interactions.

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. The parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

VUV photochemistry and nuclear spin conversion of water and water-orthohydrogen complexes in parahydrogen crystals at 4 K

Samples of H2O, HDO, and D2O were isolated in solid parahydrogen (pH2) matrices and irradiated by vacuum ultraviolet (VUV) radiation at 147 nm. Fourier-Transform Infrared (FTIR) spectra showed a clear depletion of D2O and an enrichment of both HDO and H2O by 147 nm irradiation. These irradiation-dependent changes are attributed to the production of OH and/or OD radicals through photodissociations of H2O, HDO, and D2O. The radicals subsequently react with the hydrogen matrix, leading to the observed enrichment of H2O. No trace of isolated OH or OD was detected in the FTIR spectra, indicating that the OH/OD radicals react with the surrounding matrix hydrogen molecules via quantum tunneling within our experimental timescale. The observed temporal changes in concentrations, especially the increase of HDO concentration during VUV irradiation, can be interpreted by a model with a rapid conversion from orthohydrogen (oH2) to pH2 in water-oH2 complexes upon VUV photodissociation, indicating either the acceleration of the nuclear spin conversion (NSC) of H2 due to the magnetic moment of the intermediate OH/OD radical, or the preferential reaction of the OH/OD radical with a nearby oH2 molecule over other pH2 molecules. We have also identified and quantified an anomalously slow NSC of H2O and D2O complexed with oH2 in solid pH2.

State-to-state mode specificity in H + DOH(νOH = 1) → HD + OH(ν2 = 0) reaction: vibrational non-adiabaticity or local-mode excitation?

State-of-the-art full dimensional state-to-state quantum dynamics reveal a startling observation in which the DOH(ν_OH = 1) molecule reacts with a H atom to produce a vibrationless OH product. This interesting observation is attributed to a small OD excited local-mode component in the reactant wavefunction.

Recent advances in quantum scattering calculations on polyatomic bimolecular reactions

This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.

A ring polymer molecular dynamics study of the OH + H2(D2) reaction

Using ring polymer molecular dynamics we have calculated the rate coefficients for the OH + H₂ reaction.

N-dimensional switch function for energy conservation in multiprocess reaction dynamics

The MReaDy program was designed for studying Multiprocess Reactive Dynamic systems, that is, complex chemical systems involving different and concurrent reactions. It builds a global potential energy surface integrating a variety of potential energy surfaces, each one of them representing an elementary reaction expected to play a role in the chemical process. For each elementary reaction, energy continuity problems may happen in the transition between potential energy surfaces due to differences in the functional form for each of the fragments, especially if built by different authors. A N-dimensional switch function is introduced in MReaDy in order to overcome such a problem. As an example, results of a collision trajectory calculation for H2  + OH → H3 O are presented, showing smooth transition in the potential energy, leading to conservation in the total energy. Calculations for a hydrogen combustion system from 1000 K up to 4000 K shows a variation of 0.012% when compared to the total energy of the system. © 2016 Wiley Periodicals, Inc.

A local mode picture for H atom reaction with vibrationally excited H2O: a full dimensional state-to-state quantum dynamics investigation

A local mode picture was postulated some time ago to explain the negligible OD product excitation from H reacting with D₂O in fundamental asymmetrical excitation: the D₂O asymmetric stretch state is best thought as a linear combination of local mode stretches with the H atom reacting preferentially with the excited OD bond, but has never been verified. Here we report the first full-dimensional state-to-state study for the title reaction with H₂O in the ground and the first symmetric and asymmetric stretching excited states. It is found that symmetric and asymmetric stretching excited states, which vibrate in very different ways in the normal mode picture, behave essentially identically in the reaction. More importantly, our calculations revealed that the reaction produces a small fraction of OH in the v = 1 state, with the population close to the relative reactivity between the ground and vibrationally excited states, and therefore it confirms the local mode picture for H₂O in symmetric or asymmetric stretch states when reacting with an H atom.

The READY program: Building a global potential energy surface and reactive dynamic simulations for the hydrogen combustion

READY (REActive DYnamics) is a program for studying reactive dynamic systems using a global potential energy surface (PES) built from previously existing PESs corresponding to each of the most important elementary reactions present in the system. We present an application to the combustion dynamics of a mixture of hydrogen and oxygen using accurate PESs for all the systems involving up to four oxygen and hydrogen atoms. Results at the temperature of 4000 K and pressure of 2 atm are presented and compared with model based on rate constants. Drawbacks and advantages of this approach are discussed and future directions of research are pointed out.

A classical trajectory study of the intramolecular dynamics, isomerization, and unimolecular dissociation of HO2

The classical dynamics and rates of isomerization and dissociation of HO2 have been studied using two potential energy surfaces (PESs) based on interpolative fittings of ab initio data: An interpolative moving least-squares (IMLS) surface [A. Li, D. Xie, R. Dawes, A. W. Jasper, J. Ma, and H. Guo, J. Chem. Phys. 133, 144306 (2010)] and the cubic-spline-fitted PES reported by Xu, Xie, Zhang, Lin, and Guo (XXZLG) [J. Chem. Phys. 127, 024304 (2007)]. Both PESs are based on similar, though not identical, internally contracted multi-reference configuration interaction with Davidson correction (icMRCI+Q) electronic structure calculations; the IMLS PES includes complete basis set (CBS) extrapolation. The coordinate range of the IMLS PES is limited to non-reactive processes. Surfaces-of-section show similar generally regular phase space structures for the IMLS and XXZLG PESs with increasing energy. The intramolecular vibrational energy redistribution (IVR) at energies above and below the threshold of isomerization is slow, especially for O-O stretch excitations, consistent with the regularity in the surfaces-of-section. The slow IVR rates lead to mode-specific effects that are prominent for isomerization (on both the IMLS and XXZLG) and modest for unimolecular dissociation to H + O2 (accessible only on the XXZLG PES). Even with statistical distributions of initial energy, slow IVR rates result in double exponential decay for isomerization, with the slower rate correlated with slow IVR rates for O-O vibrational excitation. The IVR and isomerization rates computed for the IMLS and XXZLG PESs are quantitatively, but not qualitatively, different from one another with the largest differences ascribed to the ~2 kcal/mol difference in the isomerization barrier heights. The IMLS and XXZLG results are compared with those obtained using the global, semi-empirical double-many-body expansion DMBE-IV PES [M. R. Pastrana, L. A. M. Quintales, J. Brandão, and A. J. C. Varandas, J. Chem. Phys. 94, 8073 (1990)], for which the surfaces-of-section display more irregular phase space structure, much faster IVR rates, and significantly less mode-specific effects in isomerization and unimolecular dissociation. The calculated IVR results for all three PESs are reasonably well represented by an analytic, coupled three-mode energy transfer model.

ROTATIONAL ALIGNMENT OF PRODUCT MOLECULES FROM THE REACTION Ca + HCl→CaCl + H

The product rotational polarization in the Ca + HCl→CaCl + H reaction at collision energy of 20 kcal/mol has been studied via the quasiclassical trajectory method on a new ab initio potential surface. The P(θ r ) distribution of angle between k and j′, and the dihedral angular distribution P(Φ r ) characterizing k - k′ - j′ correlation are discussed, the angle distribution P(θ r , Φ r ) of product rotational vectors in the form of polar plot in θ r and Φ r are shown. Furthermore, four PDDCSs (2π/σ)(dσ00/dωt), (2π/σ)(dσ20/dωt), (2π/σ)(dσ22+/dωt) and (2π/σ)(dσ21-/dωt) are also presented. The present calculations reveal that the product rotational alignment is very strong. Finally, the state distributions of the product CaCl are investigated. The results showed that the CaCl product was formed with high vibrational and rotational excitation.

Apparent failure of the Born-Oppenheimer static surface model for vibrational excitation of molecular hydrogen on copper

The accuracy of dynamical models for reactive scattering of molecular hydrogen, H(2), from metal surfaces is relevant to the validation of first principles electronic structure methods for molecules interacting with metal surfaces. The ability to validate such methods is important to progress in modeling heterogeneous catalysis. Here, we study vibrational excitation of H(2) on Cu(111) using the Born-Oppenheimer static surface model. The potential energy surface (PES) used was validated previously by calculations that reproduced experimental data on reaction and rotationally inelastic scattering in this system with chemical accuracy to within errors ≤ 1 kcal/mol ≈ 4.2 kJ/mol [Díaz C, et al. (2009) Science 326:832-834]. Using the same PES and model, our dynamics calculations underestimate the contribution of vibrational excitation to previously measured time-of-flight spectra of H(2) scattered from Cu(111) by a factor 3. Given the accuracy of the PES for the experiments for which the Born-Oppenheimer static surface model is expected to hold, we argue that modeling the effect of the surface degrees of freedom will be necessary to describe vibrational excitation with similar high accuracy.