中国激光, 2020, 47 (1): 0101005, 网络出版: 2020-01-09   

基于圆对称径向双折射的径向偏振场光纤 下载: 1121次

Radially Polarized Field Fiber Based on Circularly Symmetric Radial Birefringence
作者单位
1 深圳大学物理与光电工程学院深圳市激光工程重点实验室, 广东 深圳 518060
2 武汉安扬激光技术有限责任公司, 湖北 武汉 430000
3 国防科技大学文理学院, 湖南 长沙 410073
摘要
当面向高功率径向偏振场激光的产生和传导时,在光纤内实现高的模式区分度与大的模场面积一直是核心的技术挑战。基于此,提出一种全新的径向偏振场光纤设计方案,通过在纤芯内部引入圆对称的径向分布热应力场,使得纤芯内形成径向双折射效应,有效打破常规光纤中偏振模式之间的简并,使TM01模、TE01模和HE21模的有效折射率差为10 -4量级,从而将TM01径向偏振模区分出来。同时,此类径向偏振场光纤更易实现TM01模场的大模面积设计。
Abstract
High mode discrimination and large mode area in fibers are key technical challenges for the generation and transmission of high-power radially polarized field (RPF) laser. Based on this, a novel method for the design of a radially polarized field fiber is proposed. By introducing a circularly symmetrical radially distributed thermal stress field into the core of the fiber, a radial birefringence effect can be realized in the core, which effectively breaks the degeneracy between polarization modes in conventional optical fibers. The effective index difference among TM01, TE01, and HE21 modes is of the order of 10 -4. This allows the TM01 radial polarization mode to be separated. In addition, a large mode area design for the TM01 mode field can be achieved by using the proposed RPF fiber.

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]

2Ψ=1+ν1-ναΔT,(1)

式中:Ñ为哈密顿算子;Ψ为位移势,满足ux=Ψx,uy=Ψy, uxx方向的位移,uyy方向的位移;ν为泊松比,取0.186[13]T为室温与参考温度之差,取ΔT=20 ℃-1000 ℃=-980 ℃[13];α为热膨胀系数,是关于位置坐标的函数,可以通过掺杂来调控石英玻璃的热膨胀系数。掺杂石英玻璃的热膨胀系数计算公式为[13]

α=(1-m)αSiO2+mαdopant,(2)

式中:αSiO2为纯石英玻璃的热膨胀系数,取值5.4×10-7-1;αdopant为掺杂物质的热膨胀系数,主要涉及GeO2的热膨胀系数(7×10-6-1)和B2O3的热膨胀系数(10×10-6-1);m为掺杂物质的摩尔分数(单位为%,表示掺杂物分子的物质的量占总物质的物质的量的百分比)。由此可得x方向的应变εxy方向的应变εy之和为

εx+εy=uxx+uyy=2Ψx2+2Ψy2=1+ν1-ναΔT(3)

根据物理方程[14]

σz=E(1+ν)(1-2ν)νεx+εy)-(1+ν)αΔT]+p,(4)

式中:E为弹性模量,取78×10-9 m2· N-113;p为修正项,以确保∬σzdxdy=0[15]。由此可得z方向的应力分量为

σz=-E1-ναΔT+p,(5)

其中修正项为

p=E1-νΔTi=1nαi(ri2-ri-12)/rn2,(6)

式中:αi为从内到外第i个圆环或圆域的热膨胀系数,即α在特定位置的值;ri为从内向外第i个圆的半径,规定r0=0。

根据物理方程[15]和平衡方程[16],

σz=ν(σr+σθ)-EαΔT+p*,(7)

σrr+1rτθrθ+τzrz+σr-σθr=0,(8)

式中:r为径向坐标;θ为方位角坐标;p*为修正项;σr和σθ分别为径向和角向应力分量。由于该系统是圆对称的系统,径向即主应力方向,所以无剪切应力项τθr和τzr。联立(7)式和(8)式,消去角向应力分量,可得

r2σrr2+3σrr=0(9)

该微分方程通解为

σr=a-br2,(10)

式中:a和b为待定系数。将(10)式代入(8)式,整理得

σθ=a+br2(11)

根据条件[15]:1)纤芯中心的应力是有限的;2)σr在不同介质的交界处连续;3)σr在光纤最外表面处为零,因为无外应力,满足

p*=p(1-ν),(12)

a=12σz,(13)

b=0,rr112·E1-νΔT(α1-α2)r12,r1<rr212·E1-νΔT[(α1-α2)r12+(α2-α3)r22],r2<rr3             12·E1-νΔTi=1n-1(αi-αi+1)ri2,rn-1<rrn(14)

