潘繼環(huán)，張?jiān)?/p>

(河池學(xué)院 物理與機(jī)電工程學(xué)院，廣西 宜州 546300)

0 引言

矩形波導(dǎo)是條形波導(dǎo)中最常見的也是最基本的波導(dǎo)結(jié)構(gòu)，通常是由金屬材料(銅、鋁等)制成矩形截面、內(nèi)部填充空氣介質(zhì)的金屬波導(dǎo)[1－3]。而矩形波導(dǎo)在傳導(dǎo)電磁波的過程中具有損耗低、電磁屏蔽性能好等優(yōu)點(diǎn)，在微波通信設(shè)備領(lǐng)域有著非常廣泛而巨大的應(yīng)用價(jià)值，如在通信傳輸線平面集成化、小型化等方面具有重要的應(yīng)用[1－12]，近年來有關(guān)波導(dǎo)的研究一直是各國物理學(xué)者研究的熱點(diǎn)之一。

通過查詢相關(guān)文獻(xiàn)資料發(fā)現(xiàn)，對(duì)于一維結(jié)構(gòu)的條形平板波導(dǎo)已經(jīng)有很多的研究報(bào)道[4－12]，而有關(guān)矩形波導(dǎo)的研究還比較少見，雖然矩形波導(dǎo)也屬于一維平板結(jié)構(gòu)，但研究其芯層導(dǎo)波特性時(shí)需要考慮兩個(gè)方向在其橫截面上對(duì)光場(chǎng)的限制作用，這就涉及到一維結(jié)構(gòu)的二維問題，這類問題的處理比條形平板波導(dǎo)更為復(fù)雜。

基于此，本文主要采用馬卡梯里假設(shè)[7]，對(duì)波導(dǎo)進(jìn)行近似處理的方法，并通過Matlab編程模擬繪制芯層寬度對(duì)有效折射率的影響，找出波導(dǎo)孔徑對(duì)矩形波導(dǎo)模式特性影響的規(guī)律，使矩形波導(dǎo)具有更廣泛的應(yīng)用前景。

1 模型及計(jì)算方法

1.1 矩形波導(dǎo)結(jié)構(gòu)

矩形波導(dǎo)結(jié)構(gòu)模型[4－10]的橫截面，如圖1所示，將截面分成九個(gè)區(qū)，中間為波導(dǎo)芯層，緊鄰芯層的四個(gè)白色區(qū)域?yàn)椴▽?dǎo)包層，四個(gè)陰影區(qū)域?yàn)榻菂^(qū)。其中波導(dǎo)芯層的x方向和y方向的寬度分別為a、b，選取芯層的中心為原點(diǎn)建立直角坐標(biāo)系，z軸垂直于紙面向里，并定義電磁波沿著z軸方向傳播。其中參數(shù)λ0=1.55 μm，折射率 n1=1.5，n2=n3=n4=n5=1.2。

研究方法使用馬卡梯里法。因?yàn)轳R卡梯里法考慮到如果離截止點(diǎn)比較遠(yuǎn)，則光能量高度集中在芯區(qū)，透入到2、3、4、5四個(gè)區(qū)的光能就很少，而四個(gè)陰影的角區(qū)中光能忽略不計(jì)，這樣可以使問題大大簡化。利用這一條合理而又巧妙的假設(shè)，很容易利用已解三層平板波導(dǎo)的結(jié)果來求解傳播常數(shù)和模場(chǎng)函數(shù)。

圖1 矩形波導(dǎo)橫截面圖

1.2 特征方程的計(jì)算

從麥克斯韋電磁波方程組可以推出，矩形波導(dǎo)Eymn導(dǎo)模的磁場(chǎng)分量Hx0(x，y)的波動(dòng)方程為

則φx(x)、φy(y)滿足橫向亥姆霍茲方程，即

方程(2)，(3)中的kix，kiy分別表示5個(gè)區(qū)域的傳播系數(shù)。而方程(2)，(3)的通解分別為

式中 A、A3、A5、B、B2、B4、ζ、η 為積分常數(shù)，為各自區(qū)域的振幅因子，ζ與 η 是相位因子。

在Eymn導(dǎo)模x軸方向中，在處，利用 φx(x)連續(xù)和連續(xù)的邊界條件以及在處，利用φx(x)連續(xù)和連續(xù)的邊界條件，即可求解特征方程。

