“仪器与测量”栏目简介

（Instrumentation and Measurement）

先进的仪器设备和测量方法是推动科学进步的基石。为支持我国科研仪器设备的研发和测量方法的发展，推动物理学在科技创新中发挥更大的作用，《物理学报》开设“仪器与测量”栏目，重点关注物理学领域及其交叉学科有助于提升物理学研究水平的仪器研制与测量方法发展的创新性工作，推动科研仪器国产化，构建起前沿物理研究和工程应用之间的桥梁。本栏目接收原创性论文、研究快讯和综述，涵盖但不限于以下主题：

      1）新型仪器的设计与研发；

      2）现有仪器设备测量精度和应用范围的提升与拓展；

      3）重要仪器的自主知识产权技术方案；

      4）物理测量极限水平的探索；

      5）先进测量方法的改进、发展和应用；

      6）高通量物理测量新方法与应用；

      7）跨学科的测量问题、解决方案及测量技术的新应用；

      8）尖端仪器及测量方法国际最新动态。

基于⁷Li冷原子操控的超高真空测量

成永军，董猛^†，孙雯君，吴翔民，

张亚飞，贾文杰，冯村，张瑞芳

兰州空间技术物理研究所，真空技术与物理国家级重点实验室，兰州 73000

物理学报, 2024, 73(22): 220601

doi: 10.7498/aps.73.20241215

关键词

超高真空测量，冷原子，磁光阱，磁阱

01  引    言

7个国际单位制(SI)基本单位全部实现量子化定义促进真空测量体系向量子化转变。真空测量的研究对象是微观中性气体分子，真空参数作为导出量是最适宜进行量子测量的计量参数之一。冷原子真空测量技术属于量子计量的范畴，也是目前真空测量前瞻性和颠覆性技术研究方向之一。对照国际上对真空区域的划分，基于冷原子的量子物理特性适宜于超高真空(10^–5—10^–9 Pa)测量。超高真空测量在空间科学、真空表面分析、核聚变、高能加速器、真空微电子技术等方面具有广泛的应用价值^[1–3]。一直以来，超高真空测量主要采用电离真空计，使用前需要经过多级真空标准的校准传递才能实现量值溯源。另外，仪器自身存在的高温灯丝热辐射效应、吸放气效应、催化化学反应等固有特性会改变被测真空环境，影响测量准确性^[4–6]。而基于冷原子碰撞损失逃逸特性的真空测量，其本质是测量气体分子数密度，具有非接触、对真空环境扰动小等优点，其测量机理是基于冷原子-中性气体分子碰撞体系的本征物理特性，并可作为基础标准直接溯源至SI单位的基本量——秒和开尔文，不需要进行繁琐的校准。

利用冷原子效应测量真空的基本构想最初启发于囚禁于真空系统中冷原子的寿命由背景气体分子数密度决定，于1988年由Prentiss等^[7]首次构建理论雏形。发展至今，目前国际上最有代表性的研究机构主要有美国家标准与技术研究院(NIST)、加拿大英属哥伦比亚大学(UBC)和俄罗斯科学院。美国NIST从冷原子真空测量第一性原理理论计算和实验测量验证两个方面开展了广泛深入研究，已建立了实验室大型测量装置，占地面积约2 m²，测量下限达到了4×10^–9 Pa^[8–17]。加拿大UBC主要在磁阱中开展了⁸⁷Rb冷原子真空测量技术研究^[18–20]，俄罗斯科学院在光偶极阱中开展了6Li冷原子真空测量研究^[21]，也均已建立了实验测量验证装置。

在国内，兰州空间技术物理研究所紧随国际真空测量技术发展前沿，于2020年在国家重大项目的支持下，牵头开展了基于冷原子的超高真空测量技术研究。经过多年技术攻关，已研制出了我国首台⁷Li冷原子超高真空测量装置，如图1所示，整体技术达到国际先进水平。本文系统介绍了测量装置结构组成和工作原理，并全面开展了真空测量实验比对验证等研究工作，本项研究工作对推动我国全新跨代真空测量技术发展、创立真空量子测量标准、研制小型化量子真空测量传感器具有一定的支撑作用，并将来在空间科学探测、引力波测量、重粒子加速器等军民技术领域中具有较好的应用前景。

