Seven-dimensional quantum dynamics study of the H+NH3-->H2+NH2 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H+NH3-->H2+NH2 reaction using a seven-dimensional model and an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values. The total reaction probabilities are calculated for the initial ground and seven excited states of NH3 with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants are calculated for the temperature range 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the total reaction probabilities are overall very small, (b) the symmetric and asymmetric NH stretch excitations enhance the reaction significantly and almost all of the excited energy deposited was used to reduce the reaction threshold, (c) the excitation of the umbrella and bending motion have a smaller contribution to the enhancement of reactivity, (d) the main contribution to the thermal rate constants is thought to come from the ground state at low temperatures and from the stretch excited states at high temperatures, and (e) the calculated thermal rate constants are three to ten times smaller than the experimental data and transition state theory results.     

Differential Cross Sections for the H+D2O→HD+OD Reaction: a Full Dimensional State-to-State Quantum Dynamics Study

The time-dependent wave-packet method was employed to calculate the first full-dimensional state-to-state differential cross sections (DCS) for the title reaction with D₂O in the ground and the first symmetric (100) and asymmetric stretching (001) excited states. The calculated DCSs for these three initial states are strongly backward peaked at low collision energies. With the increase of collision energy, these DCSs become increasingly broader with the peak position shifting gradually to a smaller angle, consistent with the fact that the title reaction is a direct reaction via an abstraction mechanism. It is found that the (100) and (001) states not only have roughly the same integral cross sections, but also have essentially identical DCS, which are very close to that for the ground state at the same total energy of reaction. The reaction produces a small fraction of OD in the v=1 state, with the population close to the relative reactivity between the ground and vibrationally excited states, therefore confirming the experimental result of Zare et al. and the local mode picture[J. Phys. Chem. 97 , 2204 (1993)]. Unexpectedly, the stretching excitation reduces the rotation excitation of product HD at the same total energy.     

The investigation of spin-orbit effect for the F((2)P)+HD reaction

In this paper, we employ the time-dependent quantum wave packet method to study the reaction of F((2)P(3/2), (2)P(1/2)) with HD on the Alexander-Stark-Werner potential energy surface. The reaction probabilities and total integral cross sections of the spin-orbit ground and excited states for the two possible products of the system are calculated. Because the reaction channel of the excited spin-orbit state is closed at the resonance energy, the resonance feature does not appear in the reaction probabilities and cross section for the F((2)P(1/2))+HD(v=j=0)-->HF+D reaction, in contrast with that found for the ground spin--orbit state. We also compare the average cross sections of the two possible products with the experimental measurement. The resonance peak in the present average cross section for the HF+D product is slightly larger than the experimental result, but much smaller than that of the single-state calculations on the potential energy surface of Stark and Werner. It seems that the spin--orbit coupling would play a relatively important role in this reaction. Moreover, the isotope effects of the ground and excited spin--orbit states and the reactivity of the two product channels from the excited spin--orbit state are presented.     

Accurate quantum wave packet calculations for the F + HCl → Cl + HF reaction on the ground 1(2)A' potential energy surface

We present converged exact quantum wave packet calculations of reaction probabilities, integral cross sections, and thermal rate coefficients for the title reaction. Calculations have been carried out on the ground 1(2)A' global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged wave packet reaction probabilities at selected values of the total angular momentum up to a partial wave of J = 140 with the HCl reagent initially selected in the v = 0, j = 0-16 rovibrational states have been obtained for the collision energy range from threshold up to 0.8 eV. The present calculations confirm an important enhancement of reactivity with rotational excitation of the HCl molecule. First, accurate integral cross sections and rate constants have been calculated and compared with the available experimental data.     

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.     

Fully Coriolis-coupled quantum studies of the H + O2 (upsilon i = 0-2, j i = 0,1) --> OH + O reaction on an accurate potential energy surface: integral cross sections and rate constants

We present accurate quantum calculations of the integral cross section and rate constant for the H + O2 --> OH + O combustion reaction on a recently developed ab initio potential energy surface using parallelized time-dependent and Chebyshev wavepacket methods. Partial wave contributions up to J = 70 were computed with full Coriolis coupling, which enabled us to obtain the initial state-specified integral cross sections up to 2.0 eV of the collision energy and thermal rate constants up to 3000 K. The integral cross sections show a large reaction threshold due to the quantum endothermicity of the reaction, and they monotonically increase with the collision energy. As a result, the temperature dependence of the rate constant is of the Arrhenius type. In addition, it was found that reactivity is enhanced by reactant vibrational excitation. The calculated thermal rate constant shows a significant improvement over that obtained on the DMBE IV potential, but it still underestimates the experimental consensus.     

