Double-well structure and microscopic branching in the ()+ reaction

Theoretical study of product polarization of O(1D) + HCl(v = 0; j = 0) → ClO + H and its isotope exchange reaction

O(1D) + HCl(v = 0; j = 0) → ClO + H and its isotope exchange reaction O(1D) + DCl(v = 0; j = 0) → ClO + D are studied in the collision energy range 14.0–20.0 kcal/mol based on the potential energy surface 1[Formula: see text] state. Reaction probabilities, integral cross sections, the two angular distribution functions (concerning the initial/final velocity vector, and the product rotational momentum vector), and the product rotational alignment parameters are calculated as a function of the collision energy for the two reactions. The four generalized polarization dependent differential cross sections are presented to manifest the polarization characters. Also, the effect of the collision energy and the kinetic isotope effect are studied.

The dynamical study of O(1D) + HCl(v = 0, j = 0) reaction at hyperthermal collision energies

Backgrounds

The quasi-classical trajectory calculations for O(¹D) + HCl → OH + Cl (R1) and O(¹D) + HCl → ClO + H (R2) reactions have been performed at hyperthermal collision energies (60.0, 90.0, and 120.0 kal/mol) on the ¹A' state. Reaction probabilities and integral cross sections are calculated. The product rotational distributions for the two channels, and the product rotational alignment parameters are investigated. Also, the alignment and the orientation of the products have been predicted through the angular distribution functions (concerning the initial/final velocity vector, and the product rotational angular momentum vector). To have a deeper understanding of the natures of the vector correlation between reagent and product relative velocities, a natural generalization of the differential cross section __PDDCS₀₀, is calculated.

Results

The OH + Cl channel is the main product channel and is observed to have essentially isotropic rotational distributions. The ClO + H channel is found to be clearly rotationally polarized.

Conclusions

The dynamical, especially the stereodynamical characters are quite different for the two channels of the title reaction. Most reactions occur directly, except for R2 reaction at the collision energies of 60.0 and 120.0 kcal/mol. The alignment and orientation effects are weak/strong for R1/R2 reaction. The well structure on the potential energy surface and hyperthermal collision energies might result in the dynamical effects.

THEORETICAL STUDY OF THE STEREO-DYNAMICS OF THE REACTION O + HCl → ClO + H

Quasi-classical trajectory (QCT) method is used to study the stereo-dynamics of the title reaction on the ground 1 1A′ potential energy surface (PES). Differential cross-sections (DCSs) and alignments of the product rotational angular momentum for the reaction are reported. The influence of collision energy on the product vector properties is also studied in the present work. The distribution of angle between k and j′, P(θr), the distribution of dihedral angle denoting k-k′-j′ correlation, P(ϕr) ⋅ (2π/σ)( d σ00/ d ωt), (2π/σ)( d σ20/ d ωt), (2π/σ)( d σ22+/ d ωt) and (2π/σ)(dσ21-/dωt) have been calculated in the center of mass frame, respectively.

CHEMICAL REACTIONS IN THE O(1D) + HCl SYSTEM II

Using the accurate global potential energy surfaces for the 11A′, 11A′′, and 21A′ states reported in the previous sister Paper I, detailed quantum dynamics calculations are performed for these three adiabatic surfaces separately for J = 0 (J: total angular momentum quantum number). Overall reaction probabilities for O + HCl → OH + Cl and H + ClO, the branching ratio between the two reactions, effects of the initial rovibrational excitation, and product rovibrational distributions are evaluated in the total energy region E tot ≤ 0.9 eV. Significant contributions to the overall reaction dynamics are found from the two excited 11A′′ and 21A′ potential energy surfaces, clearly indicating the insufficiency of the dynamics only on the ground 11A′ surface. The detailed dynamics on the excited surfaces are reported in the third paper of this series.

EFFECTS OF ROTATIONAL AND VIBRATIONAL EXCITATION ON THE STEREODYNAMICS OF THE O (D-1) + HCl → OH + Cl REACTION

