德拜模型

在熱力學和固體物理學中，德拜模型（英語：Debye model）是由彼得·德拜在1912年提出的方法，用於估算聲子對固體的比熱（熱容）的貢獻。德拜模型把原子晶格的振動（熱）當作盒中的聲子處理，而與此不同的愛因斯坦模型則將固體作為許多單獨的、不相互作用的量子諧振子處理。德拜模型正確地預言了低溫時固體的熱容：與$T^3$成正比（德拜$T^3$定律）。與愛因斯坦模型一樣，它在高溫時也與杜隆-珀蒂定律相符合。但由於模型的簡單假設，它在中間的溫度不太準確。

## 推導

德拜模型是普朗克黑體輻射定律的固體版本，後者把電磁輻射視為盒中的光子氣體，而德拜模型把原子的振動視為盒中的聲子（盒子就是這塊固體）。它們大部分的計算步驟都是相同的，因為聲子氣和光子氣都是有線性色散關係的玻色氣體。

考慮一個邊長為$L$的立方體。從盒中粒子一文可知，盒中的（沿某一邊長方向傳播的）聲波由於共振效應，其波長只能為：

$$\lambda_n = {2L\over n}\,,$$

其中$n$是整數。一個聲子的能量是：

$$E_n\ =h\nu_n\,,$$

其中$h$是普朗克常數，$\nu_{n}$是聲子的頻率。近似認為頻率與波長成反比，則有：

$$E_n=h\nu_n={hc_s\over\lambda_n}={hc_sn\over 2L}\,,$$

其中$c_s$是固體中的聲速。推廣到三維的情況（沿任意方向傳播的聲子）：

$$E_n^2={p_n}^2{c_s}^2=\left({hc_s\over2L}\right)^2\left(n_x^2+n_y^2+n_z^2\right)\,.$$

其中$p_n$為聲子的動量。頻率與波長成反比的近似（即聲速恆定）對於低能量聲子是準確的，但對於高能量聲子則不然（參見聲子）。這是德拜模型的局限之一，導致其對中間溫度的預言不正確，而對低溫和高溫都很精確。

下面來計算盒中的總能量：

$$U = \sum_n E_n\,\bar{N}(E_n)\,,$$

其中$\bar{N}(E_n)$是盒中能量為$E_n$的聲子的數目。換句話說，總能量等於某能量的聲子數目乘以其能量再對各能量求和。在三維空間中，我們有：

$$U = \sum_{n_x}\sum_{n_y}\sum_{n_z}E_n\,\bar{N}(E_n)\,.$$

這是德拜模型和普朗克黑體輻射定律的不同之處。與盒中的電磁輻射不一樣，由於聲子不能有無限大的頻率，所以盒中聲子只有有限個能量狀態。它的頻率由它的傳播介質——固體的原子晶格所約束。考慮以下的橫向聲子的插圖：

可以合理假設聲子的最小波長是原子間距的兩倍，如最下一圖所示。固體中有$N$個原子，我們的固體是正方體，因此每一條邊有$\sqrt[3]{N}$個原子，原子間距為$L/\sqrt[3]{N}$，從而最小波長為：

$$\lambda_{\rm min} = {2L \over \sqrt[3]{N}}\,,$$

最大的模數$n$（對於光子是無限大）為：

$$n_{\rm max} = \sqrt[3]{N}\,.$$

這設置了三重能量求和的上限：

$$U = \sum_{n_x}^{\sqrt[3]{N}}\sum_{n_y}^{\sqrt[3]{N}}\sum_{n_z}^{\sqrt[3]{N}}E_n\,\bar{N}(E_n)\,.$$

對於緩慢變化的、表現良好的函數，求和可以用積分來代替（又稱為托馬斯-費米近似）：

$$U \approx\int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}} E(n)\,\bar{N}\left(E(n)\right)\,dn_x\, dn_y\, dn_z\,.$$

至此，我們還沒有提到$\bar{N}(E)$：能量為$E\,$的聲子的數目。聲子服從玻色-愛因斯坦統計。它們的分布由著名的玻色-愛因斯坦公式給出：

$$\langle N\rangle_{BE} = {1\over e^{E/kT}-1}\,.$$

由於一個聲子有三個可能的偏振態（一個縱向、兩個橫向，大致攜帶相同的能量），需把以上公式乘以3：

$$\bar{N}(E) = {3\over e^{E/kT}-1}\,.$$

實際上我們可以定義等效聲速$c_s:=c_{{\rm eff}}$，也就是說，德拜溫度$T_d$（見下文）與$c_{{\rm eff}}$成正比。更加精確地講，$T_D^{-3}\propto c_{{\rm eff}}^{-3}:=(1/3)c_{{\rm long}}^{-3}+(2/3)c_{{\rm trans}}^{-3}$，其中我們區分了縱向和橫向的聲波速度（貢獻分別為1/3和2/3）。德拜溫度或有效聲速是晶體的硬度的一種衡量。

把此式代入能量積分，得：

$$U = \int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}} E(n)\,{3\over e^{E(n)/kT}-1}\,dn_x\, dn_y\, dn_z\,.$$

這比光子的情況容易得多，因為（在半經典式的處理中）光子的頻率（也即這個三重積分）是無上界的。但如上面圖所示，這一點對於聲子不成立。為了近似計算這個三重積分，德拜使用了球坐標系：

$$\ (n_x,n_y,n_z)=(n\cos \theta \cos \phi,n\cos \theta \sin \phi,n\sin \theta )$$

並用八分之一球（球在一個卦限內的部分）來近似代替立方體：

