복소 분산 방정식 근을 완전 탐색하는 새로운 전산 기법

초록

**
본 논문은 복소수 형태의 분산 방정식에서 모든 복소 근을 효율적으로 찾기 위한 새로운 수치 방법을 제시한다. 격자 기반 선형·포물선 근사와 엄격한 판정 기준을 결합해 기존의 그리드 탐색과 뉴턴‑라프슨 방식이 갖는 해 누락·수렴 실패 문제를 극복한다. 방법의 유효성을 검증하기 위해 원통형 금속 나노와이어의 표면 플라스몬 폴라리톤(SPP) 분산식에 적용했으며, 복소‑ω 해와 복소‑k 해 모두에서 브루스터 각 궤도와 다중 모드( n=0~3 )의 복소 분산 곡선을 한 번의 계산으로 도출하였다. 또한 근사 해석식을 도출해 수치 결과와 일치함을 확인하였다.

**

상세 분석

**
이 연구는 복소 전파 방정식 f(x,y)=0 을 두 실변수 x, y (예: Re ω, Im ω 또는 Re k, Im k) 로 표현한 뒤, 전통적인 전역 그리드 탐색이 요구하는 초미세 격자와 허용 오차 δ 에 의존하는 방식을 탈피한다. 저자들은 먼저 격자 A 내에 근이 존재한다고 가정하고, 네 개의 코너값 f₁~f₄ 을 이용해 1차 테일러 전개를 수행한다. 이때 ∂f/∂x, ∂f/∂y 를 코너값 차이로 역산하여 Δa₁, Δb₁ 등을 구하고, “Δa₁≈Δa₂”, “Δb₁≈Δb₂”와 같은 일치 조건을 δ 이내로 만족하면 해당 격자가 근을 포함한다는 판정 기준(식 4b)을 제시한다. 이렇게 하면 격자 A 가 독립적으로 평가되므로 다중 근이 존재하더라도 누락되지 않는다.

하지만 1차 근사가 불충분한 경우(예: ∂f/∂x 또는 ∂f/∂y 가 거의 0인 경우)에는 2차 테일러 전개를 적용한다. 격자 B 의 8개 포인트를 이용해 2차 편미분 fₓₓ, f_yy, f_xy 를 계산하고, 연립식 (7a,b) 를 풀어 Δc, Δd (즉 ΔRe k, ΔIm k) 값을 얻는다. 여기서도 “fₓₓ≈f₂ₓₓ”, “f_yy≈f₂_yy”와 같은 일치 조건을 δ 이내로 검증해 근 존재 여부를 판단한다(식 8a‑c). 선형·2차 두 단계 판정 체계는 복소 함수의 급격한 변곡이나 평탄 구간에서도 안정적으로 근을 포착한다.

알고리즘 흐름은 (1) 전체 탐색 영역을 적당한 Sₓ, S_y 격자로 나누고, (2) 각 격자에 대해 위의 판정식을 적용해 후보 격자를 추출, (3) 후보 격자 내부를 재그리드하거나 최소 |f| 값을 갖는 점을 최종 근으로 선택한다. 이 과정은 병렬화가 용이하고, 격자 간 의존성이 없으므로 대규모 파라미터 스캔에 적합하다.

실제 적용 사례로 원통형 금속 나노와이어의 복소‑ω 분산식을(식 9) 사용하였다. 금속은 드루드 모델 ε₁(ω)=ε_∞