Global dynamics and transition state theories: Comparative study of reaction rate constants for gas-phase chemical reactions

Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package

A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.

A new global analytical ab initio potential energy surface for the dynamics of the C+(2P) + SH(X2Π) reaction

The global potential energy surface (PES) of HCS⁺(X¹Σ⁺) is constructed using many-body expansion (MBE) methodology. The obtained analytical function is found by fitting the 7907 ab initio energy points computed at the Davidson-corrected multi-reference configuration interaction level with the aug-cc-pV(5+d)Z basis set. The final root mean square error is 0.0419 eV, and the maximum deviation is 0.2039 eV, showing that the analytical formula agrees well with the energy points. The topological features are calculated and discussed based upon the analytical PES of HCS⁺(X¹Σ⁺). The reaction probability, integral cross sections and other details of the C⁺(²P) + SH(X²Π) → H(²S) + CS⁺(X²Σ⁺) reaction are investigated using the quasi-classical trajectory and time-dependent quantum wave packet methods.

Indirect photodegradation of fludioxonil by hydroxyl radical and singlet oxygen in aquatic environment: Mechanism, photoproducts formation and eco-toxicity assessment

Fludioxonil has been proven valuable as a broad-spectrum fungicide. However, there are concerns about its risk posed to non-target organisms in aquatic environments. In this paper, the mechanism, photoproducts transformation and eco-toxicity of fludioxonil during •OH/¹O₂-initiated process were systematically studied using quantum chemistry and computational toxicology. The results indicate that the two favorable pathways of •OH/¹O₂-initiated reactions are both occurred in pyrrole ring. It can conclude that the rate constants of •OH and ¹O₂ are 1.23 × 10¹⁰ and 3.69 × 10⁷ M^-1 s^-1 at 298K, respectively, which results in half-lives of <2 days in surface waters under sunlit near-surface conditions. Based on toxicity assessments, these photoproducts showed a decreased aquatic toxicity but the majority products are still toxic. This study gives more insight into the chemical transformation mechanism of fludioxonil in aquatic environments.

Investigation of the Abstraction and Dissociation Mechanism in the Nitrogen Trifluoride Channels: Combined Post-Hartree-Fock and Transition State Theory Approaches

The present paper concludes our series of kinetics studies on the reactions involved in the complex mechanism of nitrogen trifluoride decomposition. Two other related reactions that, along with this mechanism, take part in an efficient boron nitride growth process are also investigated. We report results concerning two abstraction reactions, namely NF2 + N ⇄ 2NF and NF3 + NF ⇄ 2NF2, and two dissociations, N2F4 ⇄ 2NF2 and N2F3 ⇄ NF2 + NF. State-of-the-art electronic structure calculations at the CCSD(T)/cc-pVTZ level of theory were considered to determine geometries and frequencies of reactants, products, and transition states. Extrapolation of the energies to the complete basis set limit was used to obtain energies of all the species. We applied transition state theory to compute thermal rate constants including Wigner, Eckart, Bell, and deformed theory corrections in order to take tunneling effects into account. The obtained results are in good agreement with the experimental data available in the literature and are expected to provide a better phenomenological understanding of the NF3 decomposition role in the boron nitride growth for a wide range of temperature values.

Classical blended quantum formulation of the parallel TDDVR method to the dynamics on furan

A major portion (∼98%) of the quantum-classical molecular dynamics approach, namely, Time-Dependent Discrete Variable Representation (TDDVR) method, is recently parallelized using shared-memory parallelization scheme with the aim of performing dynamics on relatively large molecular systems. Therefore, we have chosen the furan as a model system for dynamical simulation. Four lowest singlet excited electronic states 1A2(3s), 1B2(V), 1A1(V*), and 1B1(3p) of furan are vibronically coupled through several conical intersections in the vicinity of Frank-Condon region. The major focus of the present paper is to explore the efficiency of our newly implemented parallelized TDDVR algorithm to perform dynamics on large dimensional quantum systems in vibronically coupled electronic manifold. The present version of parallelized TDDVR algorithm show closely linear speed up with increasing number of computing processors. As a significant speed up is achieved by cycling in the correct way over arrays, all dynamical simulations are performed within a reasonable wall clock time. The photoelectron spectra calculated by the TDDVR method show peak by peak correspondence with the experimental spectra. The TDDVR calculated quantum-classical dynamical outcomes, viz., population and photoelectron spectra, etc. show good agreement with the findings of the well-established quantum dynamical method, i.e. Multi Configuration Time-Dependent Hartree (MCTDH) approach.

