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.     

Real wave packet and quasiclassical trajectory studies of the H+ + LiH reaction

Time-dependent real wave packet (RWP) and quasiclassical trajectory (QCT) calculations have been carried out to study the H(+) + LiH reaction on the ab initio potential-energy surface of Martinazzo et al. [J. Chem. Phys., 2003, 119, 11241]. Total initial state-selected and final state-resolved reaction probabilities for the two possible reaction channels, H(2)(+) + Li and LiH + H(+), have been calculated for total angular momentum J=0 at a broad range of collision energies. Integral cross sections and thermal rate coefficients have been calculated using the QCT method and from the corresponding J=0 RWP reaction probabilities by means of a capture model. The calculated thermal rate coefficients are found to be nearly independent of temperature in the 100-500 K interval with a value of approximately 10(-9) cm(3) s(-1), which is in good agreement with estimates used in evolutionary models of early-Universe lithium chemistry. The RWP results are found to be in good agreement overall with the corresponding QCT calculations.     

Time‐dependent quantum study of H(2S) + FO(2Π) → OH(2Π) + F(2P) reaction on the 13A′ and 13A″ states

The dynamics of the H(²S) + FO(²Π) → OH(²Π) + F(²P) reaction on the adiabatic potential energy surface of the 1³A′ and 1³A″ states is investigated. The initial state selected reaction probabilities for total angular momentum J = 0 have been calculated by using the quantum mechanical real wave packet method. The integral cross sections and initial state selected reaction rate constants have been obtained from the corresponding J = 0 reaction probabilities by means of the simple J‐Shifting technique. The initial state‐selected reaction probabilities and reaction cross section do not manifest any sharp oscillations and the initial state selected reaction rate constants are sensitive to the temperature. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

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.     

Quasiclassical trajectory study of the C(1D) + H2 reaction and isotopomeric variants: kinetic isotope effect and CD/CH branching ratio

The recently proposed ab initio single-sheeted double many-body expansion potential energy for the methylene molecule has been used to perform quasiclassical trajectory (QCT) calculations for the title reaction. Thermal and initial state-specific (v = 0, j = 0) rate constants for the C((1)D) + H(2)/HD/D(2) reactions have been obtained over a wide range of temperatures. Cross sections for the reaction C((1)D) + H(2) and its deuterated isotopes have also been calculated, as well as the CD/CH branching ratios for the C((1)D) + HD reaction. It is found that the CD + H product channel in the C((1)D) + HD reaction is preferred relative to the CH + D channel. The estimated rate constants are predicted to be in the order k(H2) > k(HD) > k(D2) and the calculated cross sections and rate constants compared with available theoretical and experimental data.     

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.     

Quantum wave packet study of N(2D) + H2 reactive scattering

N(²D) + H₂ → NH + H reaction at zero total angular momentum is studied by using a time dependent quantum wave packet method. State‐to‐state and state‐to‐all reactive scattering probabilities for a broad range of energy are calculated. The probabilities show many sharp peaks that ascribed to reactive scattering resonances. The probabilities for J > 0 are estimated by using the J‐shifting method. The integral cross sections and thermal rate constants are then calculated. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

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.     

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     

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.     

Quantum scattering calculation of the O(1D)+HBr reaction

Three-dimensional time-dependent quantum wave packet calculation for the O((1)D)+HBr reaction has been carried out using an accurate ab initio global potential energy surface [K. A. Peterson, J. Chem. Phys. 113, 4598 (2000)]. The calculations show that the initial state-selected reaction probabilities are dominated by resonance structures, and the lifetime of the resonance is generally in the subpicosecond time scale. The energy dependence of the reaction cross section is computed, which manifests still resonance structures, and is a decreasing function of the translational energy. The thermal rate constants are also computed, which are nearly independent on the temperature. The calculation results are discussed and compared to similar reaction with deep well.     

Quantum wave packet calculation of reaction probabilities,cross sections,and rate constants for the C(1D) + HD reaction

The time‐dependent real wave packet method has been used to study the C(¹D) + HD reaction. The state‐to‐state and state‐to‐all reactive scattering probabilities for a broad range of energies are calculated at zero total angular momentum. The probabilities for J > 0 are estimated from accurately computed J = 0 probabilities by using the J‐shifting approximation. The integral cross sections for a large energy range, and thermal rate constants are calculated. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

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.     

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.     

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

基于Horst的势能面,用SVRT(SemirigidVibratingRotorTarget)方法对D+HCN反应进行了含时波包动力学研究,计算得到了不同初始振转态的总反应几率和积分反应截面,采用UniformJ-shifting方法得到该反应的热速率常数.计算结果与H+HCN反应进行了比和讨论.     

