Kinetic and dynamic studies of the Cl(2Pu) + H2O(X̃1A1) → HCl(X̃1Σ+) + OH(X̃2Π) reaction on an ab initio based full-dimensional global potential energy surface of the ground electronic state of ClH2O

ABC+D: A time-independent coupled-channel quantum dynamics program for elastic and ro-vibrational inelastic scattering between atoms and triatomic molecules in full dimensionality

We discuss the details of a time-independent quantum mechanical method and its implementation for full-dimensional non-reactive scattering between a closed-shell triatomic molecule and a closed-shell atom. By solving the time-independent Schrödinger equation within the coupled-channel framework using a log-derivative method, the state-to-state scattering matrix (S-matrix) can be determined for inelastic scattering involving both the rotational and vibrational modes of the molecule. Various approximations are also implemented. The ABC+D code provides an important platform for understanding an array of physical phenomena involving collisions between atoms and molecules.

Full-dimensional quantum studies of vibrational energy transfer dynamics between H2O and Ar: theory assessing experiment

We report the first full-dimensional quantum mechanical calculations of the ro-vibrational inelastic scattering dynamics between water molecules and argon atoms on an accurate potential energy surface, using a recently developed time-independent quantum method based on the close-coupling approach. The state-to-state integral cross-sections and rate coefficients show strong observance of gap laws. The calculated thermal rate coefficients for the relaxation of the stretching fundamental states of H₂O are in good agreement with experimental values, while those for the bending overtone state are approximately five times smaller than the values extracted through a previous kinetic modeling of fluorescence decay data. Our state-specific quantum scattering results suggest the need to reassess the kinetic modeling of the experimental data. This work advanced our understanding of the quantum dynamics of vibrationally inelastic energy transfer processes involving polyatomic molecules.

OH + HCl Reaction on the Surface of Ice: An Ab Initio Molecular Dynamics Study

We have investigated the OH + HCl reaction on the surface of ice using Born-Oppenheimer molecular dynamics (BOMD) simulation. The present work revealed that the OH + HCl reaction becomes ∼1 order of magnitude faster on the ice surface compared to the bare reaction. The BOMD simulation also indicates that the Cl radical formed on the ice surface through the title reaction can form two hydrogen bonds at a time with the water molecules present on the ice surface; hence, the Cl radical cannot escape from the ice surface easily.

Potential energy profile for the Cl + (H2O)3 → HCl + (H2O)2OH reaction. A CCSD(T) study

Four different reaction pathways are initially located for the reaction of Cl atom plus water trimer Cl + (H₂O)₃ → HCl + (H₂O)₂OH using a standard DFT method. As found for the analogous fluorine reaction, the geometrical and energetic results for the four chlorine pathways are closely related. However, the energetics for the Cl reaction are very different from those for fluorine. In the present paper, we investigate the lowest-energy chlorine pathway using the "gold standard" CCSD(T) method in conjunction with correlation-consistent basis sets up to cc-pVQZ. Structurally, the stationary points for the water trimer reaction Cl + (H₂O)₃ may be compared to those for the water monomer reaction Cl + H₂O and water dimer reaction Cl + (H₂O)₂. Based on the CCSD(T) energies, the title reaction is endothermic by 19.3 kcal mol^-1, with a classical barrier height of 16.7 kcal mol^-1 between the reactants and the exit complex. There is no barrier for the reverse reaction. The Cl⋯(H₂O)₃ entrance complex lies 5.3 kcal mol^-1 below the separated reactants. The HCl⋯(H₂O)₂OH exit complex is bound by 8.6 kcal mol^-1 relative to the separated products. The Cl + (H₂O)₃ reaction is somewhat similar to the analogous Cl + (H₂O)₂ reaction, but qualitatively different from the Cl + H₂O reaction. It is reasonable to expect that the reactions between the chlorine atom and larger water clusters may be similar to the Cl + (H₂O)₃ reaction. The potential energy profile for the Cl + (H₂O)₃ reaction is radically different from that for the valence isoelectronic F + (H₂O)₃ system, which may be related to the different bond energies between HCl and HF.

Effect of microsolvation on the mode specificity of the OH˙(H2O) + HCl reaction

The present study investigates the mode specificity in the microsolvated OH˙(H₂O) + HCl reaction using on-the-fly direct dynamics simulation. To the best of our knowledge, this is the first study which aims to gain insights into the effect of microsolvation on the mode selectivity. Our investigation reveals that, similar to the gas phase OH˙ + HCl reaction, the microsolvated reaction is also predominantly affected by the vibrational excitation of the HCl mode, whereas the OH vibrational mode behaves as a spectator. Interestingly, in contrast to the behavior of the bare reaction, the integral cross section at the ground state of the microsolvated reaction decreases with an increase in translational energy. However, for the vibrational excited states, the reactivity of the microsolvated reaction is found to be higher than that of the bare reaction within the selected range of translational energies.

