几个分子和负离子的QCI势能曲线

用ＱＣＩＳＤ（Ｔ）／６３１１＋Ｇ（３ｄｆ，２ｐ）（简称ＱＣＩ）计算了ＢＨ－（２Π），ＢＦ－（２Π），ＢＨ（１∑），ＢｅＨ－（１∑），ＢｅＦ－（１∑）和ＢＦ（１∑）的较完整的势能曲线，探讨了这些体系的稳定性和解离过程．同时用Ｍｏｒｓｅ函数和ＥＲ函数拟合得到势能曲线，用拟合结果计算谐性频率和非谐性频率，并与已有的实验数据（ＢＨ和ＢＦ分子）进行比较，由此考察了ＱＣＩ理论方法在远离平衡构型时的可靠性．对于因单参考态近似导致的反常势能曲线进行了ＣＡＳＳＣＦ计算．还比较了ＱＣＩ和Ｇ２理论中标准的ＭＰ２（ｆｕｌ）／６３１平衡几何构型及其对ＱＣＩ总能量的影响．     

分子几何构型优化方法的系统性比较

对《ＣＲＣ物理化学手册》第７７版中收集的第三周期以前的所有已知实验构型的无机分子的构型,以ＭＰ２、Ｂ３ＬＹＰ、Ｂ３ＰＷ９１级别上进行了构型优化的系统性比较,优化采用基组为６－３１Ｇ（ｄ,ｐ）,６－３１１Ｇ（ｄ,ｐ）,６－３１１Ｇ（２ｄ,ｐ）。对大多数分子另比较了ＱＣＩＳＤ（Ｔ）方法,对多原子分子比较了ＢＰＷ９１方法,对含氢双原子分子也使用ＱＣＩＳＤ方法。结果表明,键长的平均绝对偏差（单位：ｐｐｍ）     

金属硼化物结构与稳定性的理论研究

用ＨＦ／３－２１Ｇａｂｉｎｉｔｉｏ法对金属硼化物ＭＢ２／ＭＢ２（Ｍ＝Ｌｉ，Ｎａ，Ｂｅ，Ｍｇ，Ａｌ）的７５个电子态结构进行能量梯度法优化，再用大基组二次且上互作用ＱＣＩＳＤ（Ｔ）／６－３１１Ｇ进行单点计算，得到了结构参数总能量，为了考察各原子簇的稳定性，还对２４个碎片的７０多个电子态，求得相应的ＱＣＩＳＤ（Ｔ）能量，在此基础上计算了原子化能、电离能、离解通道和碎片化能，得到了原子簇的稳定性规律。     

夹心羰基化合物（C5H5）2Ti（CO）2的电子结构及化学键研究

用量子化学从头算方法（ＳＴＯ－３Ｇ基组）对（Ｃ５Ｈ５）２Ｔｉ（ＣＯ）２进行几何优化，得到平衡几何构型，并在此基础上采用ＳＴＯ－３－２１Ｇ基组对（Ｃ５Ｈ５）２Ｔｉ（ＣＯ）２进行单点计算，计算结果表明：（Ｃ５Ｈ５）２Ｔｉ（ＣＯ）２的ＨＯＭＯ具有Ｔｉ→ＣＯ反馈π键性质，理论计算结果与实验相符。     

P14原子簇的理论研究

