高级搜索

LiCl–阴离子的光谱性质和跃迁性质

,  ,  ,  , 

downloadPDF

上一篇

下一篇

郭芮, 谭涵, 袁沁玥, 张庆, 万明杰. LiCl–阴离子的光谱性质和跃迁性质[J]. 物理学报, 2022, 71(4): 043101-1. doi: 10.7498/aps.71.20211688
引用本文: 郭芮, 谭涵, 袁沁玥, 张庆, 万明杰. LiCl–阴离子的光谱性质和跃迁性质[J]. 物理学报, 2022, 71(4): 043101-1. doi: 10.7498/aps.71.20211688
shu
Rui Guo, Han Tan, Qin-Yue Yuan, Qing Zhang, Ming-Jie Wan. Spectroscopic and transition properties of LiCl– anion[J]. Acta Physica Sinica, 2022, 71(4): 043101-1. doi: 10.7498/aps.71.20211688
Citation: Rui Guo, Han Tan, Qin-Yue Yuan, Qing Zhang, Ming-Jie Wan. Spectroscopic and transition properties of LiCl– anion[J]. Acta Physica Sinica, 2022, 71(4): 043101-1. doi: 10.7498/aps.71.20211688
shu

LiCl–阴离子的光谱性质和跃迁性质

  • 宜宾学院理学部, 宜宾　644007

    • 通讯作者: E-mail: wanmingjie1983@sina.com

  • 中图分类号: 31.15.A-, 31.15.aj, 37.10.Mn, 87.80.Cc

