赵芳婧 高峰 韩建新 周驰华 孟俊伟王叶兵 郭阳 张首刚 常宏

1)(中国科学院国家授时中心,时间频率基准重点实验室,西安 710600)

2)(中国科学院大学天文与空间科学学院,北京 100049)

(2017年12月4日收到;2017年12月21日收到修改稿)

1 引 言

冷原子锶光钟作为目前世界上性能指标最高的原子钟,其测量精度已进入10-19量级[1],比目前基准钟铯微波喷泉钟高出三个量级,是下一代时间频率标准“秒”定义的有力候选.近十几年来,国际上很多实验室都在积极地开展锶光钟的研制[2−8].

随着时间频率测量精度的不断提高,高性能原子光钟的应用越来越广泛.为了更好地开展基于光钟的高精密测量实验[9],比如里德堡常数RH测量、引力红移测量[10]、精细结构常数随时变化,以及精密探测相对论和基础理论所预言的一些“微弱现象”,并检验它们的正确性[11−13],科学家们一直在探索新的研究方向.

目前,国际上多个研究组都竞相开展空间高精度原子钟的研制,美国和欧洲均启动了利用国际空间站环境开发空间冷原子钟的研究计划[14−16].国内,中国科学院上海光学精密机械研究所从2000年开始小型化冷原子铷钟的研制,并在2016年首次成功实现了空间冷原子铷钟,地面测试中期千秒稳定度为1.5×10−14[17],已经成功搭载天宫二号进入太空,并且正常运行.最早研制空间冷原子钟的是欧洲空间局提出的空间原子钟组(Atomic Clock Ensemble in Space,ACES)计划,该计划的目的是在空间站研制冷原子微波钟.在2006年欧洲空间局开始空间光钟研制计划,目标是在2023年左右实现频率不确定度优于1×10−17、稳定度优于的空间冷原子光钟,该计划目前仍处在原理样机的地面测试阶段[18].受空间站特殊条件及光钟系统自身体积庞大和结构复杂等因素的限制,在空间光钟研制过程中,首要解决的问题就是实现光钟的小型化及低功耗设计.

近年来,关于光钟小型化的研究工作一直在进行中[19−22].2014年,意大利Tino小组对整个光钟系统进行了优化,将物理真空系统的体积缩小至120 cm×40 cm×36 cm,且系统总体积＜2 m3,最终实现了频率不确定度为7.0×10−15的可搬运88Sr光钟[19].2017年,德国联邦物理技术研究院成功实现了车载可移动锶光钟,光钟物理真空系统固定在120 cm×90 cm的平板上,最终实现了87Sr光钟,频率不确定度为7.4×10−17,稳定度优于是目前报道的不确定度和稳定度最高的可移动光钟[20].此外,中国武汉物理与数学研究所黄学人小组研制的小型化钙离子光钟目前也实现了可搬运.他们研制的小型化原子钟系统总体积＜0.54 m3,初步评估得到光钟不确定度为7.8×10−17,稳定度约为

本文主要介绍小型化锶光钟物理系统的研制实验,将构建磁光阱(magneto-optical trap,MOT)的一对反亥姆霍兹线圈置于真空腔内.线圈采用Kapton材质的真空铜线绕制,实验中线圈通电电流为2 A,在MOT中心区域产生轴向磁场梯度为43 Gs/cm,发热功率小于1.0 W.整个线圈固定在一个易于导热的铝制支架上,支架和真空腔体的紧密贴合可实现很好的散热功效,无需增加循环水冷装置以实现散热和控温.目前,实验已经成功将物理真空系统体积缩小至60 cm×20 cm×15 cm,约为实验室原锶光钟系统体积的1/10[23],并且采用真空内置磁场线圈的技术,已经成功实现了锶原子的一级多普勒冷却.实验测得俘获区锶冷原子团的直径为1.5 mm,利用短程飞行时间法[24−26]测得冷原子温度约为10.6 mK.利用荧光收集法[27]测得锶同位素88Sr和87Sr的冷原子数目分别为1.6×106和1.5×105.在一级多普勒冷却过程中,对应的原子跃迁能级结构的不封闭性造成了原子在3P2和3P0两能态上的堆积,从而影响了冷原子的俘获效率.加入707和679 nm重抽运激光,可极大地降低原子损失率,实验最终可将冷原子数目提高5倍以上.

