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.     

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.     

Full dimensional time-dependent quantum dynamics study of the H + NH3 --> H2 + NH2 reaction

A rigorous full dimensional time-dependent wave packet method has been developed for the reactive scattering between an atom and a tetra-atomic molecule. The method has been applied to the hydrogen abstraction reaction H+NH(3)-->H(2)+NH(2). Initial state-selected total reaction probabilities are investigated for the reactions from the ground vibrational state and from four excited vibrational states of ammonia. The total reaction probabilities from two lowest "tunneling doublets" due to the inversion barrier for the umbrella bending motion of NH(3) and from two pairs of doubly degenerate vibrational states of NH(3) are also inspected. Integral cross sections and rate constants are calculated for the reaction from the ground state with the centrifugal-sudden approximation. The calculated results are compared with those from the previous seven dimensional calculations [M. Yang and J. C. Corchado, J. Chem. Phys. 126, 214312 (2007)]. This work shows that the full dimensional rate constants are a factor of 3 larger than the corresponding seven dimensional calculated values at T=200 K and are overall smaller than those obtained from the variational transition state theory in the whole temperature region. The work also reveals that nonreactive NH bonds of NH(3) cannot be treated as spectators due to the fact that three NH bonds are coupled with each other during the reaction process.     

An eight-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction on a new modified potential energy surface

An eight-dimensional time-dependent quantum dynamics wave packet approach is performed for the study of the H2+C2H-->H+C2H2 reaction system on a new modified potential energy surface (PES) [L.-P. Ju et al., Chem. Phys. Lett. 409, 249 (2005)]. This new potential energy surface is obtained by modifying Wang and Bowman's old PES [J. Chem. Phys. 101, 8646 (1994)] based on the new ab initio calculation. This new modified PES has a much lower transition state barrier height at 2.29 kcal/mol than Wang and Bowman's old PES at 4.3 kcal/mol. This study shows that the reactivity for this diatom-triatom reaction system is enhanced by vibrational excitations of H2, whereas the vibrational excitations of C2H only have a small effect on the reactivity. Furthermore, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. The comparison of the rate constant between this calculation and experimental results agrees with each other very well. This comparison indicates that the new modified PES corrects the large barrier height problem in Wang and Bowman's old PES.     

A quantum reaction dynamics study of the translational, vibrational, and rotational motion effects on the HD + H3+ reaction

Time-dependent, quantum reaction dynamics wavepacket approach is employed to investigate the impacts of the translational, vibrational, and rotational motion on the HD+H(3)(+) → H(2)D(+) + H(2) reaction using the Xie-Braams-Bowman potential energy surface [Z. Xie, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 122, 224307 (2005)]. We treat this five atom reaction with a seven-degree-of-freedom model by fixing one Jacobi and one torsion angle related to H(3) (+) at the lowest saddle point geometry of the potential energy surface. The initial state selected reaction probabilities show that the rotational excitations of H(+)-H(2) greatly enhance the reactivity with the reaction probabilities increased double at high rotational states compared to the ground state. However, the vibrational excitations of H(3) (+) hinder the reactivity. The ground state reaction probability shows no reaction threshold for this exoergic reaction, and as the translational energy increases, the reaction probability decreases. Furthermore, reactive resonances and zero point energy play very important roles on the reaction dynamics. The obtained integral cross section has the character of an exoergic reaction without a threshold: it decreases with the translational energy increasing. The calculated thermal rate constants using this seven-degree-of-freedom model are in agreement with a later experiment measurement.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

Quantum dynamics of H + LiH reaction and its isotopic variants

Time-dependent quantum wave packet dynamics study is carried out to investigate the initial state selected channel specific reactivity of H + LiH collisional system on a new and more accurate ab initio potential energy surface developed by Wernli et al. [J. Phys. Chem. A 113, 1121 (2009)]. The H + LiH reaction proceeds through LiH depletion and H-exchange paths. While the former path is highly exoergic (by ～2.258 eV), the latter path is thermoneutral. State selected and energy resolved integral reaction cross sections and thermal rate constants are reported and compared with the literature data. The reactivity of the LiH depletion channel is found to be greater than the H-exchange channel. Rotational excitation of the reagent LiH molecule causes a decrease of reactivity of both the channels. On the other hand, the vibrational excitation of the reagent LiH decreases the reactivity of the LiH depletion channel and increases the reactivity of the H-exchange channel. The effect of isotopic substitution (H by D) on the reaction dynamics is also examined.     

