• Bulletin of the Korean Chemical Society
• /
• v.32 no.12
• /
• pp.4210-4214
• /
• 2011

Quasiclassical trajectory (QCT) calculations of Ne + ${H_2}^+$ reaction have been carried out on the adiabatic potential energy surface of the ground state $1^2$ A'. The reaction probability of the title reaction for J = 0 has been calculated, and the QCT result is consistent with the previous quantum mechanical wave packet result. Quasiclassical trajectory calculations of the four polarization-dependent differential cross sections have been carried out in the center of mass (CM) frame. The P(${\theta}_r$), P(${\phi}_r$) and P(${\theta}_r$, ${\phi}_r$) distributions, the k-k'-j' correlation and the angular distribution of product rotational vectors are presented in the form of polar plots. Due to the well in $1^2$ A' PES, the reagent vibrational excitation has greater influence on the polarization of the product rotational angular momentum vectors j' than the collision energy.

• Bulletin of the Korean Chemical Society
• /
• v.27 no.12
• /
• pp.1976-1980
• /
• 2006

We have calculated the reaction probability and the reaction cross-section of the $N(^4S)+O_2(X^3\sum_{g}^{-})\;\rightarrow\;NO(X^2\Pi)+O(^3P)$ reaction by the quasiclassical trajectory method with the 6th-order explicit symplectic algorithm, based on a new ground potential energy surface. The advantage of the 6th-order explicit symplectic algorithm, conserving both the total energy and the total angular momentum of the reaction system during the numerical integration of canonical equations, has firstly analyzed in this work, which make the calculation of the reaction probability more reliable. The variation of the reaction probability with the impact parameter and the influence of the relative translational energy on the reaction cross-section of the reaction have been discussed in detail. And the fact is found by the comparison that the reaction probability and the reaction cross-section of the reaction estimated in this work are more reasonable than the theoretical ones determined by Gilibert et al.

• Bulletin of the Korean Chemical Society
• /
• v.31 no.10
• /
• pp.2841-2848
• /
• 2010

Quasiclassical trajectory (QCT) method has been used to investigate stereodynamics information of the reaction $O(^1D)+H_2{\rightarrow}\;OH$+H on the DK (Dobbyn and Knowles) potential energy surface (PES) at a collision energy of 23.06 kcal/mol, with the initial quantum state of reactant $H_2$ being set for v = 0 (vibration quantum number) and j = 0-5 (rotation quantum number). The PDDCSs (polarization dependent differential cross sections) and the distributions of P($\theta_r$), P($\phi_r$), P($\theta_r$, $\phi_r$) have been presented in this work. The results demonstrate that the products are both forward and backward scattered. As j increases, the backward scattering becomes weaker while the forward scattering becomes slightly stronger. The distribution of P($\theta_r$) indicates that the product rotational angular momentum j' tends to align along the direction perpendicular to the reagent relative velocity vector k, but this kind of product alignment is found to be rather insensitive to j. Furthermore, the distribution of P($\phi_r$) indicates that the rotational angular momentum vector of the OH product is preferentially oriented along the positive direction of y-axis, and such product orientation becomes stronger with increasing j.

• Bulletin of the Korean Chemical Society
• /
• v.19 no.6
• /
• pp.634-641
• /
• 1998

Quasiclassical trajectory calculations were carried out for the reactions of $H+H_2$ (V=O, J=O) and its isotope variants on the Siegbahn-Liu-Truhlar-Horowitz potential energy surface for the relative energies E between 6 and 150 kcal/mol. The goal of the work was to understand the inter- and intramolecular isotope effects. We examine the relative motion of reactants during the collision using the method of analysis that monitors the intermolecular properties (internuclear distances, geometry of reactants, and final product). As in other works, we find that the heavier the incoming atom is, the greater the reaction cross section is at the same collision energy. Using the method of analysis we prove that the intermolecular isotope effect is contributed mainly by differences in reorientation due to the different reduced masses. We show that above E=30 kcal/mol recrossing also contributes to the intermolecular isotope effect. For the intramolecular isotope effect in the reactions of H+HD and T+HD, we reach the same conclusions as in the systems of $O(^3P)+HD$, F+HD, and Cl+HD. That is, the intramolecular isotope effect below E=150 kcal/mol is contributed by reorientation, recrossing, and knockout type reactions.

• Journal of the Korean Chemical Society
• /
• v.47 no.3
• /
• pp.191-198
• /
• 2003

The pairwise energy model (PEM) assumes that the cross section for the reaction cross section for the reaction A+BC$\{leftrightarrow}$B+C, where B and C are isotopes of hydrogen, depends on only the pairwise relative energy Es between A and B. Until now, the PEM has been used to interpret theoretically the isotope effect for the reactions such as $O(^3P)+HD,\;Ar^++(H_2,\;D_2,and\;HD)$. In this paper we carry out extensive quasiclassical trajectory calculations for the three possible reactions $Cl+H_2$ and HD and show that the PEM works very well at high energy. In particular we are able to accurately predict the intramolecular isotope effect at high energy for the reaction of Cl+HD using only the cross section data for $Cl+H_2$. To understand that the PEM works so well at high energy, the internal energy distributions for the products are examined. The distributions for three reactions are different at a fixed relative collision energy E but are approximately same at a fixed pairwise energy Es. This suggests that the PEM works very well at high energy. We believe the conclusions reached here will apply to other A+BC systems.

• Bulletin of the Korean Chemical Society
• /
• v.23 no.2
• /
• pp.229-240
• /
• 2002

Quantum mechanical calculation is performed for the $O(^1D)$ + HCl ${\rightarrow}$OH + Cl reaction using Reactive Infinite Order Sudden Approximation. Shifting approximation is also employed for the l ${\neq}$ 0 partial wave contributions. Various dynamical quantities are calculated and compared with available experimental results and quasiclassical trajectory results. Vibrational distributions agree well with experimental results i.e. product states mostly populated at $v_f$ = 3, 4. Our results also show small peak at $v_f$ = 0, which indicates bimodal vibrational distribution. The results show two significant broad peaks in ${\gamma}_i$ dependence of the cross section, one is at ${\gamma}_i$ = $15^{\circ}-35^{\circ}$ and the another is at ${\gamma}_i$= $55^{\circ}-75^{\circ}$ which can be explained as steric effects. At smaller gi, the distribution is peaked only at higher state ($v_f$ = 3, 4) while at the larger gi, both lower state ($v_f$ = 0) and higher state ($v_f$ = 3, 4) are significantly populated. Such two competing contributions (smaller and larger ${\gamma}_i$) result in the bimodal distribution. From these points we suggest two mechanisms underlying in current reaction system: one is that reaction occurs in a direct way, while the another is that reaction occurs in a indirect way.