丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

L-缬氨酸旋光异构的两种光反应可能途径

利用密度泛函理论(DFT)和从头算分子轨道理论对L-缬氨酸的旋光异构光反应机理进行了研究. 分别用B3LYP和MP2方法在6-311++G(d, p)基组级别上全优化得到了S₀和T₁态反应路径上的反应物、产物、中间体以及过渡态结构的几何构型, 给出了反应能垒, 利用含时密度泛函理论(TD-DFT)中的B3LYP/6-311++G(d, p)方法优化得到了S₁态反应路径上的平衡态几何构型. 通过分析反应途径上各个驻点的几何构型特征, 确定了L-缬氨酸在激发态可能通过手性碳上的氢原子以羰基氧或氨基氮为中转媒介发生质子迁移来完成旋光异构反应. 进一步用自洽反应场理论中的极化连续模型(PCM)方法探讨了溶剂化效应对旋光异构反应机理的影响.     

2-(2-羟基苯亚甲基胺)-4,6-二羟基-嘧啶的质子转移异构化反应

为了探索2-(2-羟基苯亚甲基胺)-4,6-二羟基-嘧啶(M1)分子醇式和酮式结构互变异构化的反应机理,利用密度泛函理论(DFT)方法,在B3LYP/6-311+G(d,p)基组水平上,对M1化合物异构化反应的势能面进行了研究,在探讨各种可能的反应途径中,发现单体至少有8种异构体和10种过渡态.结果表明:2-(2-羟基苯亚甲基胺)-6-羟基-4(3H)嘧啶酮(M6)不论是单体、与水形成的配合物,还是二聚体,比其相对应的异构体能量低,表明在通常情况下是以M6形式稳定存在的;在考察的可能反应途径中,直接进行的分子内质子转移过程需要的活化自由能为143.8 kJ· mol-1,水助催化时,反应的活化自由能为38.9 kJ· mol-1,二聚体双质子转移的活化自由能为0.6 kJ·mol-1,二聚体双质子转移所需活化自由能最低,在室温下就可以进行,由此可见氢键在降低反应活化能方面起着重要的作用.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

A theoretical study of the kinetics of the benzylperoxy radical isomerization

The rate constant of the benzylperoxy isomerization reaction has been computed using 54 different levels of theory and has been compared to the experimental value reported at 773 K. The aim of this methodology work is to demonstrate that standard theoretical methods are not adequate to obtain quantitative rate constants for the reaction under study. The use of the elaborated CASPT2 method is essential to estimate a quantitative rate constant. Geometry optimizations and vibrational frequency calculations are performed using three different methods (B3LYP, MPW1K, and MP2) and six different basis sets (6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), and cc-pVDZ). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the 6-31G(d,p) basis set, and with the CASPT2 level of theory with the ANO-L-VDZP basis set. Canonical transition-state theory with a simple Wigner tunneling correction is used to predict the high-pressure limit rate constants as a function of temperature. We recommend the use of the CASPT2/ANO-L-VDZP//B3LYP/cc-pVDZ level of theory to compute the temperature dependence of the rate constant of the four-center isomerization of the benzylperoxy radical. It is given by the following relation: k(600-2000 K) (in s (-1)) = (1.29 x 10 (10)) T (0.79) exp[(-133.1 in kJ mol (-1))/ RT]. These parameters can be used in the thermokinetic models involving aromatic compounds at high pressure. This computational procedure can be extended to predict rate constants for other similar reactions where no available experimental data exist.     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

The determination of the equilibrium structures of oxygen,ozone, and hydrogen peroxide using the ab initio and density functional theory methods

The geometries and energies of small oxygen containing molecules are studied by both the ab initio and density functional theory (DFT) methods. The RHF, MP2, and QCISD(T) ab initio methods, BHandH, BHandHLYP, BeckeSLYP, Becke3P86 DFT hybrid methods, BLYP, and the BP86 non-local DFT methods with the 3-21G¹, 6-31G(d,p), 6-311 + G(2d,2p) and 6-311 + + G(3df,3pd) basis sets were used for the computational study. The obtained results from the different methods were compared to the experimental values. The suitability of the DFT methods for reproducing experimental data were discussed.     

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters, but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.     