图1  ⁷Li冷原子真空测量装置实物图

02  测量理论模型

冷原子真空测量的基本原理是，囚禁于势阱中的冷原子与真空室的中性气体分子发生碰撞，当冷原子获得的能量大于势垒阱深时，原子发生逃逸损失，最后通过确定原子阱深、原子损失率和损失率系数等参数即可实现真空度的测量，基本理论公式为

式中，n 为分子数密度(m^–3)，k_B为玻尔兹曼常数(1.38×10^–23 J/K)，T 为气体分子温度(K)，Γ_loss为冷原子损失率(s^–1)，U 为原子阱深(K)，k_loss(U )为碰撞损失率系数(m³/s)。冷原子损失率Γ_loss通过实验测量精确获得，损失率系数k_loss(U )可采用理论计算和实验测量的方式获取，对于浅势阱，可采用ab initio第一性原理精确计算得到。由于k_loss(U )是与真空度无关的表征碰撞体系散射特性的物理常数，因此，真空度p这一宏观量与微观量k_loss紧密联系在一起，从而实现量子测量。

对于损失率系数k_loss的准确获取，在磁光阱(MOT)中，原子阱深通常在0.2—5 K范围内变化，冷原子-气体分子碰撞过程中存在量子衍射碰撞效应，即碰撞后冷原子获得的能量还不足以逃出囚禁区造成相空间原子数目的减少。因此，采用基于范德瓦耳斯长程相互作用系数C₆的半经典理论可计算损失率系数。碰撞体系的总损失率系数k_t^[9]为

式中，κ =12.88；μ =m_h⋅m_c/(m_h+m_c)为约化质量，m_h为气体分子质量；m_c为⁷Li原子质量(1.17×10^–26 kg)；T为气体分子温度(K)；为范德瓦耳斯能量；x₆=(2μC₆/ℏ²)^1/4为范德瓦耳斯长度；ℏ是约化普朗克常数(值为1.0546)。量子衍射损失率系数k_g的一阶近似为

式中，ς= 0.3755，U为阱深(K)。则⁷Li冷原子-中性气体碰撞损失率系数k_loss = k_t－k_g。图2为基于以上半经典理论计算的⁷Li-N₂，⁷Li-Ar，⁷Li-He，⁷Li-H₂碰撞体系在不同阱深条件下的损失率系数，阱深变化范围取1×10^–4—1 K，气体分子温度取293 K。

图2  碰撞体系的损失率系数k_loss与阱深U的关系曲线

结果表明，损失率系数k_loss随原子阱深的增大而逐渐减小，当阱深小于1 mK时，损失率系数为恒定值，这是由于1 mK以下的原子阱深为浅势阱，碰撞体系的量子衍射碰撞效应可以忽略。⁷Li-N₂，⁷Li-Ar，⁷Li-He和⁷Li-H₂的k_loss值分别为2.45×10^–15，2.15×10^–15，1.90×10^–15和3.93×10^–15 m³/s。可见浅势阱下k_loss与气体种类弱相关，均在10–15 m³/s 量级，除H₂外，最大偏差在30%以内，而对于电离真空计，各气体组分灵敏系数的偏差可达数倍^[22]，这也表明在真空腔中待测气体组分未知的情况下，基于冷原子的量子真空全压力测量的准确性更高。

03  实验测量装置

