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

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

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

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

序幕：直覺與洞見

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

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

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

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

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

任何宏大的理論建構(gòu)都始于堅(jiān)實(shí)的基石。對于輸運(yùn)理論，這兩塊基石是：微觀守恒定律與宏觀唯象定律。

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

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

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

2. 宏觀唯象假設(shè)（線性響應(yīng)）

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

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

3. 唯象輸運(yùn)方程

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

ρ ρ

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

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

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

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

4. 核心思想：漲落即響應(yīng)

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

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

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

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

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

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

那么我們就進(jìn)行一個(gè)關(guān)鍵的替換：在擴(kuò)散方程(3)中，用平衡態(tài)下的密度關(guān)聯(lián)函數(shù)C(x?y,t?t′)，替換掉非平衡的平均密度 ? ?。定義關(guān)聯(lián)函數(shù)為：

ρ ρ

其中下標(biāo) 0強(qiáng)調(diào)這是平衡態(tài)平均。替換后我們得到：

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

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

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

4. 傅里葉變換

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

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

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

5.走向宏觀

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

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

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

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

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

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

格林(Melville Green)

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

久保亮五(Ryogo Kubo)

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

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

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

久保公式的經(jīng)典極限，正是我們通過數(shù)學(xué)推導(dǎo)即將抵達(dá)的終點(diǎn)。它告訴我們，從“密度關(guān)聯(lián)”到“流關(guān)聯(lián)”的轉(zhuǎn)化，并非只是數(shù)學(xué)技巧，而是將問題引向最本質(zhì)、最可計(jì)算的物理量的必然一步。因?yàn)椤傲鳌敝苯訉?yīng)著微觀粒子的速度和相互作用力，是分子動(dòng)力學(xué)模擬可以直接“看到”和“記錄”的量。

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

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

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

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

8. 抵達(dá)廣義格林-久保公式

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

這個(gè)公式是一座完美的橋梁。左邊的 D是宏觀輸運(yùn)系數(shù)；右邊的積分，是微觀平衡態(tài)下“總流”自相關(guān)函數(shù)的時(shí)間積分。它告訴我們：輸運(yùn)能力的強(qiáng)弱，取決于系統(tǒng)內(nèi)部“流動(dòng)”的方向性能保持多久的記憶。如果流動(dòng)方向雜亂無章、瞬間即變（關(guān)聯(lián)衰減快），積分值就小，輸運(yùn)能力弱；如果流動(dòng)方向能長時(shí)間保持一致（關(guān)聯(lián)衰減慢），積分值就大，輸運(yùn)能力強(qiáng)。值得注意的是，公式 (8) 顯示擴(kuò)散系數(shù) D越大，密度漲落關(guān)聯(lián)衰減越快；而 Green-Kubo 公式則告訴我們，D的大小由流自相關(guān)函數(shù)的積分決定——流關(guān)聯(lián)衰減越慢，D越大。這并不矛盾，因?yàn)榍罢呙枋龅氖敲芏炔痪鶆蛐缘钠交俣?，后者揭示的是微觀流動(dòng)方向的記憶效應(yīng)。實(shí)際上，擴(kuò)散系數(shù)正是連接這兩者的橋梁：微觀流動(dòng)方向越持久（流關(guān)聯(lián)衰減慢），宏觀上物質(zhì)擴(kuò)散越快（密度關(guān)聯(lián)衰減快）

第五階段：終點(diǎn)抵達(dá)——應(yīng)用于熱導(dǎo)率

現(xiàn)在，我們將這個(gè)普適的公式應(yīng)用到最初關(guān)心的熱傳導(dǎo)問題上。

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

• 守恒量：內(nèi)能（能量）。

• 對應(yīng)的流：能量流密度，其總和即總能量流 。

• 關(guān)鍵系數(shù) C：對于能量密度，其平衡態(tài)漲落的強(qiáng)度 與系統(tǒng)的熱力學(xué)性質(zhì)直接相關(guān)。統(tǒng)計(jì)力學(xué)給出一個(gè)關(guān)鍵結(jié)果：能量漲落的方差

其中 是定容熱容。由此可以推導(dǎo)出，對于能量， 。

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

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

結(jié)語：若要知著，必先見微

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

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

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

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

參考文獻(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).

李保文課題組公眾號｜rezhi2021

文章轉(zhuǎn)載自“熱知”微信公眾號