搜索

x

留言板

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

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

锶原子光晶格钟黑体辐射频移评估

李婷 卢晓同 张强 孔德欢 王叶兵 常宏

锶原子光晶格钟黑体辐射频移评估

李婷, 卢晓同, 张强, 孔德欢, 王叶兵, 常宏
PDF
HTML
导出引用
  • 在中性原子光晶格钟的系统不确定度评估中, 通常黑体辐射引起的频移是最大的一项. 黑体辐射频移主要受周围环境温度的影响. 针对国家授时中心的锶原子光晶格钟实验系统, 通过理论分析、腔体表面温度的测量和软件模拟相结合的方法, 评估了锶原子光晶格钟黑体辐射频移的修正量和不确定度. 其中主要分析了锶原子炉、蓝宝石加热窗口、透过窗口片进入到真空腔体内的室温以及Zeeman减速装置对原子团处的热辐射引起的黑体辐射频移. 在真空腔体外表面设置了5个测温点, 利用校准过的铂电阻温度传感器监测真空腔体外表面的温度变化, 用SolidWorks绘图软件建立腔体模型, 通过有限元分析软件模拟出在真空腔体温度变化0.72 K时, 原子团所处位置温度的波动为0.34 K. 最终得到黑体辐射频移总的修正量为–2.13(1) Hz, 不确定度为2.4 × 10–17.
      通信作者: 王叶兵, wangyebing@ntsc.ac.cn ; 常宏, changhong@ntsc.ac.cn
    • 基金项目: 国家自然科学基金(批准号: 11803042, 11474282, 61775220)、中国科学院前沿科学重点研究项目(批准号: QYZDB-SSW-JSC004)和中国科学院国家授时中心青年创新人才项目资助的课题.
    [1]

    Derevianko A, Katori H 2011 Rev. Mod. Phys. 83 331

    [2]

    Bloom B J, Nicholson T L, Williams J R, Campbell S L, Bishof M, Zhang X, Zhang W, Bromley S L, Ye J 2014 Nature 506 71

    [3]

    Hinkley N, Sherman J A, Phillips N B, Schioppo M, Lemke N D, Beloy K, Pizzocaro M, Oates C W, Ludlow A D 2013 Science 341 1215

    [4]

    Nicholson T L, Campbell S L, Hutson R B, Marti G E, Bloom B J, McNally R L, Zhang W, Barrett M D, Safronova M S, Strouse G F, Tew W L, Ye J 2015 Nat. Commun. 6 6896

    [5]

    Ushijima I, Takamoto M, Das M, Ohkubo T, Katori H 2015 Nat. Photon. 9 185

    [6]

    Campbell S L, Hutson R B, Marti G E, Goban A, Oppong N D, Mcally R L, Sonderhouse L, Robinson J M, Zhang W, Bloom B J, Ye J 2017 Science 358 90

    [7]

    McGrew W F, Zhang X, Fasano1 R J, Schäffer S A , Beloy K, Nicolodi D, Brown R C, Hinkley N, Milani G, Schioppo M, Yoon T H, Ludlow A D 2018 Nature 564 87

    [8]

    Derevianko A, Pospelov M 2014 Nat. Phys. 10 933

    [9]

    Kolkowitz S, Pikovski I, Langellier N, Lukin M D, Walsworth R L, Ye J 2016 Phys. Rev. D 94 124043

    [10]

    Blatt S, Ludlow A D, Campbell G K, Thomsen J W, Zelevinsky T, Boyd M M, Ye J, Baillard X, Fouché M, Le Targat R, Brusch A, Lemonde P, Takamoto M, Hong F L, Katori H, Flambaum V V 2008 Phys. Rev. Lett. 100 140801

    [11]

    Godun R M, Nisbet-Jones P B R, Jones J M, King S A, Johnson L A M, Margolis H S, Szymaniec K, Lea S N, Bongs K, Gill P 2014 Phys. Rev. Lett. 113 210801

    [12]

    Huntemann N, Lipphardt B, Tamm C, Gerginov V, Weyers S, Peik E 2014 Phys. Rev. Lett. 113 210802

    [13]

    Ludlow A D, Boyd M M, Ye J, Peik E, Schmidt P O 2015 Rev. Mod. Phys. 87 637

    [14]

    Yudin V I, Taichenachev A V, Okhapkin M V, Bagayev S N, Tamm C, Peik E, Huntemann N, Mehlstäubler T E, Riehle F 2011 Phys. Rev. Lett. 107 030801

    [15]

    Beloy K, Hinkley N, Phillips N B, Sherman J A, Schioppo M, Lehman J, FeldmanA, Hanssen L M, Oates C W, Ludlow A D 2014 Phys. Rev. Lett. 113 260801

    [16]

    Wang Y B, Yin M J, Ren J, Xu Q F, Lu B Q, Han J X, Guo Y, Chang H 2018 Chin. Phys. B 27 023701

    [17]

    Wang Y B, Lu X T, Lu B Q, Kong D H, Chang H 2018 Appl. Sci. 8 2194

    [18]

    Itano W M, Lewis L L, Wineland D J 1982 Phys. Rev. A 25 1233

    [19]

    Hollberg L, Hall J L 1984 Phys. Rev. Lett. 53 230

    [20]

    Porsev S G, Derevianko A 2006 Phys. Rev. A 74 020502

    [21]

    Safronova M S, Porsev S G, Safronova U I, Kozlov M G, Clark C W 2013 Phys. Rev. A 87 012509

    [22]

    Middelmann T, Lisdat C, Falke S, Winfred J S R V, Riehle F, Sterr U 2011 IEEE Trans. Instrum. Meas. 60 2550

    [23]

    Middelmann T, Falke S, Lisdat C, Sterr U 2012 Phys. Rev. Lett. 109 263004

    [24]

    Falke S, Lemke N, Grebing C, Lipphardt B, Weyers S, Gerginov V, Huntemann N, Hagemann C, Al-Masoudi A, Häfner S, Vogt S, Sterr U, Lisdat C 2014 New J. Phys. 16 073023

  • 图 1  实验装置

    Fig. 1.  Experimental apparatus.

    图 2  一维光晶格

    Fig. 2.  Image of the one-dimensional optical lattice.

    图 3  87Sr原子钟跃迁谱线

    Fig. 3.  Spectra of the 87Sr clock transition.

    图 4  吸收系数对频移的影响

    Fig. 4.  Effect of the absorption coefficient on frequency shift.

    图 5  测温点的分布

    Fig. 5.  Distribution of the temperature points.

    图 6  测温点的温度波动

    Fig. 6.  Temperature fluctuations at the temperature points.

    图 7  网格划分结果

    Fig. 7.  Result of meshing.

    图 8  真空腔体内部温度分布

    Fig. 8.  Distribution of the internal temperature of the vacuum chamber.

    表 1  不同热源对应的频移修正量

    Table 1.  Frequency shift correction for different heat sources.

    热源频移修正量
    室温–0.0172(2) Hz
    原子炉–0.0031(4) Hz
    蓝宝石加热窗口–0.0013(1) Hz
    Zeeman减速窗口–0.00027(4) Hz
    下载: 导出CSV
  • [1]

    Derevianko A, Katori H 2011 Rev. Mod. Phys. 83 331

    [2]

    Bloom B J, Nicholson T L, Williams J R, Campbell S L, Bishof M, Zhang X, Zhang W, Bromley S L, Ye J 2014 Nature 506 71

    [3]

    Hinkley N, Sherman J A, Phillips N B, Schioppo M, Lemke N D, Beloy K, Pizzocaro M, Oates C W, Ludlow A D 2013 Science 341 1215

    [4]

    Nicholson T L, Campbell S L, Hutson R B, Marti G E, Bloom B J, McNally R L, Zhang W, Barrett M D, Safronova M S, Strouse G F, Tew W L, Ye J 2015 Nat. Commun. 6 6896

    [5]

    Ushijima I, Takamoto M, Das M, Ohkubo T, Katori H 2015 Nat. Photon. 9 185

    [6]

    Campbell S L, Hutson R B, Marti G E, Goban A, Oppong N D, Mcally R L, Sonderhouse L, Robinson J M, Zhang W, Bloom B J, Ye J 2017 Science 358 90

    [7]

    McGrew W F, Zhang X, Fasano1 R J, Schäffer S A , Beloy K, Nicolodi D, Brown R C, Hinkley N, Milani G, Schioppo M, Yoon T H, Ludlow A D 2018 Nature 564 87

    [8]

    Derevianko A, Pospelov M 2014 Nat. Phys. 10 933

    [9]

    Kolkowitz S, Pikovski I, Langellier N, Lukin M D, Walsworth R L, Ye J 2016 Phys. Rev. D 94 124043

    [10]

    Blatt S, Ludlow A D, Campbell G K, Thomsen J W, Zelevinsky T, Boyd M M, Ye J, Baillard X, Fouché M, Le Targat R, Brusch A, Lemonde P, Takamoto M, Hong F L, Katori H, Flambaum V V 2008 Phys. Rev. Lett. 100 140801

    [11]

    Godun R M, Nisbet-Jones P B R, Jones J M, King S A, Johnson L A M, Margolis H S, Szymaniec K, Lea S N, Bongs K, Gill P 2014 Phys. Rev. Lett. 113 210801

    [12]

    Huntemann N, Lipphardt B, Tamm C, Gerginov V, Weyers S, Peik E 2014 Phys. Rev. Lett. 113 210802

    [13]

    Ludlow A D, Boyd M M, Ye J, Peik E, Schmidt P O 2015 Rev. Mod. Phys. 87 637

    [14]

    Yudin V I, Taichenachev A V, Okhapkin M V, Bagayev S N, Tamm C, Peik E, Huntemann N, Mehlstäubler T E, Riehle F 2011 Phys. Rev. Lett. 107 030801

    [15]

    Beloy K, Hinkley N, Phillips N B, Sherman J A, Schioppo M, Lehman J, FeldmanA, Hanssen L M, Oates C W, Ludlow A D 2014 Phys. Rev. Lett. 113 260801

    [16]

    Wang Y B, Yin M J, Ren J, Xu Q F, Lu B Q, Han J X, Guo Y, Chang H 2018 Chin. Phys. B 27 023701

    [17]

    Wang Y B, Lu X T, Lu B Q, Kong D H, Chang H 2018 Appl. Sci. 8 2194

    [18]

    Itano W M, Lewis L L, Wineland D J 1982 Phys. Rev. A 25 1233

    [19]

    Hollberg L, Hall J L 1984 Phys. Rev. Lett. 53 230

    [20]

    Porsev S G, Derevianko A 2006 Phys. Rev. A 74 020502

    [21]

    Safronova M S, Porsev S G, Safronova U I, Kozlov M G, Clark C W 2013 Phys. Rev. A 87 012509

    [22]

    Middelmann T, Lisdat C, Falke S, Winfred J S R V, Riehle F, Sterr U 2011 IEEE Trans. Instrum. Meas. 60 2550

    [23]

    Middelmann T, Falke S, Lisdat C, Sterr U 2012 Phys. Rev. Lett. 109 263004

    [24]

    Falke S, Lemke N, Grebing C, Lipphardt B, Weyers S, Gerginov V, Huntemann N, Hagemann C, Al-Masoudi A, Häfner S, Vogt S, Sterr U, Lisdat C 2014 New J. Phys. 16 073023

  • [1] 卢晓同, 李婷, 孔德欢, 王叶兵, 常宏. 锶原子光晶格钟碰撞频移的测量. 物理学报, 2019, 68(23): 233401. doi: 10.7498/aps.68.20191147
    [2] 吴 翊, 荣命哲, 杨 飞, 王小华, 马 强, 王伟宗. 引入6波段P-1辐射模型的三维空气电弧等离子体数值分析. 物理学报, 2008, 57(9): 5761-5767. doi: 10.7498/aps.57.5761
    [3] 郭小云, 石才土, 张久昶, 辛洪兵. 永磁扭摆磁铁的同步辐射特性和结构分析. 物理学报, 2006, 55(4): 1731-1735. doi: 10.7498/aps.55.1731
    [4] 李树, 李刚, 田东风, 邓力. 热辐射输运问题的隐式蒙特卡罗方法求解. 物理学报, 2013, 62(24): 249501. doi: 10.7498/aps.62.249501
    [5] 史旺林, 刘兴业, 刘振兴. Vaidya-Bonner-de Sitter黑洞对狄拉克粒子的热辐射. 物理学报, 2004, 53(7): 2396-2400. doi: 10.7498/aps.53.2396
    [6] 李树, 邓力, 田东风, 李刚. 基于能量密度分布的辐射源粒子空间抽样方法研究. 物理学报, 2014, 63(23): 239501. doi: 10.7498/aps.63.239501
    [7] 郭阳, 尹默娟, 徐琴芳, 王叶兵, 卢本全, 任洁, 赵芳婧, 常宏. 锶原子光晶格钟自旋极化谱线的探测. 物理学报, 2018, 67(7): 070601. doi: 10.7498/aps.67.20172759
    [8] 孙健, 刘伟强. 内嵌定向高导热层疏导式结构热防护机理分析. 物理学报, 2012, 61(12): 124401. doi: 10.7498/aps.61.124401
    [9] 钟广明, 杜晓晴, 唐杰灵, 董向坤, 雷小华, 陈伟民. 影响倒装焊LED芯片电流分布均匀性的因素分析. 物理学报, 2012, 61(12): 127803. doi: 10.7498/aps.61.127803
    [10] 刘李辉, 邹宏新, 刘曲, 李玺. 199Hg+光频标的黑体辐射频移. 物理学报, 2012, 61(10): 103101. doi: 10.7498/aps.61.103101
    [11] 万 红, 沈仁发, 吴学忠. 对称磁电层合板磁电转换效应理论研究. 物理学报, 2005, 54(3): 1426-1430. doi: 10.7498/aps.54.1426
    [12] 万 红, 谢立强, 吴学忠, 刘希从. TbDyFe/PZT层状复合材料的磁电效应研究. 物理学报, 2005, 54(8): 3872-3877. doi: 10.7498/aps.54.3872
    [13] 赵建涛, 冯国英, 杨火木, 唐淳, 陈念江, 周寿桓. 薄片激光器热效应及其对输出功率的影响. 物理学报, 2012, 61(8): 084208. doi: 10.7498/aps.61.084208
    [14] 王存海, 郑树, 张欣欣. 非规则形状介质内辐射-导热耦合传热的间断有限元求解. 物理学报, 2020, 69(3): 034401. doi: 10.7498/aps.69.20191185
    [15] 钱治文, 商德江, 孙启航, 何元安, 翟京生. 三维浅海下弹性结构声辐射预报的有限元-抛物方程法. 物理学报, 2019, 68(2): 024301. doi: 10.7498/aps.68.20181452
    [16] 赵运进, 田锰, 黄勇刚, 王小云, 杨红, 米贤武. 基于有限元法的光子并矢格林函数重整化及其在自发辐射率和能级移动研究中的应用. 物理学报, 2018, 67(19): 193102. doi: 10.7498/aps.67.20180898
    [17] 林弋戈, 方占军. 锶原子光晶格钟. 物理学报, 2018, 67(16): 160604. doi: 10.7498/aps.67.20181097
    [18] 张登玉, 郭 萍, 高 峰. 强热辐射环境中两能级原子量子态保真度. 物理学报, 2007, 56(4): 1906-1910. doi: 10.7498/aps.56.1906
    [19] 周旺民, 蔡承宇, 王崇愚, 尹姝媛. 埋置量子点应力分布的有限元分析. 物理学报, 2009, 58(8): 5585-5590. doi: 10.7498/aps.58.5585
    [20] 张才根, 张幼文. 环境辐射对目标热辐射特性测试的影响. 物理学报, 1981, 30(7): 953-961. doi: 10.7498/aps.30.953
  • 引用本文:
    Citation:
