Poisson bracket과 관측량의 시간 변화

위치 대신 관측량의 변화를 묻기

Hamilton 방정식은 위치와 운동량이 어떻게 변하는지 알려 준다. 실제 연구에서는 그 두 변수 자체보다 곱, 이차 특징, 학습한 보존량, 시각에 따라 정의가 달라지는 관측량이 필요할 수 있다. 매번 궤적을 대입하여 미분할 수도 있지만 공통된 연산으로 정리하면 무엇이 보존되는지, 좌표를 바꾸어도 무엇이 같은지 비교하기 쉽다.

이번 단원에서는 연쇄 법칙에서 Poisson bracket을 얻는다. 물리적 계는 앞에서 검증한 자유 입자 또는 조화진동자다. 새로 배우는 것은 궤적을 생성하는 다른 수치 적분기가 아니라, 같은 실제 운동 위에서 관측량의 변화와 좌표 표현의 관계를 계산하는 방법이다.

정준 부호와 단위를 먼저 정하기

위치 q의 단위는 m, 운동량 p는 kg m/s, 질량 m은 양수이고 강성 k는 비음수라고 하자. 해밀토니안과 운동 방정식은

$$H(q,p)=\frac{p^2}{2m}+\frac{kq^2}{2},\qquad \dot q=\frac p m,\quad \dot p=-kq$$

이다. 독립 변수의 순서를 (q,p)로 고정하고 다음과 같이 정의한다.

$$\{F,G\}=F_qG_p-F_pG_q.$$

이 부호 약속에서 {q,p}=1, {p,q}=−1이고 같은 변수끼리의 bracket은 0이다. 인자의 순서를 바꾸면 부호가 바뀐다. 두 양을 곱하는 연산, 확률변수의 공분산, 확률론의 Poisson 분포와도 다르다.

q와 p의 곱은 작용의 단위 Js다. 따라서 bracket의 단위는 [F][G]/Js다. {q,H}는 m/s, {p,H}는 N이 되어 각각 속도와 힘으로 해석할 수 있다. 두 값의 숫자가 비슷하다고 단위가 다른 성분을 그대로 합쳐 벡터 크기를 만들지는 않는다. 이번 실습의 축과 잔차도 이 구분을 유지한다.

Hamilton 흐름에 연쇄 법칙 적용하기

매끄러운 관측량 A(q,p,t)가 있다고 하자. 실제 궤적을 따라가면

$$\frac{dA}{dt}=A_q\dot q+A_p\dot p+A_t =A_qH_p-A_pH_q+A_t =\{A,H\}+\partial_t A$$

이다. 편미분 ∂tA에서는 q와 p를 고정한다. 전미분 dA/dt에서는 실제 운동에 따라 둘 다 변한다. 화면에 조회 시각 하나가 있다고 해서 이 두 미분이 같은 연산이 되는 것은 아니다.

시간에 명시적으로 의존하지 않는 A의 bracket이 관련된 모든 위상점에서 0이면 A는 보존된다. 조회한 한 점에서만 0인 것은 이 항등식을 뜻하지 않는다. 명시적인 시간 항이 있는 A는 bracket과 편미분을 합한 값이 0인지 확인해야 한다. 특히 bracket만 작게 만드는 학습 손실이 언제나 올바른 보존량 조건은 아니다.

같은 작용 단위의 네 관측량 선택하기

여러 관측량의 값을 같은 축에서 비교하려면 먼저 단위를 맞춘다. 진동수가 양수이면 기준 진동수 ν=Ω, 자유 운동이면 ν=1 s⁻¹로 둔다. 선택할 기본 관측량은

$$C=qp,\qquad F=\frac{m\nu q^2}{2},\qquad G=\frac{p^2}{2m\nu},\qquad B=\frac H\nu$$

이며 모두 Js다. ν는 선택한 실험에서 고정한 기준값이다. q나 p에 따라 다시 변하는 함수로 미분하지 않는다. 직접 미분하면

$$\{C,H\}=\frac{p^2}{m}-kq^2,\qquad \{F,H\}=\nu qp,$$ $$\{G,H\}=-\frac{kqp}{m\nu},\qquad \{B,H\}=0$$

을 얻는다. 자유 입자에서는 위치에너지가 없어 G와 B가 같은 함수다. 메뉴에 네 항목이 있다고 언제나 네 개의 독립 함수가 생기는 것은 아니다.

명시적 시간 의존성을 조절하기 위해 A=A_base+βtH를 사용한다. β는 무차원이다. A는 Js, bracket과 명시적 시간 미분 및 전체 시간 변화율은 J다. 이 변화율은 작용 단위의 관측량이 변하는 속도이며 곧바로 실제 에너지 손실을 뜻하지 않는다.

유리수로 계산되는 자유 운동 예제

