KAM 정리

계산물리 물리 마지막 수정: 2026-08-02 04:29:47

1. 개요[편집]

KAM 정리 Kolmogorov–Arnold–Moser Theorem
제안	Kolmogorov (1954) → Arnold (1963) → Moser (1962)
대상	적분가능계에 가한 작은 섭동
결론	충분히 비합리적인 회전수의 불변 토러스는 살아남는다
기술적 난점	작은 분모(small divisors)
생존 측도	1 − O(√ε)
가장 질긴 토러스	황금비 회전수

섭동 급수가 발산한다는 것과 운동이 무너진다는 것은 전혀 다른 이야기다. 100년 걸려 그 차이를 구분한 정리.

KAM 정리(Kolmogorov–Arnold–Moser theorem)는 적분가능한 해밀토니안 역학 계에 충분히 작은 섭동을 가했을 때, 충분히 비합리적인 회전수를 가진 불변 토러스가 약간 일그러진 채 그대로 살아남는다는 정리다. 살아남는 토러스가 예외적인 몇 개가 아니라 위상공간에서 양의 측도를 차지한다는 것이 핵심이며, 이 한 문장이 “일반적인 역학계는 결국 에르고딕하게 뒤섞인다”는 19세기의 믿음을 정면으로 부수었다.

배경은 태양계 안정성 문제다. 행성 하나만 있으면 케플러 문제라 적분가능하고 궤도는 영원히 닫힌다. 행성끼리의 인력을 섭동으로 넣으면 라플라스·푸앵카레 시대의 섭동 급수가 나오는데, 푸앵카레는 1890년대에 이 급수가 일반적으로 발산한다는 것을 보였다. 그래서 반세기 동안 “그러므로 운동은 불안정하다”가 기본 정서였다. 콜모고로프가 1954년 국제수학자대회 강연에서 뒤집은 것이 정확히 이 추론이다 — 급수가 발산하는 것과 운동이 무너지는 것은 별개이며, 수렴 방식을 바꾸면 살아남는 궤도를 실제로 구성할 수 있다.¹

2. 작은 분모 문제[편집]

$n$ 자유도 적분가능계를 작용-각 변수로 쓰면 해밀토니안이 작용변수만의 함수 $H_0(I)$ 이고, 운동은 $n$ 차원 토러스 위의 등속 회전 $\dot\theta = \omega(I) = \partial H_0/\partial I$ 다. 여기에 섭동을 얹는다.

H(I,\theta) = H_0(I) + \epsilon H_1(I,\theta)

섭동을 없애는 정준 변환을 푸리에 급수로 찾으면, 계수마다 분모에 $k\cdot\omega$ 가 튀어나온다.

S(I,\theta) = \sum_{k \neq 0} \frac{H_{1,k}(I)}{i\,k\cdot\omega(I)}\,e^{i k\cdot\theta}

문제는 $k\cdot\omega = 0$ 이 되는 공명이 주파수 공간에 조밀하게 깔려 있다는 점이다. 정확히 0이 아니어도 임의로 작은 $k \cdot \omega$ 가 무한히 많이 나타나고, 그러면 이 급수는 폭발한다. 이것이 작은 분모(small divisors) 문제이며, 천체역학 섭동 급수가 발산하는 이유이자 KAM 이전 100년의 벽이었다.

KAM의 해법은 두 갈래다. 첫째, 문제 있는 주파수를 아예 제외한다. 분모가 너무 작아지지 않는 주파수만 상대하겠다는 것으로, 그 조건이 디오판토스 조건이다.

\lvert k\cdot\omega \rvert \;\ge\; \gamma\,\lvert k \rvert^{-\tau} \qquad \forall\, k \in \mathbb{Z}^n \setminus \{0\}

$\tau > n-1$ 이면 이 조건을 만족하는 $\omega$ 들이 완전측도를 이룬다. 즉 공명 주파수는 조밀하지만 측도 0이고, 제외해도 잃는 게 거의 없다. 둘째, 뉴턴법식 초선형 수렴을 쓴다. 섭동 차수를 $\epsilon, \epsilon^2, \epsilon^3, \dots$ 로 한 칸씩 올리는 고전 급수 대신 $\epsilon^2, \epsilon^4, \epsilon^8, \dots$ 로 제곱해 가며 잡으면, 작은 분모가 만드는 발산보다 수렴이 빨라 이긴다. 콜모고로프의 진짜 기여가 이 이차 수렴 도식이다.