Quantum wave packet dynamics of N(2D)+H2 reaction

The quantum wave packet dynamics of the title reaction within the coupled state approximation is examined here and initial state-selected reaction probabilities, integral reaction cross sections, and thermal rate constants are reported. The ab initio potential energy surface of the electronic ground state (1(2)A(")) of the system recently reported by Ho et al. [J. Chem. Phys., 119, 3063 (2003)] is employed in this investigation. All partial wave contributions up to the total angular momentum J=55 were necessary to obtain converged integral reaction cross sections up to a collision energy of 1.0 eV. Thermal rate constants are calculated from the reaction cross sections and compared with the available theoretical and experimental results. Typical resonances formed during the course of the reaction and elucidating the insertion type mechanism for the product formation are calculated. Vibrational energy levels supported by the deep well (approximately 5.5 eV) of the 1(2)A(") potential energy surface of NH(2) are also calculated for the total angular momentum J=0. A statistical analysis of the spacing between the adjacent levels of this energy spectrum is performed and the extent of irregularity in the spectral sequence is assessed.     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.     

Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces

The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.     

A reduced dimensionality quasiclassical and quantum study of the proton transfer reaction H3O(+)+H2O-->H2O+H3O(+)

We report quantum and quasiclassical calculations of proton transfer in the reaction H(3)O(+)+H(2)O in three degrees of freedom, the two OH(+) bond lengths and the OH(+)O angle. The reduced dimensional potential energy surface is obtained from the full dimensional OSS3(p) energy function of H(5)O(2) (+) [L. Ojamae, I. Shavitt, and S. J. Singer, J. Chem. Phys. 109, 5547 (1998)], with an additional long-range correction to reproduce the correct ion-molecule interaction. This surface is used to perform both quasiclassical trajectory and quantum reactive scattering calculations of the zero total angular momentum cumulative reaction probability and cross sections for initial rotational states 0, 1, and 2. Comparison of these quantities are made to assess the importance of quantum effects in this reduced dimensional reaction. Additional quasiclassical cross sections are calculated to obtain the thermal rate constant for the reaction.     

Cross sections of the O+ + H2 --> OH+ + H ion-molecule reaction and isotopic variants (D2, HD): quasiclassical trajectory study and comparison with experiments

A dynamics study [cross section and microscopic mechanism versus collision energy (E(T))] of the reaction O+ + H2 --> OH+ + H, which plays an important role in Earth's ionosphere and interstellar chemistry, was conducted using the quasiclassical trajectory method, employing an analytical potential energy surface (PES) recently derived by our group [R. Martinez et al., J. Chem. Phys. 120, 4705 (2004)]. Experimental excitation functions for the title reaction, as well as its isotopic variants with D2 and HD, were near-quantitatively reproduced in the calculations in the very broad collision energy range explored (E(T) = 0.01-6.0 eV). Intramolecular and intermolecular isotopic effects were also examined, yielding data in good agreement with experimental results. The reaction occurs via two microscopic mechanisms (direct and nondirect abstraction). The results were satisfactorily interpreted based on the reaction probability and the maximum impact parameter dependences with E(T), and considering the influence of the collinear [OHH]+ absolute minimum of the PES on the evolution from reactants to products. The agreement between theory and experiment suggests that the reaction mainly occurs through the lowest energy PES and nonadiabatic processes are not very important in the wide collision energy range analyzed. Hence, the PES used to describe this reaction is suitable for both kinetics and dynamics studies.     

State-to-state reactive differential cross sections for the H+H2-->H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.     

State-to-state reaction probabilities for the H+O2(v,j)-->O+OH(v',j') reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.     

A full dimensional, nine-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction

A full dimensional, nine-degree-of-freedom (9DOF), time-dependent quantum dynamics wave packet approach is presented for the study of the H2+C2H-->H+C2H2 reaction system. This is the first full dimensional quantum dynamics study for a diatom-triatom reaction system. The effects of the initial vibrational and rotational excitations of the reactants on the reactivity of this reaction are investigated. This study shows that vibrational excitations of H2 enhance the reactivity; whereas, the vibrational excitations of C2H only have a small effect on the reaction probability. In addition, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. Comparison of the ground state reaction probabilities of the 9DOF and 8DOF shows the reaction probability from the full dimensional calculation is larger, with more prominent resonance features.     

An eight-degree-of-freedom quantum dynamics study for the H2+C2H system

An eight-degree-of-freedom (8DOF) time-dependent wave-packet approach has been developed to study the H(2)+C(2)H-->H+C(2)H(2) reaction system. The 8DOF model is obtained by fixing one of the Jacobi torsion angle in the nine-degree-of-freedom AB+CDE reaction system. This study is an extension of the previous seven-degree-of-freedom (7DOF) computation [J. Chem. Phys. 119, 12057 (2003)] of this reaction system. This study shows that vibrational excitations of H(2) enhance the reaction probability, whereas the stretching vibrational excitations of C(2)H have only a small effect on the reactivity. Furthermore, the bending excitation of C(2)H, compared to the ground-state reaction probability, hinders the reactivity. A comparison of the rate constant between the 7DOF calculation and the present 8DOF results has been made. The theoretical and experimental results agree with each other very well when the present 8DOF results are adjusted to account for the lower transition state barrier heights found in recent ab initio calculations.     

Quantum mechanical and quasiclassical trajectory scattering calculations for the C(1D) + H2 reaction on the second excited 1 1A" potential energy surface

Time-independent quantum mechanical (QM) and quasiclassical trajectory (QCT) scattering calculations have been carried out for the C(1D) + H2 --> CH + H reaction at a collision energy of 80 meV on a newly developed ab initio potential energy surface [B. Bussery-Honvault et al., Phys. Chem. Chem. Phys. 7, 1476 (2005)] of 1 1A" symmetry, corresponding to the second singlet state 1 1B1 of CH2. A general good agreement has been found between the QM and QCT rotational distributions and differential cross sections (DCSs). In both cases, DCSs are strongly peaked in the forward direction with a small contribution in the backward direction in contrast with those obtained on the 1 1A' surface, which are nearly symmetric. Rotational distributions obtained on the 1 1A" surface are somewhat colder than those calculated on the 1 1A' surface. The specific dynamics and the contribution of the 1 1A" surface to the overall reactivity of this system are discussed.     

A time-dependent wave-packet quantum scattering study of the reaction H2+(v = 0-2,4,6;j = 1) + He--> HeH+ + H

The quantum scattering dynamics calculation was carried out for the titled reaction in the collision energy range of 0.0-2.4 eV with reactant H(2) (+) in the rotational state j = 1 and vibrational states v = 0-2, 4, and 6. The present time-dependent wave-packet calculation takes into account the Coriolis coupling (CC) and uses the accurate ab initio potential-energy surface of Palmieri et al. [Mol. Phys. 98, 1835 (2000)]. The importance of including the CC quantum scattering calculation has been revealed by the comparison between the CC calculation and the previous coupled state (CS) calculation. The CC total cross sections for the v = 2, 4, and 6 states show collision energy-dependent behaviors different from those based on the CS calculation. Furthermore, the collision energy dependence of the total cross sections obtained in the present CC calculation only exhibits minor oscillations, indicating that the chance is slim for reactive resonances in total cross sections to survive through the partial-wave averaging. The magnitude and profile of the CC total cross sections for v = 0-2 in the collision energy range of 0.0-2.5 eV are found to be consistent with experimental cross sections obtained recently by Tang et al. [J. Chem. Phys. 122, 164301 (2005)] after taking into account the experimental uncertainties.     

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.