Quantum and classical studies of the O(3P) + H2(v = 0-3,j = 0) --> OH + H reaction using benchmark potential surfaces

We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.     

Dissociative chemisorption of H2 on the Cu(110) surface: a quantum and quasiclassical dynamical study

Six-dimensional quantum dynamical and quasiclassical trajectory (QCT) calculations are reported for the reaction and vibrationally inelastic scattering of (v = 0,1,j = 0) H(2) scattering from Cu(110), and for the reaction and rovibrationally elastic and inelastic scattering of (v = 1,j = 1) H(2) scattering from Cu(110). The dynamics results were obtained using a potential energy surface obtained with density functional theory using the PW91 functional. The reaction probabilities computed with quantum dynamics for (v = 0,1,j = 0) were in excellent agreement with the QCT results obtained earlier for these states, thereby validating the QCT approach to sticking of hydrogen on Cu(110). The vibrational de-excitation probability P(v=1,j = 0 --> v = 0) computed with the QCT method is in remarkably good agreement with the quantum dynamical results for normal incidence energies E(n) between 0.2 and 0.6 eV. The QCT result for the vibrational excitation probability P(v = 0,j = 0 --> v = 1) is likewise accurate for E(n) between 0.8 and 1 eV, but the QCT method overestimates vibrational excitation for lower E(n). The QCT method gives probabilities for rovibrationally (in)elastic scattering, P(v = 1,j = 1 --> v('),j(')), which are in remarkably good agreement with quantum dynamical results. The rotationally averaged, initial vibrational state-selective reaction probability obtained with QCT agrees well with the initial vibrational state-selective reaction probability extracted from molecular beam experiments for v = 1, for the range of collision energies for which the v=1 contribution to the measured total sticking probability dominates. The quantum dynamical probabilities for rovibrationally elastic scattering of (v = 1,j = 1) H(2) from Cu(110) are in good agreement with experiment for E(n) between 0.08 and 0.25 eV.     

Ion-pair formation observed in a pulsed-field ionization photoelectron spectroscopic study of HF

The pulsed-field ionization (PFI) photoelectron (PE) spectrum of HF has been recorded at the chemical dynamics beamline of the advanced light source over the photon energy range 15.9-16.5 eV using a time-of-flight selection scheme at a resolution of 0.6 meV. Rotationally-resolved structure in the HF+(X 2 pi 3/2, 1/2, v+ = 0, 1) band systems are assigned. The spectral appearance of these systems agrees with a previous VUV laser PFI-PE study. Importantly, extensive rotationally-resolved structure between these two vibrational band systems is also observed. This is attributed to ion-pair formation via Rydberg states converging on the v+ = 1 vibrational levels of the HF+(X 2 pi 3/2, 1/2) spin-orbit states. These Rydberg states are assigned to the 1 sigma+ part of the nd-complexes (sigma, pi, and delta). Ion-pair formation is observed in this study by the detection of F- ions. Some partially rotationally-resolved structure in a previously published threshold photoelectron spectrum is similarly attributed to ion-pair formation (F- detection) through a combination of the v+ = 17 level of the (A 2 sigma+) 3s sigma Rydberg state and the (X 2 pi 3/2, 1/2, v+ = 1) 7d Rydberg states. On the basis of the present study, an accurate experimental value for the dissociation energy of the ground state of HF has been obtained, D0(HF) = 5.8650(5) eV.     

Cold and ultracold chemical reactions of F+HCl and F+DCl

We report quantum dynamics calculations of F((2)P)+HCl(v,j)-->HF(v('),j('))+Cl((2)P) and F+DCl(v,j)-->DF(v('),j('))+Cl reactions at cold and ultracold temperatures. The effect of rotational and vibrational excitations of the HCl molecule on the reactivity is investigated. It is found that, in the ultracold regime, vibrational excitation of the HCl molecule from v=0 to v=2 enhances the reactivity by four orders of magnitude. The rotational excitation from j=0 to j=1 decreases the reactivity while the rotational excitation from j=0 to j=2 increases the reactivity. The overall effect of rotational excitation was found to be much smaller than vibrational excitation. The reactivity of the F+DCl system is significantly lower than that of the F+HCl case indicating the importance of quantum tunneling at low energies. For both reactions, Feshbach resonances corresponding to Fcdots, three dots, centered HCl or Fcdots, three dots, centeredDCl triatomic states occur at low energies. We also explored the validity of the coupled-states approximation for cold collisions taking the F+HCl(v=0,j=0) reaction as an illustrative example. It is found that the coupled-states approximation is generally valid for the background scattering even at low energies but it is inadequate to accurately describe the rich resonances in the energy dependence of the cross section resulting from the decay of van der Waals complexes. It is further shown that the coupled-states approximation cannot be used for scattering in the Wigner threshold regime when the molecule is initially in a rotationally excited level.     

