A coupled states quantum reactive scattering study of H + D2 → HD + D ATErel(υ = j = 0) = 0.55 eV

In this paper we present the results of accurate quantum coupled states calculations for the H + D₂ → HD + D reaction at 0.55 eV translational energy (relative to υ = 0, j = 0) using the accurate LSTH potential surface.     

Quantum scattering calculation for reaction Br + H2 on two potential energy surfaces

Three‐dimensional time‐dependent quantum wave packet calculations have been carried out for Br + H₂ on a new global ab initio and a semi‐empirical extended London–Eyring–Polanyi–Sato potential energy surface. It is shown that on the ab initio surface, the threshold energy is much lower, and the reaction probabilities, cross sections, and rate constants are much larger. The effects of the initial rovibrational excitation have also been studied. Comparison of rate constants with experimental measurement implies that the ab initio surface is more suitable for quantum dynamic calculation. The possible reasons and mechanism for the dynamical difference on the two PES are analyzed and discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

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.     

Total scattering cross sections for ethylene by electron impact for incident electron energies from 1 to 2000 eV

We report total scattering cross sections for C₂H₄ molecule by electron impact. Calculations are performed by using two different quantum mechanical methods and they cover the energy range from 1 to 2000 eV. For low energy calculations up to 15 eV, UK molecular R‐matrix code through QUANTEMOL‐N software is used, while intermediate to high energy (15–2000 eV) calculations were carried out by applying spherical complex optical potential formalism. Comparison is made with earlier measurements and theoretical data wherever available. A shape resonance is detected around 2 eV due to the ²B_2g symmetry of an electronic state that corresponds to the temporary negative ion formation of ethylene. The differential cross sections are also calculated for the energy range from 1 to15 eV for the scattering angles between 0º and 180º. © 2013 Wiley Periodicals, Inc.     

Semi-classical calculations of rotational/vibrational transitions in Li+ − H2

Some semi-classical calculations of rot/vib transitions in Li⁺ ? H₂ at 0.869, 1.469 and 3.919 eV total energy are presented. Comparison with recent quantum mechanical and classical S-matrix calculations is made.     

Negative collision energy dependence of Br formation in the OH + HBr reaction

The reaction between HBr and OH leading to H(2)O and Br in its ground state is studied by means of a crossed molecular beam experiment for a collision energy varying from 0.05 to 0.26 eV, the initial OH being selected in the state |JOmega> = |3/2 3/2> by an electrostatic hexapole field. The reaction cross-section is found to decrease with increasing collision energy. This negative dependence suggests that there is no barrier on the potential energy surface for the formation pathway considered. The experimental results are compared with the previously reported quantum scattering calculations of Clary et al. (D. C. Clary, G. Nyman and R. Hernandez, J. Phys. Chem., 1994, 101, 3704), and briefly discussed in the light of skewed potential energy surfaces associated with heavy-light-heavy type reactions.     

A reduced dimensionality quasiclassical and quantum study of the proton transfer reaction H3O(+)+H2O-->H2O+H3O(+)

We report quantum and quasiclassical calculations of proton transfer in the reaction H(3)O(+)+H(2)O in three degrees of freedom, the two OH(+) bond lengths and the OH(+)O angle. The reduced dimensional potential energy surface is obtained from the full dimensional OSS3(p) energy function of H(5)O(2) (+) [L. Ojamae, I. Shavitt, and S. J. Singer, J. Chem. Phys. 109, 5547 (1998)], with an additional long-range correction to reproduce the correct ion-molecule interaction. This surface is used to perform both quasiclassical trajectory and quantum reactive scattering calculations of the zero total angular momentum cumulative reaction probability and cross sections for initial rotational states 0, 1, and 2. Comparison of these quantities are made to assess the importance of quantum effects in this reduced dimensional reaction. Additional quasiclassical cross sections are calculated to obtain the thermal rate constant for the reaction.     

Reduced dimensionality quantum calculations of resonances in H + H2