方程(6)、(7)中m、n為x軸和y軸方向的模階數(shù)。將式分別代入方程(6)和方程(7)，得到有效折射率與芯層寬度的關(guān)系為

2 計(jì)算結(jié)果和分析

2.1 芯層寬度a對(duì)有效折射率Nx的影響

圖2 導(dǎo)模x軸向的傳輸曲線

矩形波導(dǎo)各層介質(zhì)的參數(shù)見1.1小節(jié)，取模階數(shù)為m=0～5時(shí)，對(duì)方程(8)利用Matlab編程模擬，繪制出各模式有效折射率Nx與芯層寬度a的關(guān)系，如圖2所示。

從圖2中不難看出，隨著芯層寬度的增大各模式有效折射率均出現(xiàn)衰減現(xiàn)象，寬度a＞0.4 μm時(shí)，波導(dǎo)支持m1≥1的多模傳輸，且各模式折射率隨著芯層寬度的增大而下降，而有效折射率Nx曲線向大孔徑方向移動(dòng)。即當(dāng)a=0～8 μm增大時(shí)，零階模對(duì)應(yīng)的有效折射率由0.9衰減到0.0907左右，一階模式的有效折射率從2衰減到0.1813，二階模從2衰減至0.2718，三階模從 2衰減至 0.3620，四階模從 2衰減至0.4520，五階模從2衰減至0.5415。這些衰減現(xiàn)象主要是隨著波導(dǎo)孔徑的增大，模式出現(xiàn)襯底輻射，從而使其能量損耗所致。值得注意的是，芯層寬度a在0～8 μm范圍內(nèi)，零階模的有效折射率均小于0.9，因此該結(jié)構(gòu)波導(dǎo)在x軸向上不支持零階模的傳輸，所以零階模的波導(dǎo)中主要以襯底輻射的方式逐漸損耗能量。

電磁波能在芯層中傳導(dǎo)的條件是在有效折射率1.2＜Nx＜1.5之間，因此芯層寬度無論小到什么程度，零階模(基模)均不滿足傳導(dǎo)的條件，因此該矩形波導(dǎo)不支持基模的傳輸。高階模式中，隨著模階數(shù)的增大，導(dǎo)模的區(qū)域也會(huì)增大。

從圖2中還可以看到，每一條傳輸曲線中都存在一個(gè)凹點(diǎn)，這些凹點(diǎn)對(duì)應(yīng)的有效折射率值相同，該處為導(dǎo)模的截止點(diǎn)，也稱作導(dǎo)模和空間輻射模的交界點(diǎn)。

2.2 芯層寬度b對(duì)有效折射率Ny的影響

同樣，矩形波導(dǎo)各層介質(zhì)的參數(shù)見1.1小節(jié)，當(dāng)模階數(shù)為n=0～5時(shí)，對(duì)方程(9)利用Matlab編程模擬，繪制得到芯層寬度b和有效折射率Ny的關(guān)系，如圖3所示。

從圖3中可知，芯層寬度b＞0.4 μm時(shí)，波導(dǎo)也支持n1≥1的多模傳輸，且各模式折射率隨著芯層寬度的增大而下降，有效折射率Ny曲線同樣也向大孔徑方向移動(dòng)。即當(dāng)b=0～8 μm增大時(shí)，零階模對(duì)應(yīng)的有效折射率由2衰減到1.2處后中斷、再從0.9衰減到0.0928左右，一階模式的有效折射率從2衰減到0.1813，二階模從2衰減至0.2780，三階模從2衰減至0.3702，四階模從 2衰減至 0.4619，五階模從 2衰減至0.5528。

在有效折射率1.2＜Ny＜1.5之間的區(qū)域，波導(dǎo)同樣不支持基模的傳輸。而在同一模式階數(shù)中，y軸向的有效折射率Ny呈衰減的趨勢(shì)，并且隨著模階數(shù)的增大，導(dǎo)模的區(qū)域也會(huì)隨著增大。

圖3 導(dǎo)模y軸向的傳輸曲線

