Quasiclassical trajectory study of the F + CH4 reaction dynamics on a dual-level interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the F + CH4 reactive system has been constructed using the interpolation method of Collins and co-workers. The ab initio calculations have been performed using second-order M?ller-Plesset (MP2) perturbation theory to build the initial PES. Scaling all correlation (SAC) methodology has been employed to improve the ab initio calculations and to construct a dual-level PES. Using this PES, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations and internal energy distributions has been carried out for the F + CH4 and F + CD4 reactions and the theoretical results have been compared with the available experimental data.     

He-ThO(1Σ+) interactions at low temperatures: elastic and inelastic collisions, transport properties, and complex formation in cold 4He gas

We present an ab initio study of cold (4)He + ThO((1)Σ(+)) collisions based on an accurate potential energy surface (PES) evaluated by the coupled cluster method with single, double, and noniterative triple excitations using an extended basis set augmented by bond functions. Variational calculations of rovibrational energy levels show that the (4)He-ThO van der Waals complex has a binding energy of 10.9 cm(-1) in its ground J = 0 rotational state. The calculated energy levels are used to obtain the temperature dependence of the chemical equilibrium constant for the formation of the He-ThO complex. We find that complex formation is thermodynamically favored at temperatures below 1 K and predict the maximum abundance of free ground-state ThO(v = 0, j = 0) molecules between 2 and 3 K. The calculated cross sections for momentum transfer in elastic He + ThO collisions display a rich resonance structure below 5 cm(-1) and decline monotonically above this collision energy. The cross sections for rotational relaxation accompanied by momentum transfer decline abruptly to zero at low collision energies (<0.1 cm(-1)). We find that Stark relaxation in He + ThO collisions can be enhanced by applying an external dc electric field of less than 100 kV∕cm. Finally, we present calculations of thermally averaged diffusion cross sections for ThO in He gas, and find these to be insensitive to small variations of the PES at temperatures above 1 K.     

Quasi-classical trajectory study of the S + OH → SO + H reaction: from reaction probability to thermal rate constant

First quasi-classical trajectory calculations have been carried out for the S((3)P) + OH(X?(2)Π) → SO(X?(3)Σ(-)) + H((2)S) reaction on an ab initio global potential energy surface for the ground electronic state, X?(2)A', of HSO. Cross sections, computed for collision energies up to 1 eV, show no energy threshold and decrease with the increasing collision energy. Rate constants have been calculated in the 5-500 K temperature range. The thermal rate constant is in good agreement with approximate quantum results, while a disagreement is found at 298 K with the experimental data. Product energy distributions have also been reported at four collision energies from 0.001 to 0.5 eV. The shapes of the rovibrational and angular distributions suggest the formation of an intermediate complex that is more and more long-lived as the collision energy increases.     

Constraining the mechanism and kinetics of OH + NO2 and HO2 + NO using the multiple-well master equation

Several recent experimental studies have provided substantial new constraints for the mechanisms on the HNO3 potential energy surface. These include observations of biexponential OH decay over short time scales from OH + NO2, which constrain key properties of the short-lived HOONO intermediate, observations of both conformers of the HOONO intermediate itself, isotopic scrambling data for 18OH + NO2, and observations of HONO2 production from the HO2 + NO reaction. We combine all of these recent data in a master-equation simulation of the system. This simulation is initialized with computational values for both stable species (wells) and transition states, but parameters are then adjusted to fit the observations. All parameters are kept within limits defined by experimental and theoretical uncertainty, and all converge away from their bounds. The primary fitting is carried out on the OH kinetic data-we first fit the biexponential kinetics, then address the isotopic scrambling. Isotopic scrambling is shown to be rapid but not complete at low pressure, while at least two parameter sets are shown to be consistent with the biexponential data. Of these two parameter sets, one is far more consistent with recent observations of trans-HOONO decay, isotopic scrambling, and HONO2 production from HO2 + NO. This we regard as the most probable potential energy surface for the reaction. On this PES, cis-trans isomerization for HOONO is slow but isomerization of trans-HOONO to HONO2 is rapid. This has significant implications for observed HOONO behavior and also HONO2 formation in the atmosphere from both HO2 + NO and OH + NO2.     

