Accuracy assessment of scalar wave propagation methods for diffractive optics design: from thin elements to thick binary grating

기술적 요약: 회절 광학 설계용 스칼라 파동 전파 방법의 정확도 평가

문제 제기
회절 광학 소자 (DOE) 의 설계는 점점 더 기울기 기반 역설계에 의존하고 있으며, 여기서 전파 모델의 정확도는 결과 설계의 성능을 근본적으로 제한합니다. 푸리에 모드 방법 (FMM) 및 유한 차분 시간 영역 (FDTD) 과 같은 엄밀한 방법은 높은 정확도를 제공하지만, 반복적 최적화 루프에서는 계산 비용이 너무 많이 들어 실용적이지 않습니다. 반면, 얇은 요소 근사 (TEA) 는 계산적으로 매우 간단하지만 내부 전파 효과와 인터페이스 결합을 무시하기 때문에 두꺼운 구조나 큰 회절 각도에서는 실패합니다. 중간 스칼라 방법인 빔 전파 방법 (BPM) 과 파동 전파 방법 (WPM) 은 절충안을 제공하지만, 격자 매개변수에 따른 정확도 영역에 대한 체계적이고 정량적인 평가는 문헌에서 부족했습니다.

방법론
저자들은 엄밀한 FMM(FMMAX 를 통해 구현됨) 을 기준 표준으로 사용하여 체계적인 정확도 평가를 수행했습니다. 본 연구는 다음과 같은 실험 설계를 활용했습니다:

• 테스트 구조: $7.5 \, \mu\text{m}$ 주기 셀 위에 무작위 이진 격자 앙상블이 생성되었습니다. 격자는 공간 주파수 차단 ($f_c$) 과 물리적 두께 ($h$) 로 매개변수화되었습니다.
  • $f_c$는 $0$에서$1.33 , \mu\text{m}^{-1}$까지 샘플링되었으며, 이는 최대 회절 각도 ($\theta_{\text{max}}$) 가 약$45^\circ$까지 해당합니다.
  • 두께 ($h$) 는 $1\lambda$에서 $11\lambda$까지 샘플링되었습니다 (여기서 $\lambda = 532 \, \text{nm}$).
  • 총 1,210 개의 격자 (11 $\times$ 11 개의 매개변수 쌍 $\times$ 10 개의 무작위 실현) 가 평가되었습니다.
• 전파 모델: GPU 가속화 및 미분 가능 프레임워크인 FluxOptics.jl 에 세 가지 스칼라 방법이 구현되었습니다:
  • TEA: 무한히 얇은 위상 마스크로 모델링되었으며, 전파 단계는 유효 매질 내에서 근사화되었습니다.
  • BPM: 격자 부피를 $\Delta z = 25 \, \text{nm}$ 단계로 슬라이스하여 근축 근사를 사용했습니다.
  • WPM: 두께에 관계없이 고정된 30 개의 슬라이스를 사용하며 굴절률 불연속을 위해 매끄러운 전이 층을 활용하는 정확한 비근축 전파자 ($k_z = \sqrt{(nk_0)^2 - k_\perp^2}$) 를 적용했습니다.
• 지표: 정확도는 10 개의 실현에 대해 평균낸 스칼라 방법의 투과 필드와 FMM 기준 필드 간의 정규화된 필드 중첩 ($\eta$) 을 사용하여 정량화되었습니다. 스칼라 근사와 일관성을 유지하기 위해 교차 편광 성분은 제외되었습니다.

주요 결과
본 연구는 $(f_c, h)$ 매개변수 공간에서 정확도 맵을 생성하여 각 방법의 명확한 유효 영역을 밝혔습니다:

• TEA: 정확도는 두께가 증가함에 따라 급격히 저하됩니다. 중첩은 중간 각도에서 $h \approx 3\lambda$에서 90% 미만으로 떨어지며, $h > 6\lambda$에서는 50% 미만으로 하락하여 매우 얇은 구조로만 제한됨을 확인했습니다.
• BPM: TEA 보다 훨씬 잘 수행되지만 축적된 근축 위상 오차로 고통받습니다. 유효 영역을 중간 두께까지 확장하지만, 테스트 범위 내 큰 두께에서는 정확도가 0.4~0.5 의 중첩 수준에서 안정화됩니다.
• WPM: 전체 매개변수 공간에서 가장 높은 정확도를 유지합니다. 작은 각도 ($\theta_{\text{max}} \lesssim 12^\circ$) 의 경우, WPM 은 $11\lambda$까지의 두께에서도 $\eta > 0.90$을 유지합니다. WPM 의 BPM 대비 우위는 중간 각도 ($\theta_{\text{max}} \lesssim 30^\circ$) 에서 두꺼운 구조 ($h > 3\lambda$) 일 때 가장 두드러지는데, 이는 정확한 분산 관계가 BPM 을 제한하는 점진적인 위상 오차를 방지하기 때문입니다.

의의 및 주장
본 논문은 날카로운 굴절률 불연속으로 인해 보수적인 벤치마크인 이진 회절 격자에 특화된 이러한 스칼라 방법들에 대한 최초의 체계적인 정량적 평가를 제공한다고 주장합니다. 주요 기여는 회절 광학 역설계 파이프라인에서 전파 모델을 선택하기 위한 실용적인 지침으로 기능하는 정확도 맵을 생성한 것입니다:

1. TEA는 얇은 구조 ($h \lesssim \lambda$) 에만 적합합니다.
2. BPM은 큰 각도의 얇은 구조나 낮은 각도의 중간 두께 구조에 대한 경쟁력 있는 대안입니다.
3. WPM은 중간 각도에서 두꺼운 회절 구조의 역설계를 위한 방법으로 선정되었으며, 엄밀한 방법의 prohibitive 한 비용 없이 BPM 대비 상당한 정확도 이점을 제공합니다.

저자들은 이진 격자가 굴절률 불연속에서 스펙트럼 성분을 최대화하는 "보수적인 벤치마크"를 대표하므로, 보고된 정확도 맵은 매끄러운 연속 요철 구조에서 이러한 방법들의 성능에 대한 하한선일 가능성이 있다고 지적합니다. 본 연구는 세 가지 방법 모두 단방향 전파 및 벡터 효과의 무시와 관련된 한계를 공유하며, 분할 단계 비근축 방법으로의 확장은 향후 과제로 남겨두었다고 인정합니다.