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     

Barrier Height Effect on Cl+H2(D2) Reaction

Three-dimensional time-dependent quantum wave packet calculation was performed to study the reaction dynamics of Cl+H2(D2) on two potential energy surfaces (CW PESs). The first CW PES is with spin-orbit correction; the second is without spin-orbit correction. The integral cross-section and reaction probability as a function of collision energy are calculated in the collision energy range of 0.1 eV to 1.4 eV. For reaction of Cl with D2, the reaction section with spin-orbit correction has a shift toward the high energy because the barrier height increases. As for the reaction of Cl with H2 at low collision energy, it is more reactive on the PES with spin-orbit correction than on the low barrier height PES without spin-orbit correction, due to the tunnel effect for the reaction of the Cl with H2. When the collision energy is higher than 0.7 eV, the reactivity on the low barrier height PES is larger than that on the high barrier height PES. It is believed that the barrier height plays a very important role in the reactivity of Cl with (H2, D2). For the Cl+H2 reaction the barrier width is also very important because of the tunneling effect.     

Energy resonance and inverse population of the product vibrational states for the reaction F + H2(v = 0) → HF(v′) + H,LCAC‐SW theoretical quantum scattering study

LCAC‐SW (linear combination of arrangement channel‐scattering wavefunction) method was used to calculate collinear state‐to‐state reaction probabilities for the reaction F + H₂(v = 0) → HF(v′) + H on the 6SEC potential energy surface. The results show that reaction probabilities P₀₂ and P₀₃ [i. e., v′ = 2,3 for reaction F + H₂ (v = 0) + HF(v′) + H] are primary, the population of product vibrational states is inverse and the reaction probabilities are oscillatory with collision energies, i.e., there is energy resonance in this reaction, which agrees with a new experiment.     

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.     

三维含时量子散射研究H + CIH体系

用三维含时量子散射理论模拟了H+CIH体系在BW2,mBW2,G3势能面上的动力学行为,其计算结果表明,振动量子态对反应几率影响很大;势能面的地形对转动量子态如何影响反应几率起重要作用;反应几率表现出“黄金规则”,此外,BW2,mBW2势能面上的反应几率几乎相同,而G3势能面上的反应几率较前者低,大概由于G3的势垒高的缘故。     

Quantum wave packet and quasiclassical trajectory studies of the reaction H(2S) + CH(X2Π; v = 0, j = 1) → C(1D) + H2(X1 ): Coriolis coupling effects and stereodynamics

The quantum mechanics (QM) and quasiclassical trajectory (QCT) calculations have been carried out for the title reaction with the ground minimal allowed rotational state of CH (j = 1) on the 1 ¹A′ potential energy surface. For the reaction probability at total angular momentum J = 0, a similar trend of the QM and QCT calculations is observed, and the QM results are larger than the latter almost in the whole considered energy range (0.1–1.5 eV). The QCT integral cross sections are larger than the QM results with centrifugal sudden approximation, while smaller than those from QM method including Coriolis coupling for collision energies bigger than 0.25 eV. The quantum wave‐packet computations show that the Coriolis coupling effects get more and more pronounced with increasing of J. In addition to the scalar properties, the stereodynamical properties, such as the average rotational alignment factor <P₂( j′?k )>, the angular distributions P(θ_r), P(?_r), P(θ_r,?_r), and the polarization‐dependent generalized differential cross sections have been explored in detail by QCT approach. © 2013 Wiley Periodicals, Inc.     

Quantum mechanics and quasiclassical study of the H/D+FO → OH/OD+F,HF/DF+O reactions: Chemical stereodynamics

The time‐dependent quantum wave packet and the quasi‐classical trajectory (QCT) calculations for the title reactions are carried out using three recent‐developed accurate potential energy surfaces of the 1¹A′, 1³A′, and 1³A″ states. The two commonly used polarization‐dependent differential cross sections, dσ₀₀/dω_t, dσ₂₀/dω_t, with ω_t being the polar coordinates of the product velocity ω′, and the three angular distributions, P(θ_r), P(Φ_r), and P(θ_r,Φ_r), with θ_r, Φ_r being the polar angles of the product angular momentum, are generated in the center‐of‐mass frame using the QCT method to gain insight into the alignment and the orientation of the product molecules. Influences of the potential energy surface, the collision energy, and the isotope mass on the stereodynamics are shown and discussed. Validity of the QCT calculation has been examined and proved in the comparison with the quantum wave packet calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Exploring the reaction dynamics of O(3P)+ (X2 )→OH+(X3Σ−)+ H(2S) reaction with time‐dependent wave packet method

The O(³P)+ reaction has been investigated by employing time‐dependent quantum wave packet with split operator method on potential energy surface of the doublet ground‐state H₂O⁺(1²A″). The reaction probabilities and integral cross sections are calculated using centrifugal sudden approximation, which basically agree with the quasi‐classical results of Paniagua et al. [Phys. Chem. Chem. Phys. 2014, 16, 23594]. Moreover, the effect of vibrational and rotational excitation of reactant is investigated. The results show that the vibrational and rotational excitation effects on the integral cross section are not obvious. The little differences between Coriolis coupling results and centrifugal sudden approximation ones show that the cheaper centrifugal sudden calculations here reported are effective for this reaction.     

