統(tǒng)計(jì)物理中的一道百年小謎題及其破解

新理論不僅能給出實(shí)驗(yàn)已經(jīng)確立的結(jié)果，還能預(yù)言新的結(jié)果，而這一結(jié)果在物理上之合理之巧妙，只能歸結(jié)于大自然的安排。

撰文 | 劉全慧（理論物理博士，湖南大學(xué)物理與微電子科學(xué)學(xué)院教授）

物理難題和數(shù)學(xué)難題有一個(gè)很大的不同，世紀(jì)“高齡”的數(shù)學(xué)難題屢見(jiàn)不鮮，例如，費(fèi)馬大定理從提出到證明相隔了358年；哥德巴赫猜想已經(jīng)280年。有趣的是，數(shù)學(xué)的難題，每解決一個(gè)就相當(dāng)于殺掉了下金蛋的鵝。物理學(xué)中也有百年難題，每解決一次就相當(dāng)于發(fā)現(xiàn)了一只下金蛋的鵝。這些問(wèn)題小到冰面為什么打滑，大到宇宙的結(jié)構(gòu)及其演化；復(fù)雜如氣候變遷，具體如飛鴿回巢，等等。本文的有限尺寸效應(yīng)難題，寫在物理教材里，歷史上不斷有研究論文發(fā)表，是一個(gè)介乎數(shù)學(xué)和物理之間的小難題。

一

有限尺寸效應(yīng)難題

在所有的統(tǒng)計(jì)分布中，求最概然的一個(gè)，是統(tǒng)計(jì)學(xué)的入門知識(shí)。不過(guò)，統(tǒng)計(jì)學(xué)關(guān)心大量的樣本。稀罕的事例，不在統(tǒng)計(jì)學(xué)的范疇內(nèi)。統(tǒng)計(jì)物理也求最概然分布，但是，一個(gè)系統(tǒng)常常會(huì)和熱庫(kù)、粒子庫(kù)接觸，系統(tǒng)中就可能只有一個(gè)粒子。換言之，在統(tǒng)計(jì)物理中，即使出現(xiàn)一個(gè)粒子，也會(huì)出現(xiàn)分布。單原子熱機(jī)是近些年的一個(gè)研究熱點(diǎn)，處理的就是只有一個(gè)粒子的統(tǒng)計(jì)分布。通常的熱力學(xué)處理的是大量粒子的系統(tǒng)，數(shù)學(xué)上處理為無(wú)限多個(gè)粒子，同時(shí)系統(tǒng)體積也要取為無(wú)限大，而粒子數(shù)的密度不變，就是所謂的熱力學(xué)極限。在熱力學(xué)極限的另外一極，就是小系統(tǒng)，或者少粒子系統(tǒng)，或者有限尺寸系統(tǒng)。這里會(huì)出現(xiàn)的新效應(yīng)，可以稱之為有限尺寸效應(yīng)，或者少粒子效應(yīng)。

但是，超過(guò)一個(gè)半世紀(jì)的時(shí)間里，數(shù)以十萬(wàn)計(jì)的統(tǒng)計(jì)物理的學(xué)者們，求最概然分布，一直被如下問(wèn)題所困擾：相關(guān)理論的適用性要求粒子數(shù)很大，只有滿足這一條件才能借助成熟的數(shù)學(xué)工具給出統(tǒng)計(jì)分布。當(dāng)然，當(dāng)粒子數(shù)很大的時(shí)候，這些結(jié)果也都是正確的，受到了實(shí)驗(yàn)的嚴(yán)格檢驗(yàn)。那么，當(dāng)粒子數(shù)很少的時(shí)候，有無(wú)分布? 分布的形式如何？數(shù)學(xué)上，這一問(wèn)題歸結(jié)于如何處理一個(gè)變量x連乘積x!的對(duì)數(shù)lnx!及其微積分，其中變量x可以理解為粒子數(shù)，x=0,1,2,…。例如x在1和10之間，就是所謂的少粒子系統(tǒng)。

在統(tǒng)計(jì)物理中，常常使用lnx!的如下斯特林公式，lnx！≈xlnx-x。這一公式在x的取值非常大的時(shí)候，精度非常之高。但當(dāng)x較小的時(shí)候，例如x在1和10之間取值時(shí)，這個(gè)公式的精度非常差。但是，在物理上，只能這樣取，多一項(xiàng)少一項(xiàng)都不行！關(guān)鍵的原因在于，惟其如此，才能保證一個(gè)熱力學(xué)系統(tǒng)的熵的廣延性。如果利用更為精確的近似，必然破壞廣延性。換言之，精確的斯特林公式，給出的最后的結(jié)果如果是正確的，應(yīng)該是所謂的少粒子效應(yīng)或者有限尺寸效應(yīng)，不過(guò)大概率是錯(cuò)誤的。

