Two-flat crystal three-sided mutual test technology based on N-position image rotation method
-
摘要: 光学平面面形的绝对检验技术规避了干涉仪参考面形精度的制约,能够有效提高平晶面形的检测精度。采用N次图像旋转法的两平晶三面互检的绝对检验技术,求解待测平晶的三维绝对面形分布,结果中包含了中频波段的信息。利用递推公式构造旋转变化项的N次虚拟旋转结果,求和平均后得到旋转变化项,叠加旋转不变项结果后得到待测波面面形。推导了算法的理论误差,针对旋转角度进行优化并增加虚拟旋转次数,提高了算法精度,优化后的仿真结果的残差波面的均方根值精度为0.14 nm。对150 mm口径平晶进行两平晶三面互检实验,并将实验结果与传统三面互检法结果进行比对,均方根值偏差小于0.5 nm,验证了算法的准确性。Abstract: The absolute test technology of optical surface shape avoids the limitation of reference surface shape accuracy of interferometer and can effectively improve the detection accuracy of flat-crystal surface shape. The absolute test technology of two-flat crystal three-sided mutual test of N-position image rotation method was adopted to solve the three-dimensional absolute surface shape distribution of tested flat-crystal, and the results included the information in medium frequency band. The recursion formula was used to construct the N-position virtual rotation results of rotational change item. After summing and averaging, the rotational change item was obtained, and the tested wave surface shape was obtained after superposition of the rotational invariant item. The theoretical error of the algorithm was derived, the rotation angle was optimized and the number of virtual rotation was increased, and the accuracy of the algorithm was improved. After optimization, the simulation results show that the RMS value accuracy of residual wave surface is 0.14 nm. The two-flat crystal three-sided mutual test of flat-crystal with 150 mm caliber was carried out, and the experimental results were compared with those of the traditional three-sided mutual test method. The RMS value deviation is less than 0.5 nm, which verifies the accuracy of the algorithm.
-
Keywords:
- interferometry /
- flat absolute test /
- N-position image rotation method
-
引言
超快脉冲在超快成像、光通信、光谱学、生物医学和材料加工等各个领域都有重要的应用。被动锁模光纤激光器由于具有高光束质量、高效率、高集成度、高可靠性等优点,已被广泛用于产生飞秒级超短脉冲[1]。对于被动锁模光纤激光器,锁模器件在产生超快脉冲激光器中起着关键作用。近年来,多种纳米材料都被用来作为锁模器件,例如碳纳米管[2]、石墨烯[3]、氧化铜[4]、黑磷[5]、拓扑绝缘体[6]等,均表现出3阶非线性效应强、响应时间快、工作波长范围宽的半导体特性,可作为具有竞争力的真实可饱和吸收体应用。特别是碳纳米管可饱和吸收体,具有恢复快、饱和强度低、偏振不敏感、调制深度深、工作带宽宽、环境稳定性好、制造成本低等优点,在锁模光纤激光器中得到了广泛的研究和应用。
相比单一波长的光纤激光器,双波长锁模光纤激光器在光纤通信系统、光纤传感、太赫兹波产生和相干拉曼成像等各种领域中都有广泛潜在应用[7-12]。通过调整掺铒光纤激光器结构,可以产生双波长激光的输出。双波长激光器输出2串具有一定重频差的激光对应其不同的中心波长。因为2个脉冲串是从同一个激光腔内产生,所以具有很强的相干性和共同消噪性,由此可以产生重复频率差很小的相干频率梳。近年来,已经发展了许多技术来实现双波长锁模光纤激光器。1991年,SCHLAGER J B等人首次使用3.4 km双折射光纤在主动锁模掺铒光纤环形激光器中得到双波长锁模脉冲[13]。2009年,南洋理工大学ZHANG H等人在全正色散光纤激光器中使用SESAM锁模,调节腔内双折射系数稳定产生双波长耗散孤子[14]。2016年,四川大学WANG S等人报道了一种基于双分路非线性偏振旋转技术的波长间距可控的双波长同步锁模光纤激光器[15]。2018年,LAU K Y等人采用石墨烯和聚甲基丙烯酸甲酯饱和吸收剂的双波长被动锁模掺铒光纤激光器,通过共用增益介质和可饱和吸收体的光纤分支连接,实现了锁模光纤激光器双波长脉冲输出[16]。2020年,河北师范大学XIN Y X等人利用锥形光纤的滤波效应设计了长腔可调谐双波长掺镱光纤激光器[17]。由此可见,研究者们在双波长光纤激光器领域已经做了大量的研究。但是,目前报道的许多腔内滤波结构由于存在模式干涉不稳定,滤波带宽及滤波深度不可调谐等问题,导致双波长锁模不稳定或双波长锁模光谱带宽窄,使得双波长激光器应用受限。
本文通过将偏振控制器和在线起偏器构成Lyot滤波器引入激光器中,对腔内增益谱进行调节,从而在单个激光器中产生双波长锁模脉冲。在锁模过程中,单壁碳纳米管起到可饱和吸收作用,稳定锁模状态下,通过调节偏振控制器改变Lyot滤波器的输出特性,最终实现了双波长锁模脉冲输出。此外,本文还对腔内滤波效应导致的双波长锁模机制进行了数值模拟,数值仿真结果和实验结果相吻合,进一步验证了所引入的Loyt滤波器对于波长复用激光器的影响。
1 实验原理
图1(a)为Lyot滤波器示意图,通常Lyot滤波器由偏振控制器(polarization controller,PC)、保偏光纤(polarization-maintained fiber,PMF)和起偏器(polarizer)组成。
实验中使用的在线起偏器(inline polarizer,ILP)由PMF和起偏器组成,ILP和PC构成的Lyot滤波器的传输特性[18]可以写为
$$ T = 1 - {\sin ^2}(2\theta ){\cos ^2}(\frac{{2{\text{π }}\Delta n}}{\lambda }{L_{{\text{PMF}}}}) $$ (1) 式中:∆n为PMF的双折射系数,大小等于慢轴和快轴的折射率差;LPMF为PMF的长度;θ为输入光的偏振方向与PMF的快轴之间的夹角,可以通过调整PC来进行调谐。因此,Lyot滤波器的自由光谱范围(free spectral range,FSR),可以表示为[18]
$$ \Delta\lambda=\lambda^2/\left(2\Delta n\cdot L_{\mathrm{PMF}}\right) $$ (2) 式中PMF的双折射系数约为4×10−4。由此可知,Lyot滤波器的FSR与PMF的双折射与长度的乘积成反比。图1(b)为Lyot滤波特性,虚线和实线是LPMF分别为0.15 m和0.20 m的透过率函数。
2 实验装置
图2为双波长锁模光纤激光器实验装置。由980 nm半导体激光器(pump laser)通过波分复用器(wavelength division multiplexer,WDM)泵浦,增益介质是色散系数为15.3 ps2/km的40 cm掺铒光纤(Er-doped fiber,EDF)[19],隔离器(isolator,ISO)保证了激光器的单向工作,并且将2个偏振控制器插入腔中以优化锁模性能。锁模器件碳纳米管可饱和吸收体(saturates absorbers,SA)是通过光沉积法制备,使用泵浦光源将0.8 nm~1.2 nm管径的单壁碳纳米管无水乙醇混合液沉积在光纤跳线上,能够使激光器稳定锁模24小时以上[19]。单模光纤(single-mode fiber,SMF)的色散系数为−22.9 ps2/km,长度约为4.9 m。此外,腔中还嵌入了具有25.4 cm长的保偏光纤(polarization-maintained fiber,PMF)尾纤的在线起偏器,在线起偏器和偏振控制器形成如图1所示的Lyot滤波器。通过调节偏振控制器来改变腔内偏振,可以对中心波长和光谱带宽进行微调,最后通过30/70光纤耦合器(output coupler,OC)输出脉冲,使腔内存在更多能量,使得激光更容易产生,同时分光比也决定了输出脉冲的平均功率。
光纤激光器腔长约为5.6 m,对应约36.64 MHz的基本重复频率,其净腔色散估计为−0.106 ps2,整个光纤激光器在负色散区域工作,实验中分别使用光谱分析仪(OSA Yokogawa,AQ6370)、示波器(Tektronix MS064)、频谱分析仪(Agilent Technologies,N9000A)实时监测锁模脉冲的波形、光谱和射频频谱。
3 实验结果与分析
由于2个增益峰之间的增益竞争,当激光腔的泵浦和增益饱和条件不完全匹配时,只能获得单波长锁模脉冲,因此,通过调节PC来改变增益谱轮廓,以实现双波长锁模操作。将泵浦功率提高到70.4 mW即可实现双波长锁模操作,泵浦功率由980 nm泵浦管的P-I曲线以及双波长锁模时泵浦控制器上显示的电流来计算确定,此时脉冲输出功率为3.12 mW,脉冲宽度约为ps量级。如图3(a)所示为激光器光谱图,光谱中心波长分别为1 534.2 nm和1 546.0 nm,中心波长间隔为11.8 nm,3 dB带宽分别为1.8 nm和1.5 nm。在双波长锁模状态下工作的脉冲序列的时间波形如图3(b)所示,测量的脉冲周期为27.29 ns,与光脉冲在5.6 m腔长中的传输时间一致。射频频谱如图3(c)所示,基本重复频率为36.64 MHz。通过使用较小分辨率带宽(resolution bandwidth,RBW)的射频频谱仪测量脉冲,可以观察到这2个脉冲序列的重复频率不同,重频差约为1.3 kHz,证实了2个脉冲序列为异步脉冲,这是由于腔内的群速度色散(group velocity dispersion,GVD)导致在2个分离的激光波长上产生了不同的往返延迟。这2个脉冲的信噪比分别为56.70 dB和56.72 dB,说明激光器工作在一个稳定的状态。
4 仿真模型与结果
为了进一步探索Lyot滤波效应对双波长锁模的影响,本文采用基于非线性薛定谔方程(NLSE)的模型[20]进行了数值仿真:
$$ \frac{{\partial A}}{{\partial {\textit{z}}}} = \frac{g}{2}A - \frac{{{\mathrm{i}}{\beta _2}}}{2}\frac{{{\partial ^2}A}}{{\partial {t^2}}} + {\mathrm{i}}\gamma |A{|^2}A + \frac{g}{{2\varOmega _g^2}}\frac{{{\partial ^2}A}}{{\partial {t^2}}} $$ (3) 式中:A为归一化的慢变脉冲包络线;β2为群速度延迟参数;γ为光纤截面的三次非线性系数;变量t和z分别为时间和传播距离;Ωg是具有高斯近似的增益谱的带宽;增益由g=g0·exp(−E/Esat)表示,其中g0为小信号增益,仅对增益光纤为非零,E为脉冲能量,Esat为增益饱和能量[21]。
所设计的掺铒光纤激光器的数值模型如图4所示。它由0.4 m掺铒光纤(EDF)、光谱滤光器(Filter)、5.2 m单模光纤(SMF)、可饱和吸收体(SA)、30∶70光纤耦合器(OC)组成,其中光纤耦合器主要用于提取激光脉冲输出。
在仿真中引入了一种可饱和吸收体模拟实验中的碳纳米管,锁模技术采用一个简单的传递函数[20]:
$$ T = {R_0} + \Delta R[1 - 1/(1 + P/{P_0})] $$ (4) 式中:R0为不饱和传输;ΔR为饱和传输;P为瞬时脉冲功率;P0为饱和功率。仿真模型参数如表1所示。模拟激光腔的总长度为5.6 m,对应的净色散度为−0.114 ps2。Er掺杂光纤的增益轮廓定义为在
1540 nm处的峰值高斯曲线,带宽为40 nm。采用从输出耦合器中提取的脉冲来表示脉冲的演化。为了实现双波长脉冲输出,仿真中引入Lyot滤波函数,滤波器的自由光谱范围取决于PMF的长度和双折射,调制深度依赖于偏振入射光与PMF快轴的角度差。当θ=45°,此时Lyot滤波器的传输函数如表1所示。表 1 激光器数值模拟中各项参数Table 1. Parameters in laser numerical simulationParameter Length/m Gain coefficient/(dB·m−1) Dispersion parameters
(GVD)/ (ps2·km−1)γ/(W·km)−1 Transmittance EDF 0.4 4.2 15 3 Gaussian,Ωg=40 nm SMF 5.2 — −23 3 — Lyot Filter — — — — $ T = 1 - {\text{co}}{{\text{s}}^2}(\dfrac{{2{\text{π }}\Delta n}}{\lambda }{L_{{\text{PMF}}}}) $ SA — — — — $ T = {R_0} + \Delta R[1 - 1/(1 + P/{P_0})] $ 当SA的P0设置为40 W,增益光纤的饱和能量设置为1 nJ,g0=4.2,LPMF=0.25,Δn=4×10−4,R0=0.7,ΔR=0.3,高斯脉冲能够演化为稳定的双波长脉冲。为了观察孤子演化的稳定性,仿真中进行了400次循环,发现双波长脉冲在200圈后已经能够稳定传输。图5(a)和5(b)为400次循环后的时域图和光谱图,由于实验中产生的输出脉冲受到光纤激光器中光纤色散分布不均匀和内部损耗等影响,与仿真结果并不完全一样,存在些许差异。2个脉冲间隔约为240 ps,光谱中心波长分别为1 534 nm和1 546 nm,中心波长间隔12 nm。可以发现,中心波长间隔与实验结果基本吻合。图5(c)为初始高斯脉冲经过400圈循环后的时域演化图。可以发现,随着循环圈数增加,脉冲间隔在发生改变,这是因为中心波长不同的2个脉冲在色散的影响下,其传播速度有差别,形成异步双波长脉冲[22]。图5(d)为脉冲光谱演化图,从光谱强度上看,可以发现,光谱在200圈以后一直持续不变,说明脉冲在循环200圈后已经能够稳定传输。
5 结论
本文研究了一种基于Lyot滤波效应的双波长锁模掺铒光纤激光器。在单个光纤激光器腔内引入带有强双折射的保偏光纤尾纤的在线起偏器,突破原有脉冲在腔内的传输方式,使脉冲在2个不同中心波长下传输和锁模。通过调节腔内偏振态,当泵浦功率调节至70.4 mW,输出中心波长为1 534.2 nm和1 546.0 nm,中心波长间隔为11.8 nm,3 dB带宽分别为1.8 nm和1.5 nm,基本重复频率为36.64 MHz,重频差为1.3 kHz的稳定双波长脉冲。此外,在仿真中引入Lyot滤波函数,实现了稳定的双波长输出,与实验结果基本吻合。波长复用光纤激光器中,Lyot滤波效应研究可以进一步拓展实现可调谐双波长生成,该研究结果在双光梳光谱测量等实际应用中具有一定的参考意义。
-
表 1 两平晶三面互检中旋转角度为54°时的各波面数据
Table 1 Wavefront data at 54 ° rotation angle of two-flat crystal three-sided mutual test
nm A B C PV值 RMS值 PV值 RMS值 PV值 RMS值 原始波面 31.980 5.437 60.349 9.453 22.160 4.310 计算波面 32.223 5.440 60.235 9.456 22.020 4.307 残差波面 1.984 0.181 2.012 0.179 2.012 0.179 表 2 N次图像旋转法的不同旋转次数下的残差波面数据
Table 2 Residual wavefront data with different rotation times
nm 旋转次数 15 18 20 24 30 36 40 45 60 72 90 旋转角度 48 20 54 75 84 50 63 64 66 35 52 PV值 2.200 2.260 2.012 1.757 1.954 1.826 1.908 1.660 1.733 1.788 1.521 RMS值 0.180 0.173 0.179 0.154 0.150 0.145 0.141 0.137 0.134 0.137 0.136 表 3 绝对检验恢复的Ф150 mm平晶水平垂直轮廓线数据对比(单位:nm)
Table 3 Comparison of horizontal and vertical profiles data of Ф150 mm flat-crystal
水平方向 垂直方向 PV RMS PV RMS 传统三平晶互检法 23.222 5.950 24.014 6.485 图像旋转两平晶法 23.560 6.416 21.782 6.154 表 4 不考虑折射率非均匀性误差时仿真过程中的各波面数据(单位:nm)
Table 4 Wavefront data without refractive index heterogeneity error
A B C PV RMS PV RMS PV RMS 原始波面 22.025 3.467 34.838 4.856 22.160 4.006 计算波面 20.965 3.578 33.275 4.943 20.739 4.070 残差波面 3.358 0.539 3.348 0.538 3.348 0.538 表 5 考虑折射率非均匀性误差时仿真过程中的各波面数据(单位:nm)
Table 5 Wavefront data with refractive index heterogeneity error
A B C PV RMS PV RMS PV RMS 原始波面 22.025 3.467 34.838 4.856 22.160 4.006 计算波面 22.771 3.839 34.439 5.192 22.143 4.520 残差波面 1.855 0.134 1.873 0.133 1.873 0.133 -
[1] SCHULZ G, SCHWIDER J, HILLER C, et al. Establishing an optical flatness[J]. Applied Optics,1971,10(4):929-934. doi: 10.1364/AO.10.000929
[2] 李宇琛, 韩森, 吴泉英, 等. 基于奇偶函数法的绝对检测实验研究[J]. 应用光学,2017,38(3):469-475. LI Yuchen, HAN Sen, WU Quanying, et al. Absolute test of flats based on even or odd functions[J]. Journal of Applied Optics,2017,38(3):469-475.
[3] KEENAN P B. Pseudo-shear interferometry[J]. SPIE,1983,429:2-7.
[4] BLOEMHOF E E. Absolute surface metrology with a phase-shifting interferometer for incommensurate transverse spatial shifts[J]. Applied Optics,2014,53(5):792-797. doi: 10.1364/AO.53.000792
[5] HUANG Y, MA J, ZHU R H. Absolute measurement of optical flat surface shape based on the conjugate differential method[J]. Optics Express,2015,23(23):29687-29697. doi: 10.1364/OE.23.029687
[6] FUJIMOTO I, NISHIMURA K, TAKATSUJI T, et al. Autonomous calibration method of the reference flat surface of an interferometer without using a standard flat surface[J]. Applied Optics,2012,51(20):4754-4767. doi: 10.1364/AO.51.004754
[7] SU D Q, MIAO E L, SUI Y X, et al. Absolute surface figure testing by shift-rotation method using Zernike polynomials[J]. Optics Letters,2012,37(15):3198-3200. doi: 10.1364/OL.37.003198
[8] SONG W H, WU F, HOU X, et al. Absolute measurements of flats with the method of shift-rotation[J]. Optical Review,2013,20(5):374-377. doi: 10.1007/s10043-013-0067-5
[9] HE Y H, LI Q, WAN D M, et al. New reconstruction algorithm for absolute shape calibration in two-flat test[J]. SPIE,2016,9684:96842N.
[10] ZHAI D, CHEN S, PENG X, et al. Absolute flat test using rotated and multi-shifted maps with relative tilt measurement[J]. Optics and Lasers in Engineering,2019,114:121-128. doi: 10.1016/j.optlaseng.2018.10.020
[11] XU C, YIN J, CHEN L. Method for absolute flatness measurement of optical surface[J]. Applied Optics,2009,48(13):2536-2541. doi: 10.1364/AO.48.002536
[12] 孙文卿, 陈磊, 何勇,等. 基于频域变换的两平晶互检求解算法[J]. 光电子激光,2012,20(3):542-546. SUN Wenqin, CHEN Lei, HE Yong, et al. Two-flat test solutions based on frequency domain transform[J]. Journal of Optoelectronics. Laser,2012,20(3):542-546.
[13] KUCHEL M. A new approach to solve the three flat problem[J]. Optik,2001,112(9):381-391. doi: 10.1078/0030-4026-00076
[14] EVANS C J, KESTNER R N. Test optics error removal[J]. Applied Optics,1996,35(7):1015-1021. doi: 10.1364/AO.35.001015
[15] VANNONI M, MOLESINI G. Iterative algorithm for three flat test[J]. Optics Express,2007,15(11):6809-6816. doi: 10.1364/OE.15.006809
[16] GUBIN V B, SHARNOV V N. Algorithm for reconstructing the shape of optical surface from the results of experimental data[J]. Soviet Journal of Optical Technology,1990,57(3):147-148.
[17] SCHWIDER J, BUROW R, ELSSNER K E. Homogeneity testing by phase sampling interferometry[J]. Applied Optics,1985,24(18):3059. doi: 10.1364/AO.24.003059