Theoretical prediction of energy dependence for D+BrO→DBr+O reaction:The rate constant and product rotational polarization

A quasi-classical trajectory(QCT) calculation is used to investigate the vector and scalar properties of the D + Br O → DBr + O reaction based on an ab initio potential energy surface(X1A state) with collision energy ranging from 0.1 kcal/mol to 6 kcal/mol. The reaction probability, the cross section, and the rate constant are studied. The probability and the cross section show decreasing behaviors as the collision energy increases. The distribution of the rate constant indicates that the reaction favorably occurs in a relatively low-temperature region(T 100 K). Meanwhile, three product angular distributions P(θr), P(φr), and P(θr, φr) are presented, which reflect the positive effect on the rotational angular momentum j' polarization of the DBr product molecule. In addition, two of the polarization-dependent generalized differential cross sections(PDDCSs), PDDCS00 and PDDCS20, are computed as well. Our results demonstrate that both vector and scalar properties have strong energy dependence.     

Time-Dependent Quantum Dynamics of T＋CH4 → CH3＋HT Reaction

Reaction probability, cross section and rate constant are studied for polyatomic reaction T CH4 → CH3 HT using the semirigid vibrating rotor target (SVRT) model. The numerical calculation for the reaction system is carried out using the time-dependent wavepacket method, and the wavepacket is propagated by the splitoperator method. The calculation exhibits a variety of features that can be used for comparison with future experimental investigations. The reaction probability as a function of the translational energy shows slight oscillatory structures, similar to those observed in H abstraction reactions H H2 and H CH4. The comparisons with the H CH4 reaction are described.     

Mechanism analysis of reaction S~+(~2D)+H_2(X~1Σ_g~+) → SH~+(X~3Σ~-) + H(~2S) based on the quantum state-to-state dynamics

We present a state-to-state dynamical calculation on the reaction S~++ H_2→ SH~+ +H based on an accurate ~X2 A~″ potential surface. Some reaction properties, such as reaction probability, integral cross sections, product distribution, etc.,are found to be those with characteristics of an indirect reaction. The oscillating structures appearing in reaction probability versus collision energy are considered to be the consequence of the deep potential well in the reaction. The comparison of the present total integral cross sections with the previous quasi-classical trajectory results shows that the quantum effect is more important at low collision energies. In addition, the quantum number inversion in the rotational distribution of the product is regarded as the result of the heavy–light–light mass combination, which is not effective for the vibrational excitation. For the collision energies considered, the product differential cross sections of the title reaction are mainly concentrated in the forward and backward regions, which suggests that there is a long-life intermediate complex in the reaction process.     

Study of the H＋HS reaction on a newly built potential energy surface using the quasi-classical trajectory method

<正>The quasi-classical trajectory（QCT） method is used to study the H+HS reaction on a newly built potential energy surface（PES） of the triplet state of H₂S（³A″） in a collision energy range of 0-60 kcal/mol.Both scalar properties, such as the reaction probability and the integral cross section（ICS）,and the vector properties,such as the angular distribution between the relative velocity vector of the reactant and that of the product,etc.,are investigated using the QCT method.It is found that the ICSs obtained by the QCT method and the quantum mechanical（QM） method accord well with each other.In addition,the distribution for the product vibrational states is cold,while that for the product rotational states is hot for both reaction channels in the whole energy range studied here.     

State-to-state quantum dynamics of the N(~4S)+H_2(X~1Σ~+) → NH(X~3Σ~-)+H(~2S) reaction and its reaction mechanism analysis

Quantum state-to-state dynamics of the N(4S) + H-2(X1+Σ) → NH(X3Σ) + H(2S) reaction is reported in an accurate novel potential energy surface constructed by Zhai et al.(2011 J. Chem. Phys. 135 104314). The time-dependent wave packet method, which is implemented on graphics processing units, is used to calculate the differential cross sections. The influences of the collision energy on the product state-resolved integral cross sections and total differential cross sections are calculated and discussed. It is found that the products NH are predominated by the backward scattering due to the small impact parameter collisions, with only minor components being forward and sideways scattered, and have an inverted rotational distribution and no inversion in vibrational distributions; both rebound and stripping mechanisms exist in the case of high collision energies.     

State-to-state dynamics of F(~2P) + HO(~2Π) → O(~3P) + HF(~1Σ~+)reaction on 1~3A〞 potential energy surface

