第二章 熱力學第二定律 2.4 具有大量粒子的系統

第二章 熱力學第二定律 2.4 具有大量粒子的系統

來自專欄 Schroeder:從熱力學到統計力學4 人贊了文章

2.4 具有大量粒子的系統

By @Hey'u & @Charge & @胡大師

在上一節中,我們看到,對於兩個存在相互作用的愛因斯坦固體(每個都有一百個左右的諧振子)組成的系統,某些宏觀狀態比其他的更可能出現。然而,仍舊有大約20%的宏觀狀態可能性都不低。接下來我們將看看當系統更大時——每個固體包含10^{20}個甚至更多的諧振子——會發生什麼。在本節結束時,我們的目標是知道,在所有的宏觀狀態中,只有很小一部分是有明顯可能性的。換句話說,重數的函數變得非常尖銳(見圖 2-6)。然而,要分析這樣大的系統,我們必須首先了解非常大的數字的數學規律。

圖 2-6 兩個相互作用的愛因斯坦固體的典型的重數圖,左圖是包含幾百個諧振子和能量單位的系統,右圖則是幾千個。隨著系統尺寸的增加,峰值相對於整個橫軸變得非常窄。對於N~q~10^20的情況,峰值太尖銳而無法繪製。

非常大的數字

在統計力學中,存在三種常見的數字:小數字、大數字和非常大的數字。

小數字(Small number)指的是像6、23、42這種數字。你早就知道怎麼處理這些數字。

大數字(Large number)指的是比小數字大得多的數,通常由小數字的乘方組成。統計力學中最重要的大數字就是10^{23}數量級的阿伏伽德羅常數。對於大數字,最重要的性質是把它和一個小數字相加所得的結果其實就是它自己。舉個例子:

{10}^{32}+23={10}^{23} \(2.12)

(唯一的例外是你最後又把這個大數字減掉了:{10}^{32}+42-{10}^{23}=42。)

非常大的數字(Very large number)是比大數字還大的多的數字,往往由大數字的乘方組成。比如{10}^{{10}^{23}}。這種數字有非常漂亮的性質:你可以把它和一個大數字乘起來而不改變它的值。例如:

{10}^{{10}^{23}}	imes{10}^{23}={10}^{left({10}^{23}+23
ight)}={10}^{{10}^{23}} \(2.13)

這個性質你可能一上來難以接受,但是它對於非常大的數字十分有用。(當然,如果你最後又把這個非常大的數字除掉了,你就不能這樣寫了。)

對於處理非常大的數字,我們有一個常用的手段:取對數。這個方法把一個非常大的數字變成了一個一般的大數字,這樣我們就可以更好地對它進行運算,只要最後再乘方回去就可以了。在這一節我們就會使用這個手段。

斯特林公式

我們的重數(微觀狀態數)的公式總是和組合數有關,而組合數是由階乘組成的。為了在具有大量粒子的系統中應用這些公式,我們需要一些特別的方法來估計大數字的階乘。這個方法就是斯特林公式(Stirlings Approximation)

N!approx N^Ne^{-N}sqrt{2pi N} \(2.14)

Ngg1時,這個公式十分準確。為什麼呢?我們下面就來解釋一下。

N!是從1NN個因子的乘積,如果我們把每一個因子都替換成N,這個數就是N!approx N^N——一個十分粗略的近似。這個近似比N!實際的值大得多;事實上,平均下來每個因子都大了e倍,也就是說

N!approxleft(frac{N}{e}
ight)^N=N^Ne^{-N} \(2.15)

這還比N!相差了一個比較大的因子:sqrt{2pi N}。但若N是一個大數字,那麼N!就是一個非常大的數字,所以這個因子就可以被省略了。

如果我們只關心N!的對數,上式通常就足夠準確了:

ln{N!}approx Nln{N}-N \(2.16)

表格 2-3是在一些不太大的數上應用斯特林公式的結果。我們可以看到,即使N不太大,斯特林公式仍舊很有用——公式(2.14)N!approx N^Ne^{-N}sqrt{2pi N}甚至在N=10的時候就十分準確了;公式(2.16)ln{N!}approx Nln{N}-N則在N=100的時候就很準確了。

