적분가능계

계산물리 물리 마지막 수정: 2026-08-03 04:26:09

1. 개요[편집]

적분가능계 Integrable System
정의(리우빌)	$n$자유도에 대합(involution)인 독립 보존량이 $n$개
결과	운동이 $n$차원 원환면 위의 준주기 운동으로 환원
표준 도구	작용-각 변수, 락스 쌍, 역산란 변환
대표 예	케플러 문제, 토다 사슬, KdV, 사인-고든
현실	지극히 예외적 — 일반적인 계는 비적분

풀리는 문제는 왜 풀리는가. 답은 “우연히 똑똑해서”가 아니라 “숨은 대칭이 그만큼 많아서”다.

적분가능계(integrable system)는 $n$ 자유도 해밀토니안 역학계 중에서, 서로 포아송 괄호가 0이고 함수적으로 독립인 보존량이 자유도 수만큼( $n$ 개) 존재해 해가 유한 번의 구적(quadrature)으로 표현되는 계를 말한다.¹ 케플러 문제, 강체의 자유회전, 조화진동자 사슬, 토다 사슬, 그리고 KdV 방정식 같은 몇몇 비선형 편미분방정식이 여기 속한다.

시뮬레이션 하는 사람 입장에서 적분가능계가 중요한 이유는 두 가지다. 첫째, 정확해가 있으니 수치기법의 검증 시험대로 최상급이다 — 심플렉틱 적분기의 장기 에너지 거동을 논할 때 나오는 그림의 절반은 케플러 문제다. 둘째, 적분가능계는 예외 중의 예외라서, 거기서 살짝 벗어났을 때 무슨 일이 벌어지는가(KAM 정리 → 카오스 이론)가 곧 현대 동역학의 본론이다. 적분가능계는 목적지가 아니라 출발점이다.

2. 리우빌-아널드 정리[편집]

정리의 진술은 이렇다. $2n$ 차원 심플렉틱 다양체 위의 해밀토니안 $H$ 에 대해 함수 $F_1 = H, F_2, \dots, F_n$ 이

\{F_i, F_j\} = 0 \quad (\forall i,j), \qquad \mathrm{d}F_1, \dots, \mathrm{d}F_n \ \text{선형독립}

을 만족한다고 하자. 이때 공통 준위집합 $M_f = \{z : F_i(z) = f_i\}$ 가 콤팩트하고 연결이면, $M_f$ 는 $n$ 차원 원환면 $\mathbb{T}^n$ 에 미분동형이고, 그 근방에서 정준변환으로 작용-각 변수 $(I, \theta)$ 를 잡을 수 있어 방정식이

\dot I = 0, \qquad \dot\theta = \omega(I) = \frac{\partial H}{\partial I}

로 완전히 풀린다. 즉 모든 운동은 각기 다른 원환면 위의 준주기 운동이며, $\omega(I)$ 의 성분들이 유리수 비면 궤도가 닫히고 무리수 비면 원환면을 조밀하게 덮는다. “대합”이라는 조건이 결정적인데, 보존량이 $n$ 개 있어도 서로 괄호가 0이 아니면 이 원환면 구조가 나오지 않는다.²

여기서 푸앵카레 단면이 왜 그렇게 잘 먹히는지도 설명된다. 적분가능계의 단면에는 점들이 매끄러운 닫힌 곡선(원환면의 절단면) 위에만 찍힌다. 단면에서 곡선이 뭉개져 안개처럼 흩어지기 시작하면, 그 계는 적분가능성을 잃은 것이다. 에농-하일레스 계와 표준 사상이 그 붕괴 과정을 보여주는 표준 교보재다.

3. 적분가능성은 기적에 가깝다[편집]