2 实验系统设计

2.1 真空物理系统

锶光钟整个物理系统装置主要包括原子Dispenser源、毛细准直器、Zeeman减速器、磁光阱、离子泵及真空组件等.如图1所示,装置最右端为锶原子Dispenser源,内装约5 g固态锶,在其前端放置长15 mm、直径0.8 mm的50根毛细管堆叠而成的毛细准直器.实验中,通电后的Dispenser源释放出锶原子蒸气,经由毛细准直器准直形成热原子束喷出.装置中间部分是长29.6 cm的Zeeman减速器,经过Zeeman减速器后热原子束被减速至100 m/s以下,减速后的原子被俘获到MOT中.在锶原子Dispenser源附近接有一个10 L的离子泵,实验中其真空度约为10−5Pa,此外在MOT区还接有一个20 L的离子泵,实验中真空度约为10−6Pa.装置最前端的正八面体真空腔是最终俘获冷原子团的场所,俘获光、减速光及探测光分别从真空腔体不同的窗口导入腔中并与原子相互作用,实现一级多普勒冷却.真空腔体水平方向所有窗口尺寸均为CF16.

图1 物理系统装置示意图Fig.1.Sketch of the vacuum system for a low maintenance Sr optical lattice clock.

实验中Zeeman减速器是通过补偿原子多普勒频移来实现对原子束的持续减速,目前常用的减速器有通电线圈型和永磁体型两种.通电线圈型制作简单易调节,且产生的磁场曲线平滑,后期优化方便,缺点是需要电流源控制及循环水冷装置,系统结构复杂且体积较大.永磁体型优点是无需电流源控制及水冷装置,体积小,无功耗,但是对磁体材料及安装位置要求高,环形永磁体不易精密加工,且充磁后磁体的中心剩余磁场强度不能调节.基于对以上诸多因素及现有实验条件的综合考量,采用多匝线圈式Zeeman减速器模拟永磁体环型Zeeman减速器,通过分析实验结果选择合适的参数为下一步永磁体型Zeeman减速器的制作提供参考.减速器全长29.6 cm,共12个线圈.线圈间采用厚度为2 mm的无磁不锈钢板隔开,单个线圈均采用直径1.5 mm的铜丝绕制而成,共计15圈18层270匝.12个线圈所加电流各不相同,方便可调谐,能实现最佳减速效果.

2.2 真空内置磁场线圈设计

与单一光学阱、磁阱相比,磁光阱具有阱深较大,阱中原子温度低、囚禁时间长,冷却效果明显等特点.为了实现三维俘获,在磁光阱xyz三个方向上分别用三对光斑尺寸相同、光强相等的激光两两对射,同时调节冷却激光的偏振,使得激光偏振分别为σ+和σ−[28].此外,还需要一对反向亥姆霍兹线圈即可构建磁光阱.

通常情况下构成磁光阱的反亥姆霍兹线圈被放在MOT区真空腔之外,该结构虽然便于调节但热耗功率很大.实验室原锶光钟系统中,线圈采用1.5 mm的铜线绕220匝,线圈内径约为110 mm,外径为130 mm.要达到实验所需磁场梯度50 Gs/cm,线圈所需电流约为16.5 A,发热功率为150 W.通电一段时间后发热严重.过高的温度不仅对线圈不利,还会影响MOT区真空度以及实验系统的稳定性.因此,必须增加循环水冷装置实现散热和控温.这使得整个物理系统体积非常庞大,阻碍了物理系统向小型化方向发展,不利于未来可搬运光钟及空间光钟的实现.