Reduced dimensionality quantum reaction probabilities are reported for the H+H₂ reaction using the ab initio potential surface of Liu, Siegbahn, Truhlar and Horowitz. Resonances are found for the ground and first two excited adiabatic bending states of H₃. Comparison of the resonance energies with the new coupled states calculations of Colton and Schatz shows good agreement. Additional resonances are reported for energies greater than those considered in the coupled states studies.     

Thermal and state-selected rate constant calculations for O(3p) + H2 → OH + H and isotopic analogs

We present a new parametrization (based on ab initio calculations) of the bending potentials for the two lowest potential energy surfaces of the reaction O(³P) + H₂, and we use it for rate constant calculations by variational transition-state theory with multidimensional semiclassical tunneling corrections. We present results for the temperature range 250–2400 K for both the rate constants and the intermolecular kinetic isotope effects for the reactions of O(³P) with D₂ and HD. In general, the calculated rate constants for the thermal reactions are in excellent agreement with available experiments. We also calculate the enhancement effect for exciting H₂ to the first excited vibrational state. The calculations also provide information on which aspects of the potential energy surfaces are important for determining the predicted rate constants.     

Accurate time dependent wave packet calculations for the N + OH reaction

We present accurate quantum calculations of state-to-state cross sections for the N + OH → NO + H reaction performed on the ground (3)A' global adiabatic potential energy surface of Guadagnini et al. [J. Chem. Phys. 102, 774 (1995)]. The OH reagent is initially considered in the rovibrational state ν = 0, j = 0 and wave packet calculations have been performed for selected total angular momentum, J = 0, 10, 20, 30, 40,...,120. Converged integral state-to-state cross sections are obtained up to a collision energy of 0.5 eV, considering a maximum number of eight helicity components, Ω = 0,...,7. Reaction probabilities for J = 0 obtained as a function of collision energy, using the wave packet method, are compared with the recently published time-independent quantum mechanical one. Total reaction cross sections, state-specific rate constants, opacity functions, and product state-resolved integral cross-sections have been obtained by means of the wave packet method for several collision energies and compared with recent quasi-classical trajectory results obtained with the same potential energy surface. The rate constant for OH(ν = 0, j = 0) is in good agreement with the previous theoretical values, but in disagreement with the experimental data, except at 300 K.     

Fully Coriolis-coupled quantum studies of the H + O2 (upsilon i = 0-2, j i = 0,1) --> OH + O reaction on an accurate potential energy surface: integral cross sections and rate constants

We present accurate quantum calculations of the integral cross section and rate constant for the H + O2 --> OH + O combustion reaction on a recently developed ab initio potential energy surface using parallelized time-dependent and Chebyshev wavepacket methods. Partial wave contributions up to J = 70 were computed with full Coriolis coupling, which enabled us to obtain the initial state-specified integral cross sections up to 2.0 eV of the collision energy and thermal rate constants up to 3000 K. The integral cross sections show a large reaction threshold due to the quantum endothermicity of the reaction, and they monotonically increase with the collision energy. As a result, the temperature dependence of the rate constant is of the Arrhenius type. In addition, it was found that reactivity is enhanced by reactant vibrational excitation. The calculated thermal rate constant shows a significant improvement over that obtained on the DMBE IV potential, but it still underestimates the experimental consensus.     

A combined theoretical and experimental approach to the unimolecular loss of molecular hydrogen from protonated formaldehyde: Determination of the average internal energy of metastable [CH2OH]+ ions

Ab initio quantum chemical calculations (MP2/4–31G**) were performed for the dihydrogen elimination reaction from protonated formaldehyde. The energy difference between reactants and products and the activation energies were found to be in good agreement with the corresponding experimental quantities. Theoretical rate vs. energy curves were computed for a series of isotopic variants of the reaction using the Rice–Ramsperger–Kassel–Marcus (RRKM) method. The vibrational frequencies used in these calculations were taken from the 4–31G** geometry-optimized transition state and reactant structures. Quantum mechanical tunnelling was introduced to explain the existence of metastable CH₂OH ions, and a negative kinetic shift of about 0.1 eV was found. The intramolecular kinetic isotope effect for loss of HH/HD and DH/DD was calculated and compared with the experimental data. The result is consistent with the assumption that the average internal energy of metastable [CH₂OH]⁺ ions is very close to the critical energy for H₂ loss.     

