고품질 공진 현상에 대한 2차원 FDTD 정확도 검증

초록

본 연구는 원형 유전체 공진기에서 발생하는 고품질(Q) 위스퍼링 갤러리 모드(WGM)의 전자기 해석에 있어 표준 2차원 유한 차분 시간 영역(FDTD) 방법의 정확성을 평가한다. 원형 공진기를 선전류로 구동하고, FDTD 결과를 정확한 미에 해법과 비교하여 공진 주파수 이동 및 Q‑factor 의존성을 정량화하였다. 석영과 실리콘 재질, 다양한 반경, E·H 편광을 모두 고려한 종합적인 오류 분석을 제공한다.

상세 요약

본 논문은 고-Q 위스퍼링 갤러리 모드(WGM)의 특성을 수치적으로 재현할 때, 2차원 유한 차분 시간 영역(FDTD) 방법이 얼마나 신뢰할 수 있는지를 체계적으로 검증한다. 먼저, 원형 유전체 공진기를 선전류 소스로 직접 구동함으로써 실제 실험 환경에 근접한 전자기 흥분 조건을 설정하였다. 이때 사용된 공진기 재료는 전자기 손실이 거의 없는 석영(ε≈4.5)과 고굴절률 실리콘(ε≈11.7)이며, 반경은 파장 대비 5배에서 30배까지 다양하게 설정해 규모에 따른 수치 오차를 조사하였다.

FDTD 시뮬레이션은 표준 Yee 격자를 이용해 2차원 TM(E‑polarization) 및 TE(H‑polarization) 모드 모두에 대해 수행되었으며, 격자 간격 Δx, Δy와 시간 스텝 Δt를 각각 λ/30, λ/30, λ/(20c) 수준으로 세밀하게 조정하였다. 이러한 파라미터 설정에도 불구하고, 고-Q(>10⁴) 모드에서는 공진 주파수가 미에(Mie) 해법에 비해 몇 MHz 수준으로 이동하는 현상이 관찰되었다. 특히, 실리콘 공진기의 경우 굴절률 차이와 경계 조건의 급격한 변화가 격자화 오류를 증폭시켜, Q가 10⁵에 달하는 모드에서는 주파수 오차가 0.1%에 육박하였다.

오차 원인 분석에서는 (1) 격자화에 의한 경계면 근사, (2) 시간 스텝에 따른 수치 분산, (3) 완전 흡수 경계조건(PML)의 비정상적인 반사, (4) 선전류 소스의 공간적 이산화 등을 꼽았다. 특히, 경계면 근사는 원형 구조를 직사각형 격자에 매핑하면서 발생하는 ‘staircase’ 효과가 고-Q 모드의 전자기 에너지 저장량을 과소평가하게 만들며, 이는 공진 주파수의 체계적 하향 이동으로 이어진다. 또한, PML 두께를 충분히 확보하지 않을 경우, 고주파 성분이 반사되어 모드의 품질인자(Q)를 인위적으로 감소시키는 현상이 확인되었다.

논문은 이러한 오류를 정량화하기 위해, 공진 주파수 이동 Δf와 Q‑factor의 비례 관계를 도출하였다. Δf/f₀ ≈ α·(Δx/λ)·(1/Q) 형태의 경험적 식을 제시함으로써, 사용자가 원하는 정확도를 달성하기 위해 필요한 격자 해상도를 사전에 예측할 수 있게 하였다. 최종적으로, 고-Q 모드 해석에 있어 표준 FDTD는 충분히 작은 격자와 적절한 PML 설계가 전제될 경우, 미에 해법과 0.01% 이내의 일치도를 보이며 실용적이라고 결론지었다.