空间锶原子光钟磁场分析及主动补偿系统

空间锶原子光钟磁场分析及主动补偿系统 Magnetic field analysis and active compensation system for strontium optical lattice clock in space 1 2 first-author 任　洁（1986-），女，陕西西安人，博士研究生，工程师，2013年于中国科学院大学获得硕士学位，主要从事锶原子光频标方面的研究。E-mail： renjie@ntsc.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777772&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777756&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777774&type= 22.01332855 32.00399780 任　洁 （1986-），女，陕西西安人，博士研究生，工程师，2013年于中国科学院大学获得硕士学位，主要从事锶原子光频标方面的研究。E-mail： renjie@ntsc.ac.cn renjie@ntsc.ac.cn 1 1 2 3 corresp E-mail： liuhui_gzs@nwu.edu.cn E-mail： liuhui_gzs@nwu.edu.cn 刘　辉（1986-），男，湖南郴州人，博士后，副教授，2016年于中国科学院大学获得博士学位， 2016至2019年在法国国家科研中心进行博士后研究，主要从事激光频率合成与超快激光技术及应用的研究。E-mail：liuhui_gzs@nwu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777788&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777794&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777793&type= 22.01332855 32.00400543 刘　辉 （1986-），男，湖南郴州人，博士后，副教授，2016年于中国科学院大学获得博士学位， 2016至2019年在法国国家科研中心进行博士后研究，主要从事激光频率合成与超快激光技术及应用的研究。E-mail： liuhui_gzs@nwu.edu.cn liuhui_gzs@nwu.edu.cn 1 2 corresp changhong@ntsc.ac.cn changhong@ntsc.ac.cn 常　宏（1977-），男，山西太原人，研究员，博士生导师，2000年、2005年于山西大学分别获得学士、博士学位，2005年至2007年在法国国家科研中心进行博士后研究，主要从事锶原子光频标方面的研究。E-mail：changhong@ntsc.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777779&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777795&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777781&type= 22.01332474 32.00399017 常　宏 （1977-），男，山西太原人，研究员，博士生导师，2000年、2005年于山西大学分别获得学士、博士学位，2005年至2007年在法国国家科研中心进行博士后研究，主要从事锶原子光频标方面的研究。E-mail： changhong@ntsc.ac.cn changhong@ntsc.ac.cn 中国科学院 国家授时中心 时间频率基准重点实验室，陕西 西安 710600 中国科学院 国家授时中心 时间频率基准重点实验室，陕西 西安 710600 Key Laboratory of Time and Frequency Primary Standards， National Time Service Center， Chinese Academy of Sciences， Xi’an 710600， China Key Laboratory of Time and Frequency Primary Standards， National Time Service Center， Chinese Academy of Sciences， Xi’an 710600， China 中国科学院大学 天文与空间科学学院，北京 100049 中国科学院大学 天文与空间科学学院，北京 100049 College of Astronomy and Spatial Sciences， University of Chinese Academy of Sciences， Beijing 100049， China College of Astronomy and Spatial Sciences， University of Chinese Academy of Sciences， Beijing 100049， China 西北大学 光子学与光子技术研究所，陕西 西安710069 西北大学 光子学与光子技术研究所，陕西 西安710069 Institute of Photonics and Photon Technology， Northwest University， Xi’an 710069， China Institute of Photonics and Photon Technology， Northwest University， Xi’an 710069， China 空间锶原子光钟在轨运行过程中会感受到幅度约为80 μT的地磁场变化，通过引入同等大小的外磁场变化，实验观测了变化地磁场对空间光钟冷原子由磁光阱装载到光晶格的影响。理论分析了地磁场变化对钟跃迁探测的影响，并对磁场稳定性需求进行了分析。通过电动旋转台和永磁体模拟空间环境变化的磁场，利用三维磁通门进行实时探测，设计了主动补偿系统并对外界磁场进行实时补偿，使空间光钟磁光阱中心区域磁场波动在亚μT量级，满足了空间光钟对磁场的补偿与控制需求。 The strontium atomic optical clock in space experiences changes in the geomagnetic field during its orbit of amplitudes of approximately 80 μT. By introducing changes in the external magnetic field of the same magnitude， the influence of the changing magnetic field on the loading of cold atoms from the magneto-optical trap to the optical lattice of the space optical clock is observed. The influence of the changing geomagnetic field on the detection of clock transitions is theoretically analyzed， and the need for magnetic field stability is investigated. The magnetic field of the space environment is simulated by rotating a permanent magnet. Furthermore， a three-dimensional fluxgate probe is used for real-time detection， and an active compensation system is designed to compensate the external magnetic field synchronously. Because the magnetic field in the central area of the magnetic optical trap of the space optical clock fluctuates by a sub-μT magnitude， it can meet the compensation and control requirements of the space optical clock with respect to the magnetic field. 空间光钟 地磁补偿 锶原子光钟 时间频率 space optical clock geomagnetic compensation strontium optical lattice clock time frequency 1 引　言 1 光钟的稳定度和不确定度可达10 -18 量级，优于微波钟两个量级，有望重新定义秒 ［ 1 1 - 2 2 1-2 ］ 。实验天体物理联合研究所（JILA）的锶原子光钟实现了2×10 - 18 的不确定度，是目前最准确的原子光钟 ［ 3 3 - 4 4 3-4 ］ 。光学原子钟作为当前时间（频率）测量能力最为强大的科学与技术研究平台，有望对基础科学前沿研究产生重大影响 ［ 5 5 - 9 9 5-9 ］ 。 在空间环境中可实现较地面精度更高的原子钟，实现稳定度和不确定度在10 -18 量级甚至更高的时间频率基准。作为超高精度的时间基准平台，空间锶原子光钟（以下简称空间光钟）可以将空间天文中对大空间、大时间尺度的观测转变为一定空间尺度下对时间频率的精确测量问题，在基础研究中，可以更为主动地检验基础物理规律，探索超越现有基本物理理论的新物理规律。应用超高精度时间频率信号来同步其他原子钟、重大科技设施和工程技术设施等，可以提高国际原子时精度、测定空间站轨道等。同时，空间光钟还能为其他空间探测技术，例如空间甚长基线干涉测量（VLBI）阵列、星际导航等，提供远高于目前的时间频率基准，极大提升这些技术的性能 ［ 10 10 ］ 。中国科学院国家授时中心承担了我国空间站锶原子光钟的研制工作，原理样机已实现了闭环运行，系统分时自比对频率稳定度为2.6×10 -15 /τ 1/2 ，积分时间为4 096 s时的环内稳定度为3.3×10 -17 ，目前处于工程化阶段 ［ 11 11 ］ 。 当空间光钟在近地轨道运行时，会感受到三维地磁场的复杂变化，单个维度的磁场变化大小约为80 μT ［ 12 12 ］ 。空间光钟的磁光阱（Magneto Optical Trap， MOT）和钟跃迁探测都对外磁场较为敏感，变化地磁场势必对它产生一定程度的影响。磁场补偿一般有两种方式，即利用磁屏蔽进行被动补偿或通过改变线圈电流进行主动补偿。被动补偿系统结构较为简单，无需复杂的电子设备，且无运转功耗，但是往往需要较厚的磁屏蔽材料才能达到理想的屏蔽效果。由于空间光钟原子探测区具有复杂而紧凑的光路和光电探测设备，若只对MOT腔进行磁屏蔽，则需要留有孔洞，难以保证磁屏蔽效果，并且MOT强磁场还会造成较大的磁滞，需要较大的封闭空间才能达到全方位磁屏蔽。这会带来过高的载荷，并影响空间的开放性。因此，单靠磁屏蔽往往难以保证高精度空间光钟的需求。目前，国内外的微波空间原子钟都是采用磁屏蔽和主动补偿相互配合的方法，而空间光钟目前还没有先例。本文分析了变化地磁场对空间光钟的影响，以及对磁场稳定性和可控性的需求。根据具体需求设计并研制了地磁场主动补偿系统，通过三维磁通门探头实时采集当前磁场大小，滤波后进行PID调节，再反馈至磁光阱地磁补偿线圈，对磁光阱内磁场进行实时补偿和控制，以满足空间光钟对磁场稳定性和可控性的要求。 2 空间光钟对磁场稳定性的需求分析 1 2.1　磁光阱系统构成 2 图1 图1 （a）为磁光阱结构及其光路示意图。磁光阱是锶原子光钟进行原子冷却与俘获的场所，空间锶原子光钟系统中实现磁光阱内对锶原子的俘获，需要通过反亥姆霍兹线圈为磁光阱提供均匀稳定的磁场梯度，并通过三个维度的冷却和俘获激光将原子俘获在磁场零点附近，即磁光阱的中心位置。 图1 图1 （b）为空间光钟物理真空系统模型，尺寸为100 cm×50 cm×60 cm，其中磁光阱由14个CF16法兰和1个CF35法兰构成，用于连接Zeeman减速器、复合真空泵、笼式结构光路和探测系统等结构。其中，笼式光路由3束对打的461 nm一级冷却光、689 nm二级俘获光、679 nm和707 nm重泵浦光，813 nm晶格光和698 nm钟跃迁探测光等组成。探测系统包括光电倍增管（Photo Multiplier Tube， PMT）和电子增强型电荷耦合元件（EMCCD）等设备。该物理系统是非常紧凑和复杂的。 1;2; 2.2　外磁场对系统的影响 2 外磁场的变化会对系统产生两方面的影响：一是影响冷原子的制备；二是影响钟跃迁探测。 外磁场的变化会引起MOT的磁场零点发生变化。二级冷却中MOT的磁场梯度为10 G/cm，因此假如不加任何措施，在地磁场的影响下MOT的中心位置会发生幅度约为0.8 mm的位置漂移，而锶光钟二级窄带冷却原子团直径一般仅为0.5 mm，这势必会对后续的光晶格装载造成较大影响。因此，实验模拟测量了在与空间光钟磁场变化相当的情形下，光晶格的冷原子装载情况。实验通过永磁体的转动产生相应的磁场变化，设置磁通门距离永磁体42.5 cm，测量不同方向的磁场大小，结果如 图2 图2 和 表1 表1 所示， XYZ 的方向坐标已在图中给出。 位　置 <math id="M1"><mi>X</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777821&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777812&type= 2.20133328 2.37066650 <math id="M2"><mi>Y</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777825&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777823&type= 2.53999996 2.37066650 <math id="M3"><mi>Z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777827&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777815&type= 2.03200006 2.37066650 <math id="M4"><mroot><mrow><msup><mrow><mi>X</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>Y</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>Z</mi></mrow><mrow><mn>2</mn></mrow></msup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777833&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777828&type= 21.59000015 3.80999994 图2 图2 （a） 0 -38 -1 39 图2 图2 （b） -10 -26 80 85 图2 图2 （c） -6 -31 32 45 图2 图2 （d） 17.5 -30 27 44 图2 图2 （e） -19 -19 38 47 磁场变化幅度 37 19 81 46 μT 由 图2 图2 和 表1 表1 可知，磁场在 Z 方向的变化幅度约为81 μT，在 X 方向的变化幅度约为37 μT，在 Y 方向的变化幅度约为19 μT。可见在此距离下，永磁体磁场的变化幅度与空间光钟在轨运行时的地磁场变化相当。 将永磁体与锶原子光钟系统MOT放置在同一水平线，在距离MOT中心42.5 cm的位置转动永磁体方向，通过EMCCD观察光晶格图像，结果如 图3 图3 所示。实验表明，81 μT的磁场强度会严重影响光晶格的装载效率，空间光钟在轨运行过程中，部分位置无法实现装载，导致光钟系统实验失败或数据报错。为保证磁场零点位置偏移不超过0.1 mm，磁场波动幅值要小于10 μT。 目前，锶原子光钟的钟跃迁大多通过自旋极化谱线来探测 ［ 13 13 ］ ，即探测磁量子数 m F =±9/2所对应的两个跃迁线，如 图4 图4 （a）所示。该方法需要外加一个磁场，使磁子能级发生一定程度的塞曼分裂。在外加磁场恒定不变的情况下，由于两个跃迁线的磁量子数等大但符号相反，其产生的塞曼频移也是等大反向，通过求取平均值便可消除一阶塞曼频移的影响。而二阶塞曼频移很小，可以忽略不计。因此，可将钟激光频率锁定到自旋极化谱双峰的中心频率 <math id="M5"><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777861&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777857&type= 2.62466669 3.72533321 上 ［ 14 14 ］ 。然而，要获取双峰的中心频率，需要对钟跃迁进行4次探测，即每4个运行周期完成一次对钟激光频率的伺服反馈，每个周期运行的时间约为1 s。当外加磁场发生变化时，上述过程会受到影响，如 图4 图4 （b）所示。为简单起见，假设外加磁场发生线性漂移，相应的钟跃迁谱线的中心频率也会发生频移。在单个钟运行周期中，钟跃迁线中心频率发生的塞曼频移变化量为 <math id="M6"><mo>∆</mo><msub><mrow><mi>ν</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777866&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777863&type= 8.46666718 4.06400013 ，则4个运行周期之后获取到的双峰中心频率的变化量为： <math display="block" id="M7"><mo>∆</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">4</mn><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub><mo>=</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>∆</mo><msub><mrow><mi>ν</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777843&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777869&type= 42.07933044 7.87400007 ， 其中： <math id="M8"><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777878&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777846&type= 2.62466669 3.72533321 为外加磁场恒定时的双峰中心频率， <math id="M9"><mi>ν</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777882&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777880&type= 1.43933344 2.79399991 1 为磁场发生线性漂移时的双峰中心频率。 <math id="M10"><mo>∆</mo><msub><mrow><mi>ν</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777885&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777876&type= 8.46666718 4.06400013 为一阶塞曼频移与二阶塞曼频移之和，即： <math display="block" id="M11"><mo>∆</mo><msub><mrow><mi>ν</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub><mo>=</mo><msubsup><mrow><mo>∆</mo></mrow><mrow><mi>B</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup><mo>+</mo><msubsup><mrow><mo>∆</mo></mrow><mrow><mi>B</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777896&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777887&type= 29.12533188 4.57200003 ， 其中： <math display="block" id="M12"><msubsup><mrow><mo>∆</mo></mrow><mrow><mi>B</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn></mrow></mfenced></mrow></msubsup><mo>=</mo><mo>-</mo><mi>α</mi><mo>·</mo><msub><mrow><mi>m</mi></mrow><mrow><mi mathvariant="normal">F</mi></mrow></msub><mo>·</mo><mo>∆</mo><msub><mrow><mi>B</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777891&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777889&type= 31.49600029 4.57200003 ， <math display="block" id="M13"><msubsup><mrow><mo>∆</mo></mrow><mrow><mi>B</mi></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn></mrow></mfenced></mrow></msubsup><mo>=</mo><mo>-</mo><mi>β</mi><mo>·</mo><mo>∆</mo><msup><mrow><msub><mrow><mi>B</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777901&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777899&type= 27.77066803 5.24933338 . 对于 87 Sr， <math id="M14"><mi>α</mi><mo>=</mo><mn mathvariant="normal">1.09</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777905&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777903&type= 13.46199989 2.79399991 Hz/μT ， <math id="M15"><mi>β</mi><mo>=</mo><mn mathvariant="normal">2.33</mn><mo>×</mo><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">5</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777924&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777912&type= 24.89200020 3.89466691 Hz/μT 2［ 15 15 ］ ， <math id="M16"><mo>∆</mo><msub><mrow><mi>B</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777929&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777926&type= 9.56733322 4.06400013 为单个运行周期中所发生的磁感应强度变化，可表示为磁感应强度变化率与单个钟运行周期 <math id="M17"><msub><mrow><mi>τ</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777915&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777932&type= 5.75733376 3.72533321 之积： <math display="block" id="M18"><mo>∆</mo><msub><mrow><mi>B</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi mathvariant="normal">d</mi><mi>B</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi>t</mi></mrow></mfrac><mo>·</mo><msub><mrow><mi>τ</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777939&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777936&type= 26.67000008 7.95866632 . 由此可以得出，基于 87 Sr光钟由磁场线性漂移导致的频率不确定度为4个钟跃迁探询周期间所发生的极化谱双峰中心频率塞曼频移与无干扰下的极化谱双峰中心频率的比值，即： <math display="block" id="M19"><mfrac><mrow><mo>∆</mo><mi>ν</mi></mrow><mrow><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="normal">Δ</mi><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">4</mn><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></mrow><mrow><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mfrac><mrow><mi mathvariant="normal">Δ</mi><msub><mrow><mi>ν</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub></mrow><mrow><msub><mrow><mi>ν</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777944&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777942&type= 37.93066406 8.80533314 . 如 图5 图5 所示，假设空间光钟的目标稳定度和频率不确定度为5×10 -18 ，则钟相邻询问周期间的磁场变化约小于0.3 nT。若 <math id="M20"><msub><mrow><mi>τ</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">y</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">e</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777920&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777919&type= 12.53066730 3.72533321 s，磁场的变化率要小于0.3 nT/s。地磁场在45 min内变化约80 μT ［ 16 16 ］ ，磁场的平均变化率约为30 nT/s，而在越过南太平洋磁场反常区时，变化率远大于该值。假如不加任何磁屏蔽和主动补偿，当外磁场出现与量子化轴方向垂直的分量时，会同时激发 <math id="M21"><mi mathvariant="normal">π</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777922&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777950&type= 1.94733346 2.79399991 跃迁和 <math id="M22"><mi mathvariant="normal">σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777968&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777953&type= 1.69333339 2.79399991 跃迁，从而使谱线的形状发生复杂变化。 1;2; 综上，要保证冷原子制备不受影响，磁场的波动幅值要小于10 μT，而对于一个频率稳定度和不确定度优于5×10 -18 的空间锶原子光钟，磁场的线性漂移率要小于0.3 nT/s，这意味着要将地磁场的波动大小补偿至远小于0.8 μT。 3 磁场主动补偿系统 1 由于空间光钟物理系统紧凑而复杂，这就要求磁屏蔽结构需要对整个物理系统进行封闭包裹，面积约为2.8 m 2 ，因此会带来较大的空间载荷。磁屏蔽的另外一个缺点是具有磁滞性，当外界地磁场和MOT的磁场周期变化时都会引入相应的剩磁。因此，研制主动补偿系统对空间锶原子光钟具有重要意义。另一方面，从前面的物理系统分析可知，空间锶原子光钟本身需要3对磁场补偿线圈，以产生钟跃迁探测时需要的偏置磁场，并补偿塞曼减速器等器件带来的杂散磁场，这为主动补偿地磁场提供了天然的便利条件。 图6 图6 为自主设计的磁场实时补偿系统示意图。考虑到锶光钟的补偿线圈内为真空系统的MOT区，其内部具有MOT磁场，且为原子俘获区，无法对其内部磁场进行实时测量，同时地磁场在空间锶原子光钟空间尺度内被认为是匀强磁场，因此设计两个孪生三维磁场补偿线圈系统。其中，一个线圈系统中心放置一个高精度三维磁通门，另一个线圈系统内则为MOT区的磁场补偿线圈系统。只要将两个线圈系统的方向保持一致，则可通过实时补偿和控制其中的一个补偿线圈系统中的磁场，来稳定和控制另外一个线圈内的磁场。实验上，对其中一个线圈系统内磁场进行了实时补偿。为产生变化磁场，将一个永磁体放置在一个转速可调的电动旋转台上，并将其与磁场补偿线圈保持一定距离。三维磁通门的三路采集电压信号分别代表 X ， Y ， Z 3个方向的磁场大小，将这三路信号分别连接数据采集卡的三路模拟输入通道进行数据采集，低通滤波后进行PID运算调节和反馈，反馈电压通过模拟输出通道输出到压控恒流源模块，再反馈至磁场补偿线圈，以实现MOT磁场的实时补偿。其中，三维磁通门采用的是Bartington公司的Mag-13U100，带宽 > 3 kHz，测量范围为±100 μT，噪声水平 < 5 pTrms/ <math id="M23"><mroot><mrow><mi mathvariant="normal">H</mi><mi mathvariant="normal">z</mi></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777992&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778000&type= 7.70466709 4.06400013 @1 Hz。 设置永磁体的转动周期为120 s，对磁场进行实时补偿的测量结果如 图7 图7 所示。图中分别显示了在 X ， Y ， Z 3个方向主动补偿前后的磁场大小。自由运转时，存在较大的电路噪声，磁场波动在3个坐标轴方向均约为30 μT，在反馈环路中加入了低通滤波（截止频率为5 Hz）；当闭环锁定时，3个方向的磁场能分别稳定在设定值上（根据谱线探测需求设置 X ， Y 方向的磁场为0， Z 方向为固定值），锁定后3个方向的磁场波动均小于0.5 μT。通过优化PID参数，其稳定度还可进一步提升。 4 结　论 1 本文分析了磁场波动对空间光钟可能造成的影响，模拟测量表明，在轨运行时的变化地磁场在不加任何补偿措施的情况下，由于MOT磁场零点的移动，能对冷原子的光晶格装载造成严重干扰，甚至可能导致冷原子光晶格装载的失败；在假设闭环时，钟频率锁定于自旋极化谱双峰中心频率的情况下，分析了磁场波动对钟频率不确定度的影响；明确了空间锶原子光钟对磁场稳定性的需求。然后，提出了利用孪生线圈补偿电流对MOT区磁场进行主动补偿的方法，并搭建了实时补偿系统，实现了剩磁的主动补偿，达到的磁场稳定性能满足频率稳定度和不确定度在10 -18 量级的空间锶原子光钟的需求。 本文为空间锶原子光钟的磁场补偿系统提出了一种解决方案，建立了原理样机，后续还需进行集成化电路设计以实现在空间光钟系统中的应用，并在真实物理系统中进行测试。 本文的磁场稳定与控制系统，在与磁场相关的冷原子物理研究中也具有潜在的应用价值，如散射长度时间相关玻色-爱因斯坦凝聚中的亮孤子行为 ［ 17 17 - 18 18 17-18 ］ 、分数化涡和涡晶格的动态产生 ［ 19 19 ］ 、弱连接微腔中的光子约瑟夫森效应 ［ 20 20 ］ 等研究。 参考文献 GREBING C ， AL-MASOUDI A ， DÖRSCHER S ， et al . Realization of a timescale with an accurate optical lattice clock ［J］. Optica ， 2016 ， 3 （ 6 ）： 563 - 569 . doi: 10.1364/optica.3.000563 http://dx.doi.org/10.1364/optica.3.000563 Realization of a timescale with an accurate optical lattice clock LE TARGAT R ， LORINI L ， LE COQ Y ， et al . Experimental realization of an optical second with strontium lattice clocks ［J］. Nature Communications ， 2013 ， 4 ： 2109 . doi: 10.1038/ncomms3109 http://dx.doi.org/10.1038/ncomms3109 Experimental realization of an optical second with strontium lattice clocks OELKER E ， HUTSON R B ， KENNEDY C J ， et al . Demonstration of 4.8×10 -17 stability at 1 s for two independent optical clocks ［J］. Nature Photonics ， 2019 ， 13 （ 10 ）： 714 - 719 . doi: 10.1038/s41566-019-0493-4 http://dx.doi.org/10.1038/s41566-019-0493-4 Demonstration of 4.8×10-17 stability at 1 s for two independent optical clocks BOTHWELL T ， KEDAR D ， OELKER E ， et al . JILA SrI optical lattice clock with uncertainty of 2.0×10 -18 ［J］. Metrologia ， 2019 ， 56 （ 6 ）： 065004 . doi: 10.1088/1681-7575/ab4089 http://dx.doi.org/10.1088/1681-7575/ab4089 JILA SrI optical lattice clock with uncertainty of 2.0×10-18 BLATT S ， LUDLOW A D ， CAMPBELL G K ， et al . Optical lattice clocks ［J］. Physical Review Letters ， 2008 ， 100 （ 14 ）： 140801 . doi: 10.1103/physrevlett.100.140801 http://dx.doi.org/10.1103/physrevlett.100.140801 Optical lattice clocks GODUN R M ， NISBET-JONES P B R ， JONES J M ， et al . Frequency ratio of two optical clock transitions in 171Yb + and constraints on the time variation of fundamental constants ［J］. Physical Review Letters ， 2014 ， 113 （ 21 ）： 210801 . doi: 10.1103/physrevlett.113.210801 http://dx.doi.org/10.1103/physrevlett.113.210801 Frequency ratio of two optical clock transitions in 171Yb+ and constraints on the time variation of fundamental constants HUNTEMANN N ， LIPPHARDT B ， TAMM C ， et al . Improved limit on a temporal variation of mp/me from comparisons of Yb + and Cs atomic clocks ［J］. Physical Review Letters ， 2014 ， 113 （ 21 ）： 210802 . doi: 10.1103/physrevlett.113.210802 http://dx.doi.org/10.1103/physrevlett.113.210802 Improved limit on a temporal variation of mp/me from comparisons of Yb+ and Cs atomic clocks KOLKOWITZ S ， PIKOVSKI I ， LANGELLIER N ， et al . Gravitational wave detection with optical lattice atomic clocks ［J］. Physical Review D ， 2016 ， 94 （ 12 ）： 124043 . doi: 10.1103/physrevd.94.124043 http://dx.doi.org/10.1103/physrevd.94.124043 Gravitational wave detection with optical lattice atomic clocks DEREVIANKO A ， POSPELOV M . Hunting for topological dark matter with atomic clocks ［J］. Nature Physics ， 2014 ， 10 （ 12 ）： 933 - 936 . doi: 10.1038/nphys3137 http://dx.doi.org/10.1038/nphys3137 Hunting for topological dark matter with atomic clocks RIEHLE F . Frequency Standards ： Basics and Applications ［M］. Germany ： WILEY-VCH ， 2004 ： 387 - 415 . Frequency Standards： Basics and Applications GUO F ， TAN W ， ZHOU C H ， et al . A proof-of-concept model of compact and high-performance 87Sr optical lattice clock for space ［J］. AIP Advances ， 2021 ， 11 （ 12 ）： 125116 . doi: 10.1063/5.0064087 http://dx.doi.org/10.1063/5.0064087 A proof-of-concept model of compact and high-performance 87Sr optical lattice clock for space LI L ， JI J W ， REN W ， et al . Automatic compensation of magnetic field for a rubidium space cold atom clock ［J］. Chinese Physics B ， 2016 ， 25 （ 7 ）： 073201 . doi: 10.1088/1674-1056/25/7/073201 http://dx.doi.org/10.1088/1674-1056/25/7/073201 Automatic compensation of magnetic field for a rubidium space cold atom clock 孔德欢 . 可搬运锶光钟性能评估及空间站锶光钟预研究 ［D］. 西安 ： 中国科学院大学（中国科学院国家授时中心） ， 2021 . doi: 10.7498/aps.70.20201204 http://dx.doi.org/10.7498/aps.70.20201204 可搬运锶光钟性能评估及空间站锶光钟预研究 KONG D H . Evaluation of the Transportable Strontium Optical Clock and Pre-research of the Optical Clock for Space Station ［D］. Xi'an ： National Time Service Center， Chinese Academy of Sciences ， 2021 . （in Chinese） . doi: 10.7498/aps.70.20201204 http://dx.doi.org/10.7498/aps.70.20201204 Evaluation of the Transportable Strontium Optical Clock and Pre-research of the Optical Clock for Space Station 任洁 ， 卢晓同 ， 王叶兵 ， 等 . 锶原子光钟闭环控制系统的设计与实现 ［J］. 光学 精密工程 ， 2018 ， 26 （ 10 ）： 2546 - 2554 . doi: 10.3788/OPE.20182610.2546 http://dx.doi.org/10.3788/OPE.20182610.2546 锶原子光钟闭环控制系统的设计与实现 REN J ， LU X T ， WANG Y B ， et al . Design and realization of an auto-control system for the closed-loop operation of a strontium atomic clock ［J］. Opt. Precision Eng. ， 2018 ， 26 （ 10 ）： 2546 - 2554 . （in Chinese） . doi: 10.3788/OPE.20182610.2546 http://dx.doi.org/10.3788/OPE.20182610.2546 Design and realization of an auto-control system for the closed-loop operation of a strontium atomic clock BOYD M . M . High Precision Spectroscopy of Strontium in an Optical Lattice： Towards a New Standard for Frequency and Time［D］. B.S. ： University of Washington ， 2002 . M LIU L ， LÜ D S ， CHEN W B ， et al . In-orbit operation of an atomic clock based on laser-cooled 87Rb atoms ［J］. Nature Communications ， 2018 ， 9 ： 2760 . doi: 10.1038/s41467-018-05219-z http://dx.doi.org/10.1038/s41467-018-05219-z In-orbit operation of an atomic clock based on laser-cooled 87Rb atoms LIANG Z X ， ZHANG Z D ， LIU W M . Dynamics of a bright soliton in Bose-Einstein condensates with time-dependent atomic scattering length in an expulsive parabolic potential ［J］. Physical Review Letters ， 2005 ， 94 （ 5 ）： 050402 . doi: 10.1103/physrevlett.94.050402 http://dx.doi.org/10.1103/physrevlett.94.050402 Dynamics of a bright soliton in Bose-Einstein condensates with time-dependent atomic scattering length in an expulsive parabolic potential CHIN C ， GRIMM R ， JULIENNE P ， et al . Feshbach resonances in ultracold gases ［J］. Reviews of Modern Physics ， 2010 ， 82 （ 2 ）： 1225 - 1286 . doi: 10.1103/revmodphys.82.1225 http://dx.doi.org/10.1103/revmodphys.82.1225 Feshbach resonances in ultracold gases JI A C ， LIU W M ， SONG J L ， et al . Dynamical creation of fractionalized vortices and vortex lattices ［J］. Physical Review Letters ， 2008 ， 101 ： 010402 . doi: 10.1103/physrevlett.101.010402 http://dx.doi.org/10.1103/physrevlett.101.010402 Dynamical creation of fractionalized vortices and vortex lattices JI A C ， SUN Q ， XIE X C ， et al . Josephson effect for photons in two weakly linked microcavities ［J］. Physical Review Letters ， 2009 ， 102 （ 2 ）： 023602 . doi: 10.1103/physrevlett.102.023602 http://dx.doi.org/10.1103/physrevlett.102.023602 Josephson effect for photons in two weakly linked microcavities 10.37188/OPE.20223011.1337 0433.1 TP399 A 1004-924X（2022）11-1337-07 任洁,谭巍,郭峰等.空间锶原子光钟磁场分析及主动补偿系统[J].光学精密工程,2022,30(11):1337-1343. 任洁,谭巍,郭峰等.空间锶原子光钟磁场分析及主动补偿系统[J].光学精密工程,2022,30(11):1337-1343. REN Jie,TAN Wei,GUO Feng,et al.Magnetic field analysis and active compensation system for strontium optical lattice clock in space[J].Optics and Precision Engineering,2022,30(11):1337-1343. REN Jie,TAN Wei,GUO Feng,et al.Magnetic field analysis and active compensation system for strontium optical lattice clock in space[J].Optics and Precision Engineering,2022,30(11):1337-1343. 国家自然科学基金资助项目（No.11903042，No.61775220）；中国科学院前沿科学重点研究项目（No. QYZBD-SSW-JSC004）；中国科学院战略性先导科技专项项目（No.XDB35010202） 2021-12-22 2022-01-14 2022-06-10 2022-06-21 光学精密工程 11 30