若考虑渐变掺杂型光纤,则热膨胀系数不再局限于阶跃型分布,理论上可以是任意合理的函数分布,可以根据微元法及以上结论推得,形式上大体相同,只需要进行参数修正。

p修正为

p~=E1-νΔT0rfiber2α(r)rdr/rfiber2,(15)

式中:rfiber为整个光纤的半径(不考虑高分子涂覆层)。

b修正为

b~=12E1-νΔT0r2α(R)RdR-α(r)r2,(16)

式中:R为积分变量,表示径向坐标。

2.2 RPF光纤的应力双折射

根据光弹效应,光纤的折射率与应力的关系式为[13]

nr=n0-[B1σr+B2(σθ+σz)]nθ=n0-[B1σθ+B2(σz+σr)]nz=n0-[B1σz+B2(σr+σθ)],(17)

式中:B1B2为光弹系数,分别取0.7572448×10-12 m2·N-1和4.18775×10-12 m2·N-1;n0为光纤在无应力状态下的折射率分布;nrnθnz分别是光纤在应力状态下径向、角向和轴向的折射率分布。则阶跃型光纤关于径向和角向的双折射B和渐变掺杂型光纤关于径向和角向的双折射B~分别为

B=nr-nθ=(B1-B2)ΔTE1-ν0,rr1(α1-α2)r12/r2,r1<rr2(α1-α2)r12+(α2-α3)r22]/r2,r2<rr3           i=1n-1(αi-αi+1)ri2/r2,rn-1<rrn,(18)

B~=nr-nθ=(B1-B2)ΔTE1-ν0r2α(R)RdR/r2-α(r)(19)

关于光纤在无应力状态下的折射率分布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光纤的数值仿真方法

根据以上计算,光纤的相对介电常数张量为

ε~=nr2000nθ2000nz2(20)

以上计算都是在柱面坐标系下完成的,但有限元法数值仿真需要在直角坐标下进行,所以需要进行坐标系的变换。在直角坐标下的相对介电常数张量为

ε-=T-1ε~T,(21)

式中:T为变换矩阵,具体表达式为

T=cosθsinθ0-sinθcosθ0001(22)

位置坐标的变换为

r=x2+y2sinθ=yx2+y2z=z(23)

用有限元法对光纤进行数值仿真,将(21)式代入计算模型中计算分析光纤的波导模式,可以获得模场分布、模式有效折射率等信息。配合循环结构,可以进行相应的扫描计算,便于光纤设计。

3 RPF光纤设计

根据上述理论,提出一种全新的设计方案,纤芯包含core1和core2两个区域,如图1所示。core1中掺杂B2O3和GeO2,其中B2O3具有大的热膨胀系数,用于产生圆对称的径向分布热应力场,而GeO2主要用于调整折射率分布,从而调整模场分布,使应力双折射分布与LP11模场分布的重叠区增大。core2中掺杂GeO2用于传导RPF。包层为纯石英玻璃。

图 1. RPF光纤结构

Fig. 1. Structure of RPF fiber

下载图片 查看所有图片

图2是以图1的结构为基础,考察掺杂物质浓度和纤芯结构对模式退简并的影响,图中每一条曲线都是通过用控制变量法单独研究其中一个变量得到。在图2中,xyz三个箭头标明三个坐标轴的方向,其中x轴表示core2的直径(整个纤芯的直径,可在拉丝中改变光纤直径来控制,且core2的直径与包层直径具有固定比例1∶10),z轴表示模式间的有效折射率差Δneff,而y轴是主要考察变量。其中图2(a)考察的是core1中GeO2的摩尔分数(y轴)对模式退简并的影响,GeO2主要用于调整core1区域的折射率,从而调控模场分布;图2(b)考察的是纤芯结构对模式退简并的影响,纤芯结构被定量为core1与core2 的直径之比D'(y轴),纤芯结构的变化既可调控模场分布又可调控双折射分布。图2中方形标记曲线为HE11模与TM01模的有效折射率差,圆形标记曲线为TM01模和HE21模的有效折射率差,三角形标记曲线为HE21模与TE01模的有效折射率差。图2中点划线网格表示Δneff为1×10-4的平面,曲线在平面以上部分即为满足模式退简并的标准。预先设定core1中B2O3的摩尔分数为17.5%,用于提供较大的热应力,预先设定core2中GeO2的摩尔分数为1.18%,初步确定光纤的NA水平,RPF光纤内GeO2的摩尔分数相比于高折射率差环形纤芯光纤已是大大降低,因此其NA远低于高折射率差环形纤芯光纤。此外图2(a)预先设定D'为0.56,图2(b)预先设定core1中GeO2的摩尔分数为7%,通过图2发现,当core1中GeO2的摩尔分数被设置为7%,即图2(a)虚线所标位置,且D'为0.56时,即图2(b)虚线所标位置是最优化的取值,TM01-HE21的有效折射率差都接近极大值。此条件下,core2的直径即便在15~25 μm大范围内取值时,各模式间的有效折射率差都大于1×10-4,因此在本方案中支持实现TM01模的大模面积设计。