An experimental and quasiclassical trajectory study of the rovibrationally state-selected reactions: HD+(v=0-15,j=1)+He-->HeH+(HeD+)+D(H)

The absolute integral cross sections for the formation of HeH+ and HeD+ from the collisions of HD+(v,j=1)+He have been examined over a broad range of vibrational energy levels v=0-13 at the center-of-mass collision energies (ET) of 0.6 and 1.4 eV using the vacuum ultraviolet (VUV) pulsed field ionization photoelectron secondary ion coincidence method. The ET dependencies of the integral cross sections for products HeH+ and HeD+ from HD+(v=0-4)+He collisions in the ET range of 0-3 eV have also been measured using the VUV photoionization guided ion beam mass spectrometric technique, in which vibrationally selected HD+(v) reactant ions were prepared via excitation of selected autoionization resonances of HD. At low total energies, a pronounced isotope effect is observed in absolute integral cross sections for the HeH++D and HeD++H channels with significant favoring of the deuteron transfer channel. As v is increased in the range of v=0-9, the integral cross sections of the HeH++D channel are found to approach those of HeD++H. The observed velocity distributions of products HeD+ and HeH+ are consistent with an impulsive or spectator-stripping mechanism. Detailed quasiclassical trajectory (QCT) calculations are also presented for HD+(v,j=1)+He collisions at the same energies of the experiment. The QCT calculations were performed on the most accurate ab initio potential energy surface available. If the zero-point energy of the reaction products is taken into account, the QCT cross sections for products HeH+ and HeD+ from HD+(v)+He are found to be significantly lower than the experimental results at ET values near the reaction thresholds. The agreement between the experimental and QCT cross sections improves with translational energy. Except for prethreshold reactivity, QCT calculations ignoring the zero-point energy in the products are generally in good agreement with experimental absolute cross sections. The experimental HeH+/HeD+ branching ratios for the HD+(v=0-9)+He collisions are generally consistent with QCT predictions. The observed isotope effects can be rationalized on the basis of differences in thermochemical thresholds and angular momentum conservation constraints.     

Molecular beam scattering studies of the reaction D + H2 (v = 0) and D + H2 (v = 1) → HD + H

The reaction D + H₂ → HD + H has been investigated in two molecular beam scattering experiments. Angular and time-of-flight distributions have been measured for the initial vibrational ground state (v = 0) at a most probable collision energy of E_cm = 1.5 eV and for the first vibrational excited state (v = 1) at E_cm = 0.28 eV with the same apparatus. Results for the ground-state experiment are compared with quasiclassical trajectory calculations(QCT) on the LSTH-hypersurface transformed into the laboratory system and averaged over the apparatus distributions. The agreement isquite satisfactory. At this high collision energy the HD products are no longer scattered in a backward direction but in a wide angular region concentrated about θ = 90° in the center-of-mass system. The absolute reactive cross section has been determined and the agreement with the theoretical value from QCT calculations is within the experimental error. The high sensitivity of the experiment to different properties of the doubly differential cross section has also been demonstrated. A preliminary evaluation of the experiment with initial vibrational excitation (v = 1) shows that the HD-product molecules are preferably backward scattered and the change of internal energy is small supporting the concept of a reaction which is adiabatic with respect to the internal degrees of freedom.     

Quantum state-to-state dynamics for the quenching process of Br(2P1/2) + H2(v(i) = 0, 1, j(i) = 0)

