주기 구조 전파 파수 계산을 위한 완벽 매칭 층(PML) 방법

초록

본 논문은 주기 구조에 입사하는 전자기파의 전파 파수를 정확히 구하기 위해, 무한 영역을 PML로 제한하고 비선형 고유값 문제를 2차 고유값 문제로 변환한 뒤, 유한 요소법으로 수치 해석하는 절차를 제시한다. 이론적 수렴 증명과 다양한 수치 실험을 통해 제안 방법의 정확성과 효율성을 검증한다.

상세 분석

본 연구는 주기 구조에서 발생하는 가이드 모드(전파 파수) 문제를 ‘비선형 고유값 문제 → PML 기반 2차 고유값 문제 → 일반 고유값 문제 → FEM 이산화’라는 네 단계 흐름으로 체계화하였다. 먼저, 파동 방정식 ∆u + k²γ(x)u = 0와 레일리 방사 조건을 만족하는 α‑quasi‑periodic 해를 정의하고, 이러한 해가 존재할 경우 α를 전파 파수라 명명한다. 기존 문헌에서는 레일리 전파 조건만으로는 고유값의 유일성을 보장하지 못함을 지적했으며, 이는 전파 파수가 존재할 때 가이드 모드가 비정상적으로 발생하기 때문이다.

이를 해결하기 위해 저자들은 무한 도메인을 PML(Perfectly Matched Layer)로 트렁케이션한다. PML은 x₂ 방향에 복소 좌표 스트레칭 s(x₂)=s₁(x₂)+i s₂(x₂) (s₁≥1, s₂≥0)를 적용해 흡수 영역을 형성하고, 내부 영역(Ω_H)과 외부 PML 영역(Ω_PML±)을 연속적으로 연결한다. 이때 정의된 PML 연산자 L은 기존 라플라시안에 s(x₂) 가중치를 부여함으로써 원래 문제의 비선형성을 완전히 제거하고, 대신 경계에서 PML‑DtN(Dirichlet‑to‑Neumann) 연산자를 도입한다.

핵심 이론적 기여는 두 가지이다. 첫째, PML‑DtN 연산자가 원래 DtN 연산자에 대해 연산자 노름에서 지수적으로 수렴함을 증명하였다. 이는 PML 두께와 복소 스케일링 파라미터가 충분히 클 경우, 트렁케이션 오차가 exp(−c·δ) 형태로 급격히 감소한다는 강력한 결과다. 둘째, Gohberg‑Sigal 이론을 활용해 파라미터 α가 복소 평면에서 변할 때 고유값(전파 파수)의 연속성과 다중도(멱도) 보존성을 보였다. 특히, 전파 파수가 유한 집합 P(k)⊂