摘要: 利用Davidson修正的内收缩多参考组态相互作用(ic-MRCI +
Q )方法, 结合相对论有效芯赝势基(aug-cc-pV5Z-PP)作为Sr原子和相关一致五重基aug-cc-pV5Z为Cl原子的计算基组, 优化计算了单氯化锶(Sr
35 Cl)分子14个低激发电子态的势能曲线和跃迁偶极矩. 为了获得更加精确的光谱参数, 计算中同时引入核价电子相关和相对论效应修正势能曲线. 利用LEVEL 8.0程序拟合修正的势能曲线, 得到相应电子态的光谱常数、振动能级和分子常数等光谱性质, 结果与近来的已获得的理论计算和实验值符合得较好, 同时给出了Franck-Condon因子和辐射寿命等跃迁性质. 这些精确的光谱跃迁特性可为进一步构建Sr
35 Cl分子激光冷却方案提供理论支持.
关键词: Sr35 Cl分子 /
ic-MRCI + Q /
光谱和分子常数 /
振动能级 /
跃迁特性 English Abstract Spectrum and transition characteristics of low excited state of strontium chloride molecule Wu Dong-Lan ,Yuan Jin-Hong ,Wen Yu-Feng ,Zeng Xue-Feng ,Xie An-Dong College of Mathematic and Physical, Jinggangshan University, Ji’an 343009, China Fund Project: Project supported by the National Natural Science Foundation of China (Grant Nos. 11564019, 11147158) and the Science and Technology Project of Jiangxi Provincial Education Department, China (Grant No. GJJ170654). Received Date: 26 September 2018Accepted Date: 11 December 2018Available Online: 01 February 2019Published Online: 05 February 2019 Abstract: Sr35 Cl is a candidate system for laser cooling. The spectrum and transition characteristics are very important for constructing laser cooling schemes. In this paper, the spectral properties are analyzed by using the Davidson's modified internal contraction multi-reference interaction (ic-MRCI + Q ) method, in combination with the relativistic effective core pseudopotential group (aug-cc-pV5Z-PP) as the base group for the calculation of the Sr atom and the related consistent quintile aug-cc-pV5Z as the Cl atom. The potential energy curves and dipole moments of 14 low excited electron states of Sr35 Cl molecule are optimized. In order to obtain more accurate spectral parameters, nuclear valence electron correlation and relativistic effect correction are introduced into the calculation. Using the LEVEL 8.0 program to fit the modified potential energy curves of 5 bound states, the spectral properties such as spectral constants, vibration energy levels, and molecular constants of the corresponding electron states are obtained. The results show that there is a double potential well in ${\rm B}^2 \Sigma^+$ state and the cross phenomena are avoided in ${\rm A}^2 \Pi$ and ${\rm C}^2 \Pi$ , ${\rm B}^2 \Sigma^+$ and $3^2 \Sigma^+$ respectively. The spectrum and molecular constants are in good agreement with the recently obtained theoretical calculations and experimental values except the adiabatic excitation energy. It may be due to the fact that the effect of the interaction of electronic states is taken into account. The transition properties such as Frank-Condon factor and radiation lifetime are also given. It can be seen that the 0-0 band of ${\rm B}^2 \Sigma^+$ ?${\rm X}^2 \Sigma^+$ transition has the largest Franck-Condon factor of 0.861288, and the diagonalization is obvious, which is the condition for laser cooling. The lifetime of ${\rm B}^2 \Sigma^+$ ?${\rm X}^2 \Sigma^+$ transition is 38.89 ns, which is in accordance with the experimental value 39.6 ns ± l.6 ns. These precise spectral transition characteristics may provide theoretical support for further constructing the laser cooling scheme of Sr35 Cl molecule. Keywords: Sr35 Cl molecule /internal contraction multi-reference interaction /spectroscopic and molecular constant /vibration levels /transition characteristics 全文HTML --> --> --> 1.引 言 碱土金属卤化物在天体物理学、催化和高温化学中的应用广泛, 再加上特殊的电子组态引起了许多实验与理论光谱学家的兴趣, 一直是许多电子结构和光谱性质研究的主题[1 ?3 ] . 1928年, Walters和Barratt[4 ] 最早报道了单卤化锶(Sr35 Cl)的光学光谱; 后来由于该类化合物基态(${\rm{X}^2}\Sigma^+ $ )和第一、二激发态${\rm{A}^2}\Pi $ , ${\rm{B}^2} \Sigma^+$ )的势能曲线和核间距非常相似, 导致光谱密集不便于区分. 随着激光器的出现, 许多研究者[5 ?10 ] 精确地分析了${\rm{A}^2}\Pi $ , $ {\rm{B}^2}\Sigma^+$ 到${\rm{X}^2}\Sigma^+$ 态的光谱跃迁. 在这些研究中, 发现${\rm{B}^2}\Sigma^+$ 态的自旋相互作用常数为657 cm?1 , 表明${\rm{B}^2}\Sigma^+$ 与$ {\rm{A}^2}\Pi $ 态之间存在较大的相互作用. Zare等[11 ] 使用直接方法推导了这两个激发态到基态跃迁的微扰系数; Singh等[12 ] 对${\rm{A}^2}\Pi $ —$ {\rm{X}^2}\Sigma^+$ 和${\rm{B}^2}\Sigma^+$ —${\rm{X}^2} \Sigma^+$ 跃迁进行了振动分析; Brinkmann等[13 ] 通过分析散射光谱指出了文献[12 ]中存在的一些错误, 但是也没有成功地分析光谱; Schütze-Pahlmann等[14 ] 采用微波光谱确定了$ {\rm{X}^2}\Sigma^+$ 态的一阶转动常数; Ernst和Schr?der[15 ] 通过多普勒自由极化光谱确定了${\rm{B}^2}\Sigma^+$ —$ {\rm{X}^2}\Sigma^+$ 态跃迁的几个主要Dunham常数; 后来Schreder等[16 ] 采用同样的方法分析了${\rm{A}^2}\Pi $ —${\rm{X}^2} \Sigma^+$ 跃迁的转动结构和微扰, 获得了${\rm{A}^2}\Pi $ 态的转动常数, 同时利用微扰矩阵元分析了微扰对常数的影响; Berg等[17 ] 采用时间分辨激光共振实验获得了${\rm{B}^2}\Sigma^+$ 态的辐射寿命为39.6 ns ± l.6 ns. 通过分析发现, 大部分实验主要集中在较低的两个激发态, 对其他更高激发态的光谱参数与跃迁性质研究不多, 特别是理论计算分析光谱和分子常数的研究较少, 且部分计算结果与实验值有一定的差异, 这可能是由于计算中没有充分考虑能量大小的不一致性、核价电子关联和相对论效应等相关修正, 而这些修正都会影响光谱常数和跃迁性质的计算精度. 这些高精度的光谱性质可以提供广泛的、完整的电子结构信息, 同时对化学反应和分子动力学的理解也非常有意义. 本文利用内收缩多参考组态相互作用方法(ic-MRCI), 优化计算了Sr35 Cl分子的较低的两个解离极限14个激发态在全空间的电子结构和能量. 为了提高计算精度, 计算中考虑了Davison修正、核价电子关联效应和相对论效应修正, 最终得到14个激发态的势能曲线. 基于修正的5个较低束缚态的势能曲线和跃迁偶极矩, 利用LEVEL 8.0程序拟合计算得到相应束缚态的光谱常数、振动能级和分子常数, 同时给出了Franck-Condon因子和辐射寿命等跃迁性质. 这些精确的光谱跃迁特性可为进一步构建SrCl分子激光冷却方案提供理论支持.2.计算方法 本文所有与电子结构和能量有关的计算均使用MOLPRO2012.1程序包完成[18 ] . 经过测试不同的基组, 最终选取相对论有效芯赝势基(aug-cc-pV5Z-PP)作为Sr原子和相关一致五重基(aug-cc-pV5Z)作为Cl原子的计算基组[19 ] . Sr35 Cl分子属于C ∞v 点群. 由于MOLPRO程序自身的限制, 在实际计算中使用C ∞v 点群的阿贝尔子群C 2v , C 2v 群的四个不可约表示(A1 , B1 , B2 , A2 )和C ∞v 点群的不可约表示有如下对应关系: $\Sigma^+ $ = A1 , $\prod $ = B1 + B2 , $\Delta $ = A1 + A2 , $ \Sigma^- $ = A2 . 在计算过程中Sr35 Cl分子的8个轨道(4个a1, 2个b1和2个b2)被确定为活性空间, 它们对应Sr原子的4s4p轨道和Cl原子的3s3p轨道, Sr原子的外层电子5s2 和Cl原子的外层电子3s2 3p5 被放入活性空间进行相关能的计算, 余下的48个电子被冻结在芯轨道. 在C 2v 对称性下, 扫描计算Sr35 Cl分子在核间距R = 0.16—0.80 nm之间的单点能, 步长设为0.01 nm, 其中在平衡位置附近选取步长为0.005 nm. 扫描过程中, 首先使用Hartree-Fock方法计算Sr35 Cl分子基态的单组态波函数, 然后利用态平均全活性空间自洽场方法(SA-CASSCF)优化单组态波函数为多组态波函数[20 ,21 ] , 最后应用内收缩多参考态相互作用方法(ic-MRCI)计算Sr35 Cl分子束缚态的能量[22 ,23 ] . 为了进一步提高精度, 计算中考虑Davidson修正(+Q ), 核价电子关联效应和相对论标量效应修正, 最后得到了Sr35 Cl分子最低的14个激发态的势能曲线, 其中5个束缚态, 9个排斥态. 基于修正的5个束缚态的势能曲线和跃迁偶极矩, 利用LEVEL 8.0程序[24 ] 求解径向Schr?dinger方程获得了全部束缚态的光谱常数和分子常数, 包括离解能D e 、平衡核间距R e 、绝热激发能T e 、谐性振动常数$\omega_{\rm e}$ 和非谐性振动常数$\omega_{\rm e}$ $\chi_{\rm e}$ , 转动常数B e 和振转耦合常数$\alpha_{\rm e}$ , 以及振动能级Gv 、惯性转动常数Bv 和离心畸变常数Dv . 同时给出了Franck-Condon因子和辐射寿命等跃迁性质.3.结果与讨论 23.1.束缚态的势能曲线和光谱常数 -->3.1.束缚态的势能曲线和光谱常数 利用内收缩的多参考组态相互作用方法(ic-MRCI + Q )结合相关一致五重基aug-cc-pV5Z和相对论有效芯赝势基aug-cc-pV5Z-PP优化计算Sr35 Cl分子的电子结构. 分析对应解离极限Sr(1 Sg ) + Cl(2 P u ), Sr(3 P u ) + Cl(2 Pu )和离子对Sr+ (2 Sg ) + Cl? (1 Sg )的14个激发电子态, 同时扫描单点能获得对应电子态的势能曲线, 如图1 所示. 从图1 可看出只有${\rm{X}^2}\Sigma^+ $ , ${\rm{A}^2}\Pi$ , $ {\rm{B}^2}\Sigma^+$ , $ {\rm{C}^2}\Pi$ 和$ {\rm{3}^2}\Sigma^+$ 为束缚态, 其余都为排斥态. 其中每个束缚态在平衡键长位置处能量最低, ${\rm{B}^2}\Sigma^+$ 态出现了双势阱, 且${\rm{A}^2}\Pi$ 和${\rm{C}^2}\Pi$ 在R = 0.27 nm处, ${\rm{B}^2}\Sigma^+$ 和$ {\rm{3}^2}\Sigma^+$ 在R = 0.251 nm处出现了避免交叉, ${\rm{C}^2}\Pi$ 和${\rm{B}^2}\Sigma^+ $ 在R = 0.228 nm和0.298 nm处相交, 这些势能曲线的特征类似于相同类型的分子体系CaCl[25 ] . 利用LEVEL 8.0拟合修正的势能曲线, 得到相应束缚态的光谱常数列于表1 , 同时还列出了各个电子态在平衡位置R e 附近主要的电子组态分布及权重. 图 1 Sr35 Cl分子14个激发态的势能曲线 Figure1. Potential energy curves of 14 excited states of Sr35 Cl. Λ-S态 T e /cm?1 R e /nm $\omega $e /cm?1 $\omega $e $\chi $e /cm?1 B e /cm?1 104 α e /cm?1 D e /eV R e 附近主要电子组态(%) $ {{\rm{X}}^{2}}\Sigma^+$ 0.0 0.2575 309.78 0.8682 0.1016 4.131 3.703 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ α 9$\text{σ}$ 0 3$\text{π}$2 4$\text{π}$2 (78.8) 6${\rm{\sigma }}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 0 9$\text{σ}$ α 3$\text{π}$2 4$\text{π}$2 (7.5) 理论[26 ] 0.0 0.255 313 0.93 0.1037 — — 实验[15 ] 0.0 0.257 302[27 ] 0.95[27 ] — — — 实验[16 ] 0.0 — 302.448 ?0.9502 0.1016 4.524 — 实验[12 ] 0.0 — 302.3 0.950 — — — $ {{\rm{A}}^{2}}\Pi$ 15779.16 0.2518 330.69 — 0.1061 2.013 1.673 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 0 9$\text{σ}$ 0 3$\text{π}$ααβ 4$\text{π}$2 (85.7) 理论[26 ] 14730 0.252 323 0.95 0.1055 — — 实验[15 ] 14818 0.255 309[27 ] 0.98[27 ] — — — 实验[16 ] 14966.727 — 309.625 0.996 0.1030 4.606 — $ {{\rm{B}}^{2}}\Sigma^+$ 16612.74 0.2538 318.67 0.4874 0.1043 2.599 1.937 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ α 9$\text{σ}$ 0 3$\text{π}$2 4$\text{π}$2 (78.8) 6${\rm{\sigma }}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 0 9$\text{σ}$ α 3$\text{π}$2 4$\text{π}$2 (7.6) 理论[26 ] 15714 0.253 319 0.99 0.1055 — — 实验[15 ] 15719 0.255 306[27 ] 0.98[27 ] 0.1030[16 ] — — 实验[12 ] 15719.5 — 306.4 0.98 — — — — $ {{\rm{C}}^{2}}\Pi$ 33532.99 0.3477 425.57 15.5691 0.0554 ?11.932 1.546 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 0 9$\text{σ}$ 0 3$\text{π}$ααβ 4$\text{π}$2 (59.7) 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 2 9$\text{σ}$ 0 3$\text{π}$2 4$\text{π}$2 (13.4) 6$\text{σ}$ 2 7$\text{σ}$ α 8$\text{σ}$ β 9$\text{σ}$ 0 3$\text{π}$ααβ 4π2 (9.8) 6$\text{σ}$ 2 7$\text{σ}$ α 8$\text{σ}$ α 9$\text{σ}$ 0 3$\text{π}$αββ 4$\text{π}$2 (2.7) 理论[26 ] 26688 0.259 278 0.83 0.1095 — — 实验[27 ] 26099 — 270 0.72 — — — $ {{\rm{3}}^{2}}\Sigma^+$ 36625.55 0.3519 392.49 10.4544 0.0544 ?12.567 1.140 6$\text{σ}$ 2 7$\text{σ}$ α 8$\text{σ}$ 2 9$\text{σ}$ 0 3$\text{π}$2 4$\text{π}$2 (51.0) 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 0 9$\text{σ}$ α 3$\text{π}$2 4$\text{π}$2 (16.4) 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ α 9$\text{σ}$ 0 3$\text{π}$ααβ 4$\text{π}$2 (4.6) 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ α 9$\text{σ}$ 0 3$\text{π}$β 4$\text{π}$ααβ (4.6) 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ α 9$\text{σ}$ 0 3$\text{π}$2 4$\text{π}$2 (2.5) 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ α 9$\text{σ}$ 0 3$\text{π}$αββ 4$\text{π}$2 (1.6) 6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 0 9$\text{σ}$ α 3$\text{π}$α 4$\text{π}$αββ (1.6) 6$\text{σ}$ 2 7$\text{σ}$ α 8$\text{σ}$ 0 9$\text{σ}$ 0 3$\text{π}$4 4$\text{π}$2 (1.5) 6$\text{σ}$ 2 7$\text{σ}$ α 8$\text{σ}$ 0 9$\text{σ}$ 0 3$\text{π}$2 4$\text{π}$4 (1.5) 理论[26 ] 27979 0.248 358 1.01 0.1095 — — 实验[15 ] 28822 — 344[27 ] 1.04[27 ] 0.1030[16 ] — —
表1 5个束缚态的光谱常数Table1. Spectroscopic constants of the 5 bound states. 从图1 和表1 可看出, 基态${\rm X}^2\Sigma^+ $ 为一个深势阱态, 阱深为3.703 eV, 说明基态是一个比较稳定的电子态. 在平衡位置附近主要的电子组态为6${\text{σ}}^2$ 7${\text{σ}}$ 2 8${\text{σ}}$ α 9${\text{σ}}$ 0 3$\text{π}$ 2 4$ \text{π}$ 2 和6${\text{σ}}$ 2 7${\text{σ}}$ 2 8${\text{σ}}$ 0 9${\text{σ}}$ α 3$\text{π}$ 2 4$ \text{π}$ 2 , 表明了使用多组态方法的必要性[28 ?30 ] . 通过比较发现, 本文的光谱常数与近来获得的理论计算值[26 ] 较接近, R e , $\omega_{\rm e}$ 和B e 的偏差分别为0.0015 nm, 0.52 cm?1 和0.0112 cm?1 . 与相同类型分子体系相比, 离解能小于BeCl (94.31 kcal/mol = 4.090 eV)[28 ] 和CaCl (34283.5 cm?1 = 4.251 eV)[25 ] , 平衡键长大于BeCl (0.1792 nm)[28 ] , MgCl (0.2214 nm)[29 ] 和CaCl (0.2506 nm)[25 ] . 这些对比说明Sr35 Cl分子相比BeCl, MgCl和CaCl分子体系, 稳定性更差, 更容易离解. 对于第一激发态${\rm A}^2 \Pi$ , 与基态属于同一个解离极限Sr(1 Sg ) + Cl(2 Pu ), 但离解能为1.673 eV, 结果小于基态3.703 eV, 这说明基态更加稳定. 平衡位置附近主要的电子组态为6$\text{σ}$ 2 7$\text{σ}$ 2 8$\text{σ}$ 0 9$\text{σ}$ 0 $3{\text{π}}^{{\text{ααβ}}}4{\text{π}}^2$ , 表明有一个电子从${\rm X}^2\Sigma^+$ 激发跃迁到${\rm A}^2\Pi$ , 即8$\text{σ}$ (9$\text{σ}$ )→3$ \text{π}$ . $ {\rm A}^2\Pi$ 的绝热激发能为15779.16 cm?1 , 与实验值14966.727 cm?1 [16 ] 有一定偏差, 相对偏差为5.4%, 这可能是由于考虑了与${\rm C}^2 \Pi$ 避免交叉现象, 导致${\rm A}^2\Pi$ 的势能曲线在R = 0.348 nm处有一个势垒, 约为30951.13 cm?1 . 同时也使光谱常数$\omega_{\rm e} $ , R e 和B e 与实验结果[16 ] 有一定相差, 但与近来的理论计算[26 ] 符合得很好, 偏差分别为7.69 cm?1 , 0.0002 nm和0.0006 cm?1 . 特别值得注意的是第二激发态B2 $ \Sigma^+$ , 势能曲线出现了双势阱, 且该激发态对应离子对解离极限Sr+ (2 Sg ) + Cl? (1 Sg ), 绝热激发能为16612.74 cm?1 , 与实验值15719.5 cm?1 [12 ] 也有偏差, 相对偏差为5.7%, 这可能也是由于${\rm{B}^2}\Sigma^+$ 与${\rm{3}^2}\Sigma^+ $ 出现了避免交叉现象, 导致势能曲线发生变化, 致使拟合的光谱常数$\omega $ e , R e 和B e 与实验结果[15 ] 也有一定差异. 但与近来的理论计算[26 ] 结果也符合得较好, 偏差仅分别为0.33 cm?1 , 0.0008 nm和0.0012 cm?1 . 从图1 发现第三激发态${\rm C}^2\Pi$ , 由于与${\rm A}^2 \Pi$ 产生了避免交叉现象, 最低能量正好位于避免交叉处R = 0.348 nm, 其能量大小为31191.10 cm?1 , 与${\rm A}^2\Pi$ 能隙为239.97 cm?1 . 离解能为1.546 eV, 不同的是${\rm C}^2\Pi$ 对应的是第二解离极限Sr(3 Pu ) + Cl(2 Pu ). 从图1 还可以发现${\rm A}^2\Pi$ 与$ {\rm B}^2\Sigma^+$ 的势能曲线在0.261 nm和0.340 nm位置存在两处交叉, 本文给出了${\rm C^2}\Pi$ 和$3^2\Sigma^+ $ 激发态的光谱常数, 结果与近来获得的理论值相差较大[26 ] , 这可能是由于电子态之间相互作用的影响所致. 综上分析发现, 本文的光谱常数与实验值和近来的理论计算值有一定差异, 可能是由于本文考虑了避免交叉现象, 即考虑了态与态之间的电子关联相互作用, 这会使拟合得到的光谱常数出现一定偏差, 但当进一步考虑自旋-轨道耦合效应后, 这些势能曲线会发生变化, 从而得到更加精确的光谱常数. 限于篇幅, 后续我们将报道详细的自旋-轨道耦合效应对光谱常数和跃迁特性的影响. 23.2.束缚态的振动能级和分子常数 -->3.2.束缚态的振动能级和分子常数 利用MRCI + Q /aug-cc-pV5Z + CV + DK (aug-cc-pV5Z-PP)优化计算得到5个束缚态的势能曲线, 借助于LEVEL 8.0程序包[24 ] , 找到了J = 0时5个束缚态${\rm X}^2\Sigma^+$ , ${\rm A}^2 \Pi$ , ${\rm C}^2\Pi$ , ${\rm B}^2\Sigma^+$ 和$ 3^2\Sigma^+$ 的全部振动态. 表2 分别列出了基态和更高激发态前10个振动态的振动能级Gv , 惯性转动常数Bv 和离心畸变常数Dv . 通过前期调研发现, 文献[16 ]中${\rm A}^2\Pi$ 的0-0分子常数的实验值分别为14970.22, 0.1027584和4.579 × 108 cm?1 , 与本文的计算值符合得较好, 相对误差仅为6.2%, 3.3%和2.6%, 说明采用高精度计算方法和相关修正获得激发态的分子常数是可靠的. v 0 1 2 3 4 5 6 7 8 9 $ {{\rm{X}}^{2}}\Sigma^+$ Gv /cm?1 0 308.69 615.44 919.34 1220.80 1520.91 1820.07 2118.05 2414.53 2709.33 Bv /cm?1 0.101389 0.100964 0.100567 0.100209 0.099815 0.099369 0.098910 0.098467 0.098043 0.097629 108 Dv /cm?1 4.351839 4.364577 4.467493 4.497514 4.362226 4.242542 4.234360 4.286016 4.334124 4.353049 $ {{\rm{A}}^{2}}\Pi$ Gv /cm?1 15953.74 16281.83 16598.43 16930.62 17269.01 17605.65 17940.04 18272.59 18603.52 18933.04 Bv /cm?1 0.106234 0.106462 0.106095 0.105137 0.104604 0.104235 0.103859 0.103475 0.103082 0.102679 108 Dv /cm?1 4.461909 5.111509 3.253476 3.332790 4.007870 4.130323 4.085546 4.040997 3.976138 3.938946 $ {{\rm{B}}^{2}}\Sigma^+$ Gv /cm?1 16777.58 17097.77 17402.46 17705.63 18020.44 18339.16 18656.69 18972.44 19286.60 19599.30 Bv /cm?1 0.104260 0.104168 0.104392 0.103697 0.102831 0.102297 0.101900 0.101515 0.101130 0.100741 108 Dv /cm?1 4.416717 5.342062 4.947374 3.120903 3.500412 4.129428 4.273733 4.244199 4.206411 4.139361 $ {{\rm{C}}^{2}}\Pi$ Gv /cm?1 33822.25 34298.68 34607.21 34880.81 35124.57 35352.15 35566.01 35770.85 35971.25 36169.15 Bv /cm?1 0.056049 0.057354 0.058912 0.059995 0.060974 0.062071 0.063032 0.063843 0.064515 0.065071 108 Dv /cm?1 0.293093 1.069720 1.371596 2.128287 2.433073 3.023576 3.386636 3.382563 3.426546 3.529599 $ {{\rm{3}}^{2}}\Sigma^+$ Gv /cm?1 36879.56 37294.02 37638.56 37944.07 38221.43 38480.16 38726.99 38965.82 39199.02 39428.18 Bv /cm?1 0.055079 0.056481 0.057758 0.059028 0.060242 0.061377 0.062432 0.063414 0.064332 0.065193 108 Dv /cm?1 0.363052 0.690909 1.002658 1.373820 1.690530 1.922724 2.122355 2.300506 2.453955 2.583320
表2 Sr35 Cl分子$ {{\rm{X}}^{2}}\Sigma^+$ , $ {{\rm{A}}^{2}}\Pi$ , $ {{\rm{B}}^{2}}\Sigma^+$ , $ {{\rm{C}}^{2}}\Pi$ 和$ {{\rm{3}}^{2}}\Sigma^+$ 的Gv , Bv 和Dv 值Table2. The Gv , Bv and Dv of $ {{\rm{X}}^{2}}\Sigma^+$ , $ {{\rm{A}}^{2}}\Pi$ , $ {{\rm{B}}^{2}}\Sigma^+$ , $ {{\rm{C}}^{2}}\Pi$ and $ {{\rm{3}}^{2}}\Sigma^+$ states of Sr35 Cl. 23.3.束缚态的电偶极矩 -->3.3.束缚态的电偶极矩 由于偶极矩反映了分子的成键性质, 本文对Sr35 Cl分子的5个较低束缚态的电偶极矩进行了理论计算, 电偶极矩随着分子核间距R 的变化曲线见图2 . 图 2 Sr35 Cl分子5个束缚态的电偶极矩 Figure2. Permanent dipole moments of 5 bound states of Sr35 Cl. 从图2 可以看出, ${\rm B}^2\Sigma^+$ 的电偶极矩随着核间距的增大而变化, 当增大到0.57 nm时, $ {\rm B}^2\Sigma^+$ 的偶极矩达到最大值(4.907 a.u., 约为7.5 Debye), 表明分子离子特性较强, 与前面分析其对应的离子对解离极限的结果是一致的; 至于其他4个电子束缚态, 当核间距增加较大时, 偶极矩趋向于零, 表明对应的解离极限为中性原子. 在核间距为0.345 nm附近, 具有相同对称性的${\rm A}^2\Pi$ 和${\rm C}^2\Pi$ 出现了避免交叉现象(与势能曲线中的一致), 导致了它们的偶极矩极性反转. 同样$ {\rm B}^2\Sigma^+$ 和$3^2\Sigma^+$ 在0.353 nm附近的偶极矩也出现了极性反转现象, 这也是由于这两个态出现了避免交叉引起的. 这些变化情况与文献[29 ]中相同类型的分子MgCl变化趋势相同. 23.4.束缚态的跃迁偶极矩、Franck-Condon因子和辐射寿命 -->3.4.束缚态的跃迁偶极矩、Franck-Condon因子和辐射寿命 图3 给出了最低五个束缚态的跃迁偶极矩的绝对值随键长的变化. 从图3 发现, 在Franck-Condon区域, ${\rm C}^2\Pi$ —$ {\rm X}^2\Sigma^+$ 跃迁的偶极矩的绝对值明显大于${\rm A}^2\Pi$ —${\rm X}^2\Sigma^+$ 跃迁的数值. 当核间距越来越大, 最终五个束缚态的跃迁偶极矩都趋于0. 表3 列出了利用LEVEL软件程序包计算最低五个束缚态的Franck-Condon因子, 可看出${\rm B}^2\Sigma^+$ —$ {\rm X}^2\Sigma^+$ 跃迁的0-0带具有最大的Franck-Condon因子为0.861288, 且对角化比较明显. 这正好符合激光冷却分子体系的条件, 因此Sr35 Cl分子可以作为激光冷却的候选分子体系. 其次较大的为${\rm A}^2\Pi$ —${\rm X}^2\Sigma^+$ 跃迁的0-0带因子为0.656888, 但是没有明显的对角化; 而${\rm C}^2\Pi$ —${\rm X}^2\Sigma^+$ 和$3^2\Sigma^+$ —${\rm X}^2\Sigma^+$ 的因子非常小, 这可能是由于两者的平衡核间距较大的原因. v ′′ = 0 1 2 3 4 5 6 7 8 9 $ {{\rm{A}}^{2}}\Pi$—$ {{\rm{X}}^{2}}\Sigma^+ $ v ′ = 0 0.656888 0.266608 0.062027 0.011563 0.002170 0.000520 0.000163 0.000048 0.000007 0.000000 1 0.272308 0.192420 0.320947 0.150685 0.466479 0.012567 0.033544 0.008655 0.000174 0.000137 2 0.061741 0.365100 0.012591 0.236003 0.200935 0.086356 0.027568 0.007557 0.001791 0.000312 3 0.008378 0.145153 0.330434 0.013476 0.132392 0.200767 0.112395 0.041433 0.012093 0.002903 4 0.000641 0.027499 0.211170 0.231162 0.064589 0.062086 0.184983 0.134977 0.058263 0.018761 5 0.000040 0.003015 0.053925 0.25428 0.132467 0.115528 0.018774 0.155878 0.151198 0.077211 6 0.000004 0.000182 0.008221 0.084623 0.271397 0.058393 0.151118 0.000803 0.117380 0.157174 7 0.000000 0.000021 0.006253 0.016398 0.117920 0.263055 0.014702 0.164939 0.005113 0.077068 8 0.000000 0.000000 0.000052 0.001674 0.027547 0.151418 0.233623 0.000068 0.157208 0.024940 9 0.000000 0.000000 0.000007 0.000118 0.003627 0.041888 0.181283 0.190225 0.008630 0.133109 $ {{\rm{B}}^{2}}\Sigma^+$—$ {{\rm{X}}^{2}}\Sigma^+ $ v ′ = 0 0.861288 0.125494 0.011927 0.001065 0.000142 0.000047 0.000025 0.000091 0.000000 0.000000 1 0.129692 0.603795 0.220001 0.038332 0.006251 0.001365 0.000411 0.000120 0.000019 0.000000 2 0.008650 0.241321 0.360072 0.284163 0.081308 0.018646 0.004456 0.001109 0.000226 0.000018 3 0.000365 0.027661 0.352567 0.179385 0.288215 0.112584 0.030079 0.007227 0.001600 0.000265 4 0.000002 0.001707 0.051462 0.411903 0.083378 0.268533 0.131641 0.039281 0.009668 0.002055 5 0.000001 0.000017 0.003843 0.078423 0.421433 0.034730 0.246976 0.148790 0.049663 0.012870 6 0.000000 0.000000 0.000117 0.006341 0.108918 0.407382 0.009402 0.221296 0.163711 0.061334 7 8 0.000000 0.000000 0.000003 0.000002 0.000001 0.000007 0.000333 0.000002 0.009581 0.000763 0.141252 0.014168 0.381199 0.173321 0.000074 0.346687 0.190975 0.003801 0.174597 0.158106 9 0.000000 0.000000 0.000000 0.000000 0.000000 0.001272 0.020484 0.203605 0.306437 0.017104 $ {{\rm{C}}^{2}}\Pi$—$ {{\rm{X}}^{2}}\Sigma^+ $ v ′ = 0 0.000002 0.000057 0.000400 0.002491 0.013257 0.048381 0.121107 0.211164 0.245963 0.199742 1 0.000028 0.000755 0.003923 0.016478 0.054711 0.112499 0.129055 0.063035 0.000478 0.054393 2 0.000216 0.004536 0.018057 0.052861 0.110648 0.119527 0.040456 0.002385 0.059184 0.058927 3 0.000894 0.014294 0.042688 0.083192 0.097065 0.033697 0.003920 0.060262 0.039865 0.000817 $ {{\rm{3}}^{2}}\Sigma^+$—$ {{\rm{X}}^{2}}\Sigma^+ $ v ′ = 0 0.000005 0.000160 0.001014 0.005583 0.025973 0.082035 0.174625 0.263266 0.237744 0.147940 1 0.000072 0.001726 0.007985 0.028959 0.080861 0.133268 0.106890 0.018946 0.019348 0.137924 2 0.000444 0.008122 0.028147 0.068699 0.114064 0.085695 0.008341 0.025544 0.075168 0.024261 3 0.001887 0.025614 0.064048 0.096749 0.076464 0.007949 0.023664 0.062244 0.010184 0.021472 4 0.005823 0.057375 0.098715 0.079207 0.013459 0.015420 0.056215 0.010480 0.018687 0.044825 5 0.013104 0.091074 0.099009 0.027450 0.005603 0.051927 0.016773 0.012448 0.040966 0.001448 6 0.023467 0.109148 0.062336 0.000039 0.043234 0.033313 0.003366 0.042212 0.004752 0.023015 7 0.036552 0.106016 0.020316 0.019085 0.055807 0.001734 0.035028 0.017051 0.012699 0.031310 8 0.052050 0.085867 0.000431 0.051861 0.028088 0.011562 0.037256 0.000833 0.036176 0.001957 9 0.069548 0.057189 0.008529 0.060907 0.001791 0.038489 0.008713 0.024666 0.015229 0.013407
表3 Sr35 Cl分子$ {{\rm{A}}^{2}}\Pi$ —$ {{\rm{X}}^{2}}\Sigma^+ $ , $ {{\rm{B}}^{2}}\Sigma^+$ —$ {{\rm{X}}^{2}}\Sigma^+ $ , $ {{\rm{C}}^{2}}\Pi$ —$ {{\rm{X}}^{2}}\Sigma^+ $ 和$ {{\rm{3}}^{2}}\Sigma^+$ —$ {{\rm{X}}^{2}}\Sigma^+ $ 跃迁的Franck-Condon因子Table3. The Franck-Condon factors of the transitions $ {{\rm{A}}^{2}}\Pi$ ?$ {{\rm{X}}^{2}}\Sigma^+ $ , $ {{\rm{B}}^{2}}\Sigma^+$ ?$ {{\rm{X}}^{2}}\Sigma^+ $ , $ {{\rm{C}}^{2}}\Pi$ ?$ {{\rm{X}}^{2}}\Sigma^+ $ and $ {{\rm{3}}^{2}}\Sigma^+$ ?$ {{\rm{X}}^{2}}\Sigma^+ $ . 图 3 Sr35 Cl分子5个束缚态的跃迁偶极矩 Figure3. Transition dipole moments of 5 bound states of Sr35 Cl. 基于计算得到的激发态-基态的跃迁偶极矩、电子态的势能曲线和振动能级, 本文计算了最低五个束缚激发态振动能级的辐射寿命. 激发态振动能级v ′的辐射寿命计算公式如下[33 , 34 ] : 式中${\tau _{v^{\prime}}} $ 为激发态振动能级v ′的辐射寿命, Av′ 为激发态v ′振动能级和基态v ′′振动能级之间的爱因斯坦系数, $\Delta$ Ev′,v′′ 是激发态振动能级v ′和基态振动能级v ′′之间的能量差(单位为cm?1 ), TDM 为平均跃迁偶极矩(激发电子态的振动能级对应的经典回转点区域内的平均值, 单位为a.u.), qv′,v′′ 为Franck-Condon因子(FCFs). 根据上式计算的四个激发态到基态跃迁振动带的辐射寿命, 其中32 $\Sigma^+$ —${\rm X}^2 \Sigma^+$ 跃迁态的寿命非常大, 实验上不容易观测到, 表4 仅列出了${\rm A}^2\Pi$ —${\rm X}^2\Sigma^+$ , $ {\rm B}^2\Sigma^+$ —${\rm X}^2\Sigma^+$ 和C2 $\Pi$ —X2 $ \Sigma^+$ 跃迁态的辐射寿命. 本文给出的${\rm B}^2 \Sigma^+$ —${\rm X}^2\Sigma^+$ 跃迁的v ′ = 1寿命38.89 ns与实验研究给出的39.6 ns ± l.6 ns值符合较好[17 ] , 相对误差仅为1.8%. ${\rm A}^2\Pi$ —${\rm X}^2\Sigma^+ $ 和${\rm C}^2\Pi$ —${\rm X}^2\Sigma^+ $ 跃迁的v ′ = 1寿命分别为31.35 ns和26.01 ns, 与文献测量值31.3 ns ± 2.7 ns和26.1 ns ± l.9 ns都比较接近[35 ] . 实验研究尚未给出$3^2\Sigma^+$ —${\rm X}^2\Sigma^+$ 跃迁态振动能级的辐射寿命, 这可能是寿命较大, 不方便测量的原因. Transition Radiative lifetimes/ns v ′ = 0 v ′ = 1 v ′ = 2 $ {{\rm{A}}^{2}}\Pi $—$ {{\rm{X}}^{2}}\Sigma^+ $ 31.23 31.35 31.56 $ {{\rm{B}}^{2}}\Pi $—$ {{\rm{X}}^{2}}\Sigma^+ $ 38.83 38.89 39.12 $ {{\rm{C}}^{2}}\Pi$—$ {{\rm{X}}^{2}}\Sigma^+ $ 25.92 26.01 26.18
表4 Sr35 Cl分子$ {{\rm{A}}^{2}}\Pi $ —$ {{\rm{X}}^{2}}\Sigma^+ $ , $ {{\rm{B}}^{2}}\Sigma^+ $ —$ {{\rm{X}}^{2}}\Sigma^+ $ 和$ {{\rm{C}}^{2}}\Pi $ —$ {{\rm{X}}^{2}}\Sigma^+ $ 跃迁的辐射寿命Table4. Radiative lifetimes of the transitions $ {{\rm{A}}^{2}}\Pi $ ?$ {{\rm{X}}^{2}}\Sigma^+ $ , $ {{\rm{B}}^{2}}\Sigma^+ $ ?$ {{\rm{X}}^{2}}\Sigma^+ $ and $ {{\rm{C}}^{2}}\Pi $ ?$ {{\rm{X}}^{2}}\Sigma^+ $ . 4.结 论 利用高精度的内收缩多参考组态相互作用方法(ic-MRCI), 考虑多种相关修正研究了Sr35 Cl分子对应最低三个解离极限的激发态电子结构、势能曲线和偶极矩等. 结果表明14个低激发态分别对应三个不同解离极限, ${\rm B}^2\Sigma^+$ 对应离子对解离极限, 其中五个束缚态都存在势阱, 且以基态${\rm X}^2\Sigma^+ $ 的势阱最深, 说明这个态在实验上更稳定, 更容易在实验中观测到. ${\rm B}^2\Sigma^+$ 出现了双势阱, ${\rm A}^2\Pi$ 和${\rm C}^2 \Pi$ , ${\rm B}^2\Sigma^+ $ 和$3^2\Sigma^+$ 出现了避免交叉现象, ${\rm B}^2\Sigma^+$ 和${\rm A}^2 \Pi$ 发生了两次相交, 当考虑了自旋-轨道效应后, 这两个交叉点就会产生避免交叉现象而分开. 通过分析束缚态的电偶极矩, 发现有些地方出现了偶极矩反转, 这也正是由于避免交叉现象引起的, 与势能函数曲线得到的结论一致. 本文给出了五个束缚态的光谱常数, 与近来的理论计算和实验值都符合得较好; 同时还给出了更多激发态的振动能级和分子常数, 部分结果与实验值的百分误差都在7%以内; 最后给出了Franck-Condon因子和辐射寿命等跃迁性质, 获得了比较完善的激发态光谱性质. 这些精确的光谱跃迁特性可为进一步构建Sr35 Cl分子激光冷却方案提供理论支持.