He+H2+ HeH++H反应的从头算势能面"

运用多参考组态相互作用（MRCI）方法及aug-cc-pv5z基组构造了HeH2+体系基态（1 ）从头算势能面。采用多体展开法共拟合了8840个能量点,拟合误差为0.0677Kcal/mol。基于新势能面,采用准经典轨线方法计算了碰撞能在0-2 eV范围内的积分反应截面,计算结果与实验结果符合较好。     

用量子力学方法计算H+H2碰撞的反应几率

利用量子力学耦合通道扭曲波近似法(CCDWA)和三种势能面计算了H+H2碰撞的反应几率,结果发现在相同的势能面下利用CCDWA方法计算的反应几率和公认较好的计算结果符合很好,不同势能面共线势垒高度的差别将引起反应几率的不同.     

H+H2反应截面的全量子力学研究

利用全量子力学耦合通道扭曲波近似(CCDWA)方法和三种势能面计算了H+H2(vα=0,jα=0)→H2(vβ,jβ)+H碰撞的反应截面,并与准经典弹道计算结果及公认较好的计算结果作了比较.研究表明:在相同的势能面下利用CCDWA方法得到的截面和公认较好的计算结果符合很好,而准经典弹道计算的反应截面误差较大. 关键词： 反应散射 势能面 反应截面     

利用新构建势能面3A研究S(3P)+H2→HS+H反应含时波包量子散射

利用高精度从头算能量点构建三重态³A' 势能面. 单点能计算采用的方法是完全活性空间自洽场和多组态相互作用,计算中所用的基组是aug-cc-pV5Z,并利用多体展开方法结合AP函数形式对所有能量点进行了拟合. 在新构建的势能面基础上,在平动能0.8～2.2 eV进行了含时波包散射计算,计算中同时采用了离心突然近似和紧耦合方法. 另外还对初始振动态ν=0～3(j=0)情况下的总反应几率进行了计算,结果发现初始振动激发对该体系有很大的增强作用.     

Na2H的基态结构与势能函数

采用二次组态相关(QCISD)方法对NaH分子进行理论计算,得到它的几何结构、光谱性质,并拟合出它的Murrell-Sorbie势能函数.应用密度泛函(B3LYP)方法,在6-311G(3df,3pd)基组水平上对Na2H分子的基态结构进行优化,并用同样的基组对该分子进行了进一步的频率计算.结果发现Na2H分子的基态稳态结构为C2v构型,采用多体项展式理论导出了它的解析势能函数,其等值势能图准确再现了Na2H分子的结构特征和离解能.并报导了该分子对称伸缩振动等值势能图中存在的两个对称鞍点,对应于反应NaH Na→Na2H,活化能大约为14.56 kJ/mol.     

用神经网络方法拟合的反应体系H+CH4→H2+CH3的一个全域势能 (cited:12)

利用神经网络方法,基于47783个高精度从头算能量点构建了反应体系H+CH₄→H₂+CH₃的一个全域势能面．通过大量的准经典轨线以及量子散射计算测试了势能面的收敛性质．这个势能面对于拟合过程以及从头算点的数目都已经完全收敛,拟合误差很小且比Shepard插值的势能面计算速度更快,代表了此标志性多原子反应体系最好的势能面．     

H2S和H2S+的振动光谱和电子能谱计算

本文基于MOLPRO软件包使用从头算方法计算了星际分子H₂S及其阳离子H₂S⁺的势能面及光电子能谱.首先，在(U)CCSD/cc-pVQZ理论水平下获取了H₂S沿法线坐标展开的势能面，势能面直观描述了不同振动模式耦合对分子能量变化的影响，S-H键的非对称伸缩振动和面内弯曲振动共同作用使得系统势能变化范围明显变大.振动多组态相互作用方法被用来计算非谐振动频率和振动光谱，计算结果显示，倍频和组合频之间出现了强烈的费米共振，使得相应波段处的红外强度显著增强.最后，使用拉曼波函数和收缩不变Krylov子空间方法首次计算了H₂SΧ¹A₁→H₂S⁺ X²B₁的光电子能谱.此项研究有助于进一步理解星际分子的内部结构，并为实验研究及星际观测提供参考.     

反应体系F+H2→HF+H的一个全域从头算势能面

