해밀턴 역학 — 위상공간으로의 전환

르장드르 변환이 라그랑지언을 해밀토니언으로 바꾸고, 무대는 배위공간 MM 에서 여접다발 TMT^*M 로 옮겨가며, 운동은 2n2n 개의 1계 ODE 인 해밀턴 방정식으로 다시 쓰인다.

들어가며

7–9장의 무대는 접다발(tangent bundle) TMTM — 배위공간 MM 의 각 점 qq 에 속도 q˙\dot q 를 붙여 만든 공간 — 이었고, 라그랑지언 L(q,q˙,t)L(q, \dot q, t) 가 그 위에서 사는 함수였다.

이 장은 속도 q˙\dot q 를 켤레운동량 pp 로 갈아 끼운다. 좌표 qq 는 그대로 두고 짝꿍만 바꾸면 무대는 TMTM 에서 여접다발(cotangent bundle) TMT^*M 로 옮겨 가고, 라그랑주 쪽의 nn 개 2계 ODE 는 2n2n 개의 1계 ODE — 좌표와 운동량이 거의 대칭으로 등장하는 해밀턴 방정식 — 로 다시 쓰인다.

속도 언어에서는 잘 안 보이던 (q,p)(q, p) 의 짝 구조가 운동량 언어에서는 한눈에 들어온다. 11장 정준변환은 이 쌍을 보존하는 좌표 변환을, 12장 푸아송 괄호는 이 쌍 위의 새로운 곱셈을 다룬다.

본론 1은 켤레운동량과 르장드르 변환. 본론 2는 해밀턴 방정식과 위상공간 TMT^*M. 본론 3은 단진자의 위상초상 — 적분이 안 풀리는 비선형계의 운동 전모를 보존량의 등고선으로 읽는다.

본론 1 — 켤레운동량과 르장드르 변환

라그랑지언 L(q,q˙,t)L(q, \dot q, t) 가 주어졌을 때 ii 번째 좌표 qiq^i켤레운동량(conjugate momentum — 라그랑지언을 ii 번째 속도 슬롯으로 편미분한 양)

pi  =  Lq˙i.p_i \;=\; \frac{\partial L}{\partial \dot q^i}.

9장 뇌터 보존량 J=(L/q˙i)ξiJ = (\partial L/\partial\dot q^i)\xi^i 안의 L/q˙i\partial L/\partial\dot q^i 가 바로 이 양이다 — 거기서는 “일반화 운동량”이라 불렀다.

L=12mq˙2U(q)L = \tfrac12 m \dot q^2 - U(q)p=mq˙p = m\dot q — 익숙한 선운동량이다. 그러나 좌표가 각도면 pp 는 각운동량이 되고, 회전하는 좌표계면 회전 보정항이 끼며, 자기장 안의 하전입자면 mq˙m\dot q 에 벡터퍼텐셜 항이 더해진다. 켤레운동량은 언제나 mq˙m\dot q 가 아니라 정의식 L/q˙i\partial L/\partial\dot q^i 의 결과다.

헷갈리기 쉬운 곳 — 켤레운동량 \ne mq˙m\dot q 가 일반적이다. L=12mq˙2UL = \tfrac12 m\dot q^2 - U 라는 첫 예제가 마침 p=mq˙p = m\dot q 를 주기 때문에 그 등식이 정의처럼 머리에 박힌다. 라그랑지언에 속도가 일차로 들어가는 항이 있거나 운동에너지가 좌표에 의존하는 계수를 가지면 — 극좌표의 12m(r˙2+r2θ˙2)\tfrac12 m(\dot r^2 + r^2\dot\theta^2) 에서 θ\theta 의 켤레운동량은 mr2θ˙m r^2\dot\theta — 결과가 단순한 mq˙m\dot q 와 갈라진다. 매번 정의식 L/q˙i\partial L/\partial\dot q^i 를 직접 계산하라.

(q,q˙)(q, \dot q)(q,p)(q, p) 로 바꾸려면 pi=L/q˙ip_i = \partial L/\partial\dot q^iq˙\dot q 에 대해 거꾸로 풀 수 있어야 한다. 가능 조건이 정칙성(regularity — 헤시안 행렬이 모든 점에서 가역이라는 조건) 이다. 헤시안

