Quantum and quasiclassical dynamics of C(3P)+H2(1∑+g)→ H(2S)+CH(2Π)reaction:Coriolis coupling effects and stereodynamics

The dynamics of C+H2 → H+CH reaction is theoretically studied using the quasiclassical trajectory and quantum mechanical wave packet methods.The analysis of rea...     

Quasi-classical trajectory study of the isotope effect on the stereodynamics in the reaction H(2S)+CH(X2Π; v=0, j=1)→C(1D)+H2(X1Σg+)

The isotope effect on the stereodynamic properties in the title reaction is investigated by a quasi-classical trajectory (QCT) method on the 1¹A' potential energy surface at a collision energy of 23.06 kcal/mol. The angular distributions P(θ_r), P(ø_r), P(θ_r, ø_r), and the polarization-dependent generalized differential cross sections are calculated, which demonstrate the observable influences on the rotational polarization of the product by the isotopic substitution of H with D.     

Theoretical prediction of the optimal conditions for observing the stereodynamical vector properties of the C(~3P) + OH(X~2Π) → CO(X~1 Σ~+) + 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.     

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.     

Exact quantum dynamics study of the H(2S)+SiH+(X1Σ+) reaction on a new potential energy surface of SiH2+(X2A1)

Based on a new global potential energy surface of SiH₂⁺(X²A₁), the exact quantum dynamical calculation for the H(²S)+SiH⁺(X¹Σ⁺)→H₂+Si⁺ reaction has been carried out by using the Chebyshev wave packet method. The initial state specified (ν_i=0, j_i=0) probabilities, integral cross sections (ICS) and thermal rate constants of the title reaction are calculated. All partial wave contributions up to J=90 are calculated in exact quantum calculation including the full Coriolis coupling (CC) effect. The dynamical behaviors of probabilities, ICSs and rate constants are found to be in accord with an exothermic reaction without potential barrier. By comparing the probabilities of CC with the corresponding centrifugal sudden (CS) approximation ones, it can be concluded that neglecting CC effect will decrease the collision time, increase the amplitude of oscillation and lead to overestimation or underestimation of the reaction probability. For ICSs and rate constants, it is found that the deviation of CC and CS ICSs is small in the most of collision energy range except for the range of 0 eV-0.05 eV, while the deviation of both rate constants is considerable in the temperature range of 16 K-1000 K.     

Direct mapping of insertion reaction dynamics: S(1D)+H2→SH+H

The Doppler-selected time-of-flight method was applied to map out the differential cross sections of the title reaction at two collision energies in a crossed-beam experiment. Roughly symmetric and highly forward–backward peaking angular distributions were observed at both energies. Vibrational structures of the SH product were resolved from the velocity measurements of the counter-product H-atom. Most of the angle-integrated observables can readily be understood on statistical grounds, which suggests that statistics plays the dominant role in determining the outcomes of this prototypical insertion reaction. In terms of more detailed angle-specific reaction attributes, significant discrepancies from statistical considerations were revealed, indicative of some hidden dynamics being buried under the statistical factor. Received: 26 February 2000 / Revised version: 20 April 2000 / Published online: 13 September 2000     

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.     

反应物NO的振动激发对反应C(3P)+NO(X2Π)→CO(X1Σ+)+N(2D)立体动力学性质的影响

