선형 라그랑지안 적분기의 확률적 베르레트 분석

초록

본 논문은 베르레트 형태의 확률적 적분기를 구성‑좌표만으로 기술하고, 시간 단계에 따른 확산 상수, 기울기 이동, 조화 진동자 샘플링 세 가지 기본 통계량을 폐쇄형식으로 유도한다. 12개의 대표적 알고리즘을 분석한 결과, GJ 계열만이 모든 시간 단계에서 세 지표를 정확히 재현함을 확인하고, 이를 고품질 열역학 시뮬레이션에 권장한다.

상세 분석

논문은 먼저 연속 라그랑지 방정식 m ¨r + α ·r = f + β 에 대해, 확산 D = kBT/α, 기울기 vd = f/α, 그리고 조화 포텐셜 Ep = ½κr² 에 대한 구성‑좌표 평균 ⟨r²⟩ = kBT/κ 이라는 세 가지 기본 관계를 기준으로 삼는다. 이후 베르레트‑형식의 확률적 적분기를
rₙ₊₁ = 2c₁rₙ − c₂rₙ₋₁ + Δt²/m (c₃fₙ + c₄fₙ₋₁) + Δt²/m (c₅βₙ⁻ + c₆βₙ⁺)
형태로 일반화하고, 파라미터 cᵢ(γΔt)와 잡음 상관 ζ(γΔt)만을 변수로 두어 모든 스토캐스틱 베르레트 적분기를 포괄한다. 여기서 γ = α/m이며, 안정성 조건 |c₂|≤1, 2c₁ = 1 + c₂ 를 만족한다.

핵심은 세 개의 정규화된 지표 Γ_diff, Γ_drift, Γ_dist 을 각각 확산, 기울기, 샘플링 분포에 대해 도출한 것이다. 이들은 (A6), (B2), (C5)식에 의해 cᵢ와 ζ에 대한 함수로 명시되며, 정확한 연속‑시간 결과를 재현하려면 세 값이 모두 1이어야 한다. 논문은 이를 바탕으로 12개의 기존 적분기(SS78, EB80, MP‑80~82 등)를 분석하고, 각 파라미터 집합을 대입해 Γ값을 계산한다. 대부분의 알고리즘은 Γ_diff와 Γ_drift는 정확히 1이지만, Γ_dist가 시간 단계에 따라 크게 편차를 보이며, 특히 저감쇠(γ→0)에서 안정성 한계가 급격히 감소한다는 점을 확인한다.

반면 GJ(Grøn‑Jensen) 계열은 cᵢ와 ζ를 특별히 설계해 Γ_diff = Γ_drift = Γ_dist = 1을 모든 γΔt와 Ω₀Δt에 대해 만족한다. 이는 잡음 가중 함수 ψ(s)를 대칭적으로 선택하고, c₅·c₆·ζ = 4(γΔt)² ·