Expanding RCWA capabilities with advanced S-matrix algorithms

光は、何もない空間を移動する際には予測可能な振る舞いを見せますが、薄い膜の積層や微細な構造を持つ表面に遭遇した瞬間、そのルールは複雑になります。これは、反射防止コーティング、ミラー、そして光ファイバーネットワーク内の複雑なコンポーネントといった現代の光学デバイスのバックボーンである、層状フォトニック構造の世界です。光がこれらの層に当たると、あらゆる境界で反射、干渉、散乱を起こし、予測が困難なエネルギーの入り組んだ網を作り出します。より優れたデバイスを設計するために、エンジニアは光がこれらのスタックを通過する際にどのように振る舞うかを正確に計算する方法を必要としています。数十年にわたり、科学者たちはこの問題を解決するための強力な数学的ツールである「散乱行列法（scattering-matrix method）」に頼ってきました。この手法は、非常に厚い層や光が急速に減衰する状況においても安定しており、精度が高いことから重宝されています。こうした状況は、古い計算手法ではしばしば失敗してしまう場面です。しかし、このツールは標準的な問題にはうまく機能してきましたが、複雑な設計の探求や、光の流れの深い物理学を理解するための潜在能力は、まだ完全には引き出されていませんでした。

フランスのボルドー大学とインスティテュ・ド・オプティークの研究チームは、このツールを改良し、設計と発見の両面において、より柔軟で強力な計器へと進化させました。彼らは、散乱行列の定義方法を刷新し、焦点をつなぎ合わせる対象を「2つの層間の接続」から「単一の層自体の特性」へと移行させました。従来のアプローチでは、特定の層の計算は、その直上および直下に位置する材料に紐付けられていたため、もしエンジニアが同じ層を異なるスタックの中でテストしたいと考えた場合、すべてを一から再計算しなければなりませんでした。新しい手法は、層を、隣接する層がどうあれ、光がそこを通過する際にどのように扱うかによって記述される、独自の固有の振る舞いを持つ独立した実体として扱います。層が接する界面は、別途処理されます。この分離により、層の特性を一度計算すれば、あらゆる組み合わせにおいて際限なく再利用できるようになります。これは、新しい設計のために毎回作り直す必要のない、あらゆる壁に使用できる単一のレンガのようなものです。

この視点の転換は、単に計算速度を上げるだけでなく、これらの構造の内部における光の隠れた経路をも明らかにします。この再定式化されたアプローチを用いることで、研究者たちは、外部からの光波が、周期構造内に自然に存在する内部波（ブロッホモードとして知られる）とどのように結合するかを直接計算できるようになりました。これらのモードは、ギターの弦が特定の周波数でしか振動できないのと同様に、パターン化された材料の中を光が移動できる特定の様式です。新しい手法は、入射ビームがこれらの内部モードへとどれだけのエネルギーを運び込み、境界においてそれらがどのように相互作用するかを示す、地図のような役割を果たす一連の一般化係数を提供します。これを実証するために、チームは幅広い色に対して高い透明性を持つように設計されたシリコン・グレーティングを分析しました。従来の解説では個々のナノワイヤの共鳴に焦点を当てていましたが、新しい分析によれば、その透明性は実際には、わずか2つの特定の内部光モードの励起と伝搬によって駆動されていることが示されました。一つのモードは長い波長において支配的な挙動を示し、もう一つの非常に珍しい、有効速度が極めて低いモードが、短い波長における反射防止効果を担っていました。このレベルの詳細さにより、設計者はどの内部チャネルがエネルギーを運んでいるのかを正確に把握できるようになり、シミュレーションを、単に数値を出力する「ブラックボックス」から、デバイスの物理学への「透明な窓」へと変貌させたのです。

この再定式化による実用的なメリットは、特に設計最適化の速度において同様に重要です。新しい手法では、層の厚さを解析的な公式内の単純な変数として扱うため、以前はたった一つの計算を行うのに要した時間で、数百種類の異なる厚さに対する光学応答を計算できます。研究者たちは、上部のシリコン・グレーティングの厚さを400通りの値で変化させながら、シリカ・スペーサーで隔てられた2つのシリコン・グレーティングを持つ構造をシミュレートすることで、これをテストしました。この全スキャンは、新しいパラメータごとに完全な再実行を必要とする従来のシミュレーション手法では不可能であった時間を、ごくわずかな時間で完了しました。この能力はまた、光のコヒーレンス（干渉性）が失われる厚い材料を扱うための巧妙な解決策も提供します。太陽電池などの実世界のアプリケーションにおける一般的な問題であるランダム性を扱うために、複雑な新しい理論を構築する代わりに、多くの厚さの値にわたって結果を平均化するだけで、干渉による急速な変動を効果的に平滑化し、非コヒーレントな挙動の正確な姿を描き出せることをチームは示しました。

さらに、この新しいフレームワークにより、構造に当たる光を多くの異なる角度から同時に分析することが可能になります。従来のソフトウェアは通常、単一の光ビームの結果を計算するため、関心のある角度ごとに個別の実行を必要とします。しかし、更新されたアプローチでは、あり得るすべての入射方向および回折方向に対する散乱係数を、一度のステップで計算します。これは、人工的な周期性を用いてモデル化された、大規模で無秩序な表面を分析する際に特に有用であり、一度の計算で様々な照明条件下でのシステム全体の挙動を明らかにすることができます。これらすべての機能は、「RETICOLO」と呼ばれるフリーソフトウェアパッケージに統合され、著者らはこれを科学コミュニラティに公開しました。研究者たちは、結果を生成するために使用したコードと詳細なスクリプトを提供することで、他の科学者がこれらの高度なモード解析を行い、設計ワークフローを加速させ、複雑な層状構造を通じて光がどのように移動するかという、より深い物理学的理解を得るための実践的なツールキットを提供しています。