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.     

Full-dimensional time-dependent wave packet dynamics of H2 + D2 reaction

Collision induced dissociation (CID), four center reaction (4C), and single exchange reaction (SE) in H(2) (v(1) = high) + D(2) (v(2) = low) were studied by means of time-dependent wave packet approach within a full-dimensional model. Initial state-selected total reaction probabilities for the three competitive processes have been computed on two realistic global potential energy surfaces of Aguado-Suárez-Paniagua and Boothroyd-Martin-Keogh-Peterson (BMKP) with the total angular momentum J = 0. The role of both vibrationally excited and rotationally excited reagents was examined by varying the initial vibrational and rotational states. The vibrational excitation of the hot diatom gives an enhancement effect on the CID process, while the vibrational excitation of the cold diatom gives an inhibition effect. The rotational excitation of both reagents has a significant effect on the reaction process. The 4C and SE probabilities are at least one order of magnitude smaller than the CID probabilities over the energy range considered. Isotope substitution effects were also studied by substituting the collider D(2) by H(2) and HD on the BMKP potential energy surfaces. The CID process is most efficient for the H(2) + D(2) combination and least efficient for the H(2) + H(2) combination and is different for the 4C and SE processes.     

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     

Renner-Teller quantum dynamics of the N(2D) + H2-->NH + H reaction

We present the Born-Oppenheimer (BO) and Renner-Teller (RT) quantum dynamics of the reaction (14)N((2)D)+(1)H(2)(X (1)Sigma(g) (+))-->NH(X (3)Sigma(-))+H((2)S), considering the NH(2) electronic states X (2)B(1) and A (2)A(1). These states correlate to the same (2)Pi(u) linear species, are coupled by RT nonadiabatic effects, and give NH(X (3)Sigma(-))+H and NH(a (1)Delta)+H, respectively. We develop the Hamiltonian matrix elements in the R embedding of the Jacobi coordinates and in the adiabatic electronic representation, using the permutation-inversion symmetry, and taking into account the nuclear-spin statistics. Collision observables are calculated via the real wave-packet (WP) and flux methods, using the potential-energy surfaces of Santoro et al. [J. Phys. Chem. A 106, 8276 (2002)]. WP snapshots show that the reaction proceeds via an insertion mechanism, and that the RT-WP avoids the A (2)A(1) potential barrier, jumping from the excited to the ground surface and giving mainly the NH(X (3)Sigma(-)) products. X (2)B(1) BO probabilities and cross sections show large tunnel effects and are approximately four to ten times larger than the A (2)A(1) ones. This implies a BO rate-constant ratio k(X (2)B(1))k(A (2)A(1)) approximately 10(5) at 300 K, i.e., a negligible BO formation of NH(a (1)Delta). When H(2) is rotationally excited, RT couplings reduce slightly the X (2)B(1) reaction observables, but enhance strongly the A (2)A(1) reactivity. These couplings are important at all collision energies, reduce the collision threshold, and increase remarkably reaction probabilities and cross sections. The RT k(A (2)A(1)) is thus approximately 3.3 order of magnitude larger than the BO value, and degeneracy-averaged, initial-state-resolved rate constants increase by approximately 13% and by approximately 47% at 300 and 500 K, respectively. Owing to an overestimation of the X (2)B(1) potential barrier, the calculated thermal rate is too low with respect to that observed, but we obtain a good agreement by shifting down the calculated cross section.     

Nonadiabatic time-dependent wave packet study of the D+ + H2 reaction system

A theoretical investigation on the nonadiabatic processes of the D(+) + H(2) reaction system has been carried out by means of exact three-dimensional nonadiabatic time-dependent wave packet calculations with an extended split operator scheme (XSOS). The diabatic potential energy surface newly constructed by Kamisaka et al. (J. Chem. Phys. 2002, 116, 654) was employed in the calculations. This study provided quantum cross sections for three competing channels of the reactive charge transfer, the nonreactive charge transfer, and the reactive noncharge transfer, which contrasted markedly to many previous quantum theoretical reports on the (DH(2))(+) system restricted to the total angular momentum J = 0. These quantum theoretical cross sections derived from the ground rovibrational state of H(2) show wiggling structures and an increasing trend for both the reactive charge transfer and the nonreactive charge transfer but a decreasing trend for the reactive noncharge transfer throughout the investigated collision energy range 1.7-2.5 eV. The results also show that the channel of the reactive noncharge transfer with the largest cross section is the dominant one. A further investigation of the v-dependent behavior of the probabilities for the three channels revealed an interesting dominant trend for the reactive charge transfer and the nonreactive charge transfer at vibrational excitation v = 4 of H(2). In addition, the comparison between the centrifugal sudden (CS) and exact calculations showed the importance of the Coriolis coupling for the reactive system. The computed quantum cross sections are also compared with the experimental measurement results.     

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.     

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.     

