格林-久保(Green-Kubo)公式：深潛淵靜，回響自呈

當我們將一把金屬勺子浸入熱湯，熱量便從湯沿著勺子傳到指尖，帶來灼熱感。這一簡單現象背后，卻著深刻的物理圖景：熱量究竟如何在物質中傳導？兩個世紀前，約瑟夫·傅里葉給出了一個優美的答案：熱流正比于溫度梯度，其中的比例系數 κ，被人們之為熱導率。

故事的開端總是清晰明了，但其魅力往往在于，探其全貌卻需要數代人的努力和堅持。κ這個符號，測量它或許容易，但理解其微觀圖像卻極其困難。物體受熱后溫度的變化，可能是其內部無數原子、分子或聲子集體舞蹈的宏觀體現。若想仔細探究這舞蹈的藝術內涵，或許不能只觀看臺前的演出，而需深入幕后，觀察它們日復一日的排練——研究平衡態下物體內部自發的微小漲落。我們尋找的答案，或許就藏匿于其中，藏在那連接平衡與非平衡狀態的漲落-耗散定理(Fluctuation Dissipation Theorem)之中，并最終凝練為精妙的數學表達——格林-久保公式。

格林-久保公式是統計物理中的一座關鍵橋梁，它揭示了一個深刻的原理：系統受到微小擾動（如溫度梯度）時產生的線性響應，與其在平衡態下自發隨機漲落的統計特性，在本質上相通。它巧妙地將微觀粒子的隨機運動與宏觀可觀測的輸運現象聯系起來，為理解擴散、電導、粘滯等現象提供了統一的理論框架。

該公式常被提及，但它究竟如何演化而來？又克服了哪些思想障礙？讓我們從歷史與邏輯的雙重脈絡出發，按圖索驥。試試看能不能從微觀動力學出發，抵達計算熱導率 κ的格林-久保公式：

序幕：直覺與洞見

宏觀的耗散（系統趨向平衡）與微觀的漲落（系統自發偏離平衡），看似分屬兩個世界，但它們之間的聯系早在19世紀末便已初現端倪。

對現象背后統一規律的追尋，從來都是物理學不竭的動力。19世紀末，電化學工業的興起推動了對溶液中離子運動的理解。瓦爾特·能斯特(Walther Hermann Nernst)在研究離子擴散時，將其解釋為滲透壓梯度與粘性阻力之間的平衡，并預見了擴散系數與遷移率之間應該存在關聯（即后來的愛因斯坦關系）。雖然能斯特當時更傾向于從熱力學進行解釋，但他的工作觸及了一個核心：驅動宏觀擴散的“力”與微觀粒子受到的“阻力”，似乎同根同源。

1905年，阿爾伯特·愛因斯坦(Albert Einstein)發表了關于布朗運動的開創性論文。他洞察到，懸浮微粒的無規則運動遠非簡單的力學問題。在能斯特思想的啟發下，他將花粉顆粒的隨機熱運動（漲落）與其在流體中受粘性阻力時的定向遷移能力（耗散）聯系起來，嚴格推導出著名的愛因斯坦關系：

。其式中 D是擴散系數（描述隨機散布的快慢），μ是遷移率（描述在外力下定向運動的難易）。這是一個里程碑。它首次明確揭示：平衡態下那看似無規的分子熱運動（漲落），其強度竟決定了系統對外界驅動（耗散）的響應強弱。漲落與耗散，如同一枚硬幣的兩面。

階段一：筑底-微觀守恒定律宏觀唯象定律

任何宏大的理論建構都始于堅實的基石。對于輸運理論，這兩塊基石是：微觀守恒定律與宏觀唯象定律。

1. 微觀守恒定律（連續性方程）

守恒定律是物理學不可撼動的支柱。假設系統中存在守恒量（如粒子數、能量），其局域密度 ρ 的變化必然嚴格滿足：

此即連續性方程。它表明局部守恒量的變化率，完全由流入或流出的“流” Ja的散度決定。這是推導的出發點，精確而普適。

2. 宏觀唯象假設（線性響應）

僅有守恒律，我們尚不知“流”的具體形式。在接近平衡態時，一個可靠的經驗規律是：由梯度驅動的平均“流”，正比于梯度本身，方向相反。這就是線性響應假設：