为实现物理系统的小型化目标,实验将一对反亥姆霍兹线圈放置到真空腔内部,使得整个系统体积大大缩小.结合实验要求,对两线圈的大小及线圈之间的距离进行精确的理论计算,并模拟线圈产生的磁场.根据真空腔的体积、窗口大小以及线圈尺寸,设计加工了铝制线圈支架,可将线圈直接绕制在支架上,然后将其放入真空腔中.内置线圈装置截面图如图2(a)所示,整个真空腔长120 mm,高69 mm.线圈采用直径1 mm的铜线绕制而成,共计160匝.两线圈内径为20 mm,外径为40 mm,如图2(b)所示.

图2 内置反亥姆霍兹线圈 (a)装置截面图;(b)装置实物图Fig.2.Schematic of the built-in anti-Helmholtz coil:(a)Sectional view;(b)set-up picture.

用磁强计对MOT中心区域的轴向和径向磁场进行精细测量.数据拟合后得到轴向磁场梯度为43 Gs/cm,径向磁场梯度约为轴向的一半,为26 Gs/cm,满足实验要求.

2.3 实验主光路

图3为锶光钟一级多普勒冷却光路.实验所用激光光源型号为TA–SHG 110倍频半导体激光器(Toptica),输出为波长461 nm的线偏振光,激光线宽小于2 MHz,扫频范围约为20 GHz,输出功率为190 mW.用于一级冷却的光路主要包括三个部分:第一部分,使激光器端口输出的461 nm激光双次通过声光调制器(acoustic-optic-modulator,AOM)1,移频+570 MHz之后,将激光频率锁定到锶原子1S0→1P1跃迁线上;第二部分,锁频后,将461 nm激光分为两部分,其中一束激光直接作为减速光,与Zeeman减速器相结合,共同与原子作用,将原子从最可几速度430 m/s减速至75 m/s;第三部分,在MOT区域,冷却并俘获减速后的原子.使锁频后的另一束461 nm激光双次通过AOM2,移频+532 MHz后,分成三路.在磁光阱的xyz方向,三对偏振方向相反的圆偏振激光两两对射,形成六束激光与冷原子相互作用.同时,在MOT区,还有一对真空内置的反亥姆霍兹线圈,在中心区域产生43 Gs/cm的轴向磁场梯度,将原子俘获在MOT区中心.

图3 一级多普勒冷却光路图Fig.3.First stage Doppler cooling laser system.

3 实验结果及分析

3.1 原子速度分布曲线

Zeeman减速是冷原子制备过程中一个极其重要的技术环节.实验中Zeeman减速光采用一束功率为40 mW、失谐−570 MHz、光斑为会聚型的圆偏振461 nm激光,与相对应的Zeeman磁场共同作用使锶原子束持续减速.由于受真空腔体窗口尺寸的限制,实验采用的俘获光光斑尺寸约为13 mm,失谐−38 MHz,每束光功率均为9 mW,只有锶原子速度降至100 m/s以下时才能被MOT有效俘获.因此,先进行原子速度分布曲线探测,获得Zeeman减速效率之后再进一步优化实验参数.探测过程中,需要两台461 nm激光器,其中一台频率锁定后提供一束Zeeman减速光与原子相互作用.同时,将另一台461 nm激光器设置为扫描状态,并将其输出激光分为两部分,一束垂直原子束运动方向入射的激光作为零失谐标定;另一束激光则与原子运动方向成45°角入射作为探测光,得到原子速度分布谱线.独立调节Zeeman减速器各线圈的电流值,使得MOT中心区域俘获原子数目达到最大.

图4所示为经Zeeman减速器减速前后的原子速度分布曲线.黑色实线表示减速前,红色虚线表示减速后.从图中可以看出温度约为470°C的Dispenser喷出的锶原子束最可几速率可达到430 m/s.通过原子速度分布曲线可以看出,经过Zeeman减速后,在原子速度为220和75 m/s两处,原子数目有堆积存在,经过理论拟合可以计算出减速后速度为75 m/s的原子占原子总数的16.4%.

图4 原子的速度分布Fig.4.Atomic velocity distribution.

3.2 锶原子同位素荧光图像