A Time-Independent Quantum Approach to Ro-vibrationally Inelastic Scattering between Atoms and Triatomic Molecules

A full-dimensional time-independent quantum mechanical theory for ro-vibrationally inelastic scattering of triatomic molecules with atoms is formulated. The Jacobi-Radau coordinate system used in the calculation allows not only a near perfect description of the vibrational problem but also the adaptation of the exchange symmetry for A₂B type triatoms. The S-matrix elements are obtained by solving the close-coupling equations with contracted basis using the log-derivative method. This method is applied to the inelastic scattering of the water molecule by a chlorine atom, which sheds light on the energy gap law in energy transfer in atom-triatom collisions.

Mode- and Bond-Selected Reaction of H with Local Mode Molecule HDS

Understanding mode- and bond-selected dynamics of elementary chemical reactions is of central importance in molecular reaction dynamics. The initial state-selected time-dependent wave packet method is employed to study the mode and bond selectivity, isotopic branching ratio, and temperature dependence of rate constants of the two-channel reaction of H with local mode molecule HDS. For the abstraction channel, fundamental excitation of the HS (DS) bond of the reactant HDS significantly enhances the H-abstraction (D-abstraction) reaction, whose efficacy is higher than the same amount of translational energy except at low energies just above the energy threshold. This is in sharp contrast to the prediction of Polanyi rules: translational energy is more efficient than vibrational energy in enhancing a reaction with an early barrier. The recent sudden vector projection model is then applied to rationalize the observed mode specificity, which, however, shows that the translational mode vector has a larger coupling with the reaction coordinate than the stretching vector of the active bond, implying a reversed relative efficacy on promoting the reaction as well. In contrast, the mode and bond specificity for the exchange channel is not as strong as for the abstraction channel due to the regulation of the shallow well along the reaction path.

Probing the Partial Activation of Water by Open-Shell Interactions, Cl(H2O)1-4

The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl^•(H₂O)_n=1-4. Using an automated isomer-search procedure, dozens of unique, stable structures were computed. Among the resulting structural classes were intact, hydrated-chlorine isomers, as well as hydrogen-abstracted (HCl)(OH)(H₂O)_n-1 configurations. The latter showed increased stability as the degree of hydration increased, until n = 4, where a new class of structures was discovered with a chloride ion bound to an oxidized water network. The electronic structure of these three structural classes was investigated, and spectral signatures of this hydration-based evolution were connected to these electronic properties. An ancillary outcome of this detailed computational analysis, including coupled-cluster benchmarks, was the calibration of cost-effective quantum chemistry methods for future studies of these radical-water complexes.

Revisiting the reaction energetics of the CH3O˙ + O2 (3Σ−) reaction: the crucial role of post-CCSD(T) corrections

The CH₃O˙ + O₂ reaction has been studied by means of high level ab initio calculations to predict the reaction energy and barrier height with chemical accuracy.

Kinetics of the OH+HCl→H2 O+Cl reaction: Rate determining roles of stereodynamics and roaming and of quantum tunneling

The OH + HCl → H₂ O + Cl reaction is one of the most studied four-body systems, extensively investigated by both experimental and theoretical approaches. Here, as a continuation of our previous work on the OH + HBr and OH + HI reactions, which manifest an anti-Arrhenius behavior that was explained by stereodynamic and roaming effects, we extend the strategy to understand the transition to the sub-Arrhenius behavior occurring for the HCl case. As previously, we perform first-principles on-the-fly Born-Oppenheimer molecular dynamics calculations, thermalized at four temperatures (50, 200, 350, and 500 K), but this time we also apply a high-level transition-state-theory, modified to account for tunneling conditions. We find that the theoretical rate constants calculated with Bell tunneling corrections are in good agreement with extensive experimental data available for this reaction in the ample temperature range: (i) simulations show that the roles of molecular orientation in promoting this reaction and of roaming in finding the favorable path are minor than in the HBr and HI cases, and (ii) dominating is the effect of quantum mechanical penetration through the energy barrier along the reaction path on the potential energy surface. The discussion of these results provides clarification of the origin on different non-Arrhenius mechanisms observed along this series of reactions. © 2018 Wiley Periodicals, Inc.

Quantum dynamics of ClH2O- photodetachment: Isotope effect and impact of anion vibrational excitation

Photodetachment of the ClH₂O^- anion is investigated using full-dimensional quantum mechanics on accurate potential energy surfaces of both the anion and neutral species. Detailed analysis of the photoelectron spectrum and the corresponding wavefunctions reveals that the photodetachment leads to, in the product channel of the exothermic HCl + OH → Cl + H₂O reaction, the formation of numerous Feshbach resonances due apparently to slow energy transfer from H₂O vibrational modes to the dissociation coordinate. These long-lived resonances can be grouped into two broad peaks in the low-resolution photoelectron spectrum, which is in good agreement with available experiments, and they are assigned to the ground and first excited OH stretching vibrational manifolds of H₂O complexed with Cl. In addition, effects of isotope substitution on the photoelectron spectrum were small. Finally, photodetachment of the vibrationally excited ClH₂O^- in the ionic hydrogen bond mode is found to lead to Feshbach resonances with higher stretching vibrational excitations in H₂O.

