중첩 확산과 레비 비행: 무거운 꼬리 목표와 비정상 확산 현상

본 논문은 구속된 환경에서 확률 과정이 무거운 꼬리(heavy‑tailed) 목표 확률밀도함수(pdf)를 갖는 현상을 다각도로 탐구한다. 저자는 먼저 전통적인 Lévy‑안정분포가 독립·동일분포 점프 과정의 전형적인 결과라는 일반적인 인식을 검토하고, 이를 “역공학”(reverse engineering) 전략을 통해 목표 pdf ρ_*(x)를 사전에 지정하고 해당 pdf에 수렴하는 점프형 과정을 설계하는 방법을 제시한다. 기존의 Lévy‑Langevin 접근법에서는 드라이버가 프랙탈 차수 0<μ<2 인 분수 라플라시안을 포함한 Fokker‑Planck 방정식을 사용한다. 이 경우 구속 퍼텐셜이 존재하면, 장기적으로는 Cauchy와 같은 무거운 꼬리 pdf가 나타날 수 있으며, 순간수(특히 분산)의 존재 여부는 μ와 퍼텐셜 형태에 따라 달라진다. 그 다음 저자는 Lévy‑Schrödinger 세미그룹 접근법을 도입한다. 여기서는 확률 과정이 세미그룹 연산자 exp(−tĤ) 로 기술되며, Ĥ는 프랙탈 라플라시안과 잠재적 V(x) 로 구성된다. 목표 pdf 가 주어지면 V(x)=−λ(|Δ|^{μ/2}√ρ_*(x))/√ρ_*(x) 로 정의할 수 있고, 이때 세미그룹이 존재하면 해당 pdf는 고유 불변 상태가 된다. 중요한 점은 이 방법이 Langevin‑driven 접근법과 달리 비가우시안 형태의 포텐셜을 직접 사용함으로써 “비가스” (non‑Gibbsian) 장애물을 회피한다는 것이다. 핵심 전환점은 동일한 목표 pdf가 연속형 확산 과정에서도 달성될 수 있다는 사실이다. 저자는 로그형 퍼텐셜 V(x)=ln(1+x²) 를 갖는 전통적인 Fokker‑Planck 방정식(두 번째 미분 연산자만 사용) 을 고려한다. Shannon 엔트로피 S(ρ)=−∫ρ lnρ dx 를 최대화하면서 제약 ⟨ln(1+x²)⟩=U_α 를 부과하면, 라그랑주 승수 α에 따라 ρ_α(x)=C_α/(1+x²)^α 형태의 Cauchy 계열이 도출된다. α>½ 인 경우 정규화가 가능하고, α가 정수이면 순간수의 존재 여부가 단계적으로 변한다. 예를 들어 α=1에서는 모든 순간수가 발산하지만, α가 커질수록 2,4,… 차수의 순간수가 순차적으로 수렴한다. 이는 “모멘트 사라짐” 현상이 시간에 따라 전이되는 메커니즘을 제공한다. 시간 의존적인 분산 ⟨x²⟩(t)∼t^γ 에서 γ는 0<γ<2 로, 구속된 환경에서도 서브‑디퓨전(γ<1), 정상 디퓨전(γ=1), 슈퍼‑디퓨전(γ>1) 모두 관찰된다. 그러나 γ는 드라이버 종류(Lévy vs Wiener)와 퍼텐셜 형태에 따라 달라지며, 보편적인 시간 스케일 계층 구조는 존재하지 않는다. 저자는 수치 시뮬레이션을 통해 초기 δ‑함수 혹은 가우시안 분포가 점차 Cauchy 계열의 목표 pdf 로 수렴하는 과정을 보여준다. 특히, Wiener‑드라이버만을 사용했음에도 불구하고 무거운 꼬리 분포가 생성된다는 점은 기존의 “점프형 과정만이 무거운 꼬리를 만든다”는 통념을 뒤흔든다. 또한, 논문은 Tsallis 비가역 엔트로피와 관련된 기존 문헌을 비판한다. 기존 연구들은 무거운 꼬리 분포를 얻기 위해 비표준 제약이나 비가역 엔트로피(q‑엔트로피)를 도입했지만, 저자는 Shannon 엔트로피와 적절한 로그 제약만으로도 동일한 결과를 얻을 수 있음을 증명한다. 이는 무거운 꼬리 분포가 반드시 비가역 통계역학을 필요로 하지 않음을 시사한다. 결론적으로, 저자는 Lévy‑점프와 Wiener‑확산 두 과정이 동일한 무거운 꼬리 목표 pdf 를 공유할 수 있음을 보이며, 실험적으로 관측된 “확산 → Lévy 비행 전이” 현상이 실제로는 두 과정 사이의 구분이 모호한 전이 현상일 가능성을 제시한다. 이러한 통합적 시각은 무거운 꼬리 현상을 모델링할 때 보다 폭넓은 선택지를 제공하고, 기존의 비가스적 열역학 해석을 재고하게 만든다.