Quasi-classical trajectory study of the reaction dynamics of calcium ground state and metastable atoms with CH2Cl2

The dynamics properties of the Ca(1S0,3P) + CH2Cl2 reaction system have been calculated by means of the quasi-classical trajectory method based on the extended London–Eyring–Polanyi–Sato potential energy surface. By the calculations, the vibrational distribution, reaction cross section, rotational alignment, and reaction rate constant are obtained. The peak location of vibrational quantum numbers is at ν = 0 when the collision energy is 2.302 kcal/mol, whether the calcium atom is at the ground state or metastable state. The product vibrational distribution agrees well with the experiment value in Han, K. L.; He, G. Z.; Lou, N. Q. Chem. Phys. Lett. 1991, 178, 528. The cross section thoroughly decreases with the increase of the collision energy. The rotational alignment of the product greatly deviates from –0.5. The reaction rate constant increases with rising temperature.

Non-Born–Oppenheimer stereodynamics study for the D+ + H2 (v, j = 0) reaction using coherent switching with decay of mixing method

A non-Born–Oppenheimer (BO) investigation using coherent switching with decay of mixing method is carried out for both reactive non-charge transfer and charge transfer channels of reaction D+ + H2 (v, j = 0). For the non-charge transfer channel, it is found that the non-BO effect has a small influence on the differential cross section, but it has a large impact on the rotational polarization of the HD product, which is reflected by the joint distribution P(θr, ϕr). And it is also found that the polarization of the charge transfer products shows a quite different form. The rotational polarizations are also investigated at various collision energies and initial vibrational states for both reactive channels.

CO desorption from a catalytic surface: elucidation of the role of steps by velocity-selected residence time measurements

Directly measuring the rate of a surface chemical reaction remains a challenging problem. For example, even after more than 30 years of study, there is still no agreement on the kinetic parameters for one of the simplest surface reactions: desorption of CO from Pt(111). We present a new experimental technique for determining rates of surface reactions, the velocity-selected residence time method, and demonstrate it for thermal desorption of CO from Pt(111). We use UV−UV double resonance spectroscopy to record surface residence times at selected final velocities of the desorbing CO subsequent to dosing with a pulsed molecular beam. Velocity selection differentiates trapping-desorption from direct scattering and removes influences on the temporal profile arising from the velocity distribution of the desorbing CO. The kinetic data thus obtained are of such high quality that bi-exponential desorption kinetics of CO from Pt(111) can be clearly seen. We assign the faster of the two rate processes to desorption from (111) terraces, and the slower rate process to sequential diffusion from steps to terraces followed by desorption. The influence of steps, whose density may vary from crystal to crystal, accounts for the diversity of previously reported (single exponential) kinetics results. Using transition-state theory, we derive the binding energy of CO to Pt(111) terraces, D(0)(terr) (Pt−CO) = 34 ± 1 kcal/mol (1.47 ± 0.04 eV) for the low coverage limit (≤0.03 ML) where adsorbate−adsorbate interactions are negligible. This provides a useful benchmark for electronic structure theory of adsorbates on metal surfaces.

Study on the variation rules of the joint effects for multicomponent mixtures containing cyanogenic toxicants and aldehydes based on the transition state theory

Although the study of the variation rules of the joint effects for multicomponent mixtures has gained increasing attention, it still remains unclear how the variation occurs and what the relationships between the joint effects of multicomponent mixtures and their corresponding binary mixtures are. To explain how the variation occurs, this study first proposes a hypothesis on the variation rules of the joint effects using the well-known transition state theory. The hypothesis concluded that the joint effect of multicomponent mixtures is among the joint effects of the corresponding binary mixtures. This hypothesis was named the fishing hypothesis because there is a similarity between the action process of the joint effects and the fishing process. Next, the hypothesis was validated by use of the experimental data by evaluating the joint effects of binary, ternary and quaternary mixtures containing cyanogenic toxicants and aldehydes on Photobacterium phosphoreum. The application of the fishing hypothesis can explain the rule as to how the joint effects of a multicomponent-mixture vary with its number of components and their ratios. This study provides a good method to predict the joint effects of multicomponent mixtures using the joint effects of their corresponding binary mixtures. An improvement in the fishing hypothesis will be needed in our future studies due to the approximate assumptions used in the deduction of the hypothesis.