它連接了“驅動力”（梯度）與“響應”（流）。需要注意：傅里葉定律或菲克定律描述的是穩態非平衡過程。而我們在此處引入統計平均 ???，其深意在于假定：即便在接近平衡的、非穩態的弛豫過程中，平均而言，流與梯度之間仍滿足此線性關系。這是將穩態定律推廣至弛豫過程的關鍵一步。

3. 唯象輸運方程

將唯象關系(2)代入連續性方程(1)，我們得到一個熟悉的方程：

ρ ρ

這正是擴散方程。它描述了宏觀密度分布如何隨時間平滑演化，系數 D即是擴散系數。至此，我們有了描述近平衡宏觀行為的方程，但它仍未直接觸及微觀本質。

階段二：橋梁-漲落-耗散定理

方程(3)描述的是存在外部驅動下的非平衡弛豫。而我們希望了解的，是系統內部平衡態的微觀動力學，這道鴻溝該如何跨越？

答案就是昂薩格假設，其所蘊含的深刻原理，后經嚴格發展成為漲落-耗散定理。它告訴我們：系統對微小外力的線性響應函數，與其在平衡態下自發漲落的時間關聯函數，在數學上完全等價。

4. 核心思想：漲落即響應

愛因斯坦之后，便是站在漲落-耗散定理起點處的是拉爾斯·昂薩格(Lars Onsager)。1925年的一天，剛從挪威化學工程專業的他走進蘇黎世聯邦理工學院德拜(Peter Deby)的辦公室，質疑其的電解液中粒子傳導理論并非正確。德拜對他解釋深感震撼。他安排發表了昂薩格的成果，并邀請他到蘇黎世與自己共事數年。隨后昂薩格前往約翰霍普金斯大學工作。在1931年他了著名的昂薩格倒易關系。

昂薩格的思考更加深刻：如果多種輸運過程（如熱導與電導）耦合在一起，它們的系數矩陣為什么必須是對稱的？他的答案直指微觀世界的根本法則——微觀可逆性

昂薩格做了一個大膽而優美的假設：一個自發產生的微小漲落，其最可能的衰減方式，與系統從一個被刻意制造出來的、同等大小的非平衡狀態開始弛豫的方式，遵循完全相同的宏觀規律。換句話說，無論你是偶然看到湖面泛起漣漪（漲落），還是故意扔下一顆石子（擾動），漣漪散開的方式（耗散）都是一樣的。

這正是我們所需要的橋梁！它意味著，在研究公式(3)所描述的非平衡弛豫時，我們可以合法地將目光轉向平衡態下自發發生的、同等尺度的漲落弛豫。因為根據昂薩格的理論，它們背后的運動方式是一樣的。

從今天任何一本教材來看昂薩格的理論都舉足輕重。但是在當時，盡管昂薩格在電化學領域的研究備受推崇，但倒易關系理論因為其在統計物理上的一解釋漏洞而在十多年間幾乎未引起廣泛關注。20世紀30年代，昂薩格的研究在全球科學出版物中僅有四次被引用。直到德國物理學家約瑟夫·梅克斯納(Josef Meixner)，荷蘭物理學Sybren Ruurds de Groot和亨德里克·格哈德·卡西米爾 (Hendrik Brugt Gerhard Casimir)的研究和分析，其理論才逐漸受人認可。

回到主線，既然描述微觀狀態下量的變化，我們需要借助一種新的數學工具“自關聯函數”。他等于一個信號在不同時刻下乘積的期望。自關聯函數來自于信號科學，最初目的是描述一個信號在兩個不同時刻下的相似程度。這剛好就能衡量漲落狀態下兩個不同時刻下物體內部量的變化程度。

那么我們就進行一個關鍵的替換：在擴散方程(3)中，用平衡態下的密度關聯函數C(x?y,t?t′)，替換掉非平衡的平均密度 ? ?。定義關聯函數為：

ρ ρ

其中下標 0強調這是平衡態平均。替換后我們得到：

這是關鍵的飛躍！ 方程(5)在形式上與(3)相同，但物理內涵已然升華：它不再描述宏觀濃度的擴散，而是描述平衡態下，微觀密度漲落的關聯在時空上是如何傳播和衰減的。擴散系數 D現在同時支配著宏觀梯度的弛豫與微觀漲落的衰減。我們成功地將宏觀輸運系數與平衡態的微觀統計聯系了起來。

第三階段：數學處理——在傅里葉空間簡化問題