當(dāng)x為粒子個(gè)數(shù)的時(shí)候，lnx!是天生的離散函數(shù)。使用近似公式lnx!≈xlnx-x的一個(gè)重要目的在于把lnx!變成連續(xù)函數(shù)進(jìn)行微積分運(yùn)算。但是，當(dāng)x較小的時(shí)候，把lnx!處理為連續(xù)函數(shù)的精度非常差，只能處理為離散函數(shù)，這個(gè)時(shí)候，lnx!的微分dlnx!應(yīng)該用差分Δlnx!來(lái)代替?？墒遣罘枝nx!的定義非常多，例如，前一步差分Δlnx!= ln(x+1)!- lnx!=ln(x+1), 后一步差分Δlnx!= ln(x)!- ln(x-1)!=ln(x)，前兩步差分，后兩步差分，…，中心差分，偏心差分，…等等。定義如此之多，給出的物理結(jié)果互不相同，選擇其中的任意一個(gè)相當(dāng)于引進(jìn)了新的假設(shè)，即使這樣，給出的結(jié)果大概率也是靠不住的。

因此，統(tǒng)計(jì)物理出現(xiàn)lnx!的時(shí)候，在兩個(gè)地方需要進(jìn)行嚴(yán)格的處理，一個(gè)是lnx!精確的表達(dá)式，第二是嚴(yán)格的離散微積分。

二

尋找最穩(wěn)定的解

尋求函數(shù)lnx!的斯特林公式精確表達(dá)，物理學(xué)家玩不過(guò)數(shù)學(xué)家。把數(shù)學(xué)家發(fā)現(xiàn)的結(jié)果搬運(yùn)過(guò)來(lái)就是。所以，問(wèn)題不在是否使用精確的斯特林公式。問(wèn)題必然出現(xiàn)在函數(shù)差分上。經(jīng)過(guò)了很多次探索和失敗之后，我們終于發(fā)現(xiàn)有一條狹窄但是巧妙的解決這一問(wèn)題的途徑[1]，而且根本不需要用到斯特林公式。

統(tǒng)計(jì)物理中求最概然分布，可以理解為一個(gè)求變分的過(guò)程，走在正確的道路上，一定會(huì)碰到一些 “妖魔鬼怪”。下面通過(guò)一個(gè)簡(jiǎn)單的例子來(lái)說(shuō)明。

因此，后差解更加穩(wěn)定。給一個(gè)合理的說(shuō)法，來(lái)自數(shù)學(xué)或者物理都行，只要能挑出后差解就行了。

極少數(shù)研究人員會(huì)走到這里，然后通過(guò)一個(gè)新的物理原則，把最穩(wěn)定的解作為真實(shí)的解，把不真實(shí)的解剔除。

可是，即使到了這里，還沒(méi)有從根本上解決物理問(wèn)題。

三

一點(diǎn)數(shù)學(xué)新花樣：異步差分

經(jīng)過(guò)一點(diǎn)運(yùn)算，新理論給出的結(jié)果分為三部分。第一部分，粒子數(shù)很大的時(shí)候，結(jié)果完全回復(fù)到傳統(tǒng)的結(jié)果。這部分結(jié)果經(jīng)受了實(shí)驗(yàn)的反復(fù)檢驗(yàn)，如果新理論不能給出相同的結(jié)果，理論肯定錯(cuò)了。第二部分，如果系統(tǒng)內(nèi)只有兩個(gè)粒子，例如玻色子或者費(fèi)米子，原來(lái)的分布繼續(xù)有效。這一點(diǎn)，似乎傳統(tǒng)統(tǒng)計(jì)物理中的巨正則分布提示過(guò)這個(gè)結(jié)果。但是，提示不等于確立。巨正則分布給出的結(jié)果，是否適合于只有一兩個(gè)粒子的系統(tǒng)，統(tǒng)計(jì)物理本身給不出判斷。第三部分，如果系統(tǒng)內(nèi)只有一個(gè)粒子，新理論認(rèn)為，粒子的量子性突然消失，所有的粒子都遵守同樣的統(tǒng)計(jì)，即玻爾茲曼統(tǒng)計(jì)。這一個(gè)結(jié)果，是理論首先給出了結(jié)果，后來(lái)才理解。關(guān)于這一點(diǎn)，不妨多說(shuō)幾句。