보존량 $n$ 개는 공짜로 안 생긴다. 뇌터 정리가 주는 건 계의 명시적 대칭(시간 병진 → 에너지, 회전 → 각운동량)에 대응하는 것뿐인데, 자유도가 늘면 필요한 보존량 수도 같이 늘어나므로 어지간해서는 부족하다. 케플러 문제에서 추가로 나오는 룽게-렌츠 벡터처럼 “숨은 대칭”이 있는 경우가 있지만, 그건 행운이지 규칙이 아니다.

푸앵카레는 제한 삼체문제에 해석적 보존량이 추가로 존재하지 않음을 보였고, 이후 축적된 결과는 한 문장으로 요약된다. 적분가능계는 해밀토니안계 공간에서 측도 0에 가까운 얇은 집합이다. 교과서에 적분가능계만 잔뜩 나오는 건 그것들이 대표적이어서가 아니라, 그것만 손으로 풀리기 때문이다. 실제 계는 랴푸노프 지수가 양수인 카오스 영역과 잔존 원환면이 뒤섞인 혼합 위상공간을 갖는다.

4. 락스 쌍과 등스펙트럼 흐름[편집]

보존량을 무더기로 만들어 내는 대수적 장치가 락스 쌍(Lax pair)이다. 운동방정식을 두 행렬(또는 연산자) $L, M$ 으로

\dot L = [L, M] = LM - ML

꼴로 쓸 수 있으면, 해는 $L(t) = U(t) L(0) U(t)^{-1}$ 형태의 닮음 변환이 된다( $\dot U = MU$ ). 닮음 변환은 고유값 문제의 스펙트럼을 바꾸지 않으므로 $L$ 의 고유값이 전부 보존량이고, 동등하게 $\mathrm{tr}(L^k)$ 들이 보존된다. 이런 흐름을 등스펙트럼 흐름(isospectral flow)이라 부른다.

토다 사슬이 대표 사례다. 지수 상호작용으로 연결된 입자 사슬인데, 적당한 변수 변환에서 $L$ 이 대칭 삼중대각 행렬이 되고 흐름이 정확히 락스 형태를 띤다. 재미있는 건 그 다음이다 — 토다 흐름을 정수 시각에서 샘플링하면 QR 분해 기반 고유값 알고리즘의 반복과 일치한다. 즉 수치선형대수의 QR 알고리즘이 적분가능계의 이산 시간 흐름이라는 뜻이다.³ 란초스 알고리즘이 삼중대각화를 하는 것과 토다 사슬의 형태가 같은 것도 우연이 아니다.

5. 소리톤과 역산란 변환[편집]

무한 자유도로 가면 이야기가 더 극적이다. KdV 방정식

u_t + 6uu_x + u_{xxx} = 0

은 무한히 많은 보존량을 갖는 적분가능 편미분방정식이고, 소리톤 — 분산과 비선형성이 정확히 상쇄되어 모양을 유지한 채 전파하는 고립파 — 을 해로 갖는다. 소리톤 둘이 충돌하면 잠시 겹쳤다가 원래 모양 그대로 빠져나오고, 남는 것은 위상 이동뿐이다. 선형 파동도 아닌데 이런 일이 벌어지는 게 적분가능성의 위력이다. 사인-고든 방정식(킹크·브리더)과 비선형 슈뢰딩거 방정식(광섬유 펄스)도 같은 가족이다.

푸는 도구는 역산란 변환(IST)이다. 발상은 “비선형 푸리에 변환”에 가깝다. 초기 조건 $u(x,0)$ 을 어떤 선형 산란 문제(KdV의 경우 슈뢰딩거 연산자 $L = -\partial_x^2 + u$ )의 퍼텐셜로 보고 산란 데이터를 뽑으면, 그 데이터가 시간에 대해 선형으로 단순하게 진화한다. 진화시킨 뒤 역산란으로 $u(x,t)$ 를 복원하면 끝. 락스 쌍의 $L$ 이 여기서 산란 연산자로 재등장하고, 이산 고유값 하나가 소리톤 하나에 대응한다.