Coriolis-coupled wave packet dynamics of H + HLi reaction

We investigated the effect of Coriolis coupling (CC) on the initial state-selected dynamics of H+HLi reaction by a time-dependent wave packet (WP) approach. Exact quantum scattering calculations were obtained by a WP propagation method based on the Chebyshev polynomial scheme and ab initio potential energy surface of the reacting system. Partial wave contributions up to the total angular momentum J=30 were found to be necessary for the scattering of HLi in its vibrational and rotational ground state up to a collision energy approximately 0.75 eV. For each J value, the projection quantum number K was varied from 0 to min (J, K(max)), with K(max)=8 until J=20 and K(max)=4 for further higher J values. This is because further higher values of K do not have much effect on the dynamics and also because one wishes to maintain the large computational overhead for each calculation within the affordable limit. The initial state-selected integral reaction cross sections and thermal rate constants were calculated by summing up the contributions from all partial waves. These were compared with our previous results on the title system, obtained within the centrifugal sudden and J-shifting approximations, to demonstrate the impact of CC on the dynamics of this system.     

Reactions of C(1D) with H2 and its deuterated isotopomers, a wave packet study

Using a Chebyshev wave packet method, total and state-resolved reaction probabilities (J=0) were calculated for the reactions of C(1D) with various hydrogen isotopomers (H2, D2, and HD, nu i=0, j i=0) on a recent ab initio potential energy surface. For all the isotopic variants, it was found that the initial state specified reaction probabilities have no energy threshold and are strongly oscillatory, indicative of the involvement of long-lived resonances in this barrierless reaction. The J=0 product vibrational and rotational distributions for all three isotopic reactions, and the CH/CD branching ratio for the C+HD reaction, show strong dependence on the collision energy, further underscoring the important role played by the resonances. The generally decaying vibrational distributions and highly excited rotational distributions, which corroborate an insertion mechanism, and the dominance of the CD+H channel in the C+HD reaction are consistent with existing experimental observations. Initial state specified integral cross sections and rate constants were estimated using a capture model. The estimated rate constants were found to be close and in the order kHD>kH2>kD2. Finally, a method to calculate branching ratio in the C+HD reaction is proposed.     

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     

Wave packet and quasiclassical trajectory calculations for the N(D) + H2 reaction and its isotopic variants

The N(²D) + H₂(v = 0, j = 0) reaction and its HD and D₂ isotopic variants have been studied by means of quantum mechanical real wave packet and wave packet with split operator and quasiclassical trajectory methodologies on the potential energy surface of Ho et al. [J. Chem. Phys. 119 (2003) 6]. Total initial state-selected and final state-resolved reaction probabilities and product rotational distributions have been calculated for total angular momentum J = 0 in a broad range of collision energies. The real wave packet results are in very good agreement with the corresponding split operator wave packet calculations. A reasonable overall good agreement has been found between the wave packet and quasiclassical trajectory results. Integral cross-sections and thermal rate constants have been calculated from the wave packet reaction probabilities by means of standard J-shifting, refined J-shifting and uniform J-shifting methods in combination with the centrifugal sudden approximation for J > 0. Comparisons with available exact wave packet, quasiclassical trajectory and experimental results are made and discussed.     

Full-dimensional quantum dynamics study of the H + H2O and H + HOD exchange reactions

The initial state-selected time-dependent wave packet approach is employed to study the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H, HOH' + D exchange reactions with both OH bonds in the H(2)O reactant and OH(D) bond in the HOD reactant treated as reactive bonds. The total reaction probabilities for different partial waves, as well as the integral cross sections, which are the exact CC (coupled-channel) results, are first obtained in this study for the H(2)O(HOD) reactant initially in the ground rovibrational state. Because of the shallow C(3v) minimum along the reaction path, the reaction probabilities for the three reactions present several resonance peaks, with one dominant resonance peak just above the threshold. The cross sections for the H' + HOD → HOH' + D reaction are substantially smaller than those for the H' + H(2)O → H'OH + H and H' + HOD → H'OD + H reactions, indicating that the H'/H exchange reactions are much more favored. In the CC calculations, the resonance peaks in the reaction probabilities diminish quickly with the increase in total angular momenta J, resulting in the existence of a clear step-like feature just above the threshold in the cross sections for the title reactions, which manifests the signature of shape resonances in these reactions. In the CS calculations, the resonance peaks on reaction probabilities persist in many partial waves, and thus the resonance structures can no longer survive the partial-wave summation and are washed out completely in the CS cross sections for the title reactions.