Time-dependent quantum mechanical wave packet study of the He+H(2) (+)(v,j)-->HeH(+)+H reaction

A detailed three-dimensional time-dependent quantum dynamical study of the He+H(2) (+)(v=0-3,j=0)-->HeH(+)+H reaction is reported for different vibrational v states of H(2) (+) in its ground rotational (j=0) state over a range of translational E(trans) energies on an accurate ab initio potential energy surface published by Palmieri et al. Plots of reaction probability as a function of total energy E reveal a large number of oscillations indicating the presence of a number of reactive scattering resonances. When averaged over total angular momentum J, some of the oscillations survive, indicating that they may be amenable to experimental observation. A comparison of our present results with our earlier results on the McLaughlin-Thompson-Joseph-Sathyamurthy surface and the experimental results from different research groups reveal a good deal of agreement as well as some discrepancies between theory and experiment at the level of state-selected gas phase dynamics.     

A transition state wave packet study of the H+CH4 reaction

Transition state wave packet calculations have been carried out to obtain the thermal rate constants for the H+CH(4) reaction on the Jordan-Gilbert potential energy surface. The eight-dimensional model for the X+YCZ(3) type of reaction was employed by restricting the nonreacting CZ(3) group under a C(3V) symmetry. We calculated the cumulative reaction probability for the total angular momentum J=0, from which the thermal rate constants were obtained for the temperature range between 250 and 500 K by employing the J-K shifting approximation. It is found that the eight-dimensional rate constants agree very well with the seven dimensional ones, in which the CH bond length in the nonreacting CH(3) group is fixed, suggesting that the additional mode for the symmetry stretching in CH(3) group does not have any important effect on the reaction within the temperature range considered here. The present transition state wave packet results agree well with rate constants obtained from the previous seven dimensional initial state selected wave packet study, indicating the consistence of the treatments involved in these two studies. On the other hand, it is found that the energy threshold for the cumulative reaction probability for J=0 from the present study is higher than that from the full dimensional multiconfiguration time-dependent Hartree study by about 0.06 eV, resulting in severe discrepancy between the present rate constants and the multiconfiguration time-dependent Hartree results.     

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.     

Exact quantum dynamics of N(2D) + H2 --> NH + H reaction: cross-sections, rate constants, and dependence on reactant rotation

Using an exact Chebyshev wave packet method, initial state-specified (upsilon(i)=0, j(i)=0,2) integral cross-sections and rate constants are obtained for the title reaction on the latest ab initio potential energy surface. Reaction probabilities up to J=29 are dependent on the reactant rotation and show mild oscillations superimposed on a broad background. Due to a barrier in the entrance channel, the cross sections increase with energy with clear thresholds and the rate constants vary with temperature in the Arrhenius form. The calculated canonical rate constant is in good agreement with the experimental measurements. Our results also indicate that the quasiclassical trajectory method underestimates the rate due to the neglect of tunneling, while the quantum statistical approach overestimates because of the short lifetime of the reaction intermediate.     

A time-dependent wave packet quantum scattering study of the reaction HD+ (v = 0 - 3;j0 = 1) + He --> HeH+(HeD+) + D(H)

