Configuration interaction studies on the spectroscopic properties of PbO including spin-orbit coupling

Lead oxide(Pb O), which plays the key roles in a range of research fields, has received a great deal of attention. Owing to the large density of electronic states and heavy atom Pb including in Pb O, the excited states of the molecule have not been well studied. In this work, high level multireference configuration interaction calculations on the low-lying states of Pb O have been carried out by utilizing the relativistic effective core potential. The effects of the core-valence correlation correction, the Davidson modification, and the spin–orbital coupling on the electronic structure of the Pb O molecule are estimated. The potential energy curves of 18 Λ-S states correlated to the lowest dissociation limit(Pb(~3P_g) + O(~3P_g)) are reported. The calculated spectroscopic parameters of the electronic states below 30000 cm~(-1), for instance, X~1Σ~+, 1~3Σ~+,and 1~3Σ~-, and their spin–orbit coupling interaction, are compared with the experimental results, and good agreements are derived. The dipole moments of the 18 Λ-S states are computed with the configuration interaction method, and the calculated dipole moments of X~1Σ~+and 1~3Σ~+are consistent with the previous experimental results. The transition dipole moments from 1~1Π, 2~1Π, and 2~Σ to X~1Σ~+and other singlet excited states are estimated. The radiative lifetime of several low-lying vibrational levels of 1~1Π, 2~1Π, and 2~1Σ~+ states are evaluated.     

一氟化碳电子态的光谱性质和预解离机理的理论研究

采用考虑Davidson修正的内收缩多参考组态相互作用(icMRCI+Q)方法结合相关一致基组aug-cc-pV5Z和aug-cc-pV6Z首次计算了一氟化碳(CF)11个Λ-S 态(X²Π , a⁴Σ^-, A²Σ⁺, B²Δ, 1⁴Π, 1²Σ^-, 2⁴Π, 1^{4}Δ , 1⁴Σ⁺, 2²Σ^-和2⁴Σ^-) 所产生的25个Ω 态的势能曲线. 计算中考虑了旋轨耦合效应、核价相关和标量相对论修正以及将参考能和相关能分别外推至完全基组极限. 基于得到的势能曲线, 获得了束缚和准束缚的Λ-S态和Ω 态的光谱常数, 与已有的实验结果非常符合. 分析了束缚和准束缚的Λ-S态在各自平衡核间距R_e处的主要电子组态. 由于1⁴Π 和2⁴Π态的避免交叉, 发现准束缚态2⁴Π. 由Λ-S态势能曲线的交叉现象, 借助于计算的旋轨耦合矩阵元, 分析了a⁴Σ^-和B²Δ 态的预解离机理. 计算了25个Ω 态的离解关系, 给出了它们的离解极限. 最后研究了A²Σ⁺-X²Π 跃迁特性, 本文计算得到的A²Σ⁺-X²Π跃迁的Frank-Condon 因子和辐射寿命与已有实验值也符合得非常好.     

icMRCI+Q理论研究BF~+离子电子态的光谱性质和预解离机理