如图5所示,实验中用型号为CoolSNAP EZ的高分辨率CCD观测磁光阱中俘获的冷原子团.锶原子有四种同位素,目前用于冷原子锶光钟研究的有玻色子88Sr和费米子87Sr,88Sr在自然界中的丰度为82.6%,而87Sr在自然界中的丰度仅为7.0%,因此实验中冷原子团荧光信号的强度相差较大.利用CCD采集冷原子图像,对其荧光亮度进行高斯拟合,计算得到88Sr冷原子团直径约为1.5 mm.

实验利用自制的光电探测器加前置放大器对冷原子荧光强度进行探测.通过光电探测器收集不同冷原子团的荧光功率,可以推算出俘获冷原子的数目,88Sr和87Sr的原子数目分别为1.6×106和1.5×105,两者比值约为10.7,符合其同位素的自然丰度比例.重抽运激光707和679 nm的加入,最终可将冷原子数目提高5倍以上.

图5 锶原子同位素的冷原子荧光图像 (a)88Sr;(b)87SrFig.5.Fluorescence images of cold atoms of strontium isotopes:(a)88Sr;(b)87Sr.

3.3 短程飞行时间法测量冷原子温度

温度是表征冷原子特征的基本物理量,实验中常采用飞行时间法(time of flight,TOF)[29]测量温度.这种测量方式下,冷原子飞行距离一般较大,无需考虑冷原子团初始位置的大小及探测光斑尺寸的影响,只要探测光斑可覆盖下落的冷原子即可.但实验中由于物理系统真空腔体光学探测窗口较小,允许原子飞行的时间及路程较短,冷原子团初始位置的大小及探测光斑尺寸不能被忽略,传统的TOF法不再适用.因此,实验采用短程飞行时间法(飞行距离＜10 mm)测量冷原子温度.实验中,冷原子团被囚禁在MOT中心,半径约为0.75 mm.在距离冷原子团下方3 mm处,打一束光斑半径为0.5 mm的细圆光束作为探测光.通过时序精确控制MOT区磁场及实验中各光场的开关断,使得冷原子团得到周期性的释放与俘获,最终由探测器收集到的吸收信号推算出冷原子团温度为10.6 mK,如图6所示.

图6 短程TOF法测量的冷原子吸收信号Fig.6.Absorption curve of cold atoms by TOF.

3.4 利用Zeeman减速器研究最小减速器长度

实验中Zeeman减速器通过补偿由原子多普勒效应引起的频移来实现减速光与原子能级持续共振,从而使原子束减速.原子在减速器轴向上的运动可近似看作加速度为a的匀加速运动,因此Zeeman减速器长度l可表示为

式中r0为俘获的冷原子团半径,tb,td,ta分别为原子束离开准直器到达减速器的时间、在减速器中的时间和从减速器到达MOT的时间.由此,可由轴向最大速度径向最大速度求出长度为l的减速器作用下的俘获原子数.对原子束运动速度分布[31]积分可得[30]

式中R为任意长度的Zeeman减速器作用下俘获的原子数目,为原子束最可几速度,m为锶原子质量,k为玻尔兹曼常数,T为锶原子温度.

基于光钟物理系统小型化实验研制的考虑,在保证获得较高减速效率的情况下,应尽量缩短Zeeman减速器的尺寸.实验通过逐渐减少Zeeman减速器的线圈数来模拟不同长度的减速器.通过调节各个线圈电流值,实现不同长度下的最优MOT.最终可得到不同长度的Zeeman减速器对应的线圈最优电流值.将多螺线管Zeeman减速器的一个线圈模拟为永磁体Zeeman减速器的一对磁片,利用这些电流值可以进行磁场拟合,从中选择最优的磁场拟合结果,为下一步设计永磁体Zeeman减速器的长度及磁片的选择提供参考.因为要实现光钟的可搬运,甚至未来的空间冷原子光钟,永磁体Zeeman减速器无疑是更好的选择.永磁体Zeeman减速器不仅拥有结构简单、体积小、质量轻、无需电源供电、无功耗等显著优点,其优良的抗震性及稳定性也是实现小型化以及未来可搬运光钟和空间光钟所必需的.

