2차원 포토닉 크리스탈 밴드갭 최적화를 위한 반정밀계획법과 부분공간 기법

초록

본 논문은 2차원 정사각형 격자 포토닉 크리스탈의 밴드갭을 최적화하기 위해, 유한요소법으로 이산화한 무한 차원 고유값 문제를 부분공간으로 축소하고, 이를 반정밀계획법(SDP) 형태의 작은 볼록 최적화 문제로 변환한다. TM·TE 두 편광과 여러 주파수 밴드에 대해 수치 실험을 수행했으며, 최적 설계는 기존 직관과는 다른 복잡한 구조를 보여준다.

상세 분석

이 연구는 포토닉 크리스탈의 밴드갭을 “최대 최소” 형태의 비선형 목표함수로 정의한다. 물리적으로는 유전율 분포 ε(r)와 파동벡터 k에 따라 정의되는 무한 차원의 Hermitian 연산자 A(ε,k)의 고유값 λₙ(ε,k)가 필요하며, 밴드갭은 λₙ₊₁의 최솟값과 λₙ의 최댓값 사이의 차이로 표현된다. 이러한 목표는 비볼록이며, 고유값이 매끄럽게 변하지 않아 미분 가능성이 떨어진다.

저자들은 먼저 2차원 정사각형 격자에 대해 주기적 경계조건을 적용하고, 유한요소법(FEM)으로 공간을 이산화한다. 이 과정에서 파동벡터 k를 Brillouin zone의 고정된 샘플점들에 대해 각각 작은 차원의 고유값 문제로 변환한다. 결과적으로 전체 최적화 문제는 “다중 파라미터·다중 고유값” 형태의 대규모 비선형 프로그램이 된다.

핵심 아이디어는 관심 밴드의 하한과 상한을 결정하는 고유벡터 집합, 즉 “하위 부분공간”과 “상위 부분공간”만을 추출해 문제 차원을 크게 축소하는 것이다. 이 두 부분공간을 각각 Uₗ, Uᵤ라 하면, λₙ ≤ xᵀA(ε,k)x ∀x∈Uₗ와 λₙ₊₁ ≥ yᵀA(ε,k)y ∀y∈Uᵤ라는 반정밀제약(semi‑definite constraints)을 얻을 수 있다. 여기서 ε는 설계 변수이며, 각 물질 영역은 ε∈