[論文レビュー] 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 cm 直径で 20,000λ の高 NA メタレンズおよび 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 cm 直径のメタレンズと、1550 nm 波長での大ビーム幅グレーティングカップラー。
• 反復的設計プロセス（例：アドジョイント最適化）を加速するために、再利用可能な圧縮済みシステム行列を提供する。

提案手法

• FDTD などの微分法に内在する悪条件性と位相誤差を回避するため、電磁散乱の良好に条件付けられた積分方程式定式化を採用する。
• 高次の台形則を用いて、中程度の解像度（λ/40）で高い精度を達成し、離散化誤差を低減する。
• 階層行列圧縮とグリーン関数相互作用の低ランク圧縮を組み合わせた高速直接ソルバを実装し、メタサーフェスの直径に対して準線形スケーリングを達成する。
• 円筒型メタサーフェスの並進不変性（x–z 平面における 2D 構造）を活用して、アルゴリズムの構造を最適化し、速度とメモリ効率を向上させる。
• アドジョイント最適化ループ内で前方およびアドジョイント場励起にソルバを適用し、複雑なメタサーフェスパターンの勾配ベース設計を可能にする。
• 複数の設計反復においても再利用可能な圧縮済みシステム行列を活用し、最適化ステップあたりの計算コストを大幅に削減する。

実験結果

リサーチクエスチョン

• RQ1 高速直接積分方程式ソルバは、センチメートルスケールのデバイスにおいても、メタサーフェスの直径に対して線形スケーリングを達成しながら、高精度を維持できるか？
• RQ2 このようなソルバは、大規模シミュレーションにおいて従来の有限差分法が抱える位相蓄積と悪条件性の問題をどの程度克服できるか？
• RQ3 本手法を用いたフルウェーブ最適化によるメタサーフェス設計の最大サイズと効率はどの程度達成可能か、特に高数値孔径レンズにおいては？
• RQ4 本手法は、広いビームプロファイルを持つ高効率グレーティングカップラーのような大面積ビーム形成素子の効率的設計を可能にするか？
• RQ5 大規模メタサーフェスにおいて、従来の FDTD と比較して、本高速直接ソルバの正確性、スケーラビリティ、計算コストの観点での性能はどのようになるか？

主な発見

• 高速直接ソルバは、メタサーフェスの直径に対して解の時間計算が線形（準線形）にスケーリングし、従来の FDTD では不可能なセンチメートルスケールのデバイスシミュレーションを可能にした。
• 500 nm 波長の 1 cm 直径のメタレンズ（直径 20,000λ）が、フルウェーブアドジョイント最適化により成功裏に設計され、大規模ナノフォトニクス素子設計の実現可能性が示された。
• グレーティングカップラー設計では、1550 nm 波長で 53% の結合効率を達成し、対称的放射パターンのおかげで理論限界に近い性能を示した（FWHM ≈ 800λ の大面積ガウスビームへ）。
• システム行列の圧縮段階はわずか 30 秒で完了し、全最適化反復において再利用可能であり、設計プロセスの大幅な加速に寄与した。
• 最適化反復ごとの前方およびアドジョイントシミュレーションは 8 コアノードで 2 分をわずかに超え、約 33 時間で 750 回の最適化ステップを実行できた。
• 最適設計からの場の分布は、目標とする場とよく一致しており、対称的構造に避けがたい下方放射が生じるため、わずかなずれが見られた。