심플렉틱 적분기

수치해석 계산물리 마지막 수정: 2026-07-31 04:12:07

1. 개요[편집]

심플렉틱 적분기 Symplectic Integrator
대상	해밀토니안계 $\dot z = J\nabla H(z)$
정의 성질	한 스텝 사상의 야코비안이 $M^\top J M = J$
대표 기법	리프프로그 / 속도 베를레, 음함수 중점법, 가우스-르장드르 RK
핵심 이득	에너지 오차가 표류하지 않고 유계 진동
대표 약점	가변 시간간격과 상극

정확한 적분기는 한 걸음을 잘 걷는다. 심플렉틱 적분기는 백만 걸음 뒤에도 서 있다.

심플렉틱 적분기(symplectic integrator)는 해밀토니안 역학계를 시간 적분할 때, 한 스텝 사상 $\Phi_h:\ z_n \mapsto z_{n+1}$ 이 위상공간의 심플렉틱 구조를 정확히 보존하도록 설계된 적분기다. 여기서 “정확히”는 오차항 없이 컴퓨터의 부동소수점 연산 정밀도까지라는 뜻이다. 해를 근사하는 대신 구조를 보존하겠다는 이 발상이 기하 수치적분이라는 분야 전체의 출발점이다.

왜 이게 그렇게 대단하냐면, 궤도 자체보다 궤도가 사는 무대가 더 오래 살아남기 때문이다. 룽게-쿠타법(RK4)은 한 스텝 오차가 $O(h^5)$ 로 훨씬 작지만, 100만 스텝쯤 돌리면 에너지가 한 방향으로 슬금슬금 흘러 행성이 항성으로 빨려 들어가거나 분자동역학 계가 혼자 뜨거워진다. 반면 2차 정확도밖에 안 되는 리프프로그는 같은 시간 동안 에너지가 작은 폭으로 진동만 하고 어디로도 가지 않는다. 이 문서는 그 “왜”를 다룬다 — 구현 세부는 베를레 적분, 라그랑지안 쪽에서의 유도는 변분 적분기에 위임한다.

2. 심플렉틱 조건이란 무엇인가[편집]

정준좌표 $z = (q, p) \in \mathbb{R}^{2n}$ 와 표준 심플렉틱 행렬

J = \begin{pmatrix} 0 & I_n \\ -I_n & 0 \end{pmatrix}

을 두면 해밀턴 방정식은 $\dot z = J\,\nabla H(z)$ 한 줄로 쓰인다. 이때 정확한 흐름 $\varphi_t$ 의 야코비안 $M = \partial\varphi_t/\partial z$ 는 모든 $t$ 에서

M^\top J M = J

를 만족한다. 이 조건을 만족하는 사상을 심플렉틱 사상이라 부르고, 수치 사상 $\Phi_h$ 의 야코비안이 같은 조건을 만족하면 그 적분기가 심플렉틱이다.

기하학적으로 이건 위상공간의 2형식 $\omega = \sum_i \mathrm{d}q_i \wedge \mathrm{d}p_i$ — 즉 모든 $(q_i, p_i)$ 평면에 사영한 넓이의 합 — 이 보존된다는 뜻이다. 양변의 행렬식을 취하면 $\det M = 1$ 이 따라 나오므로 위상공간 부피도 자동 보존되고, 이것이 리우빌 정리의 이산판이다. 다만 부피 보존은 심플렉틱성의 결과이지 동치가 아니다 — 자유도가 1개( $2n=2$ )일 때만 둘이 같고, $n \ge 2$ 부터는 심플렉틱 쪽이 엄격하게 더 강한 조건이다.¹ 실제로 심플렉틱성은 포아송 괄호 구조를 통째로 보존하는 것과 같은 말이다.

3. 섀도 해밀토니안 — 에너지가 표류하지 않는 진짜 이유[편집]

심플렉틱 적분기가 에너지를 보존한다는 말은 틀렸다. 정확히는 이렇다. 후진 오차 해석(backward error analysis)에 따르면, 차수 $p$ 인 심플렉틱 적분기의 이산 궤적은 원래 $H$ 의 궤적이 아니라 수정된 해밀토니안

\tilde H(z) = H(z) + h^p H_{p+1}(z) + h^{p+1} H_{p+2}(z) + \cdots

