[논문 리뷰] Solving Three Dimensional Maxwell Eigenvalue Problem with Fourteen Bravais Lattices

Tsung‐Ming Huang, Tiexiang Li외 4명|arXiv (Cornell University)|2018. 06. 28.

Photonic Crystals and Applications참고 문헌 32인용 수 7

한 줄 요약

이 논문은 3차원 맥스웰 고유값 문제를 14개의 브라베 격자 전역에서 통합적이고 영공간이 없는 방법으로 해결하기 위해, 직교좌표계에서 블로흐 경계조건을 재정의하고 이산 컬 연산자의 교환성을 활용한다. 주요 기여는 유니타리 동치를 통한 3×3 복소 반대칭 행렬 형태의 표준형으로 이중-컬 연산자의 체계적인 고유분해를 수행함으로써, 3차원 등방성 광학 결정의 밴드 구조 계산을 효율적이고 정확하게 수행할 수 있도록 하는 것이다.

ABSTRACT

Calculation of band structure of three dimensional photonic crystals amounts to solving large-scale Maxwell eigenvalue problems, which are notoriously challenging due to high multiplicity of zero eigenvalue. In this paper, we try to address this problem in such a broad context that band structure of three dimensional isotropic photonic crystals with all 14 Bravais lattices can be efficiently computed in a unified framework. We uncover the delicate machinery behind several key results of our work and on the basis of this new understanding we drastically simplify the derivations, proofs and arguments in our framework. In this work particular effort is made on reformulating the Bloch boundary condition for all 14 Bravais lattices in the redefined orthogonal coordinate system, and establishing eigen-decomposition of discrete partial derivative operators by systematic use of commutativity among them, which has been overlooked previously, and reducing eigen-decomposition of double-curl operator to the canonical form of a 3x3 complex skew-symmetric matrix under unitary congruence. With the validity of the novel nullspace free method in the broad context, we perform some calculations on one benchmark system to demonstrate the accuracy and efficiency of our algorithm.

연구 동기 및 목표

• 14개의 브라베 격자 전역에서 3차원 맥스웰 고유값 문제를 해결하기 위한 통합적 계산 프레임워크를 개발하는 것.
• 발산 자유 제약 조건으로 인해 발생하는 일반화된 고유값 문제에서 큰 영공간 문제를 해결하는 것.
• 모든 격자를 위해 직교좌표계에서 블로흐 경계조건을 재정의하여 수치 이산화와의 일致성과 호환성을 확보하는 것.
• 서로 교환 가능한 이산 편미분 및 컬 연산자의 상호작용을 활용하여 이산 편미분 및 컬 연산자의 체계적인 고유분해를 수립하는 것.
• 유니타리 동치를 통해 이중-컬 연산자를 표준형으로 감소시켜 일반화된 고유값 문제의 효율적이고 정확한 해법을 가능하게 하는 것.

제안 방법

• 기존의 비직각 격자 벡터에서 유도된 새로운 정의된 직교좌표계에서 14개의 브라베 격자 전역에 대해 블로흐 경계조건을 재정의한다.
• 요이 유사 유한차분 기법을 사용하여 ∇×E 및 ∇×H의 컬 연산자를 행렬-벡터 곱으로 이산화함으로써 재정의된 경계조건과의 호환성을 확보한다.
• 이산 편미분 연산자 C1, C2, C3가 상호 교환 가능하다는 것을 입증하여 동시에 대각화할 수 있음을 보장한다.
• 이중-컬 연산자 A = C*C 가 유니타리 동치를 통해 3×3 복소 반대칭 행렬의 표준형으로 변환될 수 있음을 증명한다.
• C1, C2, C3의 고유분해를 활용하여 컬 연산자 C의 영공간과 치역공간의 정규직교 기저를 해석적으로 식별한다.
• 치역공간에 투영하여 영공간이 없는 일반화된 고유값 문제를 구성함으로써, 비물리적 모드가 제거된 해법을 확보한다.

실험 결과

연구 질문

• RQ1 14개의 브라베 격자 전역에서 3차원 맥스웰 고유값 문제를 해결하기 위한 단일 통합 프레임워크를 개발할 수 있는가?
• RQ2 특히 삼단계 격자에 대해, 임의의 브라베 격자에 대해 직교좌표계에서 일관성 있게 블로흐 경계조건을 정의할 수 있는가?
• RQ3 이산 편미분 연산자 간의 교환성이 이중-컬 연산자의 고유분해를 단순화하는 데 어떤 역할을 하는가?
• RQ4 유니타리 동치 하에서 이중-컬 연산자를 표준형으로 감소시킬 수 있으며, 이는 일반화된 고유값 문제의 효율적 해법을 어떻게 가능하게 하는가?
• RQ5 모든 격자에 걸쳐 컬 연산자의 영공간을 체계적으로 해석적으로 특성화하고 제거할 수 있는가?

주요 결과

• 저자들은 14개의 브라베 격자 전역에 대해 직교좌표계에서 블로흐 경계조건을 성공적으로 재정의하여 일관된 이산화를 가능하게 하였다.
• 이산 편미분 연산자 C1, C2, C3 가 상호 교환 가능하다는 것이 입증되었으며, 이는 동시에 대각화 및 해석적 고유분해를 가능하게 한다.
• 유니타리 동치를 통해 이중-컬 연산자 A = C*C 가 표준형 3×3 복소 반대칭 행렬로 감소되었으며, 고유값 문제를 단순화하였다.
• 컬 연산자의 영공간이 해석적으로 식별되고 제거되어, 비물리적 모드를 피하는 영공간이 없는 해법이 가능해졌다.
• 기준 시스템에서의 검증을 통해 3차원 광학 결정의 밴드 구조 계산에서 높은 정확도와 효율성을 입증하였다.
• 프레임워크는 일반적이며 14개의 브라베 격자 전역에 적용 가능하며, 삼단계 격자가 일반적인 경우로 기능한다.