方程(5)仍然是一個復雜的偏微分方程。為了求解并提取 D，我們需要請出強大的數學工具。

4. 傅里葉變換

利用系統的空間均勻性，對關聯函數進行空間傅里葉變換。一個美妙的結果是，變換后將微分算符 簡化為乘法因子 。在波矢 k空間方程變為：

方程在“波矢空間”中解耦了。每個波數 k（對應一個特定的空間尺度）的模式獨立演化。它的解是一個簡單的指數衰減：

這意味著，平衡態密度漲落在不同空間尺度上的關聯隨時間呈指數衰減。

5.走向宏觀

輸運系數描述的是系統在宏觀尺度（長波）和長時間下的行為。這對應數學上的雙重極限：波矢 k→0（長波極限，對應宏觀空間尺度）和頻率 ω→0（低頻極限，對應宏觀時間尺度）。對(7)式進行時間傅里葉變換，并在 k→0,ω→0的極限下分析，可以得到擴散系數的一個表達式：

雖然成功推導了擴散系數D依賴于“密度關聯函數”，但我們在微觀模擬（如分子動力學）中有一類更容易直接獲取的物理量，是粒子的速度和受力，即“流”的信息。

第四階段：關鍵轉化——從“密度關聯”到“流關聯”

審視公式(8)，它使用的是密度關聯函數。但驅動輸運的是“流”。在計算上，流（與粒子速度、力直接相關）比密度更容易從微觀模擬中獲得。因此，我們需要進行最后也是最關鍵的一步轉化：將密度關聯轉化為流關聯。

昂薩格的理論確立了漲落與耗散相關的原則，但將其轉化為一個具體、可計算的公式，我們還需要繼續等待。

20世紀50年代，所有人等來了公式的名字。

格林(Melville Green)

梅爾維爾·格林（Melville Green），他是普林斯頓大學尤金·維格納（Eugen Wigner）和埃利奧特·蒙特羅爾（Elliott Montroll）的博士生。這位格林不應與赫伯特·格林（Herbert Green，玻恩的學生）或理查德·格林（Richard Greene，卡倫的學生）相混淆。格林大量借鑒了布朗運動理論的類比，雖然他收到了多位物理學家的啟發，尤其是柯克伍德（Kirkwood）。他的方法仍是原創性的，并得出了耗散系數與漲落量的相關矩之間的更為普遍的關系。簡言之，他假設宏觀可觀察得參數 的統計力學漲落具有馬爾可夫性。基于這一假設，他推導出了系統從初始時刻的一組參數值演化到稍后時刻另一組參數值的概率。擴散（耗散）系數 因而表現為與轉移概率的二階矩相關聯，而一些柯克伍德式的操作則將這些矩轉化為自相關函數。最終結果如下：

久保亮五(Ryogo Kubo)

久保亮五(Ryogo Kubo)是20世紀日本極具影響力的理論物理學家，在統計力學和凝聚態物理領域做出了開創性貢獻。

1920年2月15日出生于東京，父親久保天隨（Tenzui Kubo）是著名的漢學家。其父親于1929年赴臺北帝國大學任教，全家遷居臺灣，直至1934年父親去世后才返回東京。受家庭影響，久保最初對文學和哲學感興趣。然而，在1936年進入第一高等學校后，他選擇了理科課程，并決心專攻物理學。

久保意識到，系統的響應可以與其內部自發漲落的時間關聯性聯系起來。他借鑒了高橋秀俊(Hidetosi Takahasi)等人的方法，運用量子力學微擾理論，得到了一個極其優美的表達式——這就是以他命名的久保公式（Kubo formula）。這個公式是昂薩格思想的定量化和普遍化。它不再限于特定模型，而是給出了一個“配方”：要計算某種輸運系數，就去計算相應“流”算符在平衡態下的時間自相關函數的積分。

久保公式的經典極限，正是我們通過數學推導即將抵達的終點。它告訴我們，從“密度關聯”到“流關聯”的轉化，并非只是數學技巧，而是將問題引向最本質、最可計算的物理量的必然一步。因為“流”直接對應著微觀粒子的速度和相互作用力，是分子動力學模擬可以直接“看到”和“記錄”的量。

7. 利用連續性方程連接二者

還記得我們的出發點——連續性方程(1)嗎？它將密度的時間導數與流聯系起來，這正是我們需要的轉換器。

