On the role of dynamical barriers in barrierless reactions at low energies: S(1D) + H2

Reaction probabilities as a function of total angular momentum (opacity functions) and the resulting reaction cross sections for the collision of open shell S((1)D) atoms with para-hydrogen have been calculated in the kinetic energy range 0.09-10 meV (1-120 K). The quantum mechanical hyperspherical reactive scattering method and quasi-classical trajectory and statistical quasi-classical trajectory approaches were used. Two different ab initio potential energy surfaces (PESs) have been considered. The widely used reproducing kernel Hilbert space (RKHS) PES by Ho et al. [T.-S. Ho, T. Hollebeek, H. Rabitz, S. D. Chao, R. T. Skodje, A. S. Zyubin, and A. M. Mebel, J. Chem. Phys 116, 4124 (2002)] and the recently published accurate double many-body expansion (DMBE)/complete basis set (CBS) PES by Song and Varandas [Y. Z. Song and A. J. C. Varandas, J. Chem. Phys. 130, 134317 (2009)]. The calculations at low collision energies reveal very different dynamical behaviors on the two PESs. The reactivity on the RKHS PES is found to be considerably larger than that on the DMBE/CBS PES as a result of larger reaction probabilities at low total (here also orbital) angular momentum values and to opacity functions which extend to significantly larger total angular momentum values. The observed differences have their origin in two major distinct topographic features. Although both PESs are essentially barrierless for equilibrium H-H distances, when the H-H bond is compressed the DMBE/CBS PES gives rise to a dynamical barrier which limits the reactivity of the system. This barrier is completely absent in the RHKS PES. In addition, the latter PES exhibits a van der Walls well in the entrance channel which reduces the height of the centrifugal barrier and is able to support resonances. As a result, a significant larger cross section is found on this PES, with marked oscillations attributable to shape resonances and/or to the opening of partial wave contributions. The comparison of the results on both PESs is illustrative of the wealth of the dynamics at low collision energy. It is also illuminating about the difficulties encountered in modeling an all-purpose global potential energy surface.     

State-to-state reaction probabilities for the H+O2(v,j)-->O+OH(v',j') reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.     

Coriolis coupling effects in the dynamics of deep well reactions: application to the H(+) + D2 reaction

We present exact and estimated quantum differential and integral cross sections as well as product state distributions for the title reaction. We employ a time-dependent wavepacket method including all Coriolis couplings and also an adapted code where the helicity quantum number and with this the Coriolis couplings have been truncated. Results from helicity truncated as well as helicity conserving (HC) calculation are presented. The HC calculations fail to reproduce the exact results due to the influence of the centrifugal barrier. While the truncated calculation overestimate the exact integral cross sections they reproduce the features of the integral cross section very well. We also find that the product rotational state distributions are well reproduced if the maximum helicity state is chosen carefully. The helicity truncated calculations fail to give a good approximation of differential cross sections.     

Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces

The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.     

Quasiclassical trajectory studies of H(D)+HBr(DBr) abstraction and exchange reactions

The abstraction and exchange reaction dynamics for H(D)+HBr(DBr) systems have been investigated on three LEPS potential-energy surfaces whose features are in accord with the surface topography suggested by recent molecular-beam and thermal experiments (abstraction barrier less than 1.0 kcal/mole, exchange reaction barriers of =5.0 kcal/mole, and no attractive wells with a depth greater than 0.209 kcal/mole). The surfaces differ primarily in the magnitude of the abstraction barrier which varies from 0.19 to 1.01 kcal/mole. Reaction cross sections have been computed on each surface as a function of relative collision energy from the results of 139000 quasiclassical trajectoris. Comparison of these results with measured relative abstraction cross sections suggests that the true abstraction barrier is very small, perhaps between 0.0 and 0.25 kcal/mole. However, thermal rate coefficients computed on the - best- surface at 300 K are about a factor of 2 larger than the most recently measured values. The calculated (H,D)/(D,H) isotope ratio at 300 K lies between the two reported experimental results. The computed thermal activation energy for abstraction is 835 cal/mole, which is in good agreement with a very early measurements but a factor of 2.5 less than the most recently reported experimental result. These results suggest that the molecular-beam and thermal rate measurements are inconsistent. The average fraction of the available energy which is partitioned into internal product modes <f_E> is found to be nearly independent of relative collision energy and the small topographical differences present in the potential surfaces used in these calculations. We find <f_E> = 0.40. In all reactions, the differential scattering cross sections are peaked in the backward direction for the molecular products, indicating a rebound mechanism.     