State-to-state time-dependent quantum dynamics calculations are carried out to study F(~2P) + HO(~2Π) → O(~3P) +HF(~1Σ~+) reaction on 1~3A〞 ground potential energy surface(PES). The vibrationally resolved reaction probabilities and the total integral cross section agree well with the previous results. Due to the heavy–light–heavy(HLH) system and the large exoergicity, the obvious vibrational inversion is found in a state-resolved integral cross section. The total differential cross section is found to be forward–backward scattering biased with strong oscillations at energy lower than a threshold of 0.10 eV, which is the indication of the indirect complex-forming mechanism. When the collision energy increases to greater than 0.10 eV, the angular distribution of the product becomes a strong forward scattering, and almost all the products are distributed at θ_t = 0°. This forward-peaked distribution can be attributed to the larger J partial waves and the property of the F atom itself, which make this reaction a direct abstraction process. The state-resolved differential cross sections are basically forward-backward symmetric for v' = 0, 1, and 2 at a collision energy of 0.07 eV; for a collision energy of 0.30 eV,it changes from backward/sideward scattering to forward peaked as v increasing from 0 to 3. These results indicate that the contribution of differential cross sections with more highly vibrational excited states to the total differential cross sections is principal, which further verifies the vibrational inversion in the products.     

Novel potential energy surface-based quantum dynamics of ion–molecule reaction O~++ D_2→OD~++ D