⁷Li冷原子的真空测量装置，其结构组成主要由⁷Li冷原子冷却囚禁系统和动态流量真空系统两部分组成。测量装置结构组成及工作原理如图3所示。基本工作原理为：原子炉加热产生的高速原子蒸气进入二维磁光阱(2D MOT)进行预冷却减速，预冷却原子通过差分管道进入三维磁光阱(3D MOT)，在磁场和光场的作用下被俘获而形成冷原子团。超高/极高真空系统通过微流量计引入气体分子，气体分子与⁷Li冷原子发生碰撞导致原子从势阱中逃逸，最后通过采用光电倍增管(PMT)荧光探测系统测量⁷Li冷原子逃逸损失特性实现真空度准确反演。

图3  ⁷Li冷原子真空测量装置原理图

3.1 ⁷Li冷原子冷却囚禁系统

⁷Li冷原子冷却囚禁系统用于在3D MOT中俘获真空测量传感⁷Li原子，主要由锂源高温原子炉、2D MOT单元、3D MOT单元、电离真空计、原子俘获与成像光路、光电倍增管原子荧光探测系统、离子泵和钛升华泵抽气单元等组成。锂源原子炉温度设置为260 ℃，对应的热Li原子饱和蒸气压仅10^–10 Pa。2D MOT的光场由波长为671 nm的4束正交红失谐圆偏振光组成，磁场由一对尺寸为60 mm×20 mm×10 mm的铷铁硼磁铁产生，磁场梯度为0.5 T/m。热原子束被2D MOT预冷却在与两组激光垂直的轴向上形成梭形原子团。3D MOT的光场由波长为671 nm的6束正交红失谐圆偏振光组成，磁场由一对反赫姆霍兹线圈组成，线圈电流15 A时径向磁场梯度为0.1 T/m。冷原子荧光测量选用日本滨松H10721-20型号光电倍增管(PMT)。差分抽气管设计为内径2 mm，长度50 mm，差分效果2个量级。两台离子泵和两台钛升华泵分别连接于2D MOT和3D MOT真空室，离子泵抽速40 L/s。图4为在3D MOT真空腔体中拍摄到的⁷Li冷原子团照片(CF40法兰中心亮点即为冷原子团)，原子个数为10⁸个，其形貌呈现尺寸约0.8 cm×0.8 cm×0.4 cm的椭球状，温度为270 μK。

图4  3D MOT真空腔体中⁷Li冷原子团拍摄照片

3.2 动态流量真空系统

动态流量真空系统用于在3D MOT中产生不同量级大小的动态平衡气体压力，主要由高精度气体微流量计、限流小孔分子流导元件、测量室和抽气室、磁悬浮涡轮分子泵和低温泵抽气机组、热阴极电离真空计和四极质谱计等组成。气体微流量计采用定容-恒压式流量计，其定容室容积为5 L，极限真空度为5×10^–8 Pa，进气小孔采用直径约30 μm的无氧铜激光烧蚀小孔，测量的其对于N₂的分子流导为5.3×10^–9 m³/s。测量室和抽气室为SUS316L加工而成的ϕ250 mm×250 mm的柱形容器，之间由直径33 mm，厚度0.2 mm的抽气小孔连接。抽气室直连抽速2000 L/s的低温泵，测量室通过DN160全金属角阀连接分子泵级联抽气机组，电离真空计和四极质谱计连接在测试室上分别用于真空度测量和气体组分分析及其真空检漏。最终在测量室可获得6×10^–9 Pa极限真空度。

04  实验结果与讨论

利用⁷Li冷原子真空测量装置在磁光阱中和 四极磁阱中开展实验验证研究，主要包括冷原子 损失率实验测量、冷原子损失率系数实验测量、 冷原子真空测量结果与电离真空计的比对验证3部分内容。

4.1 冷原子损失率实验测量

首先，在3D MOT真空室中装载冷原子，冷却光频率锁定在⁷Li原子的∣2²S_1/2,F = 1⟩→∣2²P_3/2,F ′=2⟩跃迁线下红失谐22 MHz，泵浦光频率锁定在⁷Li原子的∣2²S_1/2,F =1⟩→∣2²P_3/2,F ′=2⟩跃迁线下红失谐31 MHz，图5为⁷Li原子 的²S基态和²P激发态能级图。

