2013 年多尺度模擬諾貝爾獎：三線匯流的加冕之案#

副標題：Karplus、Levitt、Warshel 用三十七年證明—— 化學可以在電腦裡分層重演。年表定位：第四幕終章，為 QM/MM、BPTI、MD 三線匯流加冕，承前啟後。

案發現場#

2013 年 10 月 9 日，斯德哥爾摩的電話打向三地：哈佛的 Martin Karplus、史丹佛的 Michael Levitt、南加大（USC）的 Arieh Warshel。瑞典皇家科學院宣布，當年諾貝爾化學獎授予他們，理由是「發展多尺度模型研究複雜化學體系」。

外人看來，這像一樁冷案重啟：得獎工作多在七〇年代完成，為何等了三十七年？內行人卻知道，這正是計算模擬最漫長的臥底行動—— 把量子與經典、把電子與蛋白質、把方程與電腦，一層層縫進同一個案卷。

回到最初的案發現場，嫌犯是尺度本身：

第一，量子無力。薛丁格方程原則上主宰化學，但對幾千個原子的酵素求解，計算量是天文數字。全量子處理蛋白質，等於為了審一間屋子而搜查整座城市。

第二，經典無能。分子力場跑得動大蛋白，卻不懂電子重排。鍵斷鍵結之處，正是酵素最關鍵的犯罪現場，經典力場恰恰在那裡失明。

第三，動態無蹤。晶體結構是靜態照片，而催化是電影。沒有 MD 的時間軸，再精確的能量也只是死證物。

三位偵探在 1960 至 1970 年代各執一端：Karplus 懂量子化學與動力學，Levitt 懂粗粒化與結構，Warshel 懂半經驗量子與酶學。他們的匯流，花了整整十年。

1976 年，Warshel 與 Levitt 發表開創性論文：把溶菌酶（lysozyme）的反應中心用量子力學處理，外圍蛋白與溶劑用經典力學處理。這就是 QM/MM（quantum mechanics / molecular mechanics）。

其能量分割思想可寫為加法式：

\[ E_{total}=E_{QM}(QM)+E_{MM}(MM)+E_{QM-MM}(int) \]

其中 \(E_{QM}\) 是量子區能量， \(E_{MM}\) 是經典區能量， \(E_{QM-MM}\) 是兩區耦合（含靜電嵌入與邊界處理）。更精細的減法式（ONIOM 風格）則寫為：

\[ E_{total}=E_{MM}^{full}+E_{QM}^{inner}-E_{MM}^{inner} \]

其中上標 \(full\) 代表全體系， \(inner\) 代表內層量子區。兩式精神一致：好鋼用在刀刃上，量子只算反應中心。

邊界是全案最棘手的證物處理：共價鍵被 QM/MM 界面切斷時，需用氫 link 原子或 frozen orbital 封口，否則懸鍵會偽造電荷。下表是三種邊界手法：

QM/MM 的辦案哲學是分層真相：反應鍵處用量子明察秋毫，蛋白骨架用經典維持秩序。

1977 年，McCammon、Gelin 與 Karplus 發表牛胰蛋白酶抑制劑（BPTI）的真空 MD，約 9 皮秒。這是第一次讓折疊蛋白在電腦裡動起來，配上 Verlet 積分與經驗力場，原子不再是結晶照片裡的蠟像。

MD 的運動方程仍是牛頓式。對第 \(i\) 個原子有 \(F_i=m_ia_i\) 的關係，力來自勢能梯度：

\[ m_i\frac{d^2r_i}{dt^2}=-\nabla_{r_i}U(r_1,\cdots,r_N) \]

其中 \(m_i\) 是質量， \(r_i\) 是座標， \(U\) 是力場勢能。9 皮秒在今天看來是眨眼，但在 1977 年是開天闢地：它證明蛋白質的呼吸、側鏈翻轉、氫鍵開合，都可以被逐步積分出來。

Karplus 陣營在八〇年代持續添購軍火：CHARMM 力場（1983）、自由能微擾、簡正模分析、蛋白質內部水分子模擬。Levitt 則另闢蹊徑：粗粒化、摺疊路徑的簡化模型、DNA 與大複合體建模。兩條線一精一簡，恰好互補。

下表是三線匯流的時間軸：

諾獎公告特別點出三人的互補：Karplus 把量子化學帶進動力學，Levitt 把簡化與計算帶進結構生物學，Warshel 把半經驗量子帶進酶催化。三人1970年代在巴黎、劍橋、哈佛、魏茨曼之間交錯合作，師承與共事關係盤根錯節。

值得玩味的是，2002 年的元動力學、2007 年的神經勢能，乃至 2020 年的 AlphaFold，全都站在這個分層思想的肩膀上：先切分尺度，再各顯神通。這正是評審口中「多尺度模型」的真意—— 不是一個方程，而是一種斷案體制。

2013 年的加冕，結的不是一樁案，而是整整三幕：從薛丁格方程到力場，從力場到 MD，從 MD 到 QM/MM，再到後來的增強採樣與機器學習勢能。

它的遺產有四：

其一，酵素機制的計算可解。溶菌酶、光合反應中心、視紫紅質、P450，催化循環的每一步都可被 QM/MM 解剖。實驗 ’看見什麼’，計算回答 ’為何如此’。

其二，藥物設計的理性化。QM/MM 精修對接、自由能計算排先後，讓先導化合物優化不再全靠碰運氣。

其三，方法論的正統化。多尺度從旁門左道變成顯學，ONIOM、QM/MM-MD、polarizable embedding、ML/MM 相繼登場。

其四，冷案重啟的啟示。1976 年的論文在 2013 年得獎，說明計算方法的價值需要算力追上來才能被看見。AlphaFold 的 2024 年諾獎，不過是同一劇本的續集。

結案語：兇手從來不是單一尺度，而是尺度之間的勾結。三位偵探沒有發明新的物理，他們發明了讓不同物理同庭作證的法庭。

關鍵公式一：QM/MM 加法式 \(E_{total}\) ，三項分工見偵查過程。
關鍵公式二：ONIOM 減法式 \(E_{total}\) ，適配多層嵌套。
關鍵公式三：牛頓 MD 方程 \(m_i\) 與 \(r_i\) 之關係，BPTI 心跳的起點。
時間軸表：1976 至 2013 三線匯流，見上表。
工具鏈：CHARMM、AMBER、GROMACS、Gaussian、TURBOMOLE、CP2K、NAMD，QM/MM 接口以 ChemShell、Cobramm 黏合。
辦案心法：先分層，再計算。量子區寧小勿大，經典區寧大勿糙，邊界寧簡勿漏。
延伸卷宗：前案是 2007-BehlerParrinello神經勢能.md，後續第五幕見 AlphaFold 系列；源頭卷宗為 1976 年 QM/MM 與 1977 年 BPTI。