2Lq˙iq˙j\frac{\partial^2 L}{\partial \dot q^i\, \partial \dot q^j}

이 가역이면 음함수 정리가 q˙i=q˙i(q,p,t)\dot q^i = \dot q^i(q, p, t) 라는 유일한 풀이를 보장한다.

헷갈리기 쉬운 곳 — 르장드르 변환은 정보 손실 없는 가역 변환이다 (정칙성 가정 하). 헤시안이 가역이면 q˙p\dot q \to p 대응이 일대일이라 pp 에서 q˙\dot q 를 정확히 되찾는다. 직관적으로는 — 한 곡선을 “점들의 모임” 대신 “그 곡선에 닿는 접선들의 모임”으로 다시 적는 조작이다. 곡선이 볼록한 한 두 사전은 완벽하게 양방향이고, 그 “볼록함”이 헤시안 가역성이다. 정칙성이 깨지면(특이 라그랑지언) 사전이 한쪽에서 막힌다 — 이 장은 정칙인 경우만 상대한다.

해밀토니언(Hamiltonian) 의 정의는

H(q,p,t)  =  piq˙i(q,p,t)    L(q,q˙(q,p,t),t)H(q, p, t) \;=\; p_i\, \dot q^i(q, p, t) \;-\; L\bigl(q,\, \dot q(q, p, t),\, t\bigr)

이고(아인슈타인 합 규약), 이것이 LLq˙\dot q 에 대한 르장드르 변환(Legendre transform) 이다. 우변에 q˙\dot q 가 보이지만 모두 q˙i(q,p,t)\dot q^i(q, p, t) 형태로 (q,p)(q, p) 의 함수로 치환된 채 들어 있어 전체는 (q,p)(q, p) 의 함수다.

헷갈리기 쉬운 곳 — HH(q,p)(q, p) 의 함수다 — q˙\dot q 를 모두 소거해야 한다. H=pq˙LH = p\dot q - L 이라고 적어 놓고 우변에 q˙\dot q 를 그대로 남겨 두면 HH(q,q˙,p)(q, \dot q, p) 의 함수가 되어 르장드르 변환의 의미가 사라진다. 절차는 두 단계다 — 먼저 p=L/q˙p = \partial L/\partial\dot q 를 풀어 q˙=q˙(q,p)\dot q = \dot q(q, p) 를 얻고, 그것을 pq˙Lp\dot q - L 의 모든 q˙\dot q 자리에 대입해 q˙\dot q 가 식에서 완전히 사라지게 만든다. 최종 결과에 q˙\dot q 가 한 글자라도 남아 있으면 변환이 끝나지 않은 것이다.

직관. L(q˙)L(\dot q) 를 곡선으로 보면(qq 고정), 곡선 위 한 점의 기울기p=L/q˙p = \partial L/\partial\dot q 다. 기울기만으로는 곡선을 복원할 수 없으니 — 같은 기울기의 직선은 무수히 많으니 — 그 접선이 세로축을 자르는 절편 (pq˙L)=Lpq˙-(p\dot q - L) = L - p\dot q 도 함께 기록한다. 부호를 뒤집은 H=pq˙LH = p\dot q - L 이 해밀토니언이다. 곡선이 볼록한 한 점의 모임과 접선의 모임은 같은 곡선의 두 사전이다.

손계산. L=12mq˙2U(q)L = \tfrac12 m\dot q^2 - U(q) 에서 p=mq˙p = m\dot qq˙=p/m\dot q = p/m → 대입:

H=ppm(12mp2m2U)=p2mp22m+U=p22m+U(q).H = p\cdot \frac{p}{m} - \Bigl(\frac12 m \frac{p^2}{m^2} - U\Bigr) = \frac{p^2}{m} - \frac{p^2}{2m} + U = \frac{p^2}{2m} + U(q).

운동에너지와 위치에너지의 합 — 총 에너지다. 우변에 q˙\dot q 가 없다는 점이 변환의 완료 표시다. 단순한 경우 HH 가 곧 에너지지만 회전 좌표계나 자기장이 끼면 둘이 갈라진다(12장).

본론 2 — 해밀턴 방정식과 위상공간 TMT^*M

