포아송 괄호

1. 개요[편집]

포아송 괄호(Poisson bracket)는 위상공간 위에서 정의된 두 물리량 $f(\mathbf{q},\mathbf{p},t)$, $g(\mathbf{q},\mathbf{p},t)$에 대해

$$\{f,g\} = \sum_{i=1}^{N}\left(\frac{\partial f}{\partial q_i}\frac{\partial g}{\partial p_i} - \frac{\partial f}{\partial p_i}\frac{\partial g}{\partial q_i}\right)$$

로 정의되는 이항 연산이다. 해밀토니안 역학의 문법 그 자체라고 보면 된다. 운동방정식, 보존량, 대칭성, 정준변환, 그리고 양자역학의 교환자까지 전부 이 괄호 하나의 언어로 다시 쓸 수 있다.¹

시뮬레이션 하는 사람 입장에서 이 문서가 중요한 이유는 딱 하나다. 왜 베를레 적분 같은 심플렉틱 적분기가 룽게-쿠타법보다 장시간 궤도 계산에서 압도적으로 잘 버티는가에 대한 답이 포아송 괄호가 정의하는 기하 구조에 들어 있기 때문이다. 정밀도가 높아서가 아니다. 지키는 게 달라서다.

2. 정의와 대수적 성질[편집]

정의에서 곧바로 따라 나오는 성질들이 이 연산의 전부다.

• 쌍선형성: $\{af+bg,h\} = a\{f,h\} + b\{g,h\}$.
• 반대칭성: $\{f,g\} = -\{g,f\}$. 따라서 $\{f,f\}=0$.
• 라이프니츠 법칙: $\{fg,h\} = f\{g,h\} + g\{f,h\}$. 미분 연산자처럼 곱에 분배된다.
• 야코비 항등식: $\{f,\{g,h\}\} + \{g,\{h,f\}\} + \{h,\{f,g\}\} = 0$.

앞의 세 개는 계산해 보면 당연하지만, 야코비 항등식은 손으로 전개하면 항이 24개까지 불어나면서 기적처럼 상쇄된다.² 이 네 성질을 갖춘 대수 구조를 리 대수(Lie algebra)라 부르며, 위상공간 위의 매끄러운 함수 전체가 포아송 괄호를 곱으로 삼아 무한차원 리 대수를 이룬다. 리 군과 대칭성 이론이 고전역학에 들어오는 통로가 바로 여기다.

정준좌표 자신에 대해서는 아주 깔끔한 관계가 성립한다.

$$\{q_i,q_j\}=0,\qquad \{p_i,p_j\}=0,\qquad \{q_i,p_j\}=\delta_{ij}$$

이를 기본 포아송 괄호라 하고, 사실상 “정준좌표란 무엇인가”의 정의로 쓰인다.

3. 시간 발전과 보존량[편집]

포아송 괄호의 진짜 위력은 운동방정식을 한 줄로 압축한다는 데 있다. 임의의 물리량 $f$에 대해

$$\frac{df}{dt} = \{f,H\} + \frac{\partial f}{\partial t}$$

이 성립한다. $f=q_i$, $f=p_i$를 넣으면 곧장 정준방정식 $\dot q_i = \partial H/\partial p_i$, $\dot p_i = -\partial H/\partial q_i$가 튀어나온다. 즉 해밀토니안 $H$는 위상공간 위의 흐름을 생성하는 생성원(generator)이고, 포아송 괄호는 그 흐름을 따라가는 미분 연산이다.

여기서 보존량의 판정 기준이 즉시 나온다. 명시적 시간 의존성이 없는 $f$에 대해

$$\{f,H\}=0 \iff f \text{는 보존량}$$

에너지 보존은 반대칭성에서 공짜로 나온다($\{H,H\}=0$). 병진 대칭이면 총 운동량이, 회전 대칭이면 각운동량이 보존된다는 뇌터 정리의 고전역학 버전도 이 한 줄이다. 더 나아가 $f$와 $g$가 둘 다 보존량이면 야코비 항등식에 의해 $\{f,g\}$도 보존량이다. 보존량 두 개를 알면 세 번째를 공짜로 얻는 이 성질(포아송 정리)은 적분가능계 이론의 출발점이다.³

4. 정준변환과 심플렉틱 구조[편집]

좌표변환 $(\mathbf q,\mathbf p)\to(\mathbf Q,\mathbf P)$가 정준변환(canonical transformation)이라는 것은, 새 좌표가 기본 포아송 괄호를 그대로 만족한다는 뜻이다.

$$\{Q_i,Q_j\}=0,\quad \{P_i,P_j\}=0,\quad \{Q_i,P_j\}=\delta_{ij}$$

이 조건은 변환의 자코비안 행렬 $M$이

$$M^{\mathsf T} J M = J,\qquad J=\begin{pmatrix}0 & I\\ -I & 0\end{pmatrix}$$

를 만족한다는 것과 동치다. 이런 $M$을 심플렉틱 행렬이라 부르고, 위상공간에 얹힌 이 구조를 심플렉틱 구조라 한다. 유클리드 기하가 “길이를 보존하는 변환”(직교행렬, 회전행렬)의 기하라면, 해밀토니안 역학은 “심플렉틱 2형식 $\sum dq_i\wedge dp_i$를 보존하는 변환”의 기하다.