Coriolis-coupled wave packet dynamics of H + HLi reaction

We investigated the effect of Coriolis coupling (CC) on the initial state-selected dynamics of H+HLi reaction by a time-dependent wave packet (WP) approach. Exact quantum scattering calculations were obtained by a WP propagation method based on the Chebyshev polynomial scheme and ab initio potential energy surface of the reacting system. Partial wave contributions up to the total angular momentum J=30 were found to be necessary for the scattering of HLi in its vibrational and rotational ground state up to a collision energy approximately 0.75 eV. For each J value, the projection quantum number K was varied from 0 to min (J, K(max)), with K(max)=8 until J=20 and K(max)=4 for further higher J values. This is because further higher values of K do not have much effect on the dynamics and also because one wishes to maintain the large computational overhead for each calculation within the affordable limit. The initial state-selected integral reaction cross sections and thermal rate constants were calculated by summing up the contributions from all partial waves. These were compared with our previous results on the title system, obtained within the centrifugal sudden and J-shifting approximations, to demonstrate the impact of CC on the dynamics of this system.     

Adiabatic capture theory applied to N+NH-->N2+H at low temperature

The adiabatic capture centrifugal sudden approximation (ACCSA) has been applied to the ground state reaction N+NH-->N2+H over the temperature range 2-300 K using an existent potential energy surface. The resultant thermal rate constants are in agreement with available rate constants from quasi-classical trajectory calculations but are significantly larger than the available experimentally derived rate. The calculated rate constants monotonically increase with increasing temperature but could only be approximately described with a simple Arrhenius-like form. Subtle quantum effects are evident in the initial rotational state resolved cross sections and rate constants.     

Stochastic path-integral method for chemical reaction dynamics: Application to the full 3D H3 system

A stochastic path-integral (SPI) technique for chemical reaction dynamics is explored. It is shown that this technique enables the direct computation of the transition amplitude with a finite space-time range, by generating a set of classical paths subject to simultaneous stochastic differential equations. The numerical values of the Boltzmann matrix elements for a harmonic potential are in good agreement with the analytical ones. Within the quantum transition state theory, the flux-flux autocorrelation function is also evaluated at 630 K for the H + H₂ exchange reaction and is found to give a satisfactory agreement with the previous studies. To appraise the influence of the dimensionality, both one-dimensional Eckart potential and a full three-dimensional (3D) Liu-Siegbahn-Truhlar-Horowitz (LSTH) potential calculations have been performed. The calculated values of the Boltzmann matrix elements for the colinear and the full 3D cases are found to deviate slightly from each other in the lower temperature range. The 3D thermal rate constant is in very good agreement with the previous one. © 1996 John Wiley & Sons, Inc.     

A coupled states reactive scattering study of bending excited resonances in three-dimensional H + H2

Reaction probabilities from coupled states calculations on the Liu-Siegbahn-Truhlar-Horowitz surface for H+H₂ are calculated for the energy range 0.90–1.30 eV. Peaks in the vibrationally inelastic reaction probabilities near 1.10, 1.20 and 1.22–1.24 eV suggest that bending excited resonances labelled by the quantum numbers (11¹0), (12⁰0) and (12²0) exist.     

Dynamical regimes on the Cl + H2 collisions: inelastic rainbow scattering

