[論文レビュー] Ex-house 2D finite-element simulation of the whispering-gallery modes of arbitrarily shaped axisymmetric electromagnetic resonators

Mark Oxborrow|ArXiv.org|Jul 23, 2006

Photonic and Optical Devices参考文献 28被引用数 3

ひとこと要約

本論文は、商用PDEソルバを用いて、弱形式で記述されたマクスウェル方程式を用いる、追跡可能な有限要素法（FEM）に基づく手法を提示する。この手法により、軸対称な誘電体共鳴器における囁きギャラリー（WG）モードを高精度にシミュレートできる。本手法は、閉じた共鳴器および開放共鳴器の両方において、モード周波数、電場・磁場、フィリングファクター、および品質因数（Q）を高精度で計算可能であり、主な結果として、低温下におけるサファイアの誘電率の正確な推定値が得られた：ε⊥ = 9.285 (±0.010) および ε∥ = 11.366 (±0.010)。

ABSTRACT

It is described, explicitly, how a popular, commercially-available software package for solving partial-differential-equations (PDEs), as based on the finite-element method (FEM), can be configured to calculate the frequencies and fields of the whispering-gallery (WG) modes of axisymmetric dielectric resonators. The approach is traceable; it exploits the PDE-solver's ability to accept the definition of solutions to Maxwell's equations in so-called `weak form'. Associated expressions and methods for estimating a WG mode's volume, filling factor(s) and, in the case of closed(open) resonators, its wall (radiation) loss, are provided. As no transverse approxi-mation is imposed, the approach remains accurate even for so-called quasi-TM and -TE modes of low, finite azimuthal mode order. The approach's generality and utility are demonstrated by modeling several non-trivial structures: (i)two different optical microcavities [one toroidal made of silica, the other an AlGaAs microdisk]; (ii) a 3rd-order microwave Bragg cavity containing alumina layers (iii) two different cryogenic sapphire X-band microwave resonators. By fitting one of the latter to a set of measured resonance frequencies, the dielectric constants of sapphire at liquid-helium temperature have been estimated.

研究の動機と目的

• 商用有限要素PDEソルバを用いて、軸対称な誘電体共鳴器における囁きギャラリー（WG）モードを追跡可能かつ高精度にシミュレートする手法を開発すること。
• 実験的共鳴周波数からシミュレートされたモードを測定データにフィットさせることで、材料定数（例えば誘電率）の決定を可能にすること。
• 微小共鳴器、ブレグ・コアバ、低温サファイア共鳴器を含む多様な共鳴器幾何形状にわたる、本手法の汎用性と正確性を実証すること。
• 再現可能性を高め、研究コミュニティによる再利用を促進するために、自由に利用可能なソースコードおよび設定スクリプトを提供すること。
• 商用ソフトウェアを用いたシミュレーションにおけるトレーサビリティの課題に対処するため、弱形式の定式化および境界条件を明示的に文書化すること。

提案手法

• 本手法は、商用PDEソルバを介して、軸対称構造のマクスウェル方程式を弱形式で解く有限要素法（FEM）を採用する。
• 横方向の近似を回避することで、低次の準TMおよび準TEモードに対しても高精度を確保する。
• 境界条件および物性関係（例えば誘電率）をソルバのフレームワーク内で明示的に定義し、閉じた共鳴器および開放共鳴器のモデル化を可能にする。
• 放射損失または壁損失の推定を用いて、モード周波数、電場・磁場、フィリングファクター、および品質因数（Q）を計算する。
• シラス・トロイド、AlGaAsマイクロディスク、マイクロ波ブレグ共鳴器、低温サファイアXバンド共鳴器の多様な共鳴器タイプに本ソルバを適用する。
• 最小二乗フィッティング手順を用いて、誘電率を調整し、シミュレートされたモードと実験的に測定された共鳴周波数を一致させる。

実験結果

リサーチクエスチョン

• RQ1 商用FEMベースのPDEソルバを、系統的かつ追跡可能な方法で、軸対称な誘電体共鳴器における囁きギャラリーモードを高精度にシミュレートできるか？
• RQ2 本手法は、マイクロディスク、トロイド、ブレグ共鳴器といった非自明な共鳴器幾何形状に対し、どの程度正確にモード周波数および電場分布を予測できるか？
• RQ3 本手法を用いて、実験的共鳴データから材料定数（特に誘電率）をどの程度正確に抽出できるか？
• RQ4 モード同定の不確実性が、フィッティングされた誘電率の精度に与える影響は何か？ また、その影響をどのように定量化できるか？
• RQ5 メッシュ密度や幾何的不完全性といった系統的誤差が、実験測定値と比較してシミュレーションの忠実度に与える影響はどのようなものか？

主な発見

• 本手法は、シリカ・トロイド、AlGaAsマイクロディスク、3次マイクロ波ブレグ共鳴器、および2つの低温サファイアXバンド共鳴器を含む、多様な複雑な共鳴器において、囁きギャラリーモードを成功裏にシミュレートした。
• 低温サファイア共鳴器において、最良のフィットによる誘電率は ε⊥ = 9.285 (±0.010) および ε∥ = 11.366 (±0.010) であり、16個の実験的測定共鳴周波数にシミュレートされたモードをフィットさせることで得られた。
• フィット値は、メッシュ密度や幾何形状の不確実性に対して最小限の感度を示しており、モード同定の誤差が主な不確実性要因であることが示された。
• シミュレートされた中心周波数と測定された周波数との残差は、数千分の数パーセントのオーダーであり、結晶のずれがその差異の主な要因である可能性が高い。
• 独立した別の研究グループによるシミュレーション結果と照合されたことで、本手法の正確性が検証された。
• 本研究で使用したすべてのシミュレーションのソースコードおよび設定スクリプトは、公開されており、再現性を高め、今後のモデル拡張および検証を容易にする。