基于圆对称径向双折射的径向偏振场光纤 下载: 1121次
1 引言
径向偏振光是一种典型的圆柱矢量光束,具有轴对称偏振结构,其光场中心存在偏振“奇点”,光斑的中心光强为零。相比线偏振光,径向偏振光在聚焦情况下产生的焦点场呈螺旋对称不变形,可应用于粒子捕获[1]、粒子加速[2]、光学存储[3]、光学显微[4]、传感器[5]和激光加工[6-7]等领域,尤其在激光加工中,径向偏振光的聚焦焦点质量好,可有效提高加工质量和效率。少模光纤内,常见模式除了HE11模式,还有TM01模、HE21模和TE01模。其中TM01模为径向偏振场(RPF)分布,HE21模为混合偏振场分布,TE01模为角向偏振场分布。一般少模光纤中,HE21、TM01和TE01模式之间的有效折射率差在10-5量级,这种低的模式区分度使模式串扰并简并为LP11模,不利于RPF光的产生以及长距离稳定传导。已报导工作中有借助单模-少模光耦合器[8]、长周期光纤光栅[9]实现光纤内不同模式转换,以获得RPF。为打破LP11模中各矢量模的简并,2009年Ramachandran等[10]采用了一种新型的纤芯结构设计,在阶跃型纤芯之外额外增加一个更高折射率的环形区域,TM01模与其他模式的有效折射率差达到了1.8×10-4,再借助长周期光栅获得TM01模。然而,单模-少模光耦合器、长周期光纤光栅、以及大折射率差环形纤芯光纤都是在低功率应用场景下通过模式转换来实现RPF。从方案本身思考,单模-少模光耦合器与长周期光纤光栅都难以承受高功率下RPF激光的转换,高折射率差环形纤芯光纤的模场面积小,亦不利于高功率RPF激光传导。
当光纤应用于高功率RPF激光的传导与产生时,如何同时确保高的模式区分度和大的模场面积,是当前本领域内最关切的技术挑战之一。高的模式区分度可降低不同模式间的串扰,从而让RPF光束保持更高的模式纯度;大的TM01模场面积能降低高功率RPF激光运转和传导时的非线性效应,尤其对高峰值功率的超快激光脉冲的产生和传导具有重要意义。但遗憾的是,采用高折射率差环形纤芯方案不利于实现TM01大的模场面积设计,模式数量难以控制,模式纯度难以保证,而常规少模增益光纤与相位波片结合的方式需要兼顾相位波片的功率承受阈值,该方法通常在主放大级部分牺牲了RPF的模式纯度[11]。
本文提出一种全新的RPF光纤设计思路:基于圆对称的热应力双折射,创造性地在纤芯中心掺杂高热膨胀系数的B2O3,引入一个径向分布的应力场。该设计将会使RPF和角向偏振场产生模式双折射,增大相互的有效折射率差,从而打破TM01模、TE01模和HE21模之间的简并。本设计可用于实现传导高功率的低数值孔径(NA)大模场无源光纤或放大RPF的低NA大模场有源光纤。本文计算圆对称光纤(包括阶跃型光纤和渐变掺杂型光纤)的热应力分布,详细地研究了这种光纤关于径向和角向的双折射特性,获得的TM01模、TE01模和HE21模之间的有效折射率差为1×10-4量级,同时其有效模场面积较高折射率差环形纤芯光纤大幅增加。
2 基本理论
2.1 RPF光纤的热应力分布
假设光纤由n个同心圆构成,每个圆环或圆域内的热膨胀系数均匀分布,每个区域之间的热膨胀系数不尽相同,则该光纤在柱面坐标系下的热应力分布如下。
假设光纤参量不随轴向变化,则光纤为平面应变模型,轴向无应变。根据方程[12]
式中:Ñ为哈密顿算子;Ψ为位移势,满足ux=
式中:αSiO2为纯石英玻璃的热膨胀系数,取值5.4×10-7 ℃-1;αdopant为掺杂物质的热膨胀系数,主要涉及GeO2的热膨胀系数(7×10-6 ℃-1)和B2O3的热膨胀系数(10×10-6 ℃-1);m为掺杂物质的摩尔分数(单位为%,表示掺杂物分子的物质的量占总物质的物质的量的百分比)。由此可得x方向的应变εx和y方向的应变εy之和为
根据物理方程[14]
式中:E为弹性模量,取78×10-9 m2·
其中修正项为
式中:αi为从内到外第i个圆环或圆域的热膨胀系数,即α在特定位置的值;ri为从内向外第i个圆的半径,规定r0=0。
式中:r为径向坐标;θ为方位角坐标;p*为修正项;σr和σθ分别为径向和角向应力分量。由于该系统是圆对称的系统,径向即主应力方向,所以无剪切应力项τθr和τzr。联立(7)式和(8)式,消去角向应力分量,可得
该微分方程通解为
式中:a和b为待定系数。将(10)式代入(8)式,整理得
根据条件[15]:1)纤芯中心的应力是有限的;2)σr在不同介质的交界处连续;3)σr在光纤最外表面处为零,因为无外应力,满足
若考虑渐变掺杂型光纤,则热膨胀系数不再局限于阶跃型分布,理论上可以是任意合理的函数分布,可以根据微元法及以上结论推得,形式上大体相同,只需要进行参数修正。
p修正为
式中:rfiber为整个光纤的半径(不考虑高分子涂覆层)。
b修正为
式中:R为积分变量,表示径向坐标。
2.2 RPF光纤的应力双折射
根据光弹效应,光纤的折射率与应力的关系式为[13]
式中:B1和B2为光弹系数,分别取0.7572448×10-12 m2·N-1和4.18775×10-12 m2·N-1;n0为光纤在无应力状态下的折射率分布;nr、nθ和nz分别是光纤在应力状态下径向、角向和轴向的折射率分布。则阶跃型光纤关于径向和角向的双折射B和渐变掺杂型光纤关于径向和角向的双折射B~分别为
关于光纤在无应力状态下的折射率分布n0,纯石英玻璃的折射率可根据经典色散公式计算[17]。此外可以通过掺杂来改变石英玻璃的折射率,掺杂GeO2和B2O3引起的折射率变化与掺杂的摩尔分数成线性关系[18],每掺杂摩尔分数为1%的GeO2和B2O3,引起的折射率变化大约分别为+14.26×10-4和-4.2×10-4,由此可估算不同掺杂情况下的石英玻璃折射率。
从(18)式可以得出,对于阶跃型光纤:1)光纤中心没有双折射,如果纤芯只有一层结构则不能有效产生双折射,纤芯结构至少需要两层;2)应力双折射会以1/r2的衰减速度沿径向衰减,需合理设计热膨胀系数分布和折射率分布,确保应力双折射与LP11模场的分布重叠区很大,才能有效产生模式双折射(>1×10-4)。
2.3 RPF光纤的数值仿真方法
根据以上计算,光纤的相对介电常数张量为
以上计算都是在柱面坐标系下完成的,但有限元法数值仿真需要在直角坐标下进行,所以需要进行坐标系的变换。在直角坐标下的相对介电常数张量为
式中:T为变换矩阵,具体表达式为
位置坐标的变换为
用有限元法对光纤进行数值仿真,将(21)式代入计算模型中计算分析光纤的波导模式,可以获得模场分布、模式有效折射率等信息。配合循环结构,可以进行相应的扫描计算,便于光纤设计。
3 RPF光纤设计
根据上述理论,提出一种全新的设计方案,纤芯包含core1和core2两个区域,如
表 2. 本设计与几种高折射率差环形纤芯光纤的对比
Table 2. Comparison of our design with several ring-core fibers with high refractive index difference
|
表 1. 几个典型波长的有效模场面积和模式间的有效折射率差
Table 1. Effective mode field area and effective refractive index difference between modes at several typical wavelengths
|
图 2. 光纤参数对模式退简并的影响。(a) core1中GeO2浓度对模式退简并的影响;(b)纤芯结构对模式退简并的影响
Fig. 2. Effects of fiber parameters on mode degeneracy. (a) Effect of GeO2 concentration in core1 on ode degeneracy; (b) effect of fiber core structure on mode degeneracy
依据上述最优化设计条件,将RPF光纤的参数选取为:core1的直径为11.2 μm,core2的直径为20 μm(即D'为0.56),包层直径为200 μm;core1中B2O3、GeO2、SiO2的摩尔分数比例为17.5∶7∶75.5,core2中掺GeO2的摩尔分数为1.18%。下面对此参数下RPF光纤进行计算以更直观地获得光纤的特性,计算结果如
图 3. RPF光纤的各向异性折射率特性。(a) 1030 nm处RPF光纤在3个方向上的折射率分布;(b) 1030 nm处TM01模场分布与双折射分布
Fig. 3. Anisotropic refractive index properties of RPF fiber. (a) Refractive index profiles of RPF fiber in three irections at 1030 nm; (b) TM01 mode field distribution at 1030 nm and birefringence distribution
在不同波长处,TM01等模式与圆对称应力双折射区的重叠会发生相应变化,并影响到模式间的退简并。仍以相邻模式间折射率差大于1×10-4为基本要求,计算了不同波长处模式间的有效折射率差以及有效模场面积,如
图 4. RPF光纤的模式退简并情况与模场分布。(a)光纤模式之间有效折射率差随波长的变化;(b) 1030 nm处的模式分布
Fig. 4. Mode degeneracy and mode field distribution of RPF fiber. (a) Effective index difference between fiber odes versus wavelength; (b) mode distributions at 1030 nm
更进一步,core2可改为有源离子掺杂,用于RPF激光的增益放大。从
4 结论
提出一种全新的RPF光纤的方案,证明了基于圆对称的热应力场效应可打破LP11模的模式简并,相邻模式有效折射率差大于1×10-4,可降低模式之间的串扰。本设计中RPF光纤的纤芯折射率大大降低,与高折射率差环形纤芯光纤相比,纤芯与包层的折射率差可降低1个数量级。此类RPF光纤可用于设计低NA大模场无源或有源光纤,其TM01有效模场面积比一般的高折射率差环形纤芯光纤大,在高功率RPF激光的产生与传导中具有潜在的应用价值。
[1] Zhao Y Q, Zhan Q W, Zhang Y L, et al. Creation of a three-dimensional optical chain for controllable particle delivery[J]. Optics Letters, 2005, 30(8): 848-850.
[2] Salamin Y I. Electron acceleration from rest in vacuum by an axicon Gaussian laser beam[J]. Physical Review A, 2006, 73(4): 043402.
[4] KozawaY, SatoS. Enhanced detection of longitudinal field of a radially polarized beam in confocal laser microscopy[C]∥CLEO: 2015, May 10-15, 2015, San Jose, California, United States. Washington, D.C.: OSA, 2015: STh3K. 2.
[6] Venkatakrishnan K, Tan B. Interconnect microvia drilling with a radially polarized laser beam[J]. Journal of Micromechanics and Microengineering, 2006, 16(12): 2603-2607.
[7] Nesterov V, Niziev G. Laser beams with axially symmetric polarization[J]. Journal of Physics D: Applied Physics, 2000, 33(15): 1817-1822.
[8] Huang Y P, Shi F, Wang T, et al. High-order mode Yb-doped fiber lasers based on mode-selective couplers[J]. Optics Express, 2018, 26(15): 19171-19181.
[9] Zhao Y H, Liu Y Q, Zhang L, et al. Mode converter based on the long-period fiber gratings written in the two-mode fiber[J]. Optics Express, 2016, 24(6): 6186-6195.
[10] Ramachandran S, Kristensen P, Yan M F. Generation and propagation of radially polarized beams in optical fibers[J]. Optics Letters, 2009, 34(16): 2525-2527.
[11] Lin D, Baktash N, Alam S U, et al. 106 W, picosecond Yb-doped fiber MOPA system with a radially polarized output beam[J]. Optics Letters, 2018, 43(20): 4957-4960.
[12] Lu D, Ge T W, Wu J, et al. Thermal stress induced birefringence in double cladding fiber with non-circular inner cladding[J]. Journal of Modern Optics, 2009, 56(5): 638-645.
[13] Guan R F, Zhu F L, Gan Z Y, et al. Stress birefringence analysis of polarization maintaining optical fibers[J]. Optical Fiber Technology, 2005, 11(3): 240-254.
[14] 吴宇列, 郑煜, 李圣怡, 等. 熔锥型熊猫光纤耦合器的传输特性分析[J]. 光学学报, 2007, 27(6): 1111-1115.
[15] Sakai J, Kimura T. Birefringence caused by thermal stress in elliptically deformed core optical fibers[J]. IEEE Journal of Quantum Electronics, 1982, 18(11): 1899-1909.
[16] 陈昌麒, 田世兴, 陈森灿, 等. 材料学科中的固体力学[M]. 北京: 北京航空航天大学出版社, 1994: 39- 41.
Chen CQ, Tian SX, Chen SC, et al.Solid mechanics in materials science[M]. Beijing: Beihang University Press, 1994: 39- 41.
[17] Malitson I H. Interspecimen comparison of the refractive index of fused silica[J]. Journal of the Optical Society of America, 1965, 55(10): 1205-1209.
[18] 童维军. 熊猫保偏光纤结构设计、制备工艺及性能研究[D]. 武汉: 华中科技大学, 2006: 90.
Tong WJ. Structure design, fabrication process and property research of PANDA PMF[D]. Wuhan: Huazhong University of Science and Technology, 2006: 90.
[19] Yan Y, Zhang L, Wang J, et al. Fiber structure to convert a Gaussian beam to higher-order optical orbital angular momentum modes[J]. Optics Letters, 2012, 37(16): 3294-3296.
[20] PidishetyS, SrinivasP, SrinivasanB. Studies on the propagation of orbital angular momentum beams through optical fibers with different index profiles[C]∥2015 Workshop on Recent Advances in Photonics (WRAP), December 16-17, 2015, Bangalore, India. New York: IEEE, 2015: 16596551.
[21] Pidishety S, Srinivasan B, Brambilla G. All-fiber fused coupler for stable generation of radially and azimuthally polarized beams[J]. IEEE Photonics Technology Letters, 2017, 29(1): 31-34.
[22] Marciante J R. Gain filtering for single-spatial-mode operation of large-mode-area fiber amplifiers[J]. IEEE Journal of Selected Topics in Quantum Electronics, 2009, 15(1): 30-36.
朱方祥, 陈浩, 尹金德, 王金涛, 闫培光, 韦会峰, 何飞, 陈抗抗, 于洋. 基于圆对称径向双折射的径向偏振场光纤[J]. 中国激光, 2020, 47(1): 0101005. Fangxiang Zhu, Hao Chen, Jinde Yin, Jintao Wang, Peiguang Yan, Huifeng Wei, fei He, Kangkang Chen, Yang Yu. Radially Polarized Field Fiber Based on Circularly Symmetric Radial Birefringence[J]. Chinese Journal of Lasers, 2020, 47(1): 0101005.