ENZ 매질에서 비선형 광구조 시뮬레이션

초록

본 논문은 유효 유전율이 영에 가까운(ENZ) 등방성 메타물질 내부에서 약한 공간 불균일성, 케르 효과 및 외부 펌핑에 의한 선형 이득을 포함한 광파의 동역학을 기술하는 벡터형 Ginzburg‑Landau 방정식을 제안한다. 실재적인 케르 계수를 가정하고, 주된 색인 연산을 curl‑curl 형태로 두어 split‑step Fourier 방법으로 수치 시뮬레이션을 수행하였다. 결과적으로 구대칭 및 토러스형 벡터 파동 구조가 2차 포텐셜 우물에 포획되는 다양한 비선형 진화를 확인하고, 종축·횡축 파동 간의 비선형 상호작용을 관찰하였다.

상세 분석

이 연구는 ENZ(ε≈0) 매질의 전자기 파동을 기술하기 위해 기존의 스칼라 비선형 파동 방정식이 아닌, 전기장 벡터 자체를 변수로 하는 3차원 Ginzburg‑Landau 방정식을 도입한 점이 가장 큰 혁신이다. 방정식(5)는 시간 1차 미분 형태로 정리되어, 수치적으로 시간 전진이 직접 가능하도록 구성되었다. 여기서 핵심 연산자는 ‘curl‑curl’ 항으로, ENZ 매질에서 전자기 파동의 전파 속도가 거의 무한대가 되면서 전기장의 회전 성분이 지배적임을 반영한다. 또한, 비등방성 분산 파라미터 D∥와 점성 파라미터 ν⊥, ν∥을 도입해 종축(longitudinal)과 횡축(transverse) 모드의 차이를 명시적으로 구분하였다.

비선형 항은 두 개의 케르 계수 α와 β(복소수)로 구성되며, α′가 실재적인 자가위상 변조를, β가 편광 상호작용을 담당한다. 논문에서는 β/α′=η를 실험적으로 0.5 정도로 설정하고, 실수 케르 계수만을 고려해 보수적인 시뮬레이션 환경을 만든다. 선형 이득 G는 외부 펌핑에 의해 도입되며, γ(흡수)와 반대되는 부호로 작용해 파동의 성장·소멸을 조절한다.

수치 해법은 split‑step Fourier 방법을 두 단계로 나눈다. 첫 번째 단계에서는 로컬 비선형 진화(식 7)를 반시간 간격(Δt/2)만큼 적용하고, 두 번째 단계에서는 푸리에 공간에서 선형 분산(식 9)을 전 시간 간격(Δt) 동안 적용한다. 이 과정을 M번 반복해 전체 시간 스텝 M·Δt를 진행한다. 2차 정확도(Δt²)와 에너지 보존 특성을 동시에 만족시키는 알고리즘으로, 고해상도(N=192³) 격자에서도 수일간의 연산을 견딜 수 있다.

시뮬레이션 결과는 크게 세 가지 유형으로 구분된다. 첫 번째는 구대칭 ‘벌레’(central‑symmetric) 구조로, 전기장이 방사형으로 퍼지면서 Gaussian 형태의 포텐셜 우물에 포획된다. 두 번째는 토러스형 ‘도넛’(toroidal) 구조로, 전기장이 azimuthal 방향으로 회전하며 핵심부에 강한 자기 회전성을 보인다. 세 번째는 두 구조의 중첩으로, 종축과 횡축 파동이 비선형적으로 결합해 복잡한 파동 ‘패턴’이 형성된다.

특히, 이득 G⊥가 임계값(≈0.6~0.75) 이상일 때 토러스 구조는 급격한 진폭 폭발을 겪으며, 이후 dipole 모드가 우세해지는 불안정성을 보인다. 반면, 낮은 G⊥에서는 안정적인 비선형 평형 상태가 유지되며, 파동 행동이 포텐셜 우물에 의해 억제된다. 중앙 대칭 구조는 G가 증가함에 따라 진폭이 서서히 증가하다가 급격한 수축을 겪어 ‘콜랩스’ 현상에 근접한다는 점이 관찰되었다. 이는 3차 비선형과 2차 분산을 갖는 보수적인 시스템에서 알려진 불안정성(폭발적 자기집중)과 일치한다.

마지막으로, 종축·횡축 파동이 동시에 존재할 때는 작은 스케일의 ‘파동 결절’(wave‑node)들이 형성되며, 이는 ENZ 매질 내부에서 비선형 편광 상호작용이 새로운 위상 구조를 만들어낼 가능성을 시사한다. 이러한 현상은 현재 실험적으로 검증되지 않았으나, 메타물질 설계 시 편광·위상 제어를 위한 새로운 설계 지표가 될 수 있다.