搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

单氯化锶分子低激发态的光谱及跃迁特性

伍冬兰 袁金宏 温玉锋 曾学锋 谢安东

单氯化锶分子低激发态的光谱及跃迁特性

伍冬兰, 袁金宏, 温玉锋, 曾学锋, 谢安东
PDF
HTML
导出引用
  • 利用Davidson修正的内收缩多参考组态相互作用(ic-MRCI + Q)方法, 结合相对论有效芯赝势基(aug-cc-pV5Z-PP)作为Sr原子和相关一致五重基aug-cc-pV5Z为Cl原子的计算基组, 优化计算了单氯化锶(Sr35Cl)分子14个低激发电子态的势能曲线和跃迁偶极矩. 为了获得更加精确的光谱参数, 计算中同时引入核价电子相关和相对论效应修正势能曲线. 利用LEVEL 8.0程序拟合修正的势能曲线, 得到相应电子态的光谱常数、振动能级和分子常数等光谱性质, 结果与近来的已获得的理论计算和实验值符合得较好, 同时给出了Franck-Condon因子和辐射寿命等跃迁性质. 这些精确的光谱跃迁特性可为进一步构建Sr35Cl分子激光冷却方案提供理论支持.
      通信作者: 伍冬兰, wudonglan1216@sina.com
    • 基金项目: 国家自然科学基金(批准号: 11564019, 11147158)和江西省教育厅科学技术研究项目(批准号: GJJ170654)资助的课题.
    [1]

    Micheli A, Brennen G, Zoller P 2006 Nat. Phys. 2 341

    [2]

    Baron J, Campbell W C, DeMille D, Doyle J M, Gabrielse G, Gurevich Y V, Hess P W, Hutzler N R, Kirilov E, Kozyryev I, O'Leary B R, Panda C D, Parsons M F, Petrik E S, Spaun B, Vutha A C, West A D 2014 Science 343 269

    [3]

    Krems R V 2008 Phys. Chem. Chem. Phys. 10 4079

    [4]

    Walters O H, Barratt S 1928 Proc. R. Soc. London Ser. A 118 120

    [5]

    Nakagawa J, Domaille P J, Steimle T C, Harris D O 1978 J. Mol. Spectrosc. 70 374

    [6]

    Dulick M, Bernath P F, Field R W 1980 Can. J. Phys. 58 703

    [7]

    Domaille P J, Steimle T C, Wong N B, Harris D O 1977 J. Mol. Spectrosc. 65 354

    [8]

    Berg L E, Klynning L, Martin H 1980 Phys. Scr. 21 173

    [9]

    Bernath P F, Field R W, Pinchemel B, Lefebvre Y, Schamps J 1981 J. Mol. Spectrosc. 88 175

    [10]

    Reisner D E, Bernath P F, Field R W 1981 J. Mol. Spectrosc. 89 107

    [11]

    Zare R N, Schmeltekopf A L, Harrop W J, Albritton D L 1973 J. Mol. Spectrosc. 46 37

    [12]

    Singh J, Nair K P R, Upadhya K N, Rai D K 1970 Opt. Pure Appl. 3 76

    [13]

    Brinkmann U, Schmidt V H, Telle H 1982 Chern. Phys. 64 19

    [14]

    Schiitze Pahlmann H V, Ryzlewicz C H, Hoeft J, Torring T 1982 Chern. Phys. Lett. 93 74

    [15]

    Ernst W E, Schröder J O 1984 J. Chem. Phys. 81 136

    [16]

    Schruder J O, Zeller B, Ernst W E 1988 J. Mol. Spectrosc. 127 255

    [17]

    Berg L E, Royen P, Weijnitz P 1990 Mol. Phys. 69 385

    [18]

    Werner H J, Knowles P J, Knizia G, et al. 2012 MOLPRO, Version 2012.1, A Package of ab initio Programs

    [19]

    Peterson K A, Figgen D, Goll E, Stoll H, Dolg M 2003 J. Chem. Phys. 119 11099

    [20]

    Werner H J, Knowles P J 1985 J. Chem. Phys. 82 5053

    [21]

    Knowles P J, Werner H J 1985 Chem. Phys. Lett. 115 259

    [22]

    Werner H J, Knowles P J 1988 J. Chem. Phys. 89 5803

    [23]

    Knowles P J, Werner H J 1988 Chem Phys Lett. 145 514

    [24]

    Le Roy R J 2007 LEVEL 8.0: A Computer Program for Solving the Radial Schrödinger Equation for Bound and Quasi-bound Levels (Waterloo: University of Waterloo) Chemical Physics Research Report CP-663

    [25]

    Fu M K, Ma H T, Cao J W, Bian W S 2016 J. Chem. Phys. 144 184302

    [26]

    Adema Z, Makhlouf S, Taher F 2016 Comput. Theor. Chem. 1093 48

    [27]

    Huber K P, Herzberg G 1979 Constants of Diatonic Molecules, Molecular Spectra Molecular Structure (Vol. IV) (New York: Van Nostrand Reinhold)

    [28]

    Lima J C B D, Ornellas F R 2013 J. Mol. Spectrosc. 283 22

    [29]

    Wu D L, Lin C Q , Wen Y F, Xie A D, Yan B 2018 Chin. Phys. B 27 083101

    [30]

    Wu D L, Tan B, Wen Y F, Zeng X F, Xie A D, Yan B 2016 Spectrochim Acta Part A 161 101

    [31]

    魏长立, 梁桂颖, 刘晓婷, 颜培源, 闫冰 2016 物理学报 65 163101

    Wei C L, Liang G Y, Liu X T, Yan P Y, Yan B 2016 Acta Phys. Sin. 65 163101

    [32]

    Zhang X M, Liang G Y, Li R, Shi D D, Liu Y C, Liu X S, Xu H F, Yan B 2014 Chem. Phys. 443 142

    [33]

    Okabe H 1978 Photochemistry of Small Molecules (New York: Wiley-Interscience)

    [34]

    Zou W L, Liu W J 2005 J. Comput. Chem. 26 106

    [35]

    Paul J D, Howard W C, Richard N Z 1974 J. Chem. Phys. 60 2330

  • 图 1  Sr35Cl分子14个激发态的势能曲线

    Fig. 1.  Potential energy curves of 14 excited states of Sr35Cl.

    图 2  Sr35Cl分子5个束缚态的电偶极矩

    Fig. 2.  Permanent dipole moments of 5 bound states of Sr35Cl.

    图 3  Sr35Cl分子5个束缚态的跃迁偶极矩

    Fig. 3.  Transition dipole moments of 5 bound states of Sr35Cl.

    表 1  5个束缚态的光谱常数

    Table 1.  Spectroscopic constants of the 5 bound states.

    Λ-S态 Te/cm−1 Re/nm $\omega $e/cm−1 $\omega $e$\chi $e/cm−1 Be/cm−1 104αe/cm−1 De/eV Re附近主要电子组态(%)
    $ {{\rm{X}}^{2}}\Sigma^+$ 0.0 0.2575 309.78 0.8682 0.1016 4.131 3.703 6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$24$\text{π}$2(78.8)
    6${\rm{\sigma }}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$24$\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{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ 03$\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{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$24$\text{π}$2(78.8)
    6${\rm{\sigma }}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$24$\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{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ 03$\text{π}$ααβ 4$\text{π}$2(59.7)
    6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ 29$\text{σ}$ 03$\text{π}$24$\text{π}$2(13.4)
    6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ β9$\text{σ}$ 03$\text{π}$ααβ2(9.8)
    6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ α9$\text{σ}$ 03$\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{σ}$ 27$\text{σ}$ α8$\text{σ}$ 29$\text{σ}$ 03$\text{π}$24$\text{π}$2(51.0)
    6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$24$\text{π}$2(16.4)
    6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$ααβ 4$\text{π}$2(4.6)
    6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$β4$\text{π}$ααβ (4.6)
    6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$24$\text{π}$2(2.5)
    6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$αββ 4$\text{π}$2(1.6)
    6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$α4$\text{π}$αββ (1.6)
    6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ 09$\text{σ}$ 03$\text{π}$44$\text{π}$2(1.5)
    6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ 09$\text{σ}$ 03$\text{π}$24$\text{π}$4(1.5)
    理论[26] 27979 0.248 358 1.01 0.1095
    实验[15] 28822 344[27] 1.04[27] 0.1030[16]
    下载: 导出CSV

    表 2  Sr35Cl分子$ {{\rm{X}}^{2}}\Sigma^+$, $ {{\rm{A}}^{2}}\Pi$, $ {{\rm{B}}^{2}}\Sigma^+$, $ {{\rm{C}}^{2}}\Pi$$ {{\rm{3}}^{2}}\Sigma^+$Gv, BvDv

    Table 2.  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 Sr35Cl.

    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
    108Dv/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
    108Dv/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
    108Dv/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
    108Dv/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
    108Dv/cm−1 0.363052 0.690909 1.002658 1.373820 1.690530 1.922724 2.122355 2.300506 2.453955 2.583320
    下载: 导出CSV

    表 3  Sr35Cl分子$ {{\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因子

    Table 3.  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^+ $.

    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
    下载: 导出CSV

    表 4  Sr35Cl分子$ {{\rm{A}}^{2}}\Pi $$ {{\rm{X}}^{2}}\Sigma^+ $, $ {{\rm{B}}^{2}}\Sigma^+ $$ {{\rm{X}}^{2}}\Sigma^+ $$ {{\rm{C}}^{2}}\Pi $$ {{\rm{X}}^{2}}\Sigma^+ $跃迁的辐射寿命

    Table 4.  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^+ $.

    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
    下载: 导出CSV
  • [1]

    Micheli A, Brennen G, Zoller P 2006 Nat. Phys. 2 341

    [2]

    Baron J, Campbell W C, DeMille D, Doyle J M, Gabrielse G, Gurevich Y V, Hess P W, Hutzler N R, Kirilov E, Kozyryev I, O'Leary B R, Panda C D, Parsons M F, Petrik E S, Spaun B, Vutha A C, West A D 2014 Science 343 269

    [3]

    Krems R V 2008 Phys. Chem. Chem. Phys. 10 4079

    [4]

    Walters O H, Barratt S 1928 Proc. R. Soc. London Ser. A 118 120

    [5]

    Nakagawa J, Domaille P J, Steimle T C, Harris D O 1978 J. Mol. Spectrosc. 70 374

    [6]

    Dulick M, Bernath P F, Field R W 1980 Can. J. Phys. 58 703

    [7]

    Domaille P J, Steimle T C, Wong N B, Harris D O 1977 J. Mol. Spectrosc. 65 354

    [8]

    Berg L E, Klynning L, Martin H 1980 Phys. Scr. 21 173

    [9]

    Bernath P F, Field R W, Pinchemel B, Lefebvre Y, Schamps J 1981 J. Mol. Spectrosc. 88 175

    [10]

    Reisner D E, Bernath P F, Field R W 1981 J. Mol. Spectrosc. 89 107

    [11]

    Zare R N, Schmeltekopf A L, Harrop W J, Albritton D L 1973 J. Mol. Spectrosc. 46 37

    [12]

    Singh J, Nair K P R, Upadhya K N, Rai D K 1970 Opt. Pure Appl. 3 76

    [13]

    Brinkmann U, Schmidt V H, Telle H 1982 Chern. Phys. 64 19

    [14]

    Schiitze Pahlmann H V, Ryzlewicz C H, Hoeft J, Torring T 1982 Chern. Phys. Lett. 93 74

    [15]

    Ernst W E, Schröder J O 1984 J. Chem. Phys. 81 136

    [16]

    Schruder J O, Zeller B, Ernst W E 1988 J. Mol. Spectrosc. 127 255

    [17]

    Berg L E, Royen P, Weijnitz P 1990 Mol. Phys. 69 385

    [18]

    Werner H J, Knowles P J, Knizia G, et al. 2012 MOLPRO, Version 2012.1, A Package of ab initio Programs

    [19]

    Peterson K A, Figgen D, Goll E, Stoll H, Dolg M 2003 J. Chem. Phys. 119 11099

    [20]

    Werner H J, Knowles P J 1985 J. Chem. Phys. 82 5053

    [21]

    Knowles P J, Werner H J 1985 Chem. Phys. Lett. 115 259

    [22]

    Werner H J, Knowles P J 1988 J. Chem. Phys. 89 5803

    [23]

    Knowles P J, Werner H J 1988 Chem Phys Lett. 145 514

    [24]

    Le Roy R J 2007 LEVEL 8.0: A Computer Program for Solving the Radial Schrödinger Equation for Bound and Quasi-bound Levels (Waterloo: University of Waterloo) Chemical Physics Research Report CP-663

    [25]

    Fu M K, Ma H T, Cao J W, Bian W S 2016 J. Chem. Phys. 144 184302

    [26]

    Adema Z, Makhlouf S, Taher F 2016 Comput. Theor. Chem. 1093 48

    [27]

    Huber K P, Herzberg G 1979 Constants of Diatonic Molecules, Molecular Spectra Molecular Structure (Vol. IV) (New York: Van Nostrand Reinhold)

    [28]

    Lima J C B D, Ornellas F R 2013 J. Mol. Spectrosc. 283 22

    [29]

    Wu D L, Lin C Q , Wen Y F, Xie A D, Yan B 2018 Chin. Phys. B 27 083101

    [30]

    Wu D L, Tan B, Wen Y F, Zeng X F, Xie A D, Yan B 2016 Spectrochim Acta Part A 161 101

    [31]

    魏长立, 梁桂颖, 刘晓婷, 颜培源, 闫冰 2016 物理学报 65 163101

    Wei C L, Liang G Y, Liu X T, Yan P Y, Yan B 2016 Acta Phys. Sin. 65 163101

    [32]

    Zhang X M, Liang G Y, Li R, Shi D D, Liu Y C, Liu X S, Xu H F, Yan B 2014 Chem. Phys. 443 142

    [33]

    Okabe H 1978 Photochemistry of Small Molecules (New York: Wiley-Interscience)

    [34]

    Zou W L, Liu W J 2005 J. Comput. Chem. 26 106

    [35]

    Paul J D, Howard W C, Richard N Z 1974 J. Chem. Phys. 60 2330

  • [1] 陈恒杰. LiAl分子基态、激发态势能曲线和振动能级. 物理学报, 2013, 62(8): 083301. doi: 10.7498/aps.62.083301
    [2] 孙金锋, 马 恒, 朱遵略, 杨向东, 施德恒. 7Li2(X1Σ+g)分子的振动能级、转动惯量及离心畸变常数. 物理学报, 2008, 57(1): 165-171. doi: 10.7498/aps.57.165
    [3] 孙金锋, 马 恒, 朱遵略, 施德恒. 7Li2分子2 3Σ+g激发态的解析势能函数、谐振频率及振动能级. 物理学报, 2007, 56(4): 2085-2091. doi: 10.7498/aps.56.2085
    [4] 刘玉芳, 马 恒, 孙金锋, 朱遵略, 杨向东, 施德恒. 7Li2分子23Πu激发态的解析势能函数、振动能级及其转动惯量. 物理学报, 2007, 56(8): 4454-4460. doi: 10.7498/aps.56.4454
    [5] 孙卫国, 任维义. Na2分子部分电子态的完全振动能谱和离解能的精确研究. 物理学报, 2005, 54(2): 594-605. doi: 10.7498/aps.54.594
    [6] 冯 灏, 胡士德, 孙卫国, 任维义. 碱金属氢化物双原子分子部分电子态的完全振动能谱和分子离解能的精确研究. 物理学报, 2006, 55(5): 2185-2193. doi: 10.7498/aps.55.2185
    [7] 张计才, 孙金锋, 施德恒, 朱遵略. BeC分子基态和低激发态光谱性质和解析势能函数. 物理学报, 2019, 68(5): 053102. doi: 10.7498/aps.68.20181695
    [8] 李松, 韩立波, 陈善俊, 段传喜. SN-分子离子的势能函数和光谱常数研究. 物理学报, 2013, 62(11): 113102. doi: 10.7498/aps.62.113102
    [9] 黄多辉, 王藩侯, 杨俊升, 万明杰, 曹启龙, 杨明超. SnO分子的X1Σ+, a3Π和A1Π态的势能曲线与光谱性质. 物理学报, 2014, 63(8): 083102. doi: 10.7498/aps.63.083102
    [10] 孙卫国, 樊群超, 冯灏. 双原子离子振动光谱与离解能的研究. 物理学报, 2010, 59(1): 203-209. doi: 10.7498/aps.59.203
    [11] 蹇君, 雷娇, 樊群超, 范志祥, 马杰, 付佳, 李会东, 徐勇根. NO分子宏观气体热力学性质的理论研究. 物理学报, 2020, 69(5): 053301. doi: 10.7498/aps.69.20191723
    [12] 程新路, 苏欣纺, 陈恒杰, 唐海燕, 王全武. LiC分子基态及其低电子激发态的多参考组态相互作用方法研究. 物理学报, 2010, 59(7): 4556-4563. doi: 10.7498/aps.59.4556
    [13] 李会东, 樊群超, 冯灏, 孙卫国. Li2部分电子态的完全振动能谱与离解能的精确研究. 物理学报, 2010, 59(7): 4577-4583. doi: 10.7498/aps.59.4577
    [14] 高 峰, 杨传路, 张晓燕. 多参考组态相互作用方法研究ZnHg低激发态(1∏,3∏)的势能曲线和解析势能函数. 物理学报, 2007, 56(5): 2547-2552. doi: 10.7498/aps.56.2547
    [15] 邢伟, 刘慧, 施德恒, 孙金锋, 朱遵略. MRCI+Q理论研究SiSe分子X1Σ+和A1Π电子态的光谱常数和分子常数. 物理学报, 2013, 62(4): 043101. doi: 10.7498/aps.62.043101
    [16] 刘 艳, 任维义, 王阿署, 刘松红. 分子部分电子态的高阶振动能级和离解能的精确研究. 物理学报, 2008, 57(3): 1599-1607. doi: 10.7498/aps.57.1599
    [17] 刘慧, 施德恒, 孙金锋, 朱遵略. MRCI方法研究CSe(X1Σ+)自由基的光谱常数和分子常数. 物理学报, 2011, 60(6): 063101. doi: 10.7498/aps.60.063101
    [18] 向天翔, 孙胜, 龚顺生, 王嘉铭. 态—态振动能量跃迁的时间分辨研究(Ⅰ)——碘分子的自碰撞过程. 物理学报, 1990, 39(10): 1547-1554. doi: 10.7498/aps.39.1547
    [19] 刘慧, 邢伟, 施德恒, 朱遵略, 孙金锋. 用MRCI方法研究CS+同位素离子X2Σ+和A2Π态的光谱常数与分子常数. 物理学报, 2011, 60(4): 043102. doi: 10.7498/aps.60.043102
    [20] 田寅, 冯灏, 孙卫国. 碱金属双原子分子部分电子态的完全振动能谱和离解能. 物理学报, 2011, 60(2): 023301. doi: 10.7498/aps.60.023301
  • 引用本文:
    Citation:
计量
  • 文章访问数:  599
  • PDF下载量:  8
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-26
  • 修回日期:  2018-12-11
  • 上网日期:  2019-02-01
  • 刊出日期:  2019-02-05

单氯化锶分子低激发态的光谱及跃迁特性

    基金项目: 国家自然科学基金(批准号: 11564019, 11147158)和江西省教育厅科学技术研究项目(批准号: GJJ170654)资助的课题.

摘要: 利用Davidson修正的内收缩多参考组态相互作用(ic-MRCI + Q)方法, 结合相对论有效芯赝势基(aug-cc-pV5Z-PP)作为Sr原子和相关一致五重基aug-cc-pV5Z为Cl原子的计算基组, 优化计算了单氯化锶(Sr35Cl)分子14个低激发电子态的势能曲线和跃迁偶极矩. 为了获得更加精确的光谱参数, 计算中同时引入核价电子相关和相对论效应修正势能曲线. 利用LEVEL 8.0程序拟合修正的势能曲线, 得到相应电子态的光谱常数、振动能级和分子常数等光谱性质, 结果与近来的已获得的理论计算和实验值符合得较好, 同时给出了Franck-Condon因子和辐射寿命等跃迁性质. 这些精确的光谱跃迁特性可为进一步构建Sr35Cl分子激光冷却方案提供理论支持.

English Abstract

    • 碱土金属卤化物在天体物理学、催化和高温化学中的应用广泛, 再加上特殊的电子组态引起了许多实验与理论光谱学家的兴趣, 一直是许多电子结构和光谱性质研究的主题[13]. 1928年, Walters和Barratt[4]最早报道了单卤化锶(Sr35Cl)的光学光谱; 后来由于该类化合物基态(${\rm{X}^2}\Sigma^+ $)和第一、二激发态${\rm{A}^2}\Pi $, ${\rm{B}^2} \Sigma^+$)的势能曲线和核间距非常相似, 导致光谱密集不便于区分. 随着激光器的出现, 许多研究者[510]精确地分析了${\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), 优化计算了Sr35Cl分子的较低的两个解离极限14个激发态在全空间的电子结构和能量. 为了提高计算精度, 计算中考虑了Davison修正、核价电子关联效应和相对论效应修正, 最终得到14个激发态的势能曲线. 基于修正的5个较低束缚态的势能曲线和跃迁偶极矩, 利用LEVEL 8.0程序拟合计算得到相应束缚态的光谱常数、振动能级和分子常数, 同时给出了Franck-Condon因子和辐射寿命等跃迁性质. 这些精确的光谱跃迁特性可为进一步构建SrCl分子激光冷却方案提供理论支持.

    • 本文所有与电子结构和能量有关的计算均使用MOLPRO2012.1程序包完成[18]. 经过测试不同的基组, 最终选取相对论有效芯赝势基(aug-cc-pV5Z-PP)作为Sr原子和相关一致五重基(aug-cc-pV5Z)作为Cl原子的计算基组[19].

      Sr35Cl分子属于Cv点群. 由于MOLPRO程序自身的限制, 在实际计算中使用Cv点群的阿贝尔子群C2v, C2v群的四个不可约表示(A1, B1, B2, A2)和Cv点群的不可约表示有如下对应关系: $\Sigma^+ $ = A1, $\prod $ = B1 + B2, $\Delta $ = A1 + A2, $ \Sigma^- $ = A2. 在计算过程中Sr35Cl分子的8个轨道(4个a1, 2个b1和2个b2)被确定为活性空间, 它们对应Sr原子的4s4p轨道和Cl原子的3s3p轨道, Sr原子的外层电子5s2和Cl原子的外层电子3s23p5被放入活性空间进行相关能的计算, 余下的48个电子被冻结在芯轨道. 在C2v对称性下, 扫描计算Sr35Cl分子在核间距R = 0.16—0.80 nm之间的单点能, 步长设为0.01 nm, 其中在平衡位置附近选取步长为0.005 nm. 扫描过程中, 首先使用Hartree-Fock方法计算Sr35Cl分子基态的单组态波函数, 然后利用态平均全活性空间自洽场方法(SA-CASSCF)优化单组态波函数为多组态波函数[20,21], 最后应用内收缩多参考态相互作用方法(ic-MRCI)计算Sr35Cl分子束缚态的能量[22,23]. 为了进一步提高精度, 计算中考虑Davidson修正(+Q), 核价电子关联效应和相对论标量效应修正, 最后得到了Sr35Cl分子最低的14个激发态的势能曲线, 其中5个束缚态, 9个排斥态.

      基于修正的5个束缚态的势能曲线和跃迁偶极矩, 利用LEVEL 8.0程序[24]求解径向Schrödinger方程获得了全部束缚态的光谱常数和分子常数, 包括离解能De、平衡核间距Re、绝热激发能Te、谐性振动常数$\omega_{\rm e}$和非谐性振动常数$\omega_{\rm e}$$\chi_{\rm e}$, 转动常数Be和振转耦合常数$\alpha_{\rm e}$, 以及振动能级Gv、惯性转动常数Bv和离心畸变常数Dv. 同时给出了Franck-Condon因子和辐射寿命等跃迁性质.

    • 利用内收缩的多参考组态相互作用方法(ic-MRCI + Q)结合相关一致五重基aug-cc-pV5Z和相对论有效芯赝势基aug-cc-pV5Z-PP优化计算Sr35Cl分子的电子结构. 分析对应解离极限Sr(1Sg) + Cl(2Pu), Sr(3Pu) + Cl(2Pu)和离子对Sr+(2Sg) + Cl(1Sg)的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, 同时还列出了各个电子态在平衡位置Re附近主要的电子组态分布及权重.

      图  1  Sr35Cl分子14个激发态的势能曲线

      Figure 1.  Potential energy curves of 14 excited states of Sr35Cl.

      Λ-S态 Te/cm−1 Re/nm $\omega $e/cm−1 $\omega $e$\chi $e/cm−1 Be/cm−1 104αe/cm−1 De/eV Re附近主要电子组态(%)
      $ {{\rm{X}}^{2}}\Sigma^+$ 0.0 0.2575 309.78 0.8682 0.1016 4.131 3.703 6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$24$\text{π}$2(78.8)
      6${\rm{\sigma }}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$24$\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{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ 03$\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{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$24$\text{π}$2(78.8)
      6${\rm{\sigma }}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$24$\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{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ 03$\text{π}$ααβ 4$\text{π}$2(59.7)
      6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ 29$\text{σ}$ 03$\text{π}$24$\text{π}$2(13.4)
      6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ β9$\text{σ}$ 03$\text{π}$ααβ2(9.8)
      6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ α9$\text{σ}$ 03$\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{σ}$ 27$\text{σ}$ α8$\text{σ}$ 29$\text{σ}$ 03$\text{π}$24$\text{π}$2(51.0)
      6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$24$\text{π}$2(16.4)
      6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$ααβ 4$\text{π}$2(4.6)
      6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$β4$\text{π}$ααβ (4.6)
      6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$24$\text{π}$2(2.5)
      6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ α9$\text{σ}$ 03$\text{π}$αββ 4$\text{π}$2(1.6)
      6$\text{σ}$ 27$\text{σ}$ 28$\text{σ}$ 09$\text{σ}$ α3$\text{π}$α4$\text{π}$αββ (1.6)
      6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ 09$\text{σ}$ 03$\text{π}$44$\text{π}$2(1.5)
      6$\text{σ}$ 27$\text{σ}$ α8$\text{σ}$ 09$\text{σ}$ 03$\text{π}$24$\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个束缚态的光谱常数

      Table 1.  Spectroscopic constants of the 5 bound states.

      图1表1可看出, 基态${\rm X}^2\Sigma^+ $为一个深势阱态, 阱深为3.703 eV, 说明基态是一个比较稳定的电子态. 在平衡位置附近主要的电子组态为6${\text{σ}}^2$7${\text{σ}}$28${\text{σ}}$α9${\text{σ}}$03$\text{π}$24$ \text{π}$2和6${\text{σ}}$27${\text{σ}}$28${\text{σ}}$09${\text{σ}}$α3$\text{π}$24$ \text{π}$2, 表明了使用多组态方法的必要性[2830]. 通过比较发现, 本文的光谱常数与近来获得的理论计算值[26]较接近, Re, $\omega_{\rm e}$Be的偏差分别为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]. 这些对比说明Sr35Cl分子相比BeCl, MgCl和CaCl分子体系, 稳定性更差, 更容易离解.

      对于第一激发态${\rm A}^2 \Pi$, 与基态属于同一个解离极限Sr(1Sg) + Cl(2Pu), 但离解能为1.673 eV, 结果小于基态3.703 eV, 这说明基态更加稳定. 平衡位置附近主要的电子组态为6$\text{σ}$27$\text{σ}$28$\text{σ}$09$\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} $, ReBe与实验结果[16]有一定相差, 但与近来的理论计算[26]符合得很好, 偏差分别为7.69 cm−1, 0.0002 nm和0.0006 cm−1.

      特别值得注意的是第二激发态B2$ \Sigma^+$, 势能曲线出现了双势阱, 且该激发态对应离子对解离极限Sr+(2Sg) + Cl(1Sg), 绝热激发能为16612.74 cm−1, 与实验值15719.5 cm−1[12]也有偏差, 相对偏差为5.7%, 这可能也是由于${\rm{B}^2}\Sigma^+$${\rm{3}^2}\Sigma^+ $出现了避免交叉现象, 导致势能曲线发生变化, 致使拟合的光谱常数$\omega $e, ReBe与实验结果[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(3Pu) + Cl(2Pu). 从图1还可以发现${\rm A}^2\Pi$$ {\rm B}^2\Sigma^+$的势能曲线在0.261 nm和0.340 nm位置存在两处交叉, 本文给出了${\rm C^2}\Pi$$3^2\Sigma^+ $激发态的光谱常数, 结果与近来获得的理论值相差较大[26], 这可能是由于电子态之间相互作用的影响所致. 综上分析发现, 本文的光谱常数与实验值和近来的理论计算值有一定差异, 可能是由于本文考虑了避免交叉现象, 即考虑了态与态之间的电子关联相互作用, 这会使拟合得到的光谱常数出现一定偏差, 但当进一步考虑自旋-轨道耦合效应后, 这些势能曲线会发生变化, 从而得到更加精确的光谱常数. 限于篇幅, 后续我们将报道详细的自旋-轨道耦合效应对光谱常数和跃迁特性的影响.

    • 利用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
      108Dv/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
      108Dv/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
      108Dv/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
      108Dv/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
      108Dv/cm−1 0.363052 0.690909 1.002658 1.373820 1.690530 1.922724 2.122355 2.300506 2.453955 2.583320

      表 2  Sr35Cl分子$ {{\rm{X}}^{2}}\Sigma^+$, $ {{\rm{A}}^{2}}\Pi$, $ {{\rm{B}}^{2}}\Sigma^+$, $ {{\rm{C}}^{2}}\Pi$$ {{\rm{3}}^{2}}\Sigma^+$Gv, BvDv

      Table 2.  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 Sr35Cl.

    • 由于偶极矩反映了分子的成键性质, 本文对Sr35Cl分子的5个较低束缚态的电偶极矩进行了理论计算, 电偶极矩随着分子核间距R的变化曲线见图2.

      图  2  Sr35Cl分子5个束缚态的电偶极矩

      Figure 2.  Permanent dipole moments of 5 bound states of Sr35Cl.

      图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变化趋势相同.

    • 图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, 且对角化比较明显. 这正好符合激光冷却分子体系的条件, 因此Sr35Cl分子可以作为激光冷却的候选分子体系. 其次较大的为${\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  Sr35Cl分子$ {{\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因子

      Table 3.  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  Sr35Cl分子5个束缚态的跃迁偶极矩

      Figure 3.  Transition dipole moments of 5 bound states of Sr35Cl.

      基于计算得到的激发态-基态的跃迁偶极矩、电子态的势能曲线和振动能级, 本文计算了最低五个束缚激发态振动能级的辐射寿命. 激发态振动能级v′的辐射寿命计算公式如下[33, 34]:

      $\begin{split} {\tau _{v^{\prime}}} & = {({A_{v^{\prime}}})^{ - 1}} \\ &= \frac{{3h}}{{64{{\text{π}}^4}{{\left| {{a_0} \cdot e \cdot \overline {TDM} } \right|}^2}\sum\nolimits_{v^{\prime \prime}} {{q_{v^{\prime},v^{\prime \prime}}}{{(\Delta {E_{v^{\prime}, v^{\prime \prime}}})}^3}} }}\\ & = \frac{{4.936 \times {{10}^5}}}{{{{(\overline {TDM} )}^2}\sum\nolimits_{v^{\prime \prime}} {{q_{v^{\prime}, v^{\prime \prime}}}{{(\Delta {E_{v^{\prime},v^{\prime \prime}}})}^3}} }}, \end{split}$

      式中${\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  Sr35Cl分子$ {{\rm{A}}^{2}}\Pi $$ {{\rm{X}}^{2}}\Sigma^+ $, $ {{\rm{B}}^{2}}\Sigma^+ $$ {{\rm{X}}^{2}}\Sigma^+ $$ {{\rm{C}}^{2}}\Pi $$ {{\rm{X}}^{2}}\Sigma^+ $跃迁的辐射寿命

      Table 4.  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^+ $.

    • 利用高精度的内收缩多参考组态相互作用方法(ic-MRCI), 考虑多种相关修正研究了Sr35Cl分子对应最低三个解离极限的激发态电子结构、势能曲线和偶极矩等. 结果表明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因子和辐射寿命等跃迁性质, 获得了比较完善的激发态光谱性质. 这些精确的光谱跃迁特性可为进一步构建Sr35Cl分子激光冷却方案提供理论支持.

参考文献 (35)

目录

    /

    返回文章
    返回