雷前召， 何 冰， 董康軍， 張修興， 張寧寧， 吳振森

(1. 渭南師范學(xué)院 物理與電氣工程學(xué)院， 陜西 渭南 714099； 2. 西安電子科技大學(xué) 物理與光電工程學(xué)院， 陜西 西安 710071)

光纖通信的應(yīng)用已經(jīng)滲透到社會(huì)、生活的各個(gè)方面，從傳輸簡單信號(hào)到影像快速傳遞技術(shù)，從云主機(jī)、云服務(wù)器的普及應(yīng)用到數(shù)據(jù)安全預(yù)警[1]等。新型冠狀病毒感染的肺炎(新冠肺炎，novel coronavirus pneumonia，NCP)疫情暴發(fā)，更進(jìn)一步促進(jìn)了線上購物、線上教育、線上辦公以及線上醫(yī)療對(duì)傳統(tǒng)方式的加速替代。隨著世界快速進(jìn)入5G時(shí)代，社會(huì)的生產(chǎn)、生活方式，企業(yè)、國家的組織結(jié)構(gòu)、管理方法、模式都將發(fā)生深刻變化。

光纖通信系統(tǒng)容量巨大，優(yōu)勢明顯，在引領(lǐng)各國搶占科技領(lǐng)域前沿方面具有巨大潛力。目前常用的二氧化硅光纖在光纖通信中采用的低損耗工作波長為近紅外線波長，分別為1.31、1.55、1.62 μm，相比于頻率僅為幾十兆赫茲的可用通信帶寬的短波、微波，光纖通信在這3個(gè)低損耗窗口的對(duì)應(yīng)波段各自具有頻率分別為17.5、4.3、7.1 THz的帶寬儲(chǔ)備，如此豐富的帶寬資源有力地保障了光纖的大容量通信能力[2]。光纖的特殊材料決定了其具有優(yōu)良的抗干擾性[3]，即使太陽黑子暴發(fā)，電離層出現(xiàn)磁暴，光纖通信仍然能穩(wěn)定、正常運(yùn)行。此外，當(dāng)光纖與高壓輸電線平行架設(shè)或與電力導(dǎo)體構(gòu)成復(fù)合光纜時(shí)，能夠最大程度地提高電力傳輸線路的運(yùn)行效率[4]。

波導(dǎo)光學(xué)主要研究光在光纖中的產(chǎn)生、耦合、傳播和損耗等物理現(xiàn)象和性質(zhì)[5]，以及與之相應(yīng)的光波導(dǎo)器件，比如光纖耦合器、光纖傳感器[6-7]、光纖波分復(fù)用器、光開關(guān)、光放大器以及光調(diào)制器等[8]。集成光路主要研究如何把各種光波導(dǎo)器件集成在同一根光纖中，形成具有一定功能的微型光纖光學(xué)系統(tǒng)，以達(dá)到纖維集成的目的[9]。上述所有問題都離不開光纖波導(dǎo)的模式傳光研究。本文中將光纖波導(dǎo)內(nèi)矢量場進(jìn)行簡并處理，轉(zhuǎn)化為波導(dǎo)內(nèi)的線偏振模式，求解柱坐標(biāo)系下標(biāo)量場的波動(dòng)方程，得到線偏振模式的貝塞爾特征方程，并利用線偏振模式的貝塞爾函數(shù)解，模擬光纖波導(dǎo)橫截面場強(qiáng)分布，分析場強(qiáng)分布特征。

1 光纖波導(dǎo)傳光基本原理

光線在光纖中均勻同種介質(zhì)的纖芯內(nèi)沿直線傳播，在內(nèi)、外2層介質(zhì)分界面發(fā)生反射和折射。光纖結(jié)構(gòu)及其數(shù)值孔徑如圖1所示，其中，a、b為纖芯、包層的半徑，φ為入射光線與光纖軸線的夾角。

入射光線能夠通過光纖傳輸?shù)淖畲笕肷浣堑恼抑刀x為光纖的數(shù)值孔徑rNA，以2φmax為頂角的圓錐體內(nèi)的光線能通過光纖，根據(jù)光纖內(nèi)纖芯與包層間光線全反射的臨界條件，有

a、b—纖芯、包層的半徑； φ—入射光線與光纖軸線的夾角。圖1 光纖結(jié)構(gòu)及其數(shù)值孔徑

(1)

