락스 쌍

수치해석 계산물리 마지막 수정: 2026-08-10 04:47:51

1. 개요[편집]

락스 쌍 Lax Pair
락스 방정식	$L_t = [M, L] = ML - LM$
제안	1968년, 피터 락스
귀결	$L$ 의 스펙트럼이 시간 불변 → 고유값이 곧 보존량
별명	등스펙트럼 흐름(isospectral flow)
확장	영곡률 표현, AKNS 형식
주의	존재는 충분조건적 도구이지 적분가능성의 판정법이 아니다

락스 쌍(Lax pair)은 비선형 발전방정식을 두 선형 연산자 $L, M$ 에 대한 $L_t = [M, L]$ 이라는 형태로 다시 쓰는 표현이며, 이 형태가 성립하면 $L$ 의 스펙트럼이 시간에 대해 변하지 않는다. 1968년 피터 락스(Peter Lax)가 KdV 방정식의 무한 보존량과 역산란 절차를 통일적으로 설명하기 위해 도입했다.

핵심 논리는 세 줄이면 끝난다. $L\psi = \lambda\psi$ 인 고유함수를 $\psi_t = M\psi$ 로 진화시키고 $\lambda_t = 0$ 을 요구한 뒤 양변을 미분하면

L_t\psi + LM\psi = \lambda M \psi = ML\psi \quad\Longrightarrow\quad L_t = [M, L]

가 나온다. 거꾸로 읽으면, 락스 방정식이 성립하도록 $L, M$ 을 짤 수 있다면 고유값이 전부 보존량이라는 뜻이다.¹ 적분가능계 문서가 다루는 리우빌-아널드 그림에서 “보존량을 어디서 그렇게 많이 구해오느냐”에 대한 대수적 대답이 바로 이것이다.

2. 등스펙트럼 흐름[편집]

락스 방정식의 해는 다음과 같이 쓸 수 있다. $\dot U = MU$ , $U(0) = I$ 로 정의된 $U(t)$ 를 잡으면

L(t) = U(t)\,L(0)\,U(t)^{-1}

이다. 즉 $L$ 은 닮음 변환을 따라 흘러갈 뿐 스펙트럼은 그대로다. 유한차원 행렬이면 $\mathrm{tr}(L^k)$ 가 전부 보존되고(고유값의 대칭함수), 무한차원 미분연산자면 산란 데이터의 이산 고유값이 보존된다. 이런 흐름을 등스펙트럼 흐름(isospectral flow)이라 부른다. $M$ 이 반대칭이면 $U$ 가 직교행렬이 되어 흐름이 QR 분해나 슈어 분해와 같은 직교 닮음의 세계에 머무른다.

부호 규약은 문헌마다 다르다. $L_t = [M, L]$ 과 $L_t = [L, M]$ 이 모두 통용되며, $M \to -M$ 으로 서로 옮겨간다. 결론(스펙트럼 불변)은 같으니 겁먹을 필요는 없지만, 남의 논문에서 $M$ 을 그대로 베껴 오면 부호 하나 때문에 방정식이 시간 역방향으로 나오는 사고가 흔하다.

3. KdV의 락스 쌍[편집]

락스가 든 원조 예시다. 연산자를

L = -\partial_x^2 + u, \qquad M = -4\partial_x^3 + 6u\,\partial_x + 3u_x

로 놓고 $[M, L]$ 을 직접 전개하면(3계 미분 항까지 전부 상쇄되고 곱셈 연산자만 남는다)

[M, L] = 6uu_x - u_{xxx}

가 되어, $L_t = u_t$ 와 맞추면

u_t - 6uu_x + u_{xxx} = 0

이 정확히 나온다. 여기서 부호 규약을 한 번 짚고 가야 한다. 위 계산은 $u_t - 6uu_x + u_{xxx}=0$ 을 준다 — $u_t + 6uu_x + u_{xxx}=0$ 을 쓰고 싶으면 $u \to -u$ 로 바꾸거나 $L = -\partial_x^2 - u$ 를 써야 한다. 두 관례가 문헌에 혼재하며, 슈뢰딩거 작용소의 퍼텐셜 부호가 뒤집히면 속박상태(=솔리톤) 개수가 달라지므로 그냥 표기 취향이 아니다.²

