[論文レビュー] Broadband transmission properties of multilayered structures

Victor Grigoriev, Fabio Biancalana|arXiv (Cornell University)|Jun 30, 2011

Photonic Crystals and Applications参考文献 9被引用数 4

ひとこと要約

本稿では、散乱行列の極をパラメータとして用い、四分波条件下における多層構造の広帯域透過スペクトルの正確な解析的公式を提示する。レイの伝播モデルと調和反転を組み合わせることで、透過、群速度遅延、パルス再形状ダイナミクスの効率的計算が可能となり、数値的伝達行列法をはるかに超える物理的洞察を提供する。

ABSTRACT

The formalism of the scattering matrix is applied to describe the transmission properties of multilayered structures with deep variations of the refractive index and arbitrary arrangements of the layers. We show that there is an exact analytical formula for the transmission spectrum, which is valid for the full spectral range and which contains only a limited number of parameters for structures satisfying the quarter-wave condition. These parameters are related to the poles of the scattering matrix, and we present an efficient algorithm to find them, which is based on considering the ray propagation inside the structure and subsequent application of the harmonic inversion technique. These results are significant to analyze the reshaping of ultrashort pulses in multilayered structures.

研究の動機と目的

• 任意の層配置を持つ多層構造における広帯域透過の正確な解析的記述を開発すること。
• 透過スペクトルを完全に決定する物理的パラメータ（共鳴極と強度）を同定すること。
• 標準的な伝達行列法の代替として、計算的に効率的かつ物理的に直感的な手法を提供すること。
• 超短パルスに対する信号伝播の時間領域解析（群速度遅延およびパルス再形状を含む）を可能にすること。
• フォトニクスクリスタル、分散反射板、非線形光学への応用に適用可能なフレームワークを確立すること。

提案手法

• 透過スペクトルは、複素周波数平面における極の和としてモデル化される：$ T(\omega) = \sum_{p=-\infty}^{\infty} \frac{\sigma_p}{\omega - \omega_p} $、ここで$ \sigma_p $および$ \omega_p $は共鳴強度と周波数である。
• 四分波構造では、極の周期性（$ \omega_p + 4k\omega_{\rm{qw}} $）により、独立な極の数が有限集合$ P $に減少する。
• 時間領域応答関数$ G(t) $は逆フーリエ変換により導出され、指数的減衰モードの離散和として得られる：$ G(t) = \sum_{n=0}^{\infty} a_n \delta(t - n\tau) $。
• 係数$ a_n $は透過スペクトルのフーリエ成分から計算され、信号の畳み込みが効率的に行える：$ f_{\rm{out}}(t) = \sum_{n=0}^{\infty} a_n f_{\rm{in}}(t - n\tau) $。
• レイ追跡と調和反転に基づく効率的アルゴリズムが開発された：出力信号は$ \tau = n_m d_m / c $の間隔でサンプリングされ、極は行列固有値問題$ \hat{H}_1 \mathbf{V}_p = \exp(-i\omega_p \tau) \hat{H}_0 \mathbf{V}_p $から抽出される。
• 共鳴強度$ \sigma_p $は固有ベクトルの正規化により計算される：$ \sigma_p = \frac{i}{\tau} \frac{[\mathbf{V}_p^T \cdot \mathbf{A}]^2}{\mathbf{V}_p^T \cdot \hat{H}_0 \mathbf{V}_p} $。

実験結果

リサーチクエスチョン

• RQ1 任意の層配置を持つ多層構造の全周波数帯域透過スペクトルを正確な解析的公式で記述できるか？
• RQ2 広帯域透過特性を記述するために、散乱行列の極をどのように効率的に計算できるか？
• RQ3 このような構造において、時間領域応答と複素周波数の極との間にはどのような物理的関係があるか？
• RQ4 四分波条件が極の分布に周期性をもたらし、スペクトル記述をどのように簡略化するか？
• RQ5 この形式的枠組みは、群速度遅延や時間的伸長を含む超短パルス再形状をどの程度正確に予測・分析できるか？

主な発見

• 透過スペクトルは有限個の極の和として正確に記述される：$ T(\omega) = \sum_{p=1}^{P} \frac{i\sigma_p \tau}{1 - \exp[i(\omega - \omega_p)\tau]} $、これはすべての周波数で有効である。
• 四分波構造では、$ (-2\omega_{\rm{qw}}, 2\omega_{\rm{qw}}] $の範囲内にのみ$ P $個の独立した極が存在し、他のすべての極は周期性により得られる。
• 時間領域応答$ G(t) $は、$ \tau = n_m d_m / c $の間隔でサンプリングされた指数的減衰モードの離散和であり、係数$ a_n $は$ n $とともに指数的に減衰する。
• 共鳴周囲での群速度遅延$ \tau_D $は、$ \tau_D \approx -2 / (\pi \cdot \text{Im}[\omega_p / \omega_{\rm{qw}}]) $として近似可能であり、共鳴の減衰時間と一致する。
• 調和反転アルゴリズムは、$ 2P $個の信号サンプルのみで$ P $個の極とその強度を抽出でき、高い効率性を実現する。
• この手法により、パルス再形状（時間的伸長および多重共鳴系における干渉効果を含む）の正確な予測が可能になる。