F+OH reactive collisions on new excited 3A" and 3A' potential-energy surfaces

Global three-dimensional adiabatic potential-energy surfaces for the excited 2(3)A" and 1(3)A' triplet states of OHF are obtained to study the F(2P)+OH(2pi)-->O(3P)+HF(1sigma+) reaction. Highly accurate ab initio calculations are obtained for the two excited electronic states and fitted to analytical functions with small deviations. The reaction dynamics is studied using a wave-packet treatment within a centrifugal sudden approach, which is justified by the linear transition state of the two electronic states studied. The reaction efficiency presents a marked preference for perpendicular orientation of the initial relative velocity vector and the angular momentum of the OH reagent, consistent in the body-fixed frame used with an initial collinear geometry which facilitates the access to the transition state. It is also found that the reaction cross section presents a rather high threshold so that, in an adiabatic picture, the two excited triplet states do not contribute to the rate constant at room temperature. Thus, only the lowest triplet state leads to reaction under these conditions and the simulated rate constants are too low as compared with the experimental ones. Such disagreement is likely to be due to nonadiabatic transitions occurring at the conical intersections near the transition state for this reaction.     

A quasiclassical trajectory study for the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) reaction on the new 2A' and 4A' potential-energy surfaces

A quasiclassical trajectory study with the sixth-order explicit symplectic algorithm of the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) atmospheric reaction has been performed by employing the new 2A' and 4A' potential-energy surfaces reported by Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. For the translational temperature considered up to 10,000 K, the larger relative translational energy and the higher rovibrational levels of O2 molecule with respect to the previous works have been taken into account, and a clearer database about the character of the total reaction cross section has been presented in this work. The dependence of microscopic rate constants on the vibrational level of O2 molecule at T=3000, 5000, and 10,000 K has been exhibited, and we can see that the values of log10 k(T,v,J) vary almost linearly with the vibrational level of O2 molecule. The thermal rate constants at the translational temperature between 300 and 10,000 K have been evaluated and compared with the experimental and previous theoretical results. It is found that the thermal rate constants determined in this work have a better agreement with the experimental data and can provide a more valid theoretical reference at the translational temperature considered for the title reaction.     

Ab initio potential-energy surface for the reaction Ca+HCl-->CaCl+H

The potential-energy surface of the ground electronic state of CaHCl has been obtained from 6400 ab initio points calculated at the multireference configuration-interaction level and represented by a global analytical fit. The Ca+HCl-->CaCl+H reaction is endothermic by 5100 cm(-1) with a barrier of 4470 cm(-1) at bent geometry, taking the zero energy in the Ca+HCl asymptote. On both sides of this barrier are potential wells at linear geometries, a shallow one due to van der Waals interactions in the entrance channel, and a deep one attributed to the H(-)Ca(++)Cl(-) ionic configuration. The accuracy of the van der Waals well depth, approximately 200 cm(-1), was checked by means of additional calculations at the coupled-cluster singles and doubles with perturbative triples level and it was concluded that previous empirical estimates are unrealistic. Also, the electric dipole function was calculated, analytically fitted in the regions of the two wells, and used to analyze the charge shifts along the reaction path. In the insertion well, 16,800 cm(-1) deep, the electric dipole function confirmed the ionic structure of the HCaCl complex and served to estimate effective atomic charges. Finally, bound rovibrational levels were computed both in the van der Waals well and in the insertion well, and the infrared-absorption spectrum of the insertion complex was simulated in order to facilitate its detection.     

Molecular potential-energy surfaces for chemical reaction dynamics

This paper reviews the construction of molecular potential-energy surfaces by an interpolation method which has been developed over the last several years. The method uses ab initio quantum chemistry calculations of the molecular electronic energy in an automated procedure to construct global potential- energy surfaces which can be used to simulate chemical reactions with either classical or quantum dynamics. The methodology is explained and several applications are presented to illustrate the approach. Received: 22 February 2002 / Accepted: 2 May 2002 / Published online: 6 November 2002 Correspondence to: M. A. Collins e-mail: collins@rsc.anu.edu.au Acknowledgements. The methods described in this overview are the result of collaborations with former members of my group, in particular with Josef Ischtwan, Meredith Jordon, Keiran Thompson and Ryan Bettens. I am also indebted for inspiration gained from many discussions with my colleagues Leo Radom and Donghui Zhang (National University of Singapore). This work has been supported by the Supercomputer Facility of the Australian National University and the Australian Partnership for Advanced Computing.     

