[논문 리뷰] Fullwave design of cm-scale cylindrical metasurfaces via fast direct solvers

Wenjin Xue, Hanwen Zhang외 3명|arXiv (Cornell University)|2023. 08. 15.

Metamaterials and Metasurfaces Applications인용 수 7

한 줄 요약

이 논문은 센티미터 규모의 원통형 메타표면에 대한 전체파 해석 시뮬레이션을 위해 설계된 빠른 직접 적분방정식 해법을 제시한다. 이는 메타표면 지름에 대해 선형 스케일링을 달성하며, 정확하고 고해상도의 대규모 나노광학 장치 설계를 가능하게 한다. 이는 1cm 지름, 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 파장에서 1cm 지름의 금속렌즈와 1550 nm 파장에서 넓은 빔 폭을 가진 격자 커플러.
• 다양한 설계 반복 과정에서 반복적으로 사용 가능한 압축된 시스템 행렬을 제공하여 반복적 설계 프로세스(예: 인접 최적화)의 속도를 높인다.

제안 방법

• FDTD와 같은 미분 방법에서 유래하는 악조건화 및 위상 오차 문제를 피하기 위해 전자기 산란에 대한 잘 조절된 적분방정식 수식을 사용한다.
• 고차수 적분법을 적용하여 낮은 해상도(λ/40)에서도 높은 정확도를 확보하고, 이산화 오차를 감소시킨다.
• 계층적 행렬 압축과 그린 함수 상호작용의 저랭크 압축을 통해 빠른 직접 해법을 구현하여 메타표면 지름에 대해 준선형 스케일링을 달성한다.
• 원통형 메타표면의 이동 불변 기하학적 특성(2차원 x–z 평면)을 활용하여 해법의 알고리즘 구조를 최적화하여 속도와 메모리 효율을 향상시킨다.
• 인접 최적화 루프 내에서 전방 및 인접 장파면 자극에 해법을 적용하여 복잡한 메타표면 패tern의 기울기 기반 설계를 가능하게 한다.
• 다양한 설계 반복 과정에서 압축된 시스템 행렬을 재사용하여 각 최적화 단계의 계산 비용을 크게 감소시킨다.

실험 결과

연구 질문

• RQ1 빠른 직접 적분방정식 해법이 센티미터 규모의 장치에서 높은 정확도를 유지하면서도 메타표면 지름에 대해 선형 스케일링을 달성할 수 있는가?
• RQ2 이러한 해법이 대규모 시뮬레이션에서 기존 유한차분 방법의 위상 누적 및 악조건화 문제를 어느 정도 해결할 수 있는가?
• RQ3 이 방법을 통해 전체파 최적화를 통해 도달할 수 있는 최대 크기와 효율은 무엇인가? 특히 고수치 군집도 렌즈의 경우에 대해.
• RQ4 이 방법은 넓은 빔 프로파일을 가진 고효율 격자 커플러와 같은 대면적 빔 형상 조절 장치의 효율적인 설계를 가능하게 하는가?
• RQ5 대규모 메타표면에 대해 기존 FDTD 대비 이 빠른 직접 해법의 정확도, 확장성, 계산 비용 측면에서의 성능 비교는 어떠한가?

주요 결과

• 빠른 직접 해법은 해법 시간이 메타표면 지름에 대해 선형(준선형)으로 증가하여 기존 FDTD로는 비가능한 센티미터 규모 장치의 시뮬레이션을 가능하게 한다.
• 500 nm 파장에서 1cm 지름의 금속렌즈(20,000λ 지름)가 전체파 인접 최적화를 통해 성공적으로 설계되었으며, 이는 대규모 나노광학 장치 설계의 가능성과 실현 가능성을 입증한다.
• 격자 커플러 설계는 대규모 면적의 가우시안 빔(1550 nm에서 FWHM ≈ 800λ)으로 53%의 커플링 효율을 달성하였으며, 대칭적 복사 패턴으로 인해 이론적 한계에 근접하였다.
• 시스템 행렬의 압축 단계는 단지 30초가 소요되었고, 이는 모든 최적화 반복 과정에서 재사용 가능하여 설계 프로세스의 속도를 크게 향상시켰다.
• 각 최적화 반복 단계에서의 전방 및 인접 시뮬레이션은 8코어 노드에서 2분 이상 소요되었으며, 이는 약 33시간 내에 750회의 반복 최적화를 가능하게 하였다.
• 최적 설계에서 유도된 전자기장 패턴은 목표 장파면과 밀도적으로 유사하였으며, 대칭 구조에서 피할 수 없는 하향 복사로 인한 미세한 편차가 존재하였다.