무작위 매질에서 횡광 플럭스의 동역학

초록

본 논문은 상관된 무작위 굴절률을 가진 매질을 통과하는 단일 편광 단색광의 전파 과정을, 파라액시얼‑비파라액시얼 전이와 전기장 구배(‘동역학 에너지’)의 변화를 통해 분석한다. 전파 거리와 상관 길이에 따라 광소용돌이(벡터스) 수가 세제곱근 법칙으로 증가하고, 최대 섬광 지점에서 핵생성률에 ‘키크’가 나타난다. 비감쇠 파동이 필터링되면서 구동된 정상 상태에 도달하고, 비압축성 동역학 에너지 스펙트럼이 파라액시얼에서 완전 무작위 장으로 전이하는 양상을 제시한다.

상세 분석

이 연구는 맥스웰 방정식의 전기장 벡터 파동 방정식을 dyadic Green’s function을 이용해 직접 풀어, 파라액시얼 근사(∂zU≪k0∂zU)에서 비파라액시얼 영역으로의 전이를 정량적으로 기술한다. 무작위 굴절률 변동 η(r)은 Gaussian 상관 함수 C(r)로 모델링되며, 상관 길이 σ와 잡음 강도 κ가 주요 파라미터로 작용한다. 초기 단계에서는 서서히 변하는 envelope U에 대해 파라액시얼 파동 방정식(2ik0∂zU=−∇⊥²U−4πληU)을 도출하고, 전기장 구배의 제곱 적분인 ‘동역학 에너지’ Ekin(z)=∫|∇⊥Ex|² d²r⊥를 정의한다. Fokker‑Planck‑type 접근을 통해 ⟨Ekin⟩의 z‑미분을 평균화하면, ⟨Ekin⟩∝κ²/σ⁴·z라는 선형 성장 법칙을 얻는다. 그러나 수치 시뮬레이션에서는 σ가 작거나 κ가 클 때 비파라액시얼 효과가 급격히 나타나, 선형 예측과 큰 편차(Δ)가 발생한다. 이는 전기장의 y‑편광 성분이 증가하면서 벡터성(vectorial) 효과가 강화되기 때문이다.

벡터스(광소용돌이) 생성률을 정량화하기 위해 2‑D 평면에서 위상 위상점(전위의 위상 고리) 수 ρv(z)를 측정한다. 결과는 세 단계로 구분된다. (1) 파라액시얼 구간에서는 구배 에너지 축적에 의해 소수의 소용돌이가 서서히 나타난다. (2) 중간 구간에서는 ρv∝(z−zcr)¹⁄³ 형태의 급격한 증가가 관찰되며, 이때 zcr은 σ에 따라 달라지는 ‘임계 거리’이다. 특히 σ가 작을수록 이 급증 구간이 길어지고, σ가 크면 최대 섬광 지수 Sx(z)와 일치하는 지점에서 소용돌이 생성이 일시적으로 억제되는 ‘키크’가 나타난다. (3) 장거리 전파에서는 evanescent 파가 필터링되어 파수 스펙트럼이 제한되고, 소용돌이 생성과 재결합이 균형을 이루는 구동 정상 상태에 도달한다. 이 정상 상태의 ρv∞는 랜덤 필드의 스펙트럼 차단(k‑cutoff) 조건에 따라 이론적 상한과 하한 사이에 위치한다.

마지막으로, 전기장 벡터를 Helmholtz 분해하여 비압축성(회전) 성분만을 추출하고, 그 스펙트럼 Ekin(k⊥)를 계산한다. 파라액시얼 영역에서는 k⊥⁻³ 스케일(2‑D Kolmogorov‑like)과 유사한 분포를 보이지만, 비파라액시얼 영역으로 갈수록 고주파(큰 k⊥)가 억제되고, 스펙트럼이 급격히 감소한다. 이는 evanescent 파가 소멸하면서 에너지가 저주파 영역에 집중되는 현상과 일치한다. 전체적으로 논문은 파라액시얼‑비파라액시얼 전이, 동역학 에너지 성장, 소용돌이 통계, 그리고 스펙트럼 전이를 하나의 통합 프레임워크로 연결함으로써, 무작위 매질을 통한 광전파의 복합적 비선형 동역학을 정량적으로 이해하는 데 중요한 기여를 한다.