해밀토니안 역학

물리 계산물리 마지막 수정: 2026-07-21 05:03:46

1. 개요[편집]

해밀토니안 역학 Hamiltonian Mechanics
정립	윌리엄 로원 해밀턴 (1833)
핵심 개념	위상공간(phase space), 정준방정식
구조	심플렉틱 구조(symplectic structure)
관련	라그랑주 역학
주요 응용	분자동역학, N-body 시뮬레이션, 심플렉틱 적분기

해밀토니안 역학(Hamiltonian mechanics)은 계의 상태를 일반화 좌표 $q$ 와 일반화 운동량 $p$ 로 이루어진 위상공간(phase space) 위의 한 점으로 나타내고, 해밀토니안이라 불리는 스칼라 함수 $H(q, p, t)$ 하나로부터 시간에 따른 전체 운동을 기술하는 고전역학의 재정식화이다. 라그랑주 역학을 르장드르 변환(Legendre transformation)해 얻어지며 물리적으로는 완전히 동등하지만, 속도 $\dot{q}$ 대신 운동량 $p$ 를 좌표와 대등한 독립변수로 승격시킨다는 점에서 구조가 근본적으로 다르다.

라그랑주에서 해밀토니안으로 넘어가는 유도 과정 자체는 라그랑주 역학 문서에서 이미 다뤘으므로, 이 문서는 그 중복을 피해 위상공간의 기하학적 구조, 심플렉틱 성질, 그리고 분자동역학·N-body 시뮬레이션에서의 실전 활용에 집중한다. 결론부터 말하면 해밀토니안 형식이 시뮬레이션 실무에서 중요한 이유는 딱 하나다 — 에너지를 오래 보존하는 수치적분기를 설계할 수 있는 이론적 근거를 주기 때문이다.

2. 정준방정식[편집]

라그랑지안 $L(q, \dot{q}, t)$ 에서 일반화 운동량 $p_i \equiv \partial L / \partial \dot{q}_i$ 을 정의하고 $H = \sum_i p_i \dot{q}_i - L$ 로 르장드르 변환하면, 하나의 2계 미분방정식(오일러-라그랑주 방정식)이 두 개의 1계 미분방정식 — 해밀턴 정준방정식(canonical equations) — 으로 나뉜다.

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

두 방정식이 $q$ 와 $p$ 에 대해 완전히 대칭적인 모양을 하고 있다는 점이 눈에 띈다. 보존계에서는 대개 $H$ 가 계의 총 에너지 $T + V$ 와 같으며, $H$ 자체가 시간에 명시적으로 의존하지 않으면 정준방정식을 따라 움직이는 동안 $H$ 의 값, 즉 에너지가 정확히 보존된다.

3. 위상공간과 리우빌 정리[편집]

해밀토니안 역학의 무대는 좌표 $q$ 와 운동량 $p$ 가 함께 만드는 위상공간이다. 자유도가 $N$ 개인 계는 $2N$ 차원 위상공간의 한 점으로 표현되고, 시간이 흐르면 이 점은 정준방정식이 그려주는 궤적을 따라 움직인다. 뉴턴 역학에서는 위치와 속도를 따로 그리지만, 위상공간에서는 이 둘이 대등한 좌표축이 된다는 점이 이후 통계역학·양자역학과의 연결고리를 이룬다.

여기서 나오는 결정적 결과가 리우빌 정리(Liouville’s theorem)다. 위상공간 위의 한 영역을 초기 조건들의 집합이라 생각하고 이 영역이 시간에 따라 흘러가는 모습을 보면, 그 부피는 늘어나지도 줄어들지도 않는다 — 마치 비압축성 유체처럼 흐른다. 위상공간의 궤적은 서로 교차하지 않으며(결정론의 기하학적 표현), 이 비압축성은 통계역학에서 앙상블 평균을 정당화하는 근거가 된다.

4. 포아송 괄호와 정준변환[편집]

두 물리량 $A(q,p)$ , $B(q,p)$ 의 시간 변화를 다루는 도구가 포아송 괄호(Poisson bracket)다.

\{A, B\} = \sum_i \left( \frac{\partial A}{\partial q_i} \frac{\partial B}{\partial p_i} - \frac{\partial A}{\partial p_i} \frac{\partial B}{\partial q_i} \right)

임의의 물리량 $A$ 의 시간 변화율은 $\dot{A} = \{A, H\}$ 로 깔끔하게 쓰인다. 흥미롭게도 이 포아송 괄호의 대수 구조는 양자역학의 교환자(commutator) $[\hat{A}, \hat{B}]$ 와 거의 같은 규칙을 따르며, 디랙은 이 유사성을 정준양자화(canonical quantization)의 출발점으로 삼았다.¹

