Li2 K2双原子分子基态与激发态势能曲线的理论计算

本文采用STO—5G基组和闭壳层Hartree—Fock—Roothaan(HFR)从头计算法计算了Li_2,K_2双原子分子基态及许多激发态的分子势能曲线。总能量的自洽精度达到2×10~(-6)。计算结果与用其它方法计算的结果符合得很好     

用能量自洽法研究异核双原子分子的势能曲线

用研究双原子分子解析势能函数的新方法——ECM方法进一步研究了一些异核双原子分子的电子基态和激发态：CIF分子的A３Π1和B0＋(３Π)态，CH分子的X２Π态，BH分子的X1Σ＋态，XeO分子的d1Σ＋态，LaF分子的X1Σ＋态，Li7D分子的X1Σ＋态，NaRb分子的X1Σ＋态，KRb分子的(2)3Σ＋和21Π态等.获得的势能曲线表明，ECM势能很好地符合Rydberg-Klein-Rees(RKR)值，得到了比Morse势，huxley-Murrell-Sorbie(HMS)势更令人满意的结果.而且在分子渐近区和分子离解区域，ECM方法还能得到RKR可能缺乏的可靠物理数据. 关键词： 能量自洽 双原子分子 势能函数 电子态     

低激发态双原子分子的势能函数

对比分析了双原子分子的势能函数,发现用PG—高斯势函数表述低激发态双原子分子是一种可能的选择。按照微观过程原理和本文提出的“激发态选择规则”,可确定低激发态的离解极限与离解能,并导出63个激发态的PG—高斯势函数。     

一种获取双原子分子的转动常数的方法

本文给出了对于任何双原子分子或双原子自由基分子,由其振动常数结合RKR方程、结合从头计算,获得势能曲线及其转动常数的方法.并对CN的基态和激发态进行了计算,结果与实验值符合得非常一致.A     

双原子分子离子XY+部分电子态完全振动能谱的精确研究

关于双原子分子离子XY 的完全振动能谱,目前还没有实验和理论的数据报道。本文首次应用代数方法AM(Algebraic Method),获得了BeH -X1Σ 态,CO -X2Σ 态,F2 -X2Πg态,O2 -A2Πu态和Li2 -X2Σg 态的精确振动光谱常数和完全振动能谱,解决了实验方法和精确量子力学理论方法难以获得双原子分子离子XY 的包含最高振动能级在内的所有高阶振动能级的精确数值这一问题。所有研究结果表明:由部分较低的实验精确振动能级,可用AM产生双原子分子离子XY 的精确振动光谱常数和包含全部激发态的完全振动能谱;所得的AM振动能谱比其他理论方法得到的结果更好。     

双原子分子轨道能级的台劳展开与应用

在双原子分子轨道能级台劳展开的基础上,给出了2~n个同核原子体系的分子轨道能级的解析表达式,该式仅涉及双中心积分,对简化量子化学计算有一定价值。     

双原子分子离解能的精确研究

基于LeRoy与Bemstein的工作，孙卫国等最近建立了计算精确的双原子分子离解能的新解析表达式．应用该公式和最近建立的研究双原子分子精确振动能谱的代数方法（AM），我们研究了一些双原子分子部分电子态的分子离解能，并与实验值进行了比较．研究结果表明，用新解析式获得的精确分子离解能与实验值符合得非常好．该式在理论上提供了获得精确分子离解能的物理新方法．     

普适性双原子分子解析势能函数的研究

提出了一种构造解析势能函数的新方法,由此得到了一种既适用于中性双原子分子又适用于带电双原子分子离子的解析势能函数。本文用八种基本类型的双原子分子——同核中性基态双原子分子Na₂-X1Σ+_g,同核中性激发态双原子分子C₂-A1Π_u,同核带电基态双原子分子离子He+₂-X2Σ+_u,同核带电激发态双原子分子离子N+₂-B2Σ+_u,异核中性基态双原子分子NaLi-X1Σ+_g,异核中性激发态双原子分子BH-B1Σ+,异核带电基态双原子分子离子(BC)--X3Π,异核带电激发态双原子分子离子(CS)+-A2Π等共21个算例对势能函数进行了验证并与RKR (Rydberg-Klein-Rees)实验数据进行了比较,计算结果与RKR数据符合很好。     

六参数高精度双原子分子解析势能函数

提出一种构造解析势能函数的新方法,得到一种六参数解析势能函数,该势能函数适用于多种基本类型的双原子分子.用同核中性基态双原子分子Li₂-X¹Σ_g⁺、Na₂-X¹Σ_g⁺,异核带电基态双原子分子离子(BC)^--X³∏,异核中性激发态双原子分子PbS-A¹Σ⁺、BaO-A¹Σ,异核带电激发态双原子分子离子(CS)⁺-A²∏,同核中性激发态双原子分子K₂-B¹∏_u,同核带电激发态双原子分子离子N₂⁺-B²Σ_u⁺等共36个算例对该势能函数进行验证,计算结果与RKR(Rydberg-Klein-Rees)数据或振动能谱实验数据高精度符合.同时发现,在分子振动能谱计算精度方面,该势能函数总体上优于Murrell-Sorbie势能函数.     

