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.     

Energy and rotation-dependent stereodynamics of H(~2S) + NH(a~1?) → H_2(X~1Σ_g~+) + N(~2D) reaction

Quasi-classical trajectory calculations are performed to study the stereodynamics of the H(~2S) + NH(a~1?) →H_2(X~1Σ_g~+) + N(~2D) reaction based on the first excited state NH_2(1~2A') potential energy surface reported by Li et al.[Li Y Q and Varandas A J C 2010 J. Phys. Chem. A 114 9644] for the first time. We observe the changes of differential cross-sections at different collision energies and different initial reagent rotational excitations. The influence of collision energy on the k–k' distribution can be attributed to a purely impulsive effect. Initial reagent rotational excitation transforms the reaction mechanism from insertion to abstraction. The effect of initial reagent rotational excitations on k–k' distribution can be explained by the rotational excitation enlarging the rotational rate of reagent NH in the entrance channel to reduce the probability of collision between incidence H atom and H atom of target molecular. We also investigate the changes of vector correlations and find that the rotational angular momentum vector j' of the product H_2 is not only aligned, but also oriented along the y axis. The alignment parameter, the disposal of total angular momentum and the reaction mechanism are all analyzed carefully to explain the polarization behavior of the product rotational angular moment.     

Theoretical study of stereodynamics for the O（3P） ＋ H2 → OH ＋ H reaction

This paper employs the quasi-classical trajectory calculations to study the influence of collision energy on the title reaction on the potential energy surface of the ground ³A' triplet state developed by Rogers et al. (J. Phys. Chem. A 2000 104 2308). It calculates the product angular distribution of P(θ_r), P(φ_r) and P(θ_r, φ_r) which reflects vector correlation. The distribution P(θ_r) shows that product rotational angular momentum vectors j' of the products are strongly aligned along the relative velocity direction k. The distribution of P(φ_r) implies a preference for left-handed product rotation in planes parallel to the scattering plane. Four different polarisation-dependent cross-sections are also presented in the centre-of-mass frame. Results indicate that OH is sensitively affected by collision energies of H₂.     

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.     

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.     

The effect of the rotational excitation of NO on the stereodynamics for the reaction C(~3P)+NO(X~2Π)→CN(X~2Σ~+)+O(~3P)

The stereodynamic properties of the reaction C(3P)+NO(X2Π)→CN(X2Σ+)+O(3P) in different rotational states of reactant NO are studied theoretically by using the quasiclassical trajectory method on 2A' and 2A' potential energy surfaces(PESs) at a collision energy of 0.06 eV.The vector properties in different rotational states on the two surfaces are discussed in detail.The results indicate that the rotational excitation of NO has considerable influence on the stereodynamic property of the reaction occurring on the two surfaces.At the same time,the calculated polarization-dependent differential cross sections(PDDCSs) in different initial rotational states manifest that products are strongly polarized at three scattering angles.     

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.     

Dynamics of the CH4+ O(~3 P) → CH_3(ν= 0) + OH(■= 0) reaction

The dynamics of the ground-state reaction of CH_4+ O(~3P) → CH_3(ν = 0) + OH( ν= 0) have attracted a great deal of attention both theoretically and experimentally. This rapid communication represents extensive quasi-classical trajectory calculations of the vibrational distributions on a unique full-dimensional ab initio potential energy surface for the title reaction, at the collision energy of relevance to previous crossed molecular beam experiments. The surface is constructed using the all electrons coupled-cluster singles and doubles approach plus quasi-perturbative triple excitations with optimized basis sets. A modified Shepard interpolation method is also employed for the construction. Good agreement between our calculations and the available experimental results has been achieved, opening the door for accurate dynamics on this surface.     

Influence of ro-vibrational and isotope effects on the dynamics of the C(3 P)+ OD(X 2Π) → CO(X 1 Σ+) + D(2 S) reaction

The C(³ P)+OD(X ²Π) reaction has been studied by means of quantum mechanical real wave packet (RWP) and quasiclassical trajectory (QCT) methodologies on the ground potential energy surface of Zanchet et al. [J. Phys. Chem. A 110, 12017 (2006)]. Initial state selected total reaction probabilities at J?=?0 total angular momentum have been calculated for a wide range of collision energies. Product state-resolved integral cross-sections at selected collision energies and excitation functions have been determined from the RWP calculations using the J-shifting approximation and from QCT calculations. State-specific and thermal rate coefficients have been calculated using both methodologies up to 500 K. The effect of reagent rotational excitation on the dynamics for the C(³ P)+OH(X ²Π) and C(³ P)+OD(X ²Π) reactions has been investigated and interesting discrepancies between the QCT and RWP results have been found. The RWP results are found to be in an overall good agreement with the corresponding QCT results, although the QCT integral cross-section and rate coefficients are slightly smaller than those obtained from the RWP calculations.     

Vector correlations study of the reaction N(~2D) + H_2(X~1Σ_g~+) →NH(a~1?) + H(~2S) with different collision energies and reagent vibration excitations

Vector correlations of the reaction N(2D)+ H2(X1Σ+g) → NH(a1?)+ H(2S) are studied based on a recent DMBESEC PES for the first excited state of NH2[J. Phys. Chem. A 114 9644(2010)] by using a quasi-classical trajectory method.The effects of collision energy and the reagent initial vibrational excitation on cross section and product polarization are investigated for v = 0–5 and j = 0 states in a wide collision energy range(10–50 kcal/mol). The integral cross section could be increased by H2 vibration excitation remarkably based on the DMBE-SEC PES. The different phenomena of differential cross sections with different collision energies and reagent vibration excitations are explained. Particularly,the NH molecules are scattered mainly in the backward hemisphere at low vibration quantum number and evolve from backward to forward direction with increasing vibration quantum number, which could be explained by the fact that the vibrational excitation enlarges the H–H distance in the entrance channel, thus enhancing the probability of collision between N atom and H atom. A further study on product polarization demonstrates that the collision energy and vibrational excitation of the reagent remarkably influence the distributions of P(θr), P(φr), and P(θr, φr).     

O (3P) + H2 → OH+H反应的立体动力学研究

用BMS1势能面[Brandão 等, J.Chem.Phys. 121, 8861 (2004)],选取碰撞能为34.6kcal/mol,用准经典轨线方法研究了O ( ) + 反应的立体动力学性质。计算并讨论了k与j'的夹角的分布关系P ( )以及描述k-k'- j'三者关系的二面角分布。（k为反应物速度方向,k'为产物的速度方向,j'为产物的角动量方向） 的峰值在90°附近并且关于90°呈对称性分布,这表明产物角动量的方向与初速度的方向垂直。二面角分布P( )关于散射平面呈反对称性分布,这一结果表明产物的角动量具有强烈的极化效应。另外,我们还研究了振动激发对产物角动量的影响,结果表明产物的矢量性质对反应物的初始振动态非常敏感。     

Theoretical study of the aqueous medium effect on the reaction CO2+ OH− ⇄ HCO−3

A method for the theoretical study of chemical reaction mechanisms in aqueous solutions has been developed. Approximate reaction pathways are determined with the aid of a model Hamiltonian in which medium is described by a set of Langevin point dipoles. Solvation of stationary points on the potential energy surface so obtained is performed by a more exact calculation of the medium contribution using the method of minimization of the solvation shell potential energy. Usefulness of our approach is exemplified by the calculation of the reaction CO₂ + OH⁻ ⇄ HCO⁻₃. Method MIND0/13 has been employed in constructing the model Hamiltonian. A good agreement with experiment has been achieved.     

Stereodynamics study of the exchange reaction O(~3P)+CH_4 →H+OCH_3

A new London-Eyring-Polanyi-Sato potential energy surface is employed in this work to study the stereo properties of the O(3P)+CH4 →H+CH3O reaction in its rovibrationally ground state using the quasiclassical trajectory method(QCT).Our calculations are performed at a range of collision energies,Ec=1.5eV～3.5eV,and the excitation function obtained by the QCT method accords well with the experimental data.The product rotational polarization is calculated,and the product shows a strong rotational polarization in the centre-of-mass coordinate system.The orientation of the product rotational angular momenta is sensitive to the increase in collision energy,and the alignment of the product rotational angular momenta shows some of the properties of the heavy heavy-light mass combination reactions.In the isotopic substituted reaction study,when the H atoms in methane are replaced by D atoms,the rotational polarization is obviously reduced.The polarization-dependent differential cross section is also studied by this QCT calculation to provide detailed information about the rotational alignment and orientation of the product.     

Ab initio investigation of the gas phase reaction SO2 + H2O → SO2OH radical

High level ab initio calculations have been performed to investigate the reaction of the OH radical with SO₂. This reaction has been suggested as a possible first step in the atmospheric oxidation of SO₂. Results from both density functional theory (DFT) and second-order M?ller-Plesset calculations are reported. A small barrier has been located in the reaction channel and the structure of the corresponding transition state characterized. On the potential energy surface, the product, HOSO₂, occupies a double well potential that corresponds to a pair of equal energy rotamers separated by a barrier of ～12 kJ mol^?1 corresponding to a symmetric transition state.     

Kinetics and dynamics of the S(1D2) + H2 → SH + H reaction at very low temperatures and collision energies

We report combined studies on the prototypical S(1D2) + H2 insertion reaction. Kinetics and crossed-beam experiments are performed in experimental conditions approaching the cold energy regime, yielding absolute rate coefficients down to 5.8 K and relative integral cross sections to collision energies as low as 0.68 meV. They are supported by quantum calculations on a potential energy surface treating long-range interactions accurately. All results are consistent and the excitation function behavior is explained in terms of the cumulative contribution of various partial waves.     

Vibrational and rotational excitation effects of the N(2D) + D2(X1Σg +) → ND(X3Σ+) + D(2S) reaction

The effects of the rovibrational excitation of reactants in the N(²D) + D₂(X¹Σ_g⁺) → ND(X³Σ⁺) + D(²S) reaction are calculated in a collision energy range from the threshold to 1.0 eV using the time-dependent wave packet approach and a second-order split operator. The reaction probability, integral cross-section, differential cross-section and rate constant of the title reaction are calculated. The integral cross-section and rate constant of the initial states v = 0, j = 0, 1, are in good agreement with experimental data available in the literature. The rotational excitation of the D₂ molecule has little effect on reaction probability, integral cross-section and the rate constant, but it increased the sideways and forward scattering signals. The vibrational excitation of the D₂ molecule reduced the threshold and broke up the forward–backward symmetry of the differential cross-section; it also increased the forward scattering signals. This may be because the vibrational excitation of the D₂ molecule reduced the lifetime of the intermediate complex.     

The catalytic effects of H2CO3, CH3COOH,HCOOH and H2O on the addition reaction of CH2OO + H2O → CH2(OH)OOH

The addition reaction of CH₂OO + H₂O → CH₂(OH)OOH without and with X (X = H₂CO₃, CH₃COOH and HCOOH) and H₂O was studied at CCSD(T)/6-311+ G(3df,2dp)//B3LYP/6-311+G(2d,2p) level of theory. Our results show that X can catalyse CH₂OO + H₂O → CH₂(OH)OOH reaction both by increasing the number of rings, and by adding the size of the ring in which ring enlargement by COOH moiety of X inserting into CH₂OO···H₂O is favourable one. Water-assisted CH₂OO + H₂O → CH₂(OH)OOH can occur by H₂O moiety of (H₂O)₂ or the whole (H₂O)₂ forming cyclic structure with CH₂OO, where the latter form is more favourable. Because the concentration of H₂CO₃ is unknown, the influence of CH₃COOH, HCOOH and H₂O were calculated within 0–30 km altitude of the Earth's atmosphere. The results calculated within 0–5 km altitude show that H₂O and HCOOH have obvious effect on enhancing the rate with the enhancement factors are, respectively, 62.47%–77.26% and 0.04%–1.76%. Within 5–30 km altitude, HCOOH has obvious effect on enhancing the title rate with the enhancement factor of 2.69%–98.28%. However, compared with the reaction of CH₂OO + HCOOH, the rate of CH₂OO···H₂O + HCOOH is much slower.     

Force and density study of the chemical reaction process OH2+H+→OH3 +

The force and density origin of the ionic reaction process OH₂+H⁺→OH₃ ⁺ is studied on the basis of the electrostatic force theory. In order to obtain reliable results, we have used the floating AO wavefunction (with minimal basis) which satisfies the Hellmann-Feynman theorem for the approaching proton. Accuracies of the results of the ordinary nuclear-centred AO wavefunction (with minimal basis plus polarization functions on hydrogen atoms) have been also examined. It is shown that the predominant origin of the driving force is the atomic dipole force which is caused by the transfer of the lone-pair density of OH₂ into the atomic region of H⁺ and the succeeding inward polarization of the transferred electron density. The exchange force due to density accumulation in the forming O-H⁺ bond region is of secondary importance. This is in contrast with the previous results for the radical reaction processes. The present reaction may be regarded as the transformation of the lone pair into the bond pair.