1955 - Fermi–Pasta–Ulam–Tsingou 實驗——電腦上的第一個意外#

案件摘要#

1955 年夏天，Fermi、Pasta、Ulam（與程式設計師 Mary Tsingou）在 Los Alamos 的 MANIAC 電腦上模擬一維非線性晶格：64 顆質點、線性彈簧 + 微弱二次非線性。他們預期驗證能量均分定理（統計力學基石）：非線性會把能量均勻攪散。結果完全相反：能量在幾個「模式」之間近乎週期性地來回轉移，最後幾乎全部回到最初的低階模式（復發現象）。統計力學的預期在電腦上破產——這是計算物理史上第一個「意外發現」，後稱 FPUT 問題。

前因#

能量均分定理：統計力學預言非線性耦合下能量終將均分到所有模式——這是 Boltzmann–Gibbs 架構的核心。
固體比熱的懸案：若均分成立，固體比熱應為常數（Dulong–Petit），但實驗在低溫偏離——Debye 聲子理論是補救。FPUT 想用電腦直接看「非線性如何導致均分」。
MANIAC 電腦（1952）：von Neumann 架構的早期電腦，第一次讓「大量重複計算」可行。
Fermi 的直覺：大師晚年的最後一個大型研究——「我們在電腦上會看到統計力學的誕生過程」。

線索與推理#

模型：非線性晶格（β-FPUT 晶格）#

\[ m\ddot q_j = q_{j+1} - 2q_j + q_{j-1} + \beta\left[(q_{j+1}-q_j)^2 - (q_j - q_{j-1})^2\right] \]

線性部分（ \(\beta = 0\) ）是純聲子系統： \(N\) 個獨立正交模式、完全可積、能量永不分均。非線性 \(\beta\) 是「攪拌器」。預期： \(\beta \ne 0\) ⇒ 模式間耦合 ⇒ 均分。

意外的結果#

以正弦初始條件（能量集中在最低模式 \(k=1\) ）起跑，觀察各模式能量 \(E_k(t)\) ：

時間：  0 ────→ E1 幾乎 100%
時間： ~1 千周期 ──→ E1 流失到 E2..E5（但只到前 5 個模式！）
時間： ~2 千周期 ──→ 能量幾乎全部「回流」到 E1（復發！）

能量沒有均分到 64 個模式，且系統近乎週期。 均分定理在這個系統上「失蹤」。

侦查線索 1：KAM 環面的解釋（1960s）#

非線性 \(\beta\) 很弱 ⇒ 系統是「近可積」（KAM 定理的情境，見 1954-KAM定理.md）：多數不變環面存活，軌跡被環面約束在低維流形上——能量被「籠」在低階模式附近，均分被阻擋。 FPUT 復發 ≈ 軌跡在不變環面上的近週期運動。意外有了理論解釋。

侦查線索 2：孤波（soliton）的解釋（1965）#

Zabusky 與 Kruskal 把 FPUT 晶格取長波極限，得到 KdV 方程：

\[ \frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + \delta^2\frac{\partial^3 u}{\partial x^3} = 0 \]

KdV 有孤波解（soliton）：形狀不變、碰撞後恢復的非線性波。 FPUT 的「近週期性」= 一群孤波的碰撞與復位。孤子動力學（完全可積的非線性系統）誕生。

侦查線索 3：何時才會均分？#

後續研究（1970s–至今）釐清： | 情況 | 行為 | |---|---| | 弱非線性 + 低能量 | KAM 環面主導，能量受限（FPUT 復發）| | 更長時間 / 更高能量 | 混沌層擴張，緩慢均分（Chirikov 擴散）| | 強非線性 | 快速均分，統計力學恢復 |

均分沒有破產，只是時間尺度遠超預期——統計力學在更長時標上勝訴。

結案#

計算物理成為「第三種偵探方法」：理論與實驗之外，電腦模擬能「發現」意外。
三大連鎖：KAM 定理獲得具體樣本、孤子與可積非線性系統誕生（KdV、sine-Gordon、Toda 晶格）、混沌的數值研究正規化。
歷史註腳：Mary Tsingou 實際寫了程式並做了大部分計算，但 1955 年報告未署其名——2020 年代起學界補正，改稱 FPUT（含 Tsingou）。
Fermi 於 1954 年底去世，未能看到這個「意外」如何改寫非線性物理——大師留下的最後一問，答案比預期深刻得多。