H=piq˙iLH = p_i\dot q^i - L 과 오일러–라그랑주 방정식을 양변 미분해 정리하면 HH 의 편미분 두 개가 곧장 운동을 내놓는다.

q˙i=Hpi,p˙i=Hqi.\dot q^i = \frac{\partial H}{\partial p_i}, \qquad \dot p_i = -\frac{\partial H}{\partial q^i}.

이 한 쌍이 해밀턴 방정식(Hamilton’s equations) 이다. 라그랑주 쪽의 nn 개 2계 방정식이 2n2n 개의 1계 방정식 — ¨\ddot{} 가 없고 q˙,p˙\dot q, \dot p 만 들어가는 — 으로 바뀌었다. 1계 시스템은 현재 상태 (q,p)(q, p) 한 점만 알면 그 점의 (q˙,p˙)(\dot q, \dot p) 가 즉시 나와 다음 점을 가리키는 구조 — 위상공간 모든 점에 화살표가 하나씩 박힌 흐름의 그림이다.

q˙i=H/pi\dot q^i = \partial H/\partial p_ip˙i=H/qi\dot p_i = -\partial H/\partial q^i거의-대칭성 도 눈여겨볼 만하다. 차이는 부호 하나뿐이다. 라그랑주 그림에서 qq 가 주인공이고 q˙\dot q 가 종속물이었던 것이, 해밀턴 그림에서는 qqpp 가 거의 동등한 두 배우가 된다. 이 부호 하나의 대칭이 11장 정준변환과 12장 푸아송 괄호의 토대다.

라그랑주 그림에서 q˙\dot q 는 좌표를 시간으로 미분해 파생되는 양이라 — 좌표를 자유롭게 바꿔도 속도까지 마음대로 섞는 변환은 의미가 모호하다. 해밀턴 그림에서 q,pq, p 는 둘 다 위상공간의 독립 좌표라 — 새 좌표를 옛 운동량으로, 새 운동량을 옛 좌표로 잡는 섞임 변환까지 정당한 좌표 변환이 된다. 이 자유로움이 정준변환의 무대다.

헷갈리기 쉬운 곳 — q˙=H/p\dot q = \partial H/\partial p 는 정의가 아니라 운동방정식이다. 르장드르 변환 단계의 q˙(q,p)\dot q(q, p) 는 켤레운동량 정의식을 거꾸로 푼 대수적 관계 이고, 어떤 운동이 일어나든 항상 성립하는 항등식이다. 해밀턴 방정식의 q˙=H/p\dot q = \partial H/\partial p물리적 궤적 위에서만 성립하는 운동방정식이다 — 9장의 표현으로 on-shell에서만 참이다. 글자가 닮았지만 위상이 다르다. 진짜 내용물은 p˙=H/q\dot p = -\partial H/\partial q 쪽 — 이 식이 오일러–라그랑주의 해밀턴판이다.

무대도 함께 바뀐다. 라그랑주 그림에서 q˙\dot q 가 살던 곳은 접공간 TqMT_qM 이었다. 해밀턴 그림에서 pp 가 사는 곳은 여접공간(cotangent space — 접공간의 쌍대공간, 즉 접벡터를 받아 실수를 돌려주는 선형범함수들의 모임) TqMT_q^*M 이다. 켤레운동량을 미분형식으로 p=pidqip = p_i\, dq^i 로 적으면 — dqidq^i 가 접벡터의 ii 번째 성분을 뽑는 1-형식이니 — 그 선형결합이 TqMT_q^*M 의 원소다.

모든 점의 여접공간을 모은 다양체 여접다발(cotangent bundle) TMT^*M 이 해밀턴 역학의 위상공간(phase space) 이다. dimM=n\dim M = n 이면 각 점에 nn 차원 여접공간이 붙으니 TMT^*M2n2n 차원 — 해밀턴 방정식의 2n2n 개 1계 방정식과 정확히 맞는다. 한 점 (q,p)(q, p) 가 계의 완전한 순간 상태 다.

