리우빌 정리

1. 개요[편집]

리우빌 정리 Liouville's Theorem
발표	Joseph Liouville (1838)
무대	위상공간 (phase space)
내용	해밀토니안 흐름은 위상공간 부피를 보존
모태	해밀토니안 역학
주요 귀결	통계역학 앙상블, 심플렉틱 적분기

리우빌 정리(Liouville’s theorem)는 해밀토니안 방정식을 따라 흐르는 계에서 위상공간의 부피가 시간에 대해 정확히 보존된다는 정리다. 초기 조건들의 집합을 위상공간 속의 한 덩어리로 보면, 그 덩어리는 시간이 지나며 늘어나고 휘어지고 실처럼 가늘어지지만 부피만은 단 한 톨도 변하지 않는다. 위상공간의 흐름은 비압축성 유체의 흐름과 같다는 이야기다.

정리 자체는 한 줄이지만 파급은 물리학 전체를 관통한다. 통계역학의 앙상블 개념이 정당화되는 근거이고, 분자동역학에서 심플렉틱 적분기를 고집하는 이유이며, 가속기 물리의 빔 에미턴스 하한이자, 카오스가 “정보를 잃지 않으면서도 예측 불가능해지는” 방식에 대한 설명이기도 하다. 해밀토니안 역학 문서가 정준방정식과 심플렉틱 구조 전반을 다루므로, 이 문서는 위상공간 부피 보존이라는 한 가지 사실이 시뮬레이션 실무에서 무엇을 강제하는가에 집중한다.

2. 위상공간 흐름의 비압축성[편집]

자유도 $N$인 계의 상태는 $2N$차원 위상공간의 한 점 $(q_1,\dots,q_N,p_1,\dots,p_N)$이다. 이 점의 이동 속도, 즉 위상공간 속도장은 해밀턴 정준방정식이 그대로 준다.

$$\dot{q}_i = \frac{\partial H}{\partial p_i}, \qquad \dot{p}_i = -\frac{\partial H}{\partial q_i}$$

유체역학에서 비압축성의 판정 기준은 속도장의 발산이 0이냐였다. 그대로 해보자.

$$\nabla \cdot \mathbf{v} = \sum_i \left( \frac{\partial \dot{q}_i}{\partial q_i} + \frac{\partial \dot{p}_i}{\partial p_i} \right) = \sum_i \left( \frac{\partial^2 H}{\partial q_i \partial p_i} - \frac{\partial^2 H}{\partial p_i \partial q_i} \right) = 0$$

편미분의 순서 교환 하나로 두 항이 정확히 상쇄된다. 증명이 이렇게 허무하게 끝나는데, 이게 통계역학 전체를 떠받친다.¹ 결론은 위상공간의 흐름이 비압축성 유체라는 것이고, 따라서 흐름을 따라 움직이는 임의의 영역의 부피는 불변이다.

같은 사실을 분포함수 $f(q,p,t)$의 언어로 쓰면 리우빌 방정식이 된다.

$$\frac{\partial f}{\partial t} + \sum_i \left( \frac{\partial f}{\partial q_i}\dot{q}_i + \frac{\partial f}{\partial p_i}\dot{p}_i \right) = \frac{\partial f}{\partial t} + \{f, H\} = 0$$

여기 등장하는 $\{f,H\}$가 포아송 괄호다. 방정식의 뜻은 “흐름을 따라가며 보면 분포함수 값이 변하지 않는다”, 즉 $df/dt = 0$이다. 플라즈마 쪽에서 충돌 없는 블라소프 방정식이 바로 이 리우빌 방정식의 1입자 버전이고, 그래서 블라소프 솔버 품질의 척도가 위상공간 부피와 카시미르 불변량을 얼마나 잘 지키느냐가 된다.

3. 통계역학 앙상블의 근거[편집]

리우빌 정리가 없으면 통계역학은 시작조차 못 한다. 우리는 $10^{23}$개 입자의 초기 조건을 알 수 없으므로, 개별 궤적 대신 위상공간 위의 확률분포를 다루고 평균을 낸다. 이때 두 가지가 보장돼야 한다.

첫째, 평형 분포가 실제로 정지 상태여야 한다. 리우빌 방정식에서 $\partial f/\partial t = 0$이려면 $\{f, H\} = 0$, 즉 $f$가 $H$의 함수이기만 하면 된다. 미시정준 앙상블 $f \propto \delta(H - E)$, 정준 앙상블 $f \propto e^{-\beta H}$가 모두 이 조건을 만족한다. 우리가 쓰는 앙상블들이 시간에 대해 무너지지 않는 이유다.

둘째, 등에너지면 위에서 “균등하게” 센다는 말이 잘 정의돼야 한다. 부피가 보존되니 위상공간 부피 요소 $d^{3N}q\,d^{3N}p$ 자체가 자연스러운 측도가 되고, 등확률 가정(equal a priori probability)이 임의의 선택이 아니라 동역학과 양립하는 선택이 된다. 여기에 에르고드 가설을 얹으면 시간 평균과 앙상블 평균이 같아진다.²

4. 시뮬레이션이 받는 강제[편집]

여기서부터가 실무 이야기다. 해석적으로 부피가 보존된다는 것은 좋은데, 수치 적분기는 그걸 자동으로 지켜주지 않는다.

