역산란 변환

1. 개요[편집]

비선형 방정식을 푸는 가장 확실한 방법은, 그것이 사실 선형 방정식이었음을 밝히는 것이다.

역산란 변환(inverse scattering transform, IST)은 특정 비선형 편미분방정식의 초기값 문제를, 초기 파형을 어떤 선형 산란 문제의 퍼텐셜로 읽어 산란자료로 옮긴 뒤 그 산란자료를 선형으로 시간 진화시키고 다시 물리공간으로 되돌리는 풀이법이다. 1967년 가드너-그린-크루스칼-미우라(GGKM)가 KdV 방정식에 대해 발견했고, 비선형 방정식에 푸리에 변환에 해당하는 도구가 존재한다는 사실 자체가 20세기 수리물리의 사건이었다.

이 문서는 산란자료 관점에 집중한다. 솔리톤이 물리공간에서 어떻게 충돌하고 위상만 밀리는지, 락스 쌍이 어떻게 유도되는지는 각각 그 문서에 정리돼 있으니 여기서 반복하지 않는다. 여기서 볼 것은 “그 충돌 서사를 처음부터 자명하게 만드는 좌표계가 무엇인가”다.

2. 왜 「비선형 푸리에 변환」인가[편집]

선형 상수계수 PDE를 푸는 절차를 떠올려 보자.

1. 초기조건을 푸리에 변환해 $\hat u(k,0)$ 을 얻는다.
2. 각 모드가 독립적으로 $\hat u(k,t) = \hat u(k,0)e^{-i\omega(k)t}$ 로 선형·비결합 진화한다.
3. 역변환해 $u(x,t)$ 를 복원한다.

IST는 이 세 단계를 그대로 흉내 낸다. 다만 1단계와 3단계의 변환이 적분 커널이 아니라 산란 문제 풀이로 바뀐다.

선형 세계	IST
푸리에 변환	직접 산란 (direct scattering)
스펙트럼 $\hat u(k)$	산란자료 (연속 스펙트럼 + 이산 고유값)
모드별 위상 회전	산란자료의 선형 시간 전개
푸리에 역변환	겔판트-레비탄-마르첸코 방정식
정현파	솔리톤 + 분산 꼬리

결정적인 차이는 이산 스펙트럼이다. 선형 세계의 푸리에 스펙트럼은 연속뿐이지만, 산란 문제의 스펙트럼에는 속박상태에 해당하는 이산 고유값이 함께 나온다. 그리고 이산 고유값 하나가 솔리톤 하나다. 초기조건을 보고 몇 개의 솔리톤이 나올지 미리 셀 수 있다는 마법이 여기서 나온다.

3. 부호 규약 먼저[편집]

문헌에 두 관례가 혼재하므로 못 박고 시작한다. 이 문서는 GGKM 원본을 따라

$$u_t - 6uu_x + u_{xxx} = 0, \qquad L = -\partial_x^2 + u(x,t)$$

를 쓴다. 이 짝에서 솔리톤은 $u<0$ 인 퍼텐셜 우물이고, $L$ 의 속박상태가 존재한다. $u_t + 6uu_x + u_{xxx}=0$ 을 쓰고 싶으면 $u \to -u$ 로 바꾸면 되며, 이때는 $L=-\partial_x^2-u$ 를 써야 한다. 퍼텐셜 부호가 뒤집히면 속박상태 개수가 통째로 달라지므로 이건 취향 문제가 아니다.¹ 락스 쌍 문서에 이 규약 비교가 자세히 있다.

락스 형식으로는 짝

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

가 $\partial_t L = [M, L]$ 을 만족할 때 위 KdV가 나온다. 여기서 읽어야 할 것은 딱 하나다 — $L_t = [M,L]$ 은 $L$ 이 닮음 변환으로 흘러간다는 뜻이고, 따라서 $L$ 의 스펙트럼은 시간에 대해 불변이다. 그러면 스펙트럼을 좌표로 쓰지 않을 이유가 없다.

