范栋 习振华 贾文杰 成永军 李得天

(兰州空间技术物理研究所真空技术与物理重点实验室, 兰州 730000)

为进一步提高真空量值的复现性和准确性, 最新研究采用量子技术实现对真空量值的测量与表征.该方法利用Fabry-Perot 谐振腔实现腔内气体折射率的精密测量, 并反演出气体密度, 进而获得对应的真空量值,其中气体折射率的测量是影响真空量值准确性的关键.本文基于第一性原理, 利用从头计算理论计算了在已知压力和温度条件下的氦气折射率, 给出腔内气体压力与折射率关系的表达式, 并利用基于Fabry-Perot 激光谐振腔的真空测量装置, 通过双腔谐振激光拍频精确测量了充气前后谐振激光频率的变化, 测出了氦气折射率, 并分析了测量不确定度.将理论计算值与实验测量值进行了对比分析, 得出了制约准确度提高的主要因素, 并提出了修正方法.

1 引 言

2018 年, 国际单位制(International System of Units, SI)中7 个基本量全部由物理常数定义[1].面对国际计量体系的历史性变革, 诸多计量参数对应的基于经典物理学的实物标准向“自然标准”即量子标准发展, 利用量子跃迁等现象复现计量单位, 用相同装置实现任何时间、任何地点量值复现,消除宏观参数不稳定产生的影响, 将准确度提高数个量级.目前, 时间、长度、电学等方面的量子基准已逐步建立.与此同时, 采用光学干涉、冷原子碰撞损失、吸收光谱等基于物质本征物理性质复现真空量值的新方法、新概念具备突破传统技术性能极限的应用潜力, 可在多个地点及不同时间复现真空量值, 已成为国际上新的研究热点[2].其中, 基于Fabry-Perot(F-P)激光谐振腔测量非极性稀薄气体折射率, 进而反演真空度的测量方法为替代水银压力计—真空量值基准提供了可能.该方法在减小测量不确定度的同时, 可实现真空量值的扁平化传递, 为真空量值提供了新的溯源途径[3−5].近年来, 美国国家标准与技术研究院(NIST)、瑞典于默奥大学(Umeå universitet)和日本国家先进工业与技术研究所(AIST)等对基于F-P 激光谐振腔测量折射率反演真空度的真空测量技术进行了研究[6−9], 该技术主要面临的挑战是由于温度和充气导致谐振腔腔长变化, 从而影响折射率的测量准确度[10].国内此类研究尚属起步阶段[11,12].

量子动力学将非极性气体宏观介电常数与原子微观极化参数相联系.近年来, 量子力学关于单原子分子(氦气、氩气等)电磁学性质和热力学性质从头计算(ab initio)理论及方法的迅速发展, 相关参数的计算精度不断提高(例如, He 的摩尔极化率理论计算相对不确定度达到 2 × 10–7), 可进一步减小上述方法的测量不确定度[13−20].

本文采用基于F-P 激光谐振腔的真空测量装置开展非极性稀薄氦气折射率反演真空量值研究,对比理论计算折射率与测量折射率, 并评定各不确定度分量及权重系数, 分析制约准确度提高的主要因素并提出修正方法.

2 折射率理论计算

氦气是除氢气以外最简单的气体原子, 因此其微观物理参量可以达到最高精度的计算.在折射率理论值的计算过程中, 首先对氦原子的摩尔极化率、摩尔磁化率、第二维里系数和第二介电维里系数进行计算.原子极化率可以写为波长的展开式:

式中, λ 为激光波长(单位为Å), α 为原子极化率(原子单位制), A0为静极化率, Ak( k =2,4,6 )为极化率色散系数, 其值如表1 所示.

表1 He 极化率的展开系数(原子单位制)[17]Table 1.Cofficients in the expansion of the polarizability of Helium.

在激光波长λ = 780 nm 时, 计算得出He 极化 率α(780 nm) = 2.0583140283 × 10–25cm3.由于原子磁化率χ 比原子极化率α 小5 个数量级,因此可忽略χ 的色散修正, 只需考虑其静态成分,静态磁化率χ = –3.14062677086 × 10–30cm3[21].为使气体密度以mol/cm3为单位, 需将α 和χ 转化为摩尔极化率Aϵ= 4πNAα/3 和摩尔磁化率Aµ=4πNAχ/3, 其中, NA为阿伏伽德罗常数.经计算,Aϵ= 0.5192196805 cm3/mol, Aµ= –0.0000079224 cm3/mol.

根据Bich[19]和Rizzo[20]公布的各温度点下第二维里系数和第二介电维里系数的值, 对其进行拟合, 在250—350 K 温度范围内, 得出第二维里系数B(T)和第二介电维里系数Bϵ(T)关于温度的线性表达式:

