SbS电子基态及激发态势能曲线和振动能级的理论研究*

王新宇 王艺霖 石虔韩 汪庆龙 于洪洋 金园园 李松

(长江大学物理与光电工程学院,荆州 434023)

运用多参考组态相互作用(MRCI+Q)方法,对硫化锑(SbS)能量最低的3个Λ-S离解极限的所有电子态以及考虑自旋-轨道耦合效应后分裂所得的Ω态进行了计算.得到27个Λ-S电子态及能量最低的12个Ω态的电子结构、光谱常数和振动能级等信息.Sb原子和S原子能级的计算值与实验值相符很好.分析表明自旋-轨道耦合效应对光谱常数与振动能级的影响总体上并不显著.对X(3/2)→X(1/2),2(1/2)→X(1/2),4(1/2)→X(1/2),5(1/2)→X(1/2)及6(1/2)→X(1/2)跃迁的振动光谱进行了模拟与分析,其中X(3/2)→X(1/2)谱带位于中红外波段,其他谱带均位于可见光波段.此外,对氮族元素硫化物的电子态进行了验证计算,计算结果与已有实验结果吻合较好,体现了同族元素代换后相关物性的渐变规律性.

1 引言

含硫双原子体系在天体物理学、大气化学、燃烧化学、分子反应动力学等众多领域承担着重要角色,因此一直是相关领域的关注对象.硫化锑(SbS)具备优异的稳定性和丰富的元素储存,因有较大的吸收系数和1.7 eV的带隙宽度,作为良好的半导体材料和光敏材料得到了广泛应用[1−5],而且锑基硫族化合物也满足叠层太阳能电池的要求,有助于提高光电转换效率[6].Shimauchi和Nishiyama[7]于1968年对SbS自由基的电子结构与发射电子光谱进行了报道,他们确定了7个电子激发态至基态的电子跃迁谱带的带头波长,但没有对激发态进行标识.Ω基态的谐振频率分别为480 cm–1和470 cm–1,而7个激发态的谐振频率介于296—442 cm–1之间.除此以外,其他光谱常数目前仍然未知.

对氮族元素硫化物自由基的研究始于1932年对NS的 B2Σ+→X2Π和 A2Π→X2Π 谱带的实验探测[8].1951年和1954年,Zeeman等[9,10]分别对这两个谱带进行了实验光谱转动分析.此后开展的微波谱[11]、紫外与可见光[12−14]实验获得了丰富的电子激发态光谱常数与基态的精细结构常数.对PS自由基的光谱研究源于对紫外与可见光波段C2Σ−X2Π和 B2Π−X2Π 谱带的实验探测.1955—1979年[15−20],若干个工作组探测到了这两个谱带的大量谱线,通过转动分析确定了各电子态的光谱常数.随后,该体系的近红外[21]、毫米波[22]与亚毫米波[23]光谱也陆续被探测获得.Shimauchi 研究组[24−28]在1969—1973年对AsS的光谱开展了系列研究,获得了 A2Π3/2−X2Π3/2跃迁的大量数据及各电子态的分子常数.1967年Barrow等[29]首次对BiS的可见光谱进行了探测.Patiño等[30]于1984年研究了BiS的 A2Π1/2−X2Π1/2谱带高J量子数跃迁的超精细双分裂结构,确定了 X2Π1/2态的分子常数.借助傅里叶变换光谱仪,BiS的近红外[31,32]、微波谱[31]和可见光光谱[32]也被探测到,涉及 X2Π1/2和 X2Π3/2态的分子常数得以确定.

在理论研究方面,Hartree-Fock方法(HF)[33]、组态相互作用(CI)[34]、多参考双重激发组态相互作用(MRDCI)[35−37]、密度泛函理论(DFT)[38]、广义价键(GVB)[39]、二阶Møller-Plesset微扰理论(MP2)[39]、完全活性空间-平均耦合对泛函(CASACPF)[40]、单双迭代包括三重激发的耦合簇[CCSD(T)][41,42]及其显关联方法CCSD(T)-F12[43]、Gaussian-3(G3)[44]、多参考组态相互作用(MRCI)[45−52]及Rydberg-Klein-Rees (RKR)[53]等计算方法均被用于确定NS,PS,AsS和BiS等4种含硫双原子体系各电子态的结构参数、光谱常数、振动能级以及跃迁性质等,将所得计算数据分别与各自的实验数据进行了比较.并且预测了实验中没有涉及的Λ-S态及Ω激发电子态的特性,得到了这些态的跃迁偶极矩、弗兰克-康登因子、爱因斯坦系数及自发辐射寿命等数据.