Combined valence bond-molecular mechanics potential-energy surface and direct dynamics study of rate constants and kinetic isotope effects for the H + C2H6 reaction

This article presents a multifaceted study of the reaction H+C(2)H(6)-->H(2)+C(2)H(5) and three of its deuterium-substituted isotopologs. First we present high-level electronic structure calculations by the W1, G3SX, MCG3-MPWB, CBS-APNO, and MC-QCISD/3 methods that lead to a best estimate of the barrier height of 11.8+/-0.5 kcal/mol. Then we obtain a specific reaction parameter for the MPW density functional in order that it reproduces the best estimate of the barrier height; this yields the MPW54 functional. The MPW54 functional, as well as the MPW60 functional that was previously parametrized for the H+CH(4) reaction, is used with canonical variational theory with small-curvature tunneling to calculate the rate constants for all four ethane reactions from 200 to 2000 K. The final MPW54 calculations are based on curvilinear-coordinate generalized-normal-mode analysis along the reaction path, and they include scaled frequencies and an anharmonic C-C bond torsion. They agree with experiment within 31% for 467-826 K except for a 38% deviation at 748 K; the results for the isotopologs are predictions since these rate constants have never been measured. The kinetic isotope effects (KIEs) are analyzed to reveal the contributions from subsets of vibrational partition functions and from tunneling, which conspire to yield a nonmonotonic temperature dependence for one of the KIEs. The stationary points and reaction-path potential of the MPW54 potential-energy surface are then used to parametrize a new kind of analytical potential-energy surface that combines a semiempirical valence bond formalism for the reactive part of the molecule with a standard molecular mechanics force field for the rest; this may be considered to be either an extension of molecular mechanics to treat a reactive potential-energy surface or a new kind of combined quantum-mechanical/molecular mechanical (QM/MM) method in which the QM part is semiempirical valence bond theory; that is, the new potential-energy surface is a combined valence bond molecular mechanics (CVBMM) surface. Rate constants calculated with the CVBMM surface agree with the MPW54 rate constants within 12% for 534-2000 K and within 23% for 200-491 K. The full CVBMM potential-energy surface is now available for use in variety of dynamics calculations, and it provides a prototype for developing CVBMM potential-energy surfaces for other reactions.     

Mixed quantum-classical reaction path dynamics of C2H5F --> C2H4 + HF

A mixed quantum-classical method for calculating product energy partitioning based on a reaction path Hamiltonian is presented and applied to HF elimination from fluoroethane. The goal is to describe the effect of the potential energy release on the product energies using a simple model of quantized transverse vibrational modes coupled to a classical reaction path via the path curvature. Calculations of the minimum energy path were done at the B3LYP/6-311++G(2d,2p) and MP2/6-311++G** levels of theory, followed by energy-partitioning dynamics calculations. The results for the final HF vibrational state distribution were found to be in good qualitative agreement with both experimental studies and quasiclassical trajectory simulations.     