双原子分子离解能的精确研究

基于LeRoy与Bernstein的工作,孙卫国等最近建立了计算精确的双原子分子离解能的新解析表达式.应用该公式和最近建立的研究双原子分子精确振动能谱的代数方法(AM),我们研究了一些双原子分子部分电子态的分子离解能,并与实验值进行了比较.研究结果表明,用新解析式获得的精确分子离解能与实验值符合得非常好.该式在理论上提供了获得精确分子离解能的物理新方法.     

OD、OT、DT分子基态的结构和解析势能函数

本文采用密度泛函方法(B3LYP)和组态相关方法(QCISD(T))优化计算了OD、OT、DT分子基态(X2Ⅱ)的平衡结构、振动频率和离解能.根据原子分子反应静力学原理.导出了OD、OT、DT分子基态(X2 Ⅱ)的合理离解极限,采用最小二乘法拟合Murrell-Sorbie函数得到了相应的势能函数和与该基态相对应的光谱常数(Be,ae,we和wexe),计算结果与实验数据符合得相当好.     

Structure and analytical potential energy function for the ground state of the BCx (x=0, －1)

In this paper, the electronic states of the ground states and dissociation limits of BC and BC- are correctly determined based on group theory and atomic and molecular reaction statics. The equilibrium geometries, harmonic frequencies and dissociation energies of the ground state of BC and BC- are calculated by using density function theory and quadratic CI method including single and double substitutions. The analytical potential energy functions of these states have been fitted with Murrell-Sorbie potential energy function from our ab initio calculation results. The spectroscopic data （αe, ωe and ωeχe） of each state is calculated via the relation between analytical potential energy function and spectroscopic data. All the calculations are in good agreement with the experimental data.     

Ab initio calculation of accurate dissociation energy, potential energy curve and dipole moment function for the A1∑ + state 7LiH molecule

The reasonable dissociation limit of the A¹∑⁺ state $^{7}$LiH molecule is obtained. The accurate dissociation energy and the equilibrium geometry of this state are calculated using a symmetry-adapted-cluster configuration-interaction method in complete active space for the first time. The whole potential energy curve and the dipole moment function for theA¹∑⁺ state are calculated over a wide internuclear separation range from about 0.1 to 1.4\,nm. The calculated equilibrium geometry and dissociation energy of this potential energy curve are of R_{\e}=0.2487\,nm and D_{\e}=1.064\,eV, respectively. The unusual negative values of the anharmonicity constant and the vibration-rotational coupling constant are of \textit{\omega }_{\e}\textit{\chi }_{\e}=--4.7158cm^{ - 1} and \textit{\alpha }_{\e}=--0.08649cm^{ -1}, respectively. The vertical excitation energy from the ground to the A¹∑⁺ state is calculated and the value is of 3.613\,eV at 0.15875nm (the equilibrium position of the ground state). The highly anomalous shape of this potential energy curve, which is exceptionally flat over a wide radial range around the equilibrium position, is discussed in detail. The harmonic frequency value of 502.47cm¹ about this state is approximately estimated. Careful comparison of the theoretical determinations with those obtained by previous theories about the A¹∑⁺ state dissociation energy clearly shows that the present calculations are much closer to the experiments than previous theories, thus represents an improvement.     

Ab initio calculations of accurate dissociation energy and analytic potential energy function for the second excited state B1∏ of 7LiH

The reasonable dissociation limit of the second excited singlet state B¹∏ of ⁷LiH molecule is obtained. The accurate dissociation energy and equilibrium geometry of the B¹\Pi state are calculated using a symmetry-adapted-cluster configuration--interaction method in full active space. The whole potential energy curve for the B¹∏ state is obtained over the internuclear distance ranging from about 0.10nm to 0.54nm, and has a least-square fit to the analytic Murrell--Sorbie function form. The vertical excitation energy is calculated from the ground state to the B¹∏ state and compared with previous theoretical results. The equilibrium internuclear distance obtained by geometry optimization is found to be quite different from that obtained by single-point energy scanning under the same calculation condition. Based on the analytic potential energy function, the harmonic frequency value of the B¹∏ state is estimated. A comparison of the theoretical calculations of dissociation energies, equilibrium interatomic distances and the analytic potential energy function with those obtained by previous theoretical results clearly shows that the present work is more comprehensive and in better agreement with experiments than previous theories, thus it is an improvement on previous theories.     

完全活性空间自洽场方法对NO分子基态势能函数的研究

利用MOLPRO从头算程序包, 选用完全活性空间自洽场(CASSCF)方法并选取cc-pVTZ, cc-pVQZ, cc-pV5Z和cc-pV6Z基组, 计算了NO分子基态的平衡核间距与谐振频率, 从中优选出cc-pV5Z基组进行单点能扫描, 并将扫描结果用最小二乘法拟合得到4参数、5参数、6参数和7参数的Murrell-Sorbie势能函数. 通过比较由势能函数计算出的光谱数据, 最终确定6参数的结果最好. 最后, 利用拟合出的解析势能函数, 通过求解径向Schrödinger方程, 得到了NO分子基态J=0时的全部38个振动态的振动能级、经典转折点、惯性转动常数以及6个离心畸变常数.     