헷갈리기 쉬운 곳 — 위상공간 TMT^*M 과 상태공간 TMTM 의 차이. TMTM 의 한 점은 (q,q˙)(q, \dot q)(상태공간), TMT^*M 의 한 점은 (q,p)(q, p)(위상공간)다. 정칙계에서는 르장드르 변환으로 둘 사이를 오갈 수 있지만 좌표 변환에 따라 대응이 복잡하게 얽혀 두 공간을 “같다”고 말하면 안 된다. 결정적 차이는 — TMT^*M 에는 좌표와 무관하게 미리 새겨진 기하 구조가 있고 TMTM 에는 그런 구조가 자연스럽게 없다는 점이다. 운동량이 1-형식이라는 사실이 그 구조의 출발점이다. 단순히 변수 이름을 바꾼 게 아니라 더 풍부한 무대로 옮긴 것이다.

TMT^*M 위의 자연스러운 2-형식

ω  =  dqidpi\omega \;=\; dq^i \wedge dp_i

심플렉틱 형식(symplectic form) 이다. 두 가지만 알아 두자. (1) 해밀턴 방정식은 ω\omegaHH 가 만나서 만드는 벡터장의 흐름으로 다시 적을 수 있고, q˙=H/p\dot q = \partial H/\partial p, p˙=H/q\dot p = -\partial H/\partial q 의 부호 비대칭이 ω\omega 의 쐐기곱이 가진 반대칭성에서 온다. (2) ω\omega 는 좌표를 어떻게 잡든 같은 모양으로 새겨져 있다.

헷갈리기 쉬운 곳 — 심플렉틱 형식은 좌표와 무관하게 새겨져 있다. ω=dqidpi\omega = dq^i \wedge dp_i 가 특정 좌표에 묶여 보이지만, 좌표를 정준변환으로 (Q,P)(Q, P) 로 바꿔도 ω=dQidPi\omega = dQ^i \wedge dP_i 라는 똑같은 모양 으로 다시 적힌다. ω\omega 는 위상공간 다양체 자체에 박힌 기하 구조다 — 지구 표면의 곡률이 위도·경도를 어떻게 그리든 변하지 않는 것과 같다. 정준변환 장에서 왜 ω\omega 가 보존되는지를 본격적으로 본다.

본론 3 — 단진자의 위상초상

단진자의 좌표는 막대가 연직선과 이루는 각 θ\theta 하나, 좌표공간은 [π,π][-\pi, \pi] 의 양 끝을 동일시한 원. 라그랑지언 L=12m2θ˙2+mgcosθL = \tfrac12 m\ell^2 \dot\theta^2 + m g \ell \cos\theta 에서 켤레운동량은 p=L/θ˙=m2θ˙p = \partial L/\partial\dot\theta = m\ell^2 \dot\theta — 각운동량이다. 르장드르 변환을 세 단계로 굴리면 — p=m2θ˙p = m\ell^2\dot\thetaθ˙=p/(m2)\dot\theta = p/(m\ell^2) → 대입 —

H(θ,p)=p22m2mgcosθ.H(\theta, p) = \frac{p^2}{2 m \ell^2} - m g \ell \cos\theta.

HH 가 시간에 명시적으로 의존하지 않으니(9장 본론 2(c)) 궤적을 따라 보존된다 — 즉 각 궤적은 H=EH = E 등고선 위에 갇혀 머문다. 위상공간 (θ,p)(\theta, p) 위의 HH 등고선이 곧 가능한 모든 궤적의 지도다.

헷갈리기 쉬운 곳 — HH 보존 ⟹ 궤적은 등고선이다 — 적분을 못 해도 그림은 그린다. 단진자 운동방정식 θ¨=(g/)sinθ\ddot\theta = -(g/\ell)\sin\theta 는 초등함수로 풀리지 않고 — 해를 적자면 타원적분이 필요하다. 그러나 H(θ,p)=EH(\theta, p) = E 는 적분이 아니라 대수 방정식이다. θ\theta 를 주면 pp 가 풀리고, 그 점들을 이으면 궤적의 모양이 나온다. 타이밍 은 적분이 필요하지만 모양 자체는 보존량의 등고선만으로 공짜다. “풀 수 없는 문제”와 “그릴 수 없는 문제”는 다르다.

EE 의 값에 따라 등고선의 모양이 세 가지로 갈린다.

