Quantum scattering LCAC-SW theory studies on reaction probabilities of three-dimensional H + H2 (v, j) → H2 (v′, j′) + H reaction

Upon the Liu, Siegbahn, Truhlar, Horowitz (LSTH) potential energy surface, the reaction probabilities of the three-dimensional (3-D) state-to-state H + H₂ (v, j) →H ₂(v′, j′) + H reaction are calculated with the linear combination of arrangement channels-scattering wavefunction (LCAC-SW) method. In the calculation, the vibration function of H₂ and the radial propagating wave functions are expanded by the real Gauss functions. The calculated threshold energy and the resonating structure are consistent with the results of the accurate quantum scattering calculations, which shows the accuration, simplicity and practicability of the LCAC-SW method. Project supported by the National Natural Science Fondation of China and the Doctoral Foundation of the State Education Commission of China.     

Energy resonance and inverse population of the product vibrational states for the reaction F + H2(v = 0) → HF(v′) + H,LCAC‐SW theoretical quantum scattering study

LCAC‐SW (linear combination of arrangement channel‐scattering wavefunction) method was used to calculate collinear state‐to‐state reaction probabilities for the reaction F + H₂(v = 0) → HF(v′) + H on the 6SEC potential energy surface. The results show that reaction probabilities P₀₂ and P₀₃ [i. e., v′ = 2,3 for reaction F + H₂ (v = 0) + HF(v′) + H] are primary, the population of product vibrational states is inverse and the reaction probabilities are oscillatory with collision energies, i.e., there is energy resonance in this reaction, which agrees with a new experiment.     

Molecular beam scattering studies of the reaction D + H2 (v = 0) and D + H2 (v = 1) → HD + H

The reaction D + H₂ → HD + H has been investigated in two molecular beam scattering experiments. Angular and time-of-flight distributions have been measured for the initial vibrational ground state (v = 0) at a most probable collision energy of E_cm = 1.5 eV and for the first vibrational excited state (v = 1) at E_cm = 0.28 eV with the same apparatus. Results for the ground-state experiment are compared with quasiclassical trajectory calculations(QCT) on the LSTH-hypersurface transformed into the laboratory system and averaged over the apparatus distributions. The agreement isquite satisfactory. At this high collision energy the HD products are no longer scattered in a backward direction but in a wide angular region concentrated about θ = 90° in the center-of-mass system. The absolute reactive cross section has been determined and the agreement with the theoretical value from QCT calculations is within the experimental error. The high sensitivity of the experiment to different properties of the doubly differential cross section has also been demonstrated. A preliminary evaluation of the experiment with initial vibrational excitation (v = 1) shows that the HD-product molecules are preferably backward scattered and the change of internal energy is small supporting the concept of a reaction which is adiabatic with respect to the internal degrees of freedom.     

氟原子与氢分子共线反应几率的量子散射计算