与对氮族元素硫化物自由基的结构和光谱的众多研究所获得的丰富的数据及结论相比,对SbS自由基的相关研究明显不足.除了前文提到的一篇实验研究报道[7]外,到目前为止还没有相关的理论研究.因此,本文对SbS的结构和电子态进行了系统研究,以填补相关数据的空白.

2 计算方法

基于完全活性空间自洽场方法[54]构建了CI波函数,借助包含Davidson修正的多参考组态相互作用(MRCI+Q)方法[55],计算了SbS前三个离解极限27个Λ-S电子态的能量.对S原子和Sb原子分别选用aug-cc-pwCV5Z全电子基组[56]与aug-cc-pwCV5 Z-PP标量相对论基组[57],其中Sb原子的1s—3d电子用相对论有效原子实势ECP28-MDF取代.计算中将Sb原子的4s4p与S原子的1s电子作为芯电子,Sb原子的4d和S原子的2s2p原子轨道作为闭壳层分子轨道,Sb原子的5s5p与S原子的3s3p原子轨道作为活性分子轨道.以C2v群替代SbS的简并对称性C∞v群,用其不可约表示a1,b1,b2和a2表示的双占据闭壳层轨道和活性轨道分别为(4,2,2,1)与(4,2,2,0).能量点的计算范围是1.7 Å—10.0 Å (1 Å=0.1 nm),最小扫描步长为0.05 Å,每条势能曲线计算了51个数据点.

在考虑核价相关修正和标量相对论修正后,通过Murrell-Sorbie(M-S)势能函数[58]拟合单点能得到电子态的势能曲线,然后通过均方根值(RMS)来评估拟合效果.M-S函数是能够较好地反映双原子体系势能函数的解析表达式之一,本研究组也已经基于该势能函数研究了NS–[59],SCl+[60,61],SCl–[60],SF±[62],MgS+[63],SH–[64]等若干含硫双原子体系.M-S势能函数定义为:

其中,ρR−Re,R与Re分别是核间距以及平衡核间距,De是离解能,ai是拟合参量.

通过Breit-Pauli算符[65]考虑自旋-轨道耦合(SOC)效应可以计算得到Ω电子态的能量.在核价相关修正和标量相对论修正的基础上,将能量点通过最小二乘拟合法得到势能函数,并计算出各电子态的光谱常数.以上所有能量计算均基于MOLPRO软件[66]完成.利用Level程序[67]还计算得出每个电子态的振动能级与转动常数.

鉴于SbS还没有实验数据可做比较,本文还计算了PS,AsS,BiS这3种氮族元素硫化物的若干电子态.采用的方法与计算SbS的方法相同,均通过MRCI+Q方法进行.对S原子,N原子和P原子采用aug-cc-pwCV5Z基组[56],对As原子和Bi原子采用aug-cc-pwCV5Z-PP基组[57],其中分别包含ECP10MDF和ECP60MDF相对论有效原子实势.得到了每一种体系第一离解极限的Λ-S电子态及其对应的Ω电子态能量,进一步计算出每个电子态的光谱常数与振动能级.

3分析与讨论

3.1 Λ-S态的势能曲线与光谱常数

SbS前三个离解极限相对能量的计算值及对应的电子态如表1所列,本文计算值与实验值相比,符合得很好.例如S原子第一激发态1Dg相对基态3Pg的能量为9346 cm–1,与实验值[68]相差约1.2%,而Sb原子第一激发态2Du与基态4Su能量间隔为10022 cm–1,高于实验值[68]约1.7%.所有Λ-S电子态势能曲线如图1所示,其中 14Σ+,16Σ+,16Π,22Δ,34Σ+,34Δ,14Φ,34Π,44Π和 54Π 为排斥态,其余电子态均为束缚态.

图1 SbS的Λ-S态势能曲线 (a)二重、六重态;(b),(c)四重态Fig.1.Potential energy curves of Λ-S states of SbS:(a) Doublet and sextet states;(b),(c) quartet states.

表1 SbS的Λ-S态离解极限Table 1.Dissociation relationships of the Λ-S states of SbS.

基态 X2Π与第一激发态 14Π 在Re附近的能量差超过了12800 cm–1,并且没有其他电子态势能曲线与基态曲线交叉,表明与其对应的Ω态不会受到其他Ω=1/2或Ω=3/2电子态的影响,其光谱常数也不会有大的变化.在R=2.4—3.4 Å,E=25000—40000 cm–1范围内,激发态势能曲线产生了复杂的曲线(避免)交叉,预示了在此范围内对其Ω态的分析将具有很大的挑战性.