表 2. 本设计与几种高折射率差环形纤芯光纤的对比

Table 2. Comparison of our design with several ring-core fibers with high refractive index difference

ParameterMethod inRef. [10]Method inRef. [19]Method inRef. [20]Method inRef. [21]Our design
Inner/outer diameter /μm~7/~10~8/~914/18~6/~16.511.2/20
Refractive index difference~0.02300.18000.0300~0.0350~0.0017

查看所有表

表 1. 几个典型波长的有效模场面积和模式间的有效折射率差

Table 1. Effective mode field area and effective refractive index difference between modes at several typical wavelengths

Wavelength /nmEffective mode field area /μm2Effective refractive index difference /10-4
TM01HE21TE01HE11-TM01TM01-HE21HE21-TE01
7302982883101.011.151.64
10303693703862.971.201.45
10643793823983.211.201.43
155070279410436.231.021.04
155773284011726.331.001.01

查看所有表

图 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图4表1所示。

图3为在1030 nm处RPF光纤的折射率分布和圆对称热应力场所导致的应力双折射分布。其中,图3(a)为光纤的折射率分布,实线为径向折射率分布nr,虚线为角向折射率分布nθ,点划线为轴向折射率分布nz;图3(b)为TM01模的模场分布(虚线)与双折射分布(实线,即nr-nθ)。可见TM01模的模场分布与介质双折射分布的重叠区大,且位于大的双折射值范围之内,显然有利于TM01模与其他模式的退简并。

图 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

下载图片 查看所有图片

图4(a)为光纤在1 μm波段的模式有效折射率差,各模式间的有效折射率差都满足大于1×10-4的要求。在1030 nm处的模式分布如图4(b)所示,图中箭头表示电场矢量分布,其中TM01模、HE21模和TE01模的有效模场面积分别为369,370,386 μm2

在不同波长处,TM01等模式与圆对称应力双折射区的重叠会发生相应变化,并影响到模式间的退简并。仍以相邻模式间折射率差大于1×10-4为基本要求,计算了不同波长处模式间的有效折射率差以及有效模场面积,如表1所示。对此参数下的光纤在730~1570 nm波长范围内其相邻模式有效折射率差都大于1×10-4。相比文献[ 10]中环形纤芯直径约为10 μm,在通信波段其有效模场面积约为100 μm2,本设计中RPF光纤的TM01模场面积约为700 μm2,提升了7倍左右。此外,表2中第1行数据显示了几种高折射率差环形纤芯的内外直径,外直径普遍不超过20 μm;第2行数据显示了上述高折射率环形纤芯的折射率与包层的折射率之差,折射率差普遍远高于10-3量级。以上结果直观表明,采用圆对称应力双折射的方案,是实现TM01模场面积大幅度增加且能够控制模式数量的一种有效方案。

图 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中模场分布得到,通过设计可将基模的主要能量限制于无源的core1中心区域,从而降低其增益,并使其模式竞争能力大幅小于TM01模,类似于设计大模场单模光纤的限制掺杂方法[22]。通过弯曲光纤或进一步降低纤芯NA等手段,可将TE01模和HE21模充分抑制,该光纤可有效激发TM01模以实现大功率高纯度的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.

[3] Kim W C, Park N C, Yoon Y J, et al. Investigation of near-field imaging characteristics of radial polarization for application to optical data storage[J]. Optical Review, 2007, 14(4): 236-242.

[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.

[5] Yan J, Lu Y H, Wang P, et al. Improving the sensitivity of fiber-optic SPR sensor via radially polarized beam excitation[J]. Chinese Optics Letters, 2009, 7(10): 909-911.

[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.

    Wu Y L, Zheng Y, Li S Y, et al. Propagation properties analysis of fused taper PANDA fiber couplers[J]. Acta Optica Sinica, 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.

本文已被 4 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!