m=1 kg, k=0, 초기 q0=1 m, v0=1 m/s로 두자. t=1 s에서 위치는 2 m, 운동량은 1 kg m/s, 에너지는 0.5 J다. 기본 관측량 C를 고르고 β=0이면 A=2 Js이고 bracket 변화율은 p²/m=1 J다. 물리적 힘은 0이지만 위치가 변하므로 위치와 운동량의 곱도 변한다.

같은 물리적 운동에서 β=−2를 선택하면

$$A=qp-2tH=qp-tp^2/m=q_0p$$

가 된다. 이 식은 특정한 숫자 예제뿐 아니라 모든 자유 입자 궤적에서 일정하다. t=1 s에서 기본 항은 2 Js, 명시적 시간 항은 −1 Js이므로 최종 A는 1 Js다. bracket은 여전히 1 J이지만 명시적 편미분이 −1 J이어서 전체 변화율은 0이다.

따라서 시간에 명시적으로 의존하는 보존량은 H와의 bracket만 0일 필요가 없다. 시각 항을 빼먹으면 올바른 보존량을 잘못 기각한다. 반대로 잘못 학습한 관측량에 임의의 시각 항을 붙여 데이터 구간에서만 변화를 상쇄했다고 일반적인 물리적 보존 법칙을 얻은 것도 아니다.

변화율이 한 번 0이어도 보존은 아니다

m=1 kg, Ω=π/2 s⁻¹, 초기 q0=1 m, v0=0인 진동자를 생각하자. 실제 위치는 cos(Ωt), 운동량은 −Ωsin(Ωt), 에너지는 π²/8 J다.

A=qp, β=0이면 t=0과 t=1 s에서 A=0이지만 중간인 t=0.5 s에서 A=−π/4 Js다. 바로 이 중간 시각의 변화율은 0이다. 운동에너지와 위치에너지가 같아져 p²/m−kq²가 상쇄되기 때문이다. 한 점의 정지한 변화율과 구간 전체의 일정한 값을 구분해야 한다.

반면 B=H/ν는 같은 궤적 전체에서 π/4 Js로 일정하다. 조회 시각을 옮기는 일은 같은 미리 계산한 궤적에서 다른 상태를 선택하는 것이며, 궤적 자체를 다시 바꾸는 일이 아니다.

함수 전체의 보존 조건과 한 궤적의 우연한 일정함도 구분한다. 정지한 자유 입자에서는 F도 일정하지만, 자유 입자의 전체 위상 공간에서 {F,H}=νqp가 항등적으로 0인 것은 아니다. 실습은 이 두 판정을 별도로 제공한다.

정준 좌표의 올바른 크기 변경

a>0을 무차원 배율로 정하고

$$Q=aq,\qquad P=p/a$$

로 좌표를 바꾸자. 원래 bracket으로 {Q,P}를 계산하면 1이다. 새 독립 변수로 쓴 해밀토니안은

$$\widetilde H(Q,P)=\frac{a^2P^2}{2m}+\frac{kQ^2}{2a^2}$$

이다. 관측량도 같은 물리적 함수가 되도록 좌표를 대입한다. 연쇄 법칙에서 A_Q=A_q/a, A_P=a A_p를 얻으며, (Q,P)의 표준 bracket에 넣으면 배율이 상쇄된다. 따라서 좌표를 바꾸기 전후의 물리적 bracket 값은 같다. 소수 계산에서는 부동소수점 평가 차이만 남을 수 있다.

Q는 여전히 길이, P는 운동량 단위를 갖는다. 그러나 좌표의 숫자와 위상 초상은 바뀐다. 두 위상 곡선을 공간에서 움직이는 두 입자의 궤적으로 해석하지 않는다. 같은 물리적 양을 표현하는 것과 좌표 성분의 숫자가 같은 것은 다른 주장이다.

두 변수를 모두 곱하면 생기는 오류

Q=aq, Pbad=ap처럼 둘 다 같은 배율로 바꾸는 유혹을 생각하자. 원래 정준 구조에서 계산하면

$$\{Q,P_{\mathrm{bad}}\}=a^2$$

이다. 새 좌표에서도 기본 bracket이 1이라고 잘못 가정하면 관측량의 변화율은 실제 값의 1/a²가 된다. 위의 자유 운동에서 a=2이면 실제 Q 변화율은 2 m/s인데 잘못 계산한 값은 0.5 m/s다.

좌표 변환 자체는 가역이다. 잘못된 점은 변환한 뒤에도 표준 정준 bracket을 아무 조건 없이 사용한 것이다. 이 좌표가 유도하는 Poisson 행렬을 함께 사용하면 같은 운동을 표현할 수 있다.

$$J_{\mathrm{bad}}=a^2\begin{pmatrix}0&1\\-1&0\end{pmatrix}.$$

즉 잘못된 표준 bracket 계산에 a²를 곱하면 원래 값을 복원한다. a=1이면 차이가 사라진다. 여기서는 한 정준 쌍의 일정한 선형 배율만 다룬다. 이 예제를 고차원에서 부피 보존만으로 모든 정준 구조가 보존된다는 주장으로 확장하지 않는다. 일반적인 symplectic 구조는 뒤 단원에서 따로 다룬다.

