-
充有碱金属原子的气室结构是原子磁力仪的核心传感元件。一束x轴向的圆偏振光对气室内的碱金属原子进行极化,极化后的原子在y轴方向外加磁场B的作用下能级发生塞曼分裂,此时原子受磁场作用和磁矩作用,绕着y轴以一个微小的角度θ做拉莫尔进动。此时,如果用一束线偏振的探测光沿着z轴方向穿过气室,则该探测光与极化原子相互作用后的偏振面将发生一个θ角度的偏转,如图1(a)所示。通过测量θ的大小即可探测磁场的信息。基于该测量原理,设计的测量方案如图1(b)所示。该方案中采用垂直的泵浦-探测方式,圆偏振泵浦光通过线偏振激光和起偏器和四分之一波片实现。初始强度为I0的线偏振探测光从气室内出射后,经过偏振分束器进行分束,两束光以45°角分离,最终经过PC对采集到的信号进行解调得到磁场信息。由光电转换原理,测量到的两路电信号强度I1,I2=I0−I1可表示为:
图 1 (a) SERF原子磁力仪原理;(b) SERF原子磁力仪测量方案
Figure 1. (a) Mechanism of SERF atomic magnetometer; (b) Measurement of setup of SERF atomic magnetometer
$${I_1} = {I_0}{\sin ^2}\left( {\theta - \frac{\pi }{4}} \right)$$ (1) $${I_2} = {I_0}{\cos ^2}\left( {\theta - \frac{\pi }{4}} \right)$$ (2) 由于光旋角θ<<1,因此可近似得到:
$$\theta = {{{I_1} - {I_2}} / {2\left( {{I_1} + {I_2}} \right)}}$$ (3) 从光旋角的测量原理可知,气室内探测光的偏振旋转角越大,测量达到的磁场强度越精确,原子磁力仪的灵敏度也越高。因此,增加光学偏振旋转角是提高SERF原子磁力仪探测灵敏度的最直接方法。由于测量过程中需要保证钾原子的气体状态,因此需要对气室进行加热,高温增加了原子的碰撞概率。为了降低气室内碱金属气体原子的碰撞和原子的高浓度引起的淬灭现象,在气室内充入一定量的缓冲气体和惰性气体N2、He等。
根据参考文献 [12],光旋角可近似表示为:
$$\theta {\rm{ = }}\frac{{\rm{1}}}{{\rm{4}}}{\textit{z}}{r_e}cfnD\left( v \right)$$ (4) 式中:z、re、c、f、n、D(v)分别为激光有效传输距离、经典电子半径(2.8×10−13 cm)、真空光速(3×108 m/s)、激光共振强度(一般为1/3)、碱金属原子数密度和原子共振线宽。公式(4)表明,碱金属原子数密度n和有效传输距离z与光旋角成正比关系。增加碱金属原子数密度和有效传输距离可扩大光旋角的测量值。如图2所示,随着温度的上升,碱金属原子数密度增加[4, 13],自旋交换碰撞导致的共振线宽变窄(黑色线),有利于提高探测灵敏度。但是,过大的碱金属原子数密度会增加原子自旋碰撞几率(红色线),使得原本活跃的原子发生退极化,减小出射探测激光的光旋角,原子磁力仪的灵敏度下降。因此,通过增加碱金属原子浓度的途径提高探测灵敏度的效果有限,增加激光在气室内的有效传输距离有望得到较大的光旋角。
-
FP腔是许多光学设备的基本组成部分,通常由一对平行的半透半反镜组成,利用激光在腔内往返实现光束谐振 [14-16],激光在腔内传输时光路稳定,不易畸变。图3是文中理论分析和仿真计算采用的FP腔型气室结构。气室长度1 cm,气室内的碱金属选用钾原子,因为钾原子的自旋碰撞横截面最小,引起的弛豫时间最小,拥有最高的灵敏度[17]。这里,笔者考虑到了气室内充入的N2和He缓冲气体的影响作用,钾原子的线宽为30.4 GHz。由于钾原子达到SERF态需要180 ℃的高温,气室温度设定为180 ℃,根据图2曲线可得到原子数密度约为5×1013 cm−3。当激光在气室内传输一次时,根据公式(4)可达到此时的光旋角θ为0.0435 rad。在图3所示的结构中,探测激光在FP腔内由于谐振作用多次往返,有效传输距离z变长,可实现出射光旋角的倍增效果。
下面笔者通过激光谐振和传输原理[18]从理论推导方面对FP腔的弱磁倍增作用进行详细的分析。
当探测激光束进入FP腔后,在腔内激发大量激光模式,第N阶模式的激光从FP腔内出射前共反射了2N−1次。由于探测激光波长位于钾原子吸收线附近,所以激光在腔内传输时,碱金属原子的自旋碰撞会导致激光被原子吸收,同时FP腔制作工艺的精细程度等原因也会使激光腔长发生失谐。考虑到上述两个因素的影响,该出射激光模式的相位可表示为
$\left( {2N - 1} \right)\left( {\Delta \varphi - \alpha l} \right)$ ,其中,Δφ、α、l分别是FP腔失谐参数、碰撞吸收系数和气室长度。则该激光模式的输出振幅为:$${A^N} = {A_0}T{R^{\left( {N - 1} \right)}}\exp [i\left( {2N - 1} \right)\left( {\Delta \varphi - \alpha l} \right)]$$ (5) 式中:A0为入射探测激光的振幅;T、R=1-T分别为两个反射镜的功率反射率和功率透射率。这里,忽略反射镜的吸收等损耗。根据法拉第旋转效应,出射激光的偏转角为
$\left( {2N - 1} \right)\theta $ ,则第N阶模式激光x和y分量的振幅为:$$\left( {\begin{array}{*{20}{c}} {{{A}}_{{x}}^N} \\ {{{A}}_{{y}}^N} \end{array}} \right){{ = }}\left( {\begin{array}{*{20}{c}} {{{{A}}^N}{{{\rm cos}(2}}N{{ - 1)\theta }}} \\ {{{{A}}^N}{{{\rm sin}(2}}N{{ - 1)\theta }}} \end{array}} \right)$$ (6) 出射探测光两个偏振分量的琼斯矢量矩阵形式为:
$${{J}}\left( {{{\psi ,\varphi }}} \right){{ = }}\left( {\begin{array}{*{20}{c}} {{{{\rm cos}\psi }}} \\ {{{{\rm exp}( - i\varphi ){\rm sin}\psi }}} \end{array}} \right)$$ (7) 式中:ψ,φ分别为输出激光偏振旋转角和振动方向相对于传播方向的方位角。初始线偏探测激光的琼斯矩阵可表示为
$\left( {\begin{array}{*{20}{c}} {{1}} \\ {{0}} \end{array}} \right)$ 。从FP腔出射后,激光的x,y分量分别为:$${A_x} = \sum\limits_{N = 0}^\infty {A_x^N} = |{A_x}| \times \exp ( - i{\varphi _x})$$ (8) $${A_y} = \sum\limits_{N = 0}^\infty {A_y^N} = |{A_y}| \times \exp ( - i{\varphi _y})$$ (9) 为表示x,y分量与出射光旋角和相位角的关系,对公式(8)和(9)进行代数计算,发现:
$$\frac{{{A^ * }_x{A_y} + {A_x}{A^ * }_y}}{{{{\left| {{A_x}} \right|}^2} + {{\left| {{A_y}} \right|}^2}}} = \sin 2\psi \times \cos \varphi $$ (10) 根据公式(10)的代数展开项,出射探测激光旋转角和相位角分别为:
$$\psi = \frac{1}{2}\arcsin \left( {\frac{{\sin 2\theta \sqrt {{{\left( {1 - {R^2}} \right)}^2} + 4{R^2}{{\sin }^2}2\left( {\Delta \varphi - \alpha l} \right)} }}{{{R^2} + 1 - 2R\cos 2\left( {\Delta \varphi - \alpha l} \right)\cos 2\theta }}} \right)$$ (11) $$\varphi = \arccos \left( {\frac{{1 - {R^2}}}{{\sqrt {{{\left( {1 - {R^2}} \right)}^2} + 4{R^2}{{\sin }^2}2\left( {\Delta \varphi - \alpha l} \right)} }}} \right)$$ (12) 若忽略吸收系数α和失谐参数Δφ,公式(11)可简化为:
$$\psi = \arctan \left( {\frac{{1 + R}}{{1 - R}}\tan \theta } \right)$$ (13) 从公式(13)可以看出,理想状态下,出射光旋角的积累量大小主要由FP腔的参数R决定。文中,输出光旋角ψ和激光在腔内单次传输的光旋角θ的比值定义为FP腔的倍增因子。从公式(11)可以看出,影响原子磁力仪灵敏度的增益因子受FP腔的吸收和失谐影响。图4是根据理论分析结果公式(11)和(12)得到的倍增因子在吸收和失谐作用下的变化趋势。
图 4 (a)吸收和(b)失谐参数对出射光旋角的影响
Figure 4. Influence of (a) absorption and (b) off-resonance parameter on output light rotation angle
从图4可以看出,随着吸收系数α和失谐参量Δφ的增大,出射激光的光旋角ψ非线性地减小。如图4(a)所示,仅考虑吸收时,α从0增长到0.5 m−1,光旋角从0.69 rad减少到0.28 rad,倍增因子从16减少到6.5;如图4(b)所示,当只考虑FP腔失谐时,Δφ从0增长至π/16,光旋角从0.69 rad减少至0.21 rad,倍增因子从16减少到4.8。上述结果表明,PF腔的理想倍增因子为16。通过理论分析发现,碰撞吸收和腔失谐都会减弱倍增因子。在碱金属气室中,吸收主要来自碱金属原子的碰撞,包括原子间的自旋碰撞和原子与气室壁的碰撞,碰撞弛豫主要和碱金属内的原子密度和压强相关。因此需控制碱金属气室内的压强和原子数密度。除此之外,为防止失谐参数过大而影响倍增因子,FP腔的制作工艺应保持在较高精度。
-
谐振腔原理和激光传输理论表明吸收和失谐对FP腔型气室的倍增因子有削弱作用。为深入研究碰撞吸收和腔壁失谐对原子磁力仪性能的影响,这里从最基本的光场形式出发,用激光传递矩阵法对腔内光场进行数值计算。仿真以公式(8)、(9)为模型,使用的主要变量值如表1所示。数值计算中对偏振光的光场进行循环计算,假定每一次的腔长距离传输结束后激光受到气室内原子的吸收和腔失谐的作用而增加一定量的相移。经过多次迭代,数值计算结果如表1所示。
表 1 光场数值计算参数
Table 1. Numerical calculation parameters of optical field
Parameters Value Power reflectivity(R) 90% Length of cell(l)/cm 1 Alkali number density(n)/cm−3 5×1013 Radius of electron(re) 2.8×10−13 Speed of light in vacuum(c)/m·s−1 3×108 Resonance strength(f) 1/3 FWHM of resonance(D(v))/GHz 30.4 Power of laser (P)/mW 10 当探测激光在FP腔内往返传输时,激光与原子的相互作用距离会逐步增加,光旋角初始阶段近乎线性然后趋于一个最大值,如图5所示。图中3条曲线分别对应3种情况:黑色线为吸收和失谐均为0时的出射光旋角变化曲线,光旋角从ψ=0.0435 rad增至ψ=0.692 rad,最大倍增因子为15.9;当考虑碰撞吸收和腔失谐时倍增因子下降。紫红色线是不考虑失谐,吸收系数α=0.04 m−1时的光旋角变化曲线,出射激光光旋角为0.458 rad,倍增因子10.5,比最大倍增因子低了近1/3,出射激光的偏振态的琼斯矩阵形式为
$\left( {\begin{array}{*{20}{c}} {0.{\rm{896\;9}}} \\ {0.{\rm{442\;1}} + 0.00{\rm{3\;4}}i} \end{array}} \right)$ ;棕色曲线是吸收系数为0,失谐Δφ=π/32时的变化曲线,出射激光光旋角是0.373 rad,倍增因子8.6,与理想倍增因子相比,减少了近50%,此时出射激光偏振态为$\left( {\begin{array}{*{20}{c}} {0.{\rm{931\;2}}} \\ {0.{\rm{096\;9}} - 0.{\rm{35}}1\;3} \end{array}} \right)$ 。激光在腔内传输时碰撞吸收和FP腔失谐都会抑制倍增因子,所以考虑两个影响因子都不为0的情况并用控制变量的方法进行了数值仿真。如图6所示,列出了吸收系数分别为0.01 m−1,0.1 m−1和0.3 m−1以及失谐量为π/64和π/32共计6种情况的计算结果。从这些曲线变化趋势可以看出吸收对倍增因子的削弱作用明显:(1)首先,频率失谐量为π/64时,随着吸收系数从0.01 m−1增长至0.3 m−1,倍增因子从11.2减少到4.5;当频率失谐量为π/32,吸收从0.01 m−1增长至0.3 m−1,倍增因子从7.8减少到3.4;表明吸收系数越大,倍增因子越小。这里的吸收考虑了原子自旋碰撞因素,自旋碰撞会使气室内极化原子减少,降低原子磁力仪的探测灵敏度 [4,15],最终导致倍增因子的降低。(2)其次,出射光旋角会随着激光传输距离z的增长而逐渐减少。6条曲线都在初始阶段几乎线性增长然后缓慢趋于最大值,最后逐步减小,随着激光和碱金属原子相互作用的积累,光学厚效应和原子碰撞作用加剧,导致激光工作状态偏离无自旋交换弛豫状态,最终使得倍增因子降低。但是与吸收不同的是失谐对倍增因子的减弱作用与传输距离无关,如图6所示。例如,当吸收系数α=0.1 m−1而失谐量从π/64增长至π/32,倍增因子从11.2降低到7.8,但是两条曲线的变化趋势无明显差别。图 5 在3种条件下的激光相互作用距离与出射光旋角的关系图,从上到下分别是:α=0,Δφ=0;α=0.04 m−1,Δφ=0;α=0,Δφ=π/32
Figure 5. Relationship between the laser interaction distance and output light rotation angle in three conditions, from top to bottom are α=0, Δφ=0; α=0.04 m−1, Δφ=0; α=0, Δφ=π/32, respectively
图 6 倍增因子随着激光传输距离在两个集合下的变化:吸收系数不同时,失谐分别为Δφ=π/64和Δφ=π/32
Figure 6. Evolution of amplification factor along with the laser transmission distance under two sets: the off-resonance Δφ is π/64 and π/32 respectivitly for different absorption coefficients
图7(a)描述的是多维物理场特征,它更直观地展现了出射光旋角随吸收和失谐的变化趋势。随着吸收系数α和失谐参数Δφ的增大,光旋角ψ会非线性减少,当吸收和失谐都为0时,出射激光的光旋角最大。图7(b)中,随着吸收系数增加,光旋角逐步减少,减缓趋势与图4(a)中的理论分析结果一致。FP腔的最大倍增因子近似为16,和理论分析的结果保持一致。
Structure of Fabry-Perot cavity type weak magnetic sensitivity enhancement vapor cell
-
摘要: 提高小型化无自旋交换原子磁力仪的灵敏度是目前弱磁探测研究的难点问题,为解决这个难题提出一种基于法布里-珀罗腔的增敏气室结构。通过法布里-珀罗腔谐振原理和激光传输矩阵,在理论分析与数值计算方面对出射激光的光旋角倍增效果进行了研究。理论和仿真结果表明,随着激光在腔内往返传输次数的增加,倍增因子在初始阶段近线性增长然后缓慢趋于一个最大值,该最大值在理想状态下为16且由腔的结构参数决定。另外,碱金属原子自旋碰撞导致的吸收和腔体失谐会以不同的方式降低倍增因子,仿真结果表明:当失谐量为π/32时,倍增因子减少幅度接近50%。提出的这种腔增强型气室结构易于集成,为提高原子磁力仪的灵敏度和深入理解碱金属气室内的自旋碰撞提供了新的思路。Abstract: Enhancing the sensitivity of the miniaturized spin exchange relaxation free (SERF) atomic magnetometer for weak magnetic detection research is a difficult problem at present. A structure of sensitivity enhanced vapor cell based on Fabry-Perot cavity was proposed to solve this problem. With the principle of Fabry-Perot cavity resonance and theory of light transmission matrix, the amplification factor of optical rotation angle of emitted laser was studied in theoretical analysis and numerical simulation. The results of theoretical analysis and numerical simulation show that as the number of laser transmission round trips in cavity increases, the amplification factor of output rotation angle increases linearly in the initial stage then tends to a maximum value. Ideally, the maximum value is 16 and it is determined by the structural parameter of Fabry-Perot cavity. Besides, the absorption caused by alkali metal atoms spin collision and cavity off-resonance reduce the amplification factor in different ways. Numerical simulation results show that the reduction of amplification factor is close to 50%, while the off-resonance is π/32. This sensitivity enhanced vapor cell is easy to integrate, and provides a new perspective of sensitivity enhanced atomic magnetometer and thoroughly understanding of spin collision in alkali metal vapor cell.
-
表 1 光场数值计算参数
Table 1. Numerical calculation parameters of optical field
Parameters Value Power reflectivity(R) 90% Length of cell(l)/cm 1 Alkali number density(n)/cm−3 5×1013 Radius of electron(re) 2.8×10−13 Speed of light in vacuum(c)/m·s−1 3×108 Resonance strength(f) 1/3 FWHM of resonance(D(v))/GHz 30.4 Power of laser (P)/mW 10 -
[1] Chai G, Huang L, Qiao L, et al. Effect of the on-board residual magnetism on internal sensors [J]. Chinese Optics, 2019, 12(3): 515-525. (in Chinese) doi: 10.3788/co.20191203.0515 [2] Kominis I, Kornack T, Romalis M, et al. A subfemotesla multichannel atomic magnetometer [J]. Nature, 2003, 422: 596-599. doi: 10.1038/nature01484 [3] Wang F, Niu S, Yue C. et al. Design of attitude control system for ASTRU microsatellite [J]. Optics and Precision Engineering, 2020, 28(10): 2192-2202. (in Chinese) doi: 10.37188/OPE.20202810.2192 [4] Seltzer S. Developments in alkali-metal atomic magneto-metry[D]. Princeton: Princeton University, 2008. [5] Romalis M, Dang H, Mallof A. Ultrahigh sensitivity magnetic field and magnetization measurements with an atomic magnetometer [J]. Applied Physics Letter, 2010, 97(15): 151110. doi: 10.1063/1.3491215 [6] Wei K, Zhao T, Fang X J, et al. In-situ measurement of the density ratio of K-Rb hybrid vapor cell using spin-exchange collision mixing of the K and Rb light shifts [J]. Optics Express, 2019, 27(11): 16169-16183. doi: 10.1364/OE.27.016169 [7] Knappe S, Schwindt P, Gerginov V, et al. Microfabricated atomic clocks and magnetometers [J]. Journal of Optics A: Pure and Applied Optics, 2007, 8(7): 318-322. [8] Kornack T, Smullin S, Romalis M, et al. A low-noise ferrite magnetic shield [J]. Applied Physics Letters, 2007, 90(22): 223501. doi: 10.1063/1.2737357 [9] Fan W, Quan W, Zhang W, et al. Analysis on the magnetic field response for nuclear spin co-magnetometer operated in spin-exchange relaxation-free regime [J]. IEEE Access, 2019, 7: 28674-28580. [10] Wang Y, Liang M. Research on the spin-exchange-relaxation-free atomic magnetometer based on Herriott multipass cell [J]. Chinese Journal of Scientific Instrument, 2020, 41(6): 43-49. (in Chinese) [11] Sheng D, Li S, Dural N, et al. Sub-femtotesla scalar atomic magnetometer using multipass cell [J]. Physical Review Letters, 2013, 110(16): 160802. doi: 10.1103/PhysRevLett.110.160802 [12] Kornack T. A test of CPT and Lorentz symmetry using a K-3He co-magnetometer[D]. Princeton: Princeton University, 2005. [13] Xu G, Zhang Y, Jiang Q, et al. Temperature control of vapor cell based on the light absorption of nuclear magnetic resonance gyroscope [J]. Infrared and Laser Engineering, 2019, 48(S1): S106003. (in Chinese) [14] Ma G, Zhang J, Zhang H, et al. Resonant mode of Fabry-Perot microcavity regulated by metal surface platforms [J]. Chinese Optics, 2019, 12(3): 651-664. [15] Huang Y, Ma C, Hao Y, et al. Study on the lasing and thermal characteristic of square-FP coupled cavity semiconductor laser [J]. Chinese Journal of Lasers, 2020, 47(7): 0701017. (in Chinese) [16] Chen Q, Zhao H, Zhang W. External oil cavity coupled with EFPI partial discharge ultrasonic detection sensor [J]. Optics and Precision Engineering, 2020, 28(7): 1471-1479. (in Chinese) doi: 10.37188/OPE.20202807.1471 [17] Ning F, Tan R, Wang Y, et al. Investigation on threshold characteristic of laser-diode and end-pumped potassium vapor laser [J]. Infrared and Laser Engineering, 2019, 48(S1): S105002 . (in Chinese) [18] 吕百达. 激光光学: 激光束的传输变换和光束质量控制[M]. 成都: 四川大学出版社, 1992: 6-13. Lv B. Laser Optics: Laser Beam Propagation and Beam Quality Control[M]. Chengdu: Sichuan University Press, 1992: 6-13. (in Chinese)