• 심플렉틱 적분기. 베를레 적분·leapfrog·Yoshida 계열은 이산 시간 사상 자체가 심플렉틱이라 위상공간 부피를 기계 정밀도 수준까지 정확히 보존한다. 대가로 에너지가 정확히 보존되지는 않지만, 참 해밀토니안에 아주 가까운 그림자 해밀토니안(shadow Hamiltonian)을 정확히 보존하기 때문에 에너지 오차가 누적되지 않고 작은 폭으로 진동만 한다.
• 비심플렉틱 적분기. 고전적 4차 룽게-쿠타법은 국소 정확도는 훨씬 높지만 위상공간 부피를 조금씩 수축시킨다. 짧은 궤적에서는 이게 더 정확해 보이지만, 나노초·마이크로초급 장시간 분자동역학을 돌리면 에너지가 단조롭게 흘러내려 계가 조용히 얼어붙는다.³

그래서 MD 커뮤니티는 “정확도 4차” 대신 “부피 보존”을 택했다. 대부분의 MD 코드가 아직도 정확도 2차짜리 벨로시티 베를레를 기본값으로 쓰는 게 게으름이 아니라 리우빌 정리에 대한 존중인 셈이다. 이 문서 상단의 이중진자 데모가 RK4를 쓰고도 멀쩡해 보이는 건 시간 스케일이 짧기 때문이고, 같은 적분기로 며칠을 돌리면 이야기가 달라진다.

5. 서모스탯 — 일부러 부피를 안 지키기[편집]

재미있는 반전은 서모스탯에 있다. 온도를 제어하려면 계와 열욕이 에너지를 주고받아야 하고, 그러면 계 자체는 더 이상 해밀토니안이 아니다. 위상공간 부피는 당연히 보존되지 않는다. 즉 NVT 앙상블을 만들려면 리우빌 정리를 의도적으로 깨야 한다.

깨는 방식의 우아함이 서모스탯의 품질을 가른다.

• Nosé–Hoover는 열욕 자유도 $\zeta$를 상태 변수로 승격시켜 확장 위상공간을 만든다. 확장된 공간에서는 다시 해밀토니안 구조가 회복되고, 일반화된 리우빌 방정식이 성립하며, 물리 부분공간으로 투영했을 때 정확히 정준 분포가 나온다는 것을 증명할 수 있다. 그래서 “정통”으로 대접받는다.
• Berendsen은 속도를 그냥 곱해서 온도를 맞춘다. 빠르고 안정적이지만 대응하는 확장 해밀토니안이 없어 어떤 앙상블도 정확히 재현하지 못한다. 열이 계의 특정 모드로 몰리는 flying ice cube 현상이 그 부작용이다.
• **랑주뱅 동역학**은 마찰과 랜덤 힘을 넣어 아예 확률과정으로 간다. 결정론적 리우빌 방정식 대신 포커–플랑크 방정식이 무대가 되고, 정상 분포가 볼츠만 분포임을 보이는 방식으로 정당화된다.

한 줄 요약: NVE는 리우빌 정리를 지켜서 옳고, NVT는 리우빌 정리를 정확히 계산된 만큼만 깨서 옳다.

6. 카오스와 양립하는 방식[편집]

부피가 보존되는데 어떻게 예측이 불가능해지는가 — 처음 들으면 모순 같지만, 부피 보존은 형상 보존이 아니다. 위상공간의 덩어리는 어떤 방향으로 지수적으로 늘어나면서 다른 방향으로 정확히 같은 비율만큼 납작해질 수 있다. 이것이 리아푸노프 지수의 합이 0이 되는 이유이고, 초기 구가 순식간에 무한히 가는 국수 가락이 되어 등에너지면 전체에 감기는 이유다.

결과적으로 미시적 정보는 하나도 사라지지 않지만, 유한한 해상도로 관측하는 우리에게는 덩어리가 공간 전체에 균등하게 퍼진 것처럼 보인다. 되돌릴 수 없어 보이는 거시적 비가역성이 완전히 가역적인 미시 동역학에서 나오는 이 간극이 통계역학의 오래된 골칫거리이며, 로슈미트 역설과 조사 정리(fluctuation theorem)가 놀고 있는 무대다.⁴ 상단의 이중진자는 자유도가 2뿐이라 이 과정을 눈으로 따라갈 수 있는 드문 사례다.

7. 관련 문서[편집]

• 해밀토니안 역학 · 라그랑주 역학
• 분자동역학 · 서모스탯
• 베를레 적분 · 룽게-쿠타법
• 랑주뱅 동역학 · 블라소프 방정식
• 몬테카를로 방법 · 병렬 템퍼링
• 포아송 괄호
• 계산물리 · 통계

8. Footnotes[편집]

1. 증명 길이 대비 파급 효과 비율로 따지면 물리학에서 손꼽히는 가성비다. 슈바르츠 정리(편미분 교환) 한 줄로 앙상블 이론 전체가 서 있다. ↩

2. 다만 에르고드 가설은 리우빌 정리와 달리 정리가 아니라 가설이며, 실제로 성립하지 않는 계도 많다(유리 상태, 준결정, 통합가능계 등). MD에서 샘플링이 안 돼 고생하는 대부분의 상황이 사실 여기에 걸린 것이다. 병렬 템퍼링 같은 향상된 샘플링 기법이 존재하는 이유. ↩

3. 처음 겪으면 “적분기 차수를 올렸는데 결과가 더 나빠졌다”는 납득 불가능한 현상으로 보인다. 국소 오차와 장기 구조 보존은 완전히 다른 축의 품질 지표라는 것이, MD 입문자가 통과해야 하는 첫 번째 관문이다. ↩

4. 원리적으로는 모든 입자의 속도를 정확히 뒤집으면 커피에 풀린 우유가 도로 뭉친다. 문제는 “정확히”의 요구 정밀도가 우주의 나이 동안 측정해도 부족하다는 것. 리우빌 정리는 그 되돌림이 가능하다고 말하지, 우리가 할 수 있다고는 말하지 않았다. ↩