Quantum reaction dynamics of C(1D) + HD → CH(CD) + D(H) on the ground state potential energy surface

We present accurate quantum dynamic calculations of the reaction C(¹D) + HD on the latest version of the potential energy surface [Zhang et al., J. Chem. Phys. 140, 234301 (2014)]. Using a Chebyshev real wave packet method with full Coriolis coupling, we obtain the initial state‐specified ( ) reaction probabilities, integral cross sections, and rate constants. The resulting probabilities display oscillatory structures due to numerous long‐lived resonances supported by the deep potential well. The calculated rate constants and CD/CH product branching ratio at room temperature are in reasonably good agreement with the experimental measurements.     

Three-dimensional quantum dynamics study of vibrational predissociation of HeI_2 van der Waals molecule for low vibrational excitation using the time-dependent wave packet method

Three-dimensional quantum mechanical calculations for vibrational predissociation of HeI2(B) van der Waals molecules are presented using the time-dependent wave packet technique within the golden rule approxima tion.The total and partial decay widths,lifetimes,rates and their dependence on initial vibrational states were obtained for HeI2 at low initial vibrational excited levels.Our calculations show that the calculated tota decay widths,lifetimes and rates agree well with those extrapolated from experimental data available The predicted total decay widths as a function of initial vibrational states exhibit highly nonlinear behavior.The very short propagation time (less.than 1 ps) required in the golden rule wave packet calculation is determined by the duration time of the final state inter-action between the fragments on the vibrationally deexcited adiabatic potential surface.The final state interaction between the fragments is shown to play an important role in determining the final rotational distri     

Time dependent wave packet treatment of 2ΠgN2− and 3Σ−NO− shape resonances using two‐dimensional surfaces for electron‐N2 and NO interactions

The ²Π_gN and ³Σ^?NO^? resonances in electron‐N₂ and NO collisions have been treated using both nuclear and electronic degrees of freedom and a two‐dimensional (2D) time dependent wave packet approach to ascertain the importance of nonlocality in electron–nuclear interaction. The results so obtained are compared with vibrational excitation cross‐sections obtained experimentally and those from other theoretical/numerical approaches using 1D local complex potential, 2D model with a combination of the exterior complex scaling method and a finite‐element implementation of the discrete‐variable representation. The results obtained provide detailed insight into the nuclear dynamics induced by electron–molecule collision and reveal that while for resonant excitation of lower vibrational modes, the nonlocal effect may not be as critical but importance of nonlocal effects may increase with increase in quanta of resonant vibrational excitation. © 2012 Wiley Periodicals, Inc.     

Quantum dynamics of dissociative adsorption of H_2 and D_2:The time-dependent wave packet method

Ｒｅｃｅｎｔｐｒｏｇｒｅｓｓｉｎｓｕｒｆａｃｅｓｐｅｃｔｒｏｓｃｏｐｙａｎｄｍｏｌｅｃｕｌａｒｂｅａｍｓｃａｔｔｅｒｉｎｇａｎｄｄｅｔｅｃｔｉｏｎｔｅｃｈｎｉｑｕｅｓｍａｋｅｓｉｔｎｏｗｐｏｓｓｉｂｌｅａｔｔｈｅｍｉｃｒｏｓｃｏｐｉｃｌｅｖｅｌｔｏａｃｃｕｒａｔｅｌｙｍｅａｓｕｒｅｓｕｃｈｔｈｉｎｇｓａｓｄｉｓｓｏｃｉａｔｉｏｎｒａｔｅｓ,ａｄｓｏｒｂａｔｅｂｉｎｄｉｎｇａｎｄｇｅｏｍｅｔｒｙ,ａｎｄｍｏｂｉｌｉｔｉｅｓｏｆａｄｓｏｒｂｅｄｓｐｅｃｉｅｓｏｎｓｕｒｆａｃｅｓ.Ｔｈｉｓｋｉｎｄｏｆｅ…     

Highly excited vibrational levels of triatomic molecule N2O: U(4) algebraic model

In the U(4) algebraic framework, the triatomic molecules are of U₁(4) ? U₂(4) dynamical symmetry. A molecular Hamiltonian is constructed including the third‐order conbination of the invariant operators. Within this framework, the highly vibrational energy levels of the linear triatomic nitrous oxide molecule, including both bending and stretching vibrations, are studied. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Three-dimensional quantum dynamics study of vibrational predissociation of HeI2 van der Waals molecule for low vibrational excitation using the time-dependent wave packet method