Stereodynamic study of the reaction H(²S) + ClO(²Π) → HO(²Π) + Cl(²P) via quasi-classical trajectory calculations

Quasi-classical trajectory calculations were performed to investigate the stereodynamics of the reaction H(²S) + ClO(²Π) → HO(²Π) + Cl(²P) using the ground-state potential energy surface 1¹A'. The alignment and orientation of the product molecules as well as the four polarization-dependent differential cross-sections (PDDCSs) for this reaction across a wide range of collision energies (0.05-1.0 eV) and for the rovibrational state ClO(v = 0 and j = 0) were obtained and are reported here. It was found that the OH product rotational polarization is not very sensitive to the collision energy selected. We discuss this phenomenon in detail. The calculated results indicate that, for this system, the two deep wells in the potential energy surface are very likely to be powerful influences on the degree of product rotational polarization. In addition, the microscopic reaction mechanism that dictates the product angular momentum orientation was investigated. The forward peak in the PDDCS₀₀ at θ = 0° in our study showed a strong dependence on the initial collision energy.

Time-dependent wave packet quantum scattering and quasi-classical trajectory calculations of the H + FCl(v=0,j=0) → HF + Cl/HCl + F reaction

The dynamics of the title reaction are investigated using both time-dependent wave packet quantum scattering and quasi-classical trajectory (QCT) methods on adiabatic ground 1(2)A' potential energy surface (PES). Compared with the quantum results of reaction probabilities of H + FCl(J=0) → HF + Cl/HCl + F, the QCT method is proven feasible and further employed to produce integral cross sections and rate constants. Significant resonance structures are observed in the reaction probabilities using the quantum method; however, there are some undulations in the calculated QCT integral cross sections for both product channels. A comparison between the quantum mechanical coupled-channel (CC) calculation and centrifugal sudden approximation calculation reveals the very important role of Coriolis coupling effects in the quantum calculation. Comparisons between the calculated thermal rate constants for both reactions and the previous theoretical and experimental results have been done. HCl product formation is favored over the HF product in the reactive system. Finally, the HF products are found to be mainly forward scattering, and the HCl products are mainly backward scattering.

QUANTUM CHEMICAL STUDY ON THERMAL DECOMPOSITION MECHANISM OF CALCIUM CARBONATE

