Adjoint-Based Gradient Evaluation for Metasurface Inverse Design via Affine Geometric Transformations

1. 연구 배경 및 문제 제기 (Problem)

• 메타표면의 복잡성과 설계 난제: 메타표면 (Meta-surface) 은 파장 크기 이하의 산란체 (메타-아tom) 가 대량으로 배열된 구조로, 빔 조향, 광 집속 등 다양한 전자기 제어가 가능합니다. 그러나 대규모 구조를 설계할 때, 수백만 개의 자유도를 가진 전자기 역문제 (Inverse Problem) 를 해결하는 것은 계산 비용이 매우 큽니다.
• 기존 방법론의 한계:
  • 유닛 셀 (Unit-cell) 접근법: 사전에 정의된 유닛 셀 라이브러리를 기반으로 설계하지만, 고 수치 개구율 (High NA) 이나 급격한 구조 변화가 필요한 경우 성능이 제한적입니다.
  • 자유형 (Freeform) 최적화: 유전율 분포를 연속적으로 최적화하는 방식은 성능은 우수하지만, 제조 가능한 구조로 이진화 (Binarization) 하는 과정이 필요하여 계산 복잡도가 높고 대규모 메타표면으로 확장하기 어렵습니다.
  • 기울기 (Gradient) 계산의 병목: 경사하강법 (Gradient-based optimization) 을 사용할 때, 설계 변수가 많으면 유한 차분법 (Finite Difference) 으로 기울기를 계산하는 것은 비현실적입니다 (각 변수마다 직접 시뮬레이션이 필요함).

2. 제안된 방법론 (Methodology)

이 논문은 어드저인트 변수법 (Adjoint Variable Method, AVM) 을 메타표면 역설계에 적용하여, 아핀 기하학적 변환 (Affine Geometric Transformations) 을 통해 설계 변수를 효율적으로 업데이트하는 체계적인 방법론을 제시합니다.

• 매개변수화 전략 (Parametrization):

  • 자유형 최적화 대신, 각 메타-아tom 의 기하학적 매개변수 (위치, 회전, 크기 등) 를 설계 변수로 정의합니다.
  • 각 메타-아tom 은 아핀 변환 (Affine Transformation) 을 통해 변형됩니다. 이는 다음 세 가지 요소를 포함합니다:
    1. 회전 (Rotation): $\theta$
    2. 이방성 확대/축소 (Anisotropic Expansion): $x, y$ 축 방향의 스케일링 $\lambda_x, \lambda_y$
    3. 이동 (Translation): $x, y$ 좌표 이동
  • 이 접근법은 설계 변수의 수를 줄여 계산 효율성을 높이면서도, 다양한 구조 변형을 허용하여 설계 공간을 넓게 탐색할 수 있게 합니다.

• 어드저인트 기반 기울기 계산:

  • 목적 함수 (Cost Function) $I$의 기울기 $\nabla_p I$를 계산하기 위해, 한 번의 직접 문제 (Direct Problem) 와 한 번의 어드저인트 문제 (Adjoint Problem) 만을 해결하면 모든 설계 변수에 대한 기울기를 얻을 수 있음을 증명합니다.
  • 기울기 공식:
    $$\frac{\partial I}{\partial p_j} = 2 \text{Re} \left[ \beta \int_{\partial \Omega_i} \mathbf{E}_a(\mathbf{x}) \cdot \mathbf{E}_f(\mathbf{x}) w(\mathbf{x}) \, dl \right]$$
    • $\mathbf{E}_f$: 전방 (Forward) 전계 (입사파에 대한 응답)
    • $\mathbf{E}_a$: 어드저인트 전계 (출력면의 목표장에 기반한 소스에 대한 응답)
    • $w$: 기하학적 변환 (회전, 확대, 이동) 에 의해 결정되는 가중치 함수.
    • $\beta$: 목적 함수에 의존하는 스칼라 계수.
  • 핵심 이점: 목적 함수와 기하학적 변환이 분리되어 있어, 변환 유형이 바뀌더라도 추가적인 직접 시뮬레이션 없이 가중치 $w$만 재계산하면 됩니다.

