A fully ab initio potential energy surface for ClH2 reactive system

An ab initio analytical potential energy surface called BW3 for the ClH2 reactive system is presented. The fit of this surface is based on about 1 200 ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H2 and HCl, and the endothermicity of the Cl + H2 abstraction reaction. For the Cl + H2 abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H + ClH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

A fully ab initio potential energy surface for ClH2 reactive system

Anab initio analytical potential energy surface called BW3 for the ClH₂ reactive system is presented. The fit of this surface is based on about 1 200ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H₂ and HCl, and the endothermicity of the Cl + H₂ abstraction reaction. For the Cl + H₂ abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H + ClH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

A fullyab initio potential energy surface for ClH2 reactive system

Anab initio analytical potential energy surface called BW3 for the ClH₂ reactive system is presented. The fit of this surface is based on about 1 200ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H₂ and HCl, and the endothermicity of the Cl + H₂ abstraction reaction. For the Cl + H₂ abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H + ClH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

Ab initio Mechanism Study on the Reaction of Chlorine Atom with Formic Acid

The potential energy surface(PES) for the reaction of Cl atom with HCOOH is predicted using ab initio molecular orbital calculation methods at UQCIDS(T,full)6-311 G(3df,2p)//UMP2(full)/6-311 G(d,P) level of theory with zero-point vibrational energy (ZPVE) correction.The calculated results show that the reaction mechanism of Cl atom with formic acid is a C-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom with a 3.73kJ/mol reaction barrier height,leading to the formation of cis-HOCO radical which will reacts with Cl atom or other molecules in such a reaction system.Because the reaction barrier height of O-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom which leads to the formation of HCO2 radical is 67.95kJ/mol,it is a secondary reaction channel in experiment,This is in good agreement with the prediction based on the previous experiments.     

Understanding the complex dissociation dynamics of energy selected dichloroethylene ions: neutral isomerization energies and heats of formation by imaging photoelectron-photoion coincidence

The dissociative photoionization of 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2) has been investigated at high energy and mass resolution using the imaging photoelectron photoion coincidence instrument at the Swiss Light Source. The asymmetric Cl-atom loss ion time-of-flight distributions were fitted to obtain the dissociation rates in the 10(3) s(-1) < k < 10(7) s(-1) range as a function of the ion internal energy. The results, supported by ab initio calculations, show that all three ions dissociate to the same C(2v) symmetry ClC═CH(2)(+) product ion. The 0 K onset energies thus establish the relative heats of formation of the neutral isomers, that is, the isomerization energies. The experimental rate constants, k(E), as well as ab initio calculations indicate an early isomerization transition state and no overall reverse barrier to dissociation. The major high energy channels are the parallel HCl loss and the sequential ClC═CH(2)(+) → HCCH(+) + Cl process, the latter in competition with a ClC═CH(2)(+) → ClCCH(+) + H reaction. A parallel C(2)H(2)Cl(2)(+) → C(2)HCl(2)(+) + H channel also weakly asserts itself. The 0 K onset energy for the sequential Cl loss reaction suggests no barrier to the production of the most stable acetylene ion product; thus the sequential Cl-atom loss is preceded by a ClC═CH(2)(+) → HC(Cl)CH(+) reorganization step with a barrier lower than that of the second Cl-atom loss. The breakdown diagram corresponding to this sequential dissociation reveals the internal energy distribution of the first C(2)H(2)Cl(+) daughter ion, which is determined by the kinetic energy release in the first, Cl loss reaction at high excess energies. At low kinetic energy release, this distribution corresponds to the predicted two translational degrees of freedom, whereas at higher energies, the excess energy partitioning is characteristic of only one translational degree of freedom. New Δ(f)H(o)(298K) of 3.7, 2.5, and 0.2 ± 1.75 kJ mol(-1) are proposed for 1,1-C(2)H(2)Cl(2), (E)-1,2-C(2)H(2)Cl(2), and (Z)-1,2-C(2)H(2)Cl(2), respectively, and the proton affinity of ClCCH is found to be 708.6 ± 2.5 kJ mol(-1).     

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.     