4. 직접 산란 — 산란자료가 무엇인가[편집]

$u(\cdot,t)$ 를 고정된 퍼텐셜로 보고 정상상태 슈뢰딩거 방정식

$$-\psi_{xx} + u(x,t)\,\psi = \lambda\,\psi$$

을 푼다. $|x|\to\infty$ 에서 $u\to 0$ 이라 가정하면 스펙트럼이 두 부분으로 나뉜다.

연속 스펙트럼 $\lambda = k^2 > 0$. 오른쪽에서 입사하는 파동에 대해

$$\psi(x,k) \sim \begin{cases} t(k)\,e^{-ikx}, & x \to -\infty \\[2pt] e^{-ikx} + r(k)\,e^{ikx}, & x \to +\infty \end{cases}$$

로 두면 투과계수 $t(k)$ 와 반사계수 $r(k)$ 가 정의되고 $|r|^2 + |t|^2 = 1$ 이다. 물리적으로 $r(k)$ 는 “이 파형이 그 파수에서 얼마나 산란시키는가”이고, IST 안에서는 최종적으로 **분산 꼬리(radiation)**로 나타날 성분이다.

이산 스펙트럼 $\lambda_n = -\kappa_n^2 < 0$, $n=1,\dots,N$, $\kappa_1 > \kappa_2 > \cdots > \kappa_N > 0$. 대응하는 속박상태 $\psi_n$ 을 $\int_{-\infty}^{\infty}\psi_n^2\,dx = 1$ 로 규격화하고 꼬리를

$$\psi_n(x) \sim c_n\,e^{-\kappa_n x}, \qquad x \to +\infty$$

로 쓸 때의 $c_n$ 이 정규화 상수다. 이 셋을 묶은

$$S(t) = \Bigl\{\, \kappa_n,\ c_n(t)\ (n=1,\dots,N);\ \ r(k,t)\ (k\in\mathbb{R}) \,\Bigr\}$$

이 산란자료이고, 이것이 IST의 “스펙트럼”이다. 1차원에서 이 자료가 퍼텐셜을 유일하게 결정한다는 것이 역산란 이론(겔판트-레비탄, 마르첸코, 파데예프)의 내용이다. 이 방향의 문제 — 관측된 산란 정보로 안쪽을 복원하기 — 는 역문제의 원형이기도 하다.

5. 산란자료의 시간 불변성[편집]

여기가 이 이론의 심장이다. $u$ 가 KdV를 만족하며 흘러갈 때 산란자료의 진화는 다음과 같다.

$$\frac{d\kappa_n}{dt} = 0, \qquad c_n(t) = c_n(0)\,e^{4\kappa_n^3 t}, \qquad r(k,t) = r(k,0)\,e^{8ik^3 t}$$

읽는 법은 이렇다.

• $\kappa_n$ 은 상수다. 등스펙트럼 흐름의 직접적 귀결이며, 물리적으로는 “솔리톤의 진폭과 속도가 영원히 변하지 않는다”에 해당한다. $\kappa_n$ 하나가 진폭 $2\kappa_n^2$, 속도 $4\kappa_n^2$ 의 솔리톤 하나를 지정한다. 큰 솔리톤이 빠르다는 KdV의 유명한 성질이 여기서는 그냥 지수 하나다.
• $c_n$ 은 순수 지수로 자란다. 실수 지수이고, 물리공간에서는 $n$ 번째 솔리톤의 위치($x_n = 4\kappa_n^2 t + \delta_n$)를 옮기는 역할만 한다. 진폭에는 손대지 않는다.
• $r(k)$ 는 순수 위상 회전이다. $|r(k,t)| = |r(k,0)|$ 이 시간에 무관하다. 분산 꼬리의 파수별 세기는 보존되고 위상만 $8k^3 t$ 만큼 돌아간다 — 선형 KdV의 분산 관계 $\omega = -k^3$ 이 여기 그대로 박혀 있는 셈이다.

