심플렉틱 적분기

오류 제보
수치해석 계산물리 마지막 수정: 2026-09-29 04:26:55
○ 공개 검수 대기 AI 보조 작성 표시 기준 보기

1. 개요

심플렉틱 적분기
Symplectic Integrator
대상해밀토니안계 ż = J∇H(z)
정의 성질한 스텝 사상의 야코비안이 MTJM = J
대표 기법리프프로그 / 속도 베를레, 음함수 중점법, 가우스-르장드르 RK
핵심 이득에너지 오차가 표류하지 않고 유계 진동
대표 약점가변 시간간격과 상극
오해 주의단기 정확도는 같은 차수의 비심플렉틱 방법보다 나쁠 수 있다

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

심플렉틱 적분기(symplectic integrator)는 해밀토니안 역학계를 시간 적분할 때, 한 스텝 사상 Φh: zn↦zn+1\Phi_h:\ z_n \mapsto z_{n+1} 이 위상공간의 심플렉틱 구조를 정확히 보존하도록 설계된 적분기다. 여기서 “정확히”는 오차항 없이 컴퓨터의 부동소수점 연산 정밀도까지라는 뜻이다. 해를 근사하는 대신 구조를 보존하겠다는 이 발상이 기하 수치적분이라는 분야 전체의 출발점이고, 가속기 빔 동역학에서 나온 루스(Ruth, 1983)의 정준 적분 기법이 그 첫 공개 판본으로 꼽힌다.

왜 이게 그렇게 대단하냐면, 궤도 자체보다 궤도가 사는 무대가 더 오래 살아남기 때문이다. 룽게-쿠타법(RK4)은 한 스텝 오차가 O(h5)O(h^5) 로 훨씬 작지만 에너지 오차가 시간에 비례해 한 방향으로 자란다. 반면 2차 정확도밖에 안 되는 리프프로그는 같은 시간 동안 에너지가 작은 폭으로 진동만 하고 어디로도 가지 않는다. 이 문서는 그 “왜”를 다룬다 — 구현 세부는 베를레 적분, 라그랑지안 쪽에서의 유도는 변분 적분기에 위임한다.

다만 흔한 과장은 걷어내자. “심플렉틱이니까 항상 더 정확하다”는 틀렸다. 아래 검산에서 같은 스텝의 RK4가 1592 주기 뒤에도 리프프로그보다 에너지 오차가 440배 작았다. 심플렉틱의 진짜 값어치는 오차의 크기가 아니라 시간 의존성(유계 진동 대 선형 표류)에 있다.

2. 심플렉틱 조건이란 무엇인가

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

J=(0In−In0)J = \begin{pmatrix} 0 & I_n \\ -I_n & 0 \end{pmatrix}

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

M⊤JM=JM^\top J M = J

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

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

2.1. 검산: RK4 는 스텝마다 위상공간을 정확히 얼마나 줄이는가

자유도 1(2n=22n=2)에서는 심플렉틱성이 det⁡M=1\det M = 1 과 동치이므로, 선형계에서 한 스텝 사상의 행렬식을 재면 판정이 끝난다. 조화진동자 H=(p2+q2)/2H = (p^2+q^2)/2 에 대해 두 적분기의 한 스텝 행렬을 기저벡터 두 개를 통과시켜 뽑고 행렬식을 계산했다.

hh리프프로그 det⁡M−1\det M - 1RK4 det⁡M−1\det M - 1이론값 −h6/72+h8/576-h^6/72 + h^8/576
0.50 (정확히)-2.102322e-4-2.102322e-4
0.10 (정확히)-1.387153e-8-1.387153e-8
0.050 (정확히)-2.169460e-10-2.169461e-10

리프프로그는 부동소수점 0이 아니라 비트 단위로 정확히 1을 낸다(세 개의 전단 변환 곱이라 행렬식이 구조적으로 1이다). RK4는 스텝마다 위상공간 부피를 1−h6/721 - h^6/72 배로 수축시키고, 실측이 이 해석식과 일곱 자리까지 일치한다. 이 수축이 뒤에 나오는 에너지 표류의 방향까지 결정한다 — 아래 케플러 검산에서 RK4의 ΔH\Delta H 가 계속 음수인 것이 우연이 아니다.

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

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

H~(z)=H(z)+hpHp+1(z)+hp+1Hp+2(z)+⋯\tilde H(z) = H(z) + h^p H_{p+1}(z) + h^{p+1} H_{p+2}(z) + \cdots

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

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

선형계에서는 이 H~\tilde H 를 손으로 적을 수 있다. 조화진동자에 대한 리프프로그는

