CBr+O2反应在二重态势能面上的机理研究

用UB3LYP/6-311++G(d,p)和QCISD(单点能)的方法考察了CBr+O2反应在二重态势能面上的反应机理。研究发现该反应在高温过程中重要,且有两个产物通道,它们分别是BrCO+O和Br+CO2,其中前者为优势通道。为了弄清溴原子取代对次甲基与氧气反应的机理的影响,我们对CBr+O2反应与CH+O2反应的相似性和差异也作了讨论。结果表明:两反应的第一步都是CX(X=H,Br)自由基与氧气反应生成链状过氧化物XCOO,且溴原子取代对反应的活性、产物通道的数量和产物的形成过程等都有影响。     

An improved potential energy surface for the F+H2 reaction

A new ground state potential energy surface has been developed for the F+H₂ reaction. Using the UCCSD(T) method, ab initio calculations were performed for 786 geometries located mainly in the exit channel of the reaction. The new data was used to correct exit channel errors that have become apparent in the potential energy surface of Stark and Werner [J. Chem. Phys. 104 (1996) 6515]. While the entrance channel and saddlepoint properties of the Stark–Werner surface are unchanged on the new potential, the exit channel behavior is more satisfactory. The exothermicity on the new surface is much closer to the experimental value. The new surface also greatly diminishes the exit channel van der Waals well that was too pronounced on the Stark–Werner surface. Several preliminary dynamical scattering calculations were carried out using the new surface for total angular momentum equal to zero for F+H₂ and F+HD. It is found that gross features of the reaction dynamics are quite similar to those predicted by the Stark–Werner surface, in particular the reactive resonance for F+HD and F+H₂ survive. However, the most of the exit channel van der Waals resonances disappear on the new surface. It is predicted that the differential cross-sections at low collision energy for the F+H₂ reaction may be drastically modified from the predictions based on the Stark–Werner surface.     

The anharmonic potential energy surface of methyl fluoride

A large experimental spectroscopic data set sensitive to the cubic anharmonic potential energy surface (PES) of methyl fluoride has been compiled from the literature for six symmetric and asymmetric top isotopomers of methyl fluoride: 12CH3F, 13CH3FF, 12CD3F, 13CD3F, 12CHD2F and 12CH2DF. This empirical data set has been used to critically assess the best available literature ab initio cubic force field and various 'improved' theoretical force fields. A perturbation-resonance approach to the calculation of the observables from the force constants has been utilized and existing PESs were found to reproduce the data poorly. The careful treatment required for the correct theoretical reproduction of empirical spectroscopic constants arising from highly correlated least-squares fits to the original data is discussed. A new fit to the data has been performed (optimising 19 of the 38 cubic force constants) resulting in a much improved PES. The latter has been used to predict currently unknown spectroscopic constants for the least well-characterised isotopomer 12CH2DF. The prospects for a future empirical determination of the complete cubic force field of methyl fluoride are discussed and new data most likely to yield new information on the PES identified.     

Ab initio potential energy surface of CH and reaction dynamics of H + CH+

This work presents a new ground state potential energy surface (PES) for CH. The potential is tested using quasi classical trajectory (QCT) and quantum reactive scattering methods for the H + CH(+) reaction. Cross sections and rate coefficients for all reaction channels up to 300 K are calculated. The abstraction rate coefficients follow the expected slightly decreasing behaviour above 90 K, but have a positive gradient with lower temperatures. The inelastic collision and exchange reaction rate constants are increasing monotonically with temperature. The rate coefficients of the exchange reaction differ significantly between QCT and quantum reactive scattering, due to intrinsic shortcomings of the QCT final state distributions.     

Quantum mechanical study of the potential energy surface of the ClO + NO2 reaction

In the study of the reaction pathways of the ClO + NO2 reaction including reliable structures of the reactants, products, intermediates, and transition states as well as energies the MP2/6-311G(d), B3LYP/6-311G(d), and G2(MP2) methods have been employed. Chlorine nitrate, ClONO2, is formed by N-O association without an entrance barrier and is stabilized by 29.8 kcal mol(-1). It can undergo either a direct 1,3 migration of Cl or OCl rotation to yield an indistinguishable isomer. The corresponding barriers are 45.8 and 7.1 kcal mol(-1), respectively. ClONO2 can further decompose into NO3 + Cl with an endothermicity of 46.4 kcal mol(-1). The overall endothermicity of the NO2 + ClO --> NO3 + Cl reaction is calculated to be 16.6 kcal mol(-1). The formation of cis-perp and trans-perp conformer of chlorine preoxynitrite, ClOONO(cp) and ClOONO(tp), are exothermic by 5.4 and 3.8 kcal mol(-1), respectively. Calculations on the possible reaction pathways for the isomerization of ClOONO to ClONO2 showed that the activation barriers are too high to account for appreciable nitrate formation from peroxynitrite isomerization. All quoted relative energies are related to G2(MP2) calculations.     

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.     