图5  ⁷Li原子²S基态和²P激发态能级图

冷原子装载和损耗的动力学方程可表示为

式中，N为冷原子个数；R为冷原子装载率(s^–1)；Γ_loss为冷原子-气体分子碰撞损失率(s^–1)；α为冷原子间的非弹性两体碰撞损失率。当3D MOT中冷原子的损耗与装载达到动态平衡时，关闭装载光路，仅观测原子耗散过程，求解(4)式得到：

式中，N₀为稳态t = 0时刻的初始原子数目；C 为考虑到环境杂散光和本底噪声影响而引入的补偿项。

冷原子真空测量过程中需要精确控制各种操作信号的产生、稳定及撤销时间，信号间隔在μs量级，采用Labview编写的时序调控程序对损失率测量中各光路和电路的开关、冷原子团的装载时间、PMT的触发时间等进行精确操控。图6为磁光阱和磁阱中冷原子损失率测量的控制时序，其中红色和黑色线段代表磁光阱中的测量时序，绿色和黑色代表磁阱中的测量时序，V_bg表示PMT测量的背景散射荧光信号，V₀表示动态平衡后的初始装载信号，V_re表示剩余原子再捕获后的信号。

图6  磁光阱和磁阱中冷原子损失率测量控制时序

图7为在2×10^–7 Pa和1×10^–6 Pa真空度条件下，分别采用纯度为99.999%的高纯N₂和Ar在磁光阱和磁阱中测量的⁷Li冷原子的损失率曲线示意图，图中蓝色和红色数据点表示实测的冷原子荧光信号随时间的变化，蓝色和红色曲线是利用(5)式采用Mathematica软件指数拟合的结果。

图7  磁光阱和磁阱中冷原子损失率测量曲线

结果表明，在磁光阱中，2×10^–7 Pa和1×10^–6 Pa下的冷原子损失率Γ_loss分别为4.55×10^–2 s^–1和2.08×10^–1 s^–1，对应的冷原子两体碰撞损失率αα分别为5.43×10^–3 s^–1和5.30×10^–3 s^–1。同样，在磁阱中2×10^–7 Pa和1×10^–6 Pa下的损失率分别为9.39×10^–2 s^–1和3.50×10^–1 s^–1，冷原子间的两体碰撞损失率分别为2.39×10^–3 s^–1·V^–1和1.10 s^–1·V^–1。可见，在该测试状态下，磁光阱和磁阱中，初始阶 段⁷Li冷原子间两体碰撞损失占总体损失的比例αN ²/(Γ_lossN+αN ²)最大分别可达到19%和37%，占比可观，测量模型中考虑冷原子间的碰撞损失影响时才能准确获取损失率Γ_loss。

4.2 碰撞体系损失率系数实验测量

在采用基于长程范德瓦耳斯系数C₆的半经典理论计算损失率系数k_loss时，需要对磁光阱中的原子阱深U进行实验测量，目前磁光阱阱深的测量主要有“光缔合法”^[23]和“重俘获法”^[24]。本实验中，势阱深度的测量采用“光缔合法”，详细测量原理及过程见课题组前期发表文献^[25–27]，本文不再赘述，最后测量得到的阱深值为0.53 K，代入(1)式和(2)式计算得到碰撞体系的损失率系数k_loss。

另外，除了通过实验测量阱深的方法理论计算损失率系数外，在此采用磁阱中基于ab initio量子第一性原理精确计算的损失率系数值^[17]，在特定的真空压力点p下，通过分别测量磁阱和磁光阱中的损失率Γ_MT和Γ_MOT，进而可精确实验测量得到磁光阱中碰撞体系的损失率系数k_loss。测量原理如下：

则磁光阱中的损失率系数k_loss为