Three-dimensional quantum mechanical calculations for vibrational predissociation of He1₂(B) van der Waals molecules are presented using the time-dependent wave packet technique within the golden rule approximation. The total and partial decay widths, lifetimes, rates and their dependence on initial vibrational states were obtained for HeI₂ at low initial vibrational excited levels. Our calculations show that the calculated total decay widths, lifetimes and rates agree well with those extrapolated from experimental data available. The predicted total decay widths as a function of initial vibrational states exhibit highly nonlinear behavior. The very short propagation time (less than 1 ps) required in the golden rule wave packet calculation is determined by the duration time of the final state interaction between the fragments on the vibrationally deexcited adiabatic potential surface. The final state interaction between the fragments is shown to play an important role in determining the final rotational distribution. This interpretation clearly explains the dynamical effect that the final rotational distribution shifts to the lower rotational energy levels as the initial vibrational quantum numberu increases.     

The investigation of nonadiabatic effects for the N + ND → N2 + D reaction

Nonadiabatic quantum dynamical calculations have been carried out on the two coupled potential energy surfaces (1²A′ and 2²A′) (Mota et al., J Theor Comput Chem 2009, 8, 849) for the title reaction. Initial state‐resolved reaction probabilities and cross sections for ground and excited states for collision energies of 0.005–1.0 eV are determined, respectively. Nonadiabatic transition is enhanced about four times by isotopic substitution of N + NH by N + ND reaction. It turns out that the nonadiabatic effects exert no significant contribution in the N + ND → N₂ + D reaction. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Quantum dynamics of dissociative adsorption of H2 and D2: The time-dependent wave packet method

A time-dependent quantum wave packet method was used to study the dynamics of dissociative adsorption of H₂ and D₂ on a flat and static surface. The molecule-surface interaction is described using a modified London-Eyring-Polanyi-Sato (LEPS) type potential for the H₂/Ni(100) system. The three-dimensional (3-D) dissociation probabilities were calculated for different initial rovibrational states as a function of initial incident energies. Our results show that the dissociation of the diatomic rotational states whose quantum numbers satisfyj+m = odd is forbidden at low energies for the homonuclear H_z and D₂ due to the selection rule. The effect of the rotational orientation of diatoms on adsorption predicts that the in-plane rotation (m = j) is more favorable for dissociation than the out-of-plane rotation (m = 0). Enhanced dissociation for vibrationally excited molecules and the significant enhancement of the dissociation probability of H2 when compared to D₂ were explained reasonably in terms of quantum mechanical zero-point energies, the tunneling effect and the reflection from an activation barrier. Project supported by the National Natural Science Foundation of China (Grant No. 19694033) and partially by the Science Foundation for Overseas Chinese Scholars and Students, administered by the State Education Commission of China (Grant No. 1992), by the State Key Laboratory of Theoretical and Computational Chemistry of Jilin University at Changchun (Grant No. 98011, and by the Natural Science Foundation of Shandong Province (Grant No. Y96B03022)     

Global X2A′ potential energy surface of Li2H and quantum dynamics of H + Li2 (X1Σg+) → Li + LiH (X1Σ+) reaction

A global potential energy surface (PES) for the electronic ground state of Li₂H system is constructed over a large configuration space. About 30 000 ab initio energy points have been calculated by MRCI‐F12 method with aug‐cc‐pVTZ basis set. The neural network method is applied to fit the PES and the root mean square error of the current PES is only 1.296 meV. The reaction dynamics of the title reaction has been carried out by employing time‐dependent wave packet approach with second order split operator on the new PES. The reaction probability, integral cross section and thermal rate constant are obtained from the dynamics calculation. In most of the collision energy regions, the integral cross sections are in well agreement with the results reported by Gao et al. The rate constant calculated from the new PES increases in the temperature range of present investigation.     

Dynamical properties of S(3P) + HD reaction on 13A″state and their quantum wavepacket calculation

The time‐dependent wavepacket method is used to study the reaction dynamics of S(³P) + HD (v = 0, 1, 2) on the adiabatic 1³A″ potential energy surface constructed by Han and coworkers [J. Chem. Phys. 2012, 136, 094308]. The reaction probabilities and integral cross sections as a function of collision energy are obtained and discussed. The results calculated by using the CC and the CS approximation have been compared, which suggests that for this direct abstraction reaction, the cheaper CS approximation calculation is valid enough in the quantum calculation. The investigation also shows that the reaction probabilities and integral cross sections tend to increase with collision energy. By analyzing the v‐dependent behavior of the integral cross sections, the significant effect of the vibrational excitation of HD is found. Also found in the calculation is a significant resonance feature in the reaction probabilities versus collision energy. © 2014 Wiley Periodicals, Inc.     

A quantum scattering study of collinear state-testate reaction probabilities for the F + H2(v) → HF(v') + H system

A new quantum scattering approach (linear combination of arrangement channels-scattering wavefunction, LCAC-SW) proposed by Deng and his co-workers is used to calculate collinear state-to-state reaction probabilities for the F + H₂(v) → HF(v') + H system. Several interesting problems such M threshold energy, compound states and enhance by translational energy of the reactants and the vibration excitation of products are discussed and they are compared with other theoretical investigations reported in the literature. It is shown that the LCAC-SW approach is the successful one of quantum scattering methods.