实验操作如下:保证其他实验条件保持不变,从靠近锶炉的一端开始,逐渐减少Zeeman减速器线圈数,重新调节各个线圈的电流值,同时观察不同状态下光电探测器探测到的冷原子团的荧光功率,使得俘获冷原子效果最优,并计算出相应原子数目.结果如图7所示.

图7 Zeeman减速器长度与冷原子数目的关系Fig.7.Number of cold atoms under different lengths of Zeeman slower.

由图7可以看出,理论与实际拟合较好,当Zeeman减速器长度为8—12 cm时,通过调节Zeeman减速器线圈电流,仍然能俘获106数量级的冷原子,俘获效率约11%.

4 结 论

采用真空内置反亥姆霍兹线圈构建磁光阱,实现了锶原子的一级多普勒冷却.实验中观测了不同的锶同位素的冷原子团,计算得到原子数目,并用短程飞行时间法测出冷原子团温度.实验验证了真空内置反亥姆霍兹线圈构成磁光阱这一关键技术,并成功实现了锶光钟物理系统的小型化.这对于未来可搬运光钟以及空间钟的实现都具有重要意义.在已有实验的基础上,我们可以设计永磁体Zeeman减速器.永磁体Zeeman减速器体积小、重量轻、无需电源供电且具有优良的防震性能及稳定性,利用这些优点可以将Zeeman减速器的尺寸做到更小,更好地实现锶原子光钟物理系统的小型化,为可移动光钟及空间光钟的研制奠定基础.

[1]Campbell S L,Hutson R B,Marti G E,Goban A,Darkwah Oppong N,McNally R L,Sonderhouse L,Robinson J M,Zhang W,Bloom B J,Ye J 2017Science358 90

[2]Le Targat R,Lorini L,Le Coq Y,Zawada M,Guéna J,Abgrall M,Gurov M,Rosenbusch P,Rovera D G,Nagórny B,Gartman R,Westergaard P G,Tobar M E,Lours M,Santarelli G,Clairon A,Bize S,Laurent P,Lemonde P,Lodewyck J 2013Nat.Commun.4 2109

[3]Takano T,Takamoto M,Ushijima I,Ohmae N,Akatsuka T,Yamaguchi A,Kuroishi Y,Munekane H,Miyahara B,Katori H 2016Nat.Photon.10 662

[4]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 2014New J.Phys.16 073023

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

[6]Lin Y G,Wang Q,Li Y,Meng F,Lin B K,Zang E J,Sun Z,Fang F,Li T C,Fang Z J 2015Chin.Phys.Lett.32 090601

[7]Hachisu H,Ido T 2015Jpn.J.Appl.Phys.54 112401

[8]Akamatsu D,Inaba H,Hosaka K,Yasuda M,Onae A,Suzuyama T,Amemiya M,Hong F L 2014Appl.Phys.Express7 012401

[9]Schiller S,Görlitz A,Nevsky A,Koelemeij J C J,Wicht A,Gill P,Klein H A,Margolis H S,Mileti G,Sterr U,Riehle F,Peik E,Tamm C,Ertmer W,Rasel E,Klein V,Salomon C,Tino G M,Lemonde P,Holzwarth R,Hänsch T W 2007Nucl.Phys.B166 300

[10]Salomon Ch,Dimarcq N,Abgrall M,Clairon A,Laurent P,Lemonde P,Santarelli G,Uhrich P,Bernier L G,Busca G,Jornod A,Thomann P,Samain E,Wolf P,Gonzalez F,Guillemot Ph,Leon S,Nouel F,Sirmain Ch,Feltham S 2001C.R.Phys.2 1313

[11]Cacciapuoti L,Salomon C 2009Eur.Phys.J.Special Topics172 57

[12]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 2014Phys.Rev.Lett.113 210801

[13]Fortier T M,Ashby N,Bergquist J C,Delaney M J,Diddams S A,Heavner T P,Hollberg L,Itano W M,Jefferts S R,Kim K,Levi F,Lorini L,Oskay W H,Parker T E,Shirley J,Stalnaker J E 2007Phys.Rev.Lett.98 070801