计量
  • 文章访问数:  1217
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-28
  • 修回日期:  2019-02-23
  • 上网日期:  2019-05-01
  • 刊出日期:  2019-05-05

锶原子光晶格钟黑体辐射频移评估

  • 1. 中国科学院国家授时中心, 时间频率基准重点实验室, 西安 710600
  • 2. 中国科学院大学天文与空间科学学院, 北京 100049
  • 3. 河南师范大学物理与材料科学学院, 新乡 453007
  • 通信作者: 王叶兵, wangyebing@ntsc.ac.cn ; 常宏, changhong@ntsc.ac.cn
    基金项目: 国家自然科学基金(批准号: 11803042, 11474282, 61775220)、中国科学院前沿科学重点研究项目(批准号: QYZDB-SSW-JSC004)和中国科学院国家授时中心青年创新人才项目资助的课题.

摘要: 在中性原子光晶格钟的系统不确定度评估中, 通常黑体辐射引起的频移是最大的一项. 黑体辐射频移主要受周围环境温度的影响. 针对国家授时中心的锶原子光晶格钟实验系统, 通过理论分析、腔体表面温度的测量和软件模拟相结合的方法, 评估了锶原子光晶格钟黑体辐射频移的修正量和不确定度. 其中主要分析了锶原子炉、蓝宝石加热窗口、透过窗口片进入到真空腔体内的室温以及Zeeman减速装置对原子团处的热辐射引起的黑体辐射频移. 在真空腔体外表面设置了5个测温点, 利用校准过的铂电阻温度传感器监测真空腔体外表面的温度变化, 用SolidWorks绘图软件建立腔体模型, 通过有限元分析软件模拟出在真空腔体温度变化0.72 K时, 原子团所处位置温度的波动为0.34 K. 最终得到黑体辐射频移总的修正量为–2.13(1) Hz, 不确定度为2.4 × 10–17.

