페르미-파스타-울람-칭구 문제

계산물리 물리 마지막 수정: 2026-08-04 04:34:38

1. 개요[편집]

페르미-파스타-울람-칭구 문제 Fermi–Pasta–Ulam–Tsingou Problem
수행	로스앨러모스, 1953~1955 · MANIAC I
보고서	LA-1940, "Studies of Nonlinear Problems I" (1955)
계	양 끝 고정된 1차원 비선형 질점 사슬
기대	에너지가 모든 정규모드에 등분배
결과	등분배 실패 — 초기 상태로 되돌아옴(재귀)
후속	KdV 연속체 극한 · 솔리톤(1965) · 준상태 · 강한 확률화 문턱

계산이 틀린 줄 알고 며칠을 뒤졌는데, 코드가 맞았다. 물리 쪽이 틀린 거였다.

페르미-파스타-울람-칭구 문제(FPUT problem)는 1950년대 로스앨러모스에서 비선형 스프링으로 이은 1차원 질점 사슬을 컴퓨터로 적분해 에너지의 등분배를 확인하려 했으나, 에너지가 몇 개 저차 모드 사이만 오가다가 초기 상태로 거의 되돌아온 수치실험이자, 그 결과가 제기한 문제 전체를 가리킨다. 계산기로 물리 문제를 “실험”한 최초의 사례로 꼽히며, 여기서 솔리톤·카오스 전이·비선형 격자 열전도 연구가 통째로 갈라져 나왔다.

이름이 넷인 것은 뒤늦은 교정 때문이다. 원 보고서의 저자는 페르미·파스타·울람 셋이지만, MANIAC I에서 실제로 코드를 짜고 돌린 사람은 메리 칭구(Mary Tsingou)였고 보고서 본문에도 그 사실이 적혀 있다. 2008년 도수아의 글을 계기로 이름을 넷으로 부르자는 제안이 퍼졌고, 지금은 FPUT 표기가 표준에 가깝다.¹

2. 1955년, MANIAC 위에서 벌어진 일[편집]

계는 단순하다. 양 끝이 고정된 사슬에 질점 $N$ 개를 놓고, 이웃을 잇는 스프링에 약한 비선형항을 넣는다. 이차 비선형(α 모형)의 운동방정식은 다음과 같다.

\ddot{x}_i = (x_{i+1} + x_{i-1} - 2x_i) \;+\; \alpha\left[(x_{i+1}-x_i)^2 - (x_i - x_{i-1})^2\right]

세제곱 비선형을 쓰면 β 모형이 된다. $\alpha = 0$ 이면 이 계는 완전히 풀린다. 정규모드로 바꾸면 결합이 사라져 독립적인 조화진동자 다발이 되고, 각 모드의 에너지가 개별적으로 보존된다. 즉 선형 사슬은 적분가능계이고, 에너지를 한 모드에 넣으면 영원히 그 모드에 있다.

여기에 비선형항을 조금 넣으면 모드끼리 결합이 생긴다. 당시의 상식(등분배 정리와 에르고딕 가설)은 결합이 아무리 약해도 시간만 충분하면 에너지가 모든 모드에 골고루 퍼져 각 모드가 평균 $k_B T/2$ 씩 갖는 상태로 간다는 것이었다. 페르미가 확인하려던 것도 그 “완화 과정”이 어떻게 진행되는지, 얼마나 걸리는지였다. 계산은 질점 32개(보고서에는 16개·64개 계산도 있다)에 최저 모드 하나에만 에너지를 넣고 시작했다.

3. 등분배 대신 나온 것 — 재귀[편집]

초반은 예상대로였다. 에너지가 모드 1에서 2, 3, 4로 흘러 나갔다. 그런데 몇 개 저차 모드에 퍼지고 나서는 더 이상 고차 모드로 가지 않았고, 얼마 뒤 에너지가 방향을 되돌려 모드 1로 다시 모였다. 기본 모드 주기의 150여 배쯤 되는 시점에 초기 에너지의 대부분이 모드 1에 복귀했다. 이것이 FPUT 재귀(recurrence)다.

E_k(t) = \frac{1}{2}\left(\dot{A}_k^2 + \omega_k^2 A_k^2\right), \qquad \omega_k = 2\sin\frac{k\pi}{2(N+1)}

로 정의한 모드 에너지 $E_k(t)$ 를 시간에 대해 그리면, 몇 개의 저차 모드가 서로 주고받으며 준주기적으로 진동하는 그림이 나온다. 열평형의 흔적은 없다. 원 보고서의 표현이 유명한데, 요지는 “우리가 기대한 것과는 다른 결과가 나왔고, 이것이 이 계산의 주된 결론”이라는 것이다.