함수들의 대수로 읽기

Bracket은 선형성과 미분의 곱 규칙을 따른다.

$$\{F,GH\}=\{F,G\}H+G\{F,H\}.$$

선택한 작용 단위의 이차 함수에서는

$$\{F,G\}=C,\qquad \{G,C\}=-2G,\qquad \{C,F\}=-2F$$

라는 닫힌 계산 예제를 얻는다. 이 세 함수의 순환 중첩 bracket은 각각 −2C, +2C, 0으로 합이 0이다. Jacobi 항등식을 손으로 검사하기 좋은 예다.

다만 한 함수 조합에서의 수치적 상쇄가 모든 매끄러운 함수에 대한 증명은 아니다. 일반식에서는 미분을 전개하고 이계 편미분의 대칭성과 bracket의 반대칭성이 짝을 이루어 상쇄되는 것을 보인다. 필요한 미분 가능성과 정의역을 먼저 명시해야 한다.

학습한 bracket 연산자가 몇 개의 관측된 변화율과 잘 맞더라도 다른 지점에서 반대칭성, 곱 규칙, Jacobi 조건을 어길 수 있다. 따라서 궤적 적합, 보존량 조건, 대수 구조의 검증을 서로 다른 항목으로 기록한다.

계산 출력과 오차를 해석하기

모델은 두 초 동안 201개 시점을 계산한다. Ω=πf/4, k=mΩ²이므로 f를 고정하고 질량을 바꾸면 강성도 함께 바뀐다. 자유 극한은 정확한 선형 흐름을 사용하고 특별한 사분면 시각의 삼각함수 값은 앞 단원의 해석 흐름 구현을 재사용한다.

관측량 이력은 Js, 변화율은 J, 위상 공간의 축은 위치 m와 운동량 kg m/s다. 항등식 잔차는 계산한 소수 차이를 그대로 유지한다. 보고한 최대 잔차는 출력 시점의 표본 최대이며 모든 실수 시각에 대한 오차 상한이나 새 ODE 적분기의 오차가 아니다.

보존 여부의 표시는 그래프가 거의 평평한지 보고 정하지 않는다. 선택한 함수와 매개변수의 정확한 조건에서 판정한다. 원점은 여기서 직교 좌표의 정상적인 상태이므로 앞 극좌표 단원의 좌표 미정과도 구분한다.

해설이 있는 연습 문제

1. {q,H}와 {p,H}를 계산하자. 해설: q_q=1과 p_p=1인 항만 남아 각각 p/m과 −kq다. Hamilton 방정식과 비교하면 bracket 인자 순서의 부호를 점검할 수 있다.
2. 운동·위치에너지가 같은 한 시각에서 {qp,H}=0이면 왜 보존을 말할 수 없을까? 해설: 이 식은 2K−2U다. 두 에너지가 같은 시각은 고립되어 나타날 수 있고 그 전후에 qp가 변한다. 반 초 진동 예제에서 직접 확인할 수 있다.
3. 자유 운동의 명시적 시간 보존량을 검산하자. 해설: k=0에서 {qp,H}=2H이고 ∂t(−2tH)=−2H다. 합은 모든 위상점과 시각에서 0이지만 bracket 자체는 비영일 수 있다.
4. a=2의 잘못된 좌표 변경은 변화율을 얼마나 바꿀까? 해설: 새 기본 bracket은 4다. 이를 1로 취급하면 변화율은 원래의 1/4이 되며 유도된 Poisson 행렬을 사용하면 복원된다.
5. q=2 m, p=1 kg m/s에서 이차 함수의 Jacobi 예제를 계산하자. 해설: C=2 Js이므로 순환 항은 −4 Js, +4 Js, 0이다. 이 상쇄는 해당 예제의 검사이며 임의 함수 전체에 대한 증명은 아니다.
6. β=1이면 모든 궤적에서 A가 변한다고 할 수 있을까? 해설: 아니다. 에너지가 0인 정지 자유 궤적에서는 H=0이어서 시각 항이 사라진다. 선택한 궤적 밖에서 함수의 명시적 시간 의존성은 여전히 존재한다.

AI 연구와 읽을거리

물리적 보존량을 학습할 때는 명시적 시각 입력과 정준 좌표의 약속도 가설에 포함해야 한다. 궤적 손실에 더해 bracket 항등식과 독립적인 시간 미분을 대조한다. 보존량이라는 주장, 정준 좌표라는 주장, 작은 표본 잔차라는 주장을 각각 분리해서 평가하는 습관이 필요하다.

표준 부호와 정준 변환의 틀은 David Tong의 Hamiltonian Formalism 4.3·4.4절을 참고할 수 있다. 이 단원의 작용 단위 관측량, 잘못된 배율 비교, 유리수 시간 의존 예제는 직접 유도했으며 독립 기호 미분·60자리 궤적 미분·정확 분수 기준과 대조한다.