Quantum state-to-state dynamics for the quenching process Br((2)P(1/2)) + H(2)(v(i) = 0, 1, j(i) = 0) → Br((2)P(3/2)) + H(2)(v(f), j(f)) has been studied based on two-state model on the recent coupled potential energy surfaces. It was found that the quenching probabilities have some oscillatory structures due to the interference of reflected flux in the Br((2)P(1/2)) + H(2) and Br((2)P(3/2)) + H(2) channels by repulsive potential in the near-resonant electronic-to-vibrational energy transfer process. The final vibrational state resolved integral cross sections were found to be dominated by the quenching process Br((2)P(1/2)) + H(2)(v) → Br((2)P(3/2)) + H(2)(v+1) and the nonadiabatic reaction probabilities for Br((2)P(1/2)) + H(2)(v = 0, 1, j(i) = 0) are quite small, which are consistent with previous theoretical and experimental results. Our calculated total quenching rate constant for Br((2)P(1/2)) + H(2)(v(i) = 0, j(i) = 0) at room temperature is in good agreement with the available experimental data.     

Vibrational and translational enhancement of endothermic reactions: K + HCl(υ = 0,1) and K + HF(υ = 0,1)

The relative integral cross section for the two endothermic reactions K + HCl(υ = 0 and 1) → KCl + H and K + HF(υ = 0 and 1) → KF + H has been measured as a function of the collision energy E using the crossed molecular beam technique. The vibrationally excited state (υ = 1) has been populated thermally by heating the beam source to temperatures around 2000 K. The variation of the collision energy from thermal up to around 2.1 eV was achieved by seeding the K-beam with various carrier gases. The molecular reaction product was detected by surface ionization in connection with a time-of-flight method. The total energy threshold of the reactions has been found to be equal to or only slightly above the corresponding endothermicities. This suggests a vanishing or very low barrier crest on the potential energy hypersurfaces which is contradictory to recent theoretical results. The inclusion of tunneling in case of K + HF leads to a negligible rise of the barrier heights. The efficacy of translational and vibrational energy in promoting the reactive process has been directly compared over a wide range of collision energies. For K + HCl the vibrational enhancement of the reactivity descends with increasing E from approximately a factor of 10 at E = 0.08 eV to around unity for E ? 0.5 eV. The good agreement of this experimental result with phase space calculations suggests that the marked enhancements are predominantly caused by the long-range attraction between reagents in connection with an “early” barrier on the potential energy surface. In case of K + HF vibrational energy is by a factor of up to 380 more favourable in promoting the reaction than the same amount of translational energy. Again, with rising collision energy its efficacy decreases but promotes the reaction still by a factor of 70 at E = 1.7 eV. Since phase space theory fails here the reaction is certainly non-statistical and we conclude that the observed large efficacy of vibrational energy is due to a “late” barrier. The proposed barrier positions for the two systems are in accordance with theoretical results.     

Energy resonance and inverse population of the product vibrational states for the reaction F + H2(v = 0) → HF(v′) + H,LCAC‐SW theoretical quantum scattering study

LCAC‐SW (linear combination of arrangement channel‐scattering wavefunction) method was used to calculate collinear state‐to‐state reaction probabilities for the reaction F + H₂(v = 0) → HF(v′) + H on the 6SEC potential energy surface. The results show that reaction probabilities P₀₂ and P₀₃ [i. e., v′ = 2,3 for reaction F + H₂ (v = 0) + HF(v′) + H] are primary, the population of product vibrational states is inverse and the reaction probabilities are oscillatory with collision energies, i.e., there is energy resonance in this reaction, which agrees with a new experiment.     

The polarization dependent differential cross sections of the reactions:H+LiH~+(v=0,j=0)→H_2+Li~+ and H~++ LiH(v=0,j=0)→H_2~++Li

<正>Quasi-classical trajectory(QCT) calculations have been carried out to study the generalized polarization dependent differentialcross sections(PDDCSs) for the reactions H + LiH~+(v = 0,j = 0)→H_2 + Li~+ and H~+ + LiH(v = 0,j = 0)→H_2~+ + Li occurring onthe two lowest-lying electronic states of the LiH_2~+ system,using the ab initio potential energy surfaces(PESs) of Martinazzo et al.[3].Four PDDCSs,i.e.,(2π/σ)(dσ_(00)/dω_t),(2π/σ)(dσ_(20)/dω_t),(2π/σ)(dσ_(22+)/dω_t),(2π/σ)(dσ_(21-)/dω_t) have been discussed in detail.     

反应物振动激发对反应He+H+2→HeH++H影响的准经典轨线研究