According to a novel electronic ground-state potential energy surface of H_2O~+(X~4 A~'),we calculate the reaction probabilities and the integral cross section for the titled reaction O+~+ D_2→OD~++ D by the Chebyshev wave packet propagation method.The reaction probabilities in a collision-energy range of 0.0 e V–1.0 e V show an oscillatory structure for the O~++ D_2 reaction due to the existence of the potential well.Compared with the results of Mart′?nez et al.,the present integral cross section is large,which is in line with experimental data.     

Stereodynamics of the O（3P） with H2（D2） （ν=0,j=0） reaction

Quasi-classical trajectory theory is used to study the reaction of O(3 P) with H 2 (D 2) based on the ground 3 A″ potential energy surface (PES).The reaction cross section of the reaction O+H 2 →OH+H is in excellent agreement with the previous result.Vector correlations,product rotational alignment parameters P 2 (j · k) and several polarizeddependent differential cross sections are further calculated for the reaction.The product polarization distribution exhibits different characteristics that can be ascribed to different motion paths on the PES,arising from various collision energies or mass factors.     

Quasi-classical trajectory investigation on the stereodynamics of Li ＋ DF（v=1-6,j=0）→LiF＋D reaction

In this paper, the stereodynamics of Li ＋ DF → Li F ＋ D reaction is investigated by the quasi-classical trajectory（QCT）method on the ^2A＇ potential energy surface（PES） at a relatively low collision energy of 8.76 kcal/mol. The scalar properties of the title reaction such as reaction probability and cross section are studied with vibrational quantum number of v = 1–6. The product angular distributions P（θr） and P（φr） are presented in the same vibrational level range. Moreover, two polarization-dependent generalized differential cross sections（PDDCSs）, i.e., the PDDCS00 and PDDCS22＋are calculated as well. These stereodynamical results demonstrate sensitive behaviors to the vibrational quantum numbers.     

Stereodynamics of the O（^3p） with H： （D2） （v = O,j =O） reaction

Quasi-classical trajectory theory is used to study the reaction of O（3p） with H2 （D2） based on the ground 3A″ potential energy surface （PES）. The reaction cross section of the reaction O＋H2→＋OH＋H is in excellent agreement with the previous result. Vector correlations, product rotational alignment parameters （P2（j′. k）） and several polarizeddependent differential cross sections are further calculated for the reaction. The product polarization distribution exhibits different characteristics that can be ascribed to different motion paths on the PES, arising from various collision energies or mass factors.     

The effect of collision energy on the stereodynamics of the reaction H(~2S)+NH(X~3∑~-,v = 0,j = 0)→ N(~4S)+H_2

The quasi-classical trajectory(QCT) is calculated to study the stereodynamics properties of the title reaction H(2S)+NH(X3∑-) →N(4S)+H2 on the ground state 4A' potential energy surface(PES) constructed by Zhai and Han [2011 J.Chem.Phys.135 104314].The calculated QCT reaction probabilities and cross sections are in good agreement with the previous theoretical results.The effects of the collision energy on the k-k' distribution and the product polarization of H2 are studied in detail.It is found that the scattering direction of the product is strongly dependent on the collision energy.With the increase in the collision energy,the scattering directions of the products change from backward scattering to forward scattering.The distribution of P(θr) is strongly dependent on the collision energy below the lower collision energy(about 11.53 kcal/mol).In addition,the P(φr) distribution dramatically changes as the collision energy increases.The calculated QCT results indicate that the collision energy plays an important role in determining the stereodynamics of the title reaction.     

Stereodynamics study of the H′（^2S） ＋ NH（X^3∑-） 〉 N（4S） ＋ H2 reaction

The stereodynamics and reaction mechanism of the H′（^2S） ＋ NH （X^3∑^-） → N（^4S） ＋ H2 reaction are thoroughly studied at collision energies in the 0.1 eV-1.0 eV range using the quasiclassical trajectory （QCT） on the ground 4A″ potential energy surface （PES）. The distributions of vector correlations between products and reagents P（φr）, P（φr） and P（φr,φr） are presented and discussed. The results indicate that product rotational angular momentum j′ is not only aligned, but also oriented along the direction perpendicular to the scattering plane; further, the product H2 presents different rotational polarization behaviors for different collision energies. Furthermore, four polarization-dependent differential cross sections （PDDCSs） of the product He are also calculated at different collision energies. The reaction mechanism is analyzed based on the stereodynamics properties. It is found that the abstraction mechanism is appropriate for the title reaction.     

Stereodynamics in reaction O（1D） ＋ CH4 〉OH ＋ CH3

A theoretical study of the stereodynamics for reaction O(1D) + CH4→OH + CH3 has been carried out using the quasiclassical trajectory method(QCT) on a potential energy surface structured by Gonzalez et al. The integral cross sections(ICSs), differential cross sections(DCSs) and product rotational angular momentum polarization have been calculated. With the collision energy increasing, the ICS decreases. There is no threshold energy, because no barrier is found on the minimum energy path. The DCS results show that the backward and forward scatterings exist at the same time. With the collision energy increasing, the dominant rotation of the product changes from the right-handed direction to the left-handed direction in planes parallel to the scattering plane. In the isotopic effect study, the decrease of the mass factor weakens the polarization degree of the rotational angular momentum vectors of the products.     

Thermal Rate Constants of the N（^4S）＋O2（X^3∑g^-） → NO（X^2Ⅱ） ＋O（^3P） Reaction on the ^2A′ Potential Energy Surface

A quasiclassical trajectory study with the sixth-order explicit symplectic algorithm for the N（^4S）＋O2（X^3∑g^-） → NO（X^2Ⅱ） ＋O（^3P） reaction has been reported by employing a new ground potential energy surface. We have discussed the influence of the relative translational energy, the vibrational and rotational levels of O2 molecules on the total reaction cross section. Thermal rate constants at temperatures 300, 600, and 1000 K determined in this work for the reaction are 4.4 × 10^7, 1.8 × 10^10, and 3.1 × 10^11 cm^3mol^-1s^-1, respectively. It is found that they are in better agreement with the experimental data than previous theoretical values.     

Effect of reagent vibrational excitation and isotope substitution on the stereo-dynamics of the Ba+HF→BaF+H reaction

Based on an extended London-Eyring-Polanyi-Sato (LEPS) potential energy surface (PES), the Ba + HF reaction has been studied by the quasi-classical trajectory (QCT) method. The reaction integral cross section as a function of collision energy for the Ba + HF → BaF + H reaction is presented and the influence of isotope substitution on the differential cross sections (DCSs) and alignments of the product's rotational angular momentum have also been studied. The results suggest that the integral cross sections increase with increasing collision energy, and the vibrational excitation of the reagent has great influence on the DCS. In addition, the product's rotational polarization is very strong as a result of heavy-heavy-light (HHL) mass combination, and the distinct effect of isotope substitution on the stereodynamics is also revealed.     

Non-adiabatic quantum dynamical studies of Na(3p) + HD(v = 1, j5) = 0)→NaH/NaD + D/H reaction