Spectroscopic and transition properties of LiCl– anion

  • Faculty of Science, Yibin University, Yibin 644007, China

    • Corresponding author: E-mail: wanmingjie1983@sina.com

  • MSC: 31.15.A-, 31.15.aj, 37.10.Mn, 87.80.Cc

  • 摘要: 采用多参考组态相互作用方法结合全电子基组计算了LiCl- 阴离子5个电子态 (X2Σ+, A2∏, B2Σ+, 32Σ+, 22∏) 的电子结构. 为了得到精确的光谱常数, 计算中考虑了Davidson 修正、芯-价电子关联效应和自旋-轨道耦合效应. 拟合得到各电子态的光谱常数、分子常数、自发辐射速率和自发辐射寿命. 基态的光谱常数与实验值和其他理论值符合较好, 同时报道了LiCl– 阴离子激发态的光谱常数以及其到基态的跃迁性质. 计算结果表明A2∏ $\leftrightarrow $X2Σ+跃迁具有高对角分布的弗兰克-康登因子f00, 第一激发态A2∏有较短的自发辐射寿命. 构造A2∏ (ν′) $\leftrightarrow $ X2Σ+ ($\nu''$)准循环跃迁进行激光冷却LiCl– 阴离子需要一束主激光和两束抽运激光. 以上结果预测了激光冷却LiCl–阴离子是可行的.
    Abstract: The electronic structure of the X2Σ+, A2Π, B2Σ+, 32Σ+, and 22Π state of LiCl– anion are performed at an MRCI+Q level. Davison correction, core-valence correction and spin-orbit coupling effect are also considered. The ground state X2Σ+ of LiCl– anion correlates with the lowest dissociation channel Li(2Sg) + Cl–(1Sg); the A2∏ state and B2Σ+ state correlate with the second dissociation channel Li(2Pu) + Cl–(1Sg); the 32Σ+ state and 22Π state correlate with the third dissociation channel Li–(1Sg) + Cl–(2Pu).Spectroscopic parameters are calculated by solving the radial Schröedinger equation. The equilibrium internuclear distance Re of the ground state X2Σ+ is 2.1352 Å, which is a little bigger than the experimental datum, with an error being 0.5%. It is a deep potential well, and the dissociation energy De is 1.886 eV. These values are in good agreement with experimental data. The A2∏ state is at 13431.93 cm–1 above the X2Σ+ state. The Re is 2.1198 Å, which is only 0.0154 Å smaller than that of the X2Σ+ state. The values of energy level Gν and rotational constant Bν of five Λ-S states are also calculated. The values are in good agreement with available theoretical ones. The electronic structures of the excited states are also reported. The SOC effect weakly influences the spectroscopic parameters for the $ {\text{X}}{}^2\Sigma _{1/2}^ + $, $ {\text{A}}{}^2{\Pi _{1/2}} $, $ {\text{A}}{}^2{\Pi _{3/2}} $, and $ {\text{B}}{}^2\Sigma _{1/2}^ + $ state. From the analysis of the SO matrix, it can be seen that the SOC effect plays a little role in realizing the A2Π $\leftrightarrow $ X2Σ+ transition, so, it can be ignored.The scheme of laser cooling of LiCl– anion has constructed at a spin – free level. The A2∏(ν′) $\leftrightarrow $ X2Σ+($v'' $) transition has a highly diagonally distributed Franck-Condon factor f00 = 0.9898, the calculated branching ratio of the diagonal term R00 is 0.9893, and spontaneous radiative lifetime of A2∏ is 35.45 ns. A main pump laser and two repumping lasers for driving the A2∏(ν′) $\leftrightarrow $ X2Σ+($v'' $) transitions are required. The laser wavelengths are 744.10, 774.30 and 772.42 nm, respectively. Owing to the summation of R00, R01, and R02 being closer to 1, the A2∏(ν′) $\leftrightarrow $ X2Σ+($v'' $) transition is a quasicycling transition. These results imply that the LiCl– anion is a candidate for laser cooling.
  • 加载中

  • 图 1  X2Σ+, A2∏, B2Σ+, 32Σ+和22∏电子态的势能曲线

    Figure 1.  Potential energy curves of the X2Σ+, A2∏, B2Σ+, 32Σ+ and 22∏ states.

    下载: 全尺寸图片 幻灯片

    图 2  Ω态的势能曲线

    Figure 2.  Potential energy curves of the Ω states.

    下载: 全尺寸图片 幻灯片

    图 3  LiCl–阴离子的自旋-轨道矩阵元. SOi的表示见表5

    Figure 3.  Spin-orbit matrix elements of the of the LiCl– anion. The explanations of the SOi symbols are presented in Table 6.

    下载: 全尺寸图片 幻灯片

    图 4  (a) Λ-S态的电偶极矩; (b) Λ-S态的跃迁偶极矩

    Figure 4.  (a) The permanent dipole moments of the Λ-states; (b) the transition dipole moments of the Λ-states.

    下载: 全尺寸图片 幻灯片

    图 5  驱动A2∏ $\leftrightarrow $ X2Σ+跃迁进行激光冷却的途径

    Figure 5.  Proposed laser cooling scheme via A2∏ $\leftrightarrow $ X2Σ+ transition.

    下载: 全尺寸图片 幻灯片

    表 1  LiCl–阴离子Λ-S态的离解关系

    Table 1.  Calculated dissociation relationships of the Λ-S states of LiCl– anion.

    ΔE/cm–1
    原子态分子态本文工作实验值[36–38]
    Li(2Sg)+Cl–(1Sg)X2Σ+00
    Li(2Pu)+ Cl–(1Sg)A2Π, B2Σ+14903.7914253.13
    Li–(1Sg)+Cl(2Pu)32Σ+, 22Π23703.6124594.67
    下载: 导出CSV

    表 2  LiCl–阴离子Λ-S态的光谱常数

    Table 2.  Spectroscopic parameters of the Λ-S states of LiCl– anion.

    电子态Re/Åωe/cm–1ωeχe/cm–1Be/cm–1De/eVD0/eVTe/cm–1
    X2Σ+2.1352535.335.81730.72051.88861.8560
    实验值[20]2.18(4)0
    实验值[21]2.123(15)1.75(2)0
    理论值[22]2.12a0
    理论值[23]2.1354537.7b6.34b0.7203b1.810
    A2∏2.1198554.655.71570.73101.99021.955913431.93
    B2Σ+2.0282652.796.12520.79851.66531.625017491.75
    32Σ+5.859430.190.84830.09630.03620.035338607.64
    22∏7.141139.080.56160.06380.01360.011238855.32
    a采用HF方法计算得到基态的核间距.
    下载: 导出CSV

    表 3  X2Σ+, A2∏, B2Σ+, 32Σ+和22∏态的振动能级和转动常数

    Table 3.  Vibrational energy levels, rotational constants of the X2Σ+, A2∏, B2Σ+, 32Σ+ and 22∏ states.

    νX2Σ+ A2∏ B2Σ+ 32Σ+ 22∏
    Gν BνGνBνGνBνGνBνGνBν
    本文工作文献[23]本文工作文献[23]本文工作本文工作本文工作本文工作本文工作本文工作本文工作本文工作
    0266.72264.07 0.71460.7143 276.480.7252 325.330.7927 14.960.0940 6.120.0631
    1790.93791.770.70290.7023820.210.7136966.070.781143.020.089318.150.0618
    21303.271307.010.69110.69031352.130.70211594.230.769768.840.085728.920.0572
    31803.711809.920.67940.67841872.340.69072210.020.758493.220.082629.300.0095
    42292.312300.650.66770.66662380.970.67932813.560.7471116.240.079631.320.0112
    52769.152779.340.65600.65482878.130.66803404.890.7359138.170.076533.190.0125
    63234.293246.110.64440.64303363.870.65673984.090.7247158.980.072535.070.0140
    73687.833701.100.63280.63133838.290.64554551.200.7135178.260.067836.530.0310
    84129.884144.450.62120.61974301.420.63435106.170.7023195.600.062637.330.0257
    94560.614576.290.60970.60814753.370.62315648.950.6911210.770.057139.110.0188
    下载: 导出CSV

    表 4  LiCl–阴离子Ω态的离解关系

    Table 4.  Calculated dissociation relationships of the Ω states of LiCl– anion.

    ΔE/cm–1
    原子态分子态 Ω本文工作实验值[36–38]
    Li(2S1/2)+Cl–(1S0)1/200
    Li(2P1/2)+ Cl–(1S0)1/214252.7714903.62
    Li(2P3/2)+ Cl–(1S0)1/2, 3/214253.5014903.96
    Li–(1S0)+Cl(2P3/2)1/2, 3/223415.4124153.49
    Li–(1S0)+Cl(2P1/2)1/224288.4825035.84
    下载: 导出CSV

    表 5  LiCl–阴离子Ω态的光谱常数

    Table 5.  Spectroscopic parameters of the Ω states of LiCl– anion at icMRCI+Q level.

    Ω态Re/Åωe/cm–1ωeχe /cm–1 Be/cm–1De/eVTe/cm–1
    $ {\text{X}}{}^2\Sigma _{1/2}^ + $2.1352535.325.81720.72051.88860
    $ {\text{A}}{}^2{\Pi _{1/2}} $2.1196554.875.71530.73112.009413419.77
    $ {\text{A}}{}^2{\Pi _{3/2}} $2.1200554.415.71600.73082.005913444.07
    $ {\text{B}}{}^2\Sigma _{1/2}^ + $2.0282652.796.12530.79851.667917491.75
    下载: 导出CSV

    表 6  LiCl–阴离子的自旋-轨道矩阵元

    Table 6.  Notation of spin-orbit matrix element.

    SO矩阵元
    $ {\text{S}}{{\text{O}}_1} = - {\text{i}}\left\langle {{{\text{X}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_2} = \left\langle {{{\text{X}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _x}} \right\rangle $$ {\text{S}}{{\text{O}}_3} = - {\text{i}}\left\langle {{{\text{B}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_4} = \left\langle {{{\text{B}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _x}} \right\rangle $
    $ {\text{S}}{{\text{O}}_5} = - {\text{i}}\left\langle {{3^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_6} = \left\langle {{3^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _x}} \right\rangle $$ {\text{S}}{{\text{O}}_7} = - {\text{i}}\left\langle {{{\text{X}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{2^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_8} = \left\langle {{{\text{X}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{2^2}{\Pi _x}} \right\rangle $
    $ {\text{S}}{{\text{O}}_9} = - {\text{i}}\left\langle {{{\text{B}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{2^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_{10}} = \left\langle {{{\text{B}}^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{2^2}{\Pi _x}} \right\rangle $$ {\text{S}}{{\text{O}}_{11}} = - {\text{i}}\left\langle {{3^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{2^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_{12}} = \left\langle {{3^2}{\Sigma ^ + }} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{2^2}{\Pi _x}} \right\rangle $
    $ {\text{S}}{{\text{O}}_{13}} = {\text{i}}\left\langle {{{\text{A}}^2}{\Pi _x}} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_{14}} = {\text{i}}\left\langle {{2^2}{\Pi _x}} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _y}} \right\rangle $$ {\text{S}}{{\text{O}}_{15}} = {\text{i}}\left\langle {{2^2}{\Pi _y}} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{{\text{A}}^2}{\Pi _x}} \right\rangle $$ {\text{S}}{{\text{O}}_{16}} = {\text{i}}\left\langle {{2^2}{\Pi _x}} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {{2^2}{\Pi _y}} \right\rangle $
    下载: 导出CSV

    表 7  A2∏ $\leftrightarrow $ X2Σ+和B2Σ+ $\leftrightarrow $ X2Σ+跃迁的弗兰克-康登因子fν'ν'', 爱因斯坦系数Aν'ν''和自发辐射寿命 (单位: ns)

    Table 7.  Franck-Condon Factors fν'ν'', Einstein coefficients Aν'ν'', and radiative lifetimes τ of the A2∏ $\leftrightarrow $ X2Σ+ and B2Σ+ $\leftrightarrow $ X2Σ+ transitions of LiCl– anion(in ns).

    跃迁f00f01f02f03
    A00A01A02A03τ = 1/ΣA
    f10f11f12f13
    A10A11A12A13
    A2∏ $\leftrightarrow $ X2Σ+0.98980.01010.00018.70(–7)
    279042002983133848.8544.6035.45
    0.01020.96860.02090.0004
    269660273366006073351233235.43
    B2Σ+ $\leftrightarrow $ X2Σ+0.59080.29090.08940.0225
    183168007658290203873045221034.99
    0.32660.12860.28880.1671
    1227960042617807988480397939033.00
    下载: 导出CSV
  • [1] Micheli A, Brennen G, Zoller P 2006 Nat. Phys. 2 341 doi: 10.1038/nphys287
    [2] Baron J et al., (The ACME Collaboration). 2014 Science 343 269 doi: 10.1126/science.1248213
    [3] Krems R V 2008 Phys. Chem. Chem. Phys. 10 4079 doi: 10.1039/b802322k
    [4] Shuman E S, Barry J F, DeMille D 2010 Nature 467 820 doi: 10.1038/nature09443
    [5] Hummon M T, Yeo M, Stuhl B K, Collopy A L, Xia Y, Ye J 2013 Phys. Rev. Lett. 110 143001 doi: 10.1103/PhysRevLett.110.143001
    [6] You Y, Yang C L, Wang M S, Ma X G, Liu W W 2015 Phys. Rev. A 92 032502 doi: 10.1103/PhysRevA.92.032502
    [7] Wan M J, Shao J X, Huang D H, Jin C G, Yu Y, Wang F H 2015 Phys. Chem. Chem. Phys. 17 26731 doi: 10.1039/C5CP03677A
    [8] Wan M J, Shao J X, Gao Y F, Huang D H, Yang J S, Cao Q L, Jin C G, Wang F H 2015 J Chem. Phys. 143 024302 doi: 10.1063/1.4926389
    [9] Yang Q S, Li S C, Yu Y, Gao T 2018 J. Phys. Chem. A 122 3021 doi: 10.1021/acs.jpca.7b11047
    [10] Fu M K, Ma H T, Cao J W, Bian W S 2016 J. Chem. Phys. 144 184302 doi: 10.1063/1.4948631
    [11] 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 doi: 10.1063/1.4955498
    [12] Yuan X, Yin S, Shen Y, Liu Y, Lian Y, Xu H F, Yan B 2018 J. Chem. Phys. 149 094306 doi: 10.1063/1.5044387
    [13] Yzombard P, Hamamda M, Gerber S, Doser M, Comparat D 2015 Phys. Rev. Lett. 114 213001 doi: 10.1103/PhysRevLett.114.213001
    [14] Zhang Q Q, Yang C L, Wang M S, Ma X G, Liu W W 2017 Spectrochim. Acta, Part A. 182 130 doi: 10.1016/j.saa.2017.04.006
    [15] Zhang Q Q, Yang C L, Wang M S, Ma X G, Liu W W 2017 Spectrochim. Acta, Part A. 185 365 doi: 10.1016/j.saa.2017.06.001
    [16] Wan M J, Huang D H, Yu Y, Zhang Y G 2017 Phys. Chem. Chem. Phys. 17 27360
    [17] 万明杰, 李松, 金成国, 罗华锋 2019 物理学报 68 063103 doi: 10.7498/aps.68.20182039 Wan M J, Li S, Jin C G, Luo H F 2019 Acta Phys. Sin. 68 063103 doi: 10.7498/aps.68.20182039
    [18] 万明杰, 罗华锋, 袁娣, 李松 2019 物理学报 68 173102 doi: 10.7498/aps.68.20190869 Wan M J, Luo H F, Yuan D, Li S 2019 Acta Phys. Sin. 68 173102 doi: 10.7498/aps.68.20190869
    [19] 万明杰, 柳福提, 黄多辉 2021 物理学报 70 033101 doi: 10.7498/aps.70.20201413 Wan M J, Liu F T, Huang D H, 2021 Acta Phys. Sin. 70 033101 doi: 10.7498/aps.70.20201413
    [20] Carlsten J L, Peterson J R. Lineberger W C 1976 Chem. Phys. Lett. 37 5 doi: 10.1016/0009-2614(76)80148-0
    [21] Miller T M, Leopold D G, Murray K K, Lineberger W C 1986 J. Chem. Phys. 85 2368 doi: 10.1063/1.451091
    [22] Jordan K D, Luken W 1976 J Chem. Phys. 64 2760 doi: 10.1063/1.432599
    [23] Li S, Zheng R, Chen S J, Fan Q C 2015 Mol. Phys. 113 1433 doi: 10.1080/00268976.2014.1001805
    [24] Weck P F, Kirby K, Stancil P C 2004 J. Chem. Phys. 120 4216 doi: 10.1063/1.1643715
    [25] Kurosaki Y, Yokoyama K 2012 J. Chem. Phys. 137 064305 doi: 10.1063/1.4742059
    [26] Werner H J, Knowles P J, Knizia G, et al. 2010 MOLPRO, a Package of ab initio Programs (Version 2010.1)
    [27] Knowles P J, Werner H J 1985 J. Chem. Phys. 82 5053 doi: 10.1063/1.448627
    [28] Knowles P J, Werner H J 1985 Chem. Phys. Lett. 115 259 doi: 10.1016/0009-2614(85)80025-7
    [29] Werner H J, Knowles P J 1988 J. Chem. Phys. 89 5803 doi: 10.1063/1.455556
    [30] Langhoff S R, Davidson E R 1974 Int. J. Quantum Chem. 8 61 doi: 10.1002/qua.560080106
    [31] Berning A, Schweizer M, Werner H J, Knowles P J, Palmieri P 2000 Mol. Phys. 98 1283
    [32] Xiao K L, Yang C L, Wang M S, Ma X G, Liu W W 2013 J. Chem. Phys. 139 074305 doi: 10.1063/1.4818452
    [33] Weigend F 2008 J. Comput. Chem. 29 167 doi: 10.1002/jcc.20702
    [34] Peterson K A, Dunning T H 2002 J. Chem. Phys. 117 10548 doi: 10.1063/1.1520138
    [35] 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)
    [36] Haeffler G, Hanstrorp D, Kiyan I, Klinkmueller A E, Ljungblad U, Pegg D 1996 Phys. Rev. A 53 4127 doi: 10.1103/PhysRevA.53.4127
    [37] Berzinsh U, Gustafsson M, Hanstorp D, Klinkmueller A E, Ljungblad U, Maartensson-Pendrill A M 1995 Phys. Rev. A 51 231 doi: 10.1103/PhysRevA.51.231
    [38] Moore C E 1971 Atomic Energy Levels (Vol. 1) (Washington, DC: US Govt Printing Office) pp9, 195
    [39] 尹俊豪, 杨涛, 印建平 2021 物理学报 70 163302 doi: 10.7498/aps.70.20210522 Yin J H, Yang T, Yin J P 2021 Acta Phys. Sin. 70 163302 doi: 10.7498/aps.70.20210522
  • 加载中
图( 5) 表( 7)
计量
  • 文章访问数:  247
  • HTML全文浏览数:  247
  • PDF下载数:  0
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-09-10
  • 刊出日期:  2022-02-20

LiCl–阴离子的光谱性质和跃迁性质

  • 收稿日期:  2021-09-10

摘要: 采用多参考组态相互作用方法结合全电子基组计算了LiCl- 阴离子5个电子态 (X2Σ+, A2∏, B2Σ+, 32Σ+, 22∏) 的电子结构. 为了得到精确的光谱常数, 计算中考虑了Davidson 修正、芯-价电子关联效应和自旋-轨道耦合效应. 拟合得到各电子态的光谱常数、分子常数、自发辐射速率和自发辐射寿命. 基态的光谱常数与实验值和其他理论值符合较好, 同时报道了LiCl– 阴离子激发态的光谱常数以及其到基态的跃迁性质. 计算结果表明A2∏ $\leftrightarrow $X2Σ+跃迁具有高对角分布的弗兰克-康登因子f00, 第一激发态A2∏有较短的自发辐射寿命. 构造A2∏ (ν′) $\leftrightarrow $ X2Σ+ ($\nu''$)准循环跃迁进行激光冷却LiCl– 阴离子需要一束主激光和两束抽运激光. 以上结果预测了激光冷却LiCl–阴离子是可行的.

English Abstract

    全文HTML

    1.   引　言

    • 超冷分子在量子计算[1]、化学动力学[2]和控制化学[3]等领域有着广泛的应用. 由于分子内部能级结构复杂, 激光冷却分子存在巨大的挑战. Shuman等[4]于2010年首次在实验上实现了分子的激光冷却. 随后Zhelyazkova 等[5]于2014年实现了YO分子的激光冷却, 他们构造了一个闭合的三电子能级系统. 至今为止, 有很多关于激光冷却氯化物和锂化物的理论研究, 如LiBe[6], BeCl[7], MgCl[8,9], CaCl[10], AlCl[11]和TlCl[12]分子.

      2015年, Yzombard等[13]首次从理论上预测了$ {\text{C}}_2^ - $双原子分子阴离子是适合激光冷却的候选体系. 杨传路等[14,15]也分别讨论了激光冷却NH–和BH–阴离子的可能性. 对两种阴离子而言, 他们分别构造了11∏ $\leftrightarrow $ X1Σ+和12Σ+ $\leftrightarrow $ 12∏准闭合的能级系统. 近几年来, 我们讨论了激光冷却OH– [16], SH– [17], KCl– [18]和SeH– [19]阴离子的可能性. LiCl–与KCl–阴离子的结构相似. 以往的文献中很少有LiCl–阴离子和LiCl分子关于势能性质方面的研究[20-23]. 1976年, 实验上第一次报道LiCl–阴离子的光谱[20], 并得到了其基态的平衡核间距. 1986年, Miller等[21]报道了LiCl–的光电子谱, 分辨率为80 cm–1 (0.01 eV). 实验得到了其基态的平衡核间距, 离解能等光谱常数, 同时得到了LiCl分子电子亲合能为(0.593 ± 0.01) eV. 1976年, Jordan和Luken[22]采用Hartree-Fock (HF)方法计算得到了LiCl–阴离子基态X2Σ+的势能曲线, 他们得到基态的平衡核间距略大于2.12 Å. 近年来, 李松等[23]采用耦合簇方法计算了其基态的光谱常数和分子常数. 2004年, Weck等[24]采用多参考单双激发组态相互作用 (MRSDCI) 方法计算得到了LiCl分子X1Σ+和B1Σ+的势能曲线. 两个电子态的势能曲线在8 Å 左右发生避免交叉现象. 2012年, Kurosaki和Yokoyama[25]采用MRSDCI方法计算得到了LiCl分子7个Λ-S态和13个Ω态的势能曲线, 并拟合得到其基态的光谱常数.

      从以往的研究中可以发现, 只有极少数文献研究了LiCl–阴离子基态的光谱性质. 本文将采用高精度的从头算方法计算LiCl–阴离子5个较低的Λ-S态和7个Ω态的电子结构, 得到每个束缚态的光谱常数和分子常数. 同时预测A2∏ $\leftrightarrow $ X2Σ+跃迁的弗兰克-康登因子, 爱因斯坦系数和自发辐射寿命. 最后讨论通过构造A2∏ $\leftrightarrow $ X2Σ+准循环跃迁能级进行激光冷却LiCl–阴离子的可能性.

    2.   计算细节

    • 基于MOLPRO 2010程序包[26]计算了LiCl–阴离子最低的3个离解通道所对应的5个电子态 (X2Σ+, A2∏, B2Σ+, 32Σ+, 22∏) 的势能曲线. 首先采用限制性的Hartree-Fock (HF) 方法得到LiCl–阴离子的初始能量; 然后通过完全活动空间自洽场方法 (CASSCF)[27,28] 产生多参考波函数; 最后采用多参考组态相互作用方法 (MRCI)[29,30] 得到所计算Λ-S态的能量, 考虑Davidson修正 (+Q) 减小计算误差. 在MRCI+Q水平下通过Breit-Pauli哈密顿量[31]来考虑自旋-轨道耦合效应 (SOC).

      由于程序的限制, 必须把C∞V点群约化为C2V子群来计算LiCl–阴离子的电子结构. C2V群有4个不可约表示 (a1, b1, b2, a2). 在CASSCF计算中, 选择了8个分子轨道 (4, 2, 2, 0) 作为活动空间, 包含了Li 2s2p和Cl的3s3p壳层, 可以写为CAS(9, 8). Li 1s和Cl 2s2p选作闭壳层, 保持双占据. 其他的电子被冻结, 不参与能量的计算. 在MRCI+Q计算中, Li 1s和Cl 2s2p用于考虑芯-价电子 (CV) 关联效应, 有17个电子参与了相关能的计算. 和K原子一样[32], Li原子采用了def2-AQZVPP-JKFI全电子基组[33], Cl原子采用了AWCVQZ-DK全电子基组[34].

      采用LEVVEL8.0程序[35]求解径向薛定谔方程得到每个束缚态的平衡核间距 (Re), 离解能 (De), 谐振频率 (ωe), 非谐振频率 (ωeχe), 转动常数 (Be)和垂直跃迁能 (Te). 同时计算出弗兰克-康登因子, 爱因斯坦系数Aν′ν″(s–1) 和自发辐射寿命. 爱因斯坦系数可以表示为

      Λ′和Λ″分别表示高态和低态的电子角动量投影量子数, (2 – δ0,Λ′Λ″)/(2 – δ0, Λ′)表示简并因子, ΔEν′ν″表示振动能级之间的能量差, D(r)表示跃迁偶极矩, ψν′和ψν″分别表示高态和低态的振动波函数.

    3.   结果与讨论

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

    • 在MRCI+Q水平下计算得到了LiCl–阴离子X2Σ+, A2∏, B2Σ+, 32Σ+和22∏态的势能曲线. 图1描绘了5个电子态的势能曲线. 可以看到5个电子态都是束缚态. Li和Cl原子的电子亲合能的实验值分别4984.86 cm–1 [36]和29138.35 cm–1[37]. 可以得到基态X2Σ+对应最低的离解极限Li(2Sg) + Cl–(1Sg), A2∏和B2Σ+电子态对应第二离解极限Li(2Pu) + Cl–(1Sg), 而32Σ+和22∏电子态对应第三离解极限Li–(1Sg)+Cl(2Pu). Λ-S态之间的离解关系以及离解极限之间的能量差列于表1中. 从表1可以看出, 本文的计算的结果和已有实验值符合较好. 同时本文计算了LiCl分子的电子亲合能为0.4392 eV, 比Miller等[21]的实验值小0.153 eV.

      LiCl–阴离子5个束缚态的光谱常数列在表2中. 基态X2Σ+的平衡核间距为2.1352 Å, 比最新实验值[21]大0.0112 Å, 误差仅为0.5%, 与最新理论值[23]相差仅0.0002 Å. 相比与同主簇的KCl–阴离子, 其基态的平衡核间距比KCl–阴离子约小25%[18]. 本文计算的谐振频率也仅比最新理论值小2.37 cm–1. 基态在平衡核间距处的主要电子组态为1σ22σ23σ24σ21π45σ26σ22π47σ18σ03π0, 权重为74.9%, 其另外一个重要电子组态为1σ22σ23σ24σ21π45σ26σ22π47σ08σ13π0, 权重为17.2%. 基态有一个很深的势阱, 其势阱深度达到了1.8556 eV. 从表1可以看出, 本文计算的基态的光谱常数与已有实验值和理论值符合很好. 第一激发态A2∏到基态的垂直跃迁能为13431.93 cm–1. 其平衡核间距为2.1198 Å, 比基态的平衡核间距小0.0154 Å. 其在平衡核间距处的主要电子组态为1σ22σ23σ24σ21π45σ26σ22π47σ08σ03π1, 权重为92.03%, A2∏ $\leftrightarrow $ X2Σ+ 跃迁源自7σ $\leftrightarrow $ 3π跃迁. B2Σ+, 32Σ+和22∏态在平衡核间距处的主要电子组态分别为1σ22σ23σ24σ21π45σ26σ22π47σ08σ13π0, 1σ22σ23σ24σ21π45σ26σ12π47σ28σ03π0和 1σ22σ23σ24σ21π45σ26σ22π37σ28σ03π0, 权重分别为74.5%, 80.16%和80.16%, 可以看出B2Σ+ → X2Σ+, 32Σ+ → X2Σ+和22∏ → X2Σ+ 跃迁分别来源于 7σ → 8σ, 6σ → 7σ和 3π → 7σ跃迁. 32Σ+和22∏态为弱束缚态, 其谐振频率分别为30.19 cm–1和39.08 cm–1. 其离解能仅为0.00353 eV和0.0112 eV. 从图1中可以, 看出两个态的势能曲线在约3.94 Å出现交叉.

      • 3.2.   振动能级和分子常数

    • 基于精确的势能曲线, 得到了J = 0时5个Λ-S态 (X2Σ+, A2∏, B2Σ+, 32Σ+和22∏) 的振动态, 5个态分别有65, 90, 31, 30和31个振动态, 限于篇幅, 表3中分别只列出了每个态前10个振动态的振动能级Gν和转动常数Bν. 2015年, 李松等[23]计算得到了6Li35Cl–阴离子基态的分子常数[23]. 本文中G0(X2Σ+)的计算结果和已有理论值的误差为1%, 基态的其他结果与理论值[23]的误差均小于0.4%. 可见本文的计算结果是可靠的.

      • 3.3.   Ω态的势能曲线和光谱常数

    • 考虑SOC效应后, Li和Cl原子的2P原子态分裂为2P1/2和2P3/2. 5个Λ-S态分裂为7个Ω态, 包含了5个1/2和2和3/2态, 分别对应5个最低的离解极限, Ω态之间的离解关系以及离解极限之间的能量差列于表4中. Li (2P1/2) + Cl– (1S0), Li (2P3/2) + Cl– (1S0), Li– (1S0) + Cl (2P3/2)和Li– (1S0) +Cl(2P1/2)与最低离解极限Li–(1S0)+Cl(2P1/2) 之间的能量差分别为14252.77 cm–1, 14253.50 cm–1, 23415.41 cm–1 和24288.48 cm–1, 与实验值符合[36-38]的误差分别为4.37%, 4.36%, 3.06%和2.99%. 计算得到Li (2P)和Cl (2P)的分裂值分别为0.26 cm–1和873.07 cm–1, 仅比实验值[38]小0.075 cm–1和9.28 cm–1. 从表4 可以看出本文的计算结果与实验值符合较好.

      $ {\text{X}}{}^2\Sigma _{1/2}^ +, {\text{A}}{}^2{\Pi _{1/2}}, {\text{A}}{}^2{\Pi _{3/2}} , {\text{B}}{}^2\Sigma _{1/2}^ + , 3{}^2\Sigma _{1/2}^ + , 2{}^2{\Pi _{{\text{3}}/2}} $和$ 2{}^2{\Pi _{1/2}} $ 态的势能曲线绘制在图2中, 可以看出$ {\text{A}}{}^2{\Pi _{1/2}} $和$ {\text{A}}{}^2{\Pi _{3/2}} $态的曲线基本重合. Ω态的光谱常数也列于表5中. 对比表2中没有考虑SOC的Λ-S态的光谱常数可以看出, SOC效应对X2Σ+, A2∏和B2Σ+态的光谱常数几乎没有影响, Ω态和Λ-S态的光谱常数最大误差不超过1%. A2∏态的分裂不明显, 计算得到其分裂常数仅为24.3 cm–1, 而22∏态的分裂非常明显, 其分裂常数达到了804 cm–1.

      本文计算了电子态之间的SO矩阵元, 通过其大小来判断SOC相互作用的强弱, 表6中列出了所计算电子态之间可能的16个SO矩阵元, 其中$ \left\langle {\text{A}} \right|\hat H_{{\text{SO}}}^{{\text{BP}}}\left| {\text{B}} \right\rangle $表示A和B态之间的SO矩阵元. SO矩阵元随核间距变化的曲线绘制在图3中. SO1和SO2表示X2Σ+和A2∏的两个分量之间的SO矩阵元, 故两者大小相等, 同样SO3 = SO4, SO5 = SO6, SO7 = SO8, SO9 = SO10, SO11 = SO12, SO14 = SO15. 图3中可以看到, 22∏的两个分量之间以及22∏与32Σ+之间的矩阵元达到了近300 cm–1, SOC效应对22∏态的影响较大, 而其他的SO矩阵元非常小, 说明SOC效应对X2Σ+, A2∏, B2Σ+电子态的影响很小, 基本可以忽略. 因此本文在自旋无关水平下研究了LiCl–阴离子的跃迁性质及激光冷却的可能性.

      • 3.4.   偶极矩, 跃迁偶极矩, 弗朗克-康登因子, 爱因斯坦系数

    • X2Σ+, A2∏, B2Σ+, 32Σ+和22∏态的电偶极矩以及激发态到基态跃迁的跃迁偶极矩随核间距变化的曲线分别描绘于图4(a)和图4(b)中. 基态X2Σ+平衡核间距处的电偶极矩达到了2.268 a.u, 比SrF分子 (1.378 a.u) 和KCl–阴离子 (1.212 a.u) 要大很多, 所以其极性很强. 也说明激光冷却LiCl–阴离子能满足电偶极矩较大的条件. 由于其是一个强极性的分子离子, 故其电子态的电偶极矩在无穷远处不为零, 从图4(a)中可以看出, 当核间距大于4 Å时, 电偶极矩随核间距的增加基本成线性增大的关系. 由于X2Σ+, A2∏和B2Σ+态来源于Li + Cl–, 故三个态的电偶极矩有同样的变化趋势, 32Σ+和22∏态来源于Li –+ Cl, 故这两个态的电偶极矩趋近于另一极限.

      由于A2∏和B2Σ+态具有同样的离解极限, A2∏ $\leftrightarrow $ X2Σ+和B2Σ+ $\leftrightarrow $ X2Σ+跃迁在核间距无穷远处源自于Li原子2Sg到2Pu的跃迁. 它们的跃迁偶极矩随着核间距的增大趋于同一值2.561 a.u. 在平衡核间距处, A2∏ $\leftrightarrow $ X2Σ+和B2Σ+ $\leftrightarrow $ X2Σ+跃迁有很大的跃迁偶极矩, 分别达到了2.393 a.u和1.712 a.u, 22∏ $\leftrightarrow $ X2Σ+和32Σ+ $\leftrightarrow $ X2Σ+跃迁比A2∏ $\leftrightarrow $ X2Σ+和B2Σ+ $\leftrightarrow $ X2Σ+跃迁的跃迁强度弱很多.

      本文计算了A2∏ $\leftrightarrow $ X2Σ+和B2Σ+ $\leftrightarrow $ X2Σ+跃迁的弗兰克-康登因子fν'ν'', 爱因斯坦系数和自发辐射寿命, 计算结果见表7. A2∏ $\leftrightarrow $ X2Σ+跃迁具有高对角分布的弗朗克-康登因子, 而B2Σ+ $\leftrightarrow $ X2Σ+跃迁对角分布的弗朗克-康登因子比较小. 两种跃迁都有较大的爱因斯坦系数, 当ν' = 0时, 两种跃迁总的爱因斯坦系数分别为2.821 × 107 s–1和2.858 × 107 s–1, A2∏和B2Σ+态的自发辐射寿命分别为35.45 ns和34.99 ns. 当ν' = 1时, 两种跃迁总的爱因斯坦系数分别为2.823 × 107 s–1和3.03 × 107 s–1, 自发辐射寿命分别为35.43 ns和33.0 ns.

      • 3.5.   激光冷却LiCl–阴离子的可行性

    • X2Σ+, A2∏和B2Σ+的势能曲线没有交叉和避免交叉, SOC效应对着3个态的势能曲线和跃迁性质的影响也不明显, 本文将在自旋无关水平下讨论激光冷却LiCl–阴离子的可能性.

      A2∏ $\leftrightarrow $ X2Σ+跃迁的弗朗克-康登因子f00 = 0.9898, 比SrF分子和CaH分子的f00略大 [f00(SrF) = 0.98[4], f00(CaH) = 0.985[39]], 也比同主簇KCl–阴离子的f00大很多 (f00 = 0.8816)[18], 说明激光冷却LiCl–阴离子具有足够大的f00. 另一方面, A2∏态的自发辐射寿命只有35.45 ns, 足以快速的激光冷却LiCl–阴离子.

      振动分支比Rν'ν''比fν'ν''能更加准确的描述循环过程中的光子损失. 振动分支比Rν'ν''可以表示为

      A2∏ $\leftrightarrow $ X2Σ+跃迁的对角项分支比R00和R11分别为0.9893和0.9685, 非对角项分支比分别为: R01 = 0.0106, R02 = 1.36 × 10–4, R03 = 1.58 × 10–6, R10 = 0.0096, R12 = 0.0215, R13 = 4.37 × 10–4, 可以得到R00, R01和R02的和非常接近1, 说明A2∏(ν′) $\leftrightarrow $ X2Σ+($v''$)跃迁是一个准循环跃迁能级系统. 其冷却途径见图5. 红线表示泵浦激光其中主激光波长λ00为744.10 nm, 为了提高冷却效率, 增加了两束抽运激光, 其波长分别为λ10 = 774.30 nm和λ21 = 772.42 nm. 激光冷却LiCl–阴离子所需激光波长要小于冷却KCl–阴离子所需波长[18]. 蓝线表示从A2∏自发辐射到基态的振动分支比. R03+的表示ν″ ≥ 3的所有分支比之和, 计算得到R03+ < 1.60 × 10–6, 则理论上至少能散射6 × 105个光子. 故可以构造A2∏(ν') $\leftrightarrow $ X2Σ+($v''$)准闭合能级系统进行激光冷却LiCl–阴离子.

      虽然B2Σ+态的自发辐射寿命足够小, 但由于B2Σ+ $\leftrightarrow $ X2Σ+跃迁的弗兰克-康登因子太小, f00 = 0.5908, 故没有考虑构建B2Σ+ $\leftrightarrow $ X2Σ+跃迁进行激光冷却LiCl–阴离子.

    4.   结　论

    • 本文采用MRCI+Q方法计算了LiCl–阴离子X2Σ+, A2Π, B2Σ+, 32Σ+和22Π 电子态的势能曲线, 电偶极矩和跃迁偶极矩. 计算过程中考虑了CV关联和SOC效应. 计算结果表明5个Λ-S态和7个Ω态都是束缚态. 通过求解径向薛定谔方程得到束缚态的光谱常数和分子常数. 基态的光谱常数与分子常数与已有实验值和理论值符合较好. 得到了A2∏ $\leftrightarrow $ X2Σ+和B2Σ+ $\leftrightarrow $ X2Σ+跃迁的弗兰克-康登因子, 爱因斯坦系数和自发辐射寿命. 本文报道了LiCl–阴离子激发态的光谱常数, 分子常数和跃迁性质. 通过分析SO矩阵可以得到SOC效应对X2Σ+, B2Σ+, A2∏的光谱常数和A2∏ $\leftrightarrow $ X2Σ+和B2Σ+ $\leftrightarrow $ X2Σ+跃迁的跃迁性质影响很小, 可以忽略. A2∏($v' $) $\leftrightarrow $ X2Σ+($v'' $)跃迁具有高对角分布的弗兰克-康登因子以及较小的自发辐射寿命, 同时得到了激光冷却所LiCl–阴离子需的激光波长. 在自旋无关水平下预测LiCl–阴离子是适合激光冷却的候选离子.

    参考文献 (39)

目录

/

返回文章
返回