Theoretical study on the reaction mechanisms of ethylene with Cp2TiR, Cp2Ti(Cl)R, and Cp2Ti(Cl:AlH2Cl)R (R=H and CH3)

The potential energy surfaces for reactions of ethylene with Cp₂Ti⁺R, Cp₂Ti(Cl)R, and Cp₂Ti(Cl:AlH₂Cl)R (R=H and CH₃) were calculated by ab initio molecular orbital methods. These six reaction mechanisms were compared. Of the two possible reaction paths, attack of ethylene at Ti and the Cl and R ligands (path IN) and that from the opposite side of the Cl ligand (path OUT), the former is found to be more favorable, with a very low activation energy for reaction of ethylene with Cp₂Ti(Cl)H. For reaction of ethylene with Cp₂Ti(Cl)CH₃, the insertion transition states on both paths have almost the same energy barrier height above the reactants. For reaction of ethylene with Cp₂Ti(Cl:AlH₂Cl)R, the bond alternation between Ti–Cl and Cl–Al plays an important role in the mechanisms.     

CX_2(X=H,F,Cl)与甲基异丙基醚C-H键插入反应的理论研究

用从头计算方法在MP2 /6 31G(d)水平上研究了CX2 (X =H ,F ,Cl)与甲基异丙基醚的C -H键插入反应。CCl2 与甲基异丙基醚两个不同的α C的C -H键插入势垒分别为 117.2kJ/mol (甲基 )和 2 0 .6kJ/mol (异丙基 )。CF2 与异丙基α C的C -H键上插入势垒为 12 0 .0kJ/mol,在插入甲基上C -H键时会引起C -O键的断裂。CH2 的插入反应则不需要势垒。对CX2 与二甲醚、甲乙醚、甲基异丙基醚、甲基苄基醚上各种不同的C -H键插入势垒进行了比较 ,甲基和苯基都促使其毗邻的C -H键更容易被CX2 所插入     

Ab initio and DFT studies on hydrolyses of phosphorus halides

Hydrolyses of phosphorus halides, (RO)(2)POX where R = H or Me and X = F or Cl, in the gas phase and in the reaction field have been investigated theoretically with ab initio and the density functional theory (DFT). The free energy of activation in the reaction field was also estimated using the Onsager method with a correction of entropy change and basis set superposition error (BSSE). The reaction of (MeO)(2)POF proceeds through a path with bifunctional catalysis regardless of the medium, but the reaction of (MeO)(2)POCl proceeds through bifunctional and general base catalysis in the gas phase and in water, respectively. The estimated free energy barrier of 23 kcal/mol for the hydrolysis of (MeO)(2)POF is in good agreement with the experimental values of 24 kcal/mol, and relative barrier of 3 kcal/mol to the (MeO)(2)POCl is also in good agreement with the experimental values of 5 kcal/mol of diisopropyl phosphorus halides ((Pr(i)O)(2)POX, X = F and Cl).     

Modeling water exchange on an aluminum polyoxocation

For the first time, water exchange on a polymeric complex has been modeled using a combination of gas-phase ab initio calculations and molecular dynamics (MD) simulations. The GaO4Al12(OH)24(H2O)12(7+)aq ion (GaAl12) was chosen because high-quality experimental data exist, including an activation enthalpy (+63 +/- 7 kJ/mol) and an activation volume (+3 +/- 1 cm3/mol). We took a two-step approach. First, the local solvent structure and the initial states for reaction were inferred from the molecular dynamics simulations. Second, we used this information to evaluate initial-state structures in the ab initio calculations. The energy differences between the initial and transition states from the ab initio calculations varied from +59 kJ/mol to +53 kJ/mol depending upon details, closely approximating the activation enthalpy.     

环丙基硅烯的理论研究: 环丙基硅烯C3H5SiH的重排反应及其机理