The stereodynamics of the title reaction on the ground 1 1A′ potential energy surface (PES) has been studied using quasi-classical trajectory (QCT) method. Collision energy of 6.4 kcal/mol is considered, and vector properties including angular momentum alignment distributions and polarization-dependent differential cross-sections (PDDCS) of the product OH are presented. Furthermore, the influence of reagent rotational excitation and vibrational excitation on the product vector properties has also been studied in the present work. The results indicate that the distribution of the P(θr) and P(ϕr) are sensitively affected by the rotational and vibrational excitation. The rotational excitation decreases the degree of alignment and orientation, while vibrational excitation increases the degree of alignment and orientation. The PDDCS (2π/σ)(dσ20/dωt) and (2π/σ)(dσ22+/dωt) are sensitively influenced by rotational and vibrational excitations, while the PDDCS ((2π/σ)(dσ00/dωt)) and (2π/σ)(dσ21-/dωt) are not. The preference of forward scattering has been found from the results of PDDCS ((2π/σ)(dσ00/dωt)), which is in good agreement with the experimental results.

ACCURATE AND HIGHLY EFFICIENT CALCULATION OF THE O(1D)HCl VIBRATIONAL BOUND STATES, USING A COMBINATION OF METHODS

Hypochlorous acid, HOCl, is an important intermediate in the O (1D) HCl reactive system. Due in part to a large number of vibrational bound states (over 800), extremely large direct product basis sets (around 300,000) are required to compute the energy levels just below the dissociation threshold. This situation, combined with a very high density of states, results in difficult convergence for iterative methods — e.g. Lanczos requires 50,000 iterations, and filter diagonalization uses 60,000 iterations. In contrast, using new methodologies, we are able to compute the highest-lying bound states with only 271 iterations, although the CPU cost per iteration is substantially greater. Lower lying states are also computed, for a fraction of the CPU cost of the highest energy calculation.

The dynamics of the O(1D) + HCl --> OH + Cl reaction at a 0.26 eV collision energy: a comparison between theory and experiment

The dynamics of the O((1)D) + HCl(v = 0, j = 0) --> Cl + OH reaction at a 0.26 eV collision energy has been investigated by means of a quasiclassical trajectory (QCT) and statistical quantum and quasiclassical methods. State-resolved cross sections and Cl atom velocity distributions have been calculated on two different potential energy surfaces (PESs): the H2 surface (Martinez et al. Phys. Chem. Chem. Phys. 2000, 2, 589) and the latest surface by Peterson, Bowman, and co-workers (PSB2) (J. Chem. Phys. 2000, 113, 6186). The comparison with recent experimental results reveals that the PSB2 PES manages to describe correctly differential cross sections and the velocity distributions of the departing Cl atom. The calculations on the H2 PES seem to overestimate the OH scattering in the forward direction and the fraction of Cl at high recoil velocities. Although the comparison of the corresponding angular distributions is not bad, significant deviations with a statistical description are found, thus ruling out a complex-forming mechanism as the dominant reaction pathway. However, for the ClO + H product channel, the QCT and statistical predictions are found to be in good agreement.

A Crossed Molecular Beam Imaging Study of the O(1D2)+HCl→OH+Cl(2PJ=3/2, 1/2) Reaction

A crossed molecular beam study is presented for the O((1)D(2))+HCl-->OH+Cl((2)P(J)) reaction at the collision energy of 6 kcal mol(-1). State-resolved doubly differential cross sections are obtained for the Cl((2)P(J=3/2) ) and Cl*((2)P(J=1/2) ) products by velocity-map ion imaging. Both products are slightly more forward scattered, which suggests a reaction mechanism without a long-lived intermediate in the ground electronic state. A small fraction (23 %) of the energy release into the translational degree of freedom indicates strong internal excitation of the counterpart OH radical. The contribution of the electronic excited states of O--HCl to the overall reaction is also examined from the doubly differential cross sections.

Accurate and highly efficient calculation of the highly excited pure OH stretching resonances of O(1D)HCl, using a combination of methods

Accurate calculation of the energies and widths of the resonances of HOCl--an important intermediate in the O(1D)HCl reactive system--poses a challenging benchmark for computational methods. The need for very large direct product basis sets, combined with an extremely high density of states, results in difficult convergence for iterative methods. A recent calculation of the highly excited OH stretch mode resonances using the filter diagonalization method, for example, required 462,000 basis functions, and 180,000 iterations. In contrast, using a combination of new methods, we are able to compute the same resonance states to higher accuracy with a basis less than half the size, using only a few hundred iterations-although the CPU cost per iteration is substantially greater. Similar performance enhancements are observed for calculations of the high-lying bound states, as reported in a previous paper [J. Theo. Comput. Chem. 2, 583 (2003)].