①物理意義與直觀解釋 (我寫的教學說明)
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.112v 為(每種極化視為常數的)聲速;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.112N 為試樣中的原胞數;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
③推導所使用的假設 (書中)
- Debye 近似把聲速對每一種偏振都取為常數,就像古典彈性連續體那樣(p.112)
- 色散關係寫成 ω = vK,v 為該常數聲速(p.112)
- 對邊長 L 的立方體內 N³ 個原胞施加週期性邊界條件,於是 K 空間中每體積 (2π/L)³ 對應一個允許的 K,即每單位 K 空間體積有 V/8π³ 個允許值(p.111)
- 若試樣有 N 個原胞,聲學聲子模式總數即為 N;由此定出截止頻率 ωD³ = 6π²v³N/V 與截止波向量 KD = ωD/v = (6π²N/V)^{1/3},且不允許 K > KD 的模式(p.112)
- 為求簡潔,假設聲子速度與偏振無關,因而直接乘上因子 3(p.112)
- 熱能由 Planck 分布佔據數積分而得:U = ∫dω D(ω)⟨n(ω)⟩ℏω(p.112)
- Debye 溫度定義為 θ = (ℏv/kB)(6π²N/V)^{1/3}(p.112)
④適用條件與限制 (書中)
- 低溫下 T³ 近似「相當好」的條件是:只有長波長聲學模式被熱激發,而這些正是可以用宏觀彈性常數當成彈性連續體處理的模式(p.114)
- 此近似會失效的短波長模式,其能量太高,低溫下不會被顯著佔據(p.114)
- 實際晶體要看到相當純的 T³ 行為,溫度必須很低——書中說可能需要低於 θ/50(p.114)
- Fig. 7 標明 T³ 律的區域在 0.1θ 以下(p.113)
- Debye 的態密度並非實際晶體的態密度:Fig. 14 對比 Debye 固體與真實晶體,真實晶體的譜在小 ω 時同樣起於 ω²,但在奇異點處出現不連續,即 Van Hove 奇異點(p.119)
- 高溫端趨近 3NkB(p.113),但「高溫熱容為常數」屬於簡諧理論的推論,書中明說真實晶體並不精確滿足,偏差來自非簡諧項(p.119、p.120)
- 金屬在遠低於 θ 與 TF 時,熱容須寫成電子與聲子兩項之和 C = γT + AT³,Debye 的 T³ 只是其中的 AT³ 項(p.145)
⑤高溫與低溫行為 (書中)
高溫 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 章節、小節與頁碼
- Ch5 PHONON HEAT CAPACITY p.107
- Ch5 Density of States in Three Dimensions p.111
- Ch5 Debye Model for Density of States p.112
- Ch5 Debye T3 Law p.114
- Ch6 Experimental Heat Capacity of Metals p.145
書中逐字
"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。