采用考虑Davidson修正的内收缩多参考组态相互作用(icMRCI+Q)方法,结合相关一致基组aug-ccpV5Z和aug-cc-pV6Z,计算了BF+离子前两个离解极限B~+(~1S_g)+(~2P_u)和B+(~3P_u)+F(~2P_u)对应的14个Λ-S态(X~2Σ~+,1~2Π,2~2Π,2~2Σ~+,1~4Σ~+,1~4△,1~4Σ~-,1~2△,1~2Σ~-,3~2Σ~+,1~4Π,2~4Π,2~4Σ~+和3~2Π)和30个Ω态的势能曲线.在势能曲线的计算中,考虑了旋轨耦合效应、核价相关和标量相对论修正以及将参考能和相关能分别外推至完全基组极限.基于得到的势能曲线,获得了束缚和准束缚的12个A-S态和28个Ω态的光谱常数,并且X~2Σ~+态的光谱常数与已有的实验结果符合.此外,计算了BF分子X~1Σ~+态到BF+离子X~2Σ~+,1~2Π和2~2Σ~+态的垂直电离势和绝热电离势,并且BF~+(X~2Σ~+)←BF(X~1Σ~+)的垂直电离势和绝热电离势与相应的实验结果非常符合.由X~2Σ~+,2~2Π,1~4Σ~+,3~2Σ~+和3~2Π态和其他的激发A-S态势能曲线的交叉现象,借助于计算的旋轨耦合矩阵元,首次分析了X~2Σ~+和3~2Π态的预解离机理以及2~2Π(v′≥9),1~4Σ~+(v′≥4)和3~2Σ~+(v′≥4)的振动能级受到其他电子态的微扰.计算了30个Ω态离解极限处的相对能量,并且与实验结果十分符合.最后计算了2~2Π(v′=0—9)—X~2Σ~+,2~2Σ~+(v′=0—2)—X~2Σ~+,(3)1/2—(1)1/2~(势阱一)和(2)3/2(v′=0—9)—(1)1/2~(势阱一)跃迁的Franck-Condon因子、爱因斯坦自发辐射系数和辐射寿命.     

Low-lying electronic states of CuN calculated by MRCI method

The high accuracy ab initio calculation method of multi-reference configuration interaction(MRCI) is used to compute the low-lying eight electronic states of CuN.The potential energy curves(PECs) of the X~3∑~-,1~3Π,2~3∑~-,1~3△,1~1△,1~1∑~-,1~1Π,and ~5∑~- in a range of R=0.1 nm-0.5 nm are obtained and they are goodly asymptotes to the Cu(~2S_g) + N(~4S_u) and Cu(~2S_g)+N(~2D_u) dissociation limits.All the possible vibrational levels,rotational constants,and spectral constants for the six bound states of X~3∑~-,1~3Π,2~3∑~-,1~1△,1~1∑~-,and 1~1 Π are obtained by solving the radial Schrdinger equation of nuclear motion with the Le Roy provided Level 8.0 program.Also the transition dipole moments from the ground state X~3∑~- to the excited states 1~3Π and 2~3∑~- are calculated and the result indicates that the 2~3∑-X~3∑ transition has a much higher transition dipole moment than the 1~3Π-X~3∑~- transition even though the l~3Π state is much lower in energy than the 2~3∑~- state.     

SiCl自由基X2Π和A2Σ+态的光谱性质

采用内收缩多参考组态相互作用(ic MRCI)方法结合Dunning等相关一致基,计算Si Cl自由基X²Π和A²Σ⁺态的势能曲线.讨论参考能和相关能外推对X²Π和A²Σ⁺态光谱的影响.对势能进行相对论修正及核价修正计算.拟合势能曲线得到X²Π和A²Σ⁺态的光谱常数.它们与实验结果一致.利用Breit-Pauli算符,计算旋轨耦合效应,得到X²Π_1/2和X²Π_3/2的势能曲线、并计算它们的光谱常数.求解双原子分子核运动的径向SchrÖdinger方程,获得无转动SiCl自由基2个Λ-S态及X²Π态的耦合分裂态的全部振动态.得到J=0时X²Π态的自旋-轨道耦合常数、较高振动态的惯性转动常数以及X²Π_1/2和X²Π_3/2的振动能级等分子常数.     

GeO分子激发态的电子结构和跃迁性质的组态相互作用方法研究

