[论文解读] Fullwave design of cm-scale cylindrical metasurfaces via fast direct solvers

Wenjin Xue, Hanwen Zhang+3 more|arXiv (Cornell University)|Aug 15, 2023

Metamaterials and Metasurfaces Applications被引用 7

一句话总结

本文提出了一种针对厘米尺度圆柱形超表面全波仿真而定制的快速直接积分方程求解器，实现了与超表面直径的线性缩放。该方法可实现高精度、高保真度的大规模纳米光子器件设计——通过伴随优化结合全麦克斯韦方程组仿真，成功实现了1厘米直径、20,000λ高数值孔径金属透镜以及1000λ光栅耦合器的设计。

ABSTRACT

Large-scale metasurfaces promise nanophotonic performance improvements to macroscopic optics functionality, for applications from imaging to analog computing. Yet the size scale mismatch of centimeter-scale chips versus micron-scale wavelengths prohibits use of conventional full-wave simulation techniques, and has necessitated dramatic approximations. Here, we show that tailoring "fast direct" integral-equation simulation techniques to the form factor of metasurfaces offers the possibility for accurate and efficient full-wave, large-scale metasurface simulations. For cylindrical (two-dimensional) metasurfaces, we demonstrate accurate simulations whose solution time scales \emph{linearly} with the metasurface diameter. Moreover, the solver stores compressed information about the simulation domain that is reusable over many design iterations. We demonstrate the capabilities of our solver through two designs: first, a high-efficiency, high-numerical-aperture metalens that is 20,000 wavelengths in diameter. Second, a high-efficiency, large-beam-width grating coupler. The latter corresponds to millimeter-scale beam design at standard telecommunications wavelengths, while the former, at a visible wavelength of 500 nm, corresponds to a design diameter of 1 cm, created through full simulations of Maxwell's equations.

研究动机与目标

• 解决由于传统有限差分/有限元方法固有的病态性和相位误差，导致厘米尺度超表面全波麦克斯韦仿真计算不可行的问题。
• 克服微米量级波长与厘米尺度器件之间的尺寸尺度不匹配问题，该问题限制了标准仿真工具的应用。
• 通过基于良好条件积分方程的定制化快速直接求解器，实现大面积极超表面的高精度、可扩展全波设计。
• 通过高性能设计展示该方法的能力：在500 nm波长下实现1厘米直径金属透镜，以及在1550 nm波长下实现大光束宽度光栅耦合器。
• 提供可重用的压缩系统矩阵，以加速伴随优化等迭代设计过程。

提出的方法

• 采用具有良好条件的电磁散射积分方程公式，避免有限差分法（如FDTD）固有的病态性和相位误差。
• 使用高阶积分方案，在中等分辨率（λ/40）下实现高精度，减少离散化误差。
• 通过层级矩阵压缩和格林函数相互作用的低秩压缩，实现快速直接求解器，使求解时间与超表面直径呈准线性缩放。
• 利用圆柱形超表面（x–z平面二维）的平移不变几何结构，优化求解器的算法结构以提升速度和内存效率。
• 在伴随优化循环中应用求解器进行正向场和伴随场激励，实现基于梯度的复杂超表面图案设计。
• 在多个设计迭代中重用压缩后的系统矩阵，显著降低每次优化步骤的计算成本。

实验结果

研究问题

• RQ1 能否将快速直接积分方程求解器适配，实现与超表面直径的线性缩放，同时保持厘米尺度器件的高精度？
• RQ2 此类求解器在大规模仿真中，能在多大程度上克服传统有限差分方法普遍存在的相位累积和病态性问题？
• RQ3 采用该方法可实现的最大尺寸和效率的全波优化超表面设计是什么？特别是针对高数值孔径透镜？
• RQ4 该方法能否实现大面积极光束整形器件（如具有宽光束轮廓的高效率光栅耦合器）的高效设计？
• RQ5 在精度、可扩展性和计算成本方面，该快速直接求解器与传统FDTD相比表现如何，特别是在大尺寸超表面仿真中？

主要发现

• 快速直接求解器的求解时间与超表面直径呈线性（准线性）缩放，使原本在标准FDTD中不可行的厘米尺度器件仿真成为可能。
• 成功通过全波伴随优化设计出1厘米直径、500 nm波长的金属透镜（直径达20,000λ），证明了大规模纳米光子器件设计的可行性。
• 光栅耦合器设计实现了53%的耦合效率，将光耦合至大尺寸高斯光束（1550 nm下FWHM ≈ 800λ），接近理论极限，归因于对称辐射图案。
• 系统矩阵的压缩阶段仅耗时30秒，且可在所有优化迭代中重复使用，显著加速了设计流程。
• 每次优化迭代的正向与伴随仿真耗时略多于两分钟（在8核节点上），使750次迭代的优化过程总耗时约33小时。
• 最优设计的场分布与目标场高度吻合，微小偏差源于对称结构中不可避免的向下辐射。