式中，Γ_Maj为磁阱中的马约拉纳自旋翻转非理想损失率，其计算公式为Γ_Maj = ℏ/(mc⋅R²)，本实验中冷原子团的直径R为6 mm，从而得到Γ_Maj = 3.6 ×10^–4 s^–1。为消除本底背景真空(6×10^–9 Pa)对测量结果的影响，分别利用气体微流量计向真空测量室中引入1×10^–6 Pa的N₂，Ar，He，H₂四种高纯气体进行碰撞体系损失率系数的实验测量，实验温度由中央空调恒温控制在为(22±1) ℃。

磁阱中测量损失率的实验过程为：冷原子装载完毕后，关闭2D MOT 冷却光和泵浦光、推送光、3D MOT泵浦光，同时仅维持3D MOT 冷却光10 μs，将冷原子全部制备到F = 1 (m_F = –1)能态。然后，关闭所有光场，维持10 G/s梯度磁场 30 μs。最后，再次打开3D MOT冷却光和泵浦光，重新俘获磁阱中剩余冷原子，同时触发PMT观测当前3D MOT 中锂原子荧光信号，并进行数据采集分析。实验中，为避免多次装载原子数目存在的差异性，提高测量精度，在原子装载完成后同步触发PMT，将最后探测的原子荧光信号利用初始信号做归一化处理。表1为最后计算得到的磁光阱 中⁷Li-N₂，⁷Li-Ar，⁷Li-He，⁷Li-H₂的损失率系数，并和采用半经典理论的计算结果进行了对比。

表1 ab initio第一性原理实验测量的损失率系数k_loss和半经典理论计算值比对

对比结果表明，在保持磁光阱阱深不变的前提条件下，对于⁷Li-He和⁷Li-H₂，两种测量方法得到的损失率系数结果较为一致，偏差分别为24%和36%，而对于⁷Li-N₂和⁷Li-Ar，基于阱深测量的半经典理论计算值远远偏离基于ab initio第一性原理的实验测量值，出现了数量级的差异。结合4.3节表2的比对验证结果，表明该半经典理论在计算大质量数分子碰撞体系的损失率系数k_loss时，采用一阶近似计算量子衍射碰撞损失率系数k_g时存在较大的局限性，从而导致最终的碰撞系数k_loss存在近1—2个数量级的偏差，该理论只适用于Li-H₂和Li-He原子结构简单体系k_loss的预估。

表2 冷原子校准的分离规相对N₂的灵敏度因子

4.3 冷原子真空测量比对验证

当测量装置达到6×10^–9 Pa极限本底真空后，利用N₂，Ar，He和H₂四种高纯气体进样开展实验比对验证，通过动态流量真空系统在3D MOT真空室中产生10^–8—10^–5 Pa的动态平衡压力，然后同时采用电离真空计和锂冷原子进行测量，比对电离真空计选用德国Leybold公司的IE514分离规，并于测量前在国防最高真空计量标准装置上对其利用N₂进行了校准。测量过程中，实验温度由中央空调恒温控制在为(22±1) ℃，磁光阱的阱深仍然维持在0.53 K，其损失率系数采用表1中第一性原理的实验测量值，最后，按照3.1节冷原子损失率的测量方法得到不同真空度下的损失率Γ_loss后，利用(1)式计算得到冷原子真空测量结果。图8为利用4种气体测量的不同真空度p下的损失率Γ_loss(其中，内插图为3×10^–8—3×10^–7 Pa真空范围的损失率Γ_loss)，结果表明，各碰撞体系的损失率Γ_loss与真空度之间呈现良好的线性关系，线性拟合优度的皮尔森相关系数和R²值均为1，从而验证了冷原子真空测量的合理性和可靠性。

图8  损失率Γ_loss与冷原子测量真空度p的关系曲线