H~∝(1−h24)q2+p2\tilde H \propto \left(1 - \frac{h^2}{4}\right) q^2 + p^2

라는 2차형식을 정확히 보존한다. h=0.5h = 0.5 와 h=0.2h = 0.2 로 각각 20만 스텝을 돌려 이 양의 상대 변화를 재보면 6.0×10−146.0\times10^{-14} · 1.2×10−131.2\times10^{-13} 으로, 이산화 오차가 아니라 반올림만 남는다. 진짜 HH 가 아니라 hh 에 의존하는 다른 양이 보존된다는 사실 자체가 심플렉틱 적분기의 정체를 요약한다.

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

4. 스플리팅과 리프프로그

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

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

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

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

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

4.1. 검산: 케플러 궤도 1592 주기 — 유계 진동 대 선형 표류

이심률 e=0.5e = 0.5 케플러 문제(H0=−0.5H_0 = -0.5, 주기 2π2\pi)를 h=0.005h = 0.005 로 200만 스텝(t=104t = 10^4, 약 1592 주기) 적분하고 매 스텝 ΔH=H−H0\Delta H = H - H_0 를 기록했다.

적분기ΔH\Delta H 범위최종 ΔH\Delta H시간 의존성
리프프로그[0, +3.397×10−5][0,\ +3.397\times10^{-5}]+3.336e-5없음 (진폭 고정)
RK4[−7.578×10−8, 0][-7.578\times10^{-8},\ 0]-7.578e-8tt 에 선형

RK4의 표류가 선형이라는 것은 스텝 수를 두 배로 할 때마다 최종 오차가 두 배가 되는 것으로 확인된다: 25만 → -9.711e-9, 50만 → -1.919e-8, 100만 → -3.806e-8, 200만 → -7.578e-8. 리프프로그는 같은 자리에서 25만 스텝과 200만 스텝의 진폭이 둘 다 3.397e-5 로 같다.

그런데 이 표의 절대값을 보라. 1592 주기를 돈 뒤에도 RK4의 에너지 오차가 리프프로그보다 440배 작다. 흔히 인용되는 “RK4 로 장기 적분하면 행성이 별로 빨려 들어간다”는 이 파라미터에서 성립하지 않는다. 옳은 질문은 “언제 역전되는가”다. RK4 표류 기울기 7.578×10−127.578\times10^{-12}/단위시간으로 리프프로그 진폭에 도달하는 시각을 풀면 t≈4.5×106t \approx 4.5\times10^6, 약 71만 주기다.

스텝당 힘 평가 횟수(리프프로그 1회, RK4 4회)를 맞춘 등비용 비교가 더 공정하다.4 리프프로그를 h=0.00125h = 0.00125 로 800만 스텝 돌리면 진폭이 2.123×10−62.123\times10^{-6} 으로 줄고(h2h^2 스케일링을 정확히 따른다 — hh 를 4배 키운 경우 대비 실측 16.0배), 역전 시각이 t≈2.8×105t \approx 2.8\times10^5, 약 4.5만 주기로 내려온다. 태양계 수십억 년 적분이나 해밀토니안 몬테카를로의 수백만 궤적처럼 그 시간을 실제로 넘기는 계산에서만 심플렉틱성이 결정적이라는 뜻이다. 100 주기 쏘고 끝나는 계산에 심플렉틱을 들이대며 정확도를 자랑하는 것은 근거가 없다.

5. 고차화 — 요시다 합성

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

S4(h)=S2(w1h) S2(w0h) S2(w1h),w1=12−21/3,w0=−21/32−21/3S_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}}

수치로는 w1=1.3512071920w_1 = 1.3512071920, w0=−1.7024143839w_0 = -1.7024143839 이고 w0+2w1=1w_0 + 2w_1 = 1 이 배정밀도에서 정확히 성립한다(실측 1.0000000000000001.000000000000000). 같은 4차 계수가 포레스트-루스(Forest–Ruth)와 네리(Neri)에게서 독립적으로 나왔다는 점을 요시다 본인이 논문에서 밝히고 있다.

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

차수가 실제로 오르는지는 케플러 궤도를 정확히 한 주기 적분한 뒤 출발점과의 거리로 확인했다. 리프프로그는 hh 를 절반으로 줄일 때 오차가 4.00배씩(1.80e-2 → 4.49e-3 → 1.12e-3 → 2.80e-4), 요시다 4차는 16.0배씩(5.48e-5 → 3.43e-6 → 2.14e-7 → 1.34e-8) 줄어 관측 차수가 각각 2.00 · 4.00 으로 떨어진다. 스텝당 비용이 3배인 대신 오차는 h4h^4 로 줄어드니, 정확도를 조금만 높여도 4차가 이긴다.

