OH，OCI，HOCl（1A＇）的从头算与势能曲线

采用Gassian09程序包中的多种方法对OH,OCI,HOCl分子的基态结构进行优化计算,优选出QCISD／6—311G（2df）,B3P86／6—311＋G（2df）方法分别对OH（X2II）,OCI（X2Ⅱ）分子进行计算,得到平衡核间距RoH=0．09696nm,Rocl=0．1569nm,谐振频率w（OH）：3745．37cm-1,w（OCI）=892．046cm-1,与实验结果非常符合．用Murrell—Sorbie势能函数对OH和OCI分子的扫描势能点进行拟合,其扫描点都与四参数Murrell—Sorbie函数拟合曲线符合得很好．优选出QCISD（T）／D95（df,pd）方法对HOCl分子进行计算,得到基态为x1A ＇,键长ROH=0．0966nm,键角∠HOCl=102．3°,谐振频率w1（a1）：738．69cm-1,w2（b2）=1260．25cm-1,离解能De=2．24eV．通过比较发现这些结果与实验值符合得很好,并优于文献报道的结果．随后计算出了力常数,在此基础上,推导出HOCl分子的多体展式势能函数．报道了HOCl分子对称伸缩振动势能图中在H＋OCl→HOCl反应通道上有一鞍点,H原子需要越过1．74eV的能垒才能生成HOCl的稳定结构,在Cl＋OH→HOCl通道上不存在明显势垒,容易形成稳定的HOCl分子．     

OH,OCl,HOCl(~1A′)的从头算与势能曲线

采用Gassian09程序包中的多种方法对OH,OCl,HOCl分子的基态结构进行优化计算,优选出QCISD/6-311G（2df）,B3P86/6-31 1+G（2df）方法分别对OH（X²∏）,OCl（X²∏）分子进行计算,得到平衡核间距R_OH=0.09696nm,R_OCl=0.1569 nm,谐振频率w（OH）=3745.37 cm^-1,e（OCl）=892.046 cm^-1,与实验结果非常符合.用Murrell-Sorbie势能函数对OH和OCl分子的扫描势能点进行拟合,其扫描点都与四参数Murrell-Sorbie函数拟合曲线符合得很好.优选出QCISD（T）/D95（df,pd）方法对HOCl分子进行计算,得到基态为X¹A’,键长R_OH=0.0966 nm,键角∠HOCl=102.3°,谐振频率w₁（a₁）=738.69 cm^-1,w₂（b2）=1260.25 cm^-1,离解能D_e=2.24 eV.通过比较发现这些结果与实验值符合得很好,并优于文献报道的结果.随后计算出了力常数,在此基础上,推导出HOCl分子的多体展式势能函数.报道了HOCl分子对称伸缩振动势能图中在H+OCl→HOCl反应通道上有一鞍点,H原子需要越过1.74 eV的能垒才能生成HOCl的稳定结构,在Cl+OH→HOCl通道上不存在明显势垒,容易形成稳定的HOCl分子.     

HNO(1A’)自由基的从头算势能曲线

