転送行列法とは何か？

1. はじめに

転送行列法（TMM）は、多層薄膜光学構造内における電磁波の伝播をモデル化するための高速かつ計算効率に優れた手法である。屈折率の異なる多層材料から構成される光学構造において、波長依存性のある反射率、透過率、光吸収、光干渉、ならびに電界分布を算出する目的で、光学、フォトニクス、半導体デバイスシミュレーションの分野で広く利用されている。その応用例としては、薄膜太陽電池、OLEDマイクロキャビティ、誘電体ミラー、反射防止コーティング、分布ブラッグ反射器（DBR）、光学フィルター、導波路コーティング、その他の平面型フォトニックデバイスなどが挙げられる。

転送行列法はマクスウェル方程式から直接導出されるものであり、層状構造を有する媒質における電磁波方程式の一次元近似解を与える。多次元シミュレーション領域全体において完全なベクトル電磁場を解く代わりに、本手法では各層内の光を複素屈折率分布を通じて進行する進行波および後退波として表現する。

転送行列法は、光学的な問題を一連の簡潔な行列演算に還元するため、例えば完全波数数値解析手法などと比較して極めて高速である。 有限差分時間領域法転送行列法は計算的にも簡便であるにもかかわらず、定在波の形成、ファブリ・ペロー共鳴、光共振器の挙動、ならびに誘電体積層構造や多層光学コーティング内における電場の増強といったコヒーレントな薄膜干渉効果を正確に再現する。このため、多層薄膜デバイスにおけるコヒーレント光学現象のモデリング手法として最も広く用いられている手法の一つとなっている。

2. 転送行列法の物理的原理

転送行列法では、各層内の光電場を進行方向および逆方向に伝播する電磁波の重ね合わせとして表現する。

\[ E(z)=E^{+}e^{-jkz}+E^{-}e^{jkz} \]

ここで\(E^{+}\)および\(E^{-}\)はそれぞれ進行波および後退波の振幅を表し、\(z\)は多層構造内における波の伝播方向を示す。

\[ k=\frac{2\pi \tilde{n}}{\lambda}, \qquad \tilde{n}=n+j\kappa \]

これが複素光波数であり、ここで\(n\)は屈折率、\(\kappa\)は光吸収を表し、\(\lambda\)は真空中の波長である。??に示すように、前進および後退する波は各物質界面において部分的な反射および透過を通じて結合される。

この波の表現形式は、均一な媒質中におけるマクスウェル方程式から導出される一次元ヘルムホルツ波動方程式の解として得られるものである。

\[ \frac{d^2E(z)}{dz^2}+k^2E(z)=0 \]

各界面において、接線方向の電場および磁場の連続性を保つためにマクスウェルの電磁界境界条件が適用される。これらの条件からフレネル反射係数および透過係数が導出される。

\[ r=\frac{n_1-n_2}{n_1+n_2}, \qquad t=\frac{2n_1}{n_1+n_2} \]

ここで\(n_1\)および\(n_2\)は隣接する各層の複素屈折率を表す。これらの係数によって、各物質界面において前方および後方に伝播する波がどの程度反射・透過されるかが決定される。

3. 転送行列の構築

??に示すように、各層内の光電場は前方伝播波\(E^+\)および後方伝播波\(E^-\)によって表される。

\[ \mathbf{E}= \begin{bmatrix} E^{+} \\ E^{-} \end{bmatrix} \]

電磁波が厚さ\(d\)の層内を伝播する際、光電場は位相の変化および減衰を生じる。この層内での伝播過程は次のように表現できる：

\[ \begin{bmatrix} E^{+}(d) \\ E^{-}(d) \end{bmatrix} = \begin{bmatrix} e^{-jkd} & 0 \\ 0 & e^{jkd} \end{bmatrix} \begin{bmatrix} E^{+}(0) \\ E^{-}(0) \end{bmatrix} \]

二つの物質の界面において、光電場の一部は透過し、残りは反射される。この界面における関係式は行列形式でも表現可能である：

\[ \begin{bmatrix} E^{+}_2 \\ E^{-}_2 \end{bmatrix} = \frac{1}{t} \begin{bmatrix} 1 & r \\ r & 1 \end{bmatrix} \begin{bmatrix} E^{+}_1 \\ E^{-}_1 \end{bmatrix} \]

したがって、多層光学構造全体は一連の伝播行列および界面行列として表現される。

\[ \mathbf{E}_{out} = M_N I_N \cdots M_2 I_2 M_1 I_1 \mathbf{E}_{in} \]

ここで\(M\)は薄膜層内での光の伝播を表し\(I\)は界面における透過および反射を表す。こうして一連の単純な行列乗算を通じて、デバイス全体の光学的応答を算出することが可能となる。