式中T 为温度.在温度T = 299.1485 K 时, 可得B(T)=11.818 cm3/mol,Bϵ(T)=–0.053 cm6/mol2.

气体压力(单位为Pa)和密度(单位为mol/m3)的关系可通过维里状态方程来表示:

式中, p 为气体压力, NA为阿伏伽德罗常数, kB为玻尔兹曼常数, T 为温度, ρ 为气体密度, B(T)为第二维里系数.

可通过Lorentz-Lorenz 方程建立气体密度与折射率的联系, 用摩尔极化率、摩尔磁化率(以cm3/mol 为单位)及第二介电维里系数表示气体密度(以mol/cm3为单位), 有:

式中, ρ 为气体密度, n 为气体折射率, Aϵ为摩尔极化率, Aµ为摩尔磁化率, Bϵ(T)为第二介电维里系数.

结合(4)式和(5)式, 可得到压力与折射率的关系:

(6)式中, 各项系数表达式分别为:

式中, NA= 6.02214076 × 1023mol–1, kB= 1.380649 ×10–23J/K, T = 299.1485 K.计 算 得 出, c1=3.19363 × 109; c2= 4.83514 × 1010; c3= –3.32951 ×109; c4= 5.71882 × 107.在p = 100166 Pa 时, 第二项分别比第三项和第四项高5 和7 个数量级, 因此在计算折射率时, 可忽略第三项和第四项.则反写(6)式可得出折射率关于压力的简化表达式:

根据(11)式, 即可通过Aϵ、Aµ、Bϵ(T)、B(T)以及压力和温度计算出氦气折射率, 其结果如表2所示.

根据表2 所示的结果, 对折射率计算不确定度进行分析.(6)式中, 由于摩尔极化率比摩尔磁化率高5 个数量级, 且在p = 100166 Pa 时, 第一项比第二项高4 个数量级, 因此可忽略折射率非线性高阶项和摩尔磁化率对折射率计算不确定度的影响.则(6)式可简化为折射率n 的计算表达式:

式中, p 为已知气压值.

(12)式中, Aϵ的相对不确定度ur, Aϵ= 2 × 10–7,修正后的SI 中, kB和NA没有不确定度[22], 折射率计算的相对不确定度为 ur,n1.经分析可知, 当压力点为105Pa 时, 其折射率计算不确定度最大, 随着气压降低, 不确定度逐渐减小, 因此选取不确定度最大的点作为不确定度评定点.真空度p =100166 Pa 时, 利用(13)式进行折射率计算不确定度分析:

式中, uc,n1为折射率理论值n 的标准不确定度,uc,Aϵ为摩尔极化率的标准不确定度.

分析得出, 折射率理论值n 的相对不确定度ur,n1=6.27×10−12.

3 实验原理与方法

3.1 实验装置与原理

本文所用实验装置结构如图1 所示, 探测激光器为两个中心波长为780 nm 的外腔室式半导体激光器(ECDL, TOPTICA, DL pro).腔体材料为超低热膨胀玻璃(ULE), 采用双腔结构.通过调节透镜组进行模式匹配, 之后使用信号发生器(RIGOL,DG4102)发出的正弦同步信号对电光调制器(EOM,THORLIBS, CRDS)进行相位调制, 使透射信号出现边带.然后利用信号发生器产生的同频信号和腔内反射信号混频产生误差信号, 通过误差信号利用PDH(Pound Drever Hall)锁频装置将激光锁定在腔纵模上, 并利用频率计(KEYSIGHT, 53220A)测量两组激光的拍频.真空泵组包括机械泵(Leybold, SCROLLVAC)、分子泵(Leybold, TURBOVAC)和离子泵(Gamma Vacuum, TiTan lon Pump), 先用机械泵抽气实现粗真空, 再同时打开分子泵进一步抽气提高真空度, 当抽气至所需真空度时, 打开离子泵进行真空度维持.气瓶后连接一个液氮冷阱, 可去除杂气, 气体通过微漏阀控制流速进入腔体, 腔体外部连接电容薄膜真空计(INFICON Cube CDGSci)以检测充气压力.

表2 折射率计算结果Table 2.Calculation results of refractive index.

图1 基于F-P 激光谐振腔的真空测量装置结构图 (a)部分为光路图;(b)部分为气路图 (1.激光器; 2.分束光纤; 3.激光准直器; 4.1/2 波片; 5.格兰棱镜; 6.透镜组; 7.电光调制器; 8.光隔离器; 9.偏振分光棱镜; 10.1/4 波片; 11.高反镜; 12.透镜; 13.光电放大探测器; 14.PDH 锁频装置; 15.光电探测器; 16.频率计; 17.检测腔; 18.参考腔; 19.电容薄膜真空计; 20.气瓶; 21.冷阱; 22.离子泵; 23.电离规; 24.分子泵; 25.机械泵.红色实线为光路; 黑色实线为光纤; 黑色虚线表示反馈作用; 蓝色实线为气路.)Fig.1.Structure diagram of vacuum measuring device based on F-P cavity.