束缚态完整的光谱常数列于表2.需要说明的是,本文所得电子激发态的谐振频率总体上与文献[7]的数据(介于296—442 cm–1之间)是相符的,但是由于文献[7]没有标识电子态,因此无法与其数据进行比较.拟合的RMS值均较小,表明拟合质量较高.基态 X2Π 主要由电子组态构成,相比激发态其势阱最深,但Re最小.第一激发态 14Π 通过7πy→8πy电子迁移形成,虽然其Re与第二、第三激发态 22Π和 32Π 的Re非常接近,但均与基态相差超过9%,因此可预测这几个低激发态至基态跃迁的弗兰克-康登因子偏小.22Π和 32Π 均呈现多组态特征,贡献最大的电子组态分别由自旋取向不同的电子保持自旋方向性,并从7πy迁移至8πy轨道而形成.14Σ−的主要电子组态为,所占权重为83%.8π→16σ的电子迁移形成12Σ+电子态.除个别电子态以外,大多数电子态均表现出较明显的多参考特性.

表2 SbS的Λ-S态光谱常数Table 2.Spectroscopic constants of the Λ-S states of SbS.

表3汇总了氮族元素硫化物自由基的光谱常数.通过比较,发现该系列硫化物基态 X2Π 的光谱常数体现了同族元素代换后的渐变规律性.随着氮族元素核电荷数的增加,Re逐渐变大,这源于氮族元素np3价电子的弥散性渐强,而氮族元素与S原子之间的化学键减弱则导致谐振频率ωe逐渐减小.

表3 XS (X=N,P,As,Sb,Bi)自由基电子基态 X2Π 的光谱常数Table 3.Spectroscopic constants of the ground state X2Π of XS (X=N,P,As,Sb,Bi) radicals.

3.2 Ω态的势能曲线与光谱常数

由于SbS电子态数量多且曲线(避免)交叉复杂,本文仅对基态及部分低激发Ω态计算与讨论,其势能曲线如图2所示.计算涉及的Ω态离解极限,即S原子3P2,1,0原子态能级间隔与实验数据[68]吻合很好,第二、第三离解极限的计算值与实验值相差约为 3.6%和5.6% (见表4).各电子态的光谱常数列于表5.其中,仅有Ω基态的谐振频率有实验值[7]可做比较,本文计算值与实验值相符很好.各势能曲线拟合的RMS值也比较令人满意.

表4 SbS自由基Ω态的离解极限Table 4.Dissociation relationships of the Ω states of SbS.

图2 SbS的Ω态势能曲线Fig.2.Potential energy curves of Ω states of SbS.

在SOC作用下,Λ-S基态X2Π分裂为X(1/2)与X(3/2),其中前者能量更低,并且根据二者间能量差可预计其自旋-轨道耦合常数约为2025 cm–1.14Π分裂为2(1/2),3(1/2),2(3/2)和1(5/2),裂距较小,两相邻Ω态间的裂距只有基态裂距的11%.22Π分裂为裂距约700 cm–1的3(3/2)和4(1/2)电子态.与 14Σ−和 32Π 对应的Ω态势能曲线在约2.25 Å处产生了避免交叉.对比表2的Λ-S态及表5中对应Ω态的数据,发现光谱常数的变化不大,证实了本文的预测.以Re为例,Ω基态X(1/2)与X(3/2)相比,Λ-S基态 X2Π 的变化分别为0.02%和0.01%,在表中所列电子态中是最小的,因此可预测这两个态之间的跃迁会有较大的弗兰克-康登因子.Re变化率最大的态是6(1/2)和5(3/2),不过二者也分别只比 14Σ− 减小了0.46%和0.52%.从总体上看,SOC效应对这些电子态光谱常数的影响较小.此外,由于Ω基态与激发态Re相差达到了0.2 Å,表明这些Ω态间跃迁的弗兰克-康登因子均较小.

表5 SbS自由基的Ω态光谱常数Table 5.Spectroscopic constants of the Ω states of SbS.

表6列出了氮族元素硫化物自由基Ω基态的光谱常数,由图6可见,本文计算值与实验值非常相符.以BiS为例,X(1/2)与X(3/2)态Re的计算值均只与实验值[31,32]相差小于0.2%,ωe则分别高估了5%和2%左右.计算值与实验值的一致性,从侧面证明本文对SbS的计算结果具有很高的精度.

表6 XS (X=N,P,As,Sb,Bi)自由基Ω基态的光谱常数Table 6.Spectroscopic constants of the ground Ω state of XS (X=N,P,As,Sb,Bi) radicals.

表6（续） XS (X=N,P,As,Sb,Bi)自由基Ω基态的光谱常数Table 6(continued).Spectroscopic constants of the ground Ω state of XS (X=N,P,As,Sb,Bi) radicals.