Collision energy dependence of the O(1D) + HCl --> OH + Cl(2P) reaction studied by crossed beam scattering and quasiclassical trajectory calculations on ab initio potential energy surfaces

The dynamics of the O(1D) + HCl --> OH + Cl(2P) reaction are investigated by a crossed molecular beam ion-imaging method and quasiclassical trajectory calculations on the three ab initio potential energy surfaces, the ground 1(1)A' and two excited (1(1)A' and 2(1)A') states. The scattering experiment was carried out at collision energies of 4.2, 4.5, and 6.4 kcal/mol. The observed doubly differential cross sections (DCSs) for the Cl(2P) product exhibit almost no collision energy dependence over this inspected energy range. The nearly forward-backward symmetric DCS indicates that the reaction proceeds predominantly on the ground-state potential energy surface at these energies. Variation of the forward-backward asymmetry with collision energy is interpreted using an osculating complex model. Although the potential energy surfaces obtained by CASSCF-MRCI ab initio calculations exhibit relatively low potential barriers of 1.6 and 6.5 kcal/mol for 1(1)A' and 2(1)A', respectively, the dynamics calculations indicate that contributions of these excited states are small at the collision energies lower than 15.0 kcal/mol. Theoretical DCSs calculated for the ground-state reaction pathway agree well with the observed ones. These experimental and theoretical results suggest that the titled reaction at collision energies less than 6.5 kcal/mol is predominantly via the ground electronic state.     

Cross sections and rate constants for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

First quasiclassical trajectory calculations have been carried out for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a recent ab initio potential energy surface for the ground electronic state, X 2A', of HCO/COH. Total and state-specific integral cross sections have been determined for a wide range of collision energies (0.001-1 eV). Then, thermal and state-specific rate constants have been calculated in the 1-500 K temperature range. The thermal rate constant varies from 1.78x10(-10) cm3 s-1 at 1 K down to 5.96x10(-11) cm3 s-1 at 500 K with a maximum value of 3.39x10(-10) cm3 s-1 obtained at 7 K. Cross sections and rate constants are found to be almost independent of the rovibrational state of OH.     

Differential cross sections and product energy distributions for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