最后，将冷原子测量的N₂，Ar，He，H₂四种气体的真空度和分离规IE514的测量结果进行比对，得到基于冷原子反演的分离规的灵敏度因子1/CF (1/CF为分离规和冷原子测量结果的比值)，结果如图9所示，其中对于测试气体N₂，电离真空计的测量结果取校准证书对示数修正后的值。结果表明，在3×10^–8—4×10^–5 Pa整个测量范围内，除了在10^–8 Pa和10^–7 Pa两个压力点下，冷原子和分离规测量结果可能受本底极限真空和压力波动影响，灵敏度因子值有所偏离线性外(最大偏差为25%)，其他压力点下的灵敏度因子值具有良好的一致性，相对标准偏差均小于5%。

图9  分离规灵敏度因子与冷原子测量真空度的关系曲线

将冷原子测量的4种气体的灵敏度因子利用N₂的值进行归一化后，和已发表文献^[28–30]中基于传统连续膨胀法校准得到的灵敏度因子进行了比对，如表2所示。结果表明，冷原子校准得到的各测试气体的相对灵敏度因子和传统方法也具有很好的一致性，最大偏差小于8%，从而也间接验证了冷原子真空测量的准确性和可靠性。

4.4 冷原子真空测量不确定度分析评估

在以上实验研究的基础上，依据测量不确定度评定方法^[31]，对冷原子真空测量结果进行不确定分析评估。基于冷原子真空测量基本原理和实验过程，结合(1)式和(7)式，建立冷原子真空测量模型p = f(Γ_loss; Γ_MOT; k_tot; Γ_MT; Γ_Maj; k_B; T; ω)，即：

可见，冷原子真空测量的不确定度主要来源于损失率Γ_loss测量引入的不确定度；损失率系数k_loss测量时磁阱中损失率Γ_MT和Γ_Maj引入的不确定度、第一性原理计算的损失率系数k_tot引入的不确定度、磁光阱中损失率Γ_MOT引入的不确定度；玻尔兹曼常数k_B的不确定度；气体分子温度T测量引入的不确定度；测量室本底真空波动ω引入的不确定度。由于在1×10^–6 Pa真空条件下测量损失率系数k_loss时，磁阱中的Majorana非理想损失率Γ_Maj比原子损失率Γ_MT小3个量级，因此可忽略其影响。则测量模型最终简化为

由于各不确定度分量间不相关，因此相对合成标准不确定度u_c,r(p)为

原子损失率测量不确定度u_r(Γ_loss)由重复性测量、光电倍增管荧光信号测量及曲线拟合不确定度合成得到，分别采用A类和B类评定方法。其中重复性测量不确定度根据6次测量的结果，采用A类评定方法利用贝塞尔公式计算得到，光电倍增管荧光信号测量不确定度由其仪器的非线性决定，采用B类评定方法，根据产品手册和实验观测其非线性小于1%，按照均匀分布计算，则不确定度分量为1%/√3 = 0.6%，原子损失率曲线拟合引入的不确定度小于10^–4，可忽略不计；同理，损失率系数k_loss测量过程中的原子损失率测量不确定度u_r(Γ_MT)，u_r(Γ_MOT)按以上相同方法分析计算得到；第一性原理计算的损失率系数不确定度u_r(k_tot)根据文献[17]得到；玻尔兹曼常数测量不确定度u_r(k_B)可完全忽略不计；真空室内气体分子温度采用铂电阻测量，由于测量装置位于恒温恒湿洁净间内，温度波动最大为1 K，实验测量温度为295 K，则温度变化半宽度为0.7%，按照均匀分布计算，气体分子温度测量不确定度为0.7%/√3= 0.3%。对于本底真空影响引入的测量不确定度，在本底真空6×10^–9 Pa下重复6次测量，通过计算本底真空变化的实验标准偏差来评价其引入的测量不确定度，当测试气体的真空度为10^–5 Pa时，高于本底真空3个量级，不确定度可完全忽略不计。综上，冷原子真空测量不确定度各分量大小及各气体不同真空度所对应的不确定度如下表3汇总。

表 3  冷原子真空测量不确定度汇总表