A computational study of conformational interconversions in 1,4-dithiacyclohexane (1,4-dithiane)

Ab initio molecular orbital theory with the 6-31G(d), 6-31G(d,p), 6-31+G(d), 6-31+G(d,p), 6-31+G(2d,p), 6-311G(d), 6-311G(d,p), and 6-311+G(2d,p) basis sets and density functional theory (BLYP, B3LYP, B3P86, B3PW91) have been used to locate transition states involved in the conformational interconversions of 1,4-dithiacyclohexane (1,4-dithiane) and to calculate the geometry optimized structures, relative energies, enthalpies, entropies, and free energies of the chair and twist conformers. In the chair and 1,4-twist conformers the C-Hax and C-Heq bond lengths are equal at each carbon, which suggest an absence of stereoelectronic hyperconjugative interactions involving carbon-hydrogen bonds. The 1,4-boat transition state structure was 9.53 to 10.5 kcal/mol higher in energy than the chair conformer and 4.75 to 5.82 kcal/mol higher in energy than the 1,4-twist conformer. Intrinsic reaction coordinate (IRC) calculations showed that the 1,4-boat transition state structure was the energy maximum in the interconversion of the enantiomers of the 1,4-twist conformer. The energy difference between the chair conformer and the 1,4-twist conformer was 4.85 kcal/mol and the chair-1,4-twist free energy difference (deltaG degrees (c-t)) was 4.93 kcal/mol at 298.15 K. Intrinsic reaction coordinate (IRC) calculations connected the transition state between the chair conformer and the 1,4-twist conformer. This transition state is 11.7 kcal/mol higher in energy than the chair conformer. The effects of basis sets on the 1,4-dithiane calculations and the relative energies of saturated and unsaturated six-membered dithianes and dioxanes are also discussed.     

Thermal behavior of [2.1.1]propellane: a DFT/ab initio study

Density functional and ab initio molecular orbital calculations have been used to search for the low energy path of the thermal isomerization of [2.1.1]propellane 1. Three reaction modes were considered: ring opening of the bicyclo[1.1.0]butane unit in 1 to give 1,2-dimethylenecyclobutane 21, opening of the four-membered ring of 1 to afford 1,3-dimethylenecyclobutane 20, and breaking of the [2.1.1]propellane central bond and one of the bicyclo[1.1.0]butane side bonds to form carbene 17. At the CAS(12,12)PT2N/6-31G(d) level of theory, the activation barrier of the latter route was lowest in energy. Further investigation of this process at the QCISD(T)/6-311G(d,p)//QCISD/6-31G(d) and B3PW91/6-311G(d,p)// B3PW91/6-311G(d,p) level of theory indicated that the barrier of isomerization of 1 --> 17 amounts to 29 kcal/mol and that 17 is stabilized by hydrogen migration to give dienes 18 and 19.     

Investigations of silicon-nitrogen hydrides from reaction of nitrogen atoms with silane: experiments and calculations

