CH3S自由基H迁移异构化及脱H2反应的直接动力学研究

采用密度泛函方法（MPW1PW91）在6．311G（d,p）基组水平上研究了CH3S自由基H迁移反应CH3S→CH2SH（R1）,脱H2反应CH3S→HCS＋H2（R2）以及脱H2产物HCS异构化反应HCS→CSH（R3）的微观动力学机理．在QCISD（t）／6．311＋＋G（d,p）／／MPW1PW91／6．311G（d,p）＋ZPE水平上进行了单点能校正．利用经典过渡态理论（TST）与变分过渡态理论（CVT）分别计算了各反应在200-2000K温度区间内的速率常数K^TST和k^CVT,同时获得了经小曲率隧道效应模型（SCT）校正后的速率常数萨k^CVT/SCT．结果表明,反应R1,R2和R3的势垒△E^≠分别为160．69,266．61和241．63kJ／mol。R1为反应的主通道．低温下CH3S比CH2SH稳定,高温时CH2SH比CH3S更稳定．另外,速率常数计算结果显示,量子力学隧道效应在低温段对速率常数的计算有显著影响,而变分效应在计算温度段内对速率常数的影响可以忽略．     

乙硼烷离子和自由基异构体重排与体系势能面的量子化学计算研究

在B3LYP／6－311G（d,p)和QCISD（T）／6－311＋＋G（3df,2p)（单点）水平下，对B2H6^ 阳离子和B2H5自由基全优化得到9个几何异构体，B2H5^＋单态体系（D3h,C1),B2H5^ 三重态（Cs,Cs,C1),B2H5自由基（C2v,Cs,Cs,Cs)，得到势能面上与体系异构化过程相联系的5种过渡态。     

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

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

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

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组，分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化，并在B3LYP/6-311G(d,p)水平上进行了频率分析计算，在 各优化构型上，使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好；1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道，本计算预言了它们的超精细偶合 常数和最稳定构型。     

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

使用密度泛函理论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)方法的超精细偶合常数计算，得到了比以往更好的结果.     

过氧氢自由基HO_2与CH_2S的反应机理及动力学

采用密度泛函理论B3LYP方法,在6-311＋＋G（d,p）基组水平研究了HO2与CH2S的微观反应机理.在CCSD（T）/6-311＋＋G（d,p）//B3LYP/6-311＋＋G（d,p）水平上获得了势能面.利用经典过渡态理论（TST）与变分过渡态理论（CVT）并结合小曲率隧道效应模型（SCT）,计算了反应通道在2...     

BF分子X1∑+,A1∏和B1∑+电子态的势能函数

使用SAC/SAC-CI和D95＋＋,6-311＋＋g,6-311＋＋g^＊＊及D95（d）基组,分别对BF分子的基态X^1∑^＋、第一简并激发态A^1∏和第二激发态B^1∑^＋的平衡结构和谐振频率进行优化计算.对所有计算结果进行比较,得出6-311＋＋g^＊＊基组为最优基组.运用6-311＋＋g^＊＊基组和SAC方法对基态X^1∑^＋,SAC-CI方法对激发态A^1∏和B^1∑^＋进行单点能扫描计算,并用正规方程组拟合Murrell-Sorbie函数,得到相应电子态的势能函数解析式,由得到的势能函数计算了与X^1∑^＋,A^1∏和B^1∑^＋态相对应的光谱常数,结果与实验数据较为一致.     

PF(X~3∑~-)的结构、势能函数与光谱常数

运用多种方法、多种基组对PF(X3∑-)的平衡结构进行优化计算.用QCISD/6-311G(df)方法得到的平衡结构为RPF=0.158 9 nm,与实验值RPF=0.158 97 nm进行比较,最为接近,得出QCISD/6-311G(df)基组为最优基组;然后对PF(X3∑-)进行单点能扫描计算,用正规方程组拟合Murrell-Sorbie函数,得到相应电子态的势能函数解析式,由势能函数计算了与PF(X3∑-)态相对应的光谱常数,结果与实验数据较为一致.这些数据为反应动力学提供了理论依据.     

氟代乙烯阳离子的理论研究

用B3LYP和MP2方法及6-31G(d, p)、6-31＋G(d, p)、6-311G(d, p)和6-311＋G(d, p) 基组，对六种氟代乙烯阳离子做了理论研究，优化了它们的基电子态的结构，计算了对应分子的垂直电离势（VIP）和绝热电离势(AIP).结果表明，与具有非平面结构的乙烯阳离子不同，六种氟代乙烯阳离子都只具有平面结构；与分子结构相比，离子结构的C－C键增长， C－F键缩短， CCF键角变小. 自然布居分析计算表明，这些离子的正电荷主要分布在与F原子相连的C原子和各H原子上. B3LYP/6-311＋G(d, p) 级别上计算的各分子的VIP和AIP值和实验值符合得很好. 使用含弥散基函数的基集可以明显提高这类分子的电离势的计算精度.     

SiC(X~3Ⅱ)的结构与势能函数的从头计算

应用密度泛函理论B3LYP和B3P86以及组态相互作用方法 CCSD,CCSD(T),QCISD和QCISD(T),采用6-311g,6-311G(df),6-311+G(d,p),6-311++G(3df,3pd),aug-cc-pvdz和D95(d)多种基组,优化计算了SiC分子的平衡结构和能量.通过优化计算结果和实验数据R=0.171 82nm进行对比,选择B3LYP/6-311G(df),CCSD/6-311G(df)和QCISD/6-311G(df)方法对SiC(X3Ⅱ)分子进行单点能扫描,同时计算其光谱参数(Be,αe,ωe,ωexe)和力常数(f2,f3,f4),这些计算结果与实验数值相吻合,为研究SiC/SiC复合材料提供了理论数据参考.     

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