Reduced dimensionality quantum dynamics of Cl + CH4 --> HCl + CH3 on an ab initio potential

We study the reaction Cl + CH(4)--> HCl + CH(3) using a 2-D potential energy surface obtained by fitting a double Morse analytical function to high level (CCSD(T)/cc-pVTZ//MP2/cc-pVTZ)ab initio data. Dynamics simulations are performed in hyperspherical coordinates with the close-coupled equations being solved using R-matrix propagation. Quantum contributions from spectator modes are included via a harmonic zero-point correction to the ab initio data prior to fitting the potential. This is the first time this method has been applied to a heavy-light-heavy reaction and the first time it has been used to study differential cross sections. We find thermal rate constants and state-to-state differential cross sections which are in good agreement with experimental data. We discuss the applicability of our method to the study of kinetic isotope effects (KIEs), which we derive for the CH(4)/CD(4) substitution. The calculated KIE compares favourably with experiment. Finally, we discuss the sensitivity of the results of dynamics simulations on the accuracy of the fitted potential.     

Molecular beam scattering of NO+Ne: a joint theoretical and experimental study

The collision dynamics of the NO+Ne system is investigated in a molecular beam scattering experiment at a collision energy of 1055 cm(-1). Employing resonance enhanced multiphoton ionization of NO, we measured state-resolved integral and differential cross sections for the excitation to various levels of both spin-orbit manifolds. The dependence of the scattered intensity on the laser polarization is used to extract differential quadrupole moments for the collision induced angular momentum alignment. The set of cross section data is compared with results of a full quantum mechanical close coupling calculation using the set of ab initio potential energy surfaces of Alexander et al. [J. Chem. Phys. 114, 5588 (2001)]. In previous work, it was found that the positions and rotational substructures for the lowest bend-stretch vibrational states derived from these surfaces agree very well with the observed spectrum of the NO-Ne complex. For the same potential, we find that the calculated cross sections show a less satisfactory agreement with the experimental data. While the overall Jf dependence and magnitude of the integral and differential cross sections are in good agreement, noticeable discrepancies exist for the angle dependence of the differential cross sections. In general, the calculated rotational rainbow structures are shifted towards larger scattering angles indicating that the anisotropy of the potential is overestimated in the fit to the ab initio points or in the ab initio calculation itself. For most states, we find the measured alignment moments to be in excellent agreement with the results of the calculation as well as with predictions of sudden models. Significant deviations from the sudden models are observed only for those fine-structure changing collisions which are dominated by forward scattering. Results of the full quantum calculation confirm the deviations for these states.     

Adaptive temperature-accelerated dynamics

We present three adaptive methods for optimizing the high temperature T(high) on-the-fly in temperature-accelerated dynamics (TAD) simulations. In all three methods, the high temperature is adjusted periodically in order to maximize the performance. While in the first two methods the adjustment depends on the number of observed events, the third method depends on the minimum activation barrier observed so far and requires an a priori knowledge of the optimal high temperature T(high)(opt)(E(a)) as a function of the activation barrier E(a) for each accepted event. In order to determine the functional form of T(high)(opt)(E(a)), we have carried out extensive simulations of submonolayer annealing on the (100) surface for a variety of metals (Ag, Cu, Ni, Pd, and Au). While the results for all five metals are different, when they are scaled with the melting temperature T(m), we find that they all lie on a single scaling curve. Similar results have also been obtained for (111) surfaces although in this case the scaling function is slightly different. In order to test the performance of all three methods, we have also carried out adaptive TAD simulations of Ag/Ag(100) annealing and growth at T = 80 K and compared with fixed high-temperature TAD simulations for different values of T(high). We find that the performance of all three adaptive methods is typically as good as or better than that obtained in fixed high-temperature TAD simulations carried out using the effective optimal fixed high temperature. In addition, we find that the final high temperatures obtained in our adaptive TAD simulations are very close to our results for T(high)(opt)(E(a)). The applicability of the adaptive methods to a variety of TAD simulations is also briefly discussed.     

State-to-state reactive differential cross sections for the H+H2-->H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.     

Differential Cross Sections for the F + H2 Reaction: A Comparison of Classical Trajectory Results for Two Potential Energy Surfaces

Classical trajectories were calculated for the F + H₂ reaction over two potential energy surfaces: (1) the well known Muckerman 5 surface (Theoretical Chemistry: Advances and Perspectives 1981 , 6A, 1) and (2) the recent surface of Stark and Werner (J. Chem. Phys. 1996 , 104, 6515). Integral cross sections, state specified cross sections, differential cross sections and product energy distributions were calculated for the two surfaces. Since the methods for calculating the trajectories and expressing differential cross sections were identical for both surfaces, the rather substantial differences in the results are clearly due to differences in the potential surfaces. The results are discussed in terms of the special characteristics of the two surfaces.     

Accuracy of recent potential energy surfaces for the He-N2 interaction. II. Molecular beam scattering and bulk gas relaxation phenomena

A new semiempirical exchange-Coulomb model potential energy surface for the N(2)-He interaction was reported recently [A. K. Dham et al., J. Chem. Phys. 127, 054302 (2007)] and, using it, the temperature dependence of bulk gas properties of N(2)-He mixtures, such as the second virial coefficient and traditional transport phenomena, most of which depend primarily on the isotropic component of the interaction potential energy surface, was determined. Values of these properties, along with values calculated using two high-quality ab initio potential energy surfaces [C.-H. Hu and A. J. Thakkar, J. Chem. Phys. 104, 2541 (1996); K. Patel et al., ibid 119, 909 (2003)] were compared critically to available experimental data. The present paper reports on the ability of the same three potential energy surfaces to predict state-to-state and total differential cross sections, total integral cross sections, and the temperature dependence of bulk gas relaxation phenomena (including magnetic field effects on transport coefficients). While all three potential energy surfaces give total differential and higher speed integral scattering results that fall within the experimental uncertainties, integral scattering results and state-to-state differential cross section measurements consistently exceed the calculated values. All three surfaces give similar agreement with the relaxation properties of N(2)-He binary mixtures, with the semiempirical exchange-Coulomb model potential energy surface giving slightly better overall agreement with experiment than the two ab initio potential energy surfaces.     

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 hierarchical family of three-dimensional potential energy surfaces for He-CO

A hierarchical family of five three-dimensional potential energy surfaces has been developed for the benchmark He-CO system. Four surfaces were obtained at the coupled cluster singles and doubles level of theory with a perturbational estimate of triple excitations, CCSD(T), and range in quality from the doubly augmented double-zeta basis set to the complete basis set (CBS) limit. The fifth corresponds to an approximate CCSDT/CBS surface (CCSD with iterative triples/CBS, denoted CBS+corr). The CBS limit results were obtained by pointwise basis set extrapolations of the individual counterpoise-corrected interaction energies. For each surface, over 1000 interaction energies were accurately interpolated using a reproducing kernel Hilbert space approach with an R-6+R-7 asymptotic form. In each case, both three-dimensional and effective two-dimensional surfaces were developed. In standard Jacobi coordinates, the final CBS+corr surface has a global minimum at rCO=2.1322a0,R=6.418a0, and gamma=70.84 degrees with a well depth of -22.34 cm-1. The other four surfaces have well depths ranging from -14.83 cm-1 [CCSD(T)/d-aug-cc-pVDZ] to -22.02 cm-1 [CCSD(T)/CBS]. For each of these surfaces the infrared spectrum has been accurately calculated and compared to experiment, as well as to previous theoretical and empirical surfaces. The final CBS+corr surface exhibits root-mean-square and maximum errors compared to experiment (4He) of just 0.03 and 0.04 cm-1, respectively, for all 42 transitions and is the most accurate ab initio surface to date for this system. Other quantities investigated include the interaction second virial coefficient, the integral cross sections, and thermal rate coefficients for rotational relaxation of CO by He, and rate coefficients for CO vibrational relaxation by He. All the observable quantities showed a smooth convergence with respect to the quality of the underlying interaction surface.     

Correlations between dynamical properties and features of potential energy surfaces for the light atom transfer reactions O+HCl→OH+Cl AND Cl+HCl→ClH+Cl

A three-dimensional quasiclassical trajectory study of the dynamics of the light atom transfer reaction O(³P) + HCl(ν=0)→ OH + Cl was carried out employing two LEPS potential energy surfaces (I and II). Attention was focused mainly on three-dynamical properties; the oscillatory behavior of partial cross sections as a function of collision energy; the rotational excitation of the products; and the influence of reagent rotation on reactivity. Distinct differences were found between surfaces I and II with respect to these properties. The examination of individual trajectories indicated that there is a significant difference in the nature of these surfaces. While surface I is governed by weak repulsive forces, surface II is governed by strong attractive forces which tend to direct the reactants toward a collinear geometry. The present results confirm conclusions reached from an earlier study of the reaction Cl+HCl→ClH+Cl concerning correlations between dynamical properties and features of potential energy surfaces. For surfaces of the type that we termed HREP, since they are of repulsive nature and they lead to highly rotationally excited products, no significant oscillations of partial cross sections are obtained and reagent rotation promotes the reaction. On the other hand, for surfaces of the type that we termed COLD (collinearly directing), since they tend to direct the reactants toward a collinear geometry and form rotationally “cold” products, significant oscillations of partial cross sections are obtained and reagent rotation causes a decline in reactivity.     

Collision dynamics of O(3P) + DMMP using a specific reaction parameters potential form

Starting from previous benchmark CBS-QB3 electronic structure calculations (Conforti, P. F.; Braunstein, M.; Dodd, J. A. J. Phys. Chem. A 2009, 113, 13752), we develop two global potential energy surfaces for O((3)P) + DMMP collisions, using the specific reaction parameters approach. Each surface is simultaneously fit along the three major reaction pathways: hydrogen abstraction, hydrogen elimination, and methyl elimination. We then use these surfaces in classical dynamics simulations and compute reactive cross sections from 4 to 10 km s(-1) collision velocity. We examine the energy disposal and angular distributions of the reactive and nonreactive products. We find that for reactive collisions, an unusually large amount of the initial collision energy is transformed into internal energy. We analyze the nonreactive and reactive product internal energy distributions, many of which fit Boltzmann temperatures up to ~2000 K.     

Elastic/inelastic and charge transfer collisions of H++H2 at collision energies of 4.67, 6, 7.3, and 10 eV

Quantum mechanical studies of vibrational and rotational state-resolved differential cross sections, integral cross sections, and transition probabilities for both the elastic/inelastic and charge transfer processes have been carried out at collision energies of 4.67, 6, 7.3, and 10 eV using the vibrational close-coupling rotational infinite-order sudden approach. The dynamics has been performed employing our newly obtained quasidiabatic potential energy surfaces which were generated using ab initio procedures and Dunning's correlation-consistent-polarized quadrupole zeta basis set. The present theoretical results for elastic/inelastic processes provide an overall excellent agreement with the available experimental data and they are also found to be almost similar to that obtained in earlier theoretical results using the ground electronic potential energy surface, lending credence to the accuracy and reliability of the quasidiabatic potential energy surfaces. The results for the complementary charge transfer processes are also presented at these energies.     

Quantum approaches for the insertion dynamics of the H+ + D2 and D+ + H2 reactive collisions

The H(+)+D(2) and D(+)+H(2) reactive collisions are studied using a recently proposed adiabatic potential energy surface of spectroscopic accuracy. The dynamics is studied using an exact wave packet method on the adiabatic surface at energies below the curve crossing occurring at approximately 1.5 eV above the threshold. It is found that the reaction is very well described by a statistical quantum method for a zero total angular momentum (J) as compared with the exact ones, while for higher J some discrepancies are found. For J >0 different centrifugal sudden approximations are proposed and compared with the exact and statistical quantum treatments. The usual centrifugal sudden approach fails by considering too high reaction barriers and too low reaction probabilities. A new statistically modified centrifugal sudden approach is considered which corrects these two failures to a rather good extent. It is also found that an adiabatic approximation for the helicities provides results in very good agreement with the statistical method, placing the reaction barrier properly. However, both statistical and adiabatic centrifugal treatments overestimate the reaction probabilities. The reaction cross sections thus obtained with the new approaches are in rather good agreement with the exact results. In spite of these deficiencies, the quantum statistical method is well adapted for describing the insertion dynamics, and it is then used to evaluate the differential cross sections.     

Communication: Mapping water collisions for interstellar space conditions

We report a joint experimental and theoretical study that directly tests the quality of the potential energy surfaces used to calculate energy changing cross sections of water in collision with helium and molecular hydrogen, at conditions relevant for astrophysics. Fully state-to-state differential cross sections are measured for H(2)O-He and H(2)O-H(2) collisions at 429 and 575?cm(-1) collision energy, respectively. We compare these differential cross sections with theoretical ones for H(2)O+H(2) derived from state-of-the-art potential energy surfaces [P. Valiron et al., J. Chem. Phys. 129, 134306 (2008)] and quantum scattering calculations. This detailed comparison forms a stringent test of the validity of astrophysics calculations for energy changing rates in water. The agreement between theory and experiment is striking for most of the state-to-state differential cross sections measured.