采用从头算的偶合族理论和组态相互作用方法对NH, NO, HNO 自由基的基态结构进行研究, 借助多体项展式理论导出HNO 自由基的势能函数并绘制了等值势能图. 首次报道了HNO 自由基对称伸缩振动和旋转势能图中, 在O+NH→HNO, H+NO→HNO, N+HO→HNO 反应通道上都有鞍点出现, O 原子需要越过1.153 eV 的能垒, H 原子须克服1.683 eV 的能量, N 原子须克服2.126 eV 的能量才能生成稳定的HNO 自由基. 同时也首次报道了HNO(X¹A’)自由基的同分异构体HON(X³A") 在势能曲线中的位置及HNO↔HON 转变所需的能量. 关键词： HNO 从头算 势能曲线     

SeOx(x=1,2)的从头算势能曲线和光谱常数

采用从头算的多种方法和基组优化计算SeO_x(x=1,2)自由基的基态结构、谐振频率及离解能,优选出QCISD(T)/6-311+G(2df)、B3LYP/6-311G(3d2f)方法分别对SeO、SeO₂自由基进行计算,计算结果与实验结果吻合很好.对SeO自由基拟合出Murrell-Sorbie势能函数参数,计算出SeO自由基的光谱常数和力常数.计算出SeO₂自由基力常数,导出SeO₂自由基的多体展式势能函数,发现SeO₂自由基对称伸缩振动势能图中在对称的O+SeO→SeO₂反应通道上有一鞍点,其活化能约为48.24 kJ·mol^-1,O原子需要越过0.5 eV的能垒才能生成SeO₂的稳定结构.     

OH分子基态(X~2Π)的结构与势能函数

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

OH分子基态(X2Ⅱ)的结构与势能函数

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

多参考组态相互作用方法研究ZnHg低激发态(1∏，3∏)的势能曲线和解析势能函数

采用多参考组态相互作用方法计算了ZnHg二聚体两个低激发∏态(¹∏，³∏)的原子间相互作用势能曲线. 用Murrel-Sorbie函数拟合得到了相应的解析势能函数，并用其计算力常数，进而确定了光谱常数. 所得结果与仅有的理论工作进行了比较. 基于所得势能曲线，通过解分子中原子核运动的薛定谔方程预测了各电子态的振动能级. 关键词： 势能曲线 解析势能函数 光谱常数 振动能级     

OH分子基态（X2Π）的结构与势能函数

采用密度泛函理论的B3LYP方法和二次组态相互作用(QCISD(T))方法优化计算了OH分子基态(X²Π)的平衡结构、振动频率和离解能.根据原子分子反应静力学原理，导出了OH分子基态(X²Π)的合理离解极限，采用最小二乘法拟合Murrell-Sorbie函数得到了相应的势能函数和与该基态相对应的光谱常数(B_e，α_e，ω_e和ω_eχ_e)，计算结果与实验数据符合得相当好. 关键词： OH分子 2Π)')" href="#">基态(X²Π) 势能函数     

B2C(1A1)和BC2(2A′)的结构与解析势能函数

采用单双取代的二次组态相互作用方法，分别选用6-311++G(d,p)和6-311G(df,pd)基组，对B₂C和BC₂分子的结构进行了优化，得到这两个分子的基态结构为C_2v和C_s，基态电子状态为¹A₁和²A′，同时还得到了它们的平衡几何结构、离解能、谐振频率和力常数. 关键词： 碳化硼 Murrell-Sorbie函数 谐振频率 势能函数     

POX(X=1,2)的光谱常数与从头算势能曲线

采用从头算的多种方法对PO、PO₂自由基的基态结构进行优化计算,结果表明：使用密度泛函（DFT）方法计算的结果最接近实验值.对PO双原子分子优选6-311G（3df）基组进行计算、扫描并拟合.对PO₂三原子分子优选出6-311+G（3df）方法计算结构参数、谐振频率、离解能及力常数,借助多体项展式理论导出PO₂自由基的势能函数并绘制等值势能图.发现：PO₂自由基的对称伸缩振动和旋转势能图中,在O+PO→OPO反应通道上都有鞍点出现,O原子需要越过0.55 eV的能量才能生成稳定的PO₂自由基.要形成PO₂自由基只能通过两等价的通道越过势垒才能形成.     

SiF$lt;sub$gt;2$lt;/sub$gt;($lt;sup$gt;1$lt;/sup$gt;A$lt;sub$gt;1$lt;/sub$gt;)自由基的从头算及势能函数

采用从头算的单双取代的二次组态相互作用方法及耦合簇理论对SiF₂自由基的基态进行结构优化, 发现用单双取代的二次组态相互作用方法配有基组6-311G(2df)计算得到的结构参数、谐振频率、离解能及力常数与实验值最接近并优于文献值. 借助多体项展式理论导出SiF₂自由基的势能函数并绘制了等值势能图. SiF₂自由基对称伸缩振动和旋转势能图显示: 在SiF+F →SiF₂反应通道上有鞍点出现, F原子需要越过4.38 eV的能垒才能生成稳定的SiF₂自由基; 只能通过SiF+F→SiF₂两个等价的通道才能生成稳定的SiF₂ 自由基, 并且该反应是有阈能的吸热反应. 关键词： 2')" href="#">SiF₂ 从头算 势能曲线     

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.     

The structure and the analytical potential energy function of NH2 （X^2B1）

