[論文レビュー] Fast and Accurate Optimization of Metasurfaces with Gradient Descent and the Woodbury Matrix Identity

Jordan Budhu, Anthony Grbic|arXiv (Cornell University)|Jul 7, 2021

Advanced Antenna and Metasurface Technologies被引用数 4

ひとこと要約

本稿では、モーメント法（MoM）とウッドベリー行列恒等式を用いて、有限サイズのメタサーフェスの勾配ベース最適化手法を高速かつ高精度に提示する。勾配計算中に摂動を加えた行列の逆行列を効率的に計算するためにウッドベリー恒等式を適用することで、標準的な数値微分法に比べて最大26.5倍の高速化を達成し、被動的で損失のない広角反射メタサーフェスの効率的設計を可能にする。

ABSTRACT

A fast metasurface optimization strategy for finite-size metasurfaces modeled using integral equations is presented. The metasurfaces considered are constructed from finite patterned metallic claddings supported by grounded dielectric spacers. Integral equations are used to model the response of the metasurface to a known excitation and solved by Method of Moments. An accelerated gradient descent optimization algorithm is presented that enables the direct optimization of such metasurfaces. The gradient is normally calculated by solving the method of moments problem N+1 times where N is the number of homogenized elements in the metasurface. Since the calculation of each component of the N-dimensional gradient involves perturbing the moment method impedance matrix along one element of its diagonal and inverting the result, this numerical gradient calculation can be accelerated using the Woodbury Matrix Identity. The Woodbury Matrix Identity allows the inverse of the perturbed impedance matrix to be computed at a low cost by forming a rank-r correction to the inverse of the unperturbed impedance matrix. Timing diagrams show up to a 26.5 times improvement in algorithm times when the acceleration technique is applied. An example of a passive and lossless wide-angle reflecting metasurface designed using the accelerated optimization technique is reported.

研究の動機と目的

• 積分方程式を用いてモデル化された有限サイズのメタサーフェスの計算的に効率的な最適化フレームワークの開発。
• メタサーフェス設計における数値勾配計算の高い計算コストを克服すること。
• 行列恒等式を活用して、勾配評価中の摂動を加えたインピーダンス行列の逆行列を高速化すること。
• 高い精度と短い実行時間で、実用的で大規模なメタサーフェスの最適化を可能にすること。
• 被動的で損失のない広角反射メタサーフェスへの応用を示すこと。

提案手法

• メタサーフェスは、電磁応答を記述する積分方程式を用いたモーメント法（MoM）でモデル化される。
• 勾配計算は有限差分法で実行され、N個の要素に対してはMoMインピーダンス行列のN+1回の逆行列計算が必要となる。
• ウッドベリー行列恒等式を用いて、摂動を加えたインピーダンス行列の逆行列を元の逆行列に対する低ランク補正として計算する。
• これにより、各勾配成分の計算コストをN個の要素あたりO(N³)からO(N²)に削減し、最適化を著しく高速化する。
• 高速化された勾配降下法アルゴリズムを用いて、所望の波フロント制御を実現するためのメタサーフェスユニットセルの幾何形状を最適化する。
• この手法は、誘電体基板上に金属パターンを配置した被動的で損失のない広角反射メタサーフェスに対して検証された。

実験結果

リサーチクエスチョン

• RQ1 ウッドベリー行列恒等式は、メタサーフェス最適化における勾配計算の高速化に効果的に適用可能か？
• RQ2 行列恒等式に基づく逆行列計算を用いることで、勾配ベースのメタサーフェス最適化でどの程度の高速化が達成できるか？
• RQ3 高速化された手法は、複雑な波フロント制御を実現する被動的で損失のないメタサーフェスの設計においても高い精度を維持できるか？
• RQ4 計算効率と収束速度の観点から、この手法の性能限界は何か？
• RQ5 この手法は、広角応答を示す有限サイズのメタサーフェスの実用的設計を可能にするか？

主な発見

• ウッドベリー行列恒等式の使用により、各勾配成分の計算コストがO(N³)からO(N²)に削減され、最適化が高速化された。
• 高速化されたアルゴリズムは、標準的な数値勾配計算に比べて最大26.5倍の実行時間短縮を達成した。
• 高い精度を維持しており、被動的で損失のない広角反射メタサーフェスの設計が可能となった。
• タイミング図により、異なるメタサーフェスサイズや構成において一貫した性能向上が確認された。
• スケーラブルであり、複雑なユニットセル幾何形状を有する大規模メタサーフェス最適化に適している。
• ホモジゼーション近似に依存せずに、有限サイズのメタサーフェスの直接最適化が可能となった。