첫째, E<mgE < m g \ell. 진자가 아래쪽 평형점 θ=0\theta = 0 주변을 왔다 갔다 한다 — 위로 한 바퀴 넘어갈 에너지가 부족하다. 위상공간에서는 원점 (θ,p)=(0,0)(\theta, p) = (0, 0) 을 둘러싸는 닫힌 곡선진동(oscillation) — 이다. EE 가 작을수록 곡선이 원점에 가까운 작은 타원, EE 가 커지면 부풀어 오른다.

둘째, E=mgE = m g \ell. 진자가 정확히 위쪽 꼭대기 θ=±π\theta = \pm\pi(불안정 평형점)까지 도달할 수 있는 에너지다. 위상공간에서는 위쪽 평형점에 점근하는 한 쌍의 곡선 — 진동 영역과 회전 영역을 가르는 분리선(separatrix) — 이 된다.

헷갈리기 쉬운 곳 — 분리선 도달은 무한 시간이 걸린다. 진자는 위쪽 꼭대기에 영원히 도달하지 못한다. 꼭대기 θ=π\theta = \pi 가 불안정 평형점이라 가까워질수록 복원력이 약해지고 속도도 0 으로 줄어들어 — 점점 느려지며 한없이 다가가기만 한다. 분리선이 평형점에 “점근한다”는 말이 이 뜻이다. 평형점 자체는 곡선에 포함되지 않는다. 분리선은 실제 관측되는 운동이라기보다 진동과 회전이라는 두 세계의 경계선 으로 의미가 있다.

분리선의 식. H=mgH = m g \ell 에서 p22m2mgcosθ=mg\dfrac{p^2}{2m\ell^2} - mg\ell\cos\theta = mg\ell, p2=2m23g(1+cosθ)p^2 = 2m^2\ell^3 g(1 + \cos\theta), 따라서

p=±2m23g(1+cosθ).p = \pm \sqrt{2 m^2 \ell^3 g\,(1 + \cos\theta)}.

±\pm 는 시계방향·반시계방향에 대응한다. 단위계 m==g=1m = \ell = g = 1 에서는 p=±2(1+cosθ)p = \pm\sqrt{2(1 + \cos\theta)} 다.

셋째, E>mgE > m g \ell. 진자가 위쪽 꼭대기를 넘어갈 에너지를 가져 멈추지 않고 한 방향으로 계속 막대를 돌린다 — 회전(rotation). 위상공간에서는 닫히지 않고 θ\theta 축을 따라 흐르는 열린 곡선 이다. p>0p > 0 이면 θ\theta 증가, p<0p < 0 이면 감소.

헷갈리기 쉬운 곳 — θ=±π\theta = \pm\pi 동일시 — 위상공간이 원통(cylinder)이다. 회전 운동의 곡선이 θ\theta 축을 따라 끝없이 흘러가 보이는 함정은 위상공간을 평면 으로 그린 데 있다. 단진자의 θ\theta 는 원이라 θ=π\theta = \piθ=π\theta = -\pi 가 같은 상태고, 위상공간 (θ,p)(\theta, p) 는 평면이 아니라 θ\theta 방향 양 끝을 풀로 붙인 원통 이다. 평면 그림은 원통을 가위로 잘라 펼친 것이고, 오른쪽 끝 θ=π\theta = \pi 로 흘러 나간 곡선은 사실 왼쪽 끝 θ=π\theta = -\pi 로 다시 들어온다. 원통 위에서는 회전 운동도 어엿한 닫힌 곡선 — 원통을 한 바퀴 감는 고리 — 이다. 진동은 원통의 한 면에 난 작은 고리, 회전은 원통 전체를 두르는 큰 고리, 분리선은 그 둘을 가르는 임계 고리.

진동·분리선·회전 세 영역이 모두 하나의 함수 H(θ,p)H(\theta, p) 의 등고선이고, 에너지 EE 다이얼 하나로 세 세계 사이를 연속으로 옮겨 다닐 수 있다.

파이썬으로 확인

# 단진자의 위상초상: 해밀턴 방정식을 손작성 RK4 로 적분하고
# (theta, p) 평면 위에 15 개 궤적을 그린다. 분리선 E=1 은 점선으로 비교.
import numpy as np
import matplotlib.pyplot as plt