정리하면 산란자료 공간에서 KdV는 대각화되어 있다. $N+1$ 개의 결합된 비선형 상호작용이 사라지고, 서로 이야기하지 않는 지수함수들만 남는다. 솔리톤 충돌이 진폭을 보존하고 위상만 옮기는 이유를 물리공간에서 설명하려면 계산이 필요하지만, 산란자료 공간에서는 설명할 것 자체가 없다. $\kappa_n$ 이 상수니까.

이 구조를 해밀토니안 역학의 언어로 옮기면 더 선명하다. 자하로프-파데예프(1971)는 KdV를 무한차원 해밀토니안계로 보았을 때 산란자료가 정확히 작용-각 변수 역할을 함을 보였다. $\kappa_n$ 과 $|r(k)|$ 계열이 작용(불변량), $c_n$ 의 위상과 $\arg r(k)$ 가 각(선형 증가)이다. 적분가능계에서 리우빌-아널드 정리가 유한 자유도에 대해 약속하는 그림이, 무한 자유도에서 구체적으로 실현된 사례다. KdV의 무한개 보존량도 $\ln|t(k)|$ 의 $k\to\infty$ 점근 전개 계수로 한꺼번에 튀어나온다 — 보존량을 하나씩 손으로 찾던 미우라-가드너-크루스칼의 노동이 여기서는 한 줄이 된다.

6. 역산란 — 겔판트-레비탄-마르첸코[편집]

진화시킨 산란자료에서 $u(x,t)$ 를 되찾는 단계다. 커널

$$F(\xi; t) = \sum_{n=1}^{N} c_n^2(t)\,e^{-\kappa_n \xi} + \frac{1}{2\pi}\int_{-\infty}^{\infty} r(k,t)\,e^{ik\xi}\,dk$$

를 만들고, 미지함수 $K(x,y;t)$ 에 대한 선형 적분방정식

$$K(x,y) + F(x+y) + \int_x^{\infty} K(x,z)\,F(z+y)\,dz = 0, \qquad y > x$$

를 푼다. 이것이 겔판트-레비탄-마르첸코 방정식(GLM)이다. 해가 나오면 퍼텐셜은

$$u(x,t) = -2\,\frac{d}{dx}K(x,x;t)$$

로 복원된다. 중요한 것은 이 적분방정식이 선형이라는 점이다. 비선형성은 전부 1단계(직접 산란)의 스펙트럼 구조 안으로 흡수됐고, 남은 것은 프레드홀름 제2종 적분방정식 하나다.

가장 짧은 검산: $N=1$, $r\equiv 0$ 이면 $F(\xi)=c^2e^{-\kappa\xi}$ 이고 GLM이 분리형(degenerate kernel)이라 손으로 풀린다.

$$K(x,y) = -\frac{c^2 e^{-\kappa(x+y)}}{1 + \dfrac{c^2}{2\kappa}e^{-2\kappa x}} \ \Longrightarrow\ u(x,t) = -2\kappa^2 \operatorname{sech}^2\!\bigl(\kappa(x - 4\kappa^2 t) - \delta\bigr)$$

솔리톤이 나온다. $N$ 개 이산 고유값이면 커널이 $N$ 항 분리형이 되어 GLM이 $N\times N$ 선형계로 축약되고, 행렬식 꼴의 $N$-솔리톤 해가 곧장 나온다.

7. 반사 없는 퍼텐셜 — 솔리톤을 미리 세기[편집]

$r(k)\equiv 0$ 인 초기조건을 반사 없는 퍼텐셜(reflectionless potential)이라 하고, 이 경우 분산 꼬리가 전혀 없이 순수 $N$-솔리톤만 나온다. 교과서 표준 예가

$$u(x,0) = -N(N+1)\operatorname{sech}^2 x \quad \Longrightarrow \quad \kappa_n = n\ \ (n=1,\dots,N),\ \ r \equiv 0$$