[14]Sullivan D B,Ashby N,Donley E A,Heavner T P,Hollberg L W,Jefferts S R,Klipstein W M,Phillips W D,Seidel D J 2005Adv.Space Res.36 107

[15]Schiller S,Görlitz A,Nevsky A,Alighanbari S,Vasilyev S,Abou-Jaoudeh C,Mura G,Franzen T,Sterr U,Falke S,Lisdat C,Rasel E,Kulosa A,Bize S,Lodewyck J,Tino G M,Poli N,Schioppo M,Bongs K,Singh Y,Gill P,Barwood G,Ovchinnikov Y,Stuhler J,Kaenders W,Braxmaier C,Holzwarth R,Donati A,Lecomte S,Calonico D,Levi F 2012Let’s Embrace Space(Vol.II)(Luxembourg:Publications Office of the European Union)p452

[16] Świerad D,Häfner S,Vogt S,Venon B,Holleville D,Bize S,Kulosa A,Bode S,Singh Y,Bongs K,Rasel E M,Lodewyck J,Le Targat R,Lisdat C,Sterr U 2016Nat.Sci.Rep.6 33973

[17]Li L,Qu Q Z,Wang B,Li T,Zhao J B,Ji J W,Ren W,Zhao X,Ye M F,Yao Y Y,Lü D S,Liu L 2016Chin.Phys.Lett.33 063201

[18]Origlia S,Schiller S,Pramod M S,Smith L,Singh Y,He W,Viswam S,Świerad D,Hughes J,Bongs K,Sterr U,Lisdat C,Vogt S,Bize S,Lodewyck J,Le Targa R,Holleville D,Venon B,Gill P,Barwood G,Hill I R,Ovchinnikov Y,Kulosa A,Ertmer W,Rasel E M,Stuhler J,Kaenders W,the SOC2 consortium contributors 2016Quantum Opt.9900 990003

[19]Poli N,Schioppo M,Vogt S,Falke St,Sterr U,Lisdat Ch,Tino G M 2014Appl.Phys.B117 1107

[20]Koller S B,Grotti J,Vogt S,Al-Masoudi A,Dörscher S,Häfner S,Sterr U,Lisdat C 2017Phys.Rev.Lett.118 073601

[21]Cao J,Zhang P,Shang J J,Cui K F,Yuan J B,Chao S J,Wang S M,Shu H L,Huang X R 2017Appl.Phys.B123 112

[22]Vogt S,Lisdat C,Legero T,Sterr U,Ernsting I,Nevsky A,Schiller S 2011Appl.Phys.B104 741

[23]Xu Q F,Liu H,Lu B Q,Wang Y B,Yin M J,Kong D H,Ren J,Tian X,Chang H 2015Chin.Opt.Lett.13 100201

[24]Geng T,Yan S B,Wang Y H,Yang H J,Zhang T C,Wang J M 2005Acta Phys.Sin.54 5104(in Chinese)[耿涛,闫树斌,王彦华,杨海菁,张天才,王军民 2005物理学报54 5104]

[25]Fu J X,Li Y M,Chen X Z,Yang D H,Wang Y Q 2001Acta Opt.Sin.21 414(in Chinese)[付军贤,李义民,陈徐宗,杨东海,王义遒2001光学学报21 414]

[26]Brzozowski T M,Maczynska M,Zawada M,Zachorowski J,Gawlik W 2002J.Opt.B4 62

[27]Tian X 2010M.S.Thesis(Xi’an:National Time Service Center,University of Chinese Academy of Sciences)(in Chinese)[田晓 2010硕士学位论文 (西安:中国科学院大学国家授时中心)]

[28]Wang Y Q 2007Laser Cooling and Trapping of Atoms(Beijing:Peking University Press)p294(in Chinese)[王义遒2007原子的激光冷却与陷俘(北京:北京大学出版社)第294页]

[29]Lett P D,Watts R N,Westbrook C I,Phillips W D 1988Phys.Rev.Lett.61 169

[30]Savard T A 1998Ph.D.Dissertation(Durham:Duke University)

[31]Ovchinnikov Y B 2008Eur.Phys.J.Special Topics163 95