Conical intersection between the lowest spin-aligned Li3(4A') potential-energy surfaces

We have calculated new potential-energy surfaces for the lowest two spin-aligned (4)A(') states of the Li(3) trimer. This calculation shows a seam of conical intersections between these states resulting from the extra symmetry of the system when the atoms are in a collinear arrangement. This seam is especially important because of its proximity to the three-body dissociation limit of the system; ultracold scattering calculations and the bound-state energies of the system will be affected by the presence of this conical intersection. In this paper we discuss the calculation of the potential-energy surface and the location of the conical intersection seam.     

Mechanism and kinetics of the reaction NO3 + C2H4

The reaction of NO(3) radical with C(2)H(4) was characterized using the B3LYP, MP2, B97-1, CCSD(T), and CBS-QB3 methods in combination with various basis sets, followed by statistical kinetic analyses and direct dynamics trajectory calculations to predict product distributions and thermal rate constants. The results show that the first step of the reaction is electrophilic addition of an O atom from NO(3) to an olefinic C atom from C(2)H(4) to form an open-chain adduct. A concerted addition reaction mechanism forming a five-membered ring intermediate was investigated, but is not supported by the highly accurate CCSD(T) level of theory. Master-equation calculations for tropospheric conditions predict that the collisionally stabilized NO(3)-C(2)H(4) free-radical adduct constitutes 80-90% of the reaction yield and the remaining products consist mostly of NO(2) and oxirane; the other products are produced in very minor yields. By empirically reducing the barrier height for the initial addition step by 1 kcal mol(-1) from that predicted at the CBS-QB3 level of theory and treating the torsional modes explicitly as one-dimensional hindered internal rotations (instead of harmonic oscillators), the computed thermal rate constants (including quantum tunneling) can be brought into very good agreement with the experimental data for the overall reaction rate constant.     

Quantum-chemical investigation of the potential-energy surface of the T+CH3NH2 reaction

Conclusions According to the data from the calculations by the semiempirical SCF method in the MINDO/3 all-valence electron approximation and by a nonempirical method in the STO-4-31G basis, the reaction of hot T atoms in the T+CH₃NH₂ system most probably takes place according to a mechanism involving the abstraction of an H atom. In the case of hot T atoms there is also a possibility of the replacement of the NH₂ group according to a Waiden inversion mechanism or its direct abstraction.Translated from Izvestiya Akademii Nauk SSSR, Seriya Khimicheskaya, No. 9, pp. 1958–1961, September, 1982.     

Following the streambed reaction on potential-energy surfaces: a new robust method

A simple procedure with low computational efforts is proposed to follow the reaction path of the potential-energy hypersurface (PES) starting from minima or saddle points. The method uses a modification of the so-called “following the reduced gradient” [Quapp W, Hirsch M, Imig O, Heidrich D (1998) J Comput Chem 19:1087]. The original method connects points where the gradient has a constant direction. In the present article the procedure is replaced by taking iterative varying directions of the gradient controlled by the last tangent of the searched curve. The resulting minimum energy path is that valley floor gradient extremal (GE) which belongs to the smallest (absolute) eigenvalue of the Hessian and, hence, that GE which usually leads along the streambed of a chemical reaction. The new method avoids third derivatives of the PES and obtains the GE of least ascent by second-order calculations only. Nevertheless, we are able to follow the streambed GE uphill or downhill. We can connect a minimum with its saddles if the streambed leads up to a saddle, or we find a turning point or a bifurcation point. The effectiveness and the characteristic properties of the new algorithm are demonstrated by using polynomial test surfaces, an ab initio PES of H₂O, and the analytic potentials of Lennard-Jones (LJ) clusters. By tracing the streambeds we located previously identified saddle points for LJ _N with N=3, 7, 8, and 55. Saddles for LJ _N with N=15, 20, and 30 as presented here are new results. Received: 8 March 2000 / Accepted: 17 July 2000 / Published online: 24 October 2000     

An eight-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction on a new modified potential energy surface

An eight-dimensional time-dependent quantum dynamics wave packet approach is performed for the study of the H2+C2H-->H+C2H2 reaction system on a new modified potential energy surface (PES) [L.-P. Ju et al., Chem. Phys. Lett. 409, 249 (2005)]. This new potential energy surface is obtained by modifying Wang and Bowman's old PES [J. Chem. Phys. 101, 8646 (1994)] based on the new ab initio calculation. This new modified PES has a much lower transition state barrier height at 2.29 kcal/mol than Wang and Bowman's old PES at 4.3 kcal/mol. This study shows that the reactivity for this diatom-triatom reaction system is enhanced by vibrational excitations of H2, whereas the vibrational excitations of C2H only have a small effect on the reactivity. Furthermore, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. The comparison of the rate constant between this calculation and experimental results agrees with each other very well. This comparison indicates that the new modified PES corrects the large barrier height problem in Wang and Bowman's old PES.     

Modeling the H5+ potential-energy surface: a first attempt

 Ab initio molecular electronic structure calculations are performed for H₅ ⁺ at the QCISD(T) level of theory, using a correlation-consistent quadruple-zeta basis set. Structures, vibrational frequencies and thermochemical properties are evaluated for ten stationary points of the H₅ ⁺ hypersurface and are compared with previous calculations. The features of the H₃ ⁺…H₂ interaction at intermediate and large intermolecular distances are also investigated. Furthermore, an analytical functional form for the potential-energy surface of H₅ ⁺ is derived using a first-order diatomics-in-molecule perturbation theory approach. Its topology is found to be qualitatively correct for the short-range interaction region. Received: 15 March 2001 / Accepted: 5 July 2001 / Published online: 11 October 2001     

An adiabatic capture theory and quasiclassical trajectory study of C + NO and O + CN on the 2A', 2A", and 4A" potential energy surfaces

The adiabatic capture centrifugal sudden approximation (ACCSA) has been applied to the C + NO and O + CN reactions, along with quasiclassical trajectory simulations. Existing global analytic fits to the potential energy surfaces of the CNO system in the (2)A', (2)A", and (4)A" electronic states have been used. Thermal rate constants for reaction in each of the electronic states have been calculated. In all cases a strong temperature dependence is evident in the calculated rate constants. The agreement between the calculated adiabatic capture and quasiclassical trajectory rate constants is excellent in some cases, but these rate constants differ considerably in other cases. This behavior is analyzed in terms of the anisotropy of the potential energy surfaces. On the basis of this analysis, we propose a new diagnostic for the reliability of ACCSA capture calculations.     

An accurate study of the dynamics of the C+OH reaction on the second excited 1(4)A" potential energy surface

The dynamics of the C((3)P)+OH(X(2)Π) → CO(a(3)Π)+H((2)S) on its second excited potential energy surface, 1(4)A", have been investigated in detail by means of an accurate quantum mechanical (QM) time-dependent wave packet (TDWP) approach. Reaction probabilities for values of the total angular momentum J up to 50 are calculated and integral cross sections for a collision energy range which extends up to 0.1 eV are shown. The comparison with quasi-classical trajectory (QCT) and statistical methods reveals the important role played by the double well structure existing in the potential energy surface. The TDWP differential cross sections exhibit a forward-backward symmetry which could be interpreted as indicative of a complex-forming mechanism governing the dynamics of the process. The QM statistical method employed in this study, however, is not capable to reproduce the main features of the possible insertion nature in the reactive collision. The ability to stop individual trajectories selectively at specific locations inside the potential energy surface makes the QCT version of the statistical approach a better option to understand the overall dynamics of the process.     

Potential energy surface for the CCl4 + H --> CCl3 + ClH reaction: kinetics and dynamics study

An analytical potential energy surface for the gas-phase CCl4 + H --> CCl3 + ClH reaction was constructed with suitable functional forms to represent vibrational modes. This surface is completely symmetric with respect to the permutation of the four chlorine atoms and is calibrated with respect to experimental thermal rate constants available over the temperature range 297-904 K. On this surface, the thermal rate constants were calculated using variational transition-state theory with semiclassical transmission coefficients over a wider temperature range 300-2500 K, therefore obtaining kinetics information at higher temperatures than are experimentally available. This surface was also used to analyze dynamical features, such as tunneling and reaction-path curvature. In the first case, the influence of the tunneling factor is very small since a heavy chlorine atom has to pass through the barrier. In the second, it was found that vibrational excitation of the Cl-H stretching mode can be expected in the exit channel.     

State-to-state quantum dynamics calculations of the C + OH reaction on the second excited potential energy surface

Accurate three-dimensional quantum-mechanical scattering calculations using a time-indepedent hyperspherical method have been performed for the C((3)P) + OH(X(2)Π) → CO(a(3)Π) + H((2)S) reaction on the second excited potential energy surface of 1(4)A″ symmetry. State-to-state reaction probabilities at a total angular momentum J = 0 have been computed in a wide range of collision energies. Many pronounced resonances have been found, espcially at low energy. The product vibrational distributions are noninverted. The present results therefore suggest that the title reaction proceeds via a long-lived intermediate complex. An approximate quantum-mechanical rate constant has also been calculated, and large differences are observed with the quasi-classical trajectory prediction.     

Radical reaction C3H+NO: a mechanistic study

Although a number of hydrocarbon radicals including the heavier C(3)-radicals C(3)H(3) and C(3)H(5) have been experimentally shown to deplete NO effectively, no theoretical or experimental attempts have been made on the reactivity of the simplest C(3)-radical towards NO. In this article, we report our detailed mechanistic study on the C(3)H+NO reaction at the Gussian-3//B3LYP/6-31G(d) level by constructing the singlet and triplet electronic state [H,C(3),N,O] potential energy surfaces (PESs). The l-C(3)H+NO reaction is shown to barrierlessly form the entrance isomer HCCCNO followed by the direct O-elimination leading to HCCCN+(3)O on triplet PES, or by successive O-transfer, N-insertion, and CN bond-rupture to generate the product (1)HCCN+CO on singlet PES. The possible singlet-triplet intersystem crossings are also discussed. Thus, the novel reaction l-C(3)H+NO can proceed effectively even at low temperatures and is expected to play an important role in both combustion and interstellar processes. For the c-C(3)H+NO reaction, the initially formed H-cCCC-NO can most favorably isomerize to HCCCNO, and further evolution follows that of the l-C(3)H+NO reaction. Quantitatively, the c-C(3)H+NO reaction can take place barrierlessly on singlet PES, yet it faces a small barrier 2.7 kcal/mol on triplet PES. The results will enrich our understanding of the chemistry of the simplest C(3)-radical in both combustion and interstellar processes, which to date have received little attention despite their importance and available abundant studies on its structural and spectroscopic properties.     

Non-RRKM dynamics in the CH3O2 + NO reaction system

Quasi-classical trajectory (QCT) calculations on a model potential energy surface (PES) show strong deviations from statistical Rice-Ramsperger-Kassel-Marcus (RRKM) rate theory for the decomposition reaction (1) CH3OONO* --> CH3O + NO2, where the highly excited CH3OONO* was formed by (2) CH3O2 + NO --> CH3OONO*. The model PES accurately describes the vibrational frequencies, structures, and thermochemistry of the cis- and trans-CH3OONO isomers; it includes cis-trans isomerization in addition to reactions 1 and 2 but does not include nitrate formation, which is too slow to affect the decay rate of CH3OONO*. The QCT results give a strongly time-dependent rate constant for decomposition and damped oscillations in the decomposition rate, not predicted by statistical rate theory. Anharmonicity is shown to play an important role in reducing the rate constant by a factor of 10 smaller than predicted using classical harmonic RRKM theory (microcanonical variational transition state theory). Master equation simulations of organic nitrate yields published previously by two groups assumed that RRKM theory is accurate for reactions 1 and 2 but required surprising parametrizations to fit experimental nitrate yield data. In the present work, it is hypothesized that the non-RRKM rate of reaction (1) and vibrational anharmonicity are at least partly responsible for the surprising parameters.     

Theoretical study of the potential-energy surface of C2NP

 The B3LYP/6-311G(d) and CCSD(T)/6-311G(2df) (single-point) methods have been used to investigate the singlet potential energy surface of C₂NP, in which seven stationary isomers and seventeen interconversion transition states are involved. At the final CCSD(T)/6-311G(2df)//B3LYP6-311G(d) level with zero-point vibrational energy correction the lowest-lying isomer is a linear NCCP followed by two linear CNCP isomers at 23.9  and CCNP at 65.8 kcal mol⁻¹, respectively. The three isomers are kinetically very stable towards both isomerization and dissociation, and CCNP is even more kinetically stable than CNCP – by 14.3 kcal mol⁻¹ despite its high energy. Further comparative calculations were performed at the QCISD and QCISD(T) levels with the 6-311G(d) and 6-311G(2d) basis sets to obtain more reliable structures and spectroscopy for the three isomers. The calculated bond lengths, rotational constant, and dipole moment for NCCP were in close agreement with the experimentally determined values. Finally, similarities and discrepancies between the potential energy surface of C₂NP and those of the analogous species C₂N₂ and C₂P₂ were compared. The results presented in this paper might be helpful for future identification of the two still unknown yet kinetically very stable isomers CNCP and CCNP, both in the laboratory and in interstellar space. Received: 3 January 2001 / Accepted: 6 June 2001 / Published online: 30 October 2001