强离子键体系电子相关能简捷计算方案的应用

The calculation results of electron correlation energies of KF and （KF）2 were reported. The transferability of 1s^2 K , 1s^2 F and the inner core correlation effects of K and F in both K, K^＋, KF and F, F^-, KF systems were investigated respectively. The correlation energy contributions of K and F component to KF system were calculated. By applying the simple estimation scheme to the calculation of the correlation energy of the strong ionic compound KF and （KF）2, it was shown that such a powerful scheme could not only reach the chemical accuracy but also need little computational work.     

锂(100)面上H_2H+H反应动力学的理论研究

本文在H-Li(100)面吸附扩散的ab initio SCF势能面基础上构造了(H_2H+H)/Li(100)面相互作用的推广LEPS势能面,并用QCT方法研究了该体系的反应动力学行为。分析势能面特征得到:H_2在Li(100)面上的吸附无需活化能,H_2在Li(100)面上的解离吸附与吸附位及吸附模式密切相关,H_2的卧式解离比立式解离要容易得多。分析各种碰撞轨迹得到:低覆盖度下双氢原子的表面复合几率很小,H_2的表面解离几率受到H_2振动量子数的控制。本文构造了一种适合于动力学研究的气体-金属表面相互作用势能面,并且,动力学QCT计算结果能够对H_2表面活化的分子束实验作出合理的解释。     

甲基硅乙烯异构化反应的理论研究 用从头算结果讨论其热力学、动力学性质

The titled isomerizaton of methylsilene has been studied at ab initio level. Theharmonic vibrational frequencies of reactant, product and transition state, the barrier energy, reaction heat, equilibrium constant and rate constant are calculated. The results indicate that this reaction is exothermic, spontaneous and high temperature above 900 K is favourable for the isomerization. The prediction is in agreement with experiment.     

变分过渡态理论对H＋H2及其同位素选态反应的研究

将选态速度常数的计算推广到任意指定反应物、过渡态的振动激发态.用此法计算了H+H_2(v)及其同位素经不同振动激发过渡态时的速度常数,发现弯曲振动模激发所得结果与实验值更符合,并且在给定能量下,过渡态的弯曲振动模激发比其对称伸缩模激发更有利于反应进行.     

H+H2(v=0,1)→H2(v′)+H反应截面的电子非绝热三维量子散射近似计算

本文把电子非绝热一维量子散射反应几率和三维量子散射反应截面的近似公式结合起来, 对于反应物分子(H_2)不同的量子振动态(v=0, 1) 分别计算了H+H_2(v=0)→H_2(v′=0, 1)+H和H+H_2(v=1)→H_2(v′=0, 1)+H的平均反应截面σ_0和σ_1, 并同文献上用电子绝热理论计算的结果作了比较, 表明对这类中性原予-分子反应碰撞的过程, 特别是当反应物分子处于振动激发态时, 电子非绝热效应是存在的。     

OH+C2H2与OH+C2H4加成反应的热力学和动力学

OH基与乙炔(乙烯)的反应是控制大气中OH基浓度的重要化学反应。对反应OH+C_2H_2,1975年Davis等用FP-RF技术测定了反应的速率常数。1977年,Perry又用同样方法对该反应作出了研究;结果表明,此反应的速率常数强烈地依赖于压力,与Davis等人的实验结果不符。对反应OH+C_2H_4,Atkinson等人的实验研究表明,压力在30.0kPa以下,反应的速率常数随压力而改变;而在30.0—88.4 kPa之间,与压力无关.这与前人的结果不同。两个反应的产物也因温度不同而异.此外,上述反应一般是在近101kPa下发生的,而实验结果大多是在低压下获得的。且目前尚无直接的理论计算结果。为此,我们在从头算水平上用统计热力学方法,对反应     

甲酰氯异构化和光解的理论研究——从头算水平上的反应热力学、动力学分析

由于对化学反应体系解核运动Schrdinger方程困难,化学界常常求解参加反应粒子的电子运动Schrdinger方程,得到其电荷密度分布、自由价和键级等静态物理量,作为反应性指标;这虽有简明的优点,但往往缺乏传递性和统一性。如用电荷密度分布分析亲电反应,仅在同系物中有定性比较其反应性相对大小的意义,无定量价值。此外,用反应的势垒大小判断反应的难易,忽略了鞍点处曲率的影响,就动力学的Arrhenius速度常数公式看,     

两个典型氧桥三核铁(Ⅲ)配合物[Fe3O(O2CCH2OC6H5)6(3H2O)]和[Fe3O(TIEO)2(O2CPh)2Cl3]的局域自旋和磁性质的理论探讨

用从头算波函数(UHF或UDFT波函数)代替ZILSH方法中的半经验波函数得到了ABLSH方法，接着用该方法研究了两个典型氧桥三核铁(Ⅲ)配合物[Fe₃O(O₂CCH₂OC₆H₅)₆(3H₂O)]和[Fe₃O(TIEO)₂(O₂CPh)₂Cl₃]的局域自旋和磁性质。通过计算得到的局域自旋结果和前人的具有可比性，同时所得的磁交换耦合常数和实验值很吻合。该方法可作为研究海森堡型磁性系统(HM)的新工具。     

辛算法在准经典轨迹计算中的应用: H+H2体系

采用辛算法模拟了H + H₂在BKMP势能面上的长时间动力学行为, 并和Runge-Kutta方法计算的结果进行了比较.结果表明,采用Runge-Kutta方法模拟长时动力学行为时体系中出现动力学耗散行为,导致无法真实的模拟体系的物理实际;而采用辛算法模拟精确度更高,四阶辛算法误差只有Runge-Kutta法的约百分之一,六阶辛算法的误差则更小,体系中也未见动力学耗散行为.这些结果说明,在进行长时动力学行为模拟的时候,在保持体系能量守恒性上,辛算法应当是一种更好的积分方法.     

磁性体系的定域自旋研究: 模型化合物HHeH

根据Davidson和 Clark的定域自旋(Local Spin)理论的计算方程, 将它们改写为简单易行的矩阵代数形式的方程, 并根据这些方程编写了一个计算磁性体系定域自旋的程序. 为了验证程序的正确性, 在相同的计算水平下, 计算了O₂和NO的定域自旋结果并把它们和Davidson和Clark的做了比较, 结果证明程序的可行性. 进而针对简单模型磁性化合物HHeH, 计算它的定域自旋结果和磁交换耦合常数J, 探讨它们与体系几何参数之间的关系.