Differential reaction cross sections from rotationally resolved quantum scattering calculations: application to gas-phase S(N)2 reactions

Differential reaction cross sections have been computed based on previous rotationally resolved time-independent quantum-mechanical scattering calculations for the complex-forming S(N)2 reaction Cl(-) + CH(3)Br → ClCH(3) + Br(-). The results show almost isotropic cross sections for reactant molecules with high rotational quantum numbers. Backward scattering is disfavoured for reaction out of states with small rotational excitation, in particular the rovibrational ground state. This is a quantum-mechanical effect (interference of partial waves) that can partly be rationalized by simple classical arguments. In particular for higher vibrational excitations, an umbrella effect can be observed that favours the backward direction. It can be explained by the strong enhancement of the reactivity by opening a direct mechanism. The ion-dipole interaction exerts a torque onto the molecule which carries out a rotation by about 90° and then completes the reaction.     

State-to-State Quantum Dynamical Study of H+Br2→HBr+Br Reaction

The time-dependent wave packet method has been employed to calculate the state-to-state integral cross sections and differential cross sections (DCSs) for three initial states of the title reaction on the recently constructed neural network potential energy surface. It is found that the product HBr(\begin{document}$ v' $\end{document} = 2, 3, 4) states have the dominated population in the entire energy region considered here, indicating an inverted HBr vibrational state distribution. More than half of the available energy ends up as product internal motion, and most of which goes into the vibrational motion. Our calculations show that initial rotational excitation of Br\begin{document}$ _2 $\end{document} has little effect on the product ro-vibrational state distributions and DCSs of the reaction. While the initial vibrational excitation has some influences. The initial vibrational excitation to \begin{document}$ v_0 $\end{document} = 5 obviously enhance the product vibrational excitation in the low energy region. The DCSs for collision energy up to 0.5 eV at the ground and rotationally excited state are peaked in the backward direction, but the width of the angular distribution increases considerably with the increase of collision energy. For the vibrationally excited state, the DCSs are rather complicated with some strong forward scattering peaks for highly vibrationally excited products.     

Influence of collision energy and reagent rotation on the cross sections and product polarizations of the reaction F + HCl

Quasiclassical trajectory calculations have been carried out for the F+HCl reaction in three dimensions on a recent DHSN PES of the ground 1(2)A' electronic state [M. P. Deskevich, M. Y. Hayes, K. Takahashi, R. T. Skodje, and D. J. Nesbitt, J. Chem. Phys. 124, 224303 (2006)]. The effects of the collision energy and the reagent initial rotational excitation on the cross sections and product polarization are studied for the v = 0 and j ≤ 10 states of HCl over a wide collision energy range. It has been found that either the collision energy or the HCl rotational excitation increase remarkably reaction cross sections. The QCT-calculated integral cross sections are in good agreement with previous QM results. A detailed study on product polarization for the title reaction is also performed. The calculated results show that the product rotational angular momentum j' is not only aligned, but also oriented along the direction perpendicular to the scattering plane. The orientation of the HF product rotational angular momentum vector j' depends very sensitively on the collision energy and also affected by the reagent rotation. The theoretical findings and especially the roles of the collision energy and initial rotational momentum on the product polarization are discussed and reasonably explained by the HLH mass combination, the property of the PES, as well as the reactive mechanism.     

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.     

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.     

A transition state view on reactive scattering: initial state-selected reaction probabilities for the H + CH4 → H2 + CH3 reaction studied in full dimensionality

Initial state-selected reaction probabilities for the H+CH(4)→H(2)+CH(3) reaction are computed for vanishing total angular momentum by full-dimensional calculations employing the multiconfigurational time-dependent Hartree approach. An ensemble of wave packets completely describing reactivity for total energies up to 0.58 eV is constructed in the transition state region by diagonalization of the thermal flux operator. These wave packets are then propagated into the reactant asymptotic region to obtain the initial state-selected reaction probabilities. Reaction probabilities for reactants in all rotational states of the vibrational 1A(1), 1F(2), and 1E levels of methane are presented. Vibrational excitation is found to decrease reactivity when reaction probabilities at equivalent total energies are compared but to increase reaction probabilities when the comparison is done at the basis of equivalent collision energies. Only a fraction of the initial vibrational energy can be utilized to promote the reaction. The effect of rotational excitation on the reactivity differs depending on the initial vibrational state of methane. For the 1A(1) and 1F(2) vibrational states of methane, rotational excitation decreases the reaction probability even when comparing reaction probabilities at equivalent collision energies. In contrast, rotational energy is even more efficient than translational energy in increasing the reaction probability when the reaction starts from the 1E vibrational state of methane. All findings can be explained employing a transition state based interpretation of the reaction process.     

A full-dimensional time-dependent wave packet study of the OH + CO → H + CO<Subscript>2</Subscript> reaction

Full-dimensional time-dependent wave packet calculations were made to study the \(\hbox{OH}+\hbox{CO} \rightarrow \hbox{H}+\hbox{CO}_2\) reaction on the Lakin–Troya–Schatz–Harding potential energy surface. Because of the presence of deep wells supporting long-lived collision complex, one needs to propagate the wave packet up to 450,000 a.u. of time to fully converge the total reaction probabilities. Our calculation revealed that the CO bond was substantially excited vibrationally in the complex wells, making it necessary to include sufficient CO vibration basis functions to yield quantitatively accurate results for the reaction. We calculated the total reaction probabilities from the ground initial state and two vibrationally excited states for the total angular momentum J = 0. The total reaction probability for the ground initial state is quite small in magnitude with many narrow and overlapping resonances due to the small complex-formation reaction probability and small probability for complex decaying into product channel. Initial OH vibrational excitation considerably enhances the reactivity because it enhances the probability for complex decaying into product channel, while initial CO excitation has little effects on the reactivity. We also calculated the reaction probabilities for a number of J > 0 states by using the centrifugal sudden approximation. By doing some calculations with multiple K-blocks included, we found that the centrifugal sudden approximation can be employed to calculate the rate constant for the reaction rather accurately. The calculated rate constants only agree with experimental measurements qualitatively, suggesting more theoretical studies be carried out for this prototypical complex-formation four-atom reaction.     

Kinetics and dynamics of the NH3 + H → NH2 + H2 reaction using transition state methods, quasi-classical trajectories, and quantum-mechanical scattering

On a recent analytical potential energy surface developed by two of the authors, an exhaustive kinetics study, using variational transition state theory with multidimensional tunneling effect, and dynamics study, using both quasi-classical trajectory and full-dimensional quantum scattering methods, was carried out to understand the reactivity of the NH(3) + H → NH(2) + H(2) gas-phase reaction. Initial state-selected time-dependent wave packet calculations using a full-dimensional model were performed, where the total reaction probabilities were calculated for the initial ground vibrational state and for four excited vibrational states of ammonia. Thermal rate constants were calculated for the temperature range 200-2000 K using the three methods and compared with available experimental data. We found that (a) the total reaction probabilities are very small, (b) the symmetric and asymmetric N-H stretch excitations enhance the reactivity, (c) the quantum-mechanical calculated thermal rate constants are about one order of magnitude smaller than the transition state theory results, which reproduce the experimental evidence, and (d) quasi-classical trajectory calculations, which were performed with the main goal of analyzing the influence of the zero-point energy problem on the final dynamics results, reproduce the quantum scattering calculations on the same surface.     

Evidence for excited spin-orbit state reaction dynamics in F+H2: theory and experiment

We describe fully quantum, time-independent scattering calculations of the F+H2-->HF+H reaction, concentrating on the HF product rotational distributions in v'=3. The calculations involved two new sets of ab initio potential energy surfaces, based on large basis set, multireference configuration-interaction calculations, which are further scaled to reproduce the experimental exoergicity of the reaction. In addition, the spin-orbit, Coriolis, and electrostatic couplings between the three quasidiabatic F+H2 electronic states are included. The calculated integral cross sections are compared with the results of molecular beam experiments. At low collision energies, a significant fraction of the reaction is due to Born-Oppenheimer forbidden, but energetically allowed reaction of F in its excited (2P 1/2) spin-orbit state. As the collision energy increases, the Born-Oppenheimer allowed reaction of F in its ground (2P 3/2) spin-orbit state rapidly dominates. Overall, the calculations agree reasonably well with the experiment, although there remains some disagreement with respect to the degree of rotational excitation of the HF(v'=3) products as well as with the energy dependence of the reactive cross sections at the lowest collision energies.     