激光锁定在F-P 腔纵模上时, 腔内激光谐振频率可表示为

式中, q 为整数, c 为真空中的光速, n 为气体折射率, L0为谐振腔长度.

将充气前后的谐振腔内频率值代入(14)式求差, 即可推导出绝对气体折射率n 的表达式:

式中, ∆νcm为腔内谐振激光频率的变化, νf为检测腔内谐振激光末态频率.

考虑到充气过程中激光纵模模数变化和气体压力对腔体长度的影响, (15)式可写为

式中, νi为检测腔谐振激光初态频率, νf为检测腔谐振激光末态频率, Δm 为检测腔内谐振激光纵模模数变化量, νFSR是自由光谱范围, Lf(p)为充气后腔体实际物理长度, Li为充气前腔体物理长度.

(16)式中第二项长度变化因子Lf(p)和Li主要是由于充气后气体压力导致的腔长变化:

式中, ∆p 为腔内气压变化量, K 为ULE 材料的体积模量.

在腔体高真空(n = 1)状态下测量谐振激光初态频率νi, 腔体中充入一定压力的氦气时测量谐振激光末态频率νf, 即可得到充气压力点下的频率差∆νcm, 最后代入(16)式计算出每个频率差对应的气体压力.

3.2 折射率测量方法

为准确测量谐振腔内气体折射率, 利用PDH锁频模块将激光频率锁定在谐振腔内, 并通过SRS温度监测模块及Pt100 铂电阻将谐振腔内温度控制在ULE 材料的零膨胀温度点299.1485 K.首先将检测腔和参考腔分别抽气至10–5Pa, 并测量初态拍频和检测腔内谐振激光初态频率; 其后将检测腔充气至各压力点, 并测量末态拍频和检测腔内谐振激光末态频率.则结合(16)式和(17)式可得出折射率计算表达式:

式中, νr为参考激光频率, νbeati为充气前两束激光拍频, νbeatf为充气后两束激光拍频.

利用(18)式得到折射率测量结果如表3 所示.

根据表3 所示的结果, 对折射率测量不确定度进行评定分析.由于(18)式中, 分量的不确定度量级为10–15, 故可以忽略不计.经分析后发现当测量压力点为101.5 Pa 时, 气体折射率n – 1 =3.93406 × 10–8, 其测量不确定度最大, 随着腔内气体压力增大, 其测量不确定度逐渐降低, 因此选取该点作为气体折射率测量不确定度评定点.在折射率测量不确定度评定过程中, (18)式可写为

表3 折射率测量结果Table 3.Text results of refractive index.

(19)式中, 分子部分为充气前后检测腔谐振激光频率变化量, 其相对不确定度为ur,1, 分母部分为检测腔末态谐振激光频率, 其相对不确定度为ur,2, 折射率测量不确定度为 ur,n2.实验中, 参考激光频率νr和自由光谱范围νFSR用波长计测量, 其装置测量不确定度uc,r和uc,FSR均为10–6; 初态拍频νbeati和末态拍频νbeatf用频率计测量, 其装置测量不确定度uc,beati和uc,beatf均为10–6; 不确定度评定方法如下:

(19) 式中各项参数的数值分别为: νFSR=1497.9916 MHz; ∆m = 0; νr= 384366590 MHz;νbeati= 1420 MHz; νbeatf= 1434, 74 MHz.经 评定, 分子部分的测量不确定度ur, 1= 9.59 × 10–8,分母部分的测量不确定度ur, 2= 3.68 × 10–15, 各分量相互独立且不相关, 由(22)式:

可计算出当气体压强为101.5 Pa 时, 其折射率测量不确定度 ur,n2=9.59×10−8.

4 分析与讨论

结合第2 节和第3 节得到的折射率理论和测量结果, 对折射率与检测腔内真空度之间关系进行对比分析, 其结果如图2 所示.

图2(a)显示腔内氦气压力在102—103Pa 范围内, 氦气的折射率测量值大于理论值, 其主要原因是腔体材料中有余氦, 在腔内气压较低时, 腔体对氦的吸气速率小于放气速率, 因此由氦气脱附引起腔内氦气量增多, 导致了折射率测量值偏大, 这验证了ULE 材料存在渗氦现象[23].图2(b)显示腔内氦气压力在103Pa 附近时, 氦气折射率测量值开始小于理论值, 并且随着腔内气压增大, 氦气折射率测量值与理论值偏差也逐渐增大(偏差值从4.25267 × 10–9到2.38851 × 10–6), 其主要原因是随着氦气压力增大, 腔体渗氦量增加, 并且渗氦和返氦速率均增大, 但腔体对氦的吸气速率大于放气速率, 导致折射率偏差逐渐增大.图2(c)显示腔内压力在102—105Pa 范围内, 氦气折射率理论值与测量值偏差随腔内气压的总体变化情况.