여기에 하나 더 필요한 것이 비퇴화(뒤틀림) 조건이다.

\det\!\left(\frac{\partial^2 H_0}{\partial I^2}\right) = \det\!\left(\frac{\partial \omega}{\partial I}\right) \neq 0

작용변수를 조금 움직이면 주파수도 반드시 따라 움직여야 한다는 뜻이다. 이게 있어야 “좋은 주파수”를 골라 그에 맞는 토러스를 지목할 수 있다. 조화진동자 다발처럼 모든 토러스의 주파수가 같은 완전 퇴화계에는 KAM을 그대로 못 쓴다.

결론: 위 조건들 아래 $\epsilon$ 이 충분히 작으면 디오판토스 토러스들이 살아남고, 파괴되는 영역의 측도는 $O(\sqrt\epsilon)$ , 즉 생존 측도가 $1 - O(\sqrt\epsilon)$ 이다. 제곱근이 나오는 이유는 $\gamma \sim \sqrt\epsilon$ 정도까지 조건을 느슨하게 풀 수 있기 때문이며, 공명 하나가 잡아먹는 폭이 $\sqrt\epsilon$ 이라는 진자 근사 결과와도 맞는다.

3. 황금비가 마지막까지 버티는 이유[편집]

디오판토스 조건은 “유리수로 근사하기 어려운 수”를 고르는 조건이다. 그러면 가장 근사하기 어려운 수는? 연분수 전개가 답을 준다. 큰 계수가 나오는 곳에서 유리수 근사가 급격히 좋아지므로, 계수가 전부 1인 수가 가장 나쁘게 근사된다. 그게 황금비다.

\gamma_{\text{gold}} = \frac{\sqrt5 - 1}{2} = [0;1,1,1,\dots] = 0.6180\ldots

그래서 섭동을 키울 때 황금비 회전수를 가진 토러스가 가장 늦게 부서진다. 표준 사상에서 이 사실이 정량화됐다. 넓이를 보존하는 2차원 사상에서 위상공간을 가로지르는 마지막 불변곡선이 정확히 황금 회전수 곡선이며, 그것이 깨지는 값이

K_c = 0.971635406\ldots

다. 이 숫자를 뽑은 것이 그린의 잔여 판정(1979)이다. 황금비의 연분수 근사인 피보나치 비 $F_{n-1}/F_n$ 을 회전수로 갖는 주기궤도들을 찾아, 각각의 잔여

R = \frac{2 - \operatorname{tr} M}{4}

를 계산한다( $M$ 은 그 주기궤도의 단일주기 야코비안). 궤도 차수를 올릴 때 $R \to 0$ 이면 근처 토러스가 살아 있고, $R \to \infty$ 면 이미 부서졌으며, 정확히 임계에서만 $R$ 이 유한한 값( $\approx 0.25$ )으로 수렴한다. 토러스를 직접 그리는 대신 주기궤도의 안정성만 보고 판정한다는 발상이라 수치적으로 대단히 효율적이고, 지금도 사실상 표준이다.

4. 부서진 다음에 남는 것들[편집]

KAM은 살아남는 쪽 이야기고, 부서지는 쪽에는 별도의 구조가 있다.