在圖3中同樣也存在凹點(diǎn)，即導(dǎo)模和空間輻射模的交界點(diǎn)。并且還發(fā)現(xiàn)，每條傳輸曲線的導(dǎo)模區(qū)域存在一個(gè)折彎處，該處的有效折射率接近包層折射率，這說明該點(diǎn)很可能為導(dǎo)模的截止點(diǎn)。而有效折射率大于對(duì)應(yīng)的截止點(diǎn)的區(qū)域?yàn)樵撃Ｊ降慕刂箙^(qū)域，該模式的截止區(qū)域會(huì)由低一級(jí)階模式的截止區(qū)域進(jìn)行補(bǔ)償，形成一條完整的模式傳輸曲線。如五階模由四階模的截止區(qū)域“補(bǔ)償”組成新的傳輸曲線;四階模由三階模的截止區(qū)域“補(bǔ)償”組成新的傳輸曲線，以此類推。而零階模由于沒有更低一階模式的截止區(qū)域“補(bǔ)償”，即零階模只存在空間輻射模，且零階模不可導(dǎo)。

3 結(jié)論

上述是采用馬卡梯里近似法，并利用電磁波動(dòng)理論，求解矩形波導(dǎo)模式特征方程，研究了芯層寬度對(duì)矩形波導(dǎo)有效折射率的影響規(guī)律，得出以下結(jié)論:當(dāng)矩形波芯層寬度足夠大時(shí)，矩形波導(dǎo)支持多個(gè)導(dǎo)波模式的傳輸，且隨著模階數(shù)的增大，兩個(gè)獨(dú)立分量有效折射率Nx和Ny的曲線均向大孔徑方向移動(dòng);且兩個(gè)分量所有模式的有效折射率均隨芯層寬度a或b的增大而下降，這些影響規(guī)律對(duì)矩形波導(dǎo)的實(shí)際應(yīng)用具有一定的理論指導(dǎo)意義。

[1]佘守憲.導(dǎo)波光學(xué)物理基礎(chǔ)[M].北京:北方交通大學(xué)出版社，2002:124－152.

[2]詹景坤，呂天慧，陳燦輝，等.一種兩毫米波矩形波導(dǎo)帶通濾波器的研究與設(shè)計(jì)[J].電子測(cè)試，2015(6):18－22.

[3]馬春生，劉式墉.光波導(dǎo)模式理論[M].長春:吉林大學(xué)出版社，2006:172－203.

[4]孫杰，楊劍鋒，閆肅，等.等離子體輔助平板波導(dǎo)的傳輸特性及應(yīng)用研究[J].物理學(xué)報(bào)，2015，64(7):078402.

[5]陳憲鋒，沈小明，蔣美萍，等.對(duì)稱單負(fù)介質(zhì)包層平面波導(dǎo)的模式特征[J].物理學(xué)報(bào)，2008，57(6):3578－3582.

[6]張靜，何英，李春芳，等.包含左手材料的四層平板波導(dǎo)中的光導(dǎo)模[J].光學(xué)學(xué)報(bào)，2009，29(10):2673－2680.

[7]蒙成舉，盧強(qiáng)華，歐梅蓮，等.單負(fù)材料填充包層的矩形波導(dǎo)模式特性[J].廣西科學(xué)，2015，22(1):94－98.

[8]許江勇，段國華，張衛(wèi)龍，等.Metal－LHM－Metal三層對(duì)稱平板波導(dǎo)的傳輸特性研究[J].廣西科學(xué)，2012，19(3):236－238.

[9]吳美燕，羅硯濃，歐梅蓮，等.含左手材料四層平板波導(dǎo)的模式特性[J].廣西科學(xué)，2015，22(1):99－103.

[10]唐裕霞，羅成，彭家華.分別含左右手材料的三層平板波導(dǎo)中光導(dǎo)模的比較研究[J].重慶工商大學(xué)學(xué)報(bào)(自然科學(xué)版)，2014，31(7):63 －69，76.

[11]諸一琦，陳憲鋒，楊燁.特異材料對(duì)矩形波導(dǎo)導(dǎo)波模式的影響[J].常州工學(xué)院學(xué)報(bào)，2013，26(5):20－24.

[12]雷前召.矩形波導(dǎo)中電磁場(chǎng)分布特征研究[J].電子設(shè)計(jì)工程，2011，19(19):162－164.