用密度函数理论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计算值是至今与实验值符合得最好的理论计算结果。     

优化几何构型对高级别从头算能量的影响

对《ＣＲＣ物理与化学手册》（第７７版）中第三周期以前的无机双原子分子,当其理论优化构型的相对误差大于２％时,分别在实验构型和最大偏差的理论构型下,计算了ＱＣＩＳＤ（Ｔ）／６－３１１＋Ｇ（３ｄｆ,２ｐ）能量并作了比较。结果表明,大多数能量的差别在４．２ｋＪ．ｍｏｌ＾－１以内。由此说明,目前一般采用的构型优化理论方法,多数情况下不至于明显影响单点高级别从头算的计算精度。同时还发现,Ｇ２（ＱＣＩ）的高级     

An Ab Initio Study of the Complexes Formed by the Positive Acetylene Ion with the Hydrogen and Nitrogen Molecules and the Argon Atom

The complexes formed by the positive acetylene ion with the hydrogen molecule, the nitrogen molecule, and the argon atom are investigated with ab initio calculations using the 6-311G** and the 6-31+G(2df,2pd) basis sets. MP2/6-311G** energies and optimum geometries are obtained, as well as single-point MP3, MP4, and QCISD(T) energies with the MP2/6-311G** optimized geometries. Single-point calculations are performed with the 6-31+G(2df,2pd) basis set at MP2/6-311G** optimized geometries.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

An evaluation of various computational methods for the treatment of organoselenium compounds

A reliable computational method for the prediction of organoselenium geometries and bond dissociation energies (BDEs) has been determined on the basis of the performance of density functional theory (DFT: B3LYP and B3PW91) and ab initio molecular orbital procedures (Hartree-Fock (HF)) in conjunction with various Pople basis sets including (but not limited to) the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p), 6-311G(2df,p), and 6-311G(3df,3pd) sets. Predicted geometries and BDEs are compared with available experimental data and quadratic configuration interaction including single and double substitutions (QCISD) results. The B3PW91/6-311G(2df,p) level of theory is recommended for the prediction of the geometries and energetics of organoselenium compounds.     

Theoretical study for the reaction of CH3OCl with Cl atom

A direct dynamics method is employed to study the kinetics of the multiple channel reaction CH(3)OCl + Cl. The potential energy surface (PES) information is explored from ab initio calculations. Two reaction channels, Cl- and H-abstractions, have been identified. The optimized geometries and frequencies of the stationary points and the minimum-energy paths (MEPs) are calculated at the MP2 level of theory using the 6-311G(d, p) and cc-pVTZ basis sets, respectively. The single-point energies along the MEPs are further refined at the G3(MP2)//MP2/6-311G(d, p), G3//MP2/6-311G(d, p), as well as by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MP2/cc-pVTZ geometries. The enthalpies of formation for the species CH(3)OCl and CH(2)OCl are calculated via isodesmic reactions. The rate constants of the two reaction channels are evaluated by using the variational transition-state theory over a wide range of temperature, 220-2200 K. The calculated rate constants exhibit the slightly negative temperature dependence and show good agreement with the available experimental data at room temperature at the G3(MP2)//MP2/6-311G(d, p) level. The present calculations indicate that the two channels are competitive at low temperatures while H-abstraction plays a more important role with the increase of temperature. The calculated k(1a)/k(1) ratio of 0.5 at 298 K is in general agreement with the experimental one, 0.8 +/- 0.2. The high rate constant for CH(3)OCl + Cl shows that removal by reaction with Cl atom is a potentially important loss process for CH(3)OCl in the polar stratosphere.     

An ab initio MO study on the structures and electronic states of hydrogen-bonded O3- HF and SO2-HF complexes

Ab initio molecular orbital (MO) calculations have been carried out for base-hydrogen fluoride (HF) complexes (base = O3 and SO2) in order to elucidate the structures and energetics of the complexes. The ab initio calculations were performed up to the QCISD(T)/6-311++G(d,p) level of theory. In both complexes, hydrogen-bonded structures where the hydrogen of HF orients toward one of the oxygen atoms of bases were obtained as stable forms. The calculations showed that cis and trans isomers exist in both complexes. All calculations for the SO2-HF complex indicated that the cis form is more stable in energy than the trans form. On the other hand, in O3-HF complexes, the stable structures are changed by the ab initio levels of theory used, and the energies of the cis and trans forms are close to each other. From the most sophisticated calculations (QCISD(T)/6-311++G(d,p)//QCISD/6-311+G(d) level), it was predicted that the complex formation energies for cis SO2-HF, trans SO2-HF, cis O3-HF, and trans O3-HF are 6.1, 5.7, 3.4, and 3.6 kcal/mol, respectively, indicating that the binding energy of HF to SO2 is larger than that of O3. The harmonic vibrational frequencies calculated for cis O3-HF and cis SO2-HF complexes were in good agreement with the experimental values measured by Andrews et al. Also, the calculated rotation constants for cis SO2-HF agreed with the experiment.     

类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了类硅烯H2C＝SiLiBr与RH (R＝F, OH, NH2)的插入反应. 在B3LYP/6- 311＋G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C＝SiLiBr与HF, H2O或NH3发生插入反应的机理相同. QCISD/6-311＋＋G(d,p)//B3LYP/6-311＋G(d,p)计算的三个反应的势垒分别为148.62, 164.42和165.07 kJ•mol－1, 反应热分别为－69.63, －43.02和－28.27 kJ•mol－1. 相同条件下发生插入反应时, 反应活性都是H—F＞H—OH＞H—NH2.