While Cl + H(2) reactive collisions have been a subject of numerous experimental and theoretical studies, inelastic collisions leading to rotational energy transfer and/or vibrational excitation have been largely ignored. In this work, extensive quantum mechanical calculations covering the 0.5-1.5 eV total energy range and various initial rovibrational states have been carried out and used to perform a joint study of inelastic and reactive Cl + H(2) collisions. Quasiclassical trajectories calculations complement the quantum mechanical results. The analysis of the inelastic transition probabilities has revealed the existence of two distinct dynamical regimes that correlate with low and high impact parameters, b, and are neatly separated by glory scattering. It has been found that while high-b collisions are mainly responsible for |Δj| = 2 transitions which dominate the inelastic scattering, they are very inefficient in promoting higher |Δj| transitions. The effectiveness of this type of collision also drops with rotational excitation of H(2). In contrast, reactive scattering, that competes with |Δj|?> 2 inelastic transitions, is exclusively caused by low-b collisions, and it is greatly favored when the reactants get rotationally excited. Previous studies focusing on the reactivity of the Cl + H(2) system established that the van der Waals well located in the entrance channel play a key role in determining the mechanism of the collisions. Our results prove this to be also a case for inelastic processes, where the origin of the double dynamical regime can be traced back to the influence exerted by this well that shapes the topology of the entrance channel of the Cl-H(2) system.     

Coriolis coupling effects in the calculation of state-to-state integral and differential cross sections for the H+D2 reaction

The quantum wavepacket parallel computational code DIFFREALWAVE is used to calculate state-to-state integral and differential cross sections for the title reaction on the BKMP2 surface in the total energy range of 0.4-1.2 eV with D2 initially in its ground vibrational-rotational state. The role of Coriolis couplings in the state-to-state quantum calculations is examined in detail. Comparison of the results from calculations including the full Coriolis coupling and those using the centrifugal sudden approximation demonstrates that both the energy dependence and the angular dependence of the calculated cross sections are extremely sensitive to the Coriolis coupling, thus emphasizing the importance of including it correctly in an accurate state-to-state calculation.     

Fully converged sudden rotation total integral cross sections for 4He + H2 (ν = 0, j = 0) → 4He + H2 (ν′ = 1, j′)

Total integral cross sections for ⁴He + H₂ (ν = 0, j = 0) → ⁴He + H₂ (ν′ = 1, j′ = 0, 2) have been calculated in the total energy range 1.2 to 5.5 eV, according to a quantal sudden approximation for the H₂ rotational degrees of freedom and a close coupling expansion of the vibrational degree of freedom. Convergence of the above cross sections is investigated by employing four vibration basis sets in the close coupling calculations, i.e., ν = 0,1, ν = 0,1, 2, ν = 0, 1, 2, 3 and ν = 0, 1, 2, 3, 4. Between 4.2 and 5.5 eV calculations were done with three vibration basis sets; ν = 0.–4, ν = 0–5, and ν = 0–6. It is found that at least four vibrational basis functions are needed to converge (to within 5–10%) these cross sections in the above energy range. Comparison of breathing sphere calculations and summed sudden rotation results shows good agreement for the (weakly anisotropic) Mies-Krauss potential. However, as expected the former results underestimate the vibrational 0 → 1 total integral cross sections.     

Kinetics of the hydrogen abstraction R?OH + H → R?O• + H2 reaction class

This paper presents an application of the reaction class transition state theory (RC‐TST) to predict thermal rate constants for the hydrogen abstraction R? OH + H → R? O^? + H₂ reaction class, where R is an alkyl group. We have derived all parameters for the RC‐TST method for this reaction class from rate constants of 19 representative reactions, coupling with linear energy relationships (LERs) and the barrier height grouping (BHG) approach. Error analyses indicate that the RC‐TST/LER, where only reaction energy is needed, and RC‐TST/BHG, where no other information is needed, can predict rate constants for any reaction in this reaction class with satisfactory accuracy for combustion modeling. Specifically for this reaction class, the RC‐TST/LER method has less than 25% systematic errors in the predicted rate constants, whereas the RC‐TST/BHG method has less than 35% error when compared to explicit rate calculations. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 414–429, 2010