$\det M = 1$이 자동으로 따라오므로 위상공간 부피가 보존된다. 이것이 리우빌 정리이고, 통계역학에서 앙상블 평균이 성립하는 근거가 된다. 즉 리우빌 정리는 심플렉틱 구조의 부산물 중 가장 약한 결과다. 부피 하나만 지켜지는 게 아니라, 모든 짝수 차원 단면적(푸앵카레 불변량)이 전부 지켜진다.

5. 왜 심플렉틱 적분기는 에너지를 잃지 않는가[편집]

수치적분기는 결국 위상공간 위의 이산 사상 $(\mathbf q_n,\mathbf p_n)\to(\mathbf q_{n+1},\mathbf p_{n+1})$이다. 이 사상의 자코비안이 심플렉틱 조건 $M^{\mathsf T}JM=J$를 만족하면 심플렉틱 적분기, 아니면 그냥 적분기다. 룽게-쿠타법 4차는 국소 오차가 $O(h^5)$로 훨씬 작지만 심플렉틱하지 않고, velocity-Verlet은 2차밖에 안 되지만 심플렉틱하다.

후방 오차 해석(backward error analysis)이 이 차이를 설명한다. 심플렉틱 적분기의 수치해는 원래 해밀토니안 $H$의 근사해가 아니라, 섭동된 해밀토니안

$$\tilde H = H + h^2 H_2 + h^4 H_4 + \cdots$$

의 정확한 해에 지수적으로 가깝다. 수치해는 $\tilde H$를 정확히 보존하고, $\tilde H$는 $H$와 $O(h^2)$만큼만 떨어져 있으므로, 계산된 에너지는 그 폭 안에서 진동할 뿐 표류하지 않는다. 반면 비심플렉틱 적분기는 보존할 섭동 해밀토니안 자체가 존재하지 않아서, 스텝마다 생긴 미세한 부피 축소·팽창이 누적되어 에너지가 단조롭게 새거나 폭주한다. 분자동역학이 수억 스텝을 돌리면서도 2차 적분기를 고수하는 이유가 이것이다.⁴

주의할 점은 심플렉틱성이 스텝 크기를 고정했을 때만 유지된다는 것이다. 적응 스텝을 쓰는 순간 위 논리가 무너지므로, 근접조우가 잦은 계에서는 시간변수 변환(regularization) 같은 별도의 처방이 필요하다.

6. 양자역학으로의 다리[편집]

디랙이 발견한 대응 원리는 고전역학과 양자역학을 잇는 가장 유명한 사전(dictionary)이다.

$$[\hat A,\hat B] = i\hbar\,\widehat{\{A,B\}} + O(\hbar^2)$$

즉 교환자를 $i\hbar$로 나눈 것이 $\hbar\to0$ 극한에서 포아송 괄호가 된다. $\{q,p\}=1$이 $[\hat q,\hat p]=i\hbar$로, 시간 발전식 $\dot f=\{f,H\}$가 하이젠베르크 방정식 $\dot{\hat A} = \frac{1}{i\hbar}[\hat A,\hat H]$로 그대로 번역된다. 슈뢰딩거 방정식을 세울 때 고전 해밀토니안을 연산자로 바꿔 치는 “정준 양자화”가 정당화되는 근거가 여기다.

물론 이 사전은 완전하지 않다. $q^2p^2$ 같은 고차 곱은 연산자 순서 모호성 때문에 유일한 양자 대응물이 없고, 그로엔볼드-반 호프 정리는 모든 고전 물리량을 포아송 괄호 구조를 보존하며 연산자로 사상하는 것이 원리적으로 불가능함을 보였다. 그럼에도 계산물리와 양자화학 코드가 고전적 직관에 기대어 굴러가는 것은, 이 사전이 실용적인 범위에서는 거의 항상 맞기 때문이다.

7. 관련 문서[편집]

• 해밀토니안 역학 · 라그랑주 역학
• 리우빌 정리
• 베를레 적분 · 룽게-쿠타법
• 분자동역학 · 랑주뱅 동역학
• 양자역학 · 슈뢰딩거 방정식
• 자코비안 행렬 · 회전행렬
• 뇌터 정리 · 리 군

8. Footnotes[편집]

1. 시메옹 드니 푸아송이 1809년에 도입했다. 정작 본인은 이것이 100여 년 뒤 양자역학의 뼈대가 될 줄은 몰랐을 것이다. 수학자가 심심해서 만든 구조가 물리를 잡아먹는 흔한 전개. ↩

2. 학부 고전역학 과제로 나오면 대부분 “성립함을 확인했다”로 넘어간다. 좌표 미분을 벡터장으로 보면 $X_f$들의 리 괄호가 다시 $X_{\{f,g\}}$가 된다는 사실에서 세 줄로 끝나지만, 그건 미분기하를 배운 다음 이야기다. ↩

3. 자유도 $N$인 계에서 서로 포아송 괄호가 0인 독립 보존량을 $N$개 찾으면 그 계는 완전적분가능(Liouville-Arnold)하다. 실제로 그런 계는 조화진동자·케플러 문제·토다 격자 정도로 손에 꼽는다. 나머지 우주 대부분은 카오스라서 우리가 직접수치모사로 밥을 벌어먹는 것이다. ↩

4. “정확도 4차인 RK4가 2차인 Verlet보다 나쁠 리가 없다”고 믿고 100만 스텝을 돌린 뒤 총에너지가 30% 증발한 궤적을 보게 되는 것이, 수치해석 입문자의 국룰 통과의례다. 오차 차수는 한 스텝의 이야기이고, 보존은 구조의 이야기다. ↩