式中n1、n2、n0分別為光纖纖芯、包層及包層外介質(zhì)的折射率。包層外通常為空氣，光線在纖芯內(nèi)多次反射形成折線路徑，光線如果與光纖軸線相交，則入射線、反射線、法線、軸線共面，該面稱為主截面，其中光線稱為子午線，光纖子午線投影圖如圖2所示。通常，光線入射角φ滿足

(2)

φ、θ—光線通過纖芯橫截面進(jìn)入纖芯的入射角和反射角; α—纖芯內(nèi)光線在纖芯和包層分界面的入射角; n1、n2—光纖纖芯、包層的折射率；z—空間直角坐標(biāo)系的豎坐標(biāo)。圖2 光纖子午線投影圖

(3)

2 光纖傳光標(biāo)量波導(dǎo)方程及貝塞爾方程解

2.1 光纖傳光標(biāo)量波導(dǎo)方程

光纖波導(dǎo)除了可以傳輸電場矢量垂直于光線傳播方向的橫電(transverse electric，TE)波(或H波)、磁場矢量垂直于傳播方向的橫磁(transverse magnetic，TM)波(或E波)這2種線偏振波，還有傳播方向既有磁場分量又有電場分量的HE波或EH波(兼有電場、磁場矢量垂直于光線傳播方向)，直接求解TE、TM、HE、EH模的矢量精確解很復(fù)雜。通常光纖采用弱導(dǎo)光纖(纖芯、包層的折射率滿足n1≈n2)，在截止頻率附近相當(dāng)寬的頻帶內(nèi)，波的傳播系數(shù)近似相等，極化方向不隨z變化，可視為線極化波，問題轉(zhuǎn)化為求解標(biāo)量方程，則可引用線偏振模式(linearly polarized modes，LP模式)[11]。選用柱面坐標(biāo)系，使z軸與光纖軸線重合，并且視場分量A滿足標(biāo)量波動(dòng)方程，則有關(guān)系式[12]

(4)

式中：μ、ε分別為介質(zhì)的磁導(dǎo)率和電容量；r、φ、z、t分別為柱面坐標(biāo)系的3個(gè)分量以及時(shí)間因子；A(r,φ,z,t)=R(r)φ(φ)Z(z)T(t)為利用分離變量法所得的解，其中R、φ、Z、T分別為只含有r、φ、z、t的代數(shù)因子，φ(φ)=exp(±lφ)，l(l=0, 1, …)為圍繞z軸旋轉(zhuǎn)的場分量變化的周期數(shù)，Z(z)T(t)=exp[j(ωt-βz)]，ω為角頻率。

將A=R(r)φ(φ)Z(z)T(t)代入式(4)，得

(5)

式中：k為波傳播的波數(shù)；i分別取為1、2。式(5)的解為纖芯和包層中的徑向場分量[13]，即

(6)

式中：a0、ba為待定系數(shù)；V為歸一化頻率；W為貝塞爾函數(shù)的宗量；Jl為纖芯內(nèi)貝塞爾函數(shù)；Kl為包層內(nèi)貝塞爾函數(shù)。

(7)

其中

并且參數(shù)β、V、U、W隨l、ω、a等的變化而變化。

2.2 標(biāo)量近似解與傳統(tǒng)矢量精確解相對(duì)應(yīng)

根據(jù)式(7)，當(dāng)W=0時(shí)，特征方程為臨界狀態(tài)，此時(shí)V=U，由此確定各模式截止頻率。由于U≠0，因此有

Jl-1(U)=0

。

(8)

當(dāng)l=0時(shí)，J-1(U)=J1(U)=0的第1個(gè)零點(diǎn)出現(xiàn)在U=0處，傳播模式在U=0處截止，該傳播模式為LP01； 第2個(gè)零點(diǎn)出現(xiàn)在V=3.832處,傳播模式為LP02； 第3個(gè)零點(diǎn)在7.016處，傳播模式為LP03； 第4個(gè)零點(diǎn)出現(xiàn)在9.761處，傳播模式為LP04。

當(dāng)l=1時(shí)，J0(U)=0的第1、2、3個(gè)零點(diǎn)分別出現(xiàn)在U=2.405、5.52、8.654處，對(duì)應(yīng)的模式分別為LP11、LP12、LP13。

當(dāng)l=2時(shí)，J1(U)=0的第1、2、3個(gè)零點(diǎn)分別出現(xiàn)在U=3.832、7.016、9.761處，對(duì)應(yīng)的模式分別為LP21、LP22、LP23。

當(dāng)l=3時(shí)，J2(U)=0的第1、2個(gè)零點(diǎn)分別出現(xiàn)在U=5.136、8.417處，對(duì)應(yīng)的模式分別為LP31、LP32。