兩個(gè)粒子是費(fèi)米子或者玻色子，量子力學(xué)認(rèn)為是系統(tǒng)狀態(tài)滿足交換對(duì)稱性的一種后果；而量子場(chǎng)論認(rèn)為，自旋和統(tǒng)計(jì)之間有簡(jiǎn)單的一一對(duì)應(yīng)的關(guān)系。根據(jù)新的理論，如果系統(tǒng)中原本有兩個(gè)玻色子，滿足玻色分布。突然取走其中的一個(gè)，問(wèn)，剩下的一個(gè)還是玻色子嗎? 量子力學(xué)認(rèn)為，由于這個(gè)玻色子不能和其它粒子進(jìn)行交換，再說(shuō)這個(gè)粒子是玻色子或者不是，已經(jīng)失去意義。就在這里，新理論進(jìn)一步預(yù)言，這個(gè)粒子應(yīng)該滿足玻爾茲曼統(tǒng)計(jì)。這個(gè)結(jié)果一度使我非常焦慮，后來(lái)突然理解了，這才是大自然的巧妙且必然的安排。只有一個(gè)粒子的時(shí)候，整個(gè)系統(tǒng)可以作為定域區(qū)域，這個(gè)粒子就是定域粒子，當(dāng)然只能滿足玻爾茲曼統(tǒng)計(jì)。因此，新理論預(yù)言，隨著粒子從兩個(gè)變?yōu)橐粋€(gè)，出現(xiàn)兩個(gè)統(tǒng)計(jì)之間的躍變，這個(gè)時(shí)候一定會(huì)出現(xiàn)新的熱量交換。因此，文章[1]進(jìn)行了如下暢想：新理論能提供實(shí)驗(yàn)可以檢驗(yàn)的結(jié)果。

四

小 結(jié)

很明顯，這里涉及的有限尺寸效應(yīng)問(wèn)題之所以困難，是因?yàn)闆](méi)有恰當(dāng)?shù)臄?shù)學(xué)工具。我是被物理感覺(jué)所指引，覺(jué)得大概應(yīng)該有一個(gè)什么長(zhǎng)相的物理結(jié)果，然后倒逼數(shù)學(xué)，發(fā)現(xiàn)只有引入異步差分才能解決問(wèn)題。也許，在數(shù)學(xué)上有人在其它領(lǐng)域提出過(guò)異步差分的思想，但是從來(lái)沒(méi)有引進(jìn)到物理學(xué)中來(lái)。

異步差分似乎突破了數(shù)學(xué)的陳規(guī)，卻依然在數(shù)學(xué)的框架內(nèi)。異步差分就在我們身邊，距離只有一層薄紙。新理論不僅能給出實(shí)驗(yàn)已經(jīng)確立的結(jié)果，還能預(yù)言新的結(jié)果，而這一結(jié)果在物理上之合理之巧妙，只能歸結(jié)于大自然的安排。

后 記

一輩子寫文章，很多文章寫著寫著就忘記了。投稿之后，也不知道投到了哪個(gè)刊物。甚至有篇不錯(cuò)的文章，因?yàn)橐粚徱呀?jīng)通過(guò)但是沒(méi)有能及時(shí)回復(fù)而被退稿了。這篇文章[1]不同，是我二十余年講授“熱力學(xué)與統(tǒng)計(jì)物理”課程的之后給自己的一個(gè)交代，也是2020年6月份起，華盛頓大學(xué)錢纮教授召集國(guó)際納米熱力學(xué)研討會(huì) (International Seminar Nanothermodynamica series) 以來(lái)，我下決心要交付的一篇作業(yè)。文章在“五·四”青年節(jié)正式刊行，正值納米熱力學(xué)研討會(huì)兩周年之際，祝愿研討會(huì)常研常青。物理學(xué)年刊 (Annals of Physics) 的編輯迫不及待發(fā)表這篇文章，也說(shuō)明這篇論文有所原創(chuàng)。但是，之前不斷的拒稿，弄得我有點(diǎn)氣急敗壞，不但把論文的第一句改寫成一句廣告：“本文旨在解決統(tǒng)計(jì)物理中一個(gè)歷史很長(zhǎng)的基礎(chǔ)性難題”，還迫不及待地到2020-2021年中國(guó)物理學(xué)秋季年會(huì)上做了一場(chǎng)報(bào)告進(jìn)行推廣。論文接收后，居然還有點(diǎn)不敢相信；校完樣，還特地買了一盆花來(lái)紀(jì)念下自己；文章發(fā)表了，還要寫這篇文章來(lái)吹噓一下。

學(xué)術(shù)乃天下公器，誠(chéng)摯期待各位專家、老師、朋友批評(píng)指正。

參考文獻(xiàn)

[1] Q. H. Liu，Asynchronous finite differences in most probable distribution with finite numbers of particles，Annals of Physics，Vol. 441, June 2022, 168884. https://doi.org/10.1016/j.aop.2022.168884

特 別 提 示

1. 進(jìn)入『返樸』微信公眾號(hào)底部菜單“精品專欄“，可查閱不同主題系列科普文章。

2. 『返樸』提供按月檢索文章功能。關(guān)注公眾號(hào)，回復(fù)四位數(shù)組成的年份+月份，如“1903”，可獲取2019年3月的文章索引，以此類推。

版權(quán)說(shuō)明：歡迎個(gè)人轉(zhuǎn)發(fā)，任何形式的媒體或機(jī)構(gòu)未經(jīng)授權(quán)，不得轉(zhuǎn)載和摘編。轉(zhuǎn)載授權(quán)請(qǐng)?jiān)凇阜禈恪刮⑿殴娞?hào)內(nèi)聯(lián)系后臺(tái)。