Kittel·固態物理

熱容專題 / 主題 5

05

Debye 模型與 Debye 溫度

溫區定位:過渡區 T ~ θ · 熱容爬升最陡的一段

①物理意義與直觀解釋 (我寫的教學說明)

Debye(1912)修正了 Einstein 最致命的問題:他保留了「長波模式頻率趨近零」這件事。做法是把晶體當成古典彈性連續體,色散取 ω = vK,聲速 v 為常數。

這個近似在長波(低頻)時非常好——波長遠大於晶格常數時,晶體看起來確實就是連續介質,這正是 Ch3 彈性常數的適用範圍。它在短波時很差(真實色散在 Brillouin 區邊界會彎平),但低溫下短波模式本來就沒被激發,所以不影響低溫結果。

模式總數必須守恆:一個含 N 個原胞的晶體只有 3N 個模式,不能無限多。Debye 用一個截斷頻率 ω_D(對應截斷波向量 K_D)把模式數卡在 3N。這個截斷換算成溫度就是 Debye 溫度 θ = ℏω_D/k_B——它是整個理論唯一的材料參數,而且有明確的物理意義:

θ 是「最高頻的晶格振動模式,其能量量子相當於多高的溫度」。所以 T ≫ θ 表示所有模式都被充分激發(古典區,回到 Dulong–Petit);T ≪ θ 表示只有低頻的一小部分被激發(量子區,T³ 律)。θ 因此是判斷「這個溫度算高還是低」的唯一尺規——同一個 300 K,對 θ = 92 K 的固態氬是高溫,對 θ ≈ 645 K 的矽卻是低溫。這解釋了為什麼矽在室溫的比熱還沒到古典極限。

②重要公式與符號定義 (書中)

Debye 模型的態密度

ω = vK ⟹ D(ω) = Vω² / (2π²v³)

(21), (22)  Ch5 p.112

v 為(每種極化視為常數的)聲速;V 為試樣體積;由一般式 D(ω) = dN/dω = (VK²/2π²)(dK/dω)(式 (20))代入 ω = vK 得到

Debye 截斷頻率與截斷波向量

ω_D³ = 6π²v³N/V ; K_D = ω_D/v = (6π²N/V)^{1/3}

(23), (24)  Ch5 p.112

N 為試樣中的原胞數;Debye 模型不允許 K > K_D 的模態,K ≤ K_D 的模態數恰好用盡單原子晶格的自由度

Debye 溫度

θ = (ℏv/k_B) (6π²N/V)^{1/3}

(28)  Ch5 p.112

θ 為 Debye 溫度,由 ω_D 定義;並定義 x ≡ ℏω/k_BT、x_D ≡ ℏω_D/k_BT = θ/T

Debye 模型熱容(完整式)

C_V = 9Nk_B (T/θ)³ ∫₀^{x_D} dx · x⁴ eˣ / (eˣ − 1)²

(30)  Ch5 p.113

由 U = 9Nk_BT(T/θ)³∫₀^{x_D}dx x³/(eˣ−1)(p.113 明言的式 (26) 中間表達式)對溫度微分;x_D = θ/T;T ≫ θ 時趨近古典值 3Nk_B

③推導所使用的假設 (書中)

④適用條件與限制 (書中)

⑤高溫與低溫行為 (書中)

高溫 T ≫ θ

T ≫ θ 時 Debye 熱容趨近古典值 3NkB(p.113)。Fig. 7 給出 T/θ 大時的漸近值為 24.943 J mol⁻¹ deg⁻¹(p.113)。矽與鍺的實測值(Fig. 8)在室溫附近即接近此平台,並在低溫下降(p.113)。

低溫 T ≪ θ

極低溫時把 (29) 式的積分上限取為無窮,得 U ≅ 3π⁴NkBT⁴/5θ³,於是 T ≪ θ 時 CV ≅ (12π⁴/5)NkB(T/θ)³ ≅ 234 NkB(T/θ)³,即 Debye T³ 近似(p.114)。書中另給一個簡單論證:只有 ℏω < kBT 的模式會被明顯激發,其在 K 空間所佔分率約為 (KT/KD)³ = (T/θ)³,故約有 3N(T/θ)³ 個模式各帶能量 kBT,得 U ≈ 3NkBT(T/θ)³ 與 CV ≈ 12NkB(T/θ)³(p.114、Fig. 10 說明 p.115)。固態氬的低溫實驗對 T³ 作圖與 Debye T³ 律極佳吻合,θ = 92.0 K(Fig. 9,p.115)。

⑥對應的 Kittel 章節、小節與頁碼

書中逐字

"In the Debye approximation the velocity of sound is taken as constant for each polarization type, as it would be for a classical elastic continuum." (p.112)

"At sufficiently low temperature the T3 approximation is quite good; that is, when only long wavelength acoustic modes are thermally excited." (p.114)

"It may be necessary to be below T = θ/50 to get reasonably pure T3 behavior." (p.114)

【延伸知識】書中沒有處理的部分

【範圍已校正】書中有的:熱導率 K 與熱阻率(Ch5「Thermal Resistivity of Phonon Gas」p.123)、Umklapp 過程(p.125)、缺陷與晶界的尺寸效應與 Casimir 邊界散射(p.126–127)、體積比熱(p.122 式 42 明寫 C is the heat capacity per unit volume,Table 2 單位為 J cm⁻³ K⁻¹)、以及一維通道的彈道熱導 G_th(Ch18 p.561–562、564)。全書零次出現的(已逐頁搜尋確認):Kapitza/thermal boundary resistance、thermal time constant、thermal impedance。因此真正「書中沒有」的是元件層級的暫態熱分析那一整套——熱阻抗 Zth(t)/Zth(jω)、熱時間常數 τ = Rth·Cth、Foster/Cauer 熱 RC 網路、structure function,以及界面熱阻 TBR。