의 정확한 궤적 위에 (거의) 놓인다. 이 $\tilde H$ 를 섀도 해밀토니안(shadow Hamiltonian)이라 부른다. 핵심은 이 급수가 여전히 해밀토니안 형태라는 것 — 심플렉틱하지 않은 적분기에서는 수정 방정식이 해밀토니안 구조를 잃어버려 소산·발산 항이 붙는다.

결과는 극적이다. 계는 $\tilde H$ 를 정확히 보존하고 $\tilde H = H + O(h^p)$ 이므로, 실제 $H$ 의 오차는 $O(h^p)$ 크기 안에서 진동만 하고 시간에 비례해 커지지 않는다. 급수 자체는 발산하지만 최적 절단에서 잔차가 $e^{-c/h}$ 로 지수적으로 작아, 이 유계성이 $t \sim e^{c/h}$ 정도의 지수적으로 긴 시간 동안 유지된다는 정량적 결과까지 있다.² 비심플렉틱 적분기의 에너지 오차가 $t$ 에 선형으로 자라는 것(장기 표류, secular drift)과 정확히 대비되는 지점이다.

한 가지 더. 고정 스텝 심플렉틱 적분기가 $H$ 까지 정확히 보존하는 일은 (적분 가능한 예외를 빼면) 일어나지 않는다. 게-마스든(Ge–Marsden) 정리가 그걸 못 박았다.³ 즉 우리가 얻는 건 “에너지 보존”이 아니라 “에너지 오차의 유계성”이며, 대부분의 응용에서 후자면 충분하다.

4. 스플리팅과 리프프로그[편집]

실용적인 심플렉틱 적분기는 거의 전부 스플리팅에서 나온다. 해밀토니안이 가분(separable) 형태

H(q,p) = T(p) + V(q)

이면 $T$ 만 있는 계(“드리프트”: $q \mathrel{+}= h\,\partial T/\partial p$ , $p$ 불변)와 $V$ 만 있는 계(“킥”: $p \mathrel{-}= h\,\nabla V$ , $q$ 불변)를 각각 정확히 풀 수 있다. 둘 다 정확한 해밀토니안 흐름이므로 심플렉틱이고, 심플렉틱 사상의 합성은 심플렉틱이다. 그래서 아무렇게나 이어 붙여도 심플렉틱성은 공짜다.

1차: $\Phi_h = \varphi^V_h \circ \varphi^T_h$ — 심플렉틱 오일러법.
2차: $\Phi_h = \varphi^V_{h/2} \circ \varphi^T_h \circ \varphi^V_{h/2}$ — 스트랑 분할. 이게 바로 리프프로그 = 스퇴르머-베를레 = 속도 베를레이며, 셋은 같은 사상의 다른 표기다. 대칭(시간 가역)이라 차수가 자동으로 2차가 된다.

가분이 아닌 $H$ (예: 상대론적 항, 회전 강체)에서는 스플리팅이 안 먹히므로 음함수 계열을 쓴다. 음함수 중점법이 2차 심플렉틱이고, 그 일반화인 가우스-르장드르 룽게-쿠타는 $s$ 단계에서 $2s$ 차 심플렉틱이다. 대신 매 스텝 비선형 연립방정식을 풀어야 해서 비싸다.

5. 고차화 — 요시다 합성[편집]

2차 대칭 적분기 $S_2(h)$ 를 세 번 합성해 4차를 만드는 것이 요시다 합성(Yoshida, 1990)이다.

S_4(h) = S_2(w_1 h)\,S_2(w_0 h)\,S_2(w_1 h), \qquad w_1 = \frac{1}{2 - 2^{1/3}},\quad w_0 = \frac{-2^{1/3}}{2 - 2^{1/3}}

수치로는 $w_1 \approx 1.3512$ , $w_0 \approx -1.7024$ 이고 $w_0 + 2w_1 = 1$ 이다. 눈에 띄는 건 가운데 계수가 음수라는 것 — 시간을 거꾸로 감는 스텝이 반드시 끼어든다. 이건 요시다가 게을러서가 아니라, 2차를 넘는 스플리팅에서 모든 계수를 양수로 두는 것이 원리적으로 불가능하다는 정리 때문이다.⁴ 같은 방식으로 6차·8차 계수도 만들어져 있고, 천체역학 코드에서는 4차·6차가 흔하다.

6. 가변 시간간격이라는 함정[편집]