위상공간의 좌표를 바꾸되 정준방정식의 형태와 포아송 괄호 구조를 그대로 보존하는 변환을 정준변환(canonical transformation)이라 한다. 이런 변환의 집합이 만드는 기하학적 뼈대가 바로 다음 절의 심플렉틱 구조다.

5. 심플렉틱 구조와 심플렉틱 적분기[편집]

위상공간은 그냥 좌표들의 집합이 아니라 심플렉틱 구조(symplectic structure)라는 특별한 기하를 갖는다. 정준방정식이 보존하는 이 구조 덕분에, 임의의 수치적분법이 아니라 이 구조를 존중하도록 설계된 심플렉틱 적분기(symplectic integrator)를 쓰면 놀라운 일이 벌어진다. 일반적인 룽게-쿠타법이나 오일러법으로 긴 시간 적분하면 에너지가 서서히 새어나가 궤도가 나선을 그리며 무너지지만, 심플렉틱 적분기는 원래 해밀토니안과 살짝 다른 “그림자 해밀토니안”을 정확히 보존해, 에너지 오차가 특정 범위 안에서 영원히 진동할 뿐 발산하지 않는다.

가장 널리 쓰이는 심플렉틱 적분기가 베를레 적분(leapfrog) 계열이다. 계산량은 오일러법과 별 차이가 없으면서 장기 안정성은 비교가 안 될 만큼 좋아서, 굳이 안 쓸 이유가 없다는 게 업계 중론이다.²

6. 분자동역학과 N-body 시뮬레이션에서의 활용[편집]

해밀토니안 형식이 이론에만 머물지 않는 이유는, 실제 대규모 시뮬레이션이 정확히 이 정준방정식을 적분하는 문제이기 때문이다.

분자동역학: $N$ 개 입자로 이뤄진 계의 해밀토니안은 $H = \sum_i p_i^2/2m_i + V(\mathbf{r}_1, \ldots, \mathbf{r}_N)$ 꼴이다. LAMMPS나 GROMACS 같은 MD 코드는 이 방정식을 베를레 적분으로 매 스텝 적분한다. 온도를 제어하려는 서모스탯(Nosé-Hoover 등)은 아예 가상의 자유도를 추가한 확장 해밀토니안을 만들어, 정준방정식의 틀 안에서 통계역학적 앙상블(NVT, NPT)을 구현한다.
카-파리넬로 분자동역학: 전자 밀도 자유도까지 가상의 운동량을 부여해 하나의 확장된 해밀토니안으로 원자핵과 전자를 동시에 진화시키는, 해밀토니안 형식의 창의적 확장 사례다.
랑주뱅 동역학: 여기에 마찰력과 무작위 열요동을 더해, 해밀토니안 보존계를 열욕조(heat bath)와 접촉한 계로 확장한다.
N-body 중력 시뮬레이션: 천체물리에서 은하나 항성계를 시뮬레이션할 때도 해밀토니안은 $H = \sum_i p_i^2/2m_i - \sum_{i<j} Gm_im_j/|\mathbf{r}_i - \mathbf{r}_j|$ 이다. 수백만 년에 걸친 궤도를 시뮬레이션해야 하는 특성상, 에너지 오차가 누적되면 은하 전체가 서서히 부풀거나 붕괴하는 가짜 결과가 나온다. 심플렉틱 적분기가 사실상 필수인 이유다.³

7. 관련 문서[편집]

디랙은 1925년 자신의 박사 논문에서 포아송 괄호와 양자 교환자의 유사성을 알아채고, 고전역학의 $\{A,B\}$ 를 양자역학의 $\frac{1}{i\hbar}[\hat{A}, \hat{B}]$ 로 바꾸는 “정준양자화” 규칙을 제안했다. 고전과 양자가 이렇게 대놓고 같은 문법을 공유한다는 사실은 지금 봐도 소름 돋는 우연이다. ↩
베를레 적분이 일반 룽게-쿠타보다 스텝당 정확도는 오히려 낮을 수도 있다. 그런데도 장기 시뮬레이션에서 압도적으로 선호되는 이유는 정확도가 아니라 에너지 드리프트가 없다는 구조적 성질 때문이다. “정확한데 서서히 폭발하는 방법”보다 “덜 정확해도 안 폭발하는 방법”을 고르는 게 시뮬레이션 엔지니어의 현생이다. ↩
초기 우주론 시뮬레이션에서 비심플렉틱 적분기를 오래 돌렸다가 가상의 은하가 에너지를 얻어 스스로 흩어지는 광경을 본 연구자들의 심정은 “수렴은 신에게 맡긴다”는 농담으로도 커버가 안 되는 수준이었다고 한다. ↩