双原子分子激发态势能的自旋相关局域Hartree-Fock密度泛函理论方法

基于自旋相关局域Hartree-Fock (SLHF)势函数，本文提出了一种计算双原子分子激发态势能的密度泛函理论(DFT)方法，并将该方法应用于和的激发态势能曲线的计算。在只考虑交换能的情况下，本文的DFT计算结果与文献中精确方法和Hartree-Fock (HF)方法的结果符合的非常好，说明采用SLHF势函数作为交换势的DFT方法是一个很好的计算激发态势能的方法。本文还计算和探讨了电子的关联势函数和关联能，发现传统的近似方法在较大核间距的情况下大大低估了电子的关联能.     

OH+C2H3F的反应机理及产物支化比反应的理论研究

本文对HOC₂H₃F可能解离通道的势能面进行从头算CCSD(T)/CBS//B3LYP/6-311G(d，p)计算，同时对速率常数进行Rice-Ramsperger-Kassel-Marcus计算. 生成主要产物CH₂CHO+HF最有利的反应途径是OHC₂H₃F→i2→TS14→i6→TS9→i3→TS3→CH₂CHO+HF，其中速率决定步骤是HF通过TS11从CO桥接位置解离，能量比反应物高3.8 kcal/mol. 借助中间态TS14，F原子从C_α迁移到C_β位置生成CH₂O+CH₂F，然后通过中间态TS16，H从O迁移到C_α位置；通过中间态TS5，C-C键断裂生成产物，其能量比反应物低1.8 kcal/mol，比TS11低4.0 kcal/mol.     

BeF2分子基态势能函数研究

利用密度泛函理论（DFT）B3LYP/aug-cc-PVTZ方法对BeF和BeF2分子进行计算, BeF分子的基态电子态为X2Σ+,平衡核间距Re=0.1367 nm,BeF2分子最稳定构型为D∞h构型,基态电子态为X1Σg+,离解能De=13.4 eV。利用最小二乘法拟合了BeF分子基态的Murrell-Sorbie势能函数,并得到其光谱数据与力常数。结合全局多体项展式方法,导出了基态BeF2分子基态势能函数的解析表达式,该函数准确反映了BeF2分子的结构以及静态反应特征,这些结果为进一步探索BeF2的微观分子反应动力学提供了基础。     

大挠度弯曲直梁混合变量最小势能原理的应用

应用大挠度弯曲直梁混合变量最小势能原理,求解均载两端固定大挠度柱面弯曲板条的轴向挠度分布和轴向弯矩分布．实例计算表明:该方法简单实用、精度高,是一种计算大挠度柱面弯曲板条变形的有效方法．     

压电材料椭圆切口的力学分析

在线性压电本构方程框架下，用复势函数方法对椭圆切口模型进行了精确的数值计算。完整地考虑了各向异性力电耦合效应以及切口内不同电介质的介电性质。给出了切口内部不同介电性质对压电材料内部应力的影响。指出了Sosa文章里的一些计算错误。由于现在文献中很少有关于电导通边界条件下理论解的数值结果，所以本文同时提供了不同电边界条件下理论解的数值结果，所以本文同时提供了不同电边界条件下的理论解的数值结果。最后通过最小势能原理建立了8结点有限元模型，对椭圆切口问题进行了计算并与理论解进行了仔细比较。     

OCl2分子基态势能曲线的从头算研究

采用三重激发耦合族理论(CCSD),使用多种基组对OCl2分子基态结构进行几何优化,选出最优基组6-311G(2df)对OCl2分子的离解能、谐振频率和力常数等进行了计算,其结果与实验值符合很好.在此基础上推导出OCl2分子基态的多体展式势能函数,其等值势能曲线正确反应了OCl2分子的结构特征及势阱深度,进一步讨论了Cl+OCl和O+ClCl分子反应的势能面特征.这些结果可用于微观反应动力学的研究.     

ClF-分子离子的结构与势能函数

采用单、双取代包括三重激发的二次组态相互作用[QCISD(T)]方法和单、双取代包括非迭代三重激发的耦合簇理论[CCSD(T)]方法, 结合相关一致基组aug-cc-pVXZ (X=D, T, Q, 5)对基态³⁵ClF^-和³⁷ClF^- (X²Σ⁺)分子离子进行了结构优化计算. 对CCSD(T)方法的计算结果用四种方法分别外推至基组极限, 得到了体系在基组极限的平衡结构常数. 在CCSD(T)/aug-cc-pVXZ (X=D, T, Q, 5)理论水平进行了单点能扫描. 对扫描计算结果进行基组外推并用Murrell-Sorbie 势能函数拟合得到了体系的解析势能函数表达式, 并进一步得到了³⁵ClF^-和³⁷ClF^-的光谱常数. 拟合所得势能曲线准确地再现了其离解能和平衡结构特征. 对ClF 中性自由基采用完全相同的理论方法进行了计算. 所得结果与有关文献中的实验结果符合得很好, 而且在一定程度上证明了将该理论方法应用于ClF^-分子离子的计算是合适而可靠的. ClF 自由基的优化计算结果还被用于计算其电子亲和能.ClF^-的垂直解离能也同时计算得出. 基于ClF^-的结构优化和单点能扫描计算结果, 通过求解核运动的径向薛定谔方程, 得到了无转动³⁵ClF^-和³⁷ClF^-(X²Σ⁺)的全部振动态及相应的分子常数.     

基于全维高精度势能面对H+SO2→OH+SO反应的动力学研究

H+SO₂→OH+SO反应在燃烧、大气和星际化学中都扮演着重要角色. 它还是具有深势阱中间体形成的典型反应,是检验速率理论和提供有趣反应动力学现象的理想候选反应. 基于之前构建的全维高精度势能面,本文对该反应进行了准经典动力学研究. 在1400 K≤T≤2200 K的温度范围内,计算值重现了实验速率常数. 当反应物SO₂处于振-转基态,在31.0∽40.0 kcal/mol的碰撞能范围内,计算得到的积分反应截面随碰撞能增加;在40∽55 kcal/mol的碰撞能范围内,积分反应截面几乎不受碰撞能影响. 产物角度分布呈现对称的前后向双峰结构. 本文还分析了产物OH和SO的振动态分布.     

基于最新势能面的C(1D)+H2反应的环聚合分子动力学研究

本文利用环聚合分子动力学方法对C(¹D)+H₂反应开展了详细的理论研究. 计算中使用了最近构建的Zhang-Ma-Bian(ZMB)从头算势能面,该势能面对锥形交叉附近区域以及范德华区域均有精确的描述. 环聚合分子动力学计算得到的热反应速率常数与最新实验值吻合很好. 与前人计算结果比较,发现在?¹A′电子基态的ZMB-a势能面上获得的反应速率常数远大于前人构建的RKHS势能面上的结果,这是由于ZMB势能面上的范德华鞍具有与之前势能面上的范德华阱完全不同的动态学作用,表明环聚合分子动力学方法能够处理范德华作用引起的势能面拓扑结构所导致的动态学效应. 本文还揭示了b¹A′′电子激发态ZMB-b势能面以及量子效应对反应的重要性.     

在惯性系中求解落体偏东问题

在惯性系中从机械能守恒出发求解了落体偏东问题并计算了下落过程中质点的引力势能与动能的转化.