Quantum Dynamics of Li+HF/DF Reaction Investigated by a State-to-State Time-dependent Wave Packet Approach

Using the reactant coordinate based time-dependent wave packet method, on the APW potential energy surface, the differential and integral cross sections of the Li+DF/HF(v=0, j=0, 1) reactions were calculated over the collision energy range from the threshold to 0.25 eV. The initial state-specified reaction rate constants of the title reaction were also calculated. The results indicate that, compared with the Li+DF reaction, the product LiF of Li+HF reaction is a little more rotationally excited but essentially similar. The initial rotational excitation from j=0 to 1 has little effect on the Li+DF reaction. However, the rotational excitation of DF does result in a little more rotationally excited product LiF. The different cross section of both reactions is forward biased in the studied collision energy range, especially at relatively high collision energy. The resonances in the Li+HF reaction may be identifiable as the oscillations in the product ro-vibrational state-resolved integral cross sections and backward scattering as a function of collusion energy. For the Li+HF reaction, the rate constant is not sensitive to the temperature and almost has no change in the temperature range considered. For the Li+DF reaction, the rate constant increase by a factor of about 10 in the temperature range of 100?300 K. Brief comparison for the total reaction probabilities and integral cross section of the Li+HF reaction has been carried out between ours and the values reported previously. The agreement is good, and the difference should come from the better convergence of our present calculations.     