取包含因子k = 2(置信概率为95%)，则在10^–8—10^–5 Pa范围内，冷原子真空测量不确定度对于N₂为6.6%—4.2%，Ar为6.4%—4.0%，He为7.2%—6.0%，H₂为7.6%—5.4%。

05  结    论

基于超冷原子效应已实现了国际单位制中时间秒的量子化定义，在超高真空测量领域，如何基于冷原子-气体分子微观粒子间的量子化特征实现真空量值的精确反演，是目前国际前瞻性、颠覆性真空测量技术中最具代表性的研究方向，属于多学科交叉量子精密测量的范畴，在量限拓展和测量精度方面具有显著优势。本文利用自主研制的⁷Li冷原子超高真空测量装置，通过操控锂冷原子团的数目、密度、温度、能态和阱深等量子特性开展了真空测量实验和比对验证工作，主要结论如下。

1) 在3D MOT中俘获了2×10⁸ 个⁷Li原子，温度约270 μK，在确定冷原子-气体分子碰撞损失率Γ_loss时，冷原子间的两体碰撞损失率占比可达37%，为实现精确测量，测量理论模型中须考虑该影响量。

2)基于原子阱深测量的半经典理论计算方法和基于ab initio第一性原理的实验测量方法在计算⁷Li-He，⁷Li-H₂这种原子结构简单体系的损失率系数k_loss时，可认为两种方法具有一致性，最大偏差为36%。但半经典理论不适用于⁷Li-N₂，⁷Li-Ar等复杂碰撞体系k_loss的准确计算。

3)在3×10^–8—4×10^–5 Pa真空范围，除10^–8 Pa和10^–7 Pa两个压力点外，冷原子测量的4种气体的灵敏度因子的相对标准偏差均小于5%，基于⁷Li冷原子校准的分离规的相对灵敏度因子和传统方法具有良好的一致性，最大偏差小于8%，充分验证了冷原子量子真空测量的准确性和可靠性。

4)在3×10^–8—4×10^–5 Pa真空范围，对应于N₂，Ar，He，H₂四种气体，⁷Li冷原子真空测量的不确定度分别为6.6%—4.2%，6.4%—4.0%，7.2%—6.0%和7.6%—5.4%，包含因子均取k = 2。

本项研究工作为后续我国在超高极高真空测量领域建立自主知识产权的量子真空计量标准，研制小型高性能冷原子真空测量传感器具有较大促进作用，研究成果在空间科学探测、高能粒子加速器、引力波测量等军民高新技术领域具有广泛应用前景。

感谢华东师范大学武海斌教授、武跃龙研究员在锂冷原子冷却囚禁系统搭建中的支持；感谢吉林大学王德权教授团队在锂原子-气体分子碰撞体系损失率系数理论计算研究中的指导。

参考文献

[1] Li D T, Cheng Y J, Feng Y, Lu Y W, Zhao L, Xu J 2009 Chin. J. Vac. Sci. Technol. 29 522 (in Chinese) [李得天, 成永军, 冯焱, 卢耀文, 赵澜, 徐婕 2009 真空科学与技术学报 29 522]

[2] Li D T, Wang Y J, Zhang H Z, Xi Z H, Li G 2021 Space Sci. Tech. 23 7592858

[3]Liu J, Wang G, Hu Y M, Zhang T, Luo Z R, Wang Q L, Shao L J 2016 Chin. Sci. Bull 61 1502  (in Chinese) [刘见, 王刚, 胡一鸣, 张腾, 罗子人, 王晴岚, 邵立晶 2016 科学通报 61 1502]

[4] Wang Y Z, Chen X 2007 Vacuum Technology (Beijing: Beihang University Press) p6  (in Chinese) [王欲知, 陈旭 2007 真空技术 (北京: 北京航空航天大学出版社) 第6页]

[5] Jenninger B, Anderson J, Bernien M, Bundaleski M, Dimitrova H, Granovskij M, Illgen C, Setina J, Jousten K, Kucharski P, Reinhardt C, Scuderi F, Silva R A S, Stöltzel A, Teodoro O M N D, Trzpil-Jurgielewicz B, Wüest M 2021 Vacuum 183 109884

