16.12 리우빌 정리와 위상 공간 보존
1. 개요
리우빌 정리(Liouville’s theorem)는 해밀턴 역학의 중심적 결과이다. 이 정리는 해밀턴 흐름이 위상 공간의 부피를 보존한다는 것을 진술한다. 이는 통계 역학과 카오스 이론의 기반이다. 본 절에서는 리우빌 정리와 위상 공간 보존을 다룬다.
2. 리우빌 정리의 진술
2.1 부피 보존
해밀턴 흐름은 위상 공간의 부피를 보존한다. 즉, 위상 공간에서 임의의 영역 D가 시간 t 후에 영역 D_t로 변하면
\text{vol}(D_t) = \text{vol}(D)
2.2 미분 형식
미분 형식으로 표현하면, 위상 공간의 부피 형식 \Omega = dq_1 \wedge dp_1 \wedge \cdots \wedge dq_n \wedge dp_n이 해밀턴 흐름에 의해 보존된다.
\mathcal{L}_{\mathbf{X}_H}\Omega = 0
여기서 \mathcal{L}은 리 미분이고 \mathbf{X}_H는 해밀턴 벡터장이다.
2.3 발산이 0
리우빌 정리는 해밀턴 벡터장의 발산이 0이라는 것과 동등하다.
\text{div}\mathbf{X}_H = 0
3. 증명
3.1 해밀턴 벡터장
해밀턴 벡터장은 다음과 같이 정의된다.
\mathbf{X}_H = \left(\frac{\partial H}{\partial p_1}, \ldots, \frac{\partial H}{\partial p_n}, -\frac{\partial H}{\partial q_1}, \ldots, -\frac{\partial H}{\partial q_n}\right)
3.2 발산의 계산
발산은 다음과 같이 계산된다.
\text{div}\mathbf{X}_H = \sum_i\left(\frac{\partial}{\partial q_i}\frac{\partial H}{\partial p_i} + \frac{\partial}{\partial p_i}\left(-\frac{\partial H}{\partial q_i}\right)\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
(편미분의 교환성에 의해)
3.3 결론
따라서 해밀턴 벡터장의 발산이 0이며, 위상 공간의 부피가 보존된다.
4. 리우빌 방정식
4.1 분포 함수
위상 공간에서의 분포 함수 \rho(\mathbf{q}, \mathbf{p}, t)는 시스템의 가능한 상태에 대한 확률 분포이다.
4.2 보존
리우빌 정리는 분포 함수가 다음의 방정식을 만족한다는 것을 의미한다.
\frac{d\rho}{dt} = \frac{\partial \rho}{\partial t} + \{ \rho, H\} = 0
여기서 \{\rho, H\}는 푸아송 괄호이다. 이를 리우빌 방정식이라 한다.
4.3 의의
분포 함수가 해밀턴 흐름을 따라 일정하게 유지된다. 이는 통계 역학의 기반이다.
5. 통계 역학에서의 응용
5.1 통계 앙상블
통계 역학에서 시스템의 상태는 분포 함수로 표현된다. 통계 앙상블(canonical ensemble, microcanonical ensemble 등)이 위상 공간의 분포로 정의된다.
5.2 등확률 가설
같은 에너지의 모든 상태가 동일한 확률을 가진다는 등확률 가설은 리우빌 정리에 기반한다.
5.3 평형 분포
평형 상태에서 분포 함수는 시간에 의존하지 않으며, \{\rho, H\} = 0을 만족한다. 이는 \rho가 보존량의 함수임을 의미한다.
6. 카오스 이론에서의 응용
6.1 위상 공간의 변형
해밀턴 흐름은 위상 공간의 부피를 보존하지만 영역의 모양은 변할 수 있다. 카오스 시스템에서는 영역이 매우 복잡한 형태로 변한다.
6.2 혼합
카오스 시스템에서 위상 공간의 영역이 시간이 지나면서 더 복잡해지고, 결국 위상 공간 전체에 걸쳐 균일하게 분포한다. 이를 혼합(mixing)이라 한다.
6.3 에르고드성
에르고드 시스템은 시간 평균과 공간 평균이 같은 시스템이다. 이는 통계 역학의 기반이다.
7. 시뮬렉틱 적분기
7.1 동기
리우빌 정리는 해밀턴 시스템의 핵심 성질이다. 일반 수치 적분 방법은 이 부피 보존을 정확히 보존하지 않는다.
7.2 시뮬렉틱 적분기
시뮬렉틱 적분기는 시뮬렉틱 형식을 정확히 보존하는 수치 적분 방법이다. 따라서 위상 공간의 부피도 보존한다.
7.3 장기 시뮬레이션
시뮬렉틱 적분기는 해밀턴 시스템의 장기 시뮬레이션에서 우수한 정확도를 제공한다. 천체 역학, 분자 동역학 등에서 자주 사용된다.
8. 매니퓰레이터에서의 의미
8.1 보존계의 시뮬레이션
마찰이 없는 매니퓰레이터의 시뮬레이션에서 위상 공간 부피가 보존되어야 한다.
8.2 시뮬렉틱 적분기
매니퓰레이터의 정확한 시뮬레이션에 시뮬렉틱 적분기가 활용된다.
8.3 마찰의 영향
마찰이 있으면 시스템이 비보존이 되며, 위상 공간 부피가 일반적으로 감소한다(영역이 끌개로 수렴).
9. 본 절의 의의
본 절은 리우빌 정리와 위상 공간 보존을 다루었다. 리우빌 정리는 해밀턴 역학의 중심적 결과이며, 통계 역학, 카오스 이론, 정확한 수치 적분의 기반이다.
10. 참고 문헌
- Arnold, V. I. (1989). Mathematical Methods of Classical Mechanics (2nd ed.). Springer.
- Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
- Pathria, R. K., & Beale, P. D. (2011). Statistical Mechanics (3rd ed.). Academic Press.
- Hairer, E., Lubich, C., & Wanner, G. (2006). Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations (2nd ed.). Springer.
version: 1.0