푸앵카레-비르코프 정리. 회전수가 유리수 $p/q$ 인 공명 토러스는 섭동에 가장 취약해서 곧바로 깨지지만, 완전히 사라지지 않고 짝수 개( $2kq$ )의 주기점으로 부서진다. 절반은 타원점(안정), 절반은 쌍곡점(불안정)이고, 타원점 주위에 다시 섬이 생긴다. 이것이 푸앵카레 단면에서 보이는 섬 사슬이다. 각 섬 안을 확대하면 그 안에 또 섬 사슬이 있고, 자기유사 구조가 끝없이 이어진다.
캔토러스(cantorus). 임계값을 막 넘긴 토러스는 사라지는 게 아니라 칸토어 집합 꼴의 잔해로 남는다. 위상공간을 완전히 막지는 못하지만 구멍이 작아서, 궤도가 통과하는 데 걸리는 시간이 어마어마하게 길어진다. 카오스 궤도가 섬 경계에 오래 들러붙는 끈적임(stickiness) 현상의 정체이며, 통과 유량은 오브리-마더 이론의 변분 원리로 계산된다.
공명 겹침. 이웃한 두 공명의 폭이 서로 닿으면 그 사이 토러스가 남아날 수 없다. 치리코프의 이 대략적 판정은 $K_c$ 를 2.5배쯤 과대평가하지만, 손으로 임계값을 어림하는 데는 여전히 유용하다.

5. 자유도 2와 3 사이의 절벽[편집]

KAM 정리에서 가장 실용적인 따름정리이자 자주 오해받는 부분이다. 문제는 차원 계산이다. $n$ 자유도 계의 위상공간은 $2n$ 차원, 에너지 초곡면은 $2n-1$ 차원, KAM 토러스는 $n$ 차원이다.

$n = 2$ : 에너지 초곡면이 3차원인데 토러스가 2차원이다. 여분 차원이 딱 1이라 토러스가 초곡면을 두 조각으로 가른다. 카오스 궤도는 위아래 토러스 사이에 영원히 갇히고, 작용변수가 아무리 시간이 흘러도 유계로 머문다. 즉 2자유도 계는 카오스가 있어도 전역 수송이 원천 봉쇄된다.
$n \ge 3$ : 초곡면이 5차원 이상인데 토러스는 3차원이다. 여분 차원이 2 이상이면 초곡면을 가르지 못한다 — 3차원 공간에 실 한 가닥을 놔도 공간이 나뉘지 않는 것과 같은 이치다. 그래서 카오스 영역이 공명망을 따라 전부 연결되고, 궤도가 토러스 사이를 비집고 아주 느리게 흘러 다닐 수 있다. 이것이 아널드 확산(1964)이다.

아널드 확산은 실재하지만 끔찍하게 느리다. 넥호로셰프의 정리가 그 속도에 상한을 주는데, 작용변수의 변화가 $\epsilon^b$ 를 넘는 데 걸리는 시간이 $\exp(\epsilon^{-a})$ 규모로 지수적이다. “언젠가 무너지지만 우주 나이 안에는 아니다”라는 식의 결론이 나오는 게 이 동네다.

6. 수치 계산과의 연결[편집]

심플렉틱 적분기가 장시간 적분에서 왜 유리한가가 KAM의 가장 실무적인 응용이다. 립프로그 같은 심플렉틱 도식으로 적분한 궤적은 원래 $H$ 를 정확히 보존하지 않지만, $H$ 에서 $O(\Delta t^p)$ 만큼 떨어진 섀도 해밀토니안 $\tilde H$ 를 (지수적으로 작은 오차를 빼면) 정확히 보존한다. 그러면 시간 이산화 자체가 “작은 섭동”이 되고, 수치 궤적은 $\tilde H$ 라는 근처 계의 KAM 토러스 위에서 논다. 토러스는 유계이므로 에너지 오차가 표류하지 않고 진동만 한다. 반대로 룽게-쿠타법처럼 심플렉틱하지 않은 도식은 보존할 섀도 해밀토니안이 없어 에너지가 단조 표류하고, 궤도가 결국 토러스를 가로질러 나간다. 계산 정밀도의 문제가 아니라 구조의 문제라 스텝을 줄여도 표류의 방향만 늦춰질 뿐 없어지지 않는다.²

같은 논리가 분자동역학에서 베를레 적분이 40년째 표준인 이유고, 해밀토니안 몬테카를로의 수락률이 궤적 길이에 따라 붕괴하지 않는 이유이기도 하다. 다만 정직하게 말하면, 자유도가 수만인 MD 계는 근적분가능계와 거리가 멀어 KAM 토러스가 차지하는 측도가 사실상 0이다. 여기서 살아남는 것은 토러스가 아니라 섀도 해밀토니안의 유계성이라는 약한 형태다.