以Ｐ_２、Ｐ_４（Ｔｄ）、Ｐ_８（Ｃ_（２ｖ））及Ｐ_（１０）（Ｃ_（２ｖ））为结构单元设计了Ｐ_２＋Ｐ_４（Ｔ_ｄ）＋Ｐ_８（Ｃ_（２ｖ））→Ｐ_（１４）（Ｃ_ｓ）（Ａ）、Ｐｓ（Ｃ_（２ｖ）＋３Ｐ_２→Ｐ_（１４）（Ｃ_（２ｖ））（Ｉ）（Ｂ）、Ｐ_（１０）（Ｃ_（２ｖ）＋２Ｐ_２→Ｐ１４（Ｃ_（２ｖ））（Ⅱ）（Ｃ），利用从头算Ｇａｕｓｓｉａｎ－９４程序，选择６－３１Ｇ基组，对４种结构单元及Ｐ_（１４）原子簇的３种几何构型Ｃｓ、Ｃ_（２ｖ）（Ⅰ）、Ｃ_（２ｖ）（Ⅱ）进行全优化。相对能量的计算结果表明，Ｐ_（１４）（Ｃ_ｓ）与Ｐ_（１４）（Ｃ_（２ｖ））（Ⅰ）构型稳定。进一步设计（７／２）Ｐ_４（Ｔ_ｄ）→Ｐ_（１４）（Ｄ）及７Ｐ_２→Ｐ_（１４）（Ｅ），其相对能量表明，Ｐ_（１４）（Ｃ_ｓ）构型是稳定的．磷与磷连结单键键长范围为０．２２０～０．２２８ｎｍ，双键键长范围为０．２００～０．２０２ｎｍ，与实验结果大体相当。     

H＋NF_2→NF（α）＋HF反应的从头算研究

本文用量子化学从头算法对Ｈ＋ＮＦ_２→ＮＦ（α）＋ＨＦ反应进行了理论研究，采用ＭＰ_２（ＦＵＬＬ）／６－３１Ｇ（ｄ）方法计算了反应物、产物、中间体和过渡态的构型参数和振动频率，然后用Ｇａｕｓｓｉａｎ－２理论计算了各物种的精确能量，计算结果表明：反应在单重态势能面上进行，电子激发态ＮＦ（α）通过复合－消除机理产生，用Ｇ２理论计算２９８Ｋ下反应的焓变为－１１９．４ｋＪ／ｍｏｌ，过渡态比反应物能量低９８．４ｋＪ／ｍｏｌ．中间体ＨＮＦ２是势能面的最低点，比反应物的能量低３１３．５ｋＪ／ｍｏｌ．     

C70R2（R=OH,CH30的结构和电子光谱

用ＩＮＤＯ方法研究了Ｃ７０Ｒ２（Ｒ＝ＯＨ，ＣＨ３）４种异构体的结构和稳定性，表明１，９－Ｃ７０（ＯＨ）２比７，８－Ｃ７０（ＯＨ）２稳定，两者能量差为３８．５ｋＪ／ｍｏｌ，而７，８－Ｃ７０（ＣＨ３）２比１，９－Ｃ７０（ＣＨ３）２能量低２３．０ｋＪ／ｍｏｌ。以优化构型为基础，对Ｃ７０Ｒ２（Ｒ＝ＯＨ，ＣＨ３）的电子光谱进行了理论预测。     

1,2-二硒方酸的从头计算

在ＲＨＦ／ＳＴＯ－３Ｇ和ＳＴＯ－３Ｇ＊水平上用ａｂｉｎｉｔｉｏＳＣＦ方法优化得到１，２－二硒方酸（３，４－二羟基－３－环丁烯－１，２－二硒酮）三种平面异构体的平衡几何构型，发现三种平面异构体中ＺＺ型是能量最低构象。用ａｂｉｎｉｔｉｏ数值方法在ＲＨＦ／ＳＴＯ－３Ｇ＊水平上计算了三种平面异构体的谐振动频率。     

H＋NF2→NF（a）＋HF反应的从头算研究

本文用量子化学从头算法对Ｈ＋ＮＦ２→（ａ）－ＨＦ反应进行了理论研究，采用ＭＰ２（ＦＵＬＬ）／６－３１Ｇ（ｄ）方法计算了反应物，产物，中间体和过渡态的构型参数和振动频率，然后用Ｇａｕｓｓｉａｎ－２理论计算了各种物的精确能量。     

ClnAlNHn和HnAlNHn(n=1～3)的振动光谱与化学键性质的研究