应用多参考组态相互作用方法计算了GeO分子的第一解离极限(Ge(³P_g)+O(³P_g))对应的18个Λ-S电子态的电子结构. 计算中纳入了Ge原子的3d轨道电子的内壳层-价壳层电子关联效应、标量相对论效应和Davidson修正. 基于计算的电子态的电子结构, 通过求解径向Schrödinger方程获得了束缚电子态的光谱常数R_e, T_e, ω_e, ω_eχ_e, B_e, 理论计算给出的这些电子态的光谱常数与之前的实验结果符合得很好. 计算了电子态的电偶极矩随核间距的变化, 分析了电子态的组态成分的变化对电偶极矩的影响. 计算的势能曲线表明, 激发态A¹Π, 1¹Σ^-, D¹Δ, a³Π, a’³Σ⁺, d³Δ 和 e³Σ^-的绝热激发能密集地分布于26000-37000 cm^-1范围内, 这些密集分布的电子态之间的相互作用对振动波函数有明显扰动作用. 借助于激发态之间的自旋-轨道耦合矩阵元, 阐明了邻近的激发态对A¹Π和a³Π的扰动作用. 基于计算的A¹Π-X¹Σ⁺和A’¹Σ⁺-X¹Σ⁺跃迁的电偶极跃迁矩和Franck-Condon 因子, 给出了A¹Π 和A’¹Σ⁺态的最低的六个振动能级的辐射寿命.     

Theoretical study of potential energy curves,spectroscopic constants,and radiative lifetimes of low-lying states in SeO molecule

The low-lying potential energy curves of SeO molecule are computed by means of ab initio multireference configuration interaction technique, taking into account relativistic (scalar plus spin-orbit coupling) effects. The spectroscopic constants of Ω states for X³Σ^-, a¹Δ, b¹Σ⁺, A³Π, A'³Δ, and A"³Σ⁺ states are obtained, and they are in good accordance with available experimental values. The Franck-Condon factors and transition dipole moments to the ground state are computed, and the natural radiative lifetimes of low-lying Ω states are theoretically obtained. Comparisons of the natural lifetimes of Ω states with previous experimental results and those of isovalent TeO molecule are made.     

激光冷却OH分子的理论研究

采用高精度的多参考组态相互作用方法计算OH分子基态和第一激发态的势能曲线.为获得更精确的计算结果,在计算过程中考虑Davidson修正、标量相对论效应、核价相关效应和自旋轨道耦合效应.基于计算的Λ-S和?态的势能曲线,对一维径向薛定谔方程进行数值求解,得到各个电子态的光谱数据,与已报道的实验值和理论值相符合.获得OH分子的永久偶极矩、跃迁偶极矩、振动能级、Franck-Condon因子及辐射寿命,结果表明,A~2Σ~+→X~2Π跃迁具有高度对角化的Franck-Condon因子(0.9053)和短的辐射寿命(5.8363×10~(-7)s),符合激光冷却分子的条件.制定了激光冷却OH分子的具体方案,计算得到激光冷却跃迁A~2Σ~+→X~2Π所需的三束激光波长,主光束波长为307.1532 nm,两束重抽运激光波长为344.9163和349.7659 nm.计算结果为超冷OH分子的实验制备提供重要的理论依据.     

用差分收敛法研究NaLi分子部分电子态的完全振动能谱

基于双核分子振动能级的普遍表达式和差分收敛法(difference converging method, DCM), 利用微分思想将DCM应用于双核分子体系完全振动能谱的研究中. 应用DCM方法, 分别选用实验上获得的一组(10条)精确的振动能级, 对NaLi分子3¹Π, 4¹Π 和A¹Σ⁺电子态进行了研究, DCM的研究结果正确重复了已知数据并预测出了在实验上未能获得的包含高激发态在内的完全振动能谱数据, 同时计算得到了这3个电子态的振动光谱常数.     

Na_2~+离子较低电子态势能曲线和光谱常数的理论研究

通过多组态相互作用方法,结合原子有效芯势与极化势,利用非收缩的高斯基函数,计算了Na_2~+分子对应最低9个解离限的36个电子态的势能曲线.基于计算获得的束缚态势能曲线,拟合给出了相应的光谱常数,并与已有的实验和理论结果进行了比较.同时,给出了部分电子态的振动-转动能级和一些同类态避免交叉点的信息.计算获得的光谱信息对冷原子分子光谱与动力学的研究具有参考价值.     

基于ab initio的BD~+离子激光冷却理论研究

