[論文レビュー] A numerical ellipsometric analysis (NEA) of nanoscale layered systems

Giuseppe Emanuele Lio, Giovanna Palermo+2 more|arXiv (Cornell University)|Dec 12, 2019

Optical Coatings and Gratings参考文献 57被引用数 12

ひとこと要約

本論文は、COMSOL Multiphysicsを用いた数値的楕円偏光測定法（NEA）を提案し、単層および多層ナノスケール系の光学的性質を正確に予測することを目的としている。この手法は、銀単層から複雑な非周期的メタマテリアルに至るまでの3つの事例において、実験的楕円偏光測定と優れた一致を示し、製造前の高度なナノフォトニクス材料の設計における堅牢性を実証している。

ABSTRACT

A simple and robust method able to predict, with high accuracy, the optical properties of single and multi-layer nanostructures is presented. The method exploits a COMSOL Multiphysics simulation platform and it has been validated by three case studies with increasing numerical complexity: i) a single thin layer (20 nm) of Ag deposited on a glass substrate; ii) a metamaterial composed of five bi-layers of Ag/ITO (Indium Tin Oxide), with a thickness of 20 nm each; iii) a system based on a three-materials unit cell (AZO/ITO/Ag), but without any thickness periodicity (AZO stands for Al$_2$O$_3$/Zinc Oxide). Numerical results have been compared with experimental data provided by real ellipsometric measurements performed on the above mentioned nanostructures ad-hoc fabricated. The obtained agreement is excellent suggesting this research as a valid approach to design materials able to work in a broad spectrum range.

研究の動機と目的

• 実験的製造の前に、複雑なナノスケール積層系の光学的応答を予測する信頼性の高い数値的手法を開発すること。
• 高度に複雑なナノ構造において解析的手法の限界を克服するため、数値シミュレーションを活用すること。
• 高価な試行錯誤実験を減らすために、機能的ナノ材料の正確な製造前の設計を可能にすること。
• 構造的複雑性の異なる実際のサンプルに対する楕円偏光測定データとシミュレーションツールを照合すること。
• 光学的性質を調整可能なメタマテリアルおよび多層膜の設計に、包括的かつカスタマイズ可能なフレームワークを提供すること。

提案手法

• NEA手法は、COMSOL Multiphysicsを用いて、周波数領域におけるマクスウェル方程式を、積層構造内を通過する電磁波の伝播に対して解く。
• s偏光およびp偏光の両方の偏光状態に対して、適切な場の配置を設定することで、偏光に特化したシミュレーションを実施：2次元ではs偏光に対してE = (0,0,1)、p偏光に対してH = (0,0,1)を設定し、3次元ではE = (0,1,0)およびH = (0,1,0)を設定する。
• 複素数のフレネル反射係数rpとrsを計算し、ρ = tan(Ψ)e^{iΔ} = rp/rsの式を用いて、楕円偏光パラメータΨとΔを導出する。
• 計算コストを削減し、横方向に無限大の周期性を模擬するために、xおよびy方向にフロケット周期的境界条件（PBC）を適用する。
• メッシュ生成は物理的制御に基づき、重要な領域では細かく、バルク領域では粗くすることで、精度を確保しながら計算時間を最適化する。
• s偏光およびp偏光の波を同時に解く2つの周波数領域研究を並列で実行することで、シミュレーションの効率を向上させる。

実験結果

リサーチクエスチョン

• RQ1 解析的手法に依存せずに、単層および多層ナノスケール系の楕円偏光応答を正確に予測できる数値シミュレーションプラットフォームは存在するか？
• RQ2 構造的複雑性が増す系において、NEA手法は実験的反射率、透過率、および楕円偏光角（ΨおよびΔ）をどの程度正確に再現するか？
• RQ3 周期的境界条件および物理的制御メッシュの使用が、精度を損なわせることなく計算効率をどの程度向上させるか？
• RQ4 シミュレートされた光学的応答は、銀、Ag/ITO、AZO/ITO/Agナノ構造の実際の製造済みサンプルとの測定結果とどの程度一致するか？
• RQ5 NEAフレームワークは、高価な製造を伴う新しいナノフォトニクス材料の設計に信頼できるツールとして機能できるか？

主な発見

• ガラス基板上に形成された20 nmの銀膜に対して、NEAモデルは実験データと優れた一致を示し、基板サイズの縮小に起因するわずかな約20%の振幅差異を除いて良好な一致を示している。
• 5層のAg/ITOバイレイヤー（1層あたり40 nm）からなる超ひずみメタマテリアル（HMM）において、シミュレートされたR_p、R_s、T、Ψ、Δの曲線が、全スペクトルにわたり実測値と驚くほど一致している。
• 周期性の欠片な3成分のユニットセル（AZO/ITO/Ag）で、層厚が非一様な構造においても、NEAモデルは複雑な光学的特徴を正確に予測している。
• λ = 390 nmにおける電場マップは、HMMで強い場の閉じ込めと異方的挙動が確認されており、ε近零・極（εNZP）特性と整合的である。
• Ag膜におけるフェレル＝ベルマンモード（約330 nm付近）は、反射率の最小値と透過率の最大値として明確に特定されており、NEA手法により正確に捉えられている。
• PBCと並列処理された周波数領域研究の使用により、複雑な多層系のシミュレーションにおいて、計算時間を顕著に短縮しながらも高い精度を維持している。