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.     

Exact quantum scattering study of the H + HS reaction on a new ab initio potential energy surface H2S (3A")

We present an exact quantum dynamical study and quasi-classical trajectory (QCT) calculations for the exchange and abstraction processes for the H + HS reaction. These calculations were based on a newly constructed high-quality potential energy surface for the lowest triplet state of H(2)S ((3)A"). The ab initio single-point energies were computed using complete active space self-consistent field and multi-reference configuration interaction method with a basis set of aug-cc-pV5Z. The time-dependent wave packet (TDWP) method was used to calculate the total reaction probabilities and integral cross sections over the collision energy (E(col)) range of 0.0-2.0 eV for the reactant HS initially at the ground state and the first vibrationally excited state. It was found that the initial vibrational excitation of HS enhances both abstraction and exchange processes. In addition, a good agreement is found between QCT and TDWP reaction probabilities at the total momentum J = 0 as a function of collision energy for the H + HS (v = 0, j = 0) reaction.     

Time-dependent wave packet quantum and quasi-classical trajectory study of He + H₂⁺, D₂⁺ → HeH⁺ + H, HeD⁺ + D reaction on an accurate FCI potential energy surface

The quantum scattering dynamics and quasi-classical trajectory (QCT) calculations have been carried out for the title reaction on an accurate potential energy surface (PES) computed using the full configuration interaction (FCI). On the basis of the PES, the integral cross-sections of He + H?? (v = 0-3, j = 1) → HeH? + H reaction have been calculated, and the results are generally agreed with the experimental cross-sections obtained by Tang et al. [J. Chem. Phys. 2005, 122, 164301] after taking into account the experimental uncertainties, which proves the reliability of implementing dynamics calculations on the FCI PES. The reaction probability of He + D?? (v = 0-2, j = 0) → HeD? + D reactions for total angular momentum J = 0 and the integral cross-section (ICS) have been calculated. The significant quantum effect has been explored by the comparison between the QCT reaction probabilities (or ICS) and the quantum mechanical (QM) reaction probabilities (or ICS), which may be attributed to the deep well in the PES of this light atoms system. Furthermore, the role of Coriolis coupling (CC) effects has also been found not important by the comparison between the CC calculation and the centrifugal sudden (CS) approximation calculation, except that the CC total cross-sections for the v = 1 and 2 states show the collision energy-dependent behaviors in the low-energy area, which are different from those based on the CS calculation.     

Quasi-classical trajectory study of the adiabatic reactions occurring on the two lowest-lying electronic states of the LiH2+ system

Quasi-classical trajectory calculations have been performed on the adiabatically allowed reactions taking place on the two lowest-lying electronic states of the LiH2+ system, using the ab initio potential energy surfaces of Martinazzo et al. (J. Chem. Phys., 2003, 119, 11 241). These reactions comprise: (i) the exoergic H2 and H2+ formation occurring through LiH+ + H and LiH + H+ collisions in the ground and in the first electronically excited state, respectively; (ii) the endoergic (ground state) LiH+ dissociation induced by collisions with H atoms; and (iii) the endoergic (excited state) Li + H2+ --> LiH + H+ reaction. The topic is of relevance for a better understanding of the lithium chemistry in the early universe. Thermal rate constants for the above reactions have been computed in the temperature range 10-5000 K and found in reasonably good agreement with estimates based on the capture model.     

Cumulative reaction probabilities: A comparison between quasiclassical and quantum mechanical results

This article presents a quasiclassical trajectory (QCT) method for determining the cumulative reaction probability (CRP) as a function of the total energy. The method proposed is based on a discrete sampling using integer values of the total and orbital angular momentum quantum numbers for each trajectory and on the development of equations that have a clear counterpart in the quantum mechanical (QM) case. The calculations comprise cumulative reaction probabilities at a given total angular momentum J, as well as those summed over J. The latter are used to compute QCT rate constants. The method is illustrated by comparing QCT and exact QM results for the H+H2, H+D2, D+H2, and H+HD reactions. The agreement between QCT and QM results is very good, with small discrepancies between the two data sets indicating some genuine quantum effects. The most important of these involves the value of the CRP at low energies which, due to the absence of tunneling, is lower in the QCT calculations, causing the corresponding rate constants to be smaller. The second is the steplike structure that is clearly displayed in the QM CRP for J = 0, which is much smoother in the corresponding QCT results. However, when the QCT density of reactive states, i.e., the derivatives of the QCT CRP with respect to the energy, is calculated, a succession of maxima and minima is obtained which roughly resembles those found in the QM calculations, although the latter are considerably sharper. The analysis of the broad peaks in the QCT density of reactive states indicates that the distributions of collision times associated with the maxima are somewhat broader, with a tail extending to larger collision times, than those associated with the minima. In addition, the QM and QCT dynamics of the isotopic variants mentioned above are compared in the light of their CRPs. Issues such as the compliance of the QCT CRP with the law of microscopic reversibility, as well as the similarity between the CRPs for ortho and para species in the QM and QCT cases, are also addressed.     

Accurate time-dependent wave packet study of the Li + H₂⁺ reaction and its isotopic variants

The dynamics and kinetics of the Li + H?? reaction and its isotopic variants (D?? and T??) have been studied by using a time-dependent wave packet (TDWP) coupled-channel (CC) method on the ab initio potential energy surface (PES) of Martinazzo et al. [J. Chem. Phys. 2003, 119, 21]. Total initial v = 0, j = 0 state-selected reaction probabilities for the Li + H?? reaction and its isotopic variants have been calculated from the threshold up to 1 eV for total angular momenta J from 0 to 90. Integral cross sections have been evaluated from the reaction probabilities at collision energies from threshold (≈0.2 eV) up to 1.0 eV collision. The calculated rate constants as a function of temperature show an Arrhenius type behavior in the 200 ≤ T ≤ 1000 K temperature interval. It has been found to be a considerable large intermolecular kinetic isotope effect. The TDWP-CC results are in overall good agreement with those obtained applying the TDWP Centrifugal-Sudden (CS) approximation, showing that the CS approximation is rather accurate for the title reaction.     

Capture and dissociation in the complex-forming CH + H2 → CH2 + H, CH + H2 reactions

The rate coefficients for the capture process CH + H(2)→ CH(3) and the reactions CH + H(2)→ CH(2) + H (abstraction), CH + H(2) (exchange) have been calculated in the 200-800 K temperature range, using the quasiclassical trajectory (QCT) method and the most recent global potential energy surface. The reactions, which are of interest in combustion and in astrochemistry, proceed via the formation of long-lived CH(3) collision complexes, and the three H atoms become equivalent. QCT rate coefficients for capture are in quite good agreement with experiments. However, an important zero point energy (ZPE) leakage problem occurs in the QCT calculations for the abstraction, exchange and inelastic exit channels. To account for this issue, a pragmatic but accurate approach has been applied, leading to a good agreement with experimental abstraction rate coefficients. Exchange rate coefficients have also been calculated using this approach. Finally, calculations employing QCT capture/phase space theory (PST) models have been carried out, leading to similar values for the abstraction rate coefficients as the QCT and previous quantum mechanical capture/PST methods. This suggests that QCT capture/PST models are a good alternative to the QCT method for this and similar systems.     

Influence of collision energy and reagent vibrational excitation on the dynamics of the reaction H + LiH

We present a detailed quasiclassical trajectory (QCT) study of the dynamics corresponding to the reaction H + LiH proceeding via depletion and H‐exchange paths on a new potential energy surface of the electronic ground state. The effects of collision energy and reagent initial vibrational excitation on the reaction probability and cross sections are studied over a wide range of collision energies. The QCT‐calculated reaction probability and cross sections are in good agreement with previous time‐dependent wave packet results. More importantly, we found that the vibrational excitation of LiH molecule inhibits the LiH depletion reaction, whereas it promotes the H‐exchange reaction. In addition, the differential cross sections calculated for the depletion reaction at different collision energies and excitation states indicate a strong forward scattering of the product molecule H₂. © 2013 Wiley Periodicals, Inc.     

Communication: rate coefficients from quasiclassical trajectory calculations from the reverse reaction: The Mu + H2 reaction re-visited

In a previous paper [P. G. Jambrina et al., J. Chem. Phys. 135, 034310 (2011)] various calculations of the rate coefficient for the Mu + H(2) → MuH + H reaction were presented and compared to experiment. The widely used standard quasiclassical trajectory (QCT) method was shown to overestimate the rate coefficients by several orders of magnitude over the temperature range 200-1000 K. This was attributed to a major failure of that method to describe the correct threshold for the reaction owing to the large difference in zero-point energies (ZPE) of the reactant H(2) and product MuH (～0.32 eV). In this Communication we show that by performing standard QCT calculations for the reverse reaction and then applying detailed balance, the resulting rate coefficient is in very good agreement with the other computational results that respect the ZPE, (as well as with the experiment) but which are more demanding computationally.     

Quantum wave packet calculation of reaction probabilities,cross sections,and rate constants for H+ + LiH → Li + H reaction

The H⁺ + LiH → Li + H reactive scattering has been studied using a quantum real wave packet method. The state‐to‐state and state‐to‐all reaction probabilities for the entitled collision have been calculated at zero total angular momentum. The probabilities for J > 0 are estimated from the J = 0 results by using J‐shifting approximation based on the Capture model. The integral cross sections and thermal rate constants are then calculated. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Time dependent quantum dynamics study of the Ne + H2+(v=0-4)-->NeH+ + H proton transfer reaction

The Ne + H2+-->NeH+ + H proton transfer reaction was studied using the time dependent real wave packet quantum dynamics method at the helicity decoupling level, considering the H2+ molecular ion in the (v=0-4, j=0) vibrorotational states and a wide collision energy interval. The calculated reaction probabilities and reaction cross sections were in a rather good agreement with reanalyzed previous exact quantum dynamics results, where a much smaller collision energy interval was considered. Also, a quite good agreement with experimental data was found. These results suggested the adequacy of the approach used here to describe this and related systems.     

Quantum wave-packet calculation of reaction probabilities, cross sections, and rate constants for Li + H2+ reaction

The Li + H2+(upsilon,j) --> LiH(upsilon',j') + H+ reactive scattering has been studied by using quantum real wave-packet method. The state-to-state and state-to-all reaction probabilities for the entitled collision have been calculated. The probabilities show a smooth variation for all initial rotational quantum states. The J-shifting approximation has been employed to estimate the integral cross sections and thermal rate constants have been calculated.     

A quantum wave packet dynamics study of the N(2D) + H2 reaction

We report a dynamics study of the reaction N((2)D) + H(2) (v=0, j=0-5) --> NH + H using the time-dependent quantum wave packet method and a recently reported single-sheeted double many-body expansion potential energy surface for NH(2)(1(2)A' ') which has been modeled from accurate ab initio multireference configuration-interaction calculations. The calculated probabilities for (v=0, j=0-5) are shown to display resonance structures, a feature also visible to some extent in the calculated total cross sections for (v=0, j=0). A comparison between the calculated centrifugal-sudden and coupled-channel reaction probabilities validate the former approximation for the title system. Rate constants calculated using a uniform J-shifting scheme and averaged over a Boltzmann distribution of rotational states are shown to be in good agreement with the available experimental values. Comparisons with other theoretical results are also made.     

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.     

The H + Li2 bimolecular exchange reaction: dynamical and kinetical properties at J = 0

For the first time in the literature, rigorous time-independent quantum scattering formalism was applied, by means of the ABC program, to the H + Li(2) → LiH + Li reaction. The state-to-state probabilities as a function of the total energy have been computed at zero total angular momentum (J = 0) allowing us to evaluate the effect of vibrational/rotational excitation on the reaction promotion/inhibition, the energetic distribution of products, and the temperature dependence of the J-shifting thermal rate coefficients.     

Time-dependent quantum study of the kinetics of the H(2S) + FO(2II) OH(2II) + F(2P) reaction

Quantum mechanical wave packet calculations are carried out for the H((2)S) + FO((2)II) --> OH((2)II) + F((2)P) reaction on the adiabatic potential energy surface of the ground (3)A' triplet state. The state-to-state and state-to-all reaction probabilities for total angular momentum J = 0 have been calculated. The probabilities for J > 0 have been estimated from the J = 0 results by using J-shifting approximation based on a capture model. Then, the integral cross sections and initial state-selected rate constants have been calculated. The calculations show that the initial state-selected reaction probabilities are dominated by many sharp peaks. The reaction cross section does not manifest any sharp oscillations and the initial state-selected rate constants are sensitive to the temperature.     

Real wave packet and flux analysis studies of the H + F2 → HF + F reaction

The H + F₂ → HF + F reaction on ground state potential energy surface is investigated using the quantum mechanical real wave packet and Flux analysis method based on centrifugal sudden approximation. The initial state selected reaction probabilities for total angular momentum J = 0 have been calculated by both methods while the probabilities for J > 0 have been calculated by Flux analysis method. The initial state selected reaction probabilities, integral cross sections and rate coefficients have been calculated for a broad range of collision energy. The results show a large rotational enhancement of the reaction probability. Some resonances were seen in the state‐to‐state reaction probabilities while state‐to‐all reaction probabilities and the reaction cross section do not manifest any oscillations and the initial state selected reaction rate constants are sensitive to the temperature. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

A state-to-state dynamical study of the Br + H(2) reaction: comparison of quantum and classical trajectory results

We present a detailed theoretical investigation of the dynamics corresponding to the strongly endothermic Br + H(2) (v = 0-1, j = 0) → H + HBr reaction in the 0.85 to 1.9 eV total energy range. State-averaged and state-to-state results obtained through time-independent wave packet (TIWP) and time-independent quantum mechanical (TIQM) calculations and quasiclassical trajectories (QCT) are compared and analyzed. The agreement in the results obtained with both quantum mechanical results is very good overall. However, although QCT calculations reproduce the general features, their agreement with the QM results is sometimes only qualitative. The analysis of the mechanism based on state-averaged results turns out to be deceptive and conveys an oversimplified picture of the reaction consistent with a direct-rebound mechanism. Consideration of state-to-state processes, in contrast, unveils the existence of multiple mechanisms that give rise to a succession of maxima in the differential cross section (DCS). Such mechanisms correlate with different sets of partial waves and display similar collision times when analyzed through the time-dependent DCS.     

Exact quantum and TST-CEQ study of the collinear reaction O+HCl

Exact quantum calculations of reaction probabilities have been carried out using hyperspherical coordinates for the collinearr reaction O+HCl(v <1) -OH(v'<1)+Cl . A generalized LEPS potential energy surface with a barrier height of 8.12 kcal/mol has been used in the calculations. According to the calculated results we found that (1) the reaction probability oscillates with energy, (2) the reaction probability shows vibrational adiabaticity, although it is poorer than that for symmetric reaction Cl + HC1. The analysis of resonance has also been done. The reaction rate constants and average cross sections have been calculated by TST-CEQ method. The rate constants are in agreement with that by QCT and smaller than the experimental one. Finally, the threshold has been estimated and is in good agreement with that of the literature.