图2 (a)102−103 Pa 范围内折射率理论值与测量值对比图; (b) 103−105 Pa 范围内折射率理论值与测量值对比图;(c) 102−105 Pa 范围内折射率理论值与测量值总对比图Fig.2.(a) Comparison between theoretical and measured values of refractive index in the range of 102−103 Pa; (b) comparison between theoretical and measured values of refractive index in the range of 103−104 Pa; (c) total comparison between theoretical and measured values of refractive index in the range of 102−105 Pa.

分析表明, 在利用基于F-P 激光谐振腔的真空测量装置测量氦气折射率实现新的真空标准中,由于ULE 材料渗氦导致腔体长度变化对折射率测量不确定度的影响较大.为降低此影响, 需对测量结果引入必要的腔长修正.由于折射率测量不确定度在10–8量级, 此量级下无法根据已知的ULE 材料特性计算出充气前后腔长的变化, 所以可利用氦气折射率的理论计算值来量化由渗氦所致腔体形变导致的测量折射率偏差.充气后, 谐振腔实际物理长度为

式中, Lf(p)为充气后腔体实际物理长度, Li为充气前腔体物理长度, ∆LG(p)为气体压力导致的腔长变化量, ∆LHe(p)为腔体渗氦导致的腔长变化量.Lf(p), ∆LG(p), ∆LHe(p)是压力的函数.

将(23)式代入(14)式和(16)式得到加腔长修正后的折射率表达式:

式中, ∆νcm为激光谐振频率变化量, φ(p)为与腔体渗氦相关的折射率修正系数:

根据本文得到的折射率理论值和计算值, 即可计算出103—105Pa 真空度范围内本实验所用实验装置的折射率修正系数, 如表4 所示.

对表4 中数据进行拟合, 得到: 在1075—7031 Pa压力范围内, 折射率修正系数关于腔内气压的表达式为

在7031—100166 Pa 压力范围内, 折射率修正系数关于腔内气压的表达式为

为了进一步量化由气体压力导致的谐振腔腔长变化, 可以推导出气体压力所致腔体相对变化量关于折射率的表达式.气体压力所致腔体形变量关于气体压力的表达式为

式中, Li为充气前腔体的物理长度, p 为气体压力,K 为ULE 材料的体积模量.则气体压力所致腔体相对变化量为

表4 折射率修正系数Table 4.Refractive index correction coefficient.

将(30)式代入(4)式和(5)式中的压力和密度的一阶线性项, 并忽略密度一阶表达式中摩尔磁化率的影响, 即可得到气体压力所致腔体相对变化量关于折射率的表达式:

令:

式中, K 为ULE 材料的体积模量.

则可得出, 充气后气体压力所致腔体相对变化量关于折射率的表达式为

式中, ε 是无量纲腔长修正参数, 在本实验条件下,其值为3.12 × 10–2.

5 结 论

本文首先在压力和温度已知的情况下, 基于第一性原理利用摩尔极化率、摩尔极化率和维里系数等微观物理参量对氦气折射率理论值完成精确计算, 得出腔内气体压力与折射率关系的表达式并分析得出折射率理论值的相对不确定度为6.27 ×10–12.其后利用基于F-P 激光谐振腔的真空测量装置, 对102—105Pa 范围内的氦气折射率进行测量, 评定其测量不确定度为9.59 × 10–8.最后对氦气折射率理论值和测量值对比分析, 得出氦气折射率测量不确定度的主要原因是由渗氦引起的谐振腔腔体形变, 因此解决渗氦问题是建立新的真空标准的关键.文中对氦气折射率测量中由腔体渗氦导致的腔长变化进行了修正, 得到103—105Pa 真空度范围内各压力点下的折射率修正系数及其关于腔内气压的表达式, 并进一步量化气体压力导致的腔长变化量, 得到气体压力所致腔体相对变化量与折射率的关系, 计算出腔长修正参数ε = 3.12 ×10–2.在以后利用F-P 谐振腔的氦气折射率测量实验中, 可利用本文所给的各压力点下的折射率修正系数对折射率测量结果进行修正, 消除腔体渗氦对折射率测量的影响, 从而得到高准确度的气体压力.随着量子力学关于单原子分子电磁学性质和热力学性质从头计算理论及方法的迅速发展, 相关参数计算精度提高, 将来可利用其他非极性气体(如氖)实现气体折射率至气体压力的精确反演, 降低由于气体扩散导致的折射率测量不确定度, 实现真空计量新标准.