• 沒有找到結果。

第三章 人造介質材料分析方法與設計

3.4 色散關係之討論

3.4.1 等效參數萃取法的結果

首先,利用等效參數萃取法分析,其模擬環境的設定如3.1 節所述,在X-Y 平面和X-Z平面設定週期性邊界條件,在Y-Z平面設定port,由於此漁網結構之單 晶格有介質板延伸至port上,所以必須使用de-embedding的方法,其設置如圖 3-5。對於波垂直入射至單晶格,必須設定週期性邊界條件的相位差θy和θz皆為零 度,並設定port激發的波為極化方向z方向的基礎模態。

z x y

z x y

z x y

z x y

圖 3-5 漁網結構的 de-embedding 設置

3-6(a)、(b)分別是模擬所得S參數的大小和相位,等效參數萃取法將利用 此S參數,和視為等效介質的厚度dx來求得等效參數。

5 6 7 8 9 10 11 12 0

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Frequency(GHz) Magnitude of S-parameter

S11 S21

圖 3-6 (a) S 參數大小

5 6 7 8 9 10 11 12

-150 -100 -50 0 50 100 150

Frequency(GHz) Phase of S-parameter

S11 S21

3-6 (b) S 參數相位

由3.1 節的方法求得等效參數,圖 3-7 是等效折射係數,可以看出在頻率 7GHz 之前,折射係數為純虛數,並不存在電磁波傳播的模態,7GHz為此結構的截止 頻率(cutoff frequency),此截止頻率來自於此單晶格的結構。在截止頻率之後至 10.15GHz,neff為正實數,表示此頻率範圍為傳播的模態。

10.15GHz之後,neff為負實數加上虛數,此區代表帶隙出現了。在3.1 節 亦提到,在Bragg regime附近,基礎模態會受到高階模態的影響,造成模態偶合;

Notomi亦提到光子晶體的異常折射現象是來自於帶隙附近的特殊色散關係 [11],所以在帶隙附近用等效折射係數來描述材料特性是不完整的,也因此在帶 隙出現後,也就是10.15GHz之後,就不在等效係數萃取的討論範圍。

所以等效係數萃取法在7GHz~10.15GHz 是我們討論的範圍,此頻段內等效 折射係數都小於一,而不同於傳統材料的大於一或等於一,因此可以有新的運用。

5 6 7 8 9 10 11 12

-1.5 -1 -0.5 0 0.5 1 1.5

Frequency(GHz)

Ef fe ct iv e r e fr act iv e i n dex

Real(n) Imag(n)

圖 3-7 等效折射係數

在等效係數萃取的討論範圍內,非截止頻率的阻抗都為實部,因此人造介質 材料的特性就可以用等效波阻抗來描述,而視為傳輸線的特性阻抗(characteristic impedance)。此等效波阻抗對之後空間波束分歧器設計的阻抗匹配(impedance matching)是相當重要的。

5 6 7 8 9 10 11 12

-25 -20 -15 -10 -5 0 5 10 15 20 25

Frequency(GHz)

N o rm al iz ed ef fe ct iv e i m pedance

Real(z) Imag(z)

圖 3-8 正規化的等效波阻抗

圖3-9 是等效介電係數,在截止頻率介電係數為零,可以稱為此人造介質材 料的電漿頻率,在討論的頻率範圍內,介電係數皆小於一。

圖3-10 是等效導磁係數,在討論的頻率範圍內,導磁係數皆大於一。

5 6 7 8 9 10 11 12 -2

-1.5 -1 -0.5 0 0.5 1 1.5 2

Frequency(GHz)

E ffe c tiv e p e rm ittiv ity

Real(ε) Imag(ε)

圖 3-9 等效介電係數

5 6 7 8 9 10 11 12

-25 -20 -15 -10 -5 0 5 10 15 20 25

E ff e ct iv e per m eabi lit y

Real(μ) Imag(μ)

由圖3-11 的色散關係可以知道電磁波在空氣中的傳播常數是隨頻率線性增 加,但在人造介質材料中,是由截止頻率才開始傳播,其傳播常數在截止頻率為 零,由圖上可看出人造介質材料的傳播常數在7GHz 到 10GHz 都低於在空氣的 傳播常數。當電磁波從人造介質材料入射至空氣時,因為遵循切線波數連續,折 射後的電磁波將往法線偏移,之後的高指向性空間波束分歧器設計就是利用這個 特性。

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

βdx/2π Νοrmalized phase constant

5 6 7 8 9 10 11 12 13 14 15 16 17

Frequency (GHz)

Γ−X Γ−M Light line Band structure in Brillouin zone

圖 3-11 電磁波在人造介質材料和空氣的色散關係