Time-dependent wave packet quantum scattering (TWQS) calculations are presented for HD(+) (v = 0 - 3;j(0)=1) + He collisions in the center-of-mass collision energy (E(T)) range of 0.0-2.0 eV. The present TWQS approach accounts for Coriolis coupling and uses the ab initio potential energy surface of Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. For a fixed total angular momentum J, the energy dependence of reaction probabilities exhibits quantum resonance structure. The resonances are more pronounced for low J values and for the HeH(+) + D channel than for the HeD(+) + H channel and are particularly prominent near threshold. The quantum effects are no longer discernable in the integral cross sections, which compare closely to quasiclassical trajectory calculations conducted on the same potential energy surface. The integral cross sections also compare well to recent state-selected experimental values over the same reactant and translational energy range. Classical impulsive dynamics and steric arguments can account for the significant isotope effect in favor of the deuteron transfer channel observed for HD(+)(v<3) and low translational energies. At higher reactant energies, angular momentum constraints favor the proton-transfer channel, and isotopic differences in the integral cross sections are no longer significant. The integral cross sections as well as the J dependence of partial cross sections exhibit a significant alignment effect in favor of collisions with the HD(+) rotational angular momentum vector perpendicular to the Jacobi R coordinate. This effect is most pronounced for the proton-transfer channel at low vibrational and translational energies.     

The quantum dynamics of the reactions N+H2(HD,D2) and their vibrational excitation effect

The N(⁴S)+H₂ reaction and its isotopic variants have been investigated by means of time‐dependent quantum wave packet with split operator method on the ground state potential energy surface (Zhai and Han, J. Chem. Phys. 2011, 135, 104314). The reaction probabilities, integral cross sections, branching ratio of the integral cross sections, and effect of vibrational excitation of H₂, HD, and D₂ diatomic molecules are presented and discussed. The results reveal that the intramolecular isotopic effect is greater than the intermolecular one, and that the vibrational excitation of the diatomic molecules can promote the progress of this reaction. In addition, a limited number of rigorous Coriolis coupling calculations of the integral cross sections of the N(⁴S)+H₂ reaction have been carried out. Also shown is that since the Coriolis coupling plays a small part in this accurate quantum calculation, the cheaper centrifugal sudden calculations here reported are effective for this reactive system. © 2014 Wiley Periodicals, Inc.     

A time-dependent wave packet study of the H4 four-center reaction

A quantum model based on the time-dependent initial state selected wave packet approach was developed to study the four-center (4C) reaction, A₂ + B₂ → 2AB, and the competing collision induced dissociation (CID), A₂ + B₂ → A + B₂ + A, as applied to the H₂(v₁) + H₂(v₂) system important in combustion. A reduced three-dimensional model of the reaction with the atoms constrained to an isosceles trapezium and a realistic global potential energy surface of Aguado et al. [J. Chem. Phys. 101 (1994) 2742], following Hernández and Clary [J. Chem. Phys. 104 (1996) 8413], was used. A method to analyse the reaction flux for 4C and CID reaction probabilities is presented. The initial A₂ vibrational excitation is not only more efficient than translational energy in facilitating the 4C and CID processes, it also reduces the threshold energy. Both the 4C and CID processes exhibit similar threshold energy behavior. For low vibrational excitation in the A₂ diatom, the 4C process is dominant; as the A₂ diatom becomes highly excited the CID process becomes more important at low collision energies with B₂, but as the collision energy increases the 4C process is favored again.     

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.     

A quantum wave packet dynamical study of the electronic and spin-orbit coupling effects on the resonances in Cl(2P) + H2 scattering

Dynamical resonances in Cl(2P) + H2 scattering are investigated with the aid of a time-dependent wave packet approach using the Capecchi-Werner coupled ab initio potential energy surfaces [Phys. Chem. Chem. Phys. 2004, 6, 4975]. The resonances arising from the prereactive van der Waals well (approximately 0.5 kcal/mol) and the transition-state (TS) region of the 2Sigma(1/2) ground spin-orbit (SO) state of the Cl(2P) + H2 system are calculated and assigned by computing their eigenfunctions and lifetimes. The excitation of even quanta along the bending coordinate of the resonances is observed. The resonances exhibit an extended van der Waals progression, which can be attributed to the dissociative states of ClH2. Excitation of H2 vibration is also identified in the high-energy resonances. The effect of the excited 2P(1/2) SO state of Cl on these resonances is examined by considering the electronic and SO coupling in the dynamical simulations. While the electronic coupling has only a minor impact on the resonance structures, the SO coupling has significant effect on them. The nonadiabatic effect due to the SO coupling is stronger, and as a result, the spectrum becomes broad and diffuse particularly at high energies. We also report the photodetachment spectrum of ClD2- and compare the theoretical findings with the available experimental results.     

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.