基于PPA势能面(Palmieri, P.; Puzzarini, C.; Aquilanti, V. Mol. Phys., 2000, 98: 1835), 运用准经典轨线(quasi-classical trajectory)方法, 对反应He+H+2 (v=0-2, j=0)→HeH++H的立体动力学性质进行了理论研究. 结果表明, 反应物振动激发对反应的k-j'两矢量相关和k-k'-j'三矢量相关分布都产生了较大影响. 除此之外, 极化微分反应截面对振动量子数v也非常敏感.     

Isotopic branching in (He, HD+) collisions

A three-dimensional time-dependent quantum mechanical approach is used to calculate the reaction probability (P(R)) and the integral reaction cross section (sigma(R)) for both channels of the reaction He + HD+(v = 0, 1, 2, 3; j = 0) --> HeH(D)+ + D(H), over a range of translational energy (E(trans)) on two different ab initio potential energy surfaces (McLaughlin-Thompson-Joseph-Sathyamurthy and Palmieri et al.). The reaction probability plots as a function of translational energy exhibit several oscillations, which are characteristic of the system. The vibrational enhancement of the reaction probability and the integral reaction cross section values are reproduced qualitatively by our calculations, in accordance with the experimental results. The isotopic branching ratio for the reaction decreases in going from v = 0 to v = 1 and then becomes nearly v-independent in going from v = 1 to v =3 on both the surfaces.     

CH4+O(3P)→CH3+OH反应的准经典轨线研究

用准经典轨线方法研究了O（3P）与CH4的反应，计算结果表明，CH4（υ＝0，j＝0）与O（3P）的反应在低及高的碰撞参数下都是直接反应，无短寿命的碰撞复合物生成，产物OH以向后散射为主，基本上处于振转基态.CH4（υ＝1，j＝1）与O（3P）的反应在低及高的碰撞参数下反应机理不一样。在低碰撞参数下是直接反应，无短寿命的碰撞复合物生成，产物OH以向后散射为主，主要处于振动基态，转动基本上是冷的，但比高碰撞参数下的热.在高的碰撞参数下则生成短寿命的碰撞复合物，产物OH以向前散射为主，表现出明显的周边动力学反应的特征，主要处于振动激发态（υ＝1），但转动仍然是较冷的。     

Effect of reagent rotation on isotopic branching in (He, HD+) collisions

A three-dimensional time-dependent quantum mechanical wave packet approach is used to calculate reaction probability (P(R)) and integral reaction cross section (sigma(R)) values for both the channels of the reaction He + HD(+) (v = 1; j = 0, 1, 2, 3) --> HeH(D)(+) + D(H), over a range of translational energy (E(trans)) on the McLaughlin-Thompson-Joseph-Sathyamurthy (MTJS) potential energy surface using centrifugal sudden approximation for nonzero total angular momentum (J) values. The reaction probability plots as a function of translational energy for different J values exhibit several oscillations, which are characteristic of the system. It is shown that HeH(+) is preferred over HeD(+) for large J values and that HeD(+) is preferred over HeH(+) for small J values for all the rotational (j) states studied. The integral reaction cross section for both the channels and therefore the isotopic branching ratio for the reaction depend strongly on j in contrast to the marginal dependence shown by earlier QCT calculations. The computed results are in overall agreement with the available experimental results.     

Quantum dynamics of the Li + HF --> H + LiF reaction at ultralow temperatures

Quantum-mechanical calculations are reported for the Li+HF(v=0,1,j=0)-->H+LiF(v',j') bimolecular scattering process at low and ultralow temperatures. Calculations have been performed for zero total angular momentum using a recent high-accuracy potential-energy surface for the X2A' electronic ground state. For Li+HF(v=0,j=0), the reaction is dominated by resonances due to the decay of metastable states of the Li cdots,...F-H van der Waals complex. Assignment of these resonances has been carried out by calculating the eigenenergies of the quasibound states. We also find that while chemical reactivity is greatly enhanced by vibrational excitation, the resonances get mostly washed out in the reaction of vibrationally excited HF with Li atoms. In addition, we find that at low energies, the reaction is significantly suppressed due to the less-efficient tunneling of the relatively heavy fluorine atom.     

A pulsed-field ionization photoelectron secondary ion coincidence study of the H2+ (X,upsilon+=0-15,N+=1)+He proton transfer reaction