def rhs(s):                                # 해밀턴 방정식 (m=l=g=1)
    th, p = s                              # dH/dp = p, -dH/dth = -sin(th)
    return np.array([p, -np.sin(th)])

def rk4_step(s, dt):
    k1 = rhs(s)
    k2 = rhs(s + 0.5 * dt * k1)
    k3 = rhs(s + 0.5 * dt * k2)
    k4 = rhs(s + dt * k3)
    return s + dt * (k1 + 2*k2 + 2*k3 + k4) / 6

dt, T = 0.02, 20.0
N = int(T / dt)

fig, ax = plt.subplots(figsize=(7, 4))
for E in (0.2, 0.6, 1.0, 1.4, 1.8):        # 에너지를 격자처럼 훑고
    for th0 in (-2.0, 0.0, 2.0):           # 초기 각도 세 개
        val = 2 * (E + np.cos(th0))        # p0^2 = 2(E + cos th0)
        if val < 0:                        # 운동 불가 영역은 건너뛴다
            continue
        s = np.array([th0, np.sqrt(val)])
        traj = np.empty((N + 1, 2)); traj[0] = s
        for n in range(N):
            traj[n+1] = rk4_step(traj[n], dt)
        ax.plot(traj[:, 0], traj[:, 1], lw=0.7)

th = np.linspace(-np.pi, np.pi, 400)        # 분리선 E=1
ax.plot(th,  np.sqrt(2*(1 + np.cos(th))), "k--", lw=1)
ax.plot(th, -np.sqrt(2*(1 + np.cos(th))), "k--", lw=1)
ax.set_xlabel(r"$\theta$"); ax.set_ylabel(r"$p$"); plt.show()

rhs 는 단위계 m==g=1m = \ell = g = 1 의 해밀턴 방정식 — θ˙=H/p=p\dot\theta = \partial H/\partial p = p, p˙=H/θ=sinθ\dot p = -\partial H/\partial\theta = -\sin\theta — 을 두 줄로 직역한 것이다. rk4_step 은 표준 4차 룽게–쿠타 한 스텝. 1계 시스템이라 상태 (θ,p)(\theta, p) 한 점만 넣으면 다음 점이 나온다.

이중 루프. 바깥은 에너지 EE0.20.2 부터 1.81.8 까지 다섯 값으로 — E=0.2,0.6E = 0.2, 0.6 은 분리선 아래(진동), E=1.0E = 1.0 은 분리선, E=1.4,1.8E = 1.4, 1.8 은 분리선 위(회전). 안쪽은 초기 각 θ0\theta_0 세 개. val = 2 * (E + cos th0)H=p2/2cosθ=EH = p^2/2 - \cos\theta = Ep2p^2 에 대해 푼 것이고, 음수면 그 에너지로 그 위치가 금지 영역 — 위치에너지가 이미 EE 를 넘은 곳 — 이라 건너뛴다.

각 초기 조건에서 dt=0.02dt = 0.02T=20T = 20 까지 굴리고 (θ,p)(\theta, p) 자취를 그린다. 마지막 세 줄은 분리선 p=±2(1+cosθ)p = \pm\sqrt{2(1 + \cos\theta)} 를 검은 점선으로 겹친다. 결과 — 원점을 둘러싼 진동의 닫힌 곡선, 위아래로 흘러가는 회전의 열린 곡선, 그 둘을 정확히 가르며 위쪽 평형점으로 점근하는 점선. 적분으로 풀리지 않는 비선형 진자의 운동 전체가 보존량 HH 의 등고선 한 장으로 완전히 정리된다.

다음 장으로

11장: 정준변환과 심플렉틱 구조에서는 본론 2 끝에서 이름만 던졌던 심플렉틱 형식 ω=dqidpi\omega = dq^i \wedge dp_i 가 좌표 변환에 어떻게 반응하는지를 본격적으로 다룬다. 위상공간 TMT^*M 위의 좌표를 (q,p)(q, p) 에서 (Q,P)(Q, P) 로 바꿀 때 ω\omega 를 그대로 보존하는 변환이 정준변환(canonical transformation) 이고, 이 클래스만이 해밀턴 방정식의 1계 대칭 형태를 그대로 보존한다.