6. FPUT 문제 — 안 풀린 덕분에 발견된 것[편집]

1953~55년 로스앨러모스에서 페르미, 파스타, 울람이 (그리고 실제 코딩을 한 메리 칭구가⁴) MANIAC으로 돌린 계산이 FPUT 문제다. 64개 질량점을 약한 비선형 스프링으로 이은 1차원 사슬에 에너지를 최저 모드 하나에만 넣고, 등분배 정리대로 에너지가 모든 모드에 골고루 퍼져 열평형에 이르는 것을 보려는 계산이었다.

결과는 정반대였다. 에너지는 몇 개 저차 모드 사이를 왔다 갔다 하다가 거의 처음 상태로 되돌아왔다(재귀, recurrence). 등분배는 일어나지 않았다. 이 “실패”가 두 갈래로 이어진다. 하나는 1965년 자부스키와 크루스칼이 FPUT 사슬의 연속 극한이 KdV 방정식임을 보이고 그 해에 soliton이라는 이름을 붙인 것 — 소리톤 연구의 출발점이다. 다른 하나는 근적분가능계의 KAM 정리적 해석으로, 약한 비선형성에서는 불변 원환면이 대부분 살아남아 위상공간을 가둔다는 그림이다. 에르고딕 가설이 자명하지 않다는 것을 수치실험이 먼저 알려준 사건이기도 하다.

7. 근적분가능계로 가는 다리[편집]

$H(I,\theta) = H_0(I) + \varepsilon H_1(I,\theta)$ 처럼 적분가능계를 살짝 흔든 계를 근적분가능계라 한다. 여기서 벌어지는 일이 20세기 동역학의 본론이다.

KAM: 충분히 무리수적인(디오판토스) 진동수비를 가진 원환면은 작은 $\varepsilon$ 에서 살아남는다. 공명 원환면부터 부서진다.
아널드-모저 층상 구조: 부서진 자리에는 섬 사슬과 카오스 층이 생기고, 자유도 2에서는 살아남은 원환면이 위상공간을 절단해 확산을 막는다.
아널드 확산: 자유도 3 이상에서는 원환면이 위상공간을 나누지 못해, 공명망을 따라 지수적으로 느리지만 전역적인 작용 변수 표류가 가능해진다. 태양계 장기 안정성 논의가 여기 걸려 있다.

수치적으로 이 경계를 탐사할 때는 기하 수치적분 계열이 사실상 필수다. 비심플렉틱 적분기의 인공 소산은 살아있는 원환면을 죽은 것처럼 보이게 만들어, 존재하지 않는 카오스를 관측하게 한다.

8. 관련 문서[편집]

“적분가능”이라는 말은 맥락에 따라 리우빌 적분가능, 대수적 적분가능, 파인레베 성질 등 서로 다른 정의를 가리킨다. 논문에서 이 단어를 만나면 저자가 어느 정의를 쓰는지부터 확인하는 게 안전하다. ↩
예를 들어 3차원 중심력장의 각운동량 세 성분은 보존량 3개지만 서로 괄호가 0이 아니다. 그래서 대합인 조합( $L^2$ , $L_z$ )만 골라 써야 한다. 양자역학 강의에서 “동시에 측정 가능한 관측량” 얘기를 하는 그 구조와 같은 것이다. ↩
Symes (1982). 처음 이 사실을 접하면 “그러니까 내가 고유값 뽑을 때마다 적분가능계를 하나씩 시간 적분하고 있었다는 거냐”는 기분이 든다. 대략 맞다. ↩
원 보고서(Los Alamos LA-1940, 1955)는 Fermi–Pasta–Ulam으로 인용돼 왔지만, MANIAC 프로그래밍을 실제로 수행한 Mary Tsingou의 이름이 보고서에 감사로만 남아 있었다. 2008년 Dauxois의 문제제기 이후 FPUT로 표기하는 것이 표준이 되어 가는 중이다. ↩