The endothermic proton transfer reaction, H2+(upsilon+)+He-->HeH+ + H(DeltaE=0.806 eV), is investigated over a broad range of reactant vibrational levels using high-resolution vacuum ultraviolet to prepare reactant ions either through excitation of autoionization resonances, or using the pulsed-field ionization-photoelectron-secondary ion coincidence (PFI-PESICO) approach. In the former case, the translational energy dependence of the integral reaction cross sections are measured for upsilon+=0-3 with high signal-to-noise using the guided-ion beam technique. PFI-PESICO cross sections are reported for upsilon+=1-15 and upsilon+=0-12 at center-of-mass collision energies of 0.6 and 3.1 eV, respectively. All ion reactant states selected by the PFI-PESICO scheme are in the N+=1 rotational level. The experimental cross sections are complemented with quasiclassical trajectory (QCT) calculations performed on the ab initio potential energy surface provided by Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. The QCT cross sections are significantly lower than the experimental results near threshold, consistent with important contributions due to resonances observed in quantum scattering studies. At total energies above 2 eV, the QCT calculations are in excellent agreement with the present results. PFI-PESICO time-of-flight (TOF) measurements are also reported for upsilon+=3 and 4 at a collision energy of 0.6 eV. The velocity inverted TOF spectra are consistent with the prevalence of a spectator-stripping mechanism.     

A full dimensional time-dependent wave packet study for the H4 four-center, collision induced dissociation, and single exchange reactions: reaction probabilities for J=0

A time-dependent initial state selected wave packet method has been developed to study the H2(v(1)=10-11,j1=0)+H2'(v2=0,j2=0)-->HH'+HH' four-center (4C) reaction, and two other competing reactions: the H2+H2'-->H+H+H2' collision induced dissociation (CID) and the H2+H2'-->H+HH'+H' single exchange (SE) reaction, in full six dimensions. Initial state-specific total reaction probabilities for these three competing reactions are presented for total angular momentum J=0 and the effects of reagent vibration on reactions are examined. It is found that (a) the CID process is the dominant process over the whole energy range considered in this study, but the 4C and SE processes also have non-negligible probabilities; (b) the SE process has a lower threshold energy than the 4C process, but the SE probability increases slower than the 4C probability as collision energy increases; (c) the vibrational excitation of H2(v1) is much more efficient than translational motion for promoting these processes, in particular to the CID process.     

Stereodynamics of the F + HD(v = 0, j = 1) reaction: direct vs. resonant mechanisms

The stereodynamics and mechanism of the F + HD(v = 0, j = 1) → HF (DF) + D (H) reactions have been thoroughly analysed at collision energies in the 0-160 meV range. Specifically, this study is focused on (i) the comparison between the stereodynamics of the collisions leading to HF and DF formation, and (ii) the stereodynamical fingerprints of the resonance that occurs at low collision energies in the HF channel and whose manifestation in the total cross section is greatly diminished for initial j > 0. While previous studies were limited to the analysis of integral cross sections (ICS), differential cross sections (DCS) and reaction probabilities, in the present work we have included the analysis of vectorial quantities such as the direction of the initial rotational angular momentum and internuclear axis, and their effect on reactivity. In particular, polarisation parameters (PP) and polarisation dependent differential cross sections (PDDCS), quantities that describe how the intrinsic HD rotational angular momentum and molecular axis polarisations contribute to reaction, are calculated and examined. The evolution of the PPs with the collision energy differs markedly between the two reaction channels. For the DF channel, the PP values are small and change very little in the energy range in which DF formation is appreciable. In contrast, rapid fluctuations in the magnitude and sign of the PPs are observed in the HF channel at low collision energies in and around the resonance. As the collision energy increases, direct (non-resonant) scattering prevails, and the various quantities are reasonably well accounted for by the QCT calculations, as in the case of the DF channel. The intrinsic directional information has been used to access the extent of control that can be achieved through polarisation of the HD molecule prior to collision. It was found that the same extrinsic preparation leads to very different outcomes on the HF channel DCS when the collision energy is close to the resonance. It is also shown that polarisation of the HD internuclear axis along the initial relative velocity enhances the effect of the resonance and allows its clear identification. Finally, the effect of different extrinsic preparations on the angle-velocity DCS is found to be strong, thus allowing considerable control of product angular distributions.