基于相关一致基组aug-cc-pV5Z,采用内收缩多参考组态相互作用方法计算了BD~+离子两个最低解离极限B~+(~1Sg)+D(~2Sg)和B+(~3Pu)+D(~2Sg)对应的5个Λ—S态(X~2Σ~+,A~2Π,B~2Σ+,a~4Π和b~4Σ~+)的势能曲线和跃迁偶极矩.根据计算结果,求解核运动的径向薛定谔方程得到相应电子态的振-转光谱常数、Franck-Condon(F-C)因子和振动能级辐射寿命.其中A~2Π—X~2Σ~+的F-C因子(f_(00)=0.923)、辐射寿命(τ=235 ns)满足激光直接冷却的条件.因此,我们基于分子转动跃迁提出了一个可实现Doppler激光冷却的光循环方案:A~2Π_(1/2)(υ′=0)—X~2Σ+(υ′′=0,1),其中υ′=0中包含2个转动能级,υ′′=0和υ′′=1中分别包含6个和4个转动能级.根据方案,模拟了激光冷却过程中的分子布居数动力学变化过程,并计算了初速度为100 m/s的BD~+,历经5.4 ms散射1150个光子可减速到4.6 m/s、温度为13 mK.     

75As32S+和75As34S+离子的光谱常数与分子常数

采用内收缩多参考组态相互作用(MRCI)方法, 结合Dunning系列相关一致基, 分别对⁷⁵As³²S⁺和⁷⁵As³⁴S⁺离子的X³Σ^-和A¹Π电子态的势能曲线进行了计算, 进一步拟合势能曲线, 得到各电子态的光谱常数与分子常数. 首先, 采用MRCI方法结合相关一致基, aug-cc-pV5Z, 对AsS⁺离子的X³Σ^-和A¹Π 电子态进行了计算, 获得相应的势能曲线; 然后, 为进一步提高势能曲线的精度, 对其进行了三种修正计算. 采用Davidson(+Q)方法修正MRCI 方法计算过程中存在的基组大小不一致缺陷; 利用二阶Douglas-Kroll哈密顿近似, 在cc-pVQZ基组水平, 修正了相对论效应对势能曲线的影响; 利用两点能量外推法, 在aug-cc-pVQZ和aug-cc-pV5Z基组水平对各能量点的势能值进行了外推, 得到完全基组极限处的势能曲线. 最后, 利用修正(包括Davidson修正、相对论修正和基组外推)后的势能曲线, 通过Vibrot程序, 求解双原子分子核运动的径向Schrödinger方程, 并进行同位素质量识别, 得到⁷⁵As³²S⁺和⁷⁵As³⁴S⁺离子两个电子态的光谱常数(T_e, R_e, ω_e, ω_ex_e, α_e 和B_e)和分子常数(G(ϒ), B_v, D_v).     

从头计算研究BP分子的势能曲线和光谱性质

利用高精度的量子化学从头计算MRCI+Q方法结合相关一致aug-cc-pVQZ基组计算了磷化硼分子X3Π,3Σ-,5Π和5Σ-态的势能曲线,计算所得的电子态在大键长位置处收敛于同一个离解极限B(2Pu)+P(4Su).为了得到更精确的结果,计算中首次纳入了旋轨耦合(SOC)效应,使得BP分子的4个Λ-S态分裂成为15个Ω态,其中3Π0+态被确定为基态.此外,SOC效应还使两个三重态X3Π和3Σ-分裂出的0+和1态的势能曲线产生了避免交叉,表明在当前的计算中考虑SOC效应是非常必要的.利用LEVEL8.0程序对计算所得的Λ-S态和Ω态的势能曲线进行拟合,得到了相应的光谱常数,通过与其他理论和实验工作进行比较,可知我们的结果更加精确、完整,可以为实验和理论方面进一步研究BP分子的光谱性质提供可靠的参考.     

反应物NO的振动激发对反应C(3P)+NO(X2Π)→CO(X1Σ+)+N(2D)立体动力学性质的影响