이야기는 여기서 끝나지 않았다. 1972년 턱과 멘첼은 이 계산을 훨씬 오래 돌려서, 재귀가 완벽하지 않고 조금씩 어긋난 뒤 훨씬 긴 시간 후에 초기 상태에 더 가깝게 되돌아오는 초주기(superperiod)가 있음을 보였다. 재귀 위에 재귀가 얹힌 구조다.

4. 자부스키-크루스칼 — 연속체 극한과 솔리톤[편집]

1965년 자부스키와 크루스칼이 이 사슬을 연속체 극한으로 보냈다. 격자 간격이 작고 진폭이 작으며 파장이 긴 영역에서 $x_i(t) \to u(x,t)$ 로 놓고 전개하면, 이산 격자의 이차 비선형과 격자 분산이 각각 살아남아 KdV 방정식이 나온다.

u_t + u\,u_x + \delta^2 u_{xxx} = 0

여기가 결정적인 지점이다. 비선형항 $u u_x$ 은 파형의 앞을 세워 충격파로 만들려 하고, 분산항 $\delta^2 u_{xxx}$ 은 파형을 퍼뜨리려 한다. 둘이 정확히 균형을 이루면 모양을 유지한 채 달리는 국소 파동이 생긴다. 자부스키와 크루스칼은 이 해들이 서로 충돌한 뒤에도 모양과 속도를 그대로 유지하고 위상만 어긋난 채 지나간다는 것을 수치적으로 발견하고, 입자 같다는 뜻에서 솔리톤(soliton)이라 이름 붙였다.

재귀의 설명도 여기서 나온다. 초기의 매끄러운 사인파가 여러 개의 솔리톤으로 분해되어 각자 다른 속도로 주기 영역을 돌다가, 어느 순간 위상이 다시 맞아떨어지면서 원래 파형을 거의 복원한다. 재귀는 통계적 기적이 아니라 결정론적 위상 정렬이었던 것이다.

이 발견이 KdV 방정식의 완전 적분가능성, 역산란 변환(1967), 무한히 많은 보존량, 락스 쌍으로 이어지는 1970년대 가적분계 폭발의 시작점이 됐다. 적분가능계 문서의 소리톤 절이 그 계보를 다룬다.

또 하나의 설명 노선이 토다 격자다. 1967년 토다가 지수형 상호작용을 갖는 사슬이 정확히 적분가능함을 보였는데, 이 퍼텐셜을 3차까지 전개하면 정확히 α-FPUT이 된다. 즉 α-FPUT은 적분가능계인 토다 격자의 작은 섭동이며, 그 보존량들이 근사적으로 살아남아 에너지 흐름을 막는다는 그림이다. 연속체 극한을 거치지 않아도 되는 설명이라 유한 $N$ 에서도 잘 작동한다.

5. 준상태, 그리고 결국은 오는 열화[편집]

솔리톤 그림만 보면 “FPUT 사슬은 열평형에 절대 안 간다”로 오해하기 쉬운데, 사실이 아니다. 훨씬 긴 시간을 돌리면 결국 등분배로 간다. 현대적 그림은 두 단계다.

빠른 단계 — 준상태(metastable state). 저차 모드 몇 개가 빠르게 에너지를 나눠 가지며 지수적으로 감소하는 모드 에너지 분포를 만든다. 이 “자연 패킷”이 겉보기에 안정적으로 유지되는 구간이 FPUT이 본 재귀 구간이다.
느린 단계 — 진짜 등분배. 준상태에서 고차 모드로 에너지가 아주 느리게 새어 나가 결국 등분배에 도달한다. 도달 시간은 에너지 밀도 $\varepsilon = E/N$ 이 작을수록 멱함수적으로 발산한다. 정확한 지수는 모형(α인지 β인지)과 초기 조건에 따라 달라서 지금도 정밀화가 진행 중인 주제이며, α 모형에서는 최저차 공명이 3파나 4파가 아니라 6파라서 열화가 특히 느리다는 것이 파동 난류 이론으로 정리됐다.

에너지를 올리면 그림이 달라진다. 이즈라일레프와 치리코프는 1966년에 에너지가 어떤 문턱을 넘으면 공명이 겹쳐 등분배가 빠르게 일어난다고 지적했다. KAM 정리의 언어로는 낮은 에너지에서 살아남던 불변 토러스가 문턱 위에서 무너지는 것이고, 표준 사상의 공명 겹침 판정과 같은 계산이다. 그 뒤 강한 확률화 문턱(strong stochasticity threshold)이라는 이름으로, 에너지 밀도가 어떤 값을 넘으면 최대 랴푸노프 지수의 에너지 의존성이 멱지수를 바꾸며 계의 성격이 정성적으로 달라진다는 것이 확인됐다. 문턱 아래에서는 카오스가 있어도 위상공간의 좁은 영역에 갇혀 사실상 등분배가 안 일어나고, 위에서는 곧바로 뒤섞인다.

