[论文解读] 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）方法，可精确预测单层及多层纳米结构的光学性质。该方法在三个案例研究中均与实验椭偏测量结果高度一致——从单层 Ag 膜到复杂的非周期性超材料——证明了其在实际制造前设计先进纳米光子材料方面的鲁棒性。

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-极化光分别进行极化特异性模拟：在二维中，s-极化为 E = (0,0,1)，p-极化为 H = (0,0,1)；在三维中，分别为 E = (0,1,0) 和 H = (0,1,0)。
• 计算复数菲涅尔反射系数 rp 和 rs，进而通过 ρ = tan(Ψ)e^{iΔ} = rp/rs 推导出椭偏参数 Ψ 和 Δ。
• 在 x 和 y 方向应用弗洛quet周期性边界条件（PBCs），以降低计算成本，并在横向平面上模拟无限周期性。
• 网格划分采用物理控制方式，在关键区域采用更细的分辨率，而在体材料区域采用较粗的网格，以在保证精度的同时优化计算时间。
• 采用两个并行的频域研究，同时求解 s- 和 p-极化波，从而提高仿真效率。

实验结果

研究问题

• RQ1 是否能够通过数值仿真平台在不依赖解析模型的前提下，准确预测单层及多层纳米结构的椭偏响应？
• RQ2 NEA 方法在结构复杂度逐渐增加的系统中，对反射率、透射率以及椭偏角（Ψ 和 Δ）的再现能力如何？
• RQ3 使用周期性边界条件和物理控制网格划分在多大程度上提升了计算效率，同时不牺牲精度？
• RQ4 仿真得到的光学响应与实际制造的 Ag、Ag/ITO 和 AZO/ITO/Ag 纳米结构的实测结果相比如何？
• RQ5 NEA 框架是否可作为新型纳米光子材料在昂贵制造前的可靠设计工具？

主要发现

• 对于玻璃基底上的 20 nm Ag 膜，NEA 模型与实验数据高度一致，仅因模拟中基底尺寸减小导致反射率存在约 20% 的振幅偏差。
• 对于由五层 Ag/ITO 双层结构组成的超材料（每层总厚度 40 nm），模拟与实测的 R_p、R_s、T、Ψ 和 Δ 曲线在整个光谱范围内均表现出惊人的一致性。
• 在具有非均匀层厚的非周期性三组分单元结构（AZO/ITO/Ag）中，NEA 模型成功预测了复杂光学特性，尽管缺乏周期性。
• 在 λ = 390 nm 处的电场分布图证实了在超材料中存在强场局域化和各向异性行为，与介电常数接近零且具有极点（εNZP）特性一致。
• NEA 方法准确捕捉到了 Ag 膜中约 330 nm 处的 Ferrel-Berreman 模式，其特征为反射率最小值与透射率最大值。
• 使用 PBCs 和并行频域研究显著缩短了计算时间，同时在复杂多层系统仿真中保持了高精度。