采用多参考组态相互作用方法和 aug-cc-pV5z 基组并虑及 Davidson 修正计算了反应体系 F H_2→HF H 的约15,000个对称性唯一的构型点的势能值,并利用三次样条插值发展了该体系的三维全域势能面.在反应入口通道,还考虑了旋轨耦合效应.     

MgB和MgB_2(~1A_1)的结构与解析势能函数

分别采用QCISD/6-311G和QCISD/6-311++G(df)方法,对MgB和MgB2分子的微观结构进行理论计算.在此计算基础上,运用多体展式理论方法,推导出MgB2分子的解析势能函数,其等值势能面图准确再现了MgB2分子的结构特征及势阱深度,并讨论了B+MgB和Mg+BB分子反应的势能面特征.这些结果可用于微观反应动力学的研究.     

神经网络方法构建苯硫酚三维非绝热势能面（英文）

本文利用神经网络方法构建了苯硫酚三维的基态(S0)绝热势能面以及激发态~1ππ~*和~1πσ~*耦合的非绝热势能面.特别地,~1ππ~*和~1πσ~*非绝热势能面的透热化是通过神经网络拟合实现的,拟合仅基于其绝热势能,但对非绝热势能矩阵的非对角元施加了对称性的正确限制.神经网络拟合的3个态的势能面的拟合均方根误差都非常小(4 meV),体现了神经网络方法的高精度.在神经网络构建的势能面上计算得到了S0态的低位能级和S_1态0~0的寿命,结果均与早前非绝热势能面上的结果吻合,验证了神经网络方法构建的绝热和非绝热势能面的高精度和可靠性.     

基于变分Monte Carno计算的H2分子势能曲线的可靠性

基于Alexander and Coldwell报道的变分Monte Carno计算的势能曲线,计算了H2分子的振转能级.利用Murrell-Sorbie函数和最小二乘法拟合出了解析势能函数,以此为基础计算出光谱常数.通过与实验比较发现:由于该势能曲线考虑了相对论效应和非绝热效应,振转能级和光谱常数计算结果与实验值符合良好,与目前最高水平从头计算的结果相当,说明这种计算方法准确可靠,可进一步推广应用于势能函数计算之中.     

NaB2和Na2B的结构与解析势能函数

使用Gaussian03程序包,采用单双取代耦合簇(CCSD)方法,选择基组6-311+g(2df)、6-311++g(3df,3pd)分别对B2及NaB和Na2分子的基态进行优化计算,运用最小二乘法拟合得到B2、NaB和Na2分子势能函数,给出与实验值符合很好的光谱常数;选择6-31++g(3df,3pd)基组,对NaB2和Na2B分子的结构进行优化计算.在此基础上,采用多体项展式法,得到NaB2和Na2B分子基态解析势能函数.势能面静态特征表明,该势能函数准确再现了NaB2和Na2B分子基态平衡结构.     

NaB2和Na2B的结构与解析势能函数

使用Gaussian03程序包,采用单双取代耦合簇(CCSD)方法,选择基组6-311+g(2df)、6-311++g(3df,3pd)分别对B2及NaB和Na2分子的基态进行优化计算,运用最小二乘法拟合得到B2、NaB和Na2分子势能函数,给出与实验值符合很好的光谱常数;选择6-31++g（3df,3pd）基组,对NaB2和Na2B分子的结构进行优化计算.在此基础上,采用多体项展式法,得到NaB2和Na2B分子基态解析势能函数.势能面静态特征表明,该势能函数准确再现了NaB2和Na2B分子基态平衡结构.     

基态Li_n(n=2,3)结构与完全解析势能函数

使用Gaussian03程序包,采用全电子单双取代耦合簇(CCSD(full))方法,选择基组6-311+g(2df),对Li2分子的基态进行优化计算,采用十一参量Murrell-Sorbie函数,运用最小二乘法拟合得到Li2分子基态解析势能函数,给出与实验值符合很好的光谱常数;使用同样的方法和基组,对Li3分子的基态结构进行优化计算,得到Li3分子基态平衡结构.采用多体项展式法,利用Li3分子平衡结构C2v的几何参数、力常数和离解能,以及七个线性系数Ci(i=1,2,3,4,5,6,7)与两个非线性系数的函数关系,进行非线性优化拟合得到两个非线性系数,进而得到七个线性系数,得到Li3分子基态完全解析势能函数.势能面静态特征表明,该势能函数再现了Li3分子基态全部平衡结构特征.     