수치적으로 KAM 구조를 직접 확인하는 표준 도구는 푸앵카레 단면이다. 닫힌 곡선이면 살아 있는 토러스, 섬 사슬이면 공명, 흩뿌려진 점이면 카오스 바다다. 에농-하일레스 계는 에너지 하나만 올려도 이 세 가지가 순서대로 나타나는 교과서 예제이며, 랴푸노프 지수를 함께 재면 토러스 위에서 0, 카오스 바다에서 양수라는 정량 판정이 붙는다.

7. 에르고딕성과 태양계[편집]

KAM은 에르고딕 가설을 정확히 어디서 부정하는가? 에르고딕성의 정의는 “불변집합의 측도가 0 아니면 1”이다. 그런데 KAM 토러스들의 합집합은 불변집합이면서 측도가 $1 - O(\sqrt\epsilon)$ 이므로 0도 1도 아니다. 따라서 근적분가능한 해밀토니안 계는 에너지 초곡면 위에서 에르고딕하지 않다. 이것이 반례가 있다는 정도가 아니라 양의 측도를 갖는 구조적 반례라는 점이 결정타였고, 통계역학의 정당화를 에르고딕성이 아닌 다른 근거(자유도 수, 혼합 시간, 전형성)에서 찾게 만든 계기가 됐다. 반대로 자유도가 커질수록 KAM 토러스의 측도는 급격히 줄어들어, 열역학적 극한에서는 실질적으로 에르고딕한 것처럼 행동한다.

태양계 이야기는 조심해서 해야 한다. KAM이 태양계의 안정성을 증명한 것이 아니다. 정리가 요구하는 $\epsilon$ 상한은 실제 행성 질량비( $10^{-3}$ 규모)보다 몇 자릿수 작고, 케플러 문제 자체가 퇴화계라 표준 비퇴화 조건도 안 맞는다. 실제로 수치 적분에 기반한 장기 계산은 내행성계가 카오스적이며 랴푸노프 시간이 500만 년 규모라고 말하고, 50억 년 규모에서 수성의 궤도가 불안정해질 확률이 1% 정도라는 결과도 나와 있다. KAM이 준 것은 “안정하다”는 보증이 아니라 **“카오스가 있어도 규칙적 궤도가 양의 측도로 공존한다”**는 구조적 그림이다.³⁴

8. 관련 문서[편집]

콜모고로프의 1954년 발표는 4쪽짜리였고 증명이 사실상 스케치였다. 그래서 한동안 “정말 증명된 게 맞느냐”는 의심이 있었고, 아널드가 해석적 경우를, 모저가 매끄러운 경우를(원래 $C^{333}$ 이라는 전설적으로 후한 미분가능도 요구로) 각각 채워 넣으면서 세 사람 이름이 붙었다. 요즘 요구되는 미분가능도는 $C^{2n+\delta}$ 수준까지 내려왔다. ↩
그래서 태양계 장기 적분 코드에 RK4를 쓰면 수억 년 뒤 행성이 조용히 나선을 그리며 태양으로 떨어지거나 날아가 버린다. 물리가 아니라 적분기가 만든 결과다. 리뷰어가 “적분기 뭐 썼냐”부터 묻는 데는 이유가 있다. ↩
라스카가 1989년부터 돌린 태양계 장기 적분은 계산 시간 자체가 논문 단위였다. 요즘은 GPU로 며칠이면 되지만, 결론은 그때나 지금이나 “내행성계는 카오스적이고, 예측 가능 지평은 수천만 년”이다. 즉 5천만 년 뒤 지구의 궤도 위상은 원리적으로 계산 불가능하다. 다행히 궤도 반지름은 그렇지 않다. ↩
KAM 정리 이름의 순서가 K-A-M인 건 발표 연도(1954-1963-1962) 순도 아니고 알파벳 순도 아니다. 그냥 콜모고로프가 먼저고 나머지 둘은 어감이 좋은 쪽으로 붙었다는 게 정설. MAK나 KMA였으면 아무도 안 외웠을 것이다. ↩