[6] Fedchak J A, Abbott P J, Hendricks J H, Arnold P C, Peacock N T 2018 J. Vac. Sci. Technol. A 36 030802

[7] Prentiss M, Cable A, Bjorkholm J E, Chu S, Raab E L, Pritchard D E 1988 Opt. Lett. 13 452

[8] Scherschligt J, Fedchak J A, Barker D S, Eckel S, Klimov N, Makrides C, Tiesinga E 2017 Metrologia 54 125

[9] Eckel S, Barker D S, Fedchak J A, Klimov N N, Norrgard E, Scherschligt J, Makrides C, Tiesinga E 2018 Metrologia 55 182

[10] Makrides C, Barker D S, Fedchak J A, Eckel S, Tiesinga E 2019 Phys. Rev. A 99 042704

[11] Makrides C, Barker D S, Fedchak J A, Scherschligt J, Eckel S, Tiesinga E 2020 Phys. Rev. A 101 012702

[12] Makrides C, Barker D S, Fedchak J A, Scherschligt J, Eckel S, Tiesinga E 2022 Phys. Rev. A 105 029902

[13] Makrides C, Barker D S, Fedchak J A, Scherschligt J, Eckel S, Tiesinga E 2022 Phys. Rev. A 105 039903

[14] Barker D S, Klimov N N, Tiesinga E, Fedchak J A, Scherschligt J, Eckel S 2021 Measurement: Sensors 18 100229

[15] Barker D S, Acharya B P, Fedchak J A, Klimov N N, Norrgard E B, Scherschligt J, Tiesinga E, Eckel S 2022 Rev. Sci. Instrum. 93 121101

[16] Barker D S, Fedchak J A, Kłos J, Scherschligt J, Sheikh A A, Tiesinga E, Eckel S 2023 AVS Quantum Sci. 5 035001

[17] Kłos J, Tiesinga E 2023 J. Chem. Phys. 158 014308

[18] Booth J L, Shen P, Krems R V, Madison K W 2019 New J. Phys. 21 102001

[19] Shen P, Madison K W, Booth J L 2020 Metrologia 57 025015

[20] Shen P, Madison K W, Booth J L 2021 Metrologia 58 022101

[21] Makhalov V B, Martiyanov K A, Turlapov A V 2016 Metrologia 53 1287

[22] Jousten K 2016 Handbook of Vacuum Technology (Weinheim: Wiley-VCH) p634

[23] Dongen J V, Zhu C, Clement D, Dufour G, Booth J, Madison K 2011 Phys. Rev. A 84 022708

[24] Kawanaka J, Shimizu K, Takuma H 1993 Appl. Phys. B 57 113

[25] Zhang S Z, Sun W J, Dong M, Wu H B, Li R, Zhang X J, Zhang J Y, Cheng Y J 2022 Acta Phys. Sin. 71 094204  (in Chinese) [张苏钊, 孙雯君, 董猛, 武海斌, 李睿, 张雪姣, 张静怡, 成永军 2022 物理学报 71 094204]

[26] Wu X M, Cheng Y J, Dong M, Sun W J, Zhang S Z, Ma Z Y, Li Y P, Jia W J, Feng T Y, Wu C Y 2023 Vacuum 207 111561

[27] Sun W J, Wu X M, Cheng Y J, Ma Z Y, Jia W J, Zhang Y F, Zhang R F, Wu C Y, Feng C, Luo H G 2024 Vacuum 222 113079

[28] Hajime Yoshida, Kenta Arai 2018 J. Vac. Sci. Technol. A 36 031604

[29] Alper Elkatmis, Rifat Kangi 2019 Measurement 131 269

[30] Yoshida H, Arai K, Kobata T 2014 Vacuum 101 433

[31] Bich W 2014 Metrologia 51 S155

(参考文献可上下滑动查看)