$$U \approx\int_0^{\pi/2}\int_0^{\pi/2}\int_0^R E(n)\,{3\over e^{E(n)/kT}-1}n^2 \sin\theta\, dn\, d\theta\, d\phi\,,$$

其中$R$是球的半徑，通過保持立方體和八分之一球中的粒子數目相同來得出。立方體的體積是$N$：

$$N = {1\over8}{4\over3}\pi R^3\,,$$

因此：

$$R = \sqrt[3]{6N\over\pi}\,.$$

用球的積分來代替立方體積分這一近似是德拜模型不準確性的另一來源。

先積分掉角向部分：

$$U = {3\pi\over2}\int_0^R \,{hc_sn\over 2L}{n^2\over e^{hc_sn/2LkT}-1} \,dn$$

再利用變量代換$x = {hc_sn\over 2LkT}$：

$$U = {3\pi\over2} kT \left({2LkT\over hc_s}\right)^3\int_0^{hc_sR/2LkT} {x^3\over e^x-1}\, dx$$

為簡化表達式，可以定義德拜溫度$T_D$——它的量綱與溫度相同，因物質而異：

$$T_D\ \stackrel{\mathrm{def}}{=}\ {hc_sR\over2Lk} = {hc_s\over2Lk}\sqrt[3]{6N\over\pi} = {hc_s\over2k}\sqrt[3]{{6\over\pi}{N\over V}}$$

於是比內能為：

$$\frac{U}{Nk} = 9T \left({T\over T_D}\right)^3\int_0^{T_D/T} {x^3\over e^x-1}\, dx = 3T D_3 \left({T_D\over T}\right)\,,$$

其中$D_3(x)$是（第三）德拜函數。

對$T$微分，我們便得到無量綱熱容：

$$\frac{C_V}{Nk} = 9 \left({T\over T_D}\right)^3\int_0^{T_D/T} {x^4 e^x\over\left(e^x-1\right)^2}\, dx\,.$$

這個公式給出了任何溫度下德拜模型的結果。下面更簡單的公式給出了低溫和高溫極限下的漸近表現。前面已經提到，這兩個漸進表現是精確的，而中間溫度的表現則不然。其根本原因在於，低溫和高溫近似下，德拜模型都給出了其一，低頻下精確的色散關係$E(\nu )$，其二，精確的狀態密度（單位頻率間隔內的聲子數目）$(\int g(\nu ) \, {\rm d\nu}\equiv 3N)$。

## 德拜的推導

實際上，德拜用不同且更簡單的方法推出了這個公式。利用連續介質的固體力學，他發現頻率小於某個特定值的振動狀態的數目趨近於：

$$n \sim {1 \over 3} \nu^3 V F\,,$$

其中$V$是體積，$F$是一個從彈性係數和密度計算出的因子。把這與溫度T的量子諧振子的期望能量（已經由愛因斯坦在他的模型中使用）結合，便給出能量：

$$U = \int_0^\infty \,{h\nu^3 V F\over e^{h\nu/kT}-1}\, d\nu\,,$$

如果振動頻率趨於無窮大。這個形式給出了$T^3$的表現，它在低溫時是正確的。但德拜意識到N個原子不可能有超過$3N$個振動狀態。他假設在原子固體中，振動狀態的頻譜將繼續遵循以上的規則，直到一個最大頻率$\nu_m$使得總的狀態數目為$3N$為止：

$$3N = {1 \over 3} \nu_m^3 V F \,.$$

德拜知道這個假設不是真正正確的（較高的頻率比假設要更加密集），但它保證了高溫時的正確表現（杜隆-珀蒂定律）。於是，能量由以下給出：

$$U = \int_0^{\nu_m} \,{h\nu^3 V F\over e^{h\nu/kT}-1}\, d\nu\,,$$

$$= V F kT (kT/h)^3 \int_0^{T_D/T} \,{x^3 \over e^x-1}\, dx\,,$$

其中$T_D$是$h\nu_m/k$。

$$= 9 N k T (T/T_D)^3 \int_0^{T_D/T} \,{x^3 \over e^x-1}\, dx\,,$$

$$= 3 N k T D_3(T_D/T)\,,$$

其中$D_3$是一個函數，後來命名為三階德拜函數。

## 低溫極限

德拜固體的溫度稱為低的，如果$T \ll T_D$，在這個情況下：

$$\frac{C_V}{Nk} \sim 9 \left({T\over T_D}\right)^3\int_0^{\infty} {x^4 e^x\over \left(e^x-1\right)^2}\, dx$$

這個定積分可以精確計算：

$$\frac{C_V}{Nk} \sim {12\pi^4\over5} \left({T\over T_D}\right)^3$$

在低溫極限中，德拜模型的局限不適用，它給出了（聲子）熱容、溫度、彈性係數，以及每個原子的體積（後面的數量是包含在德拜溫度之中的）之間的正確關係。

## 高溫極限

德拜固體的溫度稱為高的，如果$T \gg T_D$。$e^x - 1\approx x$如果$|x|<<1$，得出：

$$\frac{C_V}{Nk} \sim 9 \left({T\over T_D}\right)^3\int_0^{T_D/T} {x^4 \over x^2}\, dx$$

$$\frac{C_V}{Nk} \sim 3\,.$$

這就是杜隆-珀蒂定律，它是相當準確的，雖然它沒有考慮非諧性，這造成了熱容進一步上升。如果固體是導體或半導體，那麼它的總熱容還可能含有電子的不可忽略的貢獻。

## 外部連結

- 用低溫恆溫器來決定石英的比熱、熱導率和電導率