여기서 자유도 수의 문제가 걸린다. KAM 토러스가 차지하는 측도는 자유도가 커질수록 급격히 줄어들므로, $N \to \infty$ 극한에서 위 그림이 그대로 살아남는지는 자명하지 않다. FPUT 문제가 60년 넘게 완전히 닫히지 않은 이유가 이것이다 — 유한 $N$ 의 관측을 열역학적 극한으로 어떻게 외삽할 것인가.²

6. 전산물리학의 출발점으로서[편집]

FPUT이 특별한 이유는 결과 자체보다 방법론에 있다. 이 계산 이전까지 물리학의 도구는 이론과 실험 둘이었다. FPUT은 종이로는 접근이 안 되고 실험실에서도 만들 수 없는 계를 컴퓨터 안에 세워 놓고 관찰한 뒤, 거기서 나온 관측이 기존 이론을 반박하게 만들었다. 시뮬레이션을 “계산 도구”가 아니라 발견의 수단으로 쓴 첫 사례이며, 계산물리와 전산과학이 독립 분야로 자리 잡는 계기가 됐다.

몇 가지 실무적 교훈도 남겼다.

예상 밖 결과가 나왔을 때 코드부터 의심하는 것은 옳지만, 거기서 멈추면 안 된다. FPUT 팀은 실제로 오랫동안 버그를 찾았다. 결과를 믿을 수 있게 만든 것은 심증이 아니라 반복·조건 변경·보존량 확인이었고, 이건 오늘날 검증 및 확인이라 부르는 절차 그 자체다.
적분기 선택이 결론을 좌우한다. 이 계는 해밀토니안 역학 계이고 관심 대상이 $10^6$ 주기 이상의 장시간 거동이라, 에너지가 단조 표류하는 적분기를 쓰면 “느린 열화”가 물리인지 수치 오차인지 구분할 수 없게 된다. FPUT 재현 계산이 베를레 적분 같은 심플렉틱 적분기를 쓰는 이유다. MANIAC 시절에는 이런 개념이 없었고, 오히려 그래서 후대 재현 연구가 필요했다.
“수렴했다”가 “평형에 도달했다”가 아니다. 준상태는 관측량이 안정적으로 보이는 구간을 만든다. 여기서 평균을 재면 오차 막대까지 예쁘게 나오지만 값은 평형값이 아니다. 분자 시뮬레이션에서 에르고딕 가설이 조용히 깨지는 방식과 똑같은 함정이다.³

FPUT은 실패한 실험으로 시작해서 하나의 분야를 여러 개 낳았다. 등분배를 확인하지 못한 대신 솔리톤·가적분계·카오스 전이·비선형 격자 열전도를 얻었으니, 투자 대비 수익률로는 아마 계산물리 역사상 최고일 것이다.⁴

7. 관련 문서[편집]

페르미는 1954년 11월에 세상을 떠나 보고서가 정식으로 나오는 것을 보지 못했다. 보고서 LA-1940은 1955년에 나왔고 널리 읽힌 것은 1965년 페르미 논문집에 실리면서였다. 즉 이 유명한 계산은 저자 한 명이 사망한 뒤에 유명해졌고, 실제 프로그래머는 반세기 뒤에 이름이 붙었다. 전산물리학의 출발점치고는 사연이 많다. ↩
그래서 이 분야 논문을 볼 때는 $N$ , 에너지 밀도 $\varepsilon = E/N$ , 적분 시간, 초기 조건(어느 모드에 넣었는지) 네 가지를 먼저 확인하는 게 습관이 된다. 이 중 하나만 달라도 “등분배가 일어난다/안 일어난다”가 뒤집히기 때문에, 조건을 안 밝힌 주장은 그냥 넘기면 된다. ↩
요즘 말로 하면 “로스가 평평해진 것과 학습이 끝난 것은 다르다”에 가깝다. 관측량이 평평한 구간은 평형의 증거가 아니라 시간 규모가 갈라졌다는 증거일 수 있다. ↩
덧붙여 이 계산은 계산 자원 대비 수익률에서도 압도적이다. MANIAC I의 성능은 지금 스마트폰의 수십억 분의 1 수준이었고, 계산 결과는 손으로 그래프를 그려 확인했다. 지금 노트북으로 몇 초면 재현되는 계산이 물리학의 한 챕터를 열었다는 사실은, 중요한 건 FLOPS가 아니라 무엇을 계산할지 고르는 안목이라는 오래된 교훈으로 남았다. ↩