ryujimiyaの日記

C#を使って数値解析したい

周波数領域FEMによるフォトニック結晶導波路CROW(coupled-resonator optical waveguide)の固有値問題、伝達問題の計算

1. はじめに

フォトニック結晶に空洞(欠陥、共振器)を1ロッドおきに並べて結合させると導波路として使えます。この構造をCROW(coupled-resonator optical waveguide)と呼んでいます。

CROWの固有モード、および反射・透過特性を去年の11月終わりころIvyFEMで実装し計算しましたが記事にしていなかったのでまとめました。

60°三角形格子、エアホール型フォトニック結晶のTEzモードを対象にしています。

エアホールの半径r = 0.30a、a:格子定数で、基盤の屈折率n1 = 3.40、エアホールの屈折率n2=1として計算しました。

 

f:id:ryujimiya:20200205141514j:plain

 

2. CROWの固有モード計算

CROWの1周期を次のように取り、左右の境界にFloquet(Bloch)の周期境界条件を適用します。

f:id:ryujimiya:20200205141839j:plain

2.1. 偶モード

f:id:ryujimiya:20200205142116j:plain

f:id:ryujimiya:20200205142127j:plain

f:id:ryujimiya:20200205142139j:plain

f:id:ryujimiya:20200205142150j:plain

2.2. 奇モード

f:id:ryujimiya:20200205142316j:plain

f:id:ryujimiya:20200205142328j:plain

f:id:ryujimiya:20200205142341j:plain

f:id:ryujimiya:20200205142406j:plain

 

3. CROW直線導波路の反射・透過特性

f:id:ryujimiya:20200205142507j:plain

f:id:ryujimiya:20200205142524j:plain

f:id:ryujimiya:20200205142534j:plain

f:id:ryujimiya:20200205142545j:plain

f:id:ryujimiya:20200205142600j:plain

偶モードを入射させたときの計算です。

ほぼ|S11|は0になっており、100%伝送できており、CROWが導波路として機能することが分かります。

4. CROW-CROW60°ベンドの反射・透過特性

f:id:ryujimiya:20200205142845j:plain

f:id:ryujimiya:20200205142858j:plain

f:id:ryujimiya:20200205142908j:plain

f:id:ryujimiya:20200205142917j:plain

f:id:ryujimiya:20200205142926j:plain

f:id:ryujimiya:20200205142935j:plain

f:id:ryujimiya:20200205142947j:plain

f:id:ryujimiya:20200205142957j:plain

異なる方向のCROWの接続としては60°ベンドが考えられます。

偶モードを入射させたときの計算です。

奇モードが存在するa/λ = 0.24~0.25あたりでは偶モードが奇モードに変換されています。奇モードが存在しない時は偶モードが出力されています。

5.CROW-欠陥導波路 90°ベンドの反射・透過特性

f:id:ryujimiya:20200205143517j:plain

f:id:ryujimiya:20200205143529j:plain

f:id:ryujimiya:20200205143541j:plain

f:id:ryujimiya:20200205143550j:plain

f:id:ryujimiya:20200205143601j:plain

f:id:ryujimiya:20200205143611j:plain

f:id:ryujimiya:20200205143621j:plain

f:id:ryujimiya:20200205143629j:plain

f:id:ryujimiya:20200205143639j:plain

f:id:ryujimiya:20200205143647j:plain

CROWと欠陥導波路の接続は90°ベンドで実現できます。

偶モードを入射させたときの計算です。

 

6. まとめ

フォトニック結晶導波路CROWの固有値問題および伝達問題を計算しました。