여기가 배신 포인트다. 오차 추정으로 스텝을 조절하는 평범한 적응 제어를 심플렉틱 적분기에 얹으면 장기 안정성이 그대로 날아간다. 이유는 명확하다 — 섀도 해밀토니안 논증은 “매 스텝 같은 $h$ 로 같은 $\tilde H$ 를 보존한다”는 전제 위에 서 있는데, $h$ 가 매번 바뀌면 보존하는 대상도 매번 바뀌어 에너지가 랜덤워크하듯 표류한다. 실제로 적응 리프프로그가 고정 스텝 리프프로그보다 장기 궤도에서 더 나빠지는 것은 잘 알려진 현상이다.

우회로는 있다. 시간 자체를 정준변환의 대상으로 끌어들이는 푸앵카레 시간 변환(확장 위상공간에서 새 해밀토니안을 정의해 균일 스텝으로 적분)이나, 스텝 크기를 대칭적으로 결정해 시간 가역성을 지키는 방식이 그것이다. 근접 조우가 잦은 계에서는 아예 스플리팅을 다시 설계한다.

7. 어디서 쓰이나[편집]

천체역학 장시간 적분. 태양계의 수십억 년 안정성 계산은 케플러 항과 상호작용 항으로 쪼갠 위즈덤-홈맨(Wisdom–Holman) 사상이 표준이다. 여기서 궁금한 건 “10억 년 뒤 지구의 좌표”가 아니라 궤도 요소의 통계적 거동이라, 심플렉틱성이 곧 신뢰도다. 이런 계의 카오스 이론적 성질과 푸앵카레 단면 구조도 심플렉틱 적분기 없이는 제대로 안 보인다.
분자동역학. 속도 베를레가 사실상 유일한 국룰. 다만 서모스탯이나 랑주뱅 동역학을 붙이는 순간 계가 해밀토니안이 아니게 되므로, 심플렉틱성 논증은 무열욕조(NVE) 구간에만 엄밀히 적용된다.
마르코프 연쇄 몬테카를로의 해밀토니안 몬테카를로(HMC). 여기선 심플렉틱성이 통계적 정확성의 필수 조건이다 — 메트로폴리스 수락비가 옳으려면 제안 사상이 시간 가역이고 부피를 보존해야 하고, 리프프로그가 정확히 그 둘을 만족한다. 게다가 에너지 오차가 유계라 수락률이 스텝 수에 따라 무너지지 않는다.
가속기 빔 동역학, 플라즈마 입자 궤도, 게임 물리. 마지막 항목은 “왜 천이 안 터지는가”의 답이기도 하다.

8. 관련 문서[편집]

$2n \times 2n$ 심플렉틱 조건은 반대칭성 때문에 $n(2n-1)$ 개의 독립 방정식을 낳는다. $n=1$ 이면 1개( $\det M = 1$ )라 부피 보존과 같지만, $n=2$ 면 6개라 $\det M = 1$ 하나로는 어림도 없다. 그로모프의 심플렉틱 낙타 정리(1985)가 이 차이를 극적으로 보여준다 — 부피만 보존하면 낙타를 바늘구멍에 통과시킬 수 있지만, 심플렉틱 사상은 그걸 못 한다. ↩
Benettin & Giorgilli (1994), Hairer & Lubich (1997). 해석적 $H$ 와 충분히 작은 $h$ 라는 전제가 붙는다. “지수적으로 긴 시간”이라 해도 무한은 아니고, 별의 근접 조우처럼 $\nabla V$ 가 폭발하는 사건이 끼면 논증의 전제가 그 자리에서 깨진다. 공짜 점심은 없다. ↩
Ge Zhong & Marsden (1988). 대략 “고정 스텝 적분기가 심플렉틱이면서 $H$ 를 정확히 보존하면, 그것은 시간 재매개변수화를 제외하고 정확한 흐름 그 자체다”라는 내용. 즉 “심플렉틱 + 에너지 정확 보존”을 동시에 요구하는 건 미분방정식을 이미 풀어놨다는 뜻이라 수치해석적으로 무의미하다. ↩
Suzuki (1991), Goldman & Kaper (1996). 3차 이상에서는 계수 중 최소 하나가 음수여야 한다. 그래서 확산 방정식처럼 시간을 거꾸로 돌리면 발산하는 계에는 3차 이상 스플리팅을 그대로 못 쓴다. 복소 계수를 허용하는 우회로가 연구되고 있는데, 이 시점에서 이미 “간단해서 좋았던” 스플리팅의 미덕은 사라진 뒤다. ↩