6. 가변 시간간격이라는 함정

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

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

7. 어디서 쓰이나

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

8. 관련 문서

9. Footnotes

  1. 2n×2n2n \times 2n 심플렉틱 조건은 반대칭성 때문에 n(2n−1)n(2n-1) 개의 독립 방정식을 낳는다. n=1n=1 이면 1개(det⁡M=1\det M = 1)라 부피 보존과 같지만, n=2n=2 면 6개라 det⁡M=1\det M = 1 하나로는 어림도 없다. 그로모프의 심플렉틱 낙타 정리(1985)가 이 차이를 극적으로 보여준다 — 부피만 보존하면 낙타를 바늘구멍에 통과시킬 수 있지만, 심플렉틱 사상은 그걸 못 한다. ↩

  2. Benettin & Giorgilli (1994), Hairer & Lubich (1997). 해석적 HH 와 충분히 작은 hh 라는 전제가 붙는다. “지수적으로 긴 시간”이라 해도 무한은 아니고, 별의 근접 조우처럼 ∇V\nabla V 가 폭발하는 사건이 끼면 논증의 전제가 그 자리에서 깨진다. 공짜 점심은 없다. ↩

  3. Ge Zhong & Marsden (1988). 대략 “고정 스텝 적분기가 심플렉틱이면서 HH 를 정확히 보존하면, 그것은 시간 재매개변수화를 제외하고 정확한 흐름 그 자체다”라는 내용. 즉 “심플렉틱 + 에너지 정확 보존”을 동시에 요구하는 건 미분방정식을 이미 풀어놨다는 뜻이라 수치해석적으로 무의미하다. ↩

  4. 그래서 벤치마크 그림을 볼 때는 가로축이 “스텝 수”인지 “힘 평가 횟수”인지를 먼저 봐야 한다. 스텝 수를 가로축으로 놓으면 스텝당 4번 힘을 계산하는 RK4가 공짜로 유리해지고, 반대로 심플렉틱 진영이 hh 를 4배 작게 잡아 비교하면 그쪽이 유리해진다. 위 등비용 표는 후자 쪽으로 맞춘 것이다. ↩

  5. Suzuki (1991)의 프랙탈 경로적분 이론이 근거를 준다. 3차 이상에서는 계수 중 최소 하나가 음수여야 한다. 그래서 확산 방정식처럼 시간을 거꾸로 돌리면 발산하는 계에는 3차 이상 스플리팅을 그대로 못 쓴다. 복소 계수를 허용하는 우회로가 연구되고 있는데, 이 시점에서 이미 “간단해서 좋았던” 스플리팅의 미덕은 사라진 뒤다. ↩

검증 자료

이 문서의 주요 출처

  1. E. Hairer, C. Lubich & G. Wanner, Geometric Numerical Integration — Structure-Preserving Algorithms for ODEs, 2nd ed. (Springer, 2006)
  2. R. D. Ruth, "A canonical integration technique", IEEE Trans. Nucl. Sci. NS-30 (1983) 2669–2671
  3. H. Yoshida, "Construction of higher order symplectic integrators", Phys. Lett. A 150 (1990) 262–268
  4. Ge Zhong & J. E. Marsden, "Lie-Poisson Hamilton-Jacobi theory and Lie-Poisson integrators", Phys. Lett. A 133 (1988) 134–139
  5. G. Benettin & A. Giorgilli, "On the Hamiltonian interpolation of near-to-the-identity symplectic mappings…", J. Stat. Phys. 74 (1994) 1117–1143
  6. E. Hairer & C. Lubich, "The life-span of backward error analysis for numerical integrators", Numer. Math. 76 (1997) 441–462
  7. M. Suzuki, "General theory of fractal path integrals with applications to many-body theories and statistical physics", J. Math. Phys. 32 (1991) 400–407
  8. J. Wisdom & M. Holman, "Symplectic maps for the N-body problem", Astron. J. 102 (1991) 1528–1538
  9. S. Duane, A. D. Kennedy, B. J. Pendleton & D. Roweth, "Hybrid Monte Carlo", Phys. Lett. B 195 (1987) 216–222
  10. M. Gromov, "Pseudo holomorphic curves in symplectic manifolds", Invent. Math. 82 (1985) 307–347
이 문서가 문제를 푸는 데 도움이 됐나요? 한 번의 클릭이 다음 검수 순서를 정합니다.

의견을 반영할게요. 감사합니다.

다음으로 읽기 해밀토니안 역학 약 5분 더 탐색 →