이 $L$ 이 바로 슈뢰딩거 방정식의 정상상태 연산자이고, $u$ 가 퍼텐셜 역할을 한다는 점이 결정적이다. 그래서 다음 순서가 자연스럽게 이어진다 — 퍼텐셜 $u(x,0)$ 의 산란 문제를 풀고, 산란 데이터를 (이제 고유값이 상수임을 아니까) 선형으로 진화시키고, 되돌린다. 이것이 역산란 변환이다. 락스 쌍은 그 프로그램의 입장권이다.

4. 영곡률 표현과 AKNS[편집]

미분연산자를 직접 짜는 락스 형식은 방정식마다 손으로 맞춰야 해서 확장성이 나쁘다. 현대적인 표준은 영곡률 표현(zero-curvature representation)이다. 스펙트럼 파라미터 $\zeta$ 를 품은 두 행렬 $U(\zeta), V(\zeta)$ 로 선형계

\psi_x = U\psi, \qquad \psi_t = V\psi

를 놓고 교차미분의 정합성 $\psi_{xt} = \psi_{tx}$ 을 요구하면

U_t - V_x + [U, V] = 0

이 되고, 이 조건이 곧 원래의 비선형 방정식이 되도록 $U, V$ 를 고른다. 기하학적으로는 접속(connection)의 곡률이 0이라는 뜻이라 게이지 이론의 언어와 그대로 이어진다.

AKNS(Ablowitz-Kaup-Newell-Segur, 1974) 형식은 $2\times2$ 행렬로 이 틀을 표준화한 것이다. 산란 문제를

\psi_x = \begin{pmatrix} -i\zeta & q \\ r & i\zeta \end{pmatrix}\psi

로 고정하고 시간 부분 $V$ 만 바꾸면, $r = -q^*$ 에서 비선형 슈뢰딩거 방정식이, $r = -q$ 에서 mKdV가, 적절한 변형에서 사인-고든 방정식이 한 가족으로 쏟아져 나온다. 한 번의 산란 문제 설정이 방정식 여러 개를 커버한다는 것이 이 형식의 위력이고, 새 적분가능 방정식을 찾는 작업이 “새 $V$ 를 찾는 작업”으로 환원된다.

5. 토다 격자와 QR 알고리즘[편집]

수치해석 하는 사람에게 락스 쌍이 남의 이야기가 아닌 이유가 여기 있다. 지수 상호작용으로 연결된 입자 사슬인 토다 격자를 적절한 변수(플라스키 변수)로 쓰면 $L$ 이 대칭 삼중대각 행렬이 되고, 흐름이

\dot L = [B, L], \qquad B = L_+ - L_-

형태가 된다( $L_\pm$ 는 상·하 삼각 부분, $B$ 는 반대칭). 여기까지는 적분가능계 이야기지만, 사임스(Symes, 1982)와 데이프트-난다-토메이(Deift-Nanda-Tomei, 1983)가 밝힌 다음 사실이 결정적이다.

e^{t L(0)} = Q(t)R(t) \quad\text{(QR 분해)}, \qquad L(t) = Q(t)^{\mathsf T} L(0)\, Q(t)

즉 토다 흐름의 시각 $t$ 에서의 해는 $e^{tL(0)}$ 의 QR 분해로 주어지고, 정수 시각에서 샘플링하면 정확히 QR 알고리즘의 반복열과 일치한다. 대칭 고유값 문제를 푸는 그 QR 알고리즘이, 알고 보니 적분가능 해밀토니안계의 시간 1 이산화였던 것이다.

이 대응은 장식이 아니라 실제 통찰을 준다.

QR 알고리즘이 삼중대각 구조를 보존하는 것은 토다 흐름이 삼중대각 다양체를 보존하는 것과 같은 사실이다(삼중대각화와 란초스 알고리즘이 이 그림에 함께 들어온다).
반복이 고유값 순서대로 정렬되며 수렴하는 것은 토다 흐름이 $t \to \infty$ 에서 대각 행렬로 가는 것과 대응한다.
시프트를 넣은 QR은 시간 재매개화된 토다 흐름에 해당하며, 수렴 가속의 기하학적 의미가 여기서 나온다.