自由基HO2与C2H2反应势能面的理论研究

应用量子化学从头计算和密度泛函理论(DFT)对HO2+C2H2反应体系的反应机理进行了研究.在B3LYP/6-311G**和CCSD(T)/6-311G**水平上计算了HO2+ C2H2反应的二重态反应势能面.计算结果表明,主要反应方式为自由基HO2的H原子和C2H2分子中的C原子结合,经过一系列异构化,最后分解得到主要产物P1 (CH2O+ HCO).此反应是放热反应,化学反应热为-321.99 kJ·mol-1.次要产物为P2 (CO2 +CH3),也是放热反应.     

Theoretical study of the Si2NO potential energy surface

The structures, spectroscopies, and stabilities of the doublet Si₂NO radical are explored at the density functional theory (DFT) and ab initio levels. Seventeen isomers are located, connected by 26 interconversion transition states. At the CCSD(T)/6‐311+G(2df)//QCISD/6‐311G(d)+ZPVE level, three low‐lying isomers are predicted, that is, one bent species SiNSiO 3 (5.1 kcal/mol) containing the important Si?N triple bonding and two four‐membered ring isomers including cyclic cSiNSiO 1 (0.0) with Si? Si cross‐bonding with C_2v symmetry and puckered cSiNSiO 1′ (11.9) with divalent carbene character. Three low‐lying isomers 1, 1′, and 3 have reasonable kinetic stabilities and might be observable either experimentally or astrophysically. The possible formation strategies of 1, 1′, and 3 in laboratory and in space are discussed in detail. The calculated vibrational frequencies and possible formation processes of 3 are consistent with recent experimental observations. In light of the fact that no cyclic nitrogen‐containing species have been detected in space, two cyclic isomers 1 and 1′ could be promising candidates. Furthermore, the bonding nature of three isomers 1, 1′, and 3 is analyzed. The calculated results are also compared with those of the analogue C₂NO radical. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

HSCCS自由基势能面的理论研究

在CCSD(T)/6-311G(d,f)//MP2/6-311G(d,f) ZPE水平下,计算得到含有8个异构体和11个过渡态的HSCCS自由基体系势能面,讨论了异构体的结构与稳定性及异构体之间的异构化过程.结果表明异构体m1的能量最低,除m1以外,异构体m2和m3的能量也比较低,在MP2水平上,过渡态TS1的能量比异构体m2仅高3.9kJ/mol,而在CCSD(T)水平上,TS1的能量比m2低14.6 kJ/mol,这说明异构体m2可以迅速转化为能量更低的m1.异构体m3的能量比异构体m1高49.99 kJ/mol,可以推断,在合适的实验条件下可以观测到异构体m1.     

Ab initio study of the potential energy surface for the OH+CO-->H+CO2 reaction

Potential energy surface for the reaction OH+CO-->H+CO2 has been calculated using the complete active space self-consistent-field and multireference configuration interaction methods with the correlation consistent triple-, quadruple-, and quintuple-zeta basis sets. A specific- reaction-parameters density functional theory has been suggested, in which the B3LYP functional is reoptimized to give the highly accurate potential energy surface with less computational efforts.     

Theoretical study on the potential energy surface and vibrational energy levels of HXeI

The potential energy surface for the electronic ground state of the HXeI molecule is constructed by using the internally contracted multi-reference configuration interaction with the Davidson correction（icMRCI＋Q）method and large basis sets. The stabilities and dissociation barriers are identified from the potential energy surfaces.The three-body dissociation channel is found to be the dominate dissociation channel for HXeI.Based on the obtained potentials,vibrational energy levels of HXeI are calculated using the Lanczos algorithm.Our theoretical results are in excellent agreement with the available observed values.     

The ground state potential energy surface of methyl fluoride dimer

Potential energy surface for methyl fluoride dimer has been studied theoretically with ab initio molecular orbital method, using a 4-31G basis set. Dimer dissociation energies, Mulliken electronic populations, and dipole moments were obtained.     

Theoretical studies of the reaction channels on the SO2/OH/NO singlet potential energy surface

Ab initio MP2/6-311++G(2d,2p) investigation of the SO2/OH/NO singlet potential energy surface (PES) has been performed with the aim to localize and describe the existing minima and transition states linking them. The systematic studies have revealed seven minima, with the trans-HONO-SO2 complex (1t) being the global minimum. Eight transition states between minima or between minima and the relevant reactant species have been described. Several available izomerization and dissociation routes have been identified and discussed. The most favorable association of HOSO2 and NO was found to be a barrierless process forming nitrososulfonic acids. Isomerizations between trans-, cis-, and gauche- nitrososulfonic acids (2t, 2c, and 2g) are possible with low-energy barriers. The HOSO2 and NO species can also react via another channels involving high-energy transition states to produce the HOSO-NO2 (3) and HNO-SO3 (4) complexes.     

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.     