이다. 즉 $u(x,0)=-6\operatorname{sech}^2x$ 는 정확히 2-솔리톤으로 갈라지고, $-12\operatorname{sech}^2x$ 는 3-솔리톤이 된다.

더 일반적으로 $u(x,0) = -U_0\operatorname{sech}^2(x/L)$ 의 속박상태 개수는

$$N = \left\lceil \frac{\sqrt{1+4U_0L^2} - 1}{2} \right\rceil$$

로 주어진다. $U_0 = N(N+1)$, $L=1$ 을 넣으면 정확히 $N$ 이 나온다. 이 공식의 실용적 함의가 크다.

• 초기 파형만 보고 최종 솔리톤 개수와 진폭을 계산으로 알 수 있다. 수치 실험을 돌리기 전에 답을 알고 시작하는 셈이라, KdV 솔버 검증 및 확인의 표준 시험 문제가 된다.
• 1차원 인력 퍼텐셜은 아무리 얕아도 속박상태를 적어도 하나 갖는다. 따라서 $u(x,0)\le 0$ 이고 어딘가 진짜 음수라면, 아무리 작은 파라도 최소 하나의 솔리톤은 반드시 나온다. 3차원이었다면 성립하지 않았을 이야기다.
• 반대로 $u(x,0) \ge 0$ (이 규약에서 척력)이면 속박상태가 없어 솔리톤이 하나도 안 나오고, 초기 파형은 전부 분산 꼬리로 흩어진다. 그 꼬리의 장시간 거동은 에어리 함수형으로 $t^{-1/3}$ 감쇠한다.

8. AKNS와 그 이후[편집]

KdV만의 요행이 아니었다는 점이 중요하다.

• 1972년 자하로프-샤바트가 $2\times2$ 산란 문제로 비선형 슈뢰딩거 방정식을 풀었다. 광섬유 펄스의 솔리톤 성분이 여기서 나온다.
• 1973~74년 AKNS(Ablowitz-Kaup-Newell-Segur)가 이를 표준화해, 하나의 산란 문제 틀에서 NLS·mKdV·사인-고든 방정식을 한 가족으로 다뤘다. 형식의 세부는 락스 쌍 문서에 있다.
• 주기 경계에서는 산란자료 대신 유한대역(finite-gap) 구조와 리만 세타 함수가 등장한다. 노비코프-두브로빈-매큔-반 뫼르베케 계열의 이론이며, 무한 구간 IST의 주기판 대응물이다.
• 장시간 점근은 리만-힐베르트 문제로 재정식화한 뒤 데이프트-저우의 비선형 최급강하법으로 다룬다. “솔리톤들이 앞서 나가고 뒤에 분산 꼬리가 남는다”는 그림이 여기서 정리로 승격된다.

이런 확장에도 불구하고 IST가 적용되는 방정식의 목록은 여전히 매우 짧다. 감쇠·강제·2차원화·변계수 중 무엇 하나만 들어와도 적분가능성이 깨진다. 적분가능계 문서의 표현대로 적분가능성은 기적에 가깝고, IST는 그 기적을 최대한 착취하는 도구다.

9. 계산 도구로서[편집]

솔직하게 적자면, IST는 KdV를 푸는 실무 수치기법이 아니다. 실제 KdV 초기값 문제는 의사스펙트럼 미분 + 지수 적분기(ETDRK4) 조합이 훨씬 빠르고 일반적이다. IST가 계산에서 갖는 값어치는 다른 데 있다.