A detailed quantum chemical calculation is performed at the MP2(full)/6-311G* level to explore the mechanism of calcium carbonate thermal decomposition. Four microscopic pathways are identified. The rate constants of rate-determining steps in four pathways are calculated over a temperature range 298–1200 K. The calculating results show that only path A (R( CaCO 3) → IM 1 → P( CaO + CO 2)) and path B (R( CaCO 3) → IM B1 → IM B2 → P( CaO + CO 2) have contributions to the CaCO 3 thermal decomposition, and path A may be more favorable than path B. The present theoretical studies may provide useful information in understanding reaction mechanism of metal carbonates at the molecular level.

Energy-dependent stereodynamics for the H(2S)+NH(X3)→H2(X1Σg+)+N(4S) reaction on the improved ZH potential energy surface

In this work, quasi-classical trajectory (QCT) calculations have been first carried out for the title reaction on a new global potential energy surface for the lowest quartet electronic state, 4A″. The average rotational alignment factor [P2(j'·k)] as a function of collision energy and the two commonly used polarization dependent generalized differential cross sections PDDCS00, PDDCS20, have been calculated in the center-of-mass (CM) frame, separately. Three angular distributions, P(θr), P(φr), and P(θr, φr) are also calculated to gain insight into the alignment and the orientation of the product molecules. Calculations show that the average rotational alignment factor on the ZH PES is almost invariant with collision energies. The distributions of P(θr) and P(φr) derived from the title reaction indicate that the product polarization is strong. Validity of the QCT calculation has been examined and proven in the comparison with the quantum-wave-packet calculation results. Comparisons with available quasi-classical trajectory results are made and discussed.

Accurate study on the quantum dynamics of the He + HeH(+) (X1Σ+) reaction on a new ab initio potential energy surface for the lowest 1(1)A' electronic singlet state

A time-dependent quantum wave packet method is used to investigate the dynamics of the He + HeH(+)(X(1)Σ(+)) reaction based on a new potential energy surface [Liang et al., J. Chem. Phys.2012, 136, 094307]. The coupled channel (CC) and centrifugal-sudden (CS) reaction probabilities as well as the total integral cross sections are calculated. A comparison of the results with and without Coriolis coupling revealed that the number of K states N(K) (K is the projection of the total angular momentum J on the body-fixed z axis) significantly influences the reaction threshold. The effective potential energy profiles of each N(K) for the He + HeH(+) reaction in a collinear geometry indicate that the barrier height gradually decreased with increased N(K). The calculated time evolution of CC and CS probability density distribution over the collision energy of 0.27-0.36 eV at total angular momentum J = 50 clearly suggests a lower reaction threshold of CC probabilities. The CC cross sections are larger than the CS results within the entire energy range, demonstrating that the Coriolis coupling effect can effectively promote the He + HeH(+) reaction.

STEREODYNAMICS, MASS EFFECT AND DOCKING FROM O(3P) + HD → OH + D AT Ecol = 0.4–1.0eV AND O(3P) + HD → OD + H AT Ecol = 0.5–1.0eV

Quasiclassical Trajectory (QCT) calculation for O(3P) + HD → OH + D and O(3P) + HD → OD + H at E col = 0.4–1.0 eV and 0.5–1.0 eV, respectively, on the lowest PES 1 3A″ of Kuppermann et al. has been done. Distribution p(ϑr) of azimuthal angles between the relative velocity k of the reactants and rotational angular momentum vector j′ of either OH or OD , p(φr) of polar as well as dihedral angles correlating k - k′ -j′, p(ϑr, φr), and PDDCS dependent upon the scattering angle ϑt of either OH , or OD between k and k′ of the reactants and products, respectively, are presented and discussed. The stereodynamics and isotopic mass effects at the smallest possible collision energies 0.4 eV and 0.5 eV for OH and OD , respectively, are significantly different. The significant mass effect with quotient 1/2 of H/D, at the corresponding collision threshold may be applied for the investigation of docking mechanism, drug modification and delivery.

QUASICLASSICAL TRAJECTORY STUDY OF PRODUCT POLARIZATIONS IN THE H(2S) + NH → N(4S) + H2 REACTION

Based on the DMBE potential surface of the 4 A ″ ground-state, the product rotational polarizations in the title reaction are studied by using quasiclassical trajectory (QCT) calculation method. Three angular distributions of P(θr), P(Φr), P(θr, Φr) and the four polarization-dependent differential cross sections (PDDCSs) were calculated for the collision energy range of 1–20 kcal/mol. The results revealed that the product is backward-scattering and the product rotational angular momentum j′ is aligned and oriented. With the increment of collision energy, the degree of the product alignment and orientation is enhanced, showing the collision energy-dependent behaviors of the product polarizations.

VECTOR CORRELATIONS AND PRODUCT POLARIZATIONS IN THE N(2D) + D2 → ND + D REACTIVE SYSTEM

The vector correlations between products and reagents of the N (2D) + D 2 reaction are investigated by employing quasi-classical trajectory (QCT) calculation on the accurate DMBE potential energy surface (PES) of the 2A″ state. Stereo-dynamic quantities, including the four generalized polarization-dependent differential cross-sections (PDDCSs), the angular distribution P(θr), the dihedral-angle distribution P(φr), as well as the product rotational angular distribution in the polar form of P(θr, φr), are calculated in the center-of-mass (CM) frame. The results indicate that the product rotational angular momentum j′ not only aligns along the y-axis, but also orients to the negative direction of the y-axis. The isotope effect in the context of chemical stereo-dynamics and influences of different versions of ground-state PESs on vector correlations are shown and discussed.

INFLUENCE OF VIBRATIONAL EXCITATION AND COLLISION ENERGY ON THE STEREO DYNAMICS OF THE REACTIONS: N(4S)+H2(v = 0-3,j = 0) → NH(X3Σ-)+H

Quasi-classical trajectory (QCT) calculations have been carried out to study the stereodynamics of the title reactions, using the double many-body expansion (DMBE) potential energy surface (PES) constructed by Poveda [Poveda LA, Varandas AJC, Phys. Chem. Chem. Phys.7:2867, 2005]. Vector correlations, such as the distributions of the polarization-dependent differential cross-sections (PDDCSs), the angular distributions of P(θr), P(ϕr), P(θr,ϕr), and the product alignment parameter P2( cos θr) are reported within the energy range of 25–140 kcal mol-1. The influences of the collision energy and the initial state-selected vibrational excitation are discussed in detail. In addition, Validity of the current QCT calculations is also examined and compared with the revelant results reported by Pascual et al. [Pascual RZ, Schatz GC, Lendvay G, Troya D, J. Phys. Chem. Aਠ106:4125, 2002].

INFLUENCE OF THE COLLISION ENERGY ON STEREODYNAMICS OF THE F + LiH (v = 0, j = 0) → LiF + H REACTION

The stereodynamics of the title reaction based on the ground 2A′ potential energy surface (PES) has been investigated using the method of the quasi-classical trajectory (QCT) at different collision energies (23 kcal/mol, 35 kcal/mol and 46 kcal/mol). The vector properties of the angular momentum (described by the distribution of K - J′P(θr), the dihedral angle distribution of K - K′ - J′P(φr) and the angular distribution P(θr, ϕr)) and the four PDDCSs [(2π/σ)(dσ00/dωt), (2π/σ)(dσ20/dωt), (2π/σ)(dσ22+/dωt), (2π/σ)(dσ21-/dωt)] of the product LiF at each collision energy have been presented, respectively. Further, the collision energy effects on the behavior of the product LiF have been discussed and studied.

Vector correlations in the F + HO → HF + O reaction and its isotopic variant

The F + H(D)O → HF(DF) + O reactions have been studied using quasi-classical trajectory (QCT) calculation method, based on the three different potential energy surfaces (PESs) of Gomez-Carrasco et al. (J Chem Phys 2004, 121:4605; J Chem Phys 2005, 123:114310; Chem Phys Lett 2007, 435:188). Facilitated with the analysis of the QCT results, the pictures for product scattering and product polarizations have been presented to investigate the vector correlations in the two reactions, with effects of isotope substitution and electronic state as well as collision energy being revealed at a chemical stereodynamical level.

COMPARATIVE STUDY OF REACTION RATE CONSTANTS FOR THE NH3 + H → NH2 + H2 REACTION WITH GLOBE DYNAMICS AND TRANSITION STATE THEORIES

The nine-dimension quasi-classical trajectory (QCT) calculations have been carried out for the title reaction with a global potential energy surface (PES) constructed by Corchado and Espinosa-García (J Chem Phys106:4013, 1997). The detailed dynamics calculations cover the specific collision energies falling in the range of 0.62–3.04 eV, which are sufficient to fit the calculated reactive cross-sections into a barrier-type excitation function and to obtain the thermal rate constants. The present QCT rate constants are in good agreement with the recent quantum dynamics (QD) results, both of which are much lower than that of the previous variational transition state theory (VTST).

QUASI-CLASSICAL TRAJECTORY STUDY OF THE REACTIONS N(2D) WITH H2, D2, AND HD

Quasi-classical trajectory (QCT) calculations are carried out for the title reactions on the potential energy surface (PES) of Ho et al.1 Our calculated integral cross-section values have been compared with the recent two quantum mechanics (QM) ones: they are close to those of one QM calculation in the high collision energy range, but they approach to another one in the low collision energy range. The product rotational alignments 〈P2 (J' ⋅ K)〉 have also been calculated.

A TRAJECTORY-FOLLOWING METHOD FOR SOLVING THE STEADY STATE OF CHEMICAL REACTION RATE EQUATIONS

This article gives a trajectory-following method for the steady state of chemical reaction rate equations. In order to avoid wasting unnecessary computing time during the steady-state phase and get roughly accurate solutions during the transient-state phase, this method is realized via adopting the semi-implicit Euler formula as the stepping direction and adaptively adjusting the time-step size by an analogous trust-region technique. Under some standard assumptions, its global convergence analysis and local superlinear convergence analysis are also given. Finally, some numerical experiments of this method, in comparison with the traditional optimization methods and ordinary differential equation methods, are reported. The numerical results show that this trajectory-following method is a promising solution for this class of problems.

THE STEREODYNAMICS STUDY OF THE C(3P) + OH (X2 Π) → CO(X1Σ+) + H(2S) REACTION

The stereodynamics of the title reaction on the ground electronic state X2A' potential energy surface (PES)1 has been studied using the quasiclassical trajectory (QCT) method. The commonly used polarization-dependent differential cross-sections (PDDCSs) of the product and the angular momentum alignment distribution, P(θr) and P(Φr), are generated in the center-of-mass frame using QCT method to gain insight of the alignment and orientation of the product molecules. Influence of collision energy on the stereodynamics is shown and discussed. The results reveal that the distribution of P(θr) and P(Φr) is sensitive to collision energy. The PDDCSs exhibit different collision energy dependency relationship at low and high collision energy ranges.

EFFECTS OF ROTATIONAL AND VIBRATIONAL EXCITATION ON THE STEREODYNAMICS OF THE O (D-1) + HCl → OH + Cl REACTION

The stereodynamics of the title reaction on the ground 1 1A′ potential energy surface (PES) has been studied using quasi-classical trajectory (QCT) method. Collision energy of 6.4 kcal/mol is considered, and vector properties including angular momentum alignment distributions and polarization-dependent differential cross-sections (PDDCS) of the product OH are presented. Furthermore, the influence of reagent rotational excitation and vibrational excitation on the product vector properties has also been studied in the present work. The results indicate that the distribution of the P(θr) and P(ϕr) are sensitively affected by the rotational and vibrational excitation. The rotational excitation decreases the degree of alignment and orientation, while vibrational excitation increases the degree of alignment and orientation. The PDDCS (2π/σ)(dσ20/dωt) and (2π/σ)(dσ22+/dωt) are sensitively influenced by rotational and vibrational excitations, while the PDDCS ((2π/σ)(dσ00/dωt)) and (2π/σ)(dσ21-/dωt) are not. The preference of forward scattering has been found from the results of PDDCS ((2π/σ)(dσ00/dωt)), which is in good agreement with the experimental results.

Kinetics and dynamics of the NH3 + H → NH2 + H2 reaction using transition state methods, quasi-classical trajectories, and quantum-mechanical scattering

On a recent analytical potential energy surface developed by two of the authors, an exhaustive kinetics study, using variational transition state theory with multidimensional tunneling effect, and dynamics study, using both quasi-classical trajectory and full-dimensional quantum scattering methods, was carried out to understand the reactivity of the NH(3) + H → NH(2) + H(2) gas-phase reaction. Initial state-selected time-dependent wave packet calculations using a full-dimensional model were performed, where the total reaction probabilities were calculated for the initial ground vibrational state and for four excited vibrational states of ammonia. Thermal rate constants were calculated for the temperature range 200-2000 K using the three methods and compared with available experimental data. We found that (a) the total reaction probabilities are very small, (b) the symmetric and asymmetric N-H stretch excitations enhance the reactivity, (c) the quantum-mechanical calculated thermal rate constants are about one order of magnitude smaller than the transition state theory results, which reproduce the experimental evidence, and (d) quasi-classical trajectory calculations, which were performed with the main goal of analyzing the influence of the zero-point energy problem on the final dynamics results, reproduce the quantum scattering calculations on the same surface.

Effect of vibrational and rotational excitation on the stereodynamics of the C(1D) + H2 (v, j) → CH + H reaction

The stereodynamics of the title reaction on the ab initio1A′ potential energy surface (PES) (B. Bussery-Honvault, P. Honvault, and J.-M. Launay. 2001. J. Chem. Phys. 115: 10701) at a collision energy of 16 kJ/mol have been studied using quasi-classical trajectory (QCT) method. Vector properties including angular momentum alignment parameters and four polarization-dependent differential cross sections (PDDCS) of the product CH are presented. Furthermore, the influence of reagent vibrational and rotational excitations on the product vector properties have also been studied in the present work. The calculated results indicate that the angle distributions of the CH product are mainly dominated by backward–forward scattering.

Quasi-classical trajectory study of the O(1D) + HF reaction dynamics on 11A′ potential energy surface

The quasi-classical trajectory (QCT) calculations for the title reaction were carried out using the recently developed, accurate potential energy surface (PES) of the [Formula: see text] singlet state of the OHF system The integral cross section and the product rotational alignment factor [Formula: see text] were calculated as a function of collision energy. In addition, I discovered the effect of isotopic substitution on stereodynamics for the title reaction, and the influence of the rotation excitation of the reagent on stereodynamics is also presented. Both the scalar and vector properties of the reaction O(1D) + HF → OH + F(2P) are studied in this paper. It was found that the reaction is mainly controlled by an indirect reaction mechanism, and that the deep noncollinear insertion HOF well has a great impact on the dynamics of the reaction. The conclusions drawn in this paper will draw from references to similar reactions, and provide a theoretical foundation for related experiments.