基态Bn(n=2,3)结构与完全解析势能函数简

使用Gaussian03程序包,采用全电子单双取代耦合簇(CCSD(full))方法,选择基组6-311+g(2df),对B2分子的基态进行优化计算,采用十一参量Murrell-Sorbie函数,运用最小二乘法拟合得到B2分子基态解析势能函数,给出与实验值符合很好的光谱常数;使用同样的方法和基组,对B3分子的基态结构进行优化计算,得到B3分子基态平衡结构.采用多体项展式法,利用B3分子平衡结构C2v的几何参数、力常数和离解能,以及七个线性系数Ci(i=1,2,3,4,5,6,7)与两个非线性系数的函数关系,进行非线性优化拟合得到两个非线性系数,进而得到七个线性系数,得到B3分子基态完全解析势能函数.势能面静态特征表明,该势能函数准确再现了B3分子基态全部平衡结构特征.     

Pb2,PdPb2分子的势能函数

采用Gaussian 98程序,运用B3LYP方法,对Pd和Pb原子采用收缩价基组LANL2DZ,对Pb2和PdPb2分子的微观结构进行了理论计算.由于Pb2分子离解后一个Pb原子处于基态,另一个Pb原子处于激发态,采用最小二乘法拟合Pb2分子的势能函数.选用的函数形式为Murrell-Sorbie势能函数加上开关函数.使用多体展式理论导出了势函数中的参数进而给出PdPb2分子基态势函数的解析表达式,其势能面准确地复现了PdPb2分子的两个稳定构型(C2V和C∞v)及其能量关系.     

Mg2Ni分子的结构和势能函数

采用密度泛函B3LYP方法在6-311+g(d,p)基组水平上对MgNi、Mg2和Mg2Ni分子的各种可能的结构进行优化,得到了它们的几何构型、平衡核间距、离解能和谐振频率.采用最小二乘法拟合MgNi和Mg2分子的Murrell-Sorbie势能函数,在此基础上导出光谱数据和力常数.通过多体项展式理论导出Mg2Ni分子的解析势能函数,并绘出了Mg2Ni分子的等值势能面,其等值势能面正确地反应了基态Mg2Ni分子的平衡构型特征.     

C (~3P)+NO (X~2Ⅱ)→CN (X~2∑+)+O (~3P) 反应的矢量相关研究

采用准经典轨线方法研究了三个不同碰撞能下,反应C(~3P)+NO (X~2∏)→CN (X~2∑+)+O(~3P) 在Nyman等拟合的两个势能面(~2A″和~2A′)上反应物与产物之间的矢量相关.结果表明在两个势能面上产物的转动角动量矢量不仅取向,而且沿着y轴垂直于散射平面定向,产物的角动量定向依赖于碰撞能和势能面.同时,计算的四个重整化极化依赖微分截面显示产物在三个散射角被强烈极化,这与在两个势能面上的三个过渡态有关.     

O2分子激发态a1△g和b1∑+g态的分析势能函数

以aug-co-pVQZ,cc-pV5Z,6-311++g(d,p)和6-311++g(3df,3pd)等为基函数,采用多组态相互作用(MRCI)方法对O2分子最低的两个激发态1△g和1∑+g的平衡结构、谐振频率和势能曲线进行了计算.并选用MurrelI-Sorbie势能函数对曲线进行拟合,利用拟合的参数值计算出了力常数和光谱数据.结果表明计算值与实验值符合较好.     

基态Bn(n=2, 3)结构与完全解析势能函数

使用Gaussian03程序包, 采用全电子单双取代耦合簇(CCSD(full))方法, 选择基组6-311+g(2df) , 对B2分子的基态进行优化计算, 采用十一参量Murrell-Sorbie函数, 运用最小二乘法拟合得到B2分子基态解析势能函数, 给出与实验值符合很好的光谱常数; 使用同样的方法和基组, 对B3分子的基态结构进行优化计算, 得到B3分子基态平衡结构. 采用多体项展式法, 利用B3分子平衡结构C2v的几何参数、力常数和离解能, 以及七个线性系数Ci(i=1, 2, 3, 4, 5, 6, 7)与两个非线性系数的函数关系, 进行非线性优化拟合得到两个非线性系数, 进而得到七个线性系数, 得到B3分子基态完全解析势能函数. 势能面静态特征表明, 该势能函数准确再现了B3分子基态全部平衡结构特征.