本文用限制的Hartree-Fock解析梯度方法在3-21G和6-31G^*水平上对环丙基硅烯的重排反应及其机理进行了从头算研究。以6-31G^*优化构型作了二级微扰计算, 并计算了各构型的频率。在此基础上得到了重排反应的热焓△H, 自由能△G和平衡常数K, 用Eyring过渡态理论计算了反应的速度常数k(T), 应用Woodward-Hoffmann规则讨论了环丙基硅烯重排反应过程中端基的旋转机理。结果表明, 环丙基硅烯经过113.4kJ/mol的势垒扩环重排为硅杂环丁烯为自发反应; 而其1,2-氢迁移重排反应热垒为190.0kJ/mol, 是非自发反应, 难于进行,不能与扩环重排相竞争。另外, 扩环重排反应可分为多步过程, 每步的端基旋转均可用Woodward-Hoffmann规则说明。     

Ab initio potential energy and dipole moment surfaces of (H2O)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with M?ller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).     

Mechanistic study of the deamination reaction of guanine: a computational study

The mechanism for the deamination of guanine with H(2)O, OH(-), H(2)O/OH(-) and for GuaH(+) with H(2)O has been investigated using ab initio calculations. Optimized geometries of the reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), B3LYP/6-31G(d), and B3LYP/6-31+G(d) levels of theory. Energies were also determined at G3MP2, G3MP2B3, G4MP2, and CBS-QB3 levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Thermodynamic properties (ΔE, ΔH, and ΔG), activation energies, enthalpies, and Gibbs free energies of activation were also calculated for each reaction investigated. All pathways yield an initial tetrahedral intermediate and an intermediate in the last step that dissociates to products via a 1,3-proton shift. At the G3MP2 level of theory, deamination with OH(-) was found to have an activation energy barrier of 155 kJ mol(-1) compared to 187 kJ mol(-1) for the reaction with H(2)O and 243 kJ mol(-1) for GuaH(+) with H(2)O. The lowest overall activation energy, 144 kJ mol(-1) at the G3MP2 level, was obtained for the deamination of guanine with H(2)O/OH(-). Due to a lack of experimental results for guanine deamination, a comparison is made with those of cytosine, whose deamination reaction parallels that of guanine.     

A quantum chemistry study of the Cl atom reaction with formaldehyde

The elementary vapor-phase reaction between Cl atoms and HCHO has been studied by ab initio methods. Calculations at the MP2, MP3, MP4(SDTQ), CCSD, CCSD(T), and MRD-CI levels of theory show that the reaction is characterized by a low electronic barrier; excluding the effects of spin-orbit splitting in Cl, our best estimate at the MRD-CI/aug-cc-pVTZ//RHF-RCCSD(T)/aug-cc-pVTZ level of theory predicts a Born-Oppenheimer barrier height of 0.7 kJ mol-1. The energies of the lowest six electronic states as resulting from MRD-CI calculations are presented at discrete points along the reaction path, and two avoided crossings are found in the transition state region. The spin-orbit splitting in Cl is also calculated along the reaction path; it is not negligible in the transition state region and is found to increase the barrier by only 1.4 kJ mol-1 at the RCCSD(T)/aug-cc-pVTZ transition state geometry. The minimum energy path of the reaction connects an energetically weakly stabilized adduct on the flat potential surface on the reactant side and an energetically strongly stabilized postreaction adduct. The reaction rate coefficient and the kinetic isotope effects were calculated using improved canonical variational theory with small curvature tunneling (ICVT/SCT), and the results were compared to experimental data. The experimental reaction rate coefficient is reproduced within its uncertainty limits by variational transition state theory with interpolated single-point energy corrections (ISPE) at the MP4(SDTQ) level of theory and by conventional transition state theory with interpolated optimized energies (IOE) at the MRD-CI//RCCSD(T) level of theory and interpolated optimized geometries at the RCCSD(T) level of theory on an MP2/aug-cc-pVTZ potential energy surface when employing scaled vibrational frequencies.     

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.     