基于最新的6SEC势能面,用邓从豪等提出的LCAC-SW方法计算得到了共线反应F+H~2(v=0)→HF(v')+H的态-态反应几率,计算结果准确地反映出势能面的特点,进一步证明LCAC-SW方法是一成功的量子散射方法。     

Quantum wave packet study of N(2D) + H2 reactive scattering

N(²D) + H₂ → NH + H reaction at zero total angular momentum is studied by using a time dependent quantum wave packet method. State‐to‐state and state‐to‐all reactive scattering probabilities for a broad range of energy are calculated. The probabilities show many sharp peaks that ascribed to reactive scattering resonances. The probabilities for J > 0 are estimated by using the J‐shifting method. The integral cross sections and thermal rate constants are then calculated. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Quasiclassical trajectory study of the reaction NH(X~3∑~-)+H→N(~4S)+H_2

The dynamics of the NH + H→N+H2 reaction has been investigated by means of the 3D quasiclassical trajectory approach by using the LEPS potential energy surface.The calculated rate coefficient is in good agreement with the experimental value.The reaction was found to occur via a direct channel.The product H2 has a cold excitation of rotational state,but has a reverse distribution of the vibrational state with a peak at v=1.Based on the potential energy surface and the trajectory analysis,the reaction mechanism has been explained successfully.     

Electronic nonadiabatic transitions in the reactive (Ar+ + H2, Ar + H+2, ArH+ + H) system. Numerical results for the collinear configuration

In this communication are presented exact quantum mechanical nonadiabatic electronic transition probabilities for the collinear reaction Ar⁺ + H₂(v_i = 0) → ArH⁺(v_f) + H. The calculations were performed using a potential surface calculated by the DIM method. It is established that large probabilities (≈ 1.0) can be obtained only if there is enough translational energy to overcome a potential barrier formed due to the crossing between v_i = 0 of the Ar⁺ + H₂ system and v_i = 2 of the Ar + H⁺₂ system. The threshold for the reaction is found to be 0.06 eV.     

Experimental study of the D + H2 (v = 1) reaction by CARS spectroscopy

The reactions D + H₂ (v = 0, 1) → HD (v = 0, 1) + H have been studiedin a discharge flow reactor by CARS-spectroscopy. For H₂(v = 0) molecules a rate constant of (4, 0 ± 1, 0) 10^?16 cm³ s^?1 is obtained at 310 K from measured HD (v = 0, 1) product yields. Keeping the degree of vibrational excitation of H₂in the microwave discharge in the range of 1% from the increase of the HD (v = 0, 1) CARS signals a rate of k_{2a, b} = (1, 0 ± 0, 4) 10^?13cm³ s^?1 is derived. The total consumption of H₂ (v = 1) in the presence of D atoms gives a rate k₂ = (1, 9 ± 0, 2) 10^?13 cm³ s^?1 at 310 K. The resultsare discussed in regard to previous measurements and theoretical treatments.     

Experimental measurement of the rate of the reaction o(3P) + H2(v = 0) → OH(v = 0) + H at T = 298 K

Measurement of the rate of the reaction is reported. The measurements were made in a flow tube apparatus. The result is based on data for the absolute density of OH(v = 0) obtained from laser-induced fluorescence measurements in the (0–0) band of the OH(A²Σ⁺ → X²II) system. The density of oxygen atoms was varied by changing the flow rate of NO which is consumed in the reaction N + NO → O + N₂. We find that k₁ (298 K) = (5.5 ± 3.0) × 10⁶ cm³/mol sec. This result was obtained with consideration and control of the effect of reaction (2): for which vibrationally excited hydrogen is created by energy transfer in the presence of active nitrogen. It was found that the addition of N₂ or CO₂ effectively suppressed the excitation of H₂(v = 1). Measurements of the density of H₂(v = 1) were made by VUV absorption in the Lyman band system of H₂. All of the reports of low-temperature measurements and some recent theoretical calculations for k₁ are discussed. The present result confirms and extends the growingevidence for significant curvature in the low-temperature end of a modified Arrhenius plot of k₁ (T).     

Competition between dissociation and exchange processes in a collinear A + BC collision. I. Exact quantum results

Exact quantum results for collision-induced dissociation on a reactive surface are presented. A modified LEPS potential-energy surface modeling the H + HD → H₂ + D system has been used. HD and H₂ bearing respectively 7 and 6 bound states. This system has been chosen because it displays significant reactive scattering for total energies above the dissociation threshold. Calculations have been performed using the time-dependent wavepacket method for two initial vibrational quantum numbers of the HD molecule (v = 0, 2). For each vibrational quantum number, two wavepackets with overlapping energy distributions have been run, covering a total energy range up to more than three times the dissociation energy. Comparison with previous collision-induced dissociation calculations shows that the dissociation is greatly enhanced by the presence of concomitant reactive scattering. A vibrational enhancement effect is also observed above the dissociation threshold; for higher energies the system exhibits a pronounced vibrational inhibition effect.     

Temperature dependence of the rate constant of ionic analogs of H + H2 → H2 + H reaction by quantum instanton method

Rate constant ratios k(T)/k(1,500K) for two symmetrical reactions H^? + H₂ → H₂ + H^? and H⁺ + H₂ → H₂ + H⁺ are reported. Direct method based on quantum instanton approximation for evaluation of the temperature dependence of the quantum‐mechanical reaction rate constant is used. Implementation of the theory involves thermodynamic integration and path integral Monte Carlo method. Results of anionic case shows resemblance to neutral case, whereas cationic case is significantly different and below 1,000K rate constant shows strong deviation form linearity of Arrhenius plot due to high activation barrier. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

A DWBA study of angular distributions for the state-to-state reactive scattering process of F + H2 → HF + H

Based on our earlier three-dimensional DWBA theory, we discuss angular distributions and the roles of various angular momenta. The theory predicts a relatively small number of partial waves, backward scattering, and broadening of angular distributions with increased energy, for the reaction F+ H₂ (v_a = 0, j_a = 0) → HF (v_b = 2 j_b = 0) + H     

Distorted-wave calculations for the three dimensional chemical reaction H + H2(v ⩽ 2, j = 0) → H2(v′ ⩽ 2, j′, mj) + H

Calculations of the vibrational—rotational product state population distributions and differential cross sections for the chemical reaction H + H₂(v ? 2, j = 0) → H₂(v′ ? 2, j′, m_j) + H have been carried out on the Porter—Karplus potential energy surface. The vibrationally-adiabatic-distorted-wave (VADW) method has been used. The relative rotational product distributions, differential cross sections and the helicity m_j, dependences of these quantities for the v = 0 reaction agree well with accurate close-coupling results. The absolute integral cross sections are considerably smaller than the accurate quantum values, however. The calculations for the v = 1 reaction agree with the findings of previous sudden quantum, limited close-coupling and quasiclassical theoretical studies and experiments that product H₂(v′ = 1) is more likely to be produced than H₂(v′ = 0). For the reaction with v = 2, it is found that at high translational energies product H₂(v′ = 2) is favoured over H₂(v′ = 1) or H₂(v′ = 0). The VADW differential cross sections for the v = 1 and v = 2 reactions have a similar shape to those of the v = 0 reaction, with backward peaking when summed over all m_j states. The relative rotational distributions for the v = 2, j = 0 → v′ = 2, j and v = 1, j = 0 → v′ = 1, j reactions are also similar to those obtained for the v = 0, j = 0 → v′ = 0, j reaction, with low rotational excitation.     

The formation of highly excited H in the reaction H2+(v) + H2 → H + H

A quasiclassical trajectory surface hopping method has been used to study H(v) + H₂ → H + H for v = 0, 3, 7, 10, 13, and 17 with an emphasis on determining the H internal energy and angular momentum distributions for high v. For v = 13 and 17, significant cross sections are found for producing H at energies above its dissociation energy. An average metastable H lifetime of 11.5 ps for v = 13 and 4.7 ps for v = 17 is found, but there is also a much longer lived component to the lifetime distributions that is more important for v = 13 than for v = 17. Some of the longer lived metastables correspond to high angular momentum orbiting states of H, but other sources of metastability are also present.     

A time-dependent wave packet quantum scattering study of the reaction HD+ (v = 0 - 3;j0 = 1) + He --> HeH+(HeD+) + D(H)

Time-dependent wave packet quantum scattering (TWQS) calculations are presented for HD(+) (v = 0 - 3;j(0)=1) + He collisions in the center-of-mass collision energy (E(T)) range of 0.0-2.0 eV. The present TWQS approach accounts for Coriolis coupling and uses the ab initio potential energy surface of Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. For a fixed total angular momentum J, the energy dependence of reaction probabilities exhibits quantum resonance structure. The resonances are more pronounced for low J values and for the HeH(+) + D channel than for the HeD(+) + H channel and are particularly prominent near threshold. The quantum effects are no longer discernable in the integral cross sections, which compare closely to quasiclassical trajectory calculations conducted on the same potential energy surface. The integral cross sections also compare well to recent state-selected experimental values over the same reactant and translational energy range. Classical impulsive dynamics and steric arguments can account for the significant isotope effect in favor of the deuteron transfer channel observed for HD(+)(v<3) and low translational energies. At higher reactant energies, angular momentum constraints favor the proton-transfer channel, and isotopic differences in the integral cross sections are no longer significant. The integral cross sections as well as the J dependence of partial cross sections exhibit a significant alignment effect in favor of collisions with the HD(+) rotational angular momentum vector perpendicular to the Jacobi R coordinate. This effect is most pronounced for the proton-transfer channel at low vibrational and translational energies.     

Accurate DMBE Potential Energy Surface For the N(2D) + H2(1Σ g + ) Reaction Using an Improved Switching Function Formalism

A single-sheeted double many-body expansion (DMBE) potential energy surface is reported for the 1 ² A′′ state of NH₂. To approximate its true multi-sheeted nature, a novel switching function that imposes the correct behavior at the H₂(X ¹Σ _g ⁺)+ N(² D) and NH(X ³Σ^-) + H(² S) dissociation limits has been suggested. The new DMBE form is shown to fit with high accuracy an extensive set of new ab initio points (calculated at the multi-reference configuration interaction level using the full valence complete active space as reference and aug-cc-pVQZ and aug-cc-pV5Z basis sets) that have been semiempirically corrected at the valence regions by scaling the n-body dynamical correlation terms such as to account for the finite basis set size and truncated configuration interaction expansion. A detailed study of the N(² D) ... H₂(X ¹Σ _g ⁺) van der Waals region has also been carried out. These calculations predict a nearly free rigid-rotor with two shallow van der Waals wells of C _2v and C _{∞ v} symmetries. Such a result contrasts with previous cc-pVTZ calculations which predict a single T-shaped van der Waals structure. Except in the vicinity of the crossing seam, which is replaced by an avoided intersection, the fit shows the correct physical behavior over the entire configurational space. The topographical features of the new DMBE potential energy surface are examined in detail and compared with those of other potential functions available in the literature. Amongst such features, we highlight the barrier for linearization (11,802 cm^-1) which is found to overestimate the most recent empirical spectroscopic estimate by only 28 cm^-1. Additionally, the T-shaped N(² D) ... H₂ van der Waals minimum is predicted to have a well depth of 90 cm^-1, being 11 cm^-1 deeper than the C _{∞ v} minimum. The title DMBE form is therefore recommendable for dynamics studies of both non-reactive and reactive N(² D)+H₂ collisions.     

Quantum wave packet and quasiclassical trajectory studies of the reaction H(2S) + CH(X2Π; v = 0, j = 1) → C(1D) + H2(X1 ): Coriolis coupling effects and stereodynamics

The quantum mechanics (QM) and quasiclassical trajectory (QCT) calculations have been carried out for the title reaction with the ground minimal allowed rotational state of CH (j = 1) on the 1 ¹A′ potential energy surface. For the reaction probability at total angular momentum J = 0, a similar trend of the QM and QCT calculations is observed, and the QM results are larger than the latter almost in the whole considered energy range (0.1–1.5 eV). The QCT integral cross sections are larger than the QM results with centrifugal sudden approximation, while smaller than those from QM method including Coriolis coupling for collision energies bigger than 0.25 eV. The quantum wave‐packet computations show that the Coriolis coupling effects get more and more pronounced with increasing of J. In addition to the scalar properties, the stereodynamical properties, such as the average rotational alignment factor <P₂( j′?k )>, the angular distributions P(θ_r), P(?_r), P(θ_r,?_r), and the polarization‐dependent generalized differential cross sections have been explored in detail by QCT approach. © 2013 Wiley Periodicals, Inc.     

Schwingungsspektren und Normalkoordinatenanalyse von N-Silylpyrrolen C4H4NSiX3 (X = Cl,CH3, H)

N-trichlorosilyl pyrrole (I) was prepared and, along with N-trimethylsilyl pyrrole (II), investigated by infrared and RAMAN spectroscopy. A normal coordinate analysis which was performed for pyrrole, I, II, and N-silyl pyrrole, suggests that no pure ?SiN stretching frequency”? occurs. This mode is mixed with ν NC₂ and δ(ring) as well as with v_s SiC₃ and v_s SiCl₃ and contributes to frequencies between ～400 and ～1100 cm^?1. The results are compared with quantum chemical calculations for pyrrole and II.     

Dynamics of H+O2 collisions on anab initio potential energy surface

The dynamics of elementary rate processes for H+O₂ collisions on an ab initio potential energy surface have been simulated by quasiclassical trajectory theory (QCT). For H+O₂ (v=0,j=1), we have obtained the reaction probabilityP _r (E,b) as a function of collision energy E and impact parameterb, the reaction cross sectionS _r as a function ofE, and the average values of the product quantum numbers of OH.For H+0₂ (v=2,j=1, 20, 40, 60, 80, 100;v=1, 3, 4, 5,j=1) atE=0.3 eV, we have found thatb _max is about 4.5a ₀ and the impact parameter at whichP _r is maximum decreases asj increases. The reaction cross section increases asj andv become large. For inelastic collisions, whenb is small andj is large, the and are both small. For reactive collisions, almost equals zero, but the probability of being larger than zero increases with increasingj; and¯v _OH even shows population inversion forj=100. Additional details of the dynamics are shown in figures of interparticle distance and stereographs.     

Collaborative control of branching ratio in the O + HO2 → OH + O2 reaction via vibrational and rotational excitation

To understand the effect of different vibrational and rotational modes of reactant to enhance the reactivity of the O + HO₂ → OH + O₂ reaction, we revisited this important atmospheric reaction. We report here a quasi-classical trajectory (QCT) study of the reaction dynamics on a recently developed full-dimensional potential energy surface (PES). Our previous work has indicated that this reaction has two pathways, the H abstraction (HA) channel and the O abstraction (OA) channel, which lead to totally different product energy distribution. In this work, we identified that the vibrational excitation of the OH stretching (v₁) mode of HO₂ is the switch of the HA channel at low collision energy; meanwhile, the rotational excitation can also greatly change the branching ratio of the two pathways. With the excitation of v₁ mode, the original negligible HA channel controlled by the tight transition state becomes quite important. This work presents an approach to control the branching ratio via collaboration between vibrational and rotational excitation and will enrich the knowledge of the O + HO₂ reaction in atmospheric chemistry and physics.