斯特林公式的推導詳見附錄B.3。(我們將在新學期到來之前更新)

表格 2-3 斯特林公式與真實值的對比(在N=1, 10, 100時)。

大型愛因斯坦固體的重數

了解了斯特林公式以後,我們現在就可以來估計具有大量諧振子和能量單位的愛因斯坦固體了。我們現在只考慮qgg N(這其實是「高溫」下的極限)——也即能量單位數遠大於諧振子數——的情況,而不考慮一般的情況。

我們先來看原先推導出的準確公式:

Omegaleft(N,q
ight)=inom{q+N-1}{q}=frac{left(q+N-1
ight)!}{q!(N-1)!}approxfrac{left(q+N
ight)!}{q!N!} \(2.17)

我們實際上乘了一個frac{q+N}{N}=frac{q}{N}+1來得到最右邊的值,Omegaleft(N,q
ight)實際上是一個非常大的數字,而frac{q}{N}最多只是一個大數字,因此這個近似其實較為準確。接下來,我們對這個公式取對數並應用對數的斯特林公式:

egin{align} ln{Omega}&=ln{left[frac{left(q+N
ight)!}{q!N!}
ight]}\ &=ln{left(q+N
ight)!}-ln{q!}-ln{N!}\ &approxleft(q+N
ight)ln{left(q+N
ight)}-left(q+N
ight)-qln{q}+q-Nln{N}+N\ &=left(q+N
ight)ln{left(q+N
ight)}-qln{q}-Nln{N} end{align} \(2.18)

目前為止我們還沒有利用qgg N,和問題2.13一樣:

ln{left(q+N
ight)}=ln{left[qleft(1+frac{N}{q}
ight)
ight]}=ln{q}+ln{left(1+frac{N}{q}
ight)}approxln{q}+frac{N}{q} \(2.19)

上式中的最後一步,我們利用了對數在x=0處的泰勒展開ln{left(1+x
ight)}approx x。(我們可以很輕鬆地證明截斷的余項在left|x
ight|<1時收斂並且是limlimits_{x	o 0}時的更高階的無窮小。)將這個結果代入ln{Omega},我們得到:

egin{align} ln{Omega} &approxleft(q+N
ight)left(ln{q}+frac{N}{q}
ight)-qln{q}-Nln{N}\ &=N+Nln{q}+frac{N^2}{q}-Nln{N}\ &=Nln{frac{q}{N}}+N+frac{N^2}{q} end{align} \(2.20)

qgg N時,最後一項其實可以忽略。所以:

Omegaleft(N,q
ight)approx e^{Nln{left({q}/{N}
ight)}}e^N=left(frac{eq}{N}
ight)^N \(2.21)

這個公式簡潔而又美觀,但是有點奇特:它的指數是一個大數字,因此Ω是一個非常大的數字;同時,即使我們略微地增加N或q,Ω都會因為這個巨大的指數N而暴漲。

重數函數的銳度

最後,我們終於回到了這一節一開始提出的問題:對於兩個有相互作用的大型愛因斯坦固體,它們的重數函數的峰到底有多尖銳?

簡單起見,我們假設每個固體都有N個諧振子。我們將總的能量單位數記為q(而非q_{mathrm{total}})。我們同樣假設qN大得多。那麼對於任意給定的宏觀狀態,這個系統的重數(微觀狀態數)是:

Omega=left(frac{eq_A}{N}
ight)^Nleft(frac{eq_B}{N}
ight)^N=left(frac{e}{N}
ight)^Nleft(q_Aq_B
ight)^N \(2.22)

其中q_Aq_B是兩個固體的能量單位數(q_A+q_B=q)。

如果我們把上式作為一個q_A的函數作圖,我們會發現它在q_A=frac{q}{2}的地方有一個非常尖銳的峰,這個峰的高度是一個非常大的數:

Omega_{mathrm{max}}=left(frac{e}{N}
ight)^Nleft(frac{q}{2}
ight)^{2N} \(2.23)

我們還對這個峰附近的形狀感興趣,所以我們令

