비대칭 광자 결정 웨이브가이드의 전계 분포 해석

최근 광집적 회로 및 나노광학 소자 설계에서 단방향 복사 특성을 갖는 구조가 주목받고 있다. 이러한 시스템은 특정 방향으로만 빛을 방출하는 특성 덕분에 고감도 센서나 일방향 광학 격리기에 응용된다. 핵심 원리는 비대칭적인 유전체 패터닝을 통해 전자기 모드의 공간 대칭성을 깨는 것이다.

예를 들어, 상부에 아파장 구조의 광자 결정 격자를 형성하고 하부는 균일한 기판으로 유지하는 경우, 전계 분포는 위쪽으로 더 강하게 확산된다. 이 현상을 수치적으로 분석하기 위해 1차원 유한차분법(FDM) 기반 파동 방정식 해법을 적용할 수 있다.

import numpy as np
from scipy.linalg import eigh

def compute_modes(grid_points=300, dx=0.015):
    # 공간 그리드 정의
    x = np.linspace(0, grid_points * dx, grid_points)
    
    # 유전율 분포 설정 (비대칭 계층 구조)
    permittivity = np.ones(grid_points)
    permittivity[40:90] = 12.25   # Si층 (n=3.5)
    permittivity[100:150] = 2.25  # 낮은 산란층 (n=1.5)
    
    # 파수 벡터 및 파장 정의
    k0 = 2 * np.pi / 1.55
    
    # 이산화된 헬름홀츠 연산자 구성
    diagonal = -(2 / dx**2) + (k0**2) * permittivity
    off_diag = 1 / dx**2 * np.ones(grid_points - 1)
    
    A = np.diag(diagonal) + np.diag(off_diag, 1) + np.diag(off_diag, -1)
    
    # 고유값 문제 해결
    eigenvalues, eigenvectors = eigh(A)
    
    return eigenvalues, eigenvectors, x

위 코드는 1D 파동 방정식을 이산화하여 전자기 모드의 고유값과 고유벡터를 계산한다. 고유값의 실수부는 모드의 전파 특성과 관련되며, 허수부는 복사 감쇠율을 반영한다. 특히, 비대칭 경계 조건에서는 특정 고조파 모드가 비대칭적으로 외부로 복사되는 것을 관찰할 수 있다.

복사 Q-인자 및 비대칭 손실 메커니즘

복사 구조의 효율성은 일반적으로 Q-인자로 평가되며, 이는 저장된 모드 에너지의 감쇠 시간과 직접적으로 관련된다. 단순한 계산식은 다음과 같다:

def calculate_quality_factor(eig_vals, mode_index=0):
    omega = eig_vals[mode_index]
    omega_r = np.real(omega)
    gamma = abs(np.imag(omega))
    
    if gamma == 0:
        return float('inf')
    return omega_r / (2 * gamma)

그러나 실제 단방향 구조에서는 상하 표면에서의 복사 손실이 불균형하게 발생한다. 이를 반영하기 위해 수정된 Q-인자 모델을 도입할 수 있다:

\[ Q_{\text{eff}} = \frac{Q}{\alpha \cdot \Gamma_{\text{top}} + (1 - \alpha) \cdot \Gamma_{\text{bottom}}} \]

여기서 \(\alpha\)는 상부 복사 비율을 조절하는 가중치이며, 실험적 데이터 피팅을 통해 약 0.65~0.75 사이의 값을 취하는 것이 일반적이다. 이 값은 격자의 주기, 심도, 그리고 재료의 굴절률 차이에 따라 달라진다.

에너지 밴드 구조 및 특이점 분석

주기적 구조의 경우, 평면파 전개법(PWE)을 사용하여 k-공간에서의 밴드 다이어그램을 계산할 수 있다. 특히, 비대칭성이 존재할 때는 비헤르미트 시스템의 특이점(EP, Exceptional Point)이 등장하며, 이는 두 개의 고유모드가 병합하는 지점으로, 방향성 복사 증폭 현상과 연결된다.

% MATLAB 기반 PWE: 푸리에 계수 계산
N = 15;  % 충분한 차수 보장을 위한 절단
G = zeros(2*N+1, 1);
period = 0.6;

for n = -N:N
    integrand = @(x) epsilon_profile(x) .* exp(-1i * 2*pi*n*x/period);
    G(n+N+1) = integral(integrand, 0, period) / period;
end

푸리에 급수의 절단 수준 \(N\)이 너무 작으면 고차 회절 성분이 무시되어 밴드 구조의 정확성이 저하된다. 특히, 금지대 경계 근처에서 비물리적 "아티팩트 피크"가 나타날 수 있으므로, \(N \geq 15\)를 권장한다.

수치 해석 시 고려사항

• 격자 크기: 구조의 특징 치수(예: 격자 폭)의 1/8 이하로 설정해야 정확한 모드 분포를 얻을 수 있다.
• 경계 조건: PML(Perfectly Matched Layer)을 사용하여 인위적인 반사를 최소화해야 한다.
• 메쉬 적응: 중요 영역(예: 유전율 급변부)에서 자동 메쉬 세분화를 적용하면 계산 효율과 정확도 모두 향상된다.

이러한 요소들은 이론적 예측과 실험 결과 간의 갭을 줄이는 데 핵심적인 역할을 하며, 고성능 단방향 광학 소자의 설계 프로세스에서 필수적인 검토 항목이다.