전역균형으로 여는 새로운 몬테카를로 시대: 이벤트 체인 알고리즘의 혁신

초록

Bernard·Krauth·Wilson(2009)이 제안한 이벤트 체인 몬테카를로(ECMC)는 상세균형을 포기하고 전역균형만을 만족함으로써, 하드 구체 시스템에서 거부 없이 결정론적인 움직임을 구현한다. 이 논평은 ECMC의 원리를 전역‑균형·리프팅 마코프 체인 이론과 연결하고, 연속적인 포텐셜과 현대 샘플링 기법으로 일반화하는 과정을 상세히 설명한다.

상세 분석

ECMC의 핵심은 “양방향 전이 확률을 동일하게 맞추는 상세균형(detailed balance)” 대신, “모든 상태에 들어오는 흐름과 나가는 흐름이 전체적으로 일치하는 전역균형(global balance)”을 이용한다는 점이다. 이를 위해 알고리즘은 입자를 고정된 방향으로 연속 이동시키고, 충돌이 발생하면 이동 권한을 다음 입자에게 전달한다(이벤트 체인). 이동 거리의 총합이 사전에 정해진 체인 길이 l에 도달하면 체인이 종료된다. 이 과정은 역방향으로 동일한 체인을 실행하면 원래 상태로 완벽히 복원되므로, 마코프 체인의 전이 행렬이 비가역적이지만 전역균형을 만족한다.

논문은 이를 “리프팅(lifted) 마코프 체인” 개념으로 일반화한다. 물리적 상태 x에 내부 변수 v(±1) — 이동 방향 혹은 속도 — 를 추가함으로써 상태 공간을 두 배로 확장한다. 리프팅된 전이에서는 충돌 시 v가 반전되며, 이는 전통적인 Metropolis 알고리즘에서의 거부(rejection)를 대체한다. 충돌 확률은 현재 위치 x와 포텐셜 기울기 ∇U(x)에 의해 결정되며, “상승 방향(에너지 감소)으로는 충돌이 없고, 하강 방향(에너지 증가)에서만 충돌이 발생한다”는 규칙을 따른다. 이 규칙은 흐름을 가능한 한 일방향으로 만들면서도 전역균형을 보장한다.

수학적으로는 전이 확률 T(x→x′) 를 “대칭 시도 커널 T_trial”과 “min(π(x),π(x′))” 형태로 분해하고, 이를 연속시간 한계(Δt→0)에서 미분 형태의 이벤트 드리븐(ED) 방정식으로 전환한다. 이때 충돌 확률 P_coll(v,x)=max