q_A=frac{q}{2}+x,quad q_B=frac{q}{2}-x \(2.24)

這裡,x可以是任意的比q小的多的數(但是可能還是挺大的)。將它代入公式(2.22)Omega=left(frac{e}{N}
ight)^Nleft(q_Aq_B
ight)^N

Omega=left(frac{e}{N}
ight)^Nleft[left(frac{q}{2}
ight)^2-x^2
ight]^N \(2.25)

為了簡化這個式子,我們還是對它取對數:

egin{align} ln{left[left(frac{q}{2}
ight)^2-x^2
ight]^N} &=Nln{left[left(frac{q}{2}
ight)^2-x^2
ight]}\ &=Nln{left{left(frac{q}{2}
ight)^2left[1-left(frac{2x}{q}
ight)^2
ight]
ight}}\ &approx Nleft[ln{left(frac{q}{2}
ight)^2}-left(frac{2x}{q}
ight)^2
ight] end{align} \(2.26)

Omega=left(frac{e}{N}
ight)^Nexp{left[Nln{left(frac{q}{2}
ight)^2}
ight]}exp{left[-Nleft(frac{2x}{q}
ight)^2
ight]}=Omega_{mathrm{max}}cdot e^{-Nleft({2x}/{q}
ight)^2} \(2.27)

這樣形式的函數被稱為高斯函數:它在x=0處有一個峰值並且在兩邊都快速衰減,就如圖 2-7所示。重數掉到這個峰值高度的1/e時,

Nleft(frac{2x}{q}
ight)^2=1Leftrightarrow x=frac{q}{sqrt N} \(2.28)

這個值其實是一個大數字,但是如果N={10}^{20}的話,它不過是整個橫軸的百億分之一罷了。

圖 2-7 兩個有相互作用的大型愛因斯坦固體組成的系統的重數圖。愛因斯坦固體的每個數諧振子具有許多份能量(高溫極限)。橫軸只畫了一小部分。

在這個圖中的比例下,峰的寬度大約一厘米,那麼整個橫軸的長度就有10^{10}厘米,也就是十萬千米——可以繞地球兩圈還多。在這個圖中最邊緣的地方,即使x只比{q}/{sqrt N}大十倍,這個x對應的重數就比最大值小了e^{100}approx{10}^{44}倍。

這個結果就告訴了我們,當兩個大型愛因斯坦固體互相到達了熱平衡時,任何偏離最可能的宏觀狀態的隨機擾動都幾乎無法觀測到。如果我們真的想要觀測到這些擾動,我們必須將能量測量到十個有效數字的準確度。因為所有的微觀狀態都是等可能的,一旦系統達到熱平衡,我們不妨假設它就處在最可能的宏觀狀態。當一個系統的大小趨於無窮的時候,任何偏離最可能的宏觀態的擾動就永遠不會出現了,這被稱作熱力學極限

總結

=========================================

  • 本節從非常大的數字的數學規律講起,說明了斯特林公式(N的階乘的近似公式)的來源。
  • 利用斯特林公式,得出了兩個存在相互作用的具有大量諧振子和更大量的能量單位數的愛因斯坦固體的重數(微觀狀態數)。
  • 接著算出了這個重數函數的峰有多麼銳利:重數掉到它的峰值高度的1/e時,峰的寬度也只有{q}/{sqrt N},是如果N的大小在阿伏伽德羅常數的量級上,寬度比上整個橫軸q也非常小。
  • 證明了在所有的宏觀狀態中,只有很小一部分是有明顯可能性的
  • 熱力學極限:當一個系統的大小趨於無窮的時候,任何偏離最可能的宏觀態的擾動就永遠不會出現了。

--------------------------------------------------------------------

Content Created: 2018年7月10日

Last updated:2018年7月10日

歡迎各位批評指正,歡迎留言或者私信。


推薦閱讀:

這位英國lady, 放著好好的攝影師不做, 偏要漂洋過海撿垃圾
論一塊好泥是怎樣煉成的
蒼蠅眼中的慢動作世界
籌備中的朱?國際論壇
假如地球人一起跳會引發地震嗎?

TAG:物理學 | 熱力學 | 自然科學 |