3.3 Λ-S态与Ω态的振动能级和转动常数

通过求解核运动的径向薛定谔方程,得到了SbS的 X2Π,14Π,22Π,32Π,14Σ−电子态及其对应Ω态的全部振动态.表7列出了v=0—5的振动能级、转动常数和离心畸变常数.受SOC效应影响,Ω态中6(1/2)和5(3/2)的振动能级相比Λ-S态均降低了15%左右,是这些低激发态中变化最大的.该现象源于这两个态分别与5(1/2)和4(3/2)在Re附近产生曲线避免交叉,使得6(1/2)和5(3/2)的势能曲线在Re处相比 14Σ− Λ-S态的绝对能量降低了约300 cm–1.总体上看,SOC效应对SbS四重态,如 14Π和 14Σ− 影响相对更为显著,而对二重态的影响则不明显.

表7 SbS的Λ-S及其对应Ω态的振动能级、转动常数和离心畸变常数(单位:cm–1)Table 7.Vibrational energy levels,rotational constants and centrifugal distortion constants for the Ω and its respective Λ-S states of SbS (in cm–1).

表7（续） SbS的Λ-S及其对应Ω态的振动能级、转动常数和离心畸变常数(单位:cm–1)Table 7(continued).Vibrational energy levels,rotational constants and centrifugal distortion constants for the Ω and its respective Λ-S states of SbS (in cm–1).

为证明本文计算结果的准确性,对验证计算并已获得光谱常数的PS,AsS和BiS进行了振动分析,计算数据及相应实验结果列于表8可见,本文计算得到的PS自由基X(1/2)与X(3/2)的转动常数及离心畸变常数与实验值[19]符合得非常好,其中偏差最大的是X(3/2)态v=7的转动常数(约0.8%),其他能级转动常数与实验值的偏差均在0.5%左右.对于AsS,X(3/2)态v=3—7能级的转动常数与实验值[25]相差不到0.1%.对于BiS来讲,只有X(1/2)的振动能级有实验值[29],计算值高估了实验值约0.6%.

表8 XS (X=P,As,Bi)自由基Ω基态的振动能级、转动常数和离心畸变常数(单位:cm–1)Table 8.Vibrational energy levels,rotational constants and centrifugal distortion constants for the ground Ω state of XS(X=P,As,Bi) radicals (in cm–1).

3.4 Ω态的振动跃迁

借助Level程序计算了若干Ω激发态至基态振动跃迁的爱因斯坦系数Av′J′v′′J′′及跃迁频率vv′J′v′′J′′.通过(2)式计算了跃迁强度Iv′J′v′′J′′,式中J′是上态转动量子数,Ev′′J′′是下态能量,Q(T)是对应温度T的配分函数,h,c,k是基本常量.

如图3(a)—(c)所示,在所有计算的谱带中,2(1/2)→X(1/2)的v′→v′′0,1,2 跃迁强度相对更大,主要分布在540—750 nm之间.4(1/2)→X(1/2)与5(1/2)→X(1/2)谱带则主要位于400—550 nm波段,其中前者除v′→v′′0 外跃迁强度均较小,而后者v′0,1,2,3,4→v′′0 跃迁则分布较密集.6(1/2)→X(1/2)跃迁的最强谱带位于430 nm附近,v′→v′′0,1,2跃迁的强度随v′′的增加而渐大.以上4个激发态至Ω基态的跃迁均位于可见光区域.

图3 SbS的振动谱带Fig.3.Vibrational transition bands of SbS.

图3(d)绘制了X(3/2)→X(1/2)跃迁的Δv−1,0,1谱带(仅绘制出v′≤15).3个谱带均位于中红外波段,其分布与位于近红外波段的BiS相应谱带分布(文献[31],图2)非常相似.强度最大的Δv0谱带位于 Δv−1与 Δv1谱带间,Δv1谱带的强度最小,在3个谱带中位于长波长一端.以上数据与结论均能够对SbS的光谱探测提供理论支持和数据支撑.

4 结论

通过MRCI+Q方法计算了SbS电子基态及低激发态的电子结构,得到了能量最低的三个Λ-S离解极限所有的电子态及部分Ω态的势能曲线.通过离解极限处的能量计算所得Sb原子与S原子的能级与实验值相符很好.计算了各电子态的光谱常数、振动能级,模拟了Ω激发态至基态的振动光谱,为后续开展光谱探测提供了参考依据.本文还对PS,AsS,BiS的电子态进行了验证计算,所得的光谱常数和振动能级均与已有的实验结果相符,间接证明本文SbS计算结果的可信度.