采用准经典轨线方法研究碰撞能为0.23 eV时,反应物分子NO在不同初始振动态(v=0～3)下发生在两个电子态(2A″和2A′)势能面上反应C(3P)+NO(X2Π)→CO(X1Σ+)+N(2D)的立体动力学性质.计算反应产物的转动角动量矢量分布(P(θr)和P(?r))以及微分散射截面(P00(ωt),(P)20(ω...     

Structural,Spectroscopic and Vibrational Properties for Low-Lying Electronic States of LiCl

A multireference configuration interaction (MRCI) study has been carried out on the LiCl molecule. The potential energy has been calculated over a wide range of internuclear separation for the 21 low-lying electronic states of the LiCl molecule dissociating into Li (²S, ²P, ³S)+Cl(²P). The (4)¹Σ⁺, (3)¹Π, 1-3³Σ⁺, 1-3³Π, ^1,3Δ, ^1,3Σ^-, (5)¹Σ⁺, (4)³Σ⁺, (4)¹Π, (4)³Π excited states are studied for the first time in theory. Molecular spectroscopic constants (R_e, D_e, ω_e, ω_eχ_e, B_e and α_e) have been derived for the 9 bound states (X¹Σ⁺, (3)¹Σ⁺, (2)¹Σ⁺, ^1,3Δ, ^1,3Σ^-, (4)¹Π, (4)³Π) with a regular shape, and the spectroscopic constants of ground states X¹Σ⁺ are in good agreement with available experimental and theoretical values. The relative differences between experimental values and our values for R_e, D_e, ω_e, ω_eχ_e, B_e and α_e are 1.02%, 0.60%, 1.72%, 9.46%, 2.0%, and 0.75%, respectively. Moreover, vibrational levels of 9 bound states, which have not been investigated experimentally, are computed.     

SnO分子的X~1∑~+,a~3Π和A~1Π态的势能曲线与光谱性质

对O原子以aug-cc-pvTZ为基组,Sn原子以aug-cc-pvTZ-PP为基组,分别采用多参考组态相互作用方法(MRCI)及运用Davidson修正的多参考组态相互作用方法对SnO分子基态X~1Σ~+及两个激发态a~3Π和A~1Π态的势能曲线进行计算,进而得到了各态的平衡键长R_e,谐振频率w_2,非谐振常数ω_ex_e,转动常数B_r,垂直跃迁能T_e以及离解能D_e,通过群论原理确定了各电子状态和离解通道,计算结果表明:3个电子态有共同的离解通道,即Sn(~3p)+O(~3p);利用Level程序对势能曲线进行拟合得到的光谱数据表明,MRCI方法的计算结果与实验值符合更好;通过求解核运动的Schrodinger方程得到了J=0时这三个电子态的前30个振动态的B_v和D_v等分子常数和振动能级E。     

Configuration interaction study on low-lying states of AlCl molecule

High-level ab initio calculations of the Λ–S states for aluminum monoiodide(Al Cl) molecule are performed by utilizing the explicitly correlated multireference configuration interaction(MRCI-F12) method. The Davidson correction and scalar relativistic correction are investigated in the calculations. Based on the calculation by the MRCI-F12 method, the spin–orbit coupling(SOC) effect is investigated with the state-interacting technique. The adiabatic potential energy curves(PECs) of the 13 Λ–S states and 24 Ω states are calculated. The spectroscopic constants of bound states are determined,which are in accordance with the results of the available experimental and theoretical studies. Finally, the transition properties of 0~+(2)–X0~+, 1(1)–X0~+, and 1(2)–X0~+ transitions are predicted, including the transition dipole moments(TDMs),Franck–Condon factors(FCFs), and the spontaneous radiative lifetimes.     

CO分子四个电子态的振转谱:两种效应修正方法的比较

基于完全活性空间自洽场方法和多参考组态相互作用（multi-reference configuration interaction method,MRCI）方法,采用MRCI+Q/CBS（TQ5）+CV+SR（方法A）和aug-cc-pwCVnZ-DK（n=T,Q,5）（方法B）方案,分别计算了包含Davidson修正（+Q）、芯-价电子关联（core-valence correlation correction,CV）效应以及标量相对论（scalar relativistic,SR）效应的CO分子的基态X¹∑⁺和激发态a³Π,a'³∑⁺和A¹Π的势能曲线.在此基础上,获得了这些电子态的振-转谱.通过与实验结果比较发现：方法A适合a'³∑⁺和A¹Π等较高激发态的振-转谱的计算,方法B更适合基态X¹∑⁺和第一激发态a³Π的振-转谱的精细计算.该研究可以为其他小分子高精度振-转谱快速计算方案选择提供参考.     

Theoretical studies on 14 Λ-S and 30 Ω states of BeBr molecule: potential energy curves,spectroscopic parameters and spin–orbit couplings

Using the complete active space self-consistent field (CASSCF) method followed by the internally contracted multi-reference configuration interaction (MRCI) approach in combination with the correlation-consistent basis sets, this paper studies the potential energy curves of X²Σ⁺, 2²Σ⁺, 3²Σ⁺, 1²Σ^?, A²Π, 2²Π, 3²Π, 1²Δ, 1⁴Σ⁺, 2⁴Σ⁺, 1⁴Σ^?, 1⁴Π, 2⁴Π and 1⁴Δ Λ-S states of BeBr molecule and the corresponding 30 Ω states for the first time. All the Λ-S states correlate to the first two dissociation channels, Be(¹S_g) + Br(²P_u) and Be(³P_u) + Br(²P_u), of BeBr molecule. Of these Λ-S states, the 3²Π and 2⁴Π are found to be repulsive without the spin–orbit coupling, whereas 1⁴Π, 2⁴Π, 3²Π and 2⁴Σ⁺ are found to be repulsive with the spin–orbit coupling included. A²Π and 2²Σ⁺ possess the double well whether the spin–orbit coupling effect is included or not. Only 1⁴Σ⁺, 1⁴Σ^?, 1²Π and 2²Π are found to be the inverted Λ-S states. The spin–orbit coupling is accounted for by the state interaction approach with Breit–Pauli Hamiltonian using the all-electron cc-pCVTZ basis set. The potential energy curves determined by the internally contracted MRCI method are corrected for size-extensivity errors by means of the Davidson correction. Core–valence correlation correction is calculated with a cc-pCVTZ basis set. Scalar relativistic correction is included using the third-order Douglas–Kroll Hamiltonian approximation at the level of cc-pVTZ basis set. The spectroscopic parameters of all the Λ-S and Ω bound states are evaluated. The spectroscopic parameters are compared with those reported in the literature. Fair agreement is found between the present results and available measurements. In particular, the energy splitting of 204.43 cm^?1 in the A²Π Λ-S state agrees well with the measurements of 201 cm^?1. Analyses demonstrate that the spectroscopic parameters reported here can be expected to be reliably predicted ones.     

BH~+离子基态及激发态的势能曲线和跃迁性质的研究

利用高精度的多组态相互作用及Davidson修正方法 (MRCI+Q),采用ACV5Z-DK全电子基组计算了BH~+离子的前4个离解通道B~+(~1S_g)+H(~2S_g),B~+(~3Pu)+H(~2S_g),B(~2P_u)+H~+(~1Sg)和B~+(~1Pu)+H(~2Sg)的9个Λ—S态的势能曲线.X~2Σ~+,A~2Π和B~2Σ~+态的光谱常数和已有实验值符合得很好,其中b~4Σ~+,3~2Σ~+,3~2Π和4~2Σ~+态的光谱常数为首次报道,3~2Π和4~2Σ~+态具有双势阱结构.预测了A~2Π和B~2Σ~+态的辐射寿命:τ(A~2Π)=239.2 ns和τ(B~2Σ~+)=431.2 ns.最后在考虑自旋轨道耦合效应下讨论了B~2Σ~+与A~2Π态的势能曲线的相交对激光冷却BH~+离子的影响.