This paper reports that the equilibrium structure of NH2 has been optimized at the QCISD/6-311＋＋G （3df, 3pd） level. The ground-state NH2 has a bent （C2v, X^2B1） structure with an angle of 103.0582°. The geometrical structure is in good agreement with the other calculational and experimental results. The harmonic frequencies and the force constants have also been calculated. Based on the group theory and the principle of microscopic reversibility, the dissociation limits of NH2（C2v, X^2B1） have been derived. The potential energy surface of NH2（X^2B1） is reasonable. The contour lines are constructed, the structure and energy of NH2 reappear on the potential energy surface.     

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.     

7Li2分子2 3Σ+g激发态的解析势能函数、谐振频率及振动能级

使用Gaussian03程序包中的“对称性匹配簇/对称性匹配簇-组态相互作用”方法, 利用多个基组对⁷Li₂分子2³Σ⁺_g态的平衡几何进行了优化计算. 同时, 在优化得到的平衡位置附近、于同一条件下通过精细的单点能扫描, 获得了相应基组下的平衡核间距. 发现两者的结果不一致，对不一致的原因进行了解释. 分析表明，由单点能扫描得到的平衡核间距应更为合理. 同时也得出了6-311++G(3df,3pd)，6-311++G(2df,2pd) 及6-311++G(2df,pd)基组均为较优基组的结论.于2.5a₀—37a₀的范围内利用6-311++G(3df,3pd)基组进行单点能扫描并使用最小二乘法拟合出了该态的解析势能函数. 利用解析势能函数的物理意义并结合Rydberg-Klein-Rees方法, 计算出了其相应的谐振频率, 进而计算了其他光谱常数. 为便于比较和分析, 对基态也进行了相应的计算.利用得到的解析势能函数, 对2³Σ⁺_g态的振动能级及振动经典转折点也进行了计算. 关键词： 势能函数 谐振频率 振动能级 Murrell-Sorbie函数     

Spin--orbit ab initio curves of 80Se2+ ion and theassignment of photoelectron spectra of 80Se2 molecule

This paper carries out ab initio calculations to study the ^80Se2（X^3Σg^-） state and ^80Se2^＋（X^2Πg）, ^80Se2^＋（a^4Πg） states by using completed active space self-consistent field and multi-reference second order perturbation theory. The electronic curves of these states including spin-orbit coupling are calculated, and then the spectroscopic parameters are obtained. The photoelectron spectra of ^80Se2 molecule in gas phase are assigned according to Franck-Condon analysis based on calculated potential energy curves. The ionization energies of ^80Se2 molecule are determined by the present calculation.     

The analytical potential energy function of flue gas SO2（X1A1）

The equilibrium structure of flue gas SO2 is optimized using the density functional theory （DFT）/B3P86 method and CC-PV5Z basis. The result shows that it has a bent （C2v, X1A1） ground state structure with an angle of 119.1184°. The vibronic frequencies and the force constants are also calculated. Based on the principles of atomic and molecular reaction statics （AMIIS）, the possible electronic states and reasonable dissociation limits for the ground state of SO2 molecule are determined. The potential functions of SO and 02 are fitted by the modified Murrell-Sorbie＋c6 （M-S＋c6） potential function and the fitted parameters, the force constants and the spectroscopic constants are obtained, which are all close to the experimental values. The analytic potential energy function of the SO2 （X1A1） molecule is derived using the many-body expansion theory. The contour liues are constructed, which show the static properties of SO2 （XIA1）, such as the equilibrium structure, the lowest energies, the most possible reaction channel, etc.     

The theoretical study on the potential energy curves for X^1 ∑^＋, A^1 ⅡI and C^1 ∑^- states of SiO molecule

This paper applies the symmetry-aziapted-cluster/symmetry-adapted-cluster configuration-interaction （SAC/SACCI） method to optimize the structures for X^1∑^＋, A^1 Ⅱ and C^1 ∑^- states of SiO molecule with the basis sets D95＋＋, 6-311＋＋G and 6-311＋＋G^＊＊. Comparing the obtained results with the experiments, it gets the conclusion that the basis set 6-311＋＋G^＊＊ is most suitable for the optimal structure calculations of X^1.∑^＋, A^Ⅱ and C^1∑^- states of SiO molecule. The whole potential energy curves for these electronic states are further scanned by using SAC/6-311＋＋G^＊＊ method for the ground state and SAC-CI/6-311＋＋G^＊＊ method for the excited states, then use a least square method to fit Murrell~Sorbie functions, at last the spectroscopic constants and force constants are calculated, which are in good agreement with the experimental data.