[논문 리뷰] Fast and Accurate Optimization of Metasurfaces with Gradient Descent and the Woodbury Matrix Identity

Jordan Budhu, Anthony Grbic|arXiv (Cornell University)|2021. 07. 07.

Advanced Antenna and Metasurface Technologies인용 수 4

한 줄 요약

이 논문은 유한 크기 메타표면의 빠르고 정확한 기울기 기반 최적화 방법을 제안한다. 여기서는 전자기 응답을 기술하기 위해 적분 방정식을 사용하는 Method of Moments(MoM)와 Woodbury 행렬 항등식을 활용한다. 기울기 계산 중에 변화된 행렬의 역행렬을 효율적으로 계산하기 위해 Woodbury 항등식을 적용함으로써, 표준 수치 미분 대비 최대 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.

연구 동기 및 목표

• 유한 크기 메타표면을 적분 방정식으로 모델링한 계산적으로 효율적인 최적화 프레임워크를 개발하기 위해.
• 메타표면 설계에서 수치 기울기 계산의 높은 계산 비용을 극복하기 위해.
• 행렬 항등식을 활용하여 기울기 평가 중 변화된 임피던스 행렬의 역행렬 계산을 가속화하기 위해.
• 높은 정확도와 단축된 런타임을 바탕으로 실용적이고 대규모의 메타표면 최적화를 가능하게 하기 위해.
• 수동적이고 손실이 없는 광역 반사 메타표면에 대해 이 방법을 적용하여 성능을 입증하기 위해.

제안 방법

• 메타표면은 전자기적 응답을 기술하기 위해 적분 방정식을 사용하는 Method of Moments(MoM)로 모델링된다.
• 기울기 계산은 유한 차분을 통해 수행되며, N개의 요소에 대해 MoM 임피던스 행렬의 역행렬을 총 N+1회 계산이 필요하다.
• 변화된 임피던스 행렬의 역행렬은 원래 역행렬에 대한 저랭크 보정으로서 Woodbury 행렬 항등식을 적용하여 계산된다.
• 이로 인해 각 기울기 성분의 계산 비용이 요소당 O(N³)에서 O(N²)로 감소하여 최적화 속도가 크게 향상된다.
• 가속화된 기울기 하강 알고리즘을 사용하여 원하는 파면 조작을 위한 메타표면 유닛 세포 기하 구조를 최적화한다.
• 이 방법은 수동적이고 손실이 없는 금속 패atters를 유전체 기판에 구현한 광역 반사 메타표면에 대해 검증되었다.

실험 결과

연구 질문

• RQ1 Woodbury 행렬 항등식은 메타표면 최적화에서 기울기 계산을 가속화하는 데 효과적으로 활용될 수 있는가?
• RQ2 행렬 항등식 기반의 역행렬 계산을 통해 기울기 기반 메타표면 최적화에서 얼마나 많은 성능 향상을 달성할 수 있는가?
• RQ3 가속화된 방법은 복잡한 파면 제어를 요구하는 수동적이고 손실이 없는 메타표면 설계에서 높은 정확도를 유지할 수 있는가?
• RQ4 계산 효율성과 수렴 속도 측면에서 이 방법의 성능 한계는 무엇인가?
• RQ5 이 방법은 유한 크기의 메타표면을 실용적으로 설계할 수 있도록 가능하게 하는가?

주요 결과

• Woodbury 행렬 항등식의 사용으로 각 기울기 성분의 계산 비용이 O(N³)에서 O(N²)로 감소하여 더 빠른 최적화가 가능해졌다.
• 가속화된 알고리즘이 표준 수치 기울기 계산 대비 최대 26.5배의 런타임 단축을 달성했다.
• 높은 정확도를 유지하여 수동적이고 손실이 없는 광역 반사 메타표면의 설계가 가능했다.
• 시간 측정도표는 다양한 메타표면 크기와 구성에서 일관된 성능 향상을 확인시켰다.
• 이 방법은 확장 가능하며, 복잡한 유닛 세포 기하 구조를 가진 대규모 메타표면 최적화에 적합하다.
• 이 방법은 동질화 근사에 의존하지 않고 직접적으로 유한 크기의 메타표면을 최적화할 수 있다.