當(dāng)l=4時(shí)，J3(U)=0的第1、2個(gè)零點(diǎn)分別出現(xiàn)在U=7.588、8.711處，對(duì)應(yīng)的模式分別為LP41、LP42。

統(tǒng)一使用LPln表示LP模式，其中n為貝賽爾函數(shù)的零點(diǎn)數(shù)，實(shí)際含義為場徑向分量變化的半波數(shù)[13]。

LPln模式遠(yuǎn)離截止時(shí)，則對(duì)應(yīng)的頻率f→∞，k0→∞，則入射角θi→90°，根據(jù)特征方程(7)，有

Jl(U)=0

。

(9)

式(9)是LPln模式遠(yuǎn)離截止的傳輸條件，l取不同值時(shí)對(duì)應(yīng)的根值標(biāo)注于貝塞爾函數(shù)曲線，如圖3所示。

V—?dú)w一化頻率；Jv—v階貝塞爾函數(shù)；J0—零階貝塞爾函數(shù)，黑線；J1—一階貝塞爾函數(shù)，紅線； J2—2階貝塞爾函數(shù)，藍(lán)線； J3—三階貝塞爾函數(shù)，紫線；LP—對(duì)應(yīng)各區(qū)間的線偏振模式場，相鄰2個(gè)根之間為一個(gè)區(qū)間；TE、TM、HE—矢量精確解，分別對(duì)應(yīng)LP01、LP11、LP02、LP12、LP03、LP13共6個(gè)線偏振模式場；LP01—HE11二重簡并，此為單模，也稱主模；LP11—HE21與TE01、TM01的簡并；LP02—HE12的二重簡并; LP12—HE22與TE02、TM02的簡并；LP03—HE13的二重簡并；LP13—HE23與TM03、TE03的簡并。圖3 貝賽爾曲線中模式的截止與遠(yuǎn)離截止

n1、n2—光纖纖芯、包層的折射率；TE、TM、HE、EH—精確場，不同的下標(biāo)代表不同的模式；LP01—HE11二重簡并，此為單模，也稱主模；LP11—HE21與TE01、TM01的簡并；LP21—EH11、HE31的簡并；LP02—HE12的二重簡并；LP31—EH21、HE41的簡并；LP12—HE22與TE02、TM02的簡并。圖4 線偏振模式(LP模式)的傳播系數(shù)與歸一化頻率V之間的關(guān)系

下面分析矢量精確解與標(biāo)量近似解的關(guān)系。

對(duì)于弱導(dǎo)光纖，橫向分量Hx、Ey、Ex、Hy的振動(dòng)模式為線偏振，矢量精確解HEl+1,m[Jl(V)]、EHl-1,m[Jl(V)](l=0, 1, …；m=1, 2, …)這2項(xiàng)本征值相同，旋向相反，2項(xiàng)進(jìn)行簡并后，合成線偏振光LPlm模式。

由表5和表6可以看出A、B、C三個(gè)影響因素中，影響因素順序?yàn)锳＞C＞B。由表5的均值結(jié)果分析，因素A的優(yōu)劣順序?yàn)?＞2＞1；因素B的優(yōu)劣順序?yàn)?＞2＞1；因素C的優(yōu)劣順序?yàn)?＞1＞2。最終選擇以A3B3C3為最佳水提工藝條件，即加10倍量水，煎煮3次，每次煎煮1.5 h。

當(dāng)l=0時(shí)，LP0m模式只有二重簡并，對(duì)應(yīng)精確解的混合波HElm。

當(dāng)l>0時(shí)，每個(gè)LPlm模式有四重簡并，徑向有2種模式，沿x方向或y方向偏振，由式(6)，徑向場分量包含角因子φ，角向有cosmφ、sinmφ之分[12]。當(dāng)l≥2時(shí)，LPlm模式由HEi+1,m與EHi-1,m模式簡并，而當(dāng)l=1時(shí)，LP1m模式由HE1+1,m與TE0m、TM0m模式簡并，簡并前、后的對(duì)應(yīng)關(guān)系如圖3、4所示。

3 纖芯截面LP模式的場強(qiáng)分布模擬及特性分析

根據(jù)LP模式的貝塞爾函數(shù)解，模擬幾個(gè)LP低階模的場強(qiáng)分布。設(shè)纖芯的內(nèi)徑為20 μm，光波波長為1.31 μm，纖芯和包層介質(zhì)折射率分別為n1=1.468，n2=1.458。