この手法により、複雑な多層光学構造における反射、透過、光吸収、光干渉、定在波の形成を極めて効率的に計算することが可能となり、かつマクスウェル方程式と直接結びついている。転送行列法は光波方程式を直接解くため、薄膜における干渉効果を自然に捉えることができる。光が複数の界面で反射する際、各層内で蓄積される光位相に応じて反射波同士が建設的または破壊的に干渉する。ここで\( \phi = \frac{2\pi n d}{\lambda} \)は屈折率\(n\)、\(d\)は層の厚さ、\(\lambda\)は波長を表す。その結果、ナノメートルスケールの厚さ変化でさえもデバイスの光学的応答に大きな影響を与え得る。このような干渉現象が、反射防止コーティング、ファブリ・ペロー共振器、分布ブラッグ反射器、光学フィルター、OLEDマイクロキャビティ、薄膜太陽電池の動作原理となっている。

4. 転送行列法では何を計算するのか？

多層構造全体の転送行列が構築された後、デバイス内部の任意の位置における電磁場を算出することが可能となる。これにより、構造全体にわたる光場の空間分布や反射率、透過率、光吸収率を求めることができる。

転送行列法は光波方程式を直接解くため、シミュレーションにおいて光干渉効果が自然に現れる。界面間で光が何度も反射することで、建設的干渉および破壊的干渉が生じ、デバイス内部に定在波パターンが形成される。

??に示されている光子密度分布は、多層構造全体における光電場強度の変化状況を表している。光電場強度が高い領域は、光干渉や共振器効果によって光子密度が増大している箇所に相当する。

局所的な光電場分布が、デバイス内部で光子がどこで吸収されるかを直接的に決定する。電磁場と材料の光吸収係数を組み合わせることで、光生成分布を算出することが可能となる。

\[ G(z)=\alpha(z)\Phi(z) \]

ここで\(G(z)\)は局所的な電荷生成率、\(\alpha(z)\)は光吸収係数、\(\Phi(z)\)は局所的な光子束密度を表す。転送行列法を用いて計算された生成プロファイルの一例を??に示す。

これらの計算は、薄膜太陽電池、OLED、光学コーティング、誘電体ミラー、光センサー、および共振型マイクロキャビティデバイスにおいて広く利用されている。転送行列法は計算コストが低いため、全波長域にわたる光学スペクトルや内部電場分布をほんの数秒で算出することが可能である。

5. 転送行列法の利点と限界

転送行列法（TMM）は、物理的精度が高い一方で計算コストが極めて低いという特長から広く利用されている。多層光学構造においては、TMMは大規模な完全波シミュレーションと同等の薄膜干渉現象を再現可能でありながら、必要なメモリ量や計算時間はそれらに比べて桁違いに少ない。

TMMの主な利点としては以下の点が挙げられる：

• 極めて高速なシミュレーションが可能であり、パラメータスイープやデバイス最適化に最適である。
• 正確な薄膜干渉のシミュレーション：コヒーレント干渉や定在波効果を厳密に再現する。
• 直接的な光学的出力として、電界強度、光子密度、および光吸収分布を容易に算出できる。
• 平面状の多層構造の解析に最適であり、特にOLEDや太陽電池、光学コーティング、および分布ブラッグ反射鏡の解析に有用である。
• 電気的モデルとの連携が容易：光生成分布をドリフト・拡散シミュレーションに直接入力可能。

しかしながら、転送行列法にも重要な制限点が存在する。この手法は根本的に、光学構造が成長方向に沿ってのみ物性が変化する近似的に一次元的な構造と見なせることを前提としている。そのため、任意の三次元散乱構造を厳密にモデル化することは困難である。

• 主に一次元的な構造に限定され、多層構造でのみ適用可能である。
• 散乱現象の取り扱いが限定的であり、拡散散乱や表面の不規則な粗さは厳密には考慮されていない。
• 任意の幾何形状には対応不可：回折格子、ナノ粒子、テクスチャー加工された表面、およびフォトニック結晶の解析には通常、全波解析ソルバーが必要となる。
• コヒーレンスに関する仮定：標準的な転送行列法では各層内における光のコヒーレントな伝播が仮定されている。

三次元的な光学的問題においては、数値計算手法としては例えば FDTD法有限要素法（FEM）や厳密結合波解析（RCWA）の方が一般的に適している。したがってTMMは、干渉効果が重要であるものの構造がほぼ一次元的である平面多層デバイスに最も適している。

8. まとめ

転送行列法は薄膜光学において最も重要な数値計算手法の一つであり、多層構造内における反射、透過、光干渉、光吸収、および光生成の現象を迅速かつ物理的に正確にモデル化する手段を提供する。

この手法は主に多層構造に限定されるものの、平面型光電子デバイスのシミュレーションにおいて利用可能な計算効率に優れた手法の一つである。

実際のデバイスシミュレーションの作業フローにおいて、TMMは電気輸送解析ソルバーと併用されることが多く、光学的現象と半導体物理現象の相互作用を解析するために用いられる。