Using a quadrupole mass filter and vacuum-ultraviolet ionization, we measured the time-of-flight spectra of species at mass-to-charge ratios of m/ z = 45-42 from the reaction of N + SiH 4 in crossed molecular beams. Species with m/ z = 44 and 43 correspond to reaction products HSiNH/SiNH 2 and HSiN/HNSi, respectively; species with m/ z = 45 and 42 are assigned to isotopic variants and daughter ions, respectively, of those two reaction products. We measured the photoionization yields and branching ratios for dissociative ionization of reaction products as a function of photoionization energy. The ionization thresholds of products HSiNH/SiNH 2 and HSiN/HNSi were determined to be 6.7 and 9.2 eV, respectively. Furthermore, we calculated the equilibrium structures, electronic energies, and vibrational wavenumbers of various silicon-nitrogen hydrides H x SiNH y ( x + y = 0-3) using quantum-chemical methods. SiNH 2 (X (2)B 2) and HNSi (X (1)Sigma (+)) are more stable than HSiNH (X (2)A') and HSiN (X (1)Sigma (+)) by 0.82 and 2.81 eV, respectively. SiNH 2 (X (2)B 2), HSiNH (X (2)A'), HNSi (X (1)Sigma (+)), and HSiN (X (1)Sigma (+)) have adiabatic ionization energies of 6.81, 8.19, 10.21, and 10.23 eV, respectively. These experimental and calculated results indicate that SiNH 2 (X (2)B 2) and HNSi (X (1)Sigma (+)) are dominant among isomeric products in the reaction of N + SiH 4. This work presents the first observation of products from the reaction of N + SiH 4 in crossed beams and extensive calculations on pertinent silicon-nitrogen hydrides.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Ab initio quantum chemical studies of reaction mechanism for CH2CO with NCO

An extensive quantum chemical study of the potential energy surface (PES) for all possible isomerization and dissociation reactions of CH₂CO with NCO is reported at the DFT (B3LYP/6-311++G(d,p)) and CCSD(T)/cc-pVDZ//B3LYP/6-311++G(d,p) levels of theory. For the CH₂CO+NCO reaction, the formation of CO+CH₂NCO via an addition–elimination mechanism is the dominant channel on the doublet surface. While the formation of CO+CH₂OCN via bimolecular substitution reaction is in the secondary. Meanwhile, the isomerization and dissociation reactions of the products, CH₂NCO and CH₂OCN, also have been investigated using the same theoretical approach. It can be concluded that these reaction channels are not feasible kinetically at low or fairly high temperatures. On the basis of the ab initio data, the total rate constants for the CH₂CO+NCO reaction in the T=296–560 K range have been computed using conventional transition state theory with Wigner tunneling correction and fitted by a rate expression as k=2.14×10⁻¹² (cm³ molecule⁻¹ s⁻¹) exp(654.29/T). The calculated total rate constants with Wigner tunneling correction for the CH₂CO+NCO reaction are in good agreement with the available experimental values.     

Vibrational frequencies and infrared intensities of the hydrogen-bonded complexes of nitrous acid with ethers: ab initio and DFT studies

The vibrational spectra of the binary complexes formed by HONO-trans and HONO-cis with dimethyl and diethyl ethers have been investigated using ab initio calculations at the SCF and MP2 levels with 6-311++G(d,p) basis set and B3LYP calculations with 6-31G(d,p) and 6-31+G(d,p) basis sets. Full geometry optimisation was made for the complexes studied. The accuracy of the ab initio calculations have been estimated by comparison between the predicted values of the vibrational characteristics (vibrational frequencies and infrared intensities) and the available experimental data. It was established, that the methods, used in this study are well adapted to the problem under examination. The predicted values with the B3LYP calculations are very near to the results, obtained with 6-311++G(d,p)/MP2. The ab initio and DFT calculations show that the changes in the vibrational characteristics (vibrational frequencies and infrared intensities) upon hydrogen bonding for the hydrogen-bonded complex (CH3)2O...HONO-trans are larger than for the complex (CH3)2O...HONO-cis.     

Ab initio study of the reaction path of the reaction HNCO+CN

HNCO is a convenient photolytic source of NCO and NH radicals for laboratory kinetics studies of elementary reaction[1] and plays an important role in the combustion and atmosphere chemistry. It can re- move deleterious compounds rapidly from exhausted ga…     

Molecular structure of the octatetranyl anion, C8H-: a computational study