我們對密度關聯函數C求二階時間導數。根據連續性方程，密度的一階時間導數就是流的散度。經過一番推導（主要涉及傅里葉變換和散度運算），我們得到：

這里 Ji(k,t)是流在傅里葉空間的分量。這個等式至關重要：密度關聯的衰減速率，直接由流-流關聯函數決定。

8. 抵達廣義格林-久保公式

將(9)式與之前的結果結合，經過一些必要的數學變換，在此不多贅述，并在 k→0極限下取流的主要部分（即系統的總流）。假設系統是各向同性的，我們最終得到簡潔而強大的廣義Green-Kubo公式：

這個公式是一座完美的橋梁。左邊的 D是宏觀輸運系數；右邊的積分，是微觀平衡態下“總流”自相關函數的時間積分。它告訴我們：輸運能力的強弱，取決于系統內部“流動”的方向性能保持多久的記憶。如果流動方向雜亂無章、瞬間即變（關聯衰減快），積分值就小，輸運能力弱；如果流動方向能長時間保持一致（關聯衰減慢），積分值就大，輸運能力強。值得注意的是，公式 (8) 顯示擴散系數 D越大，密度漲落關聯衰減越快；而 Green-Kubo 公式則告訴我們，D的大小由流自相關函數的積分決定——流關聯衰減越慢，D越大。這并不矛盾，因為前者描述的是密度不均勻性的平滑速度，后者揭示的是微觀流動方向的記憶效應。實際上，擴散系數正是連接這兩者的橋梁：微觀流動方向越持久（流關聯衰減慢），宏觀上物質擴散越快（密度關聯衰減快）

第五階段：終點抵達——應用于熱導率

現在，我們將這個普適的公式應用到最初關心的熱傳導問題上。

9. 具體化到熱傳導中的參數

• 守恒量：內能（能量）。

• 對應的流：能量流密度，其總和即總能量流 。

• 關鍵系數 C：對于能量密度，其平衡態漲落的強度 與系統的熱力學性質直接相關。統計力學給出一個關鍵結果：能量漲落的方差

其中 是定容熱容。由此可以推導出，對于能量， 。

將能量對應的 和 代入廣義公式(11)，并利用熱容與體積的關系進行化簡，我們便抵達了最終目的地：

最終詮釋：熱導率，等于系統在平衡態時，其微觀能量流方向“記憶效應”的持久程度（時間積分）。我們無需真正給材料兩端加上溫差（非平衡模擬），只需讓它“靜置”在平衡態，記錄其內部能量流的自然漲落，分析這種漲落的相關性隨時間如何衰減，就能精確預測它的導熱能力。

結語：若要知著，必先見微

回顧這段跨越世紀的旅程，我們從宏觀的傅里葉定律出發，深入微觀的守恒方程；憑借漲落-耗散定理這座由愛因斯坦、昂薩格等人構想，并由格林、久保等人精心構筑的橋梁，我們跨越了非平衡與平衡的鴻溝；再通過數學工具，最終將宏觀的熱導率 κ，錨定在了微觀能量流漲落的時間記憶之上。

格林-久保公式的建立，標志著一個范式的成熟。它不再依賴于玻爾茲曼方程等特定模型與近似，而是將計算線性輸運性質的所有復雜性，封裝在了平衡態統計力學之中——正如平衡態性質歸于配分函數的計算，線性響應性質則歸于平衡態時間關聯函數的計算。

格林-久保公式不僅僅是一個強大的計算工具，更深刻地揭示了自然內在的一種簡潔與經濟：那最終的答案，從不聲張，也無心炫耀。千百年來，它只是安然盤踞在平靜水面之下，淺唱低吟。若非心細如發之人，難以覺察其蹤跡；若非意志堅定之人，也難有勇氣循著這微弱的線索，潛入深處，去一窺其完整的樣貌。

這場從后臺排練（平衡漲落）預測前臺演出（非平衡輸運）的智慧之旅，正是理論物理思想之深邃與力量的最佳體現。

參考文獻

物理部分：

Livi R, Politi P. Nonequilibrium Statistical Physics: A Modern Perspective. 2nd ed. Cambridge University Press; 2025.

歷史部分

Darrigol, O. A history of the relation between fluctuation and dissipation. EPJ H48, 10 (2023)

Kono, H. Ryogo Kubo in his formative years as a physicist. EPJ H45, 175–204 (2020).

李保文課題組公眾號｜rezhi2021

文章轉載自“熱知”微信公眾號