Quasiclassical trajectory calculations have been carried out for the C((3)P)+OH(X (2)Pi)-->CO(X (1)Sigma(+))+H((2)S) reaction using a recent ab initio potential energy surface for the ground electronic state X (2)A(') of COH. Differential cross sections (DCSs), and product vibrational, rotational and translational distributions have been determined for a wide range of collision energies (0.001-1 eV). The role of excitations (rotation or vibration) of the OH reactant on these quantities has been investigated. Product vibrational, rotational, and translational distributions are found to be almost independent on the rovibrational state of OH, whereas DCSs show a weak dependence on the initial rotational state of OH. We also analyze the results using a study based on the lifetime of the intermediate complex and on the kinematic constraint associated with the mass combination.     

Quasiclassical trajectory simulations of OH(v) + NO2 --> HONO2* --> OH(v') + NO2: capture and vibrational deactivation rate constants

Quasiclassical trajectory calculations are used to investigate the dynamics of the OH(v) + NO(2) --> HONO(2) --> OH(v') + NO(2) recombination/dissociation reaction on an analytic potential energy surface (PES) that gives good agreement with the known structure and vibrational frequencies of nitric acid. The calculated recombination rate constants depend only weakly on temperature and on the initial vibrational energy level of OH(v). The magnitude of the recombination rate constant is sensitive to the potential function describing the newly formed bond and to the switching functions in the PES that attenuate inter-mode interactions at long range. The lifetime of the nascent excited HONO(2) depends strongly not only on its internal energy but also on the identity of the initial state, in disagreement with statistical theory. This disagreement is probably due to the effects of slow intramolecular vibrational energy redistribution (IVR) from the initially excited OH stretching mode. The vibrational energy distribution of product OH(v') radicals is different from statistical distributions, a result consistent with the effects of slow IVR. Nonetheless, the trajectory results predict that vibrational deactivation of OH(v) via the HONO(2) transient complex is approximately 90% efficient, almost independent of initial OH(v) vibrational level, in qualitative agreement with recent experiments. Tests are also carried out using the HONO(2) PES, but assuming the weaker O-O bond strength found in HOONO (peroxynitrous acid). In this case, the predicted vibrational deactivation efficiencies are significantly lower and depend strongly on the initial vibrational state of OH(v), in disagreement with experiments. This disagreement suggests that the actual HOONO PES may contain more inter-mode coupling than found in the present model PES, which is based on HONO(2). For nitric acid, the measured vibrational deactivation rate constant is a useful proxy for the recombination rate, but IVR randomization of energy is not complete, suggesting that the efficacy of the proxy method must be evaluated on a case-by-case basis.     

Theoretical Study on Mechanism of Reaction of OH with HO2NO2

The reaction of HO2NO2 (peroxynitric acid, PNA) with OH was studied by the hybrid density functional B3LYP and CBS-QB3 methods. Based on the calculated potential energy surface, five reaction channels, H2O+NO2+O2, HOOH+NO3, NO2+HO3H, HO2+HONO2 and HO2+HOONO, were examined in detail. The major reaction channel is PNA+OH→M1→TS1→H2O+NO2+O2. Taking a pre-equilibrium approximation and using the CBS-QB3 energies, the theoretical rate constant of this channel was calculated as 1.13×10-12 cm3/(molecule s) at 300 K, in agreement with the experimental result. Comparison between reactions of HOONO2+OH and HONO2+OH was carried out. For HOR+OH reactions, the total rate constants increase from R=NO2 to R=ONO2, which is consistent with experimental measurements.     

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.     

Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

气相硝酸及过氧亚硝酸的结构和光电子能谱

采用了B3LYP方法和6-311G(d) 基组，并辅以MP2方法详细地计算了硝酸(HONO2)和过氧亚硝酸(HOONO)及相应正离子的各种可能构象、能量、振动频率等，分析了它们的相对稳定性、电离势、光电子能谱，探讨了用光电子能谱去探测HOONO的可能性.结构优化结果表明， HOONO和HOONO＋均有三种稳定的构象， HOONO有三个一阶鞍点的结构.作出了HOONO和HOONO＋的能量随O－O键旋转的变化曲线图，定性地讨论了可能存在的物种以及各构象之间几何结构的变化.     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.     

Quantum dynamics of inelastic excitations and charge transfer processes in the H+ +NO collisions

State-resolved differential cross section, integral cross section, average vibrational energy transfer, and the relative transition probability are computed for the H(+)+NO system using our newly obtained ab initio potential energy surfaces (PES) at the multireference configuration interaction level of accuracy employing the correlation consistent polarized valence triple zeta basis set. The quantum dynamics is treated within the vibrational close-coupling rotational infinite-order sudden approximation using the coupled ground state and first excited state ab initio quasidiabatic PES. The computed collision attributes for the inelastic vibrational excitation are compared with the state-to-state scattering data available at E(c.m.)=9.5 eV and E(c.m.)=29.03 eV and are found to be in overall good agreement with those of the experiments. The results for the vibrational charge transfer processes at these collision energies are also presented.     

Cis-cis and trans-perp HOONO: action spectroscopy and isomerization kinetics

The weakly bound HOONO product of the OH+NO2+M reaction is studied using the vibrational predissociation that follows excitation of the first OH overtone (2nu1). We observe formation of both cis-cis and trans-perp conformers of HOONO. The trans-perp HOONO 2nu1 band is observed under thermal (223-238 K) conditions at 6971 cm(-1). We assign the previously published (warmer temperature) HOONO spectrum to the 2nu1 band at 6365 cm(-1) and 2nu1-containing combination bands of the cis-cis conformer of HOONO. The band shape of the trans-perp HOONO spectrum is in excellent agreement with the predicted rotational contour based on previous experimental and theoretical results, but the apparent origin of the cis-cis HOONO spectrum at 6365 cm(-1) is featureless and significantly broader, suggesting more rapid intramolecular vibrational redistribution or predissociation in the latter isomer. The thermally less stable trans-perp HOONO isomerizes rapidly to cis-cis HOONO with an experimentally determined lifetime of 39 ms at 233 K at 13 hPa (in a buffer gas of predominantly Ar). The temperature dependence of the trans-perp HOONO lifetime in the range 223-238 K yields an isomerization barrier of 33+/-12 kJ/mol. New ab initio calculations of the structure and vibrational mode frequencies of the transition state perp-perp HOONO are performed using the coupled cluster singles and doubles with perturbative triples [CCSD(T)] model, using a correlation consistent polarized triple zeta basis set (cc-pVTZ). The energetics of cis-cis, trans-perp, and perp-perp HOONO are also calculated at this level [CCSD(T)/cc-pVTZ] and with a quadruple zeta basis set using the structure determined at the triple zeta basis set [CCSD(T)/cc-pVQZ//CCSD(T)/cc-pVTZ]. These calculations predict that the anti form of perp-perp HOONO has an energy of DeltaE0=42.4 kJ/mol above trans-perp HOONO, corresponding to an activation enthalpy of DeltaH298 (double dagger 0)=41.1 kJ/mol. These results are in good agreement with statistical simulations based on a model developed by Golden, Barker, and Lohr. The simulated isomerization rates match the observed decay rates when modeled with a trans-perp to cis-cis HOONO isomerization barrier of 40.8 kJ/mol and a strong collision model. The quantum yield of cis-cis HOONO dissociation to OH and NO2 is also calculated as a function of photon excitation energy in the range 3500-7500 cm(-1), assuming D0=83 kJ/mol. The quantum yield is predicted to vary from 0.15 to 1 over the observed spectrum at 298 K, leading to band intensities in the action spectrum that are highly temperature dependent; however, the observed relative band strengths in the cis-cis HOONO spectrum do not change substantially with temperature over the range 193-273 K. Semiempirical calculations of the oscillator strengths for 2nu1(cis-cis HOONO) and 2nu1(trans-perp HOONO) are performed using (1) a one-dimensional anharmonic model and (2) a Morse oscillator model for the OH stretch, and ab initio dipole moment functions calculated using Becke, Lee, Yang, and Parr density functional theory (B3LYP), M?ller-Plesset pertubation theory truncated at the second and third order (MP2 and MP3), and quadratic configuration interaction theory using single and double excitations (QCISD). The QCISD level calculated ratio of 2nu1 oscillator strengths of trans-perp to cis-cis HOONO is 3.7:1. The observed intensities indicate that the concentration of trans-perp HOONO early in the OH+NO2 reaction is significantly greater than predicted by a Boltzmann distribution, consistent with statistical predictions of high initial yields of trans-perp HOONO from the OH+NO2+M reaction. In the atmosphere, trans-perp HOONO will isomerize nearly instantaneously to cis-cis HOONO. Loss of HOONO via photodissociation in the near-IR limits the lifetime of cis-cis HOONO during daylight to less than 45 h, other loss mechanisms will reduce the lifetime further.     

An ab initio spin-orbit-corrected potential energy surface and dynamics for the F + CH4 and F + CHD3 reactions

We report an analytical ab initio three degrees of freedom (3D) spin-orbit-correction surface for the entrance channel of the F + methane reaction obtained by fitting the differences between the spin-orbit (SO) and non-relativistic electronic ground state energies computed at the MRCI+Q/aug-cc-pVTZ level of theory. The 3D model surface is given in terms of the distance, R(C-F), and relative orientation, Euler angles ? and θ, of the reactants treating CH(4) as a rigid rotor. The full-dimensional (12D) "hybrid" SO-corrected potential energy surface (PES) is obtained from the 3D SO-correction surface and a 12D non-SO PES. The SO interaction has a significant effect in the entrance-channel van der Waals region, whereas the effect on the energy at the early saddle point is only ～5% of that at the reactant asymptote; thus, the SO correction increases the barrier height by ～122 cm(-1). The 12D quasiclassical trajectory calculations for the F + CH(4) and F + CHD(3) reactions show that the SO effects decrease the cross sections by a factor of 2-4 at low collision energies and the effects are less significant as the collision energy increases. The inclusion of the SO correction in the PES does not change the product state distributions.     

Evidence for excited spin-orbit state reaction dynamics in F+H2: theory and experiment

We describe fully quantum, time-independent scattering calculations of the F+H2-->HF+H reaction, concentrating on the HF product rotational distributions in v'=3. The calculations involved two new sets of ab initio potential energy surfaces, based on large basis set, multireference configuration-interaction calculations, which are further scaled to reproduce the experimental exoergicity of the reaction. In addition, the spin-orbit, Coriolis, and electrostatic couplings between the three quasidiabatic F+H2 electronic states are included. The calculated integral cross sections are compared with the results of molecular beam experiments. At low collision energies, a significant fraction of the reaction is due to Born-Oppenheimer forbidden, but energetically allowed reaction of F in its excited (2P 1/2) spin-orbit state. As the collision energy increases, the Born-Oppenheimer allowed reaction of F in its ground (2P 3/2) spin-orbit state rapidly dominates. Overall, the calculations agree reasonably well with the experiment, although there remains some disagreement with respect to the degree of rotational excitation of the HF(v'=3) products as well as with the energy dependence of the reactive cross sections at the lowest collision energies.     

Direct versus resonances mediated F+OH collisions on a new 3A" potential energy surface

A theoretical study of the F(2P) + OH(2Pi) --> HF(1Sigma+) + O(3P) reactive collisions is carried out on a new global potential energy surface (PES) of the ground 3A" adiabatic electronic state. The ab initio calculations are based on multireference configuration interaction calculations, using the aug-cc-pVTZ extended basis sets of Dunning et al. A functional representation of the PES shows no nominal barrier to reaction, contrary to previous results by others. Wave packet and quasiclassical trajectory calculations have been performed for this PES to study the F + OH(v = 0,j) reactive collision. The comparison was performed at fixed and constant values of the total angular momentum from 0 to 110 and relative translational energy up to 0.8 eV. The reaction presents a dynamical barrier, essentially due to the zero-point energy for the bending vibration near the saddle point. This determines two different reaction mechanisms. At energies higher than approximately 0.125 eV the reaction is direct, while below that value it is indirect and mediated by heavy-light-heavy resonances. Such resonances, also found in the simulations of the photodetachment spectrum of the triatomic anion, manifest themselves in the quasiclassical simulations, too, where they are associated to periodic orbits.     

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     

Theoretical study of the dynamics of Cl + O3 reaction I. Ab initio potential energy surface and quasiclassical trajectory results

We present a global full dimensional potential energy surface (PES) for the Cl + O(3)→ ClO + O(2) reaction, which is an elementary step in a catalytic cycle that leads to the destruction of ozone in the stratosphere. The PES is constructed by interpolation of quantum chemistry data using the method developed by Collins and co-workers. Ab initio data points (energy, gradients and Hessian matrix elements) have been calculated at the UQCISD/aug-cc-pVDZ (unrestricted quadratic configuration interaction with single and double excitations) level of theory. The ab initio calculations predict a markedly non-coplanar (dihedral angle of 80°) transition state for the reaction, located very early in the reactant valley and slightly below the energy of the reactants as long as the spin-orbit splitting is neglected. Quasiclassical trajectory (QCT) calculations have been carried out at several collision energies to investigate the reaction dynamics. The QCT excitation function shows no threshold, displays a minimum at a collision energy of 2.5 kcal mol(-1), and then increases monotonically at larger collision energies. This behaviour is consistent with a barrierless reaction dominated by an oxygen-abstraction mechanism. The calculated product vibrational distributions (strongly inverted for ClO) and rate constants are compared with experimental determinations. Differential cross sections (DCS) summed over all final states are found to be in fairly good agreement with those derived from crossed molecular beam experiments.