• 수치 모델:

  • 2D TEz 전파 문제를 가정하고, 경계 요소법 (Boundary Element Method, BEM) 을 사용하여 전자기장을 해석합니다.
  • 메타-아tom 의 경계에서 전계와 자계의 접선 성분을未知수로 두고, 어드저인트 필드 계산을 위해 필요한 법선 성분의 전계는 경계에서의 자계 분포로부터 유도합니다.

3. 주요 기여 (Key Contributions)

1. 일반적인 아핀 변환에 대한 기울기 유도: 기존 연구가 주로 등방성 확대나 단순한 형태 변화에 국한되었던 것과 달리, 회전, 이방성 스케일링, 이동을 포함한 일반적인 아핀 변환에 대한 기울기 식을 엄밀하게 유도했습니다.
2. 계산 효율성과 확장성: 대규모 메타표면 (수백 개 이상의 메타-아tom) 에 대해, 유한 차분법에 비해 계산 비용을 획기적으로 줄이면서도 복잡한 구조 설계를 가능하게 하는 프레임워크를 제시했습니다.
3. 다양한 목적 함수에 대한 어드저인트 소스 도출: 전계 크기, 전계 차이 노름, 전계 간 각도 등 다양한 목적 함수에 대해 필요한 어드저인트 소스 (Adjoint Source) 를 명시적으로 도출하여 실제 적용 가능성을 높였습니다.
4. 제조 가능성 고려: 매개변수화된 기하학적 접근을 통해, 최적화 과정이 자연스럽게 제조 가능한 현실적인 구조 (불규칙한 패턴 방지) 로 수렴하도록 유도했습니다.

4. 실험 결과 및 검증 (Results)

• 기울기 정확도 검증:
  • 유한 차분법 (Finite Difference) 으로 계산한 기울기와 제안된 어드저인트 방법으로 계산한 기울기를 비교했습니다.
  • 회전 각도, 메타-아tom 폭 (확대/축소), 중심 좌표 (이동) 등 다양한 설계 변수와 목적 함수 조합에서 두 방법의 결과가 높은 정확도로 일치함을 확인했습니다.
• 역설계 예시 (Light Focusing):
  • 단일 초점: 128 개의 메타-아tom 으로 구성된 메타렌즈를 설계하여, 특정 지점 ($37\lambda_0$) 에 빛을 집중시키는 데 성공했습니다. 최적화된 구조는 비대칭적이고 급격히 변화하는 회전 각도 패턴을 보였으며, 초점의 FWHM(반치폭) 은 $0.84\lambda_0$로 매우 날카로운 집속을 구현했습니다.
  • 다중 초점: 두 개 및 세 개의 초점 지점에 동시에 빛을 집중시키는 복잡한 설계도 수행되었습니다. 제안된 방법이 다중 스팟 (Multi-spot) 집속과 같은 복잡한 작업에서도 높은 정확도로 작동함을 입증했습니다.

5. 의의 및 결론 (Significance)

• 대규모 메타표면 설계의 실용화: 이 연구는 계산 비용의 병목 현상을 해결하여, 대규모 메타표면의 역설계를 실용적으로 가능하게 합니다.
• 유연한 설계 공간: 아핀 변환을 통해 단순한 크기 조절을 넘어 회전과 변형을 자유롭게 제어함으로써, 기존 유닛 셀 방식의 한계를 극복하고 더 높은 성능 (고 NA, 넓은 대역폭 등) 을 달성할 수 있는 가능성을 열었습니다.
• 이론적 엄밀성: 기하학적 변환과 목적 함수를 분리하여 기울기를 계산하는 체계적인 수학적 프레임워크를 제공함으로써, 향후 다양한 메타물질 설계 문제에 적용 가능한 표준적인 방법론을 제시했습니다.

결론적으로, 이 논문은 어드저인트 변수법과 기하학적 매개변수화를 결합하여 대규모 메타표면의 역설계 문제를 효율적이고 정확하게 해결할 수 있는 강력한 도구를 제공했다는 점에서 의의가 큽니다.