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

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

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

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

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

序幕：直覺與洞見

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

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

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

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

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

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

1. 微觀守恒定律（連續(xù)性方程）

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

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

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

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

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

3. 唯象輸運方程

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

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

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

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

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

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

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

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

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

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

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

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

那么我們就進行一個關鍵的替換：在擴散方程(3)中，用平衡態(tài)下的密度關聯(lián)函數(shù) ，替換掉非平衡的平均密度 。定義關聯(lián)函數(shù)為：

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

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

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

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

5.傅里葉變換

利用系統(tǒng)的空間均勻性，對關聯(lián)函數(shù)進行空間傅里葉變換。一個美妙的結果是，變換后將微分算符 簡化為乘法因子 。在波矢 空間方程變?yōu)椋?/p>

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

這意味著，平衡態(tài)密度漲落在不同空間尺度上的關聯(lián)隨時間呈指數(shù)衰減。

6.走向宏觀

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

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

第四階段：關鍵轉化——從“密度關聯(lián)”到“流關聯(lián)”

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

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

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

格林(Melville Green)

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

久保亮五(Ryogo Kubo)

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

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

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

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

7. 利用連續(xù)性方程連接二者

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

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

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

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

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

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

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

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

9. 具體化到熱傳導中的參數(shù)

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

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

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

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

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

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

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

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

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

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

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

參考文獻

物理部分：

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).

來源：熱知

編輯：LogicMoriaty

轉載內容僅代表作者觀點

不代表中科院物理所立場

如需轉載請聯(lián)系原公眾號