D CH_4反应的SVRT含时波包理论研究

基于Jordan和Gilbert势能面,用SVRT(semirigidvibratingrotortarget)模型,对D CH4反应进行了含时波包动力学研究,计算得到了不同初始振动转动态的总反应几率、积分散射截面和热速率常数.计算结果与H CH4反应进行了比较和讨论.反应几率随平动能的变化曲线呈现出显著的量子共振特性.通过对j=0时,v=0、1、2的反应几率的计算,看出H-CH3的振动激发极大地提高了反应几率,而反应阈能明显降低,说明反应分子的振动能对分子的碰撞反应有重要贡献.对v=0时,j=0、1、2、3的反应几率的计算,得出转动量子数j的增大也会使反应几率有较大的提高,但反应阈能基本不变.     

Quantum wave packet study of the H+ + D2 reaction on diabatic potential energy surfaces

The exact three-dimensional nonadiabatic quantum dynamics calculations were carried out for the title reaction by a time-dependent wave packet approach based on a newly constructed diabatic potential energy surface (Kamisaka et al. J. Chem. Phys. 2002, 116, 654). Three processes including those of reactive charge transfer, nonreactive charge transfer, and reactive noncharge transfer were investigated to determine the initial state-resolved probabilities and reactive cross sections. The results show that a large number of resonances can be observed in the calculated probabilities due to the deep well on adiabatic ground surface and the dominant process is the reactive noncharge-transfer process. Some interesting dynamical features such as v-dependent and j-dependent behaviors of the probabilities are also revealed. In addition, a good agreement has been achieved in the comparison between the calculated quantum cross sections from the ground rovibrational initial state and the experimental measurement data.     

Fully converged integral cross sections of collision induced dissociation, four-center, and single exchange reactions, and accuracy of the centrifugal sudden approximation in H2 + D2 reaction

The initial state selected time-dependent wave packet method was employed to calculate the integral cross sections for the H(2) + D(2) reaction with and without the centrifugal sudden (CS) approximation by including all important K (the projection of the total angular momentum on the body-fixed axis) blocks. With a full-dimensional model, the first fully converged coupled-channel (CC) cross sections for different competitive processes from the ground rotational state were obtained: collision induced dissociation (CID), four-center (4C) reaction and single exchange (SE) reaction. The effect of the total angular momentum J on the reaction dynamics of H(2) + D(2) and the accuracy of the CS approximation have also been studied. It was found that the CID and SE processes occur in a wide range of J values while the 4C process can only take place in a narrow window of J values. For this reason, the CC cross section for the 4C channel is merely comparable to the SE channel. A comparison of the integral cross sections from CC and CS calculations showed that the CS approximation works well for the CID process but not for the 4C and SE processes, and the discrepancy between the CC and CS cross sections grows larger as the translational energy and/or the vibrational energy increase(s).     

Resonant vibrational excitation and de-excitation of N2(v) by low-energy electrons

We have calculated cross sections and rate coefficients for low-energy electron impact excitation of the nitrogen molecule from vibrationally excited levels N2(v) 1-8. Calculations are performed in the 2Pig shape resonance energy region, from 0 to 5 eV. The cross sections are determined by using our recent integral cross section measurements of the ground level vibrational excitation and the most recent cross sections for elastic electron scattering, applying the principle of detailed balance. The rate coefficient calculations are performed for the Maxwellian electron energy distribution. By using extended Monte Carlo simulations, the electron energy distribution functions (EEDF) and the rate coefficients are also determined for the nonequilibrium conditions, in the presence of the homogeneous external electric field for the typical, moderate values of the electric field over gas number density ratios, E/N.