English Abstract

    • 原子钟的性能主要包括稳定度和不确定度两个方面[1]. 基于中性原子的光晶格钟[2-5], 有大量的原子样品作参考, 可以在更短的时间内达到较高的稳定度. 2017年, 美国天体物理联合实验室的锶原子光晶格钟[6]测量精度达到了5 × 10–19. 2018年, 美国标准技术研究院的镱原子光钟[7]系统的不确定度达到了1.4 × 10–18. 随着光钟性能指标的不断提高, 其在暗物质探测[8]、引力波测量[9]以及基本物理常数的测量[10-12]等方面有广泛的应用前景.

      在中性原子光晶格钟不确定度性能评估中, 通常情况下黑体辐射引起的频移为最大的一项[13]. 它的大小与动/静态修正因数、环境温度以及原子的极化率等密切相关[14]. 国际上对黑体辐射频移不确定度的评估一般有两种方法. 第一种是控制真空腔体周围环境温度的波动, 通过对腔体内部温度的精确测量, 计算出频率偏移修正量和不确定度. 如2014年美国天体物理联合实验室的锶原子光晶格钟[2], 他们把真空腔体放入一个黑体箱中并进行精确的主动控温, 阻挡了大部分的环境温度波动对真空腔体的影响. 然后使用两个内置在真空腔体中的硅二极管温度传感器测温, 描绘了原子团周围的环境温度分布图, 黑体辐射频移不确定度为2 × 10–18. 2014年, 美国标准技术研究院的镱原子光晶格钟[15]在真空腔外安装了一个镀膜的黑体辐射屏蔽腔, 阻挡了周围的黑体辐射, 同时减小了环境温度对原子团处温度的影响, 黑体辐射频移不确定度达到了1 × 10–18. 第二种方法是通过低温冷冻的方法降低原子团处的温度, 根据黑体辐射频移与环境温度的四次方成反比的关系[14], 降低原子团所处位置的温度, 可以大幅度地降低黑体辐射频移和不确定度, 然后通过测量计算得到频率偏移修正量和不确定度. 如2015年日本东京大学的研究组采用低温冷冻的方法[5], 通过运动光晶格把冷原子团转移到温度为95 K的低温冷冻腔里进行探测, 根据测量计算得到黑体辐射的频移不确定度为0.9 × 10–18. 相比之下, 第二种方法比第一种方法更为复杂, 操作难度更大, 但是第二种方法可以从根本上降低黑体辐射频移, 并且容易降低黑体辐射频移的不确定度.

      国家授时中心于2017年实现了锶冷原子光晶格钟的闭环锁定[16]. 目前, 系统闭环不稳定度为1.6 × 10–15/$\sqrt \tau $, 在2000 s的平均时间内环内不稳定度[17]为2.8 × 10–17. 在对系统的黑体辐射频移评估时, 由于在真空腔体内部没有安装内置的测温装置, 无法精确测量出原子团所处位置的温度, 只能通过测量腔体表面的温度, 再根据腔体外表面的温度推测原子团所处位置的温度. 然而, 在锶原子光晶格钟正常运行过程中, 实验环境和真空腔体周围其他的热源使真空腔体周围的环境变得复杂, 真空腔体外表面的温度不能真实反映出真空腔体内的温度分布. 我们采用测量腔体表面温度与软件模拟相结合的方法, 得到原子团处的温度. 本文就黑体辐射频移不确定度评估实验和理论部分进行详细说明. 通过监测系统稳定运行过程中真空腔外表面的温度, 根据真空腔体外表面的平均温度利用ANSYS有限元分析软件模拟出真空腔体内原子团所处位置的温度, 由原子团处的温度及温度波动计算相应的频移修正量和不确定度.

    • 任何非绝对零度的物体都会向周围辐射电磁波[13]. 在光钟正常运行过程中, 原子所处的环境温度并不是绝对零度, 辐射能量使得原子钟跃迁上下能级发生不同程度的位移, 最终导致钟跃迁绝对频率的移动, 该频移为黑体辐射频移[18-21].

    • 锶原子光晶格钟的整个物理系统的构成如图1所示. 图1从左至右: 1为锶炉, 常温下锶原子为固体, 在超冷原子87Sr光钟实验装置中, 锶炉被加热到500 ℃, 锶炉与磁光阱真空腔在水平方向上有一个夹角, 减小了热辐射源对磁光阱真空腔体内原子团的热辐射; 2为准直器; 3为二维准直窗口; 4和6分别为40 L/s和60 L/s的离子泵; 5为Zeeman减速器; 7为磁光阱真空腔, 该部分共有十四个窗口, 其中最左端的窗口与Zeeman 减速器相连, 最右端的窗口8为蓝宝石窗口, 为避免锶原子光晶格钟系统长时间运行时锶原子堆积在窗口片上, 窗口外缠绕了加热丝, 在实验时将其加热到120 ℃. 这个窗口和真空磁光阱之间用一个三通连接, 使加热窗口对磁光阱真空腔内原子团的热辐射减小. 磁光阱真空腔上下是一对反亥姆霍兹线圈, 为了减小线圈通电产生的热量, 避免线圈发热对光钟系统产生影响, 线圈结构采用浸润式, 将两个线圈接入水冷装置中, 实现对线圈的温度控制, 线圈整体温度控制为与室温一致.

      图  1  实验装置

      Figure 1.  Experimental apparatus.

      图2为通过电荷耦合器件(CCD)探测到真空腔体中的一维光晶格囚禁原子的图像, 光晶格的长度约为1 mm, 势阱深度为56Er (Er是晶格光子反冲能量), 囚禁的原子数目大约为3.6 × 104. 在模拟时, 把光晶格中的原子团视为一个直径为1 mm的球体. 图3为边带可分辨的钟跃迁谱线, 通过红蓝边带的面积比可以估算势阱中原子的温度约为4.2 μK. 图3右上角的插图是在可分辨边带谱的基础上加了补偿磁场, 使得原子团周围的磁场几乎为零, 降低钟激光的功率至150 nW, 在载波跃迁附近扫谱得到的窄线宽钟跃迁简并谱线.

      图  2  一维光晶格

      Figure 2.  Image of the one-dimensional optical lattice.

      图  3  87Sr原子钟跃迁谱线

      Figure 3.  Spectra of the 87Sr clock transition.

    • 在锶原子光晶格钟正常运行阶段, 主要热源为锶炉、Zeeman减速器、透过窗口片进入真空腔体内的室温和蓝宝石窗口, 它们对原子团的辐射可以根据立体角计算得出. 立体角计算公式为

      $S = 2{\text{π}}{R^2}\left( {1 - {\rm{cos}}\theta } \right),\;\varOmega = \frac{S}{{{R^2}}} = 2{\text{π}}\left( {1 - {\rm{cos}}\theta } \right).$

      原子团在水平方向上到窗口的距离为R1, 在竖直方向上到窗口的距离为R2, 原子团到锶炉的距离为R3, 蓝宝石窗口到原子团的距离为R4, Zeeman窗口到原子团的距离为R5, 每个普通窗口的半径为${r_1}$, 蓝宝石窗口的半径为${r_2}$, 原子炉处的半径为${r_3}$, Zeeman减速窗口的半径为${r_5}$. 根据(1)式可知, 一个普通窗口在水平方向上对应原子团的立体角为$0.01027{\text{π}}\;{\rm{sr}}$, 竖直方向上对应原子团的立体角为$0.02994{\text{π}}\;{\rm{sr}}$, 真空腔体内壁对应原子团的立体角Ω1$3.83604{\text{π}}\;{\rm{sr}}$, 普通窗口对应原子团总的立体角Ω2$0.16258{\text{π}}\;{\rm{sr}}$, 锶炉对应原子团的立体角Ω3$0.00009{\text{π}}\;{\rm{sr}}$, 蓝宝石窗口对应原子团的立体角Ω4$0.00078{\text{π}}\;{\rm{sr}}$, Zeeman减速窗口对应原子团的立体角Ω5$0.00051{\text{π}}\;{\rm{sr}}$.

    • 锶原子光晶格钟的黑体辐射频移可以表示为[14]

      $\Delta {V_{{\rm{BBR}}}}(T) = {V_{{\rm{stat}}}}{\left( {\frac{T}{{{T_0}}}} \right)^4} + {V_{{\rm{dyn}}}}{\left( {\frac{T}{{{T_0}}}} \right)^6},$

      式中T为环境温度, ${T_0} = 300$ K; ${V_{{\rm{stat}}}}$为静态频移系数, ${V_{{\rm{dyn}}}}$为动态频移系数. 其中德国联邦物理技术研究院测量得到[22,23]${V_{{\rm{stat}}}} = - $2.13023(6) Hz, 美国天体物理联合实验室测量得到[4]${V_{{\rm{dyn}}}}= $$ -0.1487(7)$ Hz. 由于辐射源有多个, 每处的温度也不相同, (2)式按照加权平均的方法得到适用于我们实验装置的黑体辐射计算公式为

      $\begin{split} \Delta {V_{{\rm{BBR}}}}(T) = &\; \frac{{{\varOmega _1}}}{{4{\text{π}}}} \times \Delta {V_{{\rm{BBR}}}}({T_1}) + \frac{{{\varOmega _2}}}{{4{\text{π}}}} \times \Delta {V_{{\rm{BBR}}}}({T_2}) \\ & + \frac{{{\varOmega _3}}}{{4{\text{π}}}}\times \Delta {V_{{\rm{BBR}}}}({T_3}) + \frac{{{\varOmega _4}}}{{4{\text{π}}}} \times \Delta {V_{{\rm{BBR}}}}({T_4})\\ & + \frac{{{\varOmega _5}}}{{4{\text{π}}}} \times \Delta {V_{{\rm{BBR}}}}({T_5}), \end{split}$

      其中, 真空腔体温度为T1, 室温为T2, 原子炉窗口温度为T3, 蓝宝石窗口温度为T4, Zeeman减速窗口温度为T5.

      根据(3)式计算出光钟在稳定运行时各个辐射源处的频移修正量(如表1所列, 括号内的数字是通过温度和相应的温度波动代入(2)式计算得到的): 1)室温可以透过玻璃窗口进入真空腔体内, 窗口的玻璃片上有镀膜, 一部分能量能够经过透射进入真空腔体[24], 在模拟时取窗口片的吸收系数为0.8, 透射系数小于0.2, 室温为294.89 K, 其波动值为0.72 K, 则对原子团的辐射引起的频移为–0.0172(2) Hz; 2)原子炉的温度为773.15 K, 其波动值为20 K, 则对原子团的辐射引起的频移为–0.0031(4) Hz; 3)蓝宝石加热窗口的温度为393.15 K, 其波动值为5 K, 则对原子团的辐射引起的频移为–0.0013(1) Hz; 4)把测温点1的温度视为Zeeman减速窗口的温度, 其波动值为10 K, 则对原子团的辐射引起的频移为–0.00027(4) Hz.

      热源频移修正量
      室温–0.0172(2) Hz
      原子炉–0.0031(4) Hz
      蓝宝石加热窗口–0.0013(1) Hz
      Zeeman减速窗口–0.00027(4) Hz

      表 1  不同热源对应的频移修正量

      Table 1.  Frequency shift correction for different heat sources.

    • 在分析真空腔体时, 先假设真空腔体内表面的吸收系数为a, 玻璃窗口的吸收系数为b, 假设真空腔体内壁的漫反射和能量密度在空间上是均匀且各向同性的, 根据能量守恒方程, 真空腔体吸收的能量等于辐射所散发出的能量[22,23], 得到:

      ${\varPhi _1} = {\text{π}}r_3^2c\rho({T_{{\rm{room}}}}),\tag{4a}$

      ${\varPhi _2} ={\text{π}}(12r_1^2 + r_2^2)bc\rho ({T_{{\rm{room}}}}),\tag{4b}$

      ${\varPhi _3} = {\text{π}}(2r_5^2 + 2{r_5}h - r_3^2 - 12r_1^2 - r_2^2)ac\rho ({T_{{\rm{cav}}}}),\tag{4c}$

      ${\varPhi _{\rm{4}}} = - {\rm{{\text{π}}(}}12r_1^2 + r_2^2{\rm{)}}bc{\rho _{{\rm{cav}}}},\tag{4d}$

      ${\varPhi _5} = - {\text{π}}(2r_5^2 + 2{r_5}h - r_3^2 - 12r_1^2 - r_2^2)ac{\rho _{{\rm{cav}}}},\tag{4e}$

      ${\varPhi _{\rm{6}}} = - {\text{π}}r_3^2c{\rho _{{\rm{cav}}}},\tag{4f}$

      ${\varPhi _1}$为通过Zeeman减速器与真空腔体连接处进入真空腔体内的辐射散发的能量, ${\varPhi _2}$为玻璃窗口内表面辐射散发的能量, ${\varPhi _3}$为真空腔体内表面辐射散发的能量, ${\varPhi _4}$为窗口内表面吸收的能量, ${\varPhi _5}$为真空腔体内表面吸收的能量, ${\varPhi _6}$为Zeeman减速器与真空腔体连接处真空腔体泄漏的能量, c为光速, ${r_5}$为圆柱型真空腔体的直径, h为圆柱型腔体的高. 当达到平衡状态时,

      $({\varPhi _1}+{\varPhi _2}+{\varPhi _3})+({\varPhi _4} +{\varPhi _5}+{\varPhi _6}) = 0,$

      能量密度为

      ${\rho _{{\rm{cav}}}} = \frac{{\left( {2r_5^2 + 2{r_{\rm{5}}}h - r_3^2 - 12r_1^2 - r_2^2} \right)a\rho ({T_{{\rm{cav}}}}) + \left( {12r_1^2 + r_2^2} \right)b\rho ({T_{{\rm{room}}}}) + r_3^2\rho ({T_{{\rm{room}}}})}}{{\left( {2r_5^2 + 2{r_{\rm{5}}}h - r_3^2 - 12r_1^2 - r_2^2} \right)a + \left( {12r_1^2 + r_2^2} \right)b + r_3^2}}.$

      测量的腔体温度${T_{{\rm{cav}}}}$为295.07 K, $ {T_{{\rm{room}}}}$为294.89 K, 当$a = 0.9,\;b = 0.8$时, 由$\Delta v = - \Delta \alpha {\rho _{{\rm{cav}}}}/2h{\varepsilon _0}$可知, 真空腔体内表面吸收系数和玻璃窗口内表面吸收系数对频移的影响如图4所示.

      图  4  吸收系数对频移的影响

      Figure 4.  Effect of the absorption coefficient on frequency shift.

      图4(a)中, 红色曲线表示b = 0.8时吸收系数a与频移的关系, 即当玻璃的吸收系数为0.8时, 由腔体吸收系数a的取值引起的频移变化量为0.0035 Hz; 图4(b)中蓝色曲线表示a = 0.9时, 吸收系数b与频移的关系, 即当腔体的吸收系数为0.9时, 由玻璃吸收系数b引起的频移变化量为0.00035 Hz. 由于吸收系数引起的频移变化量所占的比例不足总频移量的千分之二, 所以由组成腔体的材料不同所造成的频移可以忽略, 即腔体内表面各部分吸收系数的大小对频移的影响可以忽略.

    • 利用经过中国计量院校准的铂电阻(PT100)和数字万用表(Agilent 34411A)进行腔体表面温度的测量. 在锶原子光晶格钟正常运行时, 在真空腔体主腔周围设置了5个温度测温点, 见图5, 测温点1位于靠近Zeeman减速器方向, 测温点3在蓝宝石窗口方向, 测温点2和测温点4位于腔体两侧, 测温点5悬挂真空腔体周围, 用于监测室温.

      图  5  测温点的分布

      Figure 5.  Distribution of the temperature points.

      在锶原子光晶格钟运行时, 真空腔体温度达到平衡时各个测温点的温度有差异, 具体如图6所示. 测温点1的平均温度为22.06 ℃, 测温点2的平均温度为21.74 ℃, 测温点3的平均温度为22.23 ℃, 测温点4的平均温度为21.65 ℃, 测温点5的温度波动为0.72 ℃. 由多次温度测量结果可知, 靠近锶炉和Zeeman减速装置方向的窗口的温度和靠近蓝宝石的窗口的温度较高.

      图  6  测温点的温度波动

      Figure 6.  Temperature fluctuations at the temperature points.

      通过有限元分析法得到真空腔内的温度分布. 温度模拟的目的是为了得到环境温度变化时原子团处的温度波动情况. 对真空腔体进行分析时, 需要建立一个封闭的腔体, 把Zeeman减速装置连接处和蓝宝石窗口连接处用两个普通的玻璃窗口来替代, 然后用SolidWorks绘图软件建立了封闭的真空腔体的模型.

      我们选用了热稳态分析的分析类型, 将利用SolidWorks绘图软件建立的腔体模型导入到ANSYS软件中. 首先, 定义材料属性, 腔体的材料是316L不锈钢, 普通玻璃窗口的材料为K9玻璃, Zeeman减速窗口的材料为蓝宝石. 在对材料进行热分析过程中, 需要定义与温度有关的量, 如密度、比热和导热系数等. 分析真空腔体内部的温度分布时, 先将腔体内部的真空部分填充为实体, 通过调整其导热系数、热值和密度等参数使填充实体特性接近真空. 原子团部分用一个直径为1 mm的球体填充(即球体为黑体), 此时, 球体的温度就等同于原子团处的温度. 为了保证在各个实体间有热传递, 对模型实体进行布尔运算操作, 使整个腔体为一个整体. 然后, 对整个腔体进行有限元网格划分. 在划分网格时, 较为复杂的模型通常采用自动网格划分的方法. 由于本文实验模型是不规则、且较为复杂的, 所以选用了自动网格划分. 通过设置各部分实体之间的关联度和网格尺寸, 得到的网格划分结果如图7所示.

      图  7  网格划分结果

      Figure 7.  Result of meshing.

      根据对热源的分析, 我们的模型除了室温辐射外其他辐射暂不考虑, 锶炉、Zeeman减速装置和蓝宝石窗口的辐射单独计算. 由于测温点1靠近Zeeman减速装置, 测温点3靠近蓝宝石加热窗口, 选取测温点1和测温点3正对的法兰窗口作为恒定温度载荷. 在相同时间里温度波动较大的地方表明与空气发生热交换的速度快, 表现为对流换热系数大. 根据多次不同位置的测温结果, 在腔体侧表面和上下表面施加不同的对流系数. 由多次实际测温点温度与对应的模拟点温度的对比, 最终确定施加在内侧表面的对流系数为4 W/(m2·K), 外侧表面的对流系数为2.7 W/(m2·K), 施加在腔体上下表面施加对流系数为1.5 W/(m2·K). 在真空腔体内表面和玻璃窗口施加辐射载荷, 辐射系数分别取为0.9和0.8, 黑体小球表面设置的辐射系数为1, 初始温度设置为293.15 K. 当我们得到的模拟值和测量值的温度差不超过0.1 K时, 认为选取的模型是合理的. 我们取小球的平均温度作为原子团所处位置的温度${T_{{\rm{eff}}}}$, 通过ANSYS Workbench的求解及后处理, 得到的模拟结果为: ${T_{{\rm{eff}}}}$ = 294.57 K, 相应的频移值为–2.11(1) Hz. 真空腔体内部温度分布如图8所示, 其中左下角的插图为腔体中心处的原子团的温度分布. 从图8内插图可以看出原子团处的温度是不均匀的, 这是由于真空腔体周围的环境温度不是恒温引起的.

      图  8  真空腔体内部温度分布

      Figure 8.  Distribution of the internal temperature of the vacuum chamber.

      改变环境温度得到的原子团处的温度差值记为$\Delta {T_{{\rm{eff}}}}$. 当环境温度变化0.72 K时, $\Delta {T_{{\rm{eff}}}}$ = 0.34 K. 其引起的不确定度为2.3 × 10–17, 加上原子炉和Zeeman减速窗口, 蓝宝石窗口以及室温透过玻璃窗口的影响, 得到的黑体辐射频移总的修正量为–2.13(1) Hz, 总的不确定度为2.4 × 10–17.

      在以后的实验改进方面, 我们会将真空腔体置于恒温环境中, 或者置于黑体辐射腔里, 可以提高腔体周围环境温度的均匀性, 并且进一步提高软件模拟的精确度.

    • 黑体辐射频移系统不确定度主要受原子团处温度波动的影响. 在真空腔体内没有安装内置测温装置的情况下, 为了更精确地测量黑体辐射频移的修正量及不确定度, 选择了理论分析计算、腔体表面温度的精确测量和软件模拟相结合的方法. 首先对热源和辐射系数进行了分析, 然后通过对真空腔体外表面温度的监测得到测温点温度的平均值, 最后利用ANSYS有限元分析软件对真空腔体内部的温度分布进行了模拟. 得到原子团周围的温度为294.57 K, 在腔体周围环境温度改变0.72 K时, 原子团所处的位置的温度变化值为0.34 K. 最终得出总的黑体辐射频移修正量为–2.13(1) Hz, 不确定度为2.4 × 10–17.

参考文献 (24)

目录

    /

    返回文章
    返回