Modeling the global potential energy surface of the N + N2 reaction from ab initio data

A new global potential energy surface for the N + N2 exchange reaction has been built from ab initio data. To overcome the difficulty of carrying out ab initio calculations for a large set of geometries the alternative strategy of fitting the minimum energy paths of the surface and their angular dependence using a modified LAGROBO functional form has been adopted. In this way we have been able to reproduce all the main features of the potential using a fairly small set of ab initio values.     

Theoretical study on the potential energy surface of SiC2O

The DFT/B3LYP/6-311G(d) and CCSD(T)/6-311G(2df) (single-point) methods have been used to investigate the potential energy surface of SiC₂O. A total of 11 isomers are located, which are connected by 13 interconversion transition states. The global minimum is the linear SiCCO with ¹Σ⁺ electronic state whose structure can be described as cumulenic double bonding form :SiCCÖ:. The isomer O-cSiCC possessing SiCC three-membered ring with an exocyclic OSi bond also shows considerable stability. Both SiCCO and O-cSiCC are kinetically stable towards isomerization and dissociation. Further calculations are performed at the CCSD(T)/6-311G(2df)//QCISD/6-311G(2d) level to obtain more reliable energies, geometrical parameters, harmonic vibrational frequencies and infrared intensities of the isomers. The energies and geometrical parameters are also calculated with the CASPT2//CASSCF methods.     

Exact quantum mechanical reaction probabilities for the collinear H + H2 reaction on a porter-karplus potential energy surface

Exact quantum mechanical calculation of the reaction probability for the collinear H + H₂ reaction on a Porter-Karplus potential energy surface are carried out by the finite-difference boundary value method at 6 energes in the threshold region and compared to close coupling, distorted wave, classical S matrix, transition state theory, and vibrational adiabatic calculations.     

Quasiclassical trajectory calculations of the OH+NO2 association reaction on a global potential energy surface

We report a full-dimensional potential energy surface (PES) for the OH+NO(2) reaction based on fitting more than 55,000 energies obtained with density functional theory-B3LYP6-311G(d,p) calculations. The PES is invariant with respect to permutation of like nuclei and describes all isomers of HOONO, HONO(2), and the fragments OH+NO(2) and HO(2)+NO. Detailed comparison of the structures, energies, and harmonic frequencies of various stationary points on the PES are made with previous and present high-level ab initio calculations. Two hydrogen-bond complexes are found on the PES and confirmed by new ab initio CASPT2 calculations. Quasiclassical trajectory calculations of the cross sections for ground rovibrational OH+NO(2) association reactions to form HOONO and HONO(2) are done using this PES. The cross section to form HOONO is larger than the one to form HONO(2) at low collision energies but the reverse is found at higher energies. The enhancement of the HOONO complex at low collision energies is shown to be due, in large part, to the transient formation of a H-bond complex, which decays preferentially to HOONO. The association cross sections are used to obtain rate constants for formation of HOONO and HONO(2) for the ground rovibrational states in the high-pressure limit.     

Accurate potential energy surface and quantum reaction rate calculations for the H+CH4-->H2+CH3 reaction

Calculations for the cumulative reaction probability N(E) (for J=0) and the thermal rate constant k(T) of the H+CH(4)-->H(2)+CH(3) reaction are presented. Accurate electronic structure calculations and a converged Shepard-interpolation approach are used to construct a potential energy surface which is specifically designed to allow the precise calculation of k(T) and N(E). Accurate quantum dynamics calculations employing flux correlation functions and multiconfigurational time-dependent Hartree wave packet propagation compute N(E) and k(T) based on this potential energy surface. The present work describes in detail the various convergence test performed to investigate the accuracy of the calculations at each step. These tests demonstrate the predictive power of the present calculations. In addition, approximate approaches for reaction rate calculations are discussed. A quite accurate approximation can be obtained from a potential energy surface which includes only interpolation points on the minimum energy path.     

Kinetic study on the H + SiH4 abstraction reaction using an ab initio potential energy surface

Variational transition state theory calculations with the correction of multidimensional tunneling are performed on a 12-dimensional ab initio potential energy surface for the H + SiH(4) abstraction reaction. The surface is constructed using a dual-level strategy. For the temperature range 200-1600 K, thermal rate constants are calculated and kinetic isotope effects for various isotopic species of the title reaction are investigated. The results are in very good agreement with available experimental data.