거꾸로 흐름 쪽에서 알고리즘 쪽으로 배우는 것도 있다. $\dot L = [B, L]$ 을 일반 룽게-쿠타법으로 적분하면 고유값이 표류한다 — 닮음 구조가 이산 수준에서 깨지기 때문이다. 처방은 $L$ 을 직접 전진시키지 말고 $\dot U = BU$ 를 직교군 위에서 적분한 뒤 $L = U^{\mathsf T}L_0U$ 로 복원하는 것이고, 이것이 리 군 적분기와 기하 수치적분이 다루는 등스펙트럼 적분기다.³

6. 락스 쌍이 있으면 적분가능한가[편집]

여기서 뉘앙스를 지켜야 한다. 흔히 “락스 쌍이 있으면 적분가능”이라고 요약되지만, 엄밀하게는 그렇게 단순하지 않다.

자명한 락스 쌍은 언제나 만들 수 있다. 스펙트럼 파라미터가 들어 있지 않은 락스 쌍은 정보를 거의 담지 못하며, 실제로 쓸모 있는 것은 $\zeta$ 에 의존하는 락스 쌍이다. $\zeta$ 의 각 차수에서 보존량이 계열로 뽑혀 나오기 때문이다.
보존량이 많다고 곧 리우빌 적분가능은 아니다. 고유값이 보존된다는 것과, 그것들이 함수적으로 독립이며 서로 포아송 괄호가 0(대합)이라는 것은 별개의 주장이다. 대합성은 고전 $r$ -행렬 구조가 있을 때 따라 나오며, 락스 쌍을 찾은 다음 별도로 확인해야 하는 항목이다.
없다고 비적분이라는 증명도 아니다. 주어진 방정식에 락스 쌍이 존재하는지를 판정하는 일반 알고리즘은 없다. 실무에서 쓰는 것은 파인레베 시험, 대칭·보존법칙 계산, 3-솔리톤 조건 같은 필요조건 검사이고, 그것들을 통과한 뒤에야 손으로 락스 쌍을 찾아 나선다.

정리하면 락스 쌍은 적분가능성을 확인하고 활용하기 위한 매우 강력한 도구이지, 보편적 판정 기준이 아니다. 논문에서 “락스 쌍을 찾았으므로 이 계는 적분가능하다”는 문장을 보면 스펙트럼 파라미터가 실제로 들어 있는지부터 확인하는 것이 습관이 된다.⁴

7. 관련 문서[편집]

이 유도에서 진짜 어려운 부분은 논리가 아니라 $M$ 을 찾는 것이다. $L$ 은 대개 물리에서 자연스럽게 주어지지만( $u$ 를 퍼텐셜로 갖는 산란 연산자 등), 짝이 되는 $M$ 은 차수를 하나씩 올려가며 미정계수법으로 두들겨 맞추는 수작업이다. 컴퓨터 대수의 도움을 받아도 여전히 노가다. ↩
심위키 안에서도 KdV 방정식 문서는 $u_t + 6uu_x + u_{xxx}=0$ 관례를 쓴다. 두 문서를 나란히 놓고 부호가 다르다고 당황하지 말 것 — $u \to -u$ 한 번이면 같은 이야기다. 이런 사소한 관례 차이가 코드 이식에서 반나절씩 잡아먹는다는 게 문제일 뿐. ↩
“고유값 뽑는 알고리즘을 돌리는 것이 사실은 적분가능계를 시간적분하는 것”이라는 문장은 처음 들으면 그냥 신기한 얘기 같지만, 등스펙트럼 적분기 설계에서는 실전 지침이 된다. 보존해야 할 구조가 무엇인지 알면 어떤 적분기를 쓰면 안 되는지도 알게 되니까. ↩
사실 “적분가능”이라는 단어 자체가 리우빌 적분가능·대수적 적분가능·파인레베 성질 등 여러 정의로 쓰인다는 점도 같이 봐야 한다. 이 바닥에서 용어 정의부터 확인하는 것은 트집이 아니라 생존 기술이다. ↩