• 정확해 생성기. $N$-솔리톤 해석해를 임의 정밀도로 만들어 솔버의 수렴차수를 재고, 충돌 후 진폭 복원과 위상 이동을 이론값과 대조한다. 수치 소산과 분산이 있는 도식은 이 시험을 절대 통과하지 못한다.
• 진단 도구 — 비선형 스펙트럼 해석. 측정된 파형에 직접 산란만 적용해 “이 데이터에 솔리톤 성분이 몇 개, 어떤 진폭으로 들어 있는가”를 뽑는다. 해양 파랑 시계열 분석, 광섬유 펄스 진단이 대표 용도이며, 이 방향을 비선형 푸리에 해석이라 부른다. 광통신에서는 아예 정보를 산란자료에 실어 보내자는 발상(비선형 주파수 분할 다중화)까지 나왔다 — 채널의 비선형성이 산란자료 공간에서는 대각이니까.
• 수치 IST 자체. 직접 산란은 결국 고유값 문제이고, 역산란은 프레드홀름 적분방정식 또는 리만-힐베르트 문제의 수치해다. 후자를 안정적으로 이산화하면 시간 $t$ 에 무관한 비용으로 임의의 $t$ 에서 해를 평가할 수 있다는 특이한 장점이 생긴다. 시간 적분기는 $t$ 가 커지면 그만큼 스텝을 밟아야 하지만, IST는 $c_n(t)$ 와 $r(k,t)$ 에 지수를 곱하기만 하면 된다.²

10. 여담[편집]

• 이 모든 것의 출발점은 페르미-파스타-울람-칭구 문제의 “실패한” 수치실험이었다. 등분배가 안 일어난 이유를 캐다가 자부스키-크루스칼이 연속체 극한에서 KdV를 얻었고, 그 KdV의 이상한 안정성을 설명하려다 GGKM이 IST에 도달했다. 수치실험 하나가 새로운 수학 분야를 낳은 몇 안 되는 사례다.³
• GGKM 논문은 4페이지짜리 Physical Review Letters다. 결과가 얼마나 이상했던지, 당시에는 “KdV라는 특수한 방정식의 요행”으로 여겨졌다는 회고가 많다. 락스가 이듬해 락스 쌍으로 구조를 드러내고, 자하로프-샤바트가 NLS로 확장하고 나서야 이것이 방법론임이 인정됐다.
• 산란자료 관점이 주는 부수적 교훈 하나 — 어떤 문제가 도저히 안 풀린다면, 그 문제가 어려운 것이 아니라 좌표계가 나쁜 것일 수 있다.⁴

11. 관련 문서[편집]

• KdV 방정식 · 솔리톤 · 락스 쌍 · 적분가능계
• 비선형 슈뢰딩거 방정식 · 사인-고든 방정식 · 페르미-파스타-울람-칭구 문제
• 슈뢰딩거 방정식 · 고유값 문제 · 역문제 · 푸리에 변환
• 겔판트-레비탄-마르첸코 방정식 · 리만-힐베르트 문제 · 비선형 푸리에 해석
• 해밀토니안 역학 · 포아송 괄호 · 의사스펙트럼 · 지수 적분기

12. Footnotes[편집]

1. 이 부호를 잘못 짚으면 “솔리톤이 하나도 안 나오는데요?”라는 상황이 벌어진다. 코드가 아니라 규약이 틀린 것이므로 디버거로는 절대 못 찾는다. 논문 식을 옮길 때 KdV의 부호부터 확인하는 습관을 들이자. ↩

2. 물론 공짜는 아니다. 직접 산란 단계에서 이산 고유값을 놓치거나 $r(k)$ 를 부정확하게 잡으면 그 오차가 모든 $t$ 에 그대로 실린다. 그리고 초기조건이 충분히 빨리 감쇠하지 않으면 이론의 전제가 무너진다. ↩

3. 페르미는 결과를 보지 못하고 1954년에 세상을 떠났다. 본인 이름이 붙은 문제가 솔리톤과 IST와 KAM 정리 논의를 동시에 열어젖힐 줄은 몰랐을 것이다. ↩

4. 이 교훈의 수치해석판이 전처리(preconditioning)다. 좌표계를 바꿔 문제를 쉽게 만든다는 점에서 정신은 같은데, 아쉽게도 IST처럼 문제를 완전히 대각화해 주는 전처리자는 아직 없다. ↩