GaN和GaN2分子基态的结构和解析势能函数

本文利用密度泛函方法,优化了CaN2基态(X2B2)的几何平衡结构,其对称性属于C2v,它的键角、平衡核间距和离解能分别为26.428°,0.2526 nm和9.574 eV.基于量子化学计算得出GaN和GaN2基电子态及其离解极限之后,利用计算数据和最小二乘法得到了GaN的M-S解析势能函数,并计算出各态的谐振频率、力常数和光谱数据.利用多体项展式理论首次得到了GaN2基态的解析势能函数,它准确表达了GaN2的平衡几何结构,在此基础上讨论了Ga N2和N GaN反应系统的反应动力学.     

7LiH分子基态及若干激发态的平衡几何及垂直激发能的从头算

使用“对称性匹配簇-组态相互作用方法”(SAC/SAC-CI),在多种基组下计算了7LiH分子X1Σ 、A1Σ 、B1Π及b3Π态的平衡几何,并将由“几何优化”得到的相应各态的平衡几何与“单点能扫描”得到的结果进行了比较。比较的结果表明,由“几何优化”得到的相应各态的平衡几何,与“单点能扫描”得到的结果存在着差异。分析的结果表明,在完全活性空间中,由“SAC/SAC-CI”方法进行“单点能扫描”得到的结果,比“几何优化”得到的结果更加合理。首次报告了使用6-311G(3df,3pd)基组进行“单点等扫描”得到的相应各态的平衡几何。其值分别是:基态(X1Σ )为0.1588 nm、单重态的第一激发态(A1Σ )为0.2487 nm、单重态的第二激发态(B1Π)为0.2434 nm、三重态的第二激发态(b3Π)为0.1958 nm。这一计算结果与实验值非常接近。还研究了从基态到上述相应各态的垂直激发能,在基态的平衡位置处,其值分别为(A1Σ ←X1Σ )3.613 eV(、B1Π←X1Σ )4.612 eV和(b3Π←X1Σ )4.233 eV。与其它理论计算结果进行比较后得出,本文的计算结果非常接近于使用很复杂的计算方法获得的结果。     

Accurate double many-body expansion potential energy surface of HS_2(A~2A') by scaling the external correlation

A globally accurate single-sheeted double many-body expansion potential energy surface is reported for the first excited state of HS_2 by fitting the accurate ab initio energies, which are calculated at the multireference configuration interaction level with the aug-cc-pVQZ basis set. By using the double many-body expansion-scaled external correlation method,such calculated ab initio energies are then slightly corrected by scaling their dynamical correlation. A grid of 2767 ab initio energies is used in the least-square fitting procedure with the total root-mean square deviation being 1.406 kcal · mol~(-1).The topographical features of the HS_2(A_2A') global potential energy surface are examined in detail. The attributes of the stationary points are presented and compared with the corresponding ab initio results as well as experimental and other theoretical data, showing good agreement. The resulting potential energy surface of HS_2(A_2A') can be used as a building block for constructing the global potential energy surfaces of larger S/H molecular systems and recommended for dynamic studies on the title molecular system.     

CO基态的偶极矩的量化计算研究

运用MOLPRO软件,采用基组6-31G,6-31G(d),6-31G(2d,3f),aug-cc-pv5z,方法HF(Hartree-Fock)、CASSCF(完全活性空间自洽场)和MRCI(多参考组态相互作用),计算了CO的基态X1Σ+的偶极矩随核间距变化的12条曲线,与基于实验光谱数据的偶极矩在平衡位置的值以及在平衡位置的斜率进行了比较,结果表明,采用MRCI方法,aug-cc-pv5z基组的偶极矩曲线最为准确.平衡位置斜率的偏差小于1%.接着采用MRCI方法,aug-cc-pv5z基组计算得到了CO基态的势能曲线,进而通过求解核的薛定谔方程得到振动能级以及分子力常数,分子力常数与实验值符合的很好.     

用耦合簇理论及系列相关一致基研究基态SH+ 和 SD+ 的光谱常数*

利用耦合簇理论CCSD(T)和相关一致基组cc-pVXZ和aug-co-pVXZ(X=2,3,4,5)计算了SH 和SD 的基态平衡几何re、谐振频率ωe和离解能De.计算结果表明,所用基组越大,得到的结果与实验值之间的差别就越小,因此我们选用了大基组aug-cc-pV5Z在0.07～2.50 nm的核间距范围内对SH 和SD 的基态进行单点能计算,并将计算结果拟合成了Murrell-Sorbie函数.利用得到的解析势能函数,计算了SH 和SD 的其余3个光谱常数(ωexe,αe和Be),其结果与实验值符合得相当好.