用从头计算分子轨道法和密度泛函理论，在ＨＦ／６－３１Ｇ＊和Ｂ３ＬＹＰ／６－３１Ｇ＊水平上对ＣｌｎＡｌＮＨｎ和ＨａＡｌＮＨｎ（ｎ＝１～３）及其碎片分子的几何构型、电子结构、振动光谱和化学热力学性质进行了理论研究。结果表明，优化几何参数与实验值相。ＣｌＡｌＮＨ２和Ｈ２ＡｌＮＨ２分子中，Ａｌ－Ｎ键为由一个σ键和一个π键组成的双重键，旋转势垒分别为３４．１０和５４．３５ＫＪ．ｍｏｌ＾－１，而Ｃｌ３ＡｌＮＨ     

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.     

低硝酸盐离子N2O22-和亚硝酸HONO的异构体及其异构化反应的量化研究

采用Gaussian-98程序进行，在HF／6—31G(d)，B3LYP／6—31G(d)和MP2／6—31G(d)水平下优化分子结构并寻找过渡态，对于MP2／6—31G(d)结果在QCISD(T，E4T)，MP4／6—311 G(d，p)，MP4／6—311 G(2df，p)水平下重新计算能级．并用内禀反应坐标(IRC)法研究了N2O2^2-和亚硝酸HONO的异构化反应机理。     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

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.     

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.     

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.     

CH3 CH2 +N(4S)反应机理的理论研究

The reaction for CH3CH2+N(4S) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single point calculations for all the stationary points were carried out at the QCISD(T)/ 6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2CH2+3NH and H2CN＋CH3, and the minor products are the CH3CHN+H in the reaction. The majority of the products CH2CH2+3NH are formed via a direct hydrogen abstraction channel. The products H2CN＋CH3 are produced via an addition/dissociation channel. The products CH3CHN+H are produced via an addition/dissociation channel.     

Theoretical studies on dynamics and thermochemistry of the reactions CF(3)CHCl(2)+Cl-->CF(3)CCl(2)+HCl and CF(3)CHFCl+Cl-->CF(3)CFCl+HCl

A dual-level direct dynamics study has been carried out for the two hydrogen abstraction reactions CF(3)CHCl(2)+Cl and CF(3)CHFCl+Cl. The geometries and frequencies of the stationary points are optimized at the BHLYP/6-311G(d,p), B3LYP/6-311G(d,p), and MP2/6-31G(d) levels, respectively, with single-point calculations for energy at the BHLYP/6-311++G(3df,2p), G3(MP2), and QCISD(T)/6-311G(d,p) levels. The enthalpies of formation for the species CF(3)CHCl(2), CF(3)CHFCl, CF(3)CCl(2), and CF(3)CFCl are evaluated at higher levels. With the information of the potential energy surface at BHLYP/6-311++G(3df,2p)//6-311G(d,p) level, we employ canonical variational transition-state theory with small-curvature tunneling correction to calculate the rate constants. The agreement between theoretical and experimental rate constants is good in the measured temperature range 276-382 K. The effect of fluorine substitution on reactivity of the C-H bond is discussed.     

N(4S)+CH3X(X=Cl、Br)反应机理的理论研究

The reaction of N(4S)+CH3X(X=Cl、Br) was studied by the ab initio method. The geometries of the reactants, transition states and products were optimized at the MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the MP2/6-311++G(3df,2p) and the QCISD(T)/6-311+G(d,p) levels using the MP2/6-311+G(d,p) optimized geometries. The energies of all the stationary points were calculated by the G2MP2 method. The results of this theoretical study indicate that the reaction has three reaction channels: H abstraction reaction channel a, Cl or Br abstraction reaction channel b and substitution reaction channel c. For the N(4S)+CH3Cl reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have a slight contribution in the reaction. For the N(4S)+CH3Br reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have some contribution in the reaction.