The equilibrium molecular structure of the octatetranyl anion, C8H(-), which has been recently detected in two astronomical environments, is investigated with the aid of both ab initio post-Hartree-Fock and density functional theory (DFT) calculations. The model chemistry adopted in this study was selected after a series of benchmark calculations performed on molecular acetylene for which accurate gas-phase structural data are available. Geometry optimizations performed at the CCSD/6-311+G(2d,p), QCISD/6-311+G(2d,p), and MP4(SDQ)/6-311+G(2d,p) levels of theory yield for C8H(-) an interesting polyyne-type structure that defies the chemical formula displaying a simple alternation of triple and single carbon-carbon bonds, [:C[triple bond]C-C[triple bond]C-C[triple bond]C-C[triple bond]CH](1-). In the optimized geometry of C8H(-), as one proceeds from the naked carbon atom on one side of the chain to the CH unit on the opposite side of the chain, the short (formally triple) carbon-carbon bonds decrease in length from 1.255 to 1.213 A whereas the long (formally single) carbon-carbon bonds increase (albeit only slightly) in length from 1.362 to 1.378 A (CCSD results). In striking contrast, both MP2 and DFT (B3LYP and PBE0) calculations fail in reproducing the pattern of the carbon-carbon bond lengths obtained with the CCSD, QCISD, and MP4 methods. The structures of three shorter n-even chains, C(n)H(-) (n = 2, 4, and 6), along with those of four n-odd compounds (n = 3, 5, 7, and 9) are also investigated at the CCSD/6-311+G(2d,p) level of theory.     

Theoretical study of the regioselectivity of the cycloaddition reaction between cyclopentadiene and methyleneketene

Three possible reaction schemes for the cycloaddition reaction between methyleneketene and cyclopentadiene were studied by DFT (density functional theory) and ab initio calculations. All of these cycloaddition reactions are exothermic, concerted but nonsynchronous processes. The computed activation energies indicate that the norbornene product yielded from a 1,2-addition of methyleneketene with cyclopentadiene (reaction (1)) is the primary product. The performance of various computational methodologies, MP2, MP4, and CCSD(T), in conjunction with a wide array of basis sets, 6-31G(d), 6-311+G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ, in obtaining reliable activation and reaction energies of the reactions under investigation has been critically analyzed.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.     

An assessment of theoretical procedures for predicting the thermochemistry and kinetics of hydrogen abstraction by methyl radical from benzene

The reaction enthalpy (298 K), barrier (0 K), and activation energy and preexponential factor (600-800 K) have been examined computationally for the abstraction of hydrogen from benzene by the methyl radical, to assess their sensitivity to the applied level of theory. The computational methods considered include high-level composite procedures, including W1, G3-RAD, G3(MP2)-RAD, and CBS-QB3, as well as conventional ab initio and density functional theory (DFT) methods, with the latter two classes employing the 6-31G(d), 6-31+G(d,p) and/or 6-311+G(3df,2p) basis sets, and including ZPVE/thermal corrections obtained from 6-31G(d) or 6-31+G(d,p) calculations. Virtually all the theoretical procedures except UMP2 are found to give geometries that are suitable for subsequent calculation of the reaction enthalpy and barrier. For the reaction enthalpy, W1, G3-RAD, and URCCSD(T) give best agreement with experiment, while the large-basis-set DFT procedures slightly underestimate the endothermicity. The reaction barrier is slightly more sensitive to the choice of basis set and/or correlation level, with URCCSD(T) and the low-cost BMK method providing values in close agreement with the benchmark G3-RAD value. Inspection of the theoretically calculated rate parameters reveals a minor dependence on the level of theory for the preexponential factor. There is more sensitivity for the activation energy, with a reasonable agreement with experiment being obtained for the G3 methods and the hybrid functionals BMK, BB1K, and MPW1K, especially in combination with the 6-311+G(3df,2p) basis set. Overall, the high-level G3-RAD composite procedure, URCCSD(T), and the cost-effective DFT methods BMK, BB1K, and MPW1K give the best results among the methods assessed for calculating the thermochemistry and kinetics of hydrogen abstraction by the methyl radical from benzene.