LP01模式的場強(qiáng)為

ELP01(x,y)=CJ0(U01r/a)/J0(U01),

(10)

式中：C為常數(shù)；U01為纖芯內(nèi)貝塞爾函數(shù)宗量。

由貝塞爾函數(shù)解的特性可知，J0(U01r/a)與J0(U01)在取值范圍(0，2.405)內(nèi)均為正值，因此導(dǎo)模的場沿徑向無零點(diǎn)。由于l=0，因此導(dǎo)模的場無角向分量。利用式(10)模擬LP01模式場在纖芯截面的場強(qiáng)分布，結(jié)果如圖5(a)所示。從圖中可以看出，光以纖芯軸心為中心分布，越靠近軸心時(shí)場強(qiáng)越大。

LP02模式的場強(qiáng)為

ELP02(x,y)=CJ0(U02r/a)/J0(U02), 3.823<

U02<5.520, 0

(11)

式中U02為纖芯內(nèi)貝塞爾函數(shù)宗量。

由貝塞爾函數(shù)解特性可知, 在宗量取值范圍內(nèi)，J0(U02r/a)有零點(diǎn)r=2.405a/U02，在該零點(diǎn)兩側(cè)值變號(hào)，而分母J0(U02)不變號(hào)。由于l=0，因此導(dǎo)模的場無角向分量，纖芯截面LP02模式的場強(qiáng)分布模擬結(jié)果如圖5(b)所示。由圖可知，纖芯中心為一個(gè)亮斑； 以纖芯軸心為中心分布一個(gè)亮圓環(huán)，以r=2.405a/U02環(huán)線為界，兩側(cè)場強(qiáng)變號(hào)。

LP11模式的場強(qiáng)為

ELP11(x,y)=C[J1(U11r/a)/J1(U11)]cosφ, 2.405<

U11<3.823, 0

(12)

(a)LP01模

(b)LP02模

(c)LP11模

J1(U11r/a)與J1(U11)在其宗量范圍內(nèi)均為正，角因子φ有零點(diǎn)π/2、3π/2，這2個(gè)角的邊線對(duì)稱且構(gòu)成連線，稱為零線，零線兩側(cè)場強(qiáng)符號(hào)相反，模擬結(jié)果如圖5(c)所示。

LP21模式的場強(qiáng)為

ELP21(x,y)=C[J2(U21r/a)/J2(U21)]cos 2φ, 3.823<

U21<5.136, 0

(13)

J2(U21r/a)與J2(U21)在其宗量范圍內(nèi)均為正，角因子φ有零點(diǎn)π/4、3π/4、5π/4和7π/4，且兩兩構(gòu)成對(duì)稱對(duì)角線，為2條零線，零線兩側(cè)場強(qiáng)符號(hào)相反，模擬結(jié)果如圖5(d)所示。

通過幾種低階模的模擬及其結(jié)果分析可知，與利用模式的截止、遠(yuǎn)離截止進(jìn)行的分析相比，兩者得出的結(jié)果一致，但是橫截面場強(qiáng)分布模擬更形象地顯示了各模式場強(qiáng)分布特點(diǎn)。

4 結(jié)語

本文中通過對(duì)光纖纖芯軸線橫截面的場強(qiáng)分布進(jìn)行模擬，并分析模擬結(jié)果，得到以下結(jié)論及啟示：

1)當(dāng)歸一化頻率V小于2.405時(shí)，光纖只有主模LP01，即為單模傳輸，模擬結(jié)果顯示，光線集中在光纖內(nèi)芯中心部分。單模光纖能夠完全避免模間色散，可以進(jìn)行大容量、長距離傳輸。

2)當(dāng)V大于2.405時(shí)，光纖內(nèi)的傳輸方式為多模傳輸，會(huì)受到色散影響，將無法實(shí)現(xiàn)大容量、遠(yuǎn)距離傳輸。

3)在設(shè)計(jì)制造光纖時(shí)，應(yīng)通盤考慮纖芯折射率、半徑、相對(duì)折射率差以及光頻率參數(shù)，使結(jié)構(gòu)參量滿足單模傳輸條件，保證光纖中只有LP01模式傳輸。

4)本文中主要研究階躍光纖傳播模式，在實(shí)際應(yīng)用中還有漸變光纖，階躍光纖是基礎(chǔ)。如果選擇漸變光纖的折射率按雙曲正割規(guī)律分布，能實(shí)現(xiàn)完善的子午線自聚焦效果，從而有效地消除模式色散的影響。