-
Potential energy curves (PECs), permanent dipole moments (PDMs) and transition dipole moments (TMDs) of five Λ-S states of SeH− anion are calculated by the MRCI + Q method with ACVQZ-DK basis set. The core-valence corrections, Davidson corrections, scalar relativistic corrections, and spin-orbit coupling (SOC) effects are also considered. In the CASSCF step, Se(1s2s2p3s3p) shells are put into the frozen orbitals, which are not optimized. Six molecular orbitals are chosen as active space, including H(1s) and Se(4s4p5s) shells, and eight electrons are distributed in a (4, 1, 1, 0) active space, which is referred to as CAS (8, 6), and the Se(3d) shell is selected as a closed-shell, which keeps doubly occupation. In the MRCI step, the remaining Se(3d) shell is used for core-valence calculations of SeH− anion. The SOC effects are taken into account in the one- and two- electron Breit-Pauli operators. The b3Σ+ state is a repulsive state. Other excited states are bound, and all states possess two potential wells. The b3Σ+0− andb3Σ+1 both turn into bound states when the SOC effect is considered. All spectroscopic parameters of Λ-S states and Ω states are reported for the first time. The TDMs of theA1Π1↔X1Σ+0+ ,a3Π1↔X1Σ+0+ ,a3Π0+↔X1Σ+0+ ,A1Π1↔a3Π1 , andA1Π1↔a3Π0+ transitions are also calculated. The TDMs of theA1Π1↔X1Σ+0+ anda3Π1↔X1Σ+0+ transitions are large in the Franck-Condon region, which are about –2.05 Debye (D) and 1.45 D at Re. Notably, the TDMs of thea3Π0+↔X1Σ+0+ transition cannot be ignored. The value of TDM at Re equals –0.15 D.Based on the accurately PECs and PDMs, the values of Franck-Condon factor fυ′υ″, vibrational branching ratio Rυ′υ″ and radiative coefficient of the a3Π1↔X1Σ+0+ ,a3Π0+↔X1Σ+0+ , andA1Π1↔X1Σ+0+ transitions are also calculated. Highly diagonally distributed Franck-Condon factor f00 and the values of vibrational branching ratio R00 of thea3Π1(υ′)↔X1Σ+0+(υ″) ,a3Π0+(υ′)↔X1Σ+0+(υ″) , andA1Π1(υ′)↔X1Σ+0+(υ″) transitions are obtained, respectively. Spontaneous radiation lifetimes of thea3Π1 ,a3Π0+ , andA1Π1 excited states are all short for rapid laser cooling. The influences of intervening states of theA1Π1(υ′)↔X1Σ+0+(υ″) transition can be ignored. The proposed cooling wavelengths using thea3Π1(υ′)↔X1Σ+0+(υ″) ,a3Π0+(υ′)↔X1Σ+0+(υ″) , andA1Π1(υ′)↔X1Σ+0+(υ″) transitions are all in the visible region.1. 引 言
由于分子内部能级的复杂性, 首次激光冷却分子的实验于2010年实现, Shuman等[1]采用横向激光冷却技术对SrF分子进行了冷却. 2013年, Hummon等[2]采用同样的技术实现了YO分子的冷却. 这是1个三能级体系, 其存在1个中间态A′2Δ3/2. 2014年, Zhelyazkova等[3]采用纵向激光冷却技术对CaF分子进行了冷却. 近十年, 有很多学者对一系列的双原子分子进行了激光冷却的理论研究, 如MgH[4], LiBe[5], CH[6], AlCl[7], OH[8], BaF[9]等.
相比于中性分子和分子阳离子, 双原子分子阴离子具有更加复杂的能级, 光谱数据非常缺乏, 直到2015年, Yzombard等[10]才首次从理论上预测了激光冷却双原子分子阴离子的可能性. 随后激光冷却阴离子吸引了越来越多学者的关注[11-16]. 鲁东大学的杨传路教授等预测了激光冷却NH–[11]和BH–[12]阴离子的冷却途径. Zeid等[13]在自旋无关水平下计算了XH– (X = Mg, Ca, Sr, Ba)体系的势能曲线和跃迁偶极矩. XH–体系的11Π
↔ X1Σ+跃迁均具有很大的跃迁偶极矩和对角化的弗兰克-康登因子, 他们认为这些阴离子可能是激光冷却的潜在候选体系. 在我们前期的工作中, 研究了激光冷却VI簇氢化物的阴离子(OH–, SH–, TeH–)的可能性[14-16]. SeH在实验上和理论上都是一种重要的自由基, 也是一种比较简单的重原子分子. 本文将主要研究在考虑自旋-轨道耦合效应(SOC)下SeH–阴离子的光谱和跃迁性质, 并讨论其激光冷却的可行性.迄今为止, 对SeH–阴离子的光谱常数的实验和理论研究很少. 1987年, Stoneman和Larson[17]在磁场中采用光致分离技术观察SeH–阴离子, 他们得到了SeH–基态的平衡核间距Re是SeH分子基态1.0038(3)倍. 同时采用最小二乘法得到了转动常数B为7.7289 cm–1. 1982年, Brown和Fackerell[18]采用CO激光磁共振观察了SeD和SeH分子的振转光谱, 得到了SeH分子基态X2Π的平衡核间距为1.4640715(25) Å. 1987年, Balasubramanian等[19]采用二阶组态相互作用方法计算了SeH和SeH+较低电子态的势能曲线, 并得到了SeH和SeH+基态的平衡核间距分别为1.467 Å和1.58 Å. 1990年, Binning Jr和Curtiss[20]采用Hartree-Fock方法计算得到SeH基态的平衡核间距为1.463 Å. 到目前为止, 没有发现针对SeH–阴离子的跃迁性质的理论或实验研究.
本文将在SOC水平下研究SeH–阴离子的基态(X1Σ+)和低激发(a3Π, A1Π, b3Σ+, 21Σ+)的势能曲线和跃迁偶极矩, 拟合得到各束缚态的光谱常数, 预测其弗兰克-康登因子、振动分支比和自发辐射系数, 最后讨论其激光冷却的可能性, 并给出其冷却途径.
2. 计算细节
本文采用从头算方法计算了SeH–阴离子4个Λ-S态(X1Σ+, a3Π, 11Π, b3Σ+)的电子结构. 在计算中增加了21Σ+态, 以便比较11Π和21Σ+态的能级高低. 所有计算都在MOLPRO 2010程序包[21]下完成. X1Σ+, a3Π, 11Π, b3Σ+态对应于最低离解极限Se–(2Pu)+H(2Sg), 第二离解极限为Se(3Pg) + H–(1Sg), 其对应的电子态只能是3重态, 所以第二个1Σ+态对应于第三离解通道Se(1Du) + H–(1Sg). 首先, 采用限制性的Hartree-Fock方法对SeH–阴离子进行波函数的初始猜测, 得到初始波函数; 然后, 采用完全活动空间自洽场方法(CASSCF)[22,23]产生多参考的波函数; 最后, 采用多参考组态相互作用方法(MRCI)[24,25]得到各Λ-S态的能量, 同时也考虑了Davidson修正[26]. 通过二阶Douglas-Kroll哈密顿量[27,28]来考虑标量相对论效应.
在CASSCF中, Se原子的1s2s壳层为冻结轨道, 6个分子轨道被选为活性轨道, 包括H原子1s和Se原子4s4p5s壳层, 8个电子占据了(4, 1, 1, 0)活性轨道, 写为CAS(8, 6), Se原子的1s2s2p3s3p为冻结轨道, Se原子的3d轨道为双占据的闭壳层. 在MRCI计算中, 价-芯(CV)电子关联被考虑, CASSCF过程中的闭壳层Se(3d)的10个电子参与CV关联计算. 而Se原子的1s2s2p3s3p仍被冻结. 此外, 本文在MRCI+Q水平下通过Breit-Pauli算符[29]考虑了SOC效应. 在Λ-S和Ω态的所有计算中, Se原子选择了aug-cc-pCVQZ-DK (ACVQZ-DK)全电子基组[30], H原子选择了aug-cc-pVQZ-DK(AVQZ-DK)全电子基组[31].
采用LEVEL8.0程序[32], 通过求解径向薛定谔方程得到SeH–阴离子低电子态的光谱常数、弗兰克-康登因子、自发辐射系数和自发辐射寿命.
3. 结果与讨论
3.1 势能曲线与光谱常数
计算了SeH–阴离子X1Σ+, a3Π, 11Π和b3Σ+态的势能曲线, 这4个电子态有1个共同的离解通道Se–(2Pu) + H(2Sg). 为了确定11Π态的能级顺序, 同时计算了第二个1Σ+态. 它对应于第三离解通道Se(1Du) + H–(1Sg). 所计算的电子态的势能曲线如图1所示. 可以看出11Π是单重态的第一激发态, 写为A1Π.
SeH–阴离子的基态为X1Σ+, 具有1个很深的势阱, 其离解能De为3.729 eV. X1Σ+态在平衡核间距Re处的主要组态为(1σ22σ23σ21π44σ25σ21π46σ22π41δ4)7σ28σ23π4, 简写为(core)7σ28σ23π4, 其权重为90.8%. Stoneman和Larson[17]观察到SeH–基态的平衡核间距Re(X1Σ+)是SeH分子Re(X2Π)的1.0038倍. Brown和Larson[18]测得SeH分子基态X2Π的平衡核间距Re为1.464 Å. 虽然Huber和Herzberg[33]收集了SeH分子基态的实验值为1.475 Å, 但结果不准确. 本文的计算值为1.4694 Å (表1), 只比实验值[17,18]小0.0002 Å, 相对误差仅为0.01%. 本文计算的转动常数Be为7.8507 cm–1, 比采用最小二乘法拟合得到的实验值[17]大1.58%. 可以看出本文计算结果与已有的实验值符合很好. 同时基态的谐振频率(ωe)和非谐振频率(ωeχe)也列于表1中, 分别为2300.77和46.10 cm–1. 同时在MRCI+Q/AVQZ-DK水平下计算了SeH–阴离子低激发态的势能曲线, 在计算中没有考虑CV关联效应, 基态的光谱常数也列于表1中, 可以看出考虑CV关联后, 基态的光谱常数(Re, Be)和已有实验值符合的更好.
表 1 Λ-S的光谱常数Table 1. Spectroscopic parameters of the Λ-S states.Λ-S态 来源 Re/Å ωe/cm–1 ωeχe/cm–1 Be/cm–1 De/eV Te/cm–1 X1Σ+ ACVQZ-DK 1.4694 2300.77 46.10 7.8507 3.487 0 AVQZ-DK 1.4614 2380.32 45.57 7.9326 3.711 实验[17] 1.4696 a 7.7289 c 1.4806 b a3Π 本文工作 第一势阱 1.4778 2206.52 123.45 7.8428 0.519 20642.90 第二势阱 2.1787 839.87 49.66 3.44016 0.450 24549.11 A1Π 本文工作 第一势阱 1.4726 2373.65 127.14 7.8391 0.734 21240.75 第二势阱 2.2780 437.62 44.07 3.0932 0.147 26997.57 b3Σ+ 本文工作 repulsive 21Σ+ 本文工作 第一势阱 1.6188 1336.45 — 6.1955 0.228 51684.73 第二势阱 4.0808 198.90 9.96 1.0190 0.135 46349.30 注: a 为SeH分子基态的平衡核间距的实验值来源于文献[18]; b为SeH分子基态的平衡核间距的实验值来源于文献[33], 结果不准确; c 为采用最小二乘法得到转动惯量B. 从图1可以看到b3Σ+是1个排斥态. 其他3个激发态都是具有双势阱结构的束缚态. a3Π, A1Π和21Σ+态的第一势阱分别位于比基态高20642.90, 21240.75, 51684.73 cm–1处, 势阱深度分别为0.518, 0.737, 0.205 eV; 其平衡核间距Re分别为1.4679, 1.4627, 1.6113 Å. 在平衡核间距位置处的主要组态分别为(core)7σ28σ29σα3
παx 3παβy , (core)7σ28σ29σα3πβx 3παβy 和(core)7σ28σα9σβ3π4, 权重分别为91.32%, 91.36%和88.28%. X1Σ+→A1Π和X1Σ+ → 21Σ+跃迁主要是由3πx → 9σ和8σ → 9σ跃迁引起的. 由于a3Π, A1Π和21Σ+态电偶极矩的曲线非常平滑, 没有发生突变, 可以推测这3个电子态都可能由预解离所致. 从图1可以看出这3个态分别约在1.98, 2.04, 2.12 Å处发生预解离. a3Π, A1Π和21Σ+态的第二势阱的平衡核间距Re分别为2.1787, 2.2780, 4.0808 Å. 其势阱深度分别为0.450, 0.147, 0.135 eV. 由于阴离子体系的光谱常数的实验值非常少, 本文对比了第VI主簇氢化物阴离子基态的光谱常数, 分别列于表2中, 可以看出随着第VI主簇原子质量的增加, 其氢化物阴离子基态的平衡核间距Re逐渐地增大, 但其氢化物阴离子基态的谐振频率ωe、非谐振频率ωeχe和势阱深度De都逐渐地减小.表 2 第VI主簇氢化物阴离子基态的光谱常数Table 2. Spectroscopic parameters of the ground state X1Σ+ of the Group VI-hydride anions.考虑SOC效应后, Se–离子的2P原子态分裂为2P1/2和2P3/2态.
X1Σ+0+ ,a3Π2 ,a3Π1 ,a3Π0− 和A1Π1 态对应于Se–(2P3/2) + H(2S1/2)离解通道,a3Π0+,b3Σ+0−,b3Σ+1 态对应于Se–(2P1/2) + H(2S1/2)离解通道. 本文预测了Se–(2P)原子态的分裂常数, ASO(2P1/2 – 2P3/2) = 2303.77 cm–1. 另外,21Σ+0+ 对应于较高的离解通道Se(1D2) + H–(1S0), 其与最低离解通道的相对能量为20032.24 cm–1, 比实验值[34-36]大241.36 cm–1, 相对误差仅为1.22%, 本文计算值与已有实验值符合很好. Ω态之间的离解关系列于表3. 表3也同样可以看出考虑CV关联效应后, 第三离解极限与最低离解极限的相对能量和已有实验值符合得更好.9个Ω态的势能曲线如图2所示, 所计算的电子态都是束缚态, 各Ω态的光谱常数列于表4. 对比表2, 基态
X1Σ+0+ 的势阱深度比X1Σ+态的要小约0.1 eV. 从图2可看出a3Π2 ,a3Π1 ,a3Π0− ,a3Π0+ ,A1Π1 和21Σ+0+ 态具有2个势阱.a3Π2 态是第一激发态, 其第一势阱到基态的垂直跃迁能为19787.17 cm–1. a3Π态的第一势阱的分裂常数为: ASO(a3Π1 –a3Π2 ) = 249.1 cm–1, ASO(a3Π0− –a3Π1 ) = 1436.25 cm–1, ASO(a3Π0+ –a3Π0− ) = 4.60 cm–1. 而a3Π态的第二势阱具有更大的分裂常数,a3Π0+ 和a3Π0− 态之间的分裂常数达到了192.26 cm–1. 从图2还可以看出a3Π0+ 和A1Π1 在2.89 Å处交叉. 考虑SOC效应后, b3Σ+态的2个分裂态b3Σ+0− 和b3Σ+1 态都具有了1个很浅的势阱, 势阱深度分别为0.096和0.066 eV, 平衡核间距分别为3.1807和3.2046 Å, 这两个态都有较大的谐振频率, 均超过200 cm–1, 说明这2个分裂态都是弱束缚态. 同时可以预测处b3Σ+态的分裂常数ASO(b3Σ+1 –b3Σ+0− )为239.22 cm–1. 同样21Σ+0+ 态的第二势阱也是弱束缚态. 由于A1Π1 的势阱太小, 可以认为其由弱束缚态变为排斥态. 由于a3Π和b3Σ+态的分裂常数很大, 可以看出SOC效应对SeH–阴离子的势能曲线和光谱常数的影响很大.表 4 Ω电子态的光谱常数Table 4. Spectroscopic parameters of the Ω states.Ω态 Re/Å ωe/cm–1 ωeχe/cm–1 Be/cm–1 De/eV Te/cm–1 X1Σ+0+ 1.4694 2301.31 47.01 7.8499 3.395 0 a3Π2 第一势阱 1.4777 2207.22 122.39 7.8416 0.523 19787.17 第二势阱 2.1739 861.02 52.10 3.4081 0.454 23751.54 a3Π1 第一势阱 1.4759 2232.16 111.70 7.8434 0.560 20036.27 第二势阱 2.1822 818.11 55.64 3.3929 0.386 24301.10 a3Π0− 第一势阱 1.4778 2205.83 124.97 7.8485 0.513 21472.52 第二势阱 2.1986 778.28 73.70 3.4048 0.267 25261.96 a3Π0+ 第一势阱 1.4777 2208.03 122.90 7.8422 0.522 21477.12 第二势阱 2.1619 904.15 49.02 3.4355 0.527 25454.22 A1Π1 第一势阱 1.4744 2368.50 144.22 7.8262 0.686 21821.04 第二势阱 排斥态 b3Σ+0− 第一势阱 3.1807 318.89 35.64 1.6988 0.096 28945.41 b3Σ+1 第二势阱 3.2046 239.13 30.94 1.6662 0.066 29184.63 21Σ+0+ 第一势阱 1.6190 1332.50 — 6.1895 0.225 51714.58 第二势阱 4.0800 190.57 8.84 1.0190 0.135 46351.94 3.2 Ω态的跃迁偶极矩
考虑SOC效应后,
a3Π2↔X1Σ+0+ 和a3Π0−↔ X1Σ+0+ 跃迁不被允许. 本文计算了A1Π1↔X1Σ+0+ ,a3Π1↔X1Σ+0+ ,a3Π0+↔X1Σ+0+ ,A1Π1↔a3Π1 和A1Π1↔a3Π0+ 跃迁的跃迁偶极矩, 其与核间距的函数关系如图3所示. 可以看出A1Π1↔X1Σ+0+ 和a3Π1↔X1Σ+0+ 跃迁在弗兰克-康登区域具有很大的跃迁偶极矩, 在平衡核间距处分别达到了–2.05 Debye (D)和1.45 D. 值得注意的是,a3Π0+↔ X1Σ+0+ 跃迁也具有较大的偶极矩, 在平衡核间距处为–0.15 D. 在讨论激光冷却SeH–阴离子的过程中必须要考虑a3Π0+↔X1Σ+0+ 跃迁. 由于X1Σ+0+ ,a3Π1 和A1Π1 来源于同一离解通道, 故当核间距大于6 Å后A1Π1↔X1Σ+0+ ,a3Π1↔X1Σ+0+ 和A1Π1↔a3Π1 跃迁的跃迁偶极矩均趋于零.a3Π0+↔X1Σ+0+ 和A1Π1↔a3Π0+ 跃迁在核间距趋于无穷大时源于Se–离子的2P3/2↔ 2P1/2跃迁. 由于其跃迁强度很小, 故核间距大于8 Å后a3Π0+↔X1Σ+0+ 和A1Π1↔ a3Π0+ 跃迁的跃迁偶极矩也趋于零.3.3 振动分支比及激光冷却SeH–阴离子的可能性
基于精确的势能曲线和跃迁偶极矩, 本文计算了
A1Π1↔X1Σ+0+ ,a3Π1↔X1Σ+0+ 和a3Π0+↔X1Σ+0+ 跃迁的弗兰克-康登因子fυ′υ″ 和自发辐射系数Aυ′υ″ . 由于在激光冷却循环过程中光子损失路径的相对强度与振动分支比Rυ′υ″ 有直接关系[37]. 振动分支比可以表示为Rυ′υ″=Aυ′υ″/∑υAυ′υ″, (1) 其中,
Aυ′υ″ 表示A1Π1↔X1Σ+0+ ,a3Π1↔X1Σ+0+ 和a3Π0+↔X1Σ+0+ 跃迁的自发辐射系数. 表5列出了以上3种跃迁的弗兰克-康登因子、自发辐射系数和振动分支比. 本文分别讨论自旋阻禁和直接跃迁2种情况来进行激光冷却SeH–阴离子的可能性.表 5 ,a3Π1↔X1Σ+0+ 和a3Π0+↔X1Σ+0+ A1Π1↔ 跃迁的辐射系数X1Σ+0+ 、弗兰克-康登因子Aυ′υ″ 和振动分支比fυ′υ″ Rυ′υ″ Table 5. Emission rates , Franck-Condon FactorsAυ′υ″ , branching ratiosfυ′υ″ of theRυ′υ″ ,a3Π1↔X1Σ+0+ , anda3Π0+↔X1Σ+0+ A1Π1↔ transitions.X1Σ+0+ Index a3Π1↔X1Σ+0+ a3Π0+↔X1Σ+0+ A1Π1↔X1Σ+0+ fυ′υ″ f00 0.9949 0.9922 0.9974 f01 0.0047 0.0072 0.0025 f02 0.0004 0.0006 0.0001 f10 0.0051 0.0079 0.0026 f11 0.9541 0.9324 0.9792 f12 0.0337 0.0486 0.0159 Aυ′υ″/s A00 5.02×106 8.02×104 1.36×107 A01 1.88×102 4.28×103 1.87×104 A02 2.81×101 7.48×101 2.00×103 A10 1.10×105 6.50×102 5.79×104 A11 4.13×106 9.13×104 1.32×107 A12 1.32×104 1.57×104 1.45×105 Rυ′υ″ R00 0.99996 0.9484 0.9985 R01 3.7×10–5 0.0506 0.0014 R02 5.6×10–6 0.0009 0.0001 R10 0.02592 0.0060 0.0043 R11 0.9707 0.8394 0.9836 R12 0.0031 0.1446 0.0108 3.3.1 自旋阻禁跃迁
a3Π1(υ′)↔X1Σ+0+(υ″) ,a3Π0+(υ′)↔X1Σ+0+(υ″) 阻禁跃迁都有高对角分布的弗兰克-康登因子f00, 分别为0.9949和0.9922, 相应的振动分支比R00分别为0.99996和0.9484. 高对角分布的f00和R00使SeH–阴离子满足了激光冷却的首要条件. 同时本文也计算了两种跃迁的非对角项振动分支比. 对于两种跃迁, R00, R01和R02的代数和非常接近1, 这可以保证在激光冷却SeH–阴离子过程中的准循环跃迁.a3Π1(υ′)↔X1Σ+0+(υ″) 和a3Π0+(υ′)↔ X1Σ+0+(υ″) 跃迁也有很大的自发辐射系数, 当υ′ = 0时, 这两种跃迁总的自发辐射系数分别为5.18 × 106和8.45 × 104 s–1, 因此a3Π1 和a3Π0+ 态的自发辐射寿命分别为0.199和11.8 μs.a3Π0+ 态的辐射寿命比杨传路教授等[38]预测的LiRb分子的辐射寿命(τ = 10.26 μs)要略大. 说明这两种跃迁构建的准闭合循环都能快速的冷却SeH–阴离子.a3Π1(υ′)↔X1Σ+0+(υ″) 跃迁具有非常大的振动分支比R00, 和TeH–阴离子一样, 采用a3Π1(υ′)↔ X1Σ+0+(υ″) 阻禁跃迁对SeH–阴离子进行激光冷却时只需要选取一束主激光来驱动就可以保证跃迁循环的准闭合性, 其散射的光子数目为Nscat = 1/R01+, 理论上可以散射光子数大于2×104个. 所需的主激光波长为501.01 nm. 而a3Π0+(υ′)↔ X1Σ+0+(υ″) 跃迁的振动分支比R00只有0.9484, 故在一束主激光的基础上增加了两束抽运激光来保证其跃迁循环的准闭合性, 其理论上散射光子数可以达到2.5 × 104个. 所需的主激光波长为λ00 = 467.59 nm, 抽运激光波长分别为λ10 = 521.91和λ21 = 526.15 nm. 其冷却方案绘于图4(a)和图4(b)中. 图中的实线表示驱动激光的波长, 虚线表示自发辐射.3.3.2 三能级跃迁
本文计算得到
A1Π1(υ′)↔X1Σ+0+(υ″) 跃迁的f00和R00分别为0.9974和0.9985, 同样其具有很大的自发辐射系数, 为1.36 × 107 s–1, 因此A1Π1 态的自发辐射寿命为73.5 ns, 其比a3Π1 和a3Π0+ 态的辐射寿命要短得多. 采用A1Π1(υ′)↔X1Σ+0+(υ″) 跃迁能更快速地冷却SeH–阴离子. 在构建A1Π1(υ′)↔X1Σ+0+(υ″) 准闭合的循环跃迁时存在中间态a3Π1 和a3Π0+ , 必须考虑中间态是否对跃迁循环准闭合性有影响. 本文采用振动分支损失比(η)来评估. 振动分支损失比可以表示为η=γ/γΣ, (2) 其中, 对于SeH–阴离子, γ1和γ2分别表示
A1Π1↔a3Π1 和A1Π1↔a3Π0+ 跃迁的自发辐射系数. 本文的计算值分别为2.31和1.21 × 10–2 s–1. γΣ表示A1Π1↔X1Σ+0+ 态的自发辐射速率. 这样可以得到η1 = γ1/γΣ < 2.0 × 10–7, η2 = γ2/γΣ < 9.0 × 10–10, 本文的计算结果远远小于YO分子的振动分支损失比(< 4 × 10–4)[2]. 结果表明中间态的存在对激光冷却SeH–阴离子几乎没有影响. 可以选取一束主激光来驱动A1Π1↔X1Σ+0+ 跃迁时, 理论上散射光子数仅为600个. 其主激光波长为459.52 nm. 为了提高冷却效率, 增加了两束抽运激光λ10 = 511.88和λ21 = 512.83 nm, 其理论上散射光子数可达到3.3 × 105个. 相应的三能级跃迁进行激光冷却SeH–阴离子的冷却途径如图5所示.4. 结 论
在MRCI+Q/AVQZ-DK水平下, 本文计算了SeH–阴离子5个Λ-S态和9个Ω态的势能曲线. 计算中考虑了标量相对论效应和CV关联效应. 考虑SOC效应后b3Σ+态的2个分裂态变为了弱束缚态. 所有电子态的光谱常数都是首次报道.
a3Π1(υ′)↔X1Σ+0+(υ″) ,a3Π0+(υ′)↔X1Σ+0+(υ″) 阻禁跃迁和A1Π1(υ′)↔X1Σ+0+(υ″) 三能级跃迁都具有高对角分布的弗兰克-康登因子和振动分支比. 并且a3Π1 ,a3Π0+ 和A1Π1 态的自发辐射寿命都很短, 能快速地冷却SeH–阴离子. 中间态a3Π1 和a3Π0+ 的存在不影响构造A1Π1(υ′)↔X1Σ+0+(υ″) 准闭合的循环跃迁. 驱动a3Π1(υ′)↔X1Σ+0+(υ″) 跃迁只需要一束主激光, 由于振动分支比不够大, 驱动a3Π0+(υ′) ↔X1Σ+0+(υ″) 和A1Π1(υ′)↔ X1Σ+0+(υ″) 跃迁需要增加两束抽运激光. 所有的激光波长都在可见光的范围内. 本文预测SeH–阴离子是适合激光冷却的潜在体系.[1] Shuman E S, Barry J F, DeMille D 2010 Nature 467 820
Google Scholar
[2] Hummon M T, Yeo M, Stuhl B K, Collopy A L, Xia Y, Ye J 2013 Phys. Rev. Lett. 110 143001
Google Scholar
[3] Zhelyazkova V, Cournol A, Wall T E, Matsushima A, Hudson J J, Hinds E A, Tarbutt M R, Sauer B E 2014 Phys. Rev. A 89 053416
Google Scholar
[4] Gao Y, Gao T 2014 Phys. Rev. A 90 052506
Google Scholar
[5] You Y, Yang C L, Wang M S, Ma X G, Liu W W 2015 Phys. Rev. A 92 032502
Google Scholar
[6] Cui J, Xu J G, Qi J X, Dou G, Zhang Y G 2018 Chin. Phys. B 27 103101
Google Scholar
[7] Wan M J, Yuan D, Jin C G, Wang F H, Yang Y J, Yu Y, Shao J X 2016 J. Chem. Phys. 145 024309
Google Scholar
[8] 张云光, 张华, 窦戈, 徐建刚 2017 物理学报 66 233101
Google Scholar
Zhang Y G, Zhang H, Dou G, Xu J G 2017 Acta Phys. Sin. 66 233101
Google Scholar
[9] Xu L, Wei W, Xia Y, Deng L Z, Yin J P 2017 Chin. Phys. B 26 033702
Google Scholar
[10] Yzombard P, Hamamda M, Gerber S, Doser M, Comparat D 2015 Phys. Rev. Lett. 114 213001
Google Scholar
[11] Zhang Q Q, Yang C L, Wang M S, Ma X G, Liu W W 2017 Spectrochim. Acta. Part A 185 365
Google Scholar
[12] Zhang Q Q, Yang C L, Wang M S, Ma X G, Liu W W 2017 Spectrochim. Acta. Part A 182 130
Google Scholar
[13] Zeid I, Abdallah R A, El Kork N, Korek M 2020 Spectrochim. Acta. Part A 224 117461
Google Scholar
[14] Wan M J, Huang D H, Yu Y, Zhang Y G 2017 Phys. Chem. Chem. Phys. 19 27360
Google Scholar
[15] 万明杰, 李松, 金成国, 罗华锋 2019 物理学报 68 063103
Google Scholar
Wan M J, Li S, Jin C G, Luo H F 2019 Acta. Phys. Sin. 68 063103
Google Scholar
[16] Deng B L, Wan M J, Zhao X F, Tang K, Zhang X Q 2020 Spectrochim. Acta. Part A 227 117684
Google Scholar
[17] Stoneman R C, Larson D J 1987 Phys. Rev. A 35 2928
Google Scholar
[18] Brown J M, Fackerell A D 1982 Physica. Scripta. 25 351
Google Scholar
[19] Balasubramanian K, Liao M Z, Han M 1987 Chem. Phys. Lett. 139 551
Google Scholar
[20] Binning Jr R C, Curtiss L A 1990 J. Chem. Phys. 92 1860
Google Scholar
[21] Werner H J, Knowles P J, Knizia G, et al. 2010 MOLPRO, a Package of ab initio Programs (version 2010.1)
[22] Knowles P J, Werner H J 1985 J. Chem. Phys. 82 5053
Google Scholar
[23] Werner H J, Meyer W 1980 J. Chem. Phys. 73 2342
Google Scholar
[24] Werner H J, Knowles P J 1988 J. Chem. Phys. 89 5803
Google Scholar
[25] Knowles P J, Werner H J 1988 Chem. Phys. Lett. 145 514
Google Scholar
[26] Laughoff S R, Davidson E R 1974 Int. J Quantum. Chem. 8 61
Google Scholar
[27] Douglas N, Kroll N M 1974 Ann. Phys. 82 89
Google Scholar
[28] Hess B A 1986 Phys. Rev. A 33 3742
Google Scholar
[29] Berning A, Schweizer M, Werner H J, Knowles P J, Palmieri P 2000 Mol. Phys. 98 1823
Google Scholar
[30] DeYonker N J, Peterson K A, Wilson A K 2007 J. Phys Chem A 111 11383
Google Scholar
[31] Dunning Jr T H 1989 J. Chem. Phys. 90 1007
Google Scholar
[32] Le Roy R J 2007 LEVEL 8.0: a Computer Program for Solving the Radial Schröinger Equation for Bound and Quasibound Levels (Waterloo: University of Waterloo) Chemical Physics Research Report CP-663
[33] Huber K, Herzberg G 1979 Molecular Spectra and Molecular Structure Vol. 4. Constants of Diatomic Molecules (New York: Van Nostrand Reinhold) p586
[34] Moore C E 1971 Atomic Energy Levels (Vol. 1) (Washington, DC: US Govt Printing Office) pp2, 150
[35] Lykke K R, Murray K K, Lineberger W C 1991 Rhys. Rev. A 43 6104
Google Scholar
[36] Hotop H, Lineberger W C 1985 J. Phys. Chem. Ref. Data 14 731
Google Scholar
[37] Lane I C 2015 Phys. Rev. A 92 022511
Google Scholar
[38] You Y, Yang C L, Zhang Q Q, Wang M S, Ma X G, Wang W W 2016 Phys. Chem. Chem. Phys. 18 19838
Google Scholar
-
表 1 Λ-S的光谱常数
Table 1. Spectroscopic parameters of the Λ-S states.
Λ-S态 来源 Re/Å ωe/cm–1 ωeχe/cm–1 Be/cm–1 De/eV Te/cm–1 X1Σ+ ACVQZ-DK 1.4694 2300.77 46.10 7.8507 3.487 0 AVQZ-DK 1.4614 2380.32 45.57 7.9326 3.711 实验[17] 1.4696 a 7.7289 c 1.4806 b a3Π 本文工作 第一势阱 1.4778 2206.52 123.45 7.8428 0.519 20642.90 第二势阱 2.1787 839.87 49.66 3.44016 0.450 24549.11 A1Π 本文工作 第一势阱 1.4726 2373.65 127.14 7.8391 0.734 21240.75 第二势阱 2.2780 437.62 44.07 3.0932 0.147 26997.57 b3Σ+ 本文工作 repulsive 21Σ+ 本文工作 第一势阱 1.6188 1336.45 — 6.1955 0.228 51684.73 第二势阱 4.0808 198.90 9.96 1.0190 0.135 46349.30 注: a 为SeH分子基态的平衡核间距的实验值来源于文献[18]; b为SeH分子基态的平衡核间距的实验值来源于文献[33], 结果不准确; c 为采用最小二乘法得到转动惯量B. 表 2 第VI主簇氢化物阴离子基态的光谱常数
Table 2. Spectroscopic parameters of the ground state X1Σ+ of the Group VI-hydride anions.
表 3 Ω态的离解关系
Table 3. Calculated dissociation relationships of the Ω states.
表 4 Ω电子态的光谱常数
Table 4. Spectroscopic parameters of the Ω states.
Ω态 Re/Å ωe/cm–1 ωeχe/cm–1 Be/cm–1 De/eV Te/cm–1 X1Σ+0+ 1.4694 2301.31 47.01 7.8499 3.395 0 a3Π2 第一势阱 1.4777 2207.22 122.39 7.8416 0.523 19787.17 第二势阱 2.1739 861.02 52.10 3.4081 0.454 23751.54 a3Π1 第一势阱 1.4759 2232.16 111.70 7.8434 0.560 20036.27 第二势阱 2.1822 818.11 55.64 3.3929 0.386 24301.10 a3Π0− 第一势阱 1.4778 2205.83 124.97 7.8485 0.513 21472.52 第二势阱 2.1986 778.28 73.70 3.4048 0.267 25261.96 a3Π0+ 第一势阱 1.4777 2208.03 122.90 7.8422 0.522 21477.12 第二势阱 2.1619 904.15 49.02 3.4355 0.527 25454.22 A1Π1 第一势阱 1.4744 2368.50 144.22 7.8262 0.686 21821.04 第二势阱 排斥态 b3Σ+0− 第一势阱 3.1807 318.89 35.64 1.6988 0.096 28945.41 b3Σ+1 第二势阱 3.2046 239.13 30.94 1.6662 0.066 29184.63 21Σ+0+ 第一势阱 1.6190 1332.50 — 6.1895 0.225 51714.58 第二势阱 4.0800 190.57 8.84 1.0190 0.135 46351.94 表 5
a3Π1↔X1Σ+0+ ,a3Π0+↔X1Σ+0+ 和A1Π1↔ X1Σ+0+ 跃迁的辐射系数Aυ′υ″ 、弗兰克-康登因子fυ′υ″ 和振动分支比Rυ′υ″ Table 5. Emission rates
Aυ′υ″ , Franck-Condon Factorsfυ′υ″ , branching ratiosRυ′υ″ of thea3Π1↔X1Σ+0+ ,a3Π0+↔X1Σ+0+ , andA1Π1↔ X1Σ+0+ transitions.Index a3Π1↔X1Σ+0+ a3Π0+↔X1Σ+0+ A1Π1↔X1Σ+0+ fυ′υ″ f00 0.9949 0.9922 0.9974 f01 0.0047 0.0072 0.0025 f02 0.0004 0.0006 0.0001 f10 0.0051 0.0079 0.0026 f11 0.9541 0.9324 0.9792 f12 0.0337 0.0486 0.0159 Aυ′υ″/s A00 5.02×106 8.02×104 1.36×107 A01 1.88×102 4.28×103 1.87×104 A02 2.81×101 7.48×101 2.00×103 A10 1.10×105 6.50×102 5.79×104 A11 4.13×106 9.13×104 1.32×107 A12 1.32×104 1.57×104 1.45×105 Rυ′υ″ R00 0.99996 0.9484 0.9985 R01 3.7×10–5 0.0506 0.0014 R02 5.6×10–6 0.0009 0.0001 R10 0.02592 0.0060 0.0043 R11 0.9707 0.8394 0.9836 R12 0.0031 0.1446 0.0108 -
[1] Shuman E S, Barry J F, DeMille D 2010 Nature 467 820
Google Scholar
[2] Hummon M T, Yeo M, Stuhl B K, Collopy A L, Xia Y, Ye J 2013 Phys. Rev. Lett. 110 143001
Google Scholar
[3] Zhelyazkova V, Cournol A, Wall T E, Matsushima A, Hudson J J, Hinds E A, Tarbutt M R, Sauer B E 2014 Phys. Rev. A 89 053416
Google Scholar
[4] Gao Y, Gao T 2014 Phys. Rev. A 90 052506
Google Scholar
[5] You Y, Yang C L, Wang M S, Ma X G, Liu W W 2015 Phys. Rev. A 92 032502
Google Scholar
[6] Cui J, Xu J G, Qi J X, Dou G, Zhang Y G 2018 Chin. Phys. B 27 103101
Google Scholar
[7] Wan M J, Yuan D, Jin C G, Wang F H, Yang Y J, Yu Y, Shao J X 2016 J. Chem. Phys. 145 024309
Google Scholar
[8] 张云光, 张华, 窦戈, 徐建刚 2017 物理学报 66 233101
Google Scholar
Zhang Y G, Zhang H, Dou G, Xu J G 2017 Acta Phys. Sin. 66 233101
Google Scholar
[9] Xu L, Wei W, Xia Y, Deng L Z, Yin J P 2017 Chin. Phys. B 26 033702
Google Scholar
[10] Yzombard P, Hamamda M, Gerber S, Doser M, Comparat D 2015 Phys. Rev. Lett. 114 213001
Google Scholar
[11] Zhang Q Q, Yang C L, Wang M S, Ma X G, Liu W W 2017 Spectrochim. Acta. Part A 185 365
Google Scholar
[12] Zhang Q Q, Yang C L, Wang M S, Ma X G, Liu W W 2017 Spectrochim. Acta. Part A 182 130
Google Scholar
[13] Zeid I, Abdallah R A, El Kork N, Korek M 2020 Spectrochim. Acta. Part A 224 117461
Google Scholar
[14] Wan M J, Huang D H, Yu Y, Zhang Y G 2017 Phys. Chem. Chem. Phys. 19 27360
Google Scholar
[15] 万明杰, 李松, 金成国, 罗华锋 2019 物理学报 68 063103
Google Scholar
Wan M J, Li S, Jin C G, Luo H F 2019 Acta. Phys. Sin. 68 063103
Google Scholar
[16] Deng B L, Wan M J, Zhao X F, Tang K, Zhang X Q 2020 Spectrochim. Acta. Part A 227 117684
Google Scholar
[17] Stoneman R C, Larson D J 1987 Phys. Rev. A 35 2928
Google Scholar
[18] Brown J M, Fackerell A D 1982 Physica. Scripta. 25 351
Google Scholar
[19] Balasubramanian K, Liao M Z, Han M 1987 Chem. Phys. Lett. 139 551
Google Scholar
[20] Binning Jr R C, Curtiss L A 1990 J. Chem. Phys. 92 1860
Google Scholar
[21] Werner H J, Knowles P J, Knizia G, et al. 2010 MOLPRO, a Package of ab initio Programs (version 2010.1)
[22] Knowles P J, Werner H J 1985 J. Chem. Phys. 82 5053
Google Scholar
[23] Werner H J, Meyer W 1980 J. Chem. Phys. 73 2342
Google Scholar
[24] Werner H J, Knowles P J 1988 J. Chem. Phys. 89 5803
Google Scholar
[25] Knowles P J, Werner H J 1988 Chem. Phys. Lett. 145 514
Google Scholar
[26] Laughoff S R, Davidson E R 1974 Int. J Quantum. Chem. 8 61
Google Scholar
[27] Douglas N, Kroll N M 1974 Ann. Phys. 82 89
Google Scholar
[28] Hess B A 1986 Phys. Rev. A 33 3742
Google Scholar
[29] Berning A, Schweizer M, Werner H J, Knowles P J, Palmieri P 2000 Mol. Phys. 98 1823
Google Scholar
[30] DeYonker N J, Peterson K A, Wilson A K 2007 J. Phys Chem A 111 11383
Google Scholar
[31] Dunning Jr T H 1989 J. Chem. Phys. 90 1007
Google Scholar
[32] Le Roy R J 2007 LEVEL 8.0: a Computer Program for Solving the Radial Schröinger Equation for Bound and Quasibound Levels (Waterloo: University of Waterloo) Chemical Physics Research Report CP-663
[33] Huber K, Herzberg G 1979 Molecular Spectra and Molecular Structure Vol. 4. Constants of Diatomic Molecules (New York: Van Nostrand Reinhold) p586
[34] Moore C E 1971 Atomic Energy Levels (Vol. 1) (Washington, DC: US Govt Printing Office) pp2, 150
[35] Lykke K R, Murray K K, Lineberger W C 1991 Rhys. Rev. A 43 6104
Google Scholar
[36] Hotop H, Lineberger W C 1985 J. Phys. Chem. Ref. Data 14 731
Google Scholar
[37] Lane I C 2015 Phys. Rev. A 92 022511
Google Scholar
[38] You Y, Yang C L, Zhang Q Q, Wang M S, Ma X G, Wang W W 2016 Phys. Chem. Chem. Phys. 18 19838
Google Scholar
期刊类型引用(1)
1. 郭芮,谭涵,袁沁玥,张庆,万明杰. LiCl~–阴离子的光谱性质和跃迁性质. 物理学报. 2022(04): 40-47 . 百度学术
其他类型引用(0)
Catalog
Metrics
- Abstract views: 6568
- PDF Downloads: 90
- Cited By: 1