Decomposition of protonated noradrenaline and normetanephrine assisted by NH2 migration studied by electrospray tandem mass spectrometry and molecular orbital calculations

As part of a research program on neurotransmitters in a biological fluid, the fragmentations characterising catecholamines protonated under electrospray ionisation (ESI) conditions, under low collision energy in a triple-quadrupole mass spectrometer, were investigated. The decompositions of protonated noradrenaline (VH) and normetanephrine (VIH) were studied. Both precursor ions eliminate first H2O at very low collision energy, and the fragmentations of [MH-H2O]+ occur at higher collision energy. The breakdown graphs of [MH-H2O]+ ions, with collision energy varying from 0-40 eV in the laboratory frame, are presented. [VIH-H2O]+ ions lose competitively NH3 and CH3OH. For [VH-H2O]+ the loss of NH3 is dominant while H2O is eliminated at very low abundance at all collision energies. All of these secondary fragmentations are followed at higher collision energies by elimination of CO. These fragmentations are interpreted by means of ab initio calculations up to the B3LYP/6-311+G(2d,2p) level of theory. The elimination of H2O requires first the isomerisation of N-protonated forms, chosen as energy references, to O-protonated forms. The isomerisation barriers are calculated to be lower than 81 kJ/mol above the N-protonated forms. The elimination of NH3 from [MH-H2O]+ requires first the migration, via a cyclisation, of the amine function from the linear chain to the aromatic ring in order to prevent the formation of unstable disubstituted carbocations in the ring. The barriers associated with the loss of NH3 are located 220 and 233 kJ/mol above VH and 219 kJ/mol above VIH. The energy barrier for the loss of ROH is located 236 and 228 kJ/mol above VH and VIH, respectively. The absence of ions corresponding to [VH-2H2O]+ is due to a parasitic mechanism with an activation barrier lower than 236 kJ/mol that leads to a stable species unable to fragment, thus preventing the second loss of H2O. Losses of CO following the secondary fragmentations involve activation barriers higher than 330 kJ/mol.     

氧原子和甲基自由基反应机理的理论研究

用分子轨道从头算和密度泛函理论（DFT）中的B3LYP方法以及适中基组6-311+G(2df,2p)对氧原子与甲基CH3反应进行了系统的研究。计算给出了反应通道上各驻点物种的构型参数、振动频率和能量。结果表明： CH2OH比CH3O稳定,能量约低26.63 kJ/mol,且生成氢和甲醛为其最主要反应通道。     

亚甲基环丙烷卡宾的开环与1,2-氢迁移的从头算研究

用从头计算-解析梯度法研究了亚甲基环丙烷卡宾的重排反应。开环重排经Hückel型过渡态,在过渡态前后分别属2电子和4电子反应,取对旋和顺旋方式,遵从4n+2和4n规则。1,2氢重排经迁移氢进攻卡宾碳的空p轨道发生。用RHF/6-31-G^**∥STO-3G方法算得开环势垒为24kJ/mol。1,2氢重排不是开环的竞争反应,其势垒为118kJ/mol。与环丙烷卡宾比较,亚甲基取代未改变其开环反应的类型与机理,但增大开环反应活性。     

H+NH3反应势能面的量子化学研究

本文采用UCCSD(T)/aug-cc-pVTZ方法研究了H NH3反应势能面,获得了夺氢反应和交换反应过渡态的的几何结构和振动频率。夺氢反应的过渡态具有Cs对称性,其能垒为61.92 kJ/mol。交换反应的过渡态具有C3v对称性,其能垒为39.69 kJ/mol。H NH3发生形成Td对称性的反应中间体NH4里德堡自由基。与夺氢反应相比,交换反应具有更低的反应能垒,并且NH4自由基在反应中可形成长寿命的共振态,和夺氢反应形成竞争关系,因此在H NH3反应的量子动力学研究中必须同时考虑这两类反应。本文还采用更大的基组aug-cc-pVQZ和aug-cc-pV5Z研究了势能面对基组的收敛行为。