采用准经典轨线方法研究碰撞能为0.23 eV时,反应物分子NO在不同初始振动态(v=0～3)下发生在两个电子态(2A″和2A′)势能面上反应C(3P)+NO(X2Π)→CO(X1Σ+)+N(2D)的立体动力学性质.计算反应产物的转动角动量矢量分布(P(θr)和P(?r))以及微分散射截面(P00(ωt),(P)20(ω...     

Direct ab initio dynamics study of rate constants and kinetic isotope effects for C2(A3Π u ) + CH3OH reaction

The hydrogen abstraction of CH₃OH by C₂ (A³Π _u ) has been investigated by direct ab initio dynamics over a wide temperature range 200–3000?K. The potential energy surfaces (PESs) have been constructed at the UCCSD(T)/aug-cc-pVTZ//UMP2/6-311++G(d,p) levels of theory. Two different hydrogen abstractions on the methyl and hydroxyl sites of methanol are considered. For the methyl H-abstraction, it is essentially a hydrogen atom transfer (HAT), whereas the hydroxyl site H-abstraction is better described as a proton coupled electron transfer (PCET) according to the Natural Bond Orbital (NBO) analysis. The results suggest that the methyl site reaction is dominant, and the calculated rate constants are roughly consistent with available experimental values. On the other hand, the temperature dependence of deuterium kinetic isotope effects (KIEs) analysis reveals a substantial normal isotope effect in the methyl H-abstraction process, while normal and inverse KIEs coexist in the hydroxyl H-abstraction channel. Furthermore, the three and four–parameter expressions of Arrhenius rate constants are also provided within 200–3000?K.     

气相反应N(4S)+O2(X3∑g-)→NO(X2Π)+O(3P)的准经典轨线研究

基于从头算数据拟合得到的新的基态势能面,对气相反应N(4S) O2(X3∑g-)→NO(X2Π) O(3P)进行了系统的准经典轨线研究.文中计算了该反应的反应几率、反应截面以及在T=300、500、700 K时的热速率常数.通过与以前的理论值和实验值比较发现,我们算得的结果与实验值符合得很好.     

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.     

基于新CH2(X-3A″)势能面的C(3P)+H2(X1∑g^+)→H(2S)+CH(2Π)反应量子波包动力学

基于一个最新的CH2(X-3A″)势能面,运用切比雪夫波包方法对初始态为(v=0,j=0)的C(3P)+H2(X1∑g^+)→H(2S)+CH(2Π)反应体系在1.0-2.0 eV的碰撞能量范围内进行了动力学研究.通过对角动量量子数J=60以下的所有分波进行计算,得到了反应几率、积分散射截面和速率常数.计算中用到了耦合态近似方法和考虑科里奥利耦合效应的精确量子方法.通过对比发现,随着角动量量子数以及能量的增加,科里奥利耦合效应的影响越发显著,因而对于该反应体系,科里奥利耦合效应不可忽略.本文计算所得的积分散射截面和速率常数尚无实验数据可以比较,对该反应的后续研究有一定的参考价值.     

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.     

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.     

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.     

Surface corrugation effects on the adsorption dynamics of xenon on Pt(110)−(1 × 2)

The adsorption probability of xenon on Pt(110)(1 × 2) shows a strong dependence on both the polar (θ_i) and aximuthal angles of incidence. For atoms incident along the (11̄0) direction the adsorption probability increases slowly with increasing angle of incidence, scaling as E_i, cos^0.5θ, whereas along the (100) azimuth it decreases with increasing angle of incidence, exhibiting “negative” energy scaling. The hard cube washboard model does not account quantitatively for these differences with a single consistent set of parameters.     

基于新CH~2(■)势能面的C(~3P)+H_2(■)→H(~2S)+CH(~2Π)反应量子波包动力学

基于一个最新的CH■_2势能面,运用切比雪夫波包方法对初始态为(v=0,j=0)的■反应体系在1.0—2.0 eV的碰撞能量范围内进行了动力学研究.通过对角动量量子数J=60以下的所有分波进行计算,得到了反应几率、积分散射截面和速率常数.计算中用到了耦合态近似方法和考虑科里奥利耦合效应的精确量子方法.通过对比发现,随着角动量量子数以及能量的增加,科里奥利耦合效应的影响越发显著,因而对于该反应体系,科里奥利耦合效应不可忽略.本文计算所得的积分散射截面和速率常数尚无实验数据可以比较,对该反应的后续研究有一定的参考价值.     

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.