Non-adiabatic dynamical calculations are carried out for the Na(3 p)+HD(ν = 1, j = 0)→NaH/NaD+D/H reaction on the diabatic potential energy surfaces of Wang et al.(Sci. Rep. 2018, 8, 17960) by using the time-dependent wave packet method. The state-to-state integral cross sections and differential cross sections of two reaction channels(NaH/NaD+D/H)are calculated for collision energy up to 0.4 eV. The cross section branching ratio indicates that the dominant reaction channel changes from NaD+H to NaH+D when the collision energy is larger than 0.227 eV. The products from two reaction channels both prefer to form in vibrationally cold but rotationally hot states, and they both tend to forward scattering.     

Photodetachment of H- Near a Dielectric Surface

Using the closed orbit theory, the photodetachment cross section of H- near a dielectric surface has been derived and calculated. The results show that the dielectric surface has great influence on the photodetachment process of negative ion near the ionization threshold. Above the ionization threshold, the photodetachment cross section starts to oscillate. With the increase of the energy, the oscillating amplitude decreases and the oscillating frequency increases. The oscillation in the photodetachment cross section of H- in the presence of a dielectric surface is either larger or smaller than the photodetachment of H- without the surface. As the photon energy is larger than the critical value Epc, the oscillatory structure disappeared and the cross section approaches to the case of the photodetachment of H- without any external fields. For a given detached-electron energy, the photodetachment cross section becomes decreased with the increase of the ion-surface distance. Besides, the dielectric constant has great influence on the photodetachment of H-. With the increase of the dielectric constant, the oscillation in the cross section becomes increased. As the dielectric constant increases to infinity, the cross section is the same as the photodetachment of H- near a metal surface. This study provides a new understanding on the photodetachment process of H- in the presence of a dielectric surface.     

Effects of a reagent＇s rotational and vibrational excitations on reactionO（3p）＋ H2（v=0, 3, j = 0, 3, 5, 7, 9, 12, 15） → OH ＋ H

To investigate the effect of a reagent’s rotational and vibrational excitations on the stereo-dynamics of the reaction product, the title reaction is theoretically simulated using the quasi-classical trajectory (QCT) method on the 3 A and 3 A potential energy surfaces (PESs). The reaction cross section is considered as the only scalar property in this work at four different collision energies. Furthermore the vector properties including two polarization-dependent differential cross sections (PDDCSs), the angular distributions of product’ rotational momentum are discussed at one fixed collision energy. Effects of reagents’ rotational excitation on the reaction do exist regularly.     

Theoretical prediction of the optimal conditions for observing the stereodynamical vector properties of the C（3P）＋OH （X2∏）→CO（X1S＋）＋H（2S） reaction

The best optimal initial reactant state and collision energy for observing the stereodynamical vector properties of the title reaction in the ground electronic state X2A’ potential energy surface （PES）[Zanchet et al. 2006 J. Phys. Chem. A 110 12017] are theoretically predicted using the quasi-classical trajectory （QCT） method for the first time. The calculated results reveal that the smallest value of the rotational quantum number j, larger vibrational quantum number v, and the lower strength of collision energy should be selected for offering the most obvious picture about the stereodynamical vector properties. Polarization-dependent differential cross sections and the angular momentum alignment distribution, P（θr） and P（Φr） in the center-of-mass frame, are obtained to gain an insight into the alignment and orientation of the product molecules. The rotational angular momentum vector j’ of CO is aligned to be perpendicular to reagent relative velocity k. The product polarizations align along the y axis, pointing to the positive direction of the y axis. A new method is developed to investigate massive reactions with various initial states and to further study the vector properties of the fundamental reactions in detail.     

Quasi-classical trajectory study of H+LiH(v=0,1,2,j=0)→Li+H_2 reaction on a new global potential energy surface

Quasi-classical trajectory(QCT) calculations are reported for the H+LiH(v = 0–2, j = 0)→Li+H_2 reaction on a new ground electronic state global potential energy surface(PES) of the LiH_2 system. Reaction probability and integral cross sections(ICSs) are calculated for collision energies in the range of 0 eV–0.5 eV. Reasonable agreement is found in the comparison between present results and previous available theoretical results. We carried out statistical analyses with all the trajectories and found two main distinct reaction mechanisms in the collision process, in which the stripping mechanism(i.e., without roaming process) is dominated over the collision energy range. The polarization dependent differential cross sections(PDDCSs) indicate that forward scattering dominates the reaction due to the dominated mechanism. Furthermore,the reactant vibration leads to a reduction of the reactivity because of the barrierless and attractive features of PES and mass combination of the system.