Accurate Determination of Tunneling-Affected Rate Coefficients: Theory Assessing Experiment

The thermal rate coefficients of a prototypical bimolecular reaction are determined on an accurate ab initio potential energy surface (PES) using ring polymer molecular dynamics (RPMD). It is shown that quantum effects such as tunneling and zero-point energy (ZPE) are of critical importance for the HCl + OH reaction at low temperatures, while the heavier deuterium substitution renders tunneling less facile in the DCl + OH reaction. The calculated RPMD rate coefficients are in excellent agreement with experimental data for the HCl + OH reaction in the entire temperature range of 200-1000 K, confirming the accuracy of the PES. On the other hand, the RPMD rate coefficients for the DCl + OH reaction agree with some, but not all, experimental values. The self-consistency of the theoretical results thus allows a quality assessment of the experimental data.

A global coupled cluster potential energy surface for HCl + OH ↔ Cl + H2O

A new and more accurate full-dimensional global potential energy surface (PES) for the ground electronic state of the ClH₂O system is developed by using the permutation invariant polynomial-neural network (PIP-NN) method to fit 15 777 points obtained using an explicitly correlated unrestricted coupled-cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b).

Quasi-classical trajectory studies on the full-dimensional accurate potential energy surface for the OH + H2O = H2O + OH reaction

The first accurate PES for the OH + H₂O reaction is developed by using the permutation invariant polynomial-neural network method to fit ∼48 000 CCSD(T)-F12a/AVTZ calculated points.

Insights into the bond-selective reaction of Cl + HOD(nOH) → HCl + OD

The theoretical study of the title reaction reveals strong bond selectivity, vibrational enhancement, and unique product distributions, in agreement with the experiment.

Energy disposal and thermal rate constants for the OH + HBr and OH + DBr reactions: quasiclassical trajectory calculations on an accurate potential energy surface

We report reaction cross sections, energy disposal, and rate constants for the OH + HBr → Br + H2O and OH + DBr → Br + HDO reactions from quasiclassical trajectory calculations using an ab initio potential energy surface [ de Oliveira-Filho , A. G. S. ; Ornellas , F. R. ; Bowman , J. M. J. Phys. Chem. Lett. 2014 , 5 , 706 - 712 ]. Comparison with available experiments are made and generally show good agreement.

The exothermic HCl + OH·(H2O) reaction: removal of the HCl + OH barrier by a single water molecule

The entrance complex, transition state, and exit complex for the title reaction have been investigated using the CCSD(T) method with correlation consistent basis sets up to cc-pVQZ. The stationary point geometries for the reaction are related to but different from those for the water monomer reaction HCl + OH → Cl + H2O. Our most important conclusion is that the hydrogen-bonded water molecule removes the classical barrier entirely. For the endothermic reverse reaction Cl + (H2O)2, the second water molecule lowers the relative energies of the entrance complex, transition state, and exit complex by about 4 kcal/mol. The title reaction is exothermic by 17.7 kcal/mol. The entrance complex HCl⋯OH·(H2O) is bound by 6.9 kcal/mol relative to the separated reactants. The classical barrier height for the reverse reaction is predicted to be 16.5 kcal/mol. The exit complex Cl⋯(H2O)2 is found to lie 6.8 kcal/mol below the separated products. The potential energy surface for the Cl + (H2O)2 reaction is radically different from that for the valence isoelectronic F + (H2O)2 system.

Quasiclassical Trajectory Calculations of the Rate Constant of the OH + HBr → Br + H2O Reaction Using a Full-Dimensional Ab Initio Potential Energy Surface Over the Temperature Range 5 to 500 K

We report a permutationally invariant, ab initio potential energy surface (PES) for the OH + HBr → Br + H2O reaction. The PES is a fit to roughly 26 000 spin-free UCCSD(T)/cc-pVDZ-F12a energies and has no classical barrier to reaction. It is used in quasiclassical trajectory calculations with a focus on the thermal rate constant, k(T), over the temperature range 5 to 500 K. Comparisons with available experimental data over the temperature range 23 to 416 K are made using three approaches to treat the OH rotational and associated electronic partition function. All display an inverse temperature dependence of k(T) below roughly 160 K and a nearly constant temperature dependence above 160 K, in agreement with experiment. The calculated rate constant with no treatment of spin-orbit coupling is overall in the best agreement with experiment, being (probably fortuitously) within 20% of it.

A nine-dimensional global potential energy surface for NH4(X2A1) and kinetics studies on the H + NH3 ↔ H2 + NH2 reaction

A nine-dimensional global potential energy surface (PES) for the NH₄ system is developed from ∼10⁵ high-level ab initio points and the hydrogen abstraction kinetics on the PES agree with experiment.