.4.2 特徵模態分析法的結果與討論

,使用特徵模態分析法來分析漁網結構的單晶格。圖3-4(a)是漁網結

構的單晶格,邊界條件的設定是在 和 平面設定週期性邊界條件,相位差

分別為 xy, 平面設為 ;模擬軟體計算的模態設為 個,計算特徵頻 3

接下來

X-Z Y-Z

θ θ X-Y PEC 5

率起始點為1GHz。

其計算得到的頻帶結構用布里淵區來表示,如圖3-12。由於單晶格結構的複 雜,可以看出特徵頻率點相當的混亂,整個布里淵區在7.5GHz 附近都有兩個模 態,而影響到人造介質材料想要的基礎模態之判定。所以,為了精準的得到基礎 模態的色散關係,必須對不同週期性邊界條件解得的電場結果作判斷,將正確的 電場對應的特徵頻率連成一線即為圖3-12 的基礎模態。

7 8 9 10 11 12 13 14 15 16 17

Fr eq ue nc y (G H z)

Γ X M Γ

圖 3-12 布里淵區的頻帶結構

觀察由特徵模態求得的電場和Poynting’s vector,如圖 3-13 和圖 3-14。圖 3-13 是週期性邊界條件的相位差為0 度的電場分佈,可以看出在基礎模態下,電場的 方向為z 方向,和等效介質參數萃取法一致。圖 3-14 是週期性邊界條件的相位 差為10 度的 Poynting’s vector,可以看出能量傳播的方向和相位差的方向相同。

圖 3-13 單晶格的電場分佈

Phase difference=10°

Phase difference=10°

圖 3-14 相位差為 10°的 Poynting’s vector

將頻帶關係用三維的圖來表示,如圖3-15。x軸和y軸分別是單晶格x方向和y 方向的相位差,z軸是特徵頻率,對此圖每隔 0.5GHz做頻率的等高線,並將x軸 和y軸的相位差利用kxx/dx和kyy/dy的關係轉成kx和ky,得到等頻率面(equal frequency surface)的相位關係,如圖 3-16。

圖 3-15 三維頻帶結構

3-16,在頻率為 9.7GHz 之下,等頻率面為圓形,表示此漁網結構在 x 方 向和y 方向具備等向性。在 X-Y 平面上,不管波向量的方向為何,其傳播常數 大小皆為等頻率圓的半徑,這正是等效介質參數萃取法中只設定垂直入射,而不 需要設定斜向入射的原因。因此,等效介質參數萃取法和特徵模態分析法的結果 是否一致就是十分重要的問題。

在圖3-17 繪製出等效介質參數萃取法和特徵模態分析法的結果,比較等效 介質參數萃取法的等效折射係數,以及特徵模態的正規化相位常數作比較,對於 實部而言,可以看出在截止頻率至10.15GHz 是一致的,但是當參數萃取中帶隙 出現後,實部大小就不一致了,這可以對應至第3.4.1 節提到的在帶隙附近不能 夠用等效係數來描述等效介質的特性。以上兩種方法對於人造介質材料在基礎模 態的分析都有相當一致的結果。

6 8 10 12

5 7 9 11

Frequency(GHz) -2

-1 0 1 2

-1.5 -0.5 0.5 1.5

n

Real(n) from extraction Imag(n) from extraction Real(n) from Eigenmode

圖 3-17 等效介質參數萃取法和特徵模態分析法的結果比較

3.5 漁網結構等效電路模型

對於人造介質材料的特性,除了用等效參數來描述外,也希望能夠找出等效 電路,以加快分析速度。

在Marcuvitz 的 Waveguide handbook[12]裡,5.18 和 5.19 節給出了沒有厚度 而且上下無窮延伸的平面金屬柵(metallic strip)之等效電路,金屬柵方向分別有和 電場與磁場平行的兩種平面波入射形態,如圖3-18 和圖 3-19,等效電路數值如 Appendix A。

E

圖 3-18 (a) 金屬柵之正面與電場方向 (b) 等效電路

H

圖 3-19 (a) 金屬柵之正面與磁場方向 (b) 等效電路

本論文的漁網結構可以看成圖3-18 之金屬柵和圖 3-19 之金屬柵轉 90 度重 疊而成,水平和垂直的金屬柵結合成金屬網格,如圖3-20。此金屬網格的等效電 路由Ulrich[13]將Marcuvitz的兩種金屬柵延伸,提出了當波長小於金屬網格的週 期時,其穿透和反射的特性與金屬柵相近。以漁網結構而言,稱為電感性網格,

其行為與圖3-18 的電感性金屬柵類似,但在入射波的波長接近金屬柵週期時(通

網格的等效電路之電感再並聯電容,電感和電容如同Marcuvitz所提供(圖 3-18 和 圖3-19),並再加入在λ=dx時完全穿透的修正項,此電容修正成為

B ⎟⎟⎠ ⋅B

⎜⎜ ⎞

⎛ −

=

−1 0

0

' 0

ω ω ω

ω ω where = , ω0 =1

ω dλy (3.4.1) 電感為原始的數值

X'=X (3.4.2)

dy dx

az ax

/2

dy dx

az ax

/2

圖 3-20 金屬網格的單晶格

而金屬網格週期性在y方向的堆疊,可用傳輸線和金屬網格的等效電路作堆 疊,如圖3-21。其傳輸線的長度為dy的一半,傳輸線的特性阻抗為平行板波導的 波阻抗。

圖 3-21 金屬網格單晶格的等效電路模型

金屬網格在z方向和x方向週期性延伸的特性已包含至LC等效電路裡,現在 只要考慮在y方向的並聯LC和傳輸線的堆疊,把此單晶格表示成ABCD矩陣,以

兩端接特性阻抗為Z0的傳輸線,轉換為S參數,如圖 3-22。圖 3-22 是單晶格的S 參數,電磁波在18GHz附近完全穿透,此漁網結構具有濾波的特性;當串接三 個單晶格時,如圖3-23。在 7.4GHz和 9.3GHz完全穿透,15~20GHz是通帶(pass band),並且有兩個峰值(peak),這是來自於三個單晶格串接時,相鄰的漁網結構 形成兩個共振腔,而有兩個相鄰的峰值。

5 10 15 20

Frequency (GHz) 0

0.2 0.4 0.6 0.8 1

S-parameter

S11 simulation S11 model S21 simulation S21 model

圖 3-22 單晶格的散射參數

5 10 15 20 Frequency (GHz)

0

0 0.1 0.2 0.3 0.4 0.5 kxdx/2π

0 4 8 12 16 20

Freq uenc y (GHz)

Eigenmode solver Model

圖 3-24 正向入射的色散關係

等效電路的存在更是為了加快分析與設計,當改變金屬網格的結構參數時,

等效電路所得的色散關係也要因此而改變。因為漁網結構是屬於電感性金屬網格 中,在相同晶格常數下,對截止頻率影響最大的參數是ax的大小,越大的ax對應 至等效電路越大的並聯電感值,而使截止頻率降低。圖3-25 是不同ax由特徵模 態和等效電路所求得的截止頻率,截止頻率的差距在4%以內,因此,由等效電 路模型得到的截止頻率具有相當的準確性。

7 8 9 10 11 12 13

Cutoff Frequency (GHz)

dx=13.8mm dz=12.8mm az=9mm Model

Eigenmode solver

圖 3-25 金屬網格在不同結構參數下的截止頻率

由上式可以知道當考慮介質後,電容值加大,而導致諧振頻率下降,使得單 晶格的截止頻率下降。在本漁網結構中,介質基板相對於晶格常數而言是相當的 小,因此電容值的估算和(3.4.4)式的情形有所不同,對截止頻率的降低的程度較 小。

dy

dx h

dy

dx h

圖 3-26 兩層金屬網格堆疊的單晶格

圖 3-27 兩層金屬網格的電路模型

以上電路模型的討論都是多層堆疊的漁網結構,並未討論本論文提出的漁網 結構。電磁波垂直入射,X-Z 平面的等效電路是由水平和垂直的金屬柵構成,但 對Y-Z 平面而言, x 方向磁場正好對金屬網格上的環形成電感,因此對於漁網

3-28。此電感值和矩形環內面積的大小相

圖 3-28 漁網結構的電路模型

3-28 是漁網結構在基礎模態的電路模型,但在本論文裡還尚未找到Xy和k,

因此漁網結構的設計,還必須藉助圖3-21 的等效電路,精確的算出截止頻率和 色散關係。由以上的討論,知道當漁網結構的介質基板,相對於晶格常數很小時,

其上的雙面漁網結構的偶合以及基板的介電係數大小,對截止頻率的影響很小,

而可以用圖3-21 的等效電路來估計截止頻率。又圖 3-21 的等效電路,因為未考 慮y方向的電感,算出的截止頻率較漁網結構為低。

儘管目前尚未找出對應至漁網結構的等效電路參數,但仍可以用已知的金屬 網格堆疊的等效電路模型,估計漁網結構基礎模態的大致範圍和截止頻率,這對 單晶格尺寸上的設計有相當大的幫助。

值得一提的是,以上等效電路的討論都是基於多層堆疊週期性結構的觀點,

值得一提的是,以上等效電路的討論都是基於多層堆疊週期性結構的觀點,

相關文件