02聲子熱容：三維態密度與 Debye 模型

這一段在解決什麼問題我的判讀

第 5 章從 p.107 開場就把整條路線圖畫好了。式(1) 說晶格的總能量是「把所有聲子模態的能量加起來」：U_lat = Σ_K Σ_p ⟨n_K,p⟩ℏω_K,p。式(7) 的 Planck 分布給了每個模態在溫度 T 下的平均佔據數 ⟨n⟩ = 1/(e^{ℏω/τ} − 1)。到這裡，每一個模態該貢獻多少能量，已經完全解決了。

剩下的只有一個問題，而且是純粹的計數問題：到底有多少個模態，它們的頻率怎麼分布？

Kittel 在 p.108 把這件事講得很白：「The central problem is to find D(ω), the number of modes per unit frequency range. This function is called the density of modes or, more often, density of states.」把式(1) 的離散求和換成積分之後，式(9) 變成 U = Σ_p ∫dω D_p(ω) ℏω/(e^{ℏω/τ} − 1)，式(10) 給出 C_lat = k_B Σ_p ∫dω D_p(ω) x² e^x/(e^x − 1)²。這兩條式子裡唯一未知的就是 D(ω)。

p.109-110 已經在一維線鏈上把 D₁(ω) 算完了（式(15)：D₁(ω)dω = (L/π)(dω)/(dω/dK)）。但一維只是暖身——真實晶體是三維的，而且三維的 D(ω) 有一個一維完全看不出來的結構：D(ω) ∝ ω²。這個 ω² 就是低溫 T³ 定律的全部來源。

所以這一段（p.111 下半到 p.113）要做兩件事：

第一件（Density of States in Three Dimensions, p.111-112）：把「一個立方體裡有多少個允許的 K」這個純幾何問題算乾淨，得到式(20) D(ω) = (VK²/2π²)(dK/dω)。這一步完全不需要任何物理假設——它只用到週期性邊界條件，對任何色散關係都成立。

第二件（Debye Model for Density of States, p.112-113）：引入 Debye 近似 ω = vK（式(21)），把 D(ω) 具體化成 Vω²/2π²v³（式(22)），再用「N 個原胞只能有 N 個聲學模態」把頻率截在 ω_D（式(23)），最後一路積到熱容式(30)，並在高溫回到古典的 3Nk_B。

這兩件事之間有一個常被忽略的分界：第一件是恆等式，第二件是近似。讀這一段時如果把兩者混在一起，後面遇到「為什麼 T³ 定律在低溫是精確的、但 Debye 模型在中溫是錯的」就會答不出來。

還有一件事值得先講：這一段裡冒出來的常數 8π³ 與 2π²，是整章最常被抄錯的兩個數字。它們都不是查表來的，而是從 (2π/L)³ 與 (4π/3)/(8π³) 這兩個幾何量硬算出來的。下面會把每一步都攤開，讀完之後你應該永遠不需要再背它們。

原文的起手式：從一維的奇異點跨到三維

p.111 的上半頁還在收尾一維的討論。Kittel 寫下一句在後面（p.118-119 與習題 1）會回來咬人的話：

「We can obtain the group velocity dω/dK from the dispersion relation ω versus K. There is a singularity in D₁(ω) whenever the dispersion relation ω(K) is horizontal; that is, whenever the group velocity is zero.」

這句話的分量要看懂。一維的式(15) 是 D₁(ω) = (L/π)/(dω/dK)——群速度在分母。凡是色散曲線出現水平切線的地方（單原子鏈的區邊界 K = π/a、光學支的 K = 0），群速度歸零，D₁(ω) 就發散。這是 van Hove 奇異點的一維版本。

Kittel 並沒有把這個觀察擱著——他在同一章把它做完了。p.118 他把式(19) 的球形計數推廣成任意等頻率面的一般式(37)：「We divide both sides by dω and write V = L³ for the volume of the crystal: the result for the density of states is (37) D(ω) = V/(2π)³ ∫ dS_ω/v_g.」（同頁已定義 v_g = |∇_Kω|。）p.119 接著替這些點正式命名：「There is a special interest in the contribution to D(ω) from points at which the group velocity is zero. Such critical points produce singularities (known as Van Hove singularities) in the distribution function (Fig. 14).」而 Fig. 14 就是把兩種態密度並排畫給你看，圖說逐字為：「Density of states as a function of frequency for (a) the Debye solid and (b) an actual crystal structure. The spectrum for the crystal starts as ω² for small ω, but discontinuities develop at singular points.」——Debye 的光滑曲線與真實晶體的不連續，是書上明講、明畫的對比。

p.128 的習題 1（「Singularity in density of states」）則是要你把 Fig. 14 的定性圖算成定量式子：一維單原子鏈的 D(ω) = (2N/π)·1/(ω_m² − ω²)^{1/2}，以及三維光學支 ω(K) = ω₀ − AK² 附近的 D(ω) = (L/2π)³(2π/A^{3/2})(ω₀ − ω)^{1/2}。

接著 Kittel 就直接開新小節 Density of States in Three Dimensions，一句廢話都沒有。

這個奇異點預告了 Debye 近似要付的代價我的判讀

（本節是本講解把 p.111 的觀察與 p.119 的 Fig. 14 接起來的判讀框架，Kittel 在 p.111-113 的行文裡並沒有這樣整理。）

為什麼開場要提奇異點？因為它預告了接下來整段最大的代價：Debye 模型把色散關係硬寫成 ω = vK，也就是群速度處處等於 v、永遠不為零。換句話說，Debye 模型在構造上就消滅了所有 van Hove 奇異點。真實晶體的 D(ω) 在區邊界附近有尖銳的峰與扭折，Debye 的 D(ω) 卻是一條光滑的拋物線 Vω²/2π²v³，然後在 ω_D 突然垂直落到零。

這不是臆測：p.119 的 Fig. 14 就是把 (a) Debye 固體與 (b) 真實晶體的 D(ω) 並排，圖說明寫真實晶體同樣以 ω² 起步、但在臨界點處發展出不連續。讀 p.111-113 時把這張圖記在腦子裡，後面每次看到「Debye 的拋物線」就會自動想到它抹掉了什麼；做完 p.128 的習題 1，才會真正明白那條拋物線有多粗糙。

為什麼上界恰好是 Nπ/L：截止不是人為塞進去的我的判讀

這一點在課堂上很少被講清楚，但它決定了你能不能理解後面式(23) 的 ω_D。

式(17) 的允許值從 0 開始、以 2π/L 為間距往兩邊長，每個方向共 N 個（含 K = 0 的均勻模態）。三維合起來 N³ 個。而樣品裡有 N³ 個原胞。所以：

允許的 K 點數 = 原胞數。

這不是巧合，而是離散晶格的一條硬定理。它的物理內容是：一個由 N³ 個原胞組成的晶格，每個原胞的每個原子有 3 個自由度；若每個原胞只有 1 個原子，總自由度是 3N³，正好分配給 N³ 個 K 點 × 3 個偏振。自由度守恆。

Kittel 在 p.112 把這句話寫成：「The number of modes with K ≤ K_D exhausts the number of degrees of freedom of a monatomic lattice.」——注意他明確限定 monatomic（單原子晶格，一個原胞一個原子）。若每個原胞有 s 個原子，總自由度是 3sN³，其中 3N³ 個屬於聲學支（Debye 模型處理的部分），剩下 3(s−1)N³ 個屬於光學支（Debye 模型完全不管，要用第 5 章後面的 Einstein 模型 p.114-116 去接）。

矽和鍺就是這種情況：金剛石結構，每個原胞 2 個原子，6 支聲子（3 聲學 + 3 光學）。這件事在讀 Fig. 8 的時候會回來——待會兒會算給你看。

所以請把這句話記牢：Debye 的截止 ω_D 不是「為了讓積分收斂而人為加上去的」，而是「K 空間本來就只有這麼多格點」這個事實的頻率表述。很多教科書（不含 Kittel）把 ω_D 描述成一個 cutoff parameter，那是本末倒置。

Fig. 6 在說什麼（網頁無圖，用文字還原）我的判讀

p.111 上半頁的 Fig. 6 是整段最值得看懂的一張圖，但它畫的是二維的情況（因為三維畫不出來）。圖說逐字如下：

「Allowed values in Fourier space of the phonon wavevector K for a square lattice of lattice constant a, with periodic boundary conditions applied over a square of side L = 10a. The uniform mode is marked with a cross. There is one allowed value of K per area (2π/10a)² = (2π/L)², so that within the circle of area πK² the smoothed number of allowed points is πK²(L/2π)².」

圖面內容：一個灰色的正方形，裡面規則排列著黑點。正方形代表第一布里淵區，圖頂端有一條標註 π/a 的線段，量的是從中心到右邊界的一半——實際上灰色方塊橫跨 −π/a 到 +π/a，邊長 2π/a。黑點就是式(17) 的二維版本允許的 K 值，相鄰兩點的間距是 2π/L = 2π/10a，所以沿每個方向恰好有 10 個點（對應 10 個原胞），整個方塊裡約 10 × 10 = 100 個點。正中央那個點（K = 0，均勻模態）依圖說標了一個叉。方塊中間疊了一個白色虛線圓，圓心在原點，一條白色箭頭從圓心指向圓周，箭頭旁標著 K——這個圓就是「所有波向量大小小於 K 的模態」的邊界。

這張圖要你看懂的三件事：

第一，格點是均勻的。不論靠近中心還是靠近邊界，黑點的疏密完全一樣。這是「數格點 = 算體積」這個換算能成立的唯一理由。（對照一下：如果用 p.109 Fig. 3 的固定端邊界條件，K 只能取正值，格點密度會差一倍，而且沒有 K = 0。週期性邊界條件讓 ±K 都存在、間距加倍，總數不變——這就是 p.110 說的「gives the same number of modes (one per mobile atom)」。）

第二，圓面積除以單格面積 = 圓內的點數。圖說最後那句 πK²(L/2π)² 就是二維版的式(19)。三維時圓變成球、面積變成體積，(L/2π)² 變成 (L/2π)³，就得到式(19)。Fig. 6 是式(19) 的二維預演。

第三，圓遲早會頂到方塊的邊。當 K 增大到接近 π/a，虛線圓就會切出灰色方塊之外——這時「圓內的點數」這個算法就開始失真，因為方塊外面根本沒有點。Debye 模型的最大破綻就在這裡：它假設一路到 K_D 都還是完整的球，但真實的第一布里淵區是多面體（矽是截角八面體），球和多面體在角落處對不上。Kittel 選 K_D 讓球的體積等於區的體積（下面式(24) 的推導會證明這一點），所以總模態數是對的，但分布是錯的——球把區角落的模態挪到了球面上。這正是 Debye 模型在中溫（T ~ 0.1θ 到 θ）誤差最大的根源。

D(ω) 到底是什麼：三個檢查點我的判讀

在往下走之前，先把 D(ω) 這個物件釘死，否則後面的積分只會變成符號操作。

檢查點一：它是密度，不是數量。D(ω) 本身沒有物理意義，D(ω)dω 才有——它是「頻率落在 ω 到 ω + dω 之間的模態有幾個」。所以 D(ω) 的單位是「個/(rad·s⁻¹)」= 秒。任何時候看到 D(ω) 單獨出現而沒有 dω 跟著，就要警覺。

檢查點二：它對偏振的計數要講清楚。Kittel 在式(20) 前明寫「for each polarization type」——式(20)、(22) 都是單一偏振的態密度。三維晶體每個 K 有 3 支聲學模態（1 縱 + 2 橫），所以總態密度是 3 倍。Kittel 在式(26) 前才補上這個因子，而且明說是為了省事：「For brevity we assume that the phonon velocity is independent of the polarization, so that we multiply by the factor 3 to obtain」。這是 Debye 模型的第二個近似（第一個是 ω = vK），而且比第一個更粗糙——真實晶體的縱波速度典型是橫波的 1.4 到 2 倍（矽約 1.44、銅約 2.05、鋁約 2.11；標準手冊值，非本頁讀到）。嚴謹的做法是對每個偏振用各自的 v_p 積分，最後定義一個平均速度 v_D 滿足 3/v_D³ = 1/v_L³ + 2/v_T³（這條式子 Kittel 沒寫，是本講解補充；注意是三次方倒數的平均，因為 D(ω) ∝ 1/v³）。

檢查點三：它必須歸一化到正確的模態總數。∫₀^{ω_D} D(ω)dω = N（單一偏振）。這一條在上面已經驗算過了。順帶可以算出 Debye 態密度在截止點的值：

D(ω_D) = Vω_D²/2π²v³，而 V/2π²v³ = 3N/ω_D³（由式(23)），所以 D(ω_D) = 3N/ω_D（單一偏振）。

也就是說，如果 D(ω) 是一條「高度固定」的水平線而不是拋物線，它只需要高度 N/ω_D 就能裝下 N 個模態；實際的拋物線在截止點的高度是那條水平線的 3 倍。這個「3」再一次提醒你：模態全都擠在高頻端。

Debye 近似的物理動機：為什麼可以把晶體當成連續彈性體

p.112 開新小節 Debye Model for Density of States，第一句話就是整個模型的全部內容：

「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.」

請把重點放在後半句：as it would be for a classical elastic continuum（就像古典彈性連續體那樣）。Debye 近似的本質不是「隨便取個線性色散」，而是把晶體徹底當成一塊連續的彈性介質——一塊沒有原子、只有密度 ρ 與彈性常數的果凍。在果凍裡，聲波的色散關係嚴格是 ω = vK，v 由彈性常數與密度決定，與波長無關。

這個近似什麼時候好？當波長遠大於原子間距時。因為波長很長的時候，波「看不見」原子的離散性，晶體對它而言確實就是連續介質。用波向量講：K ≪ π/a 時，Debye 近似是精確的，不是近似。第 4 章算單原子鏈時已經證明了這一點——p.92 的式(9) ω = 2(C/M)^{1/2}|sin(Ka/2)| 在 Ka ≪ 1 時展開成 p.94 的式(15) ω ≈ (C/M)^{1/2}aK，確實線性。

Kittel 在 p.114（下一節 Debye T³ Law）把這層意思講得最清楚：

「At sufficiently low temperature the T³ approximation is quite good; that is, when only long wavelength acoustic modes are thermally excited. These are just the modes that may be treated as an elastic continuum with macroscopic elastic constants.」

以及緊接著那句：「The energy of the short wavelength modes (for which this approximation fails) is too high for them to be populated significantly at low temperatures.」

這兩句話合起來就是理解整個 Debye 模型的鑰匙：低溫時只有長波模態被激發，而長波模態恰好是 Debye 近似精確的那些。

三個溫區：Debye 模型什麼時候準、什麼時候不準我的判讀

（以下三條分區是本講解由 p.112 與 p.114 的原文整理出來的判讀框架，不是 Kittel 的原話。）

• T ≪ θ：Debye 模型不只是「好」，它是漸近精確的（T³ 定律的係數 234Nk_B 是真的，不是擬合出來的）；
• T ≫ θ：Debye 模型也對，但理由完全不同——高溫時每個模態都拿到 k_BT，根本不在乎頻率分布長什麼樣，所以任何把模態數算對的模型都會給出 3Nk_B；
• T ~ 0.1θ 到 θ：Debye 模型最差，因為這時被激發的正是中高頻模態，而那裡的真實 D(ω) 與拋物線差最多。這也是為什麼實驗上「Debye 溫度」會隨溫度變化——你若在不同溫度區間用式(30) 去擬合同一塊材料，會得到不一樣的 θ。Table 1（p.116）的表頭因此特別註明取的是 low temperature limit of θ。

Debye 截止的邏輯：三式一氣呵成

p.112 中段是整節的樞紐，原文只有三句話加三條式子，但邏輯鏈非常緊。逐字引用：

「If there are N primitive cells in the specimen, the total number of acoustic phonon modes is N.」（式(23) 前）

「To this frequency there corresponds a cutoff wavevector in K space:」（式(24) 前）

「On the Debye model we do not allow modes of wavevector larger than K_D. The number of modes with K ≤ K_D exhausts the number of degrees of freedom of a monatomic lattice.」（式(24) 後）

把這三句翻成一條推理鏈：

前提：樣品有 N 個原胞 → 每個偏振恰有 N 個聲學模態（見前面「為什麼上界恰好是 Nπ/L」一節的自由度守恆論證）。

做法：既然 Debye 的 D(ω) = Vω²/2π²v³ 一路往上不會自己停，就人為地在 ω_D 處截斷，並要求截斷之後剛好裝下 N 個模態。用式子寫就是 ∫₀^{ω_D}D(ω)dω = N，或等價地把式(19) 的 K 換成 K_D、N 換成原胞數。

結果：式(23) ω_D³ = 6π²v³N/V。

換算：除以 v³ 再開立方，得式(24) K_D = ω_D/v = (6π²N/V)^{1/3}。

再換算：乘上 ℏ/k_B 換成溫度單位，得式(28) θ = (ℏv/k_B)(6π²N/V)^{1/3}。

三條式子是同一個內容的三種單位：ω_D 是頻率、K_D 是波向量、θ 是溫度。ℏ 與 k_B 只是換算因子。如果你覺得要記三條式子，那就是還沒看穿它們是同一條。

（上面「三條式子是同一個內容的三種單位」這個整理是本講解的說法，Kittel 只是依序寫出三式而未點破。）

最後留意 Kittel 那個小心翼翼的限定詞：monatomic lattice。他很清楚「N 個聲學模態用光全部自由度」這句話只對一個原胞一個原子的晶格成立。對矽、鍺這種一個原胞兩個原子的晶體，3N 個聲學模態只用掉了 6N 個自由度裡的一半，另外 3N 個屬於光學支。Debye 模型對矽鍺是不完整的——這一點在讀 Fig. 8 時會直接反映在數字上。

校註：Kittel 的 N 在這三頁裡有三個意思我的判讀

這是讀 p.111-113 最容易出錯的地方，而且錯了會一路錯到算不出正確數值。把三個 N 攤開：

第一個和第三個差一個立方（每邊 N 個 vs 總共 N³ 個）。第三個和第四個差一個「每原胞原子數 s」。第二個則根本是另一回事——它是一個隨 K 變動的函數 N(K)，只是碰巧在 K = K_D 時等於第三個 N。

為什麼 Kittel 在式(29) 突然改口說 N 是原子數？因為他從式(23) 到式(30) 一路都隱含著 monatomic lattice（他在式(24) 後明講了）。對單原子晶格，原胞數 = 原子數，兩個說法一致，所以他就換了個比較直觀的說法。這在教學上可以理解，但對矽鍺這種一原胞兩原子的材料會造成實質的計算歧義。

實務上怎麼辦？建議固定用以下版本，並在自己的筆記裡註明：

• 式(23)(24)(28) 裡的 N/V 永遠是原胞數密度（Debye 只管 3 支聲學模態）。
• 式(29)(30) 高溫極限給出的 3Nk_B，這個 N 應該是原子數（Dulong–Petit 定律數的是原子自由度）。

對矽這兩者差 2 倍。所以嚴格來說，Debye 模型對矽最多只能解釋 3N_cell k_B = 1.5 N_atom k_B 的熱容，剩下一半要靠光學支（Kittel 隨後在 p.114-116 引入 Einstein 模型正是為此）。

⚠️ 另一個校註（關於式(30) 的印刷）：本 PDF 的式(30) 中間式，積分變數與指數裡印的是斜體 w（∫₀^{w_D} dw ω⁴e^{ℏw/τ}/(e^{ℏw/τ} − 1)²），只有 ω⁴ 那一項用了 ω。經把該頁算繪成影像放大目視確認：這是書本第 8 版本身的排版瑕疵，不是 PDF 文字層的抽取錯誤。全式應一律讀作 ω。本講解在引用式(30) 時已統一改回 ω，特此註明依據。

θ 沒有可調參數：這決定了它能不能當材料指紋我的判讀

這一節要把式(28) 的意義講透，因為它是整章最有實用價值的一條式子。

先看清楚式(28) 裡有什麼、沒有什麼：

有：ℏ、k_B（普世常數）、v（聲速）、N/V（數密度）。

沒有：溫度、樣品尺寸、樣品形狀、雜質濃度、缺陷密度、任何擬合參數。

這意味著三件事。

第一，θ 是可以「算」出來的，不是只能「擬」出來的。給我一塊材料的彈性常數與密度，我就能在不做任何熱容量測的情況下預測它的 θ，進而預測整條 C(T) 曲線。這是 Debye 模型相對於 Einstein 模型的決定性優勢——Einstein 模型的 θ_E 沒有獨立的計算管道，只能回頭去 fit 熱容資料，所以它解釋數據但不預測數據。

第二，θ 是一個真正的材料常數，可以拿來比對與辨識。Table 1（p.116，標題「Debye temperature and thermal conductivity」，θ 那一列的表頭註明是 low temperature limit of θ, in Kelvin）就是拿它當元素週期表的一個欄位在用。θ 高代表「硬且密」，θ 低代表「軟且疏」。

第三，也是最要緊的：正因為沒有可調參數，θ 對不上就是模型錯了。實務上大家確實會發現「Debye 溫度隨溫度變化」——同一塊矽，用 10 K 附近的數據 fit 得到一個 θ，用 300 K 附近的數據 fit 得到另一個。這不是 θ 真的在變，而是式(30) 這個曲線形狀本身不對，你硬要用一個參數去 fit 一條形狀不對的曲線，得到的參數自然隨你 fit 的區間而變。Table 1 特別標明取「low temperature limit」，正是因為只有在 T ≪ θ 時 Debye 的形狀才是漸近正確的，那裡 fit 出來的 θ 才有唯一而確定的意義。

θ(T) 的變化本身是有用的診斷訊號。如果你把量到的 C(T) 逐點反解成 θ(T)（每個溫度點都解一次式(30)），畫出來的 θ(T) 曲線就等於在告訴你真實 D(ω) 偏離拋物線的方式：θ(T) 在某溫度下凹，代表那個頻段的真實模態數比拋物線多（真實 D(ω) 有 van Hove 峰）。這是一種不需要中子散射就能窺見聲子譜結構的窮人做法，在文獻上叫做 θ(T) plot，Kittel 在這幾頁沒有介紹（這一段是本講解補充）。

最後一句實用的話：做熱模擬時，別直接從文獻抄一個 θ 就用。先確認那個 θ 是從哪個溫區、用什麼 N 的約定（原胞還是原子）擬合出來的。前面「N 有三個意思」一節已經指出，矽用原子密度或原胞密度會讓反解的聲速差 26%。

Fig. 7 在說什麼（網頁無圖，用文字還原並補上數值）我的判讀

p.113 上半頁的 Fig. 7 是式(30) 的圖解。圖說逐字如下：

「Heat capacity C_V of a solid, according to the Debye approximation. The vertical scale is in J mol⁻¹ K⁻¹. The horizontal scale is the temperature normalized to the Debye temperature θ. The region of the T³ law is below 0.1θ. The asymptotic value at high values of T/θ is 24.943 J mol⁻¹ deg⁻¹.」

圖面內容：橫軸標 T/θ（帶一個向右的箭頭），從 0 到 1.6，主要刻度每 0.2 一格，格線是直的細線；縱軸標 Heat capacity, in J mol⁻¹ K⁻¹，從 0 到 25，刻度每 5 一格。畫面上只有一條粗實線。曲線從原點出發，一開始幾乎貼著橫軸（這是 T³ 區，掉得極快），約在 T/θ = 0.1 開始離開軸線，在 0.15 到 0.35 之間以近乎直線的陡峭斜率上衝，之後開始明顯彎折，過了 0.6 就趨於平緩，在 1.0 以後幾乎水平，緩慢逼近 24.943 這條漸近線但在圖的右端（T/θ = 1.6）仍未真正到達。

關鍵讀值（本講解對式(30) 數值積分所得，可與圖上曲線逐點對照）：

三個要記住的數字：

24.943 J mol⁻¹deg⁻¹ = 3R。驗算：3 × 8.3145 = 24.943 ✔。這是每莫耳原子的古典熱容，與材料完全無關——這就是 Dulong–Petit 定律。（圖說用的是 deg⁻¹ 這個舊寫法，與 K⁻¹ 同義。）

0.1θ 是 T³ 區的上界。驗算精度：T³ 公式 C = 234Nk_B(T/θ)³ 在 T = 0.1θ 給 1.946，精確值 1.891，高估 2.9%；在 T = 0.05θ 給 0.2432 vs 精確 0.2430，誤差 0.1%；但在 T = 0.2θ 給 15.56 vs 精確 9.20，高估 69%。圖說寫 0.1θ 是相當寬鬆的說法——要拿 T³ 定律做定量分析，實務上應該用 0.05θ 以下。Kittel 自己在 p.114 講得更嚴：對實際晶體，可能必須低到 T = θ/50（= 0.02θ）才能得到相當純粹的 T³ 行為。（此句原文中的「T = θ/50」在本 PDF 的文字層漂到句外，故此處採轉述而非逐字引用。）

T = θ 只達到 95%。這一點反直覺。很多人以為 T = θ 就「飽和了」，其實還差 5%，要到 T ≈ 3θ 才到 99.4%。「高溫極限」在 Debye 模型裡來得比想像中慢。

Fig. 8 在說什麼：矽與鍺的實測熱容（並用 θ 驗證 Debye 的普適曲線）我的判讀

p.113 下半頁的 Fig. 8 是本段唯一的實驗圖，也是與你的 Cth 研究最直接相關的一張。圖說逐字如下：

「Heat capacity of silicon and germanium. Note the decrease at low temperatures. To convert a value in cal/mol-K to J/mol-K, multiply by 4.186.」

圖面內容：橫軸標 Temperature, K，從 0 到 300 出頭，主要格線在 100、200、300；縱軸標 C_p, in cal/mol-K（注意是 C_p，定壓熱容，不是理論的 C_V），從 0 到 6，每 1 一格。圖上兩條實線，各自標了名字：上面那條標 Germanium，下面那條標 Silicon。

兩條曲線都從原點出發，起步時幾乎貼著橫軸（這就是圖說要你注意的 "the decrease at low temperatures"），然後向上彎起。鍺的曲線始終在矽的上方，而且上升得早得多：鍺約在 50 K 就明顯抬頭，在 100 K 已經到 3.3 左右，200 K 約 5.0，300 K 約 5.5 並已明顯趨於平緩；矽則要到約 40 K 起才明顯抬頭、70 K 左右才通過 1 cal/mol-K，100 K 時約 1.7，200 K 約 3.8，300 K 約 4.7，在圖的右端仍在上升、尚未飽和。

這張圖的第一個重點：兩條曲線是同一條曲線。Fig. 7 說 Debye 熱容只是 T/θ 的普適函數。把 Table 1（p.116）的 θ 代進去驗算——

兩者都對到小數點後一位。（換算用圖說給的 1 cal/mol-K = 4.186 J/mol-K；3R = 24.943 J/mol-K = 5.958 cal/mol-K，這就是兩條曲線最終要趨近的水平線，但圖的縱軸只到 6，所以那條漸近線幾乎就貼在圖框頂端。）

這件事本身就是 Debye 模型最漂亮的驗證：矽和鍺是同一種晶體結構、化學性質相近，唯一的重大差別是鍺的原子重約 2.6 倍、鍵較軟，所以聲速低、θ 低。θ 低代表「用同一個絕對溫度去量，鍺已經走到曲線更右邊的位置」，所以鍺的熱容看起來總是比矽高。不是鍺「比較會吸熱」，是鍺的聲子模態比較容易被激發。

第二個重點：為什麼是 C_p 不是 C_V？Fig. 8 畫的是實驗量到的定壓熱容，理論式(30) 算的是定容熱容。Kittel 在 p.107 的註腳 1 處理了這個落差：「A thermodynamic relation gives C_p − C_V = 9α²BVT, where α is the temperature coefficient of linear expansion, V the volume, and B the bulk modulus. The fractional difference between C_p and C_V is usually small in solids and often may be neglected.」對矽而言，α ≈ 2.6×10⁻⁶ K⁻¹ 很小，而且 C_p − C_V ∝ α²T，在 300 K 以下只有百分之零點幾——所以上表能對得那麼準。但對熱膨脹係數大的材料（例如某些聚合物或軟金屬），這個差別不能忽略。

第三個重點（這是本段對 Cth 研究最重要的一句話）：把 Fig. 8 的縱軸從 cal/mol-K 換成工程上用的 J cm⁻³K⁻¹，就得到元件熱模擬需要的體積熱容——下面有專門的一節處理這個換算。

「低溫下降」到底在下降什麼：模態凍結的機制我的判讀

Fig. 8 圖說那句 「Note the decrease at low temperatures.」只有六個字，但它是整章最需要講透的一句。

古典物理預測熱容是常數。能量均分定理說每個自由度拿 ½k_BT，晶體有 3N 個振動自由度（每個 ½k_BT 動能 + ½k_BT 位能 = k_BT），所以 U = 3Nk_BT，C = 3Nk_B，與溫度無關。這就是 Dulong–Petit 定律，它在 19 世紀初就被實驗證實——在室溫附近。

但低溫時熱容明明就掉下去了，而且掉得比線性還快（∝ T³）。古典物理對此完全無解。這是十九世紀末「古典物理的烏雲」之一，也是 Einstein 在 1907 年、Debye 在 1912 年用量子論解決的第一批固態問題。

量子解答的機制，可以完全用這一段的式子講清楚。

關鍵在 Planck 因子 ⟨n(ω)⟩ = 1/(e^{ℏω/k_BT} − 1)。一個頻率 ω 的模態，平均能量是 ⟨n⟩ℏω：

• 當 ℏω ≪ k_BT（即 x ≪ 1，「軟模態」）：展開 e^x − 1 ≈ x，得 ⟨n⟩ ≈ k_BT/ℏω，所以 ⟨n⟩ℏω ≈ k_BT——完全古典，這個模態貢獻 k_B 的熱容。
• 當 ℏω ≫ k_BT（即 x ≫ 1，「硬模態」）：⟨n⟩ ≈ e^{−ℏω/k_BT}，能量 ≈ ℏω·e^{−ℏω/k_BT}，指數地趨近於零。這個模態幾乎完全不參與熱交換——它被凍結（frozen out）了。

凍結的物理圖像：量子諧振子的能階是離散的，間距 ℏω。要激發它，至少要一整份 ℏω 的能量。如果環境能提供的典型能量 k_BT 遠小於 ℏω，那麼「湊足一整份」的機率是 Boltzmann 抑制的 e^{−ℏω/k_BT}，小得可憐。這個模態就像被鎖住一樣，溫度再怎麼變它都不吸熱，所以對熱容沒有貢獻。

把這件事和 D(ω) ∝ ω² 合起來，就得到完整的低溫圖像（Kittel 在 p.114 用 Fig. 10 的圖像論證講了同一件事）：

1. 溫度 T 定義了一個「熱波向量」K_T，滿足 ℏvK_T = k_BT（p.114 原文的定義）。只有 K < K_T 的模態是活的，其餘全部凍結。
2. 活著的模態佔 K 空間的比例是 (K_T/K_D)³ = (T/θ)³——這就是三維的功勞，球體積 ∝ 半徑三次方。
3. 活著的模態共 3N(T/θ)³ 個，每個貢獻 k_BT 的能量，所以 U ≈ 3Nk_BT(T/θ)³。
4. 微分得 C_V ≈ 12Nk_B(T/θ)³ ∝ T³。（p.114 本文給的正是這個 12Nk_B，Fig. 10 的圖說——印在 p.115——也寫同一條；精確係數是 234Nk_B，差別來自「活/死」的階梯近似太粗糙，但冪次完全正確。）

所以「低溫下降」= 高頻聲子模態被逐一凍結。溫度每降一半，活著的模態數就掉到 1/8。

這對熱容量測意味著什麼：熱容不是材料的固有常數，而是「在這個溫度下有多少個振動自由度還醒著」的計數器。你在 300 K 量到的矽熱容（4.7 cal/mol-K，只有古典值的 79%）並不是「矽的熱容」，而是「矽在 T/θ = 0.465 這個位置上的熱容」。如果元件在 400 K 工作，T/θ = 0.62，熱容會升到 21.9 J/mol-K（88%）——同一塊矽，熱容差了 10%。做暫態熱分析時把 Cth 當常數，在這個溫區會引入這個量級的誤差。

五個常見誤解（含每一個的反例或驗算）我的判讀

誤解一：「Debye 溫度是一個相變溫度，材料在 T = θ 時會發生什麼事。」

不會。θ 只是 ℏω_D/k_B，是把「最高聲子頻率」換算成溫度單位。Fig. 7 在 T/θ = 1 附近完全光滑，沒有任何奇異性。θ 的唯一作用是分隔量子區與古典區的尺度，就像 Fermi 溫度 T_F 之於電子氣（Kittel 在 p.145 特別警告過同一件事：「Recall that although T_F is called the Fermi temperature, it is not the electron temperature, but only a convenient reference notation.」）。驗算：T = θ 時 C_V = 23.74 J/mol·K（95% of 3R），T = 1.1θ 時 23.97——毫無異常。

誤解二：「T³ 定律是 Debye 近似的產物，所以是近似的。」

正好相反。T³ 定律在 T → 0 是漸近精確的，與 Debye 近似的粗糙無關。理由在前面「Debye 近似的物理動機」一節已經給過：低溫時只有長波聲學模態被激發，而長波模態的色散本來就是線性的、D(ω) 本來就是 ∝ ω²。Kittel 在 p.114 講得很清楚：「These are just the modes that may be treated as an elastic continuum with macroscopic elastic constants.」所以 T³ 定律不需要 Debye 模型背書；反而是 Debye 模型因為在低溫恰好正確而獲得信任。

⚠️ 但要加一個限定：T³ 定律的前提是「三維、有長程週期性的晶體」。Kittel 在 p.580（第 19 章 NONCRYSTALLINE SOLIDS，不是第 18 章）的 Fig. 8 給了反例——熔融石英（vitreous silica）這類非晶固體，低溫熱容在 1 K 以下大致線性於 T，該圖並用虛線畫出計算的 Debye 熱容作對照，兩者明顯分岔。非晶固體沒有本段賴以成立的週期性，低頻端多出額外的激發態。所以嚴格的說法是：T³ 定律對三維晶體在 T → 0 漸近精確，對非晶固體不成立。

誤解三：「ω_D 是為了讓積分收斂而人為引入的截斷參數，可以當作 fitting parameter。」

不是。ω_D 由式(23) 完全決定，內容是「N 個原胞恰好有 N 個聲學模態」這條自由度守恆律。它沒有調整餘地。若你把 ω_D 當可調參數去 fit，你就不是在做 Debye 模型，而是在做一個兩參數的唯象擬合。（真正的自由度在於：你要用哪個 v？縱波、橫波，還是 3/v_D³ = 1/v_L³ + 2/v_T³ 的 Debye 平均？這個選擇會影響 θ，但它是物理選擇，不是擬合。）

誤解四：「Debye 模型把 D(ω) 畫成拋物線，所以低頻模態最多。」

完全相反。D(ω) ∝ ω² 代表高頻模態最多：低於 ω_D/2 的模態只佔 12.5%，高於 0.8ω_D 的佔 48.8%（前面已算過）。拋物線的形狀容易讓人誤讀——它在低頻貼近橫軸，在高頻才衝高。這個「頭重腳輕」的分布正是熱容上升得那麼慢、要到 T ≈ 3θ 才飽和的原因。

誤解五：「Debye 模型對矽、鍺的熱容給出完整的描述。」

不完整。矽鍺是金剛石結構，一個原胞兩個原子，所以有 6 支聲子：3 支聲學 + 3 支光學。Debye 模型只處理聲學支。Kittel 在式(24) 後的限定詞已經預告了這件事：「exhausts the number of degrees of freedom of a monatomic lattice」。嚴格說，矽的 Debye 部分最多只能解釋 3N_cell k_B = 1.5 N_atom k_B，另一半要靠光學支（Kittel 隨後在 p.114 引入 Einstein 模型 D(ω) = Nδ(ω − ω₀) 正是為此，所以標準做法是 Debye + Einstein 混合）。那為什麼 Fig. 8 的數字用單一 θ 就對得那麼準？是因為約定一致，不是因為光學支被 θ 吸收了。Table 1 的 θ 依慣例以「每原子 3 個模態」定義（N = 原子數），而前面式(29)(30) 的所有數值也一律用 N = 原子數（3R = 24.943 J mol⁻¹K⁻¹ 是每莫耳原子的古典上限），兩邊約定相同，所以對得上。若改用原胞數密度去算 θ，同一組彈性常數給出的 θ 會從約 696 K 掉到約 553 K（取矽的 Debye 平均聲速 v_D ≈ 6.35 km/s；彈性常數為標準手冊值，非本書讀到），Fig. 8 就完全對不上。

光學支的影響必須另外談。Table 1 的表頭明寫取的是 low temperature limit of θ，而低溫下矽的光學支（Γ 點約 15.5 THz，ℏω/k_B ≈ 745 K）是完全凍結的，對 T³ 係數毫無貢獻——低溫極限的 θ 由聲學支的聲速決定，不可能「吸收」一個在那個溫區根本沒醒的分支（這一點與前面「θ 沒有可調參數」一節強調 Table 1 取低溫極限的理由是同一件事）。光學支真正的後果，是讓真實的 C(T) 在中高溫偏離 Debye 形狀，也是 θ(T) 會隨擬合溫區漂移的原因之一。用 N = 原子數的 Debye 模型，實質上是把 3s 支全部塞進同一條拋物線的單參數近似：它保證高溫回到正確的 3N_atom k_B，但不代表光學支的形狀被描述對了。

實務上仍要區分兩種 θ：熱容擬合值與由彈性常數計算的值即使約定相同也未必相等（矽：擬合 645 K vs 彈性常數算出的約 696 K），引用前務必看清楚是哪一種、用的又是原子還是原胞約定。

接到 Cth：從 Fig. 8 的 cal/mol-K 到元件熱容延伸知識

⚠️ 本節是本講解補充的跨領域連結。經對 Kittel 第 8 版全書逐詞檢索確認：Kapitza resistance、thermal boundary resistance、thermal time constant、thermal impedance、self-heating、Cauer/Foster 網路等元件層級暫態熱分析的概念，在全書出現次數為零。所以以下「從材料熱容到元件 Cth」的橋接，書上完全沒有，但用到的每一個材料量都來自這幾頁。

第一步：把單位換成工程用的體積熱容。

Fig. 8 的縱軸是 cal/mol-K，熱模擬需要的是 J cm⁻³K⁻¹（因為 Cth = C_vol × 體積）。換算需要莫耳體積 V_m = N_A/n：

（原子濃度取自 Table 4，書本 p.21，影像頁目視判讀；莫耳熱容取自前面對式(30) 的數值積分，並已用 Fig. 8 的圖上讀值交叉驗證。矽的 1.66 J cm⁻³K⁻¹ 與工程手冊常用的 1.63-1.66 一致。）

值得注意的巧合：矽和鍺的莫耳熱容差很多（20.0 vs 23.1，差 15%），但體積熱容幾乎一樣（1.66 vs 1.70，差 2%）。原因是鍺的原子較大（莫耳體積大 13%），恰好抵消了它較高的莫耳熱容。這其實是 Dulong-Petit 的一個推論：在古典極限下，體積熱容 = 3nk_B，只依賴原子數密度。取矽：3 × 5.00×10²² × 1.381×10⁻²³ = 2.07 J cm⁻³K⁻¹（這是古典上限；300 K 只走到 1.66，即 80%）。

第二步：接到 Kittel 自己的熱導率式子。

第 5 章 p.122 的式(42) 是 K = ⅓Cvℓ，Kittel 明講「where C is the heat capacity per unit volume, v is the average particle velocity, and ℓ is the mean free path of a particle between collisions」——這裡的 C 就是上表的 C_vol，正是本段算出來的量。（順帶一提，p.122 的 Table 2 也把 C 列成 J cm⁻³K⁻¹，單位與上表一致。）所以本段的態密度機器，在書內的直接下游就是熱導率。

第三步：元件層級的 Cth（書上完全沒有的部分）。

若把 FinFET 的鰭視為一塊矽，Cth = C_vol × V_fin。取一個示意的鰭尺寸 10 nm（寬）× 30 nm（高）× 100 nm（長）= 3×10⁴ nm³ = 3×10⁻¹⁷ cm³：

Cth ≈ 1.66 × 3×10⁻¹⁷ ≈ 5×10⁻¹⁷ J/K。

這個數量級提供了一個有用的下界：τ = R_th C_th，若把 R_th 取在 10³–10⁵ K/W 的範圍，單就鰭本體而言 τ 落在 5×10⁻¹⁴ 到 5×10⁻¹² s——也就是次皮秒到皮秒量級。實測的 FinFET 自熱時間常數之所以落在 ns 量級，是因為散熱路徑上串著體積大得多的熱容（STI、埋氧層、基板、金屬內連線），那些才是暫態響應的主導項。這正好說明為什麼 Cth 萃取必須看整條熱路徑，而不是只算通道矽的體積熱容。

第四步：本段的物理對 Cth 萃取的四個實際提醒。

（1）Cth 隨溫度變，而且在元件工作溫區變得最兇。由式(30) 算出的矽體積熱容：

300 → 400 K 之間 C_vol 上升約 10%，而這正是 Fig. 7 曲線最彎的那一段。把 Cth 當常數做暫態擬合，在 300-400 K 的溫升範圍內會引入約 10% 的系統性誤差，而且這個誤差與溫升相關（溫升愈大偏差愈大），在資料上會偽裝成「非線性熱阻」。

（2）Cth 是「醒著的自由度」的計數器，不是材料的固有常數。這是前面「低溫下降」那一節的結論。低溫量測拿到的 Cth 會遠小於室溫值——上表的 77 K 那一列：T/θ = 0.119，C_vol 只有 0.254 J cm⁻³K⁻¹，是室溫值的 15%。若某個特性化流程在低溫下做、卻拿室溫的 C_vol 去反推熱阻，誤差會是數倍而不是數個百分點。

（3）奈米尺度可能改變 D(ω) 本身，不只是改變體積。式(17) 的 K 空間格點間距是 2π/L。當鰭寬 L = 10 nm 時，2π/L = 6.3×10⁸ m⁻¹，只比 K_D = 1.44×10¹⁰ m⁻¹ 小 23 倍——橫向每個方向的允許 K 點只剩 L/a ≈ 18 個，離散性開始顯現，而且表面會引入額外的聲子模態。此時「用塊材 C_vol 乘體積」這個做法的前提就要打問號。Kittel 第 18 章（p.557 起）談的正是這個尺度的熱性質，其中 p.561-562 給了一維彈道熱導 G_th 的量子極限。

（4）Debye 只管聲學支，但矽有光學支。前面誤解五已經指出：矽一個原胞兩個原子，光學支承擔一半的自由度。室溫下矽的光學聲子（Γ 點約 15.5 THz，換算 ℏω/k_B ≈ 745 K；標準值，非本書讀到）大部分仍是凍結的，所以它們對 300 K 的 Cth 貢獻有限；但在 400 K 以上就開始醒來，這是 C_vol(T) 在高溫仍持續上升的部分原因。

這一段在全書的位置：後面哪幾頁會回來用它我的判讀

把 p.111-113 的成果放回全書的脈絡，可以看出它是一個被反覆調用的工具庫，而不是一個孤立的章節。以下每一條都是本講解為這幾頁查證過原文之後列出的下游（頁碼均為書本頁碼）。

（一）緊接著的 Debye T³ Law（p.114-115） ——直接用式(29)。把積分上限取到無窮，用 ∫₀^∞x³/(e^x − 1)dx = π⁴/15（式(31)，該頁給了完整推導），得到 U ≈ 3π⁴Nk_BT⁴/5θ³，微分得式(32) C_V = (12π⁴/5)Nk_B(T/θ)³ = 234Nk_B(T/θ)³。Fig. 9（p.115）用固態氬的數據驗證，圖說註明 θ = 92.0 K（與 Table 1 的 Ar 欄位 92 一致），並把縱軸的熱容對橫軸的 T³ 畫成直線。Fig. 10（p.115）則用「K_T 球 vs K_D 球」的體積比 (K_T/K_D)³ = (T/θ)³ 給出本段前面複述過的直觀論證。p.114 還有一句實務上比 Fig. 7 圖說更嚴格的提醒：可能必須低到 T = θ/50 才能得到相當純粹的 T³ 行為。

（二）Einstein Model（p.114 起） ——本段的對照組。Einstein 取 D(ω) = Nδ(ω − ω₀)，把 Debye 的拋物線換成一根 δ 函數。兩者在高溫都給 3Nk_B（因為模態總數相同），但低溫差天差地遠（Einstein 給指數衰減、Debye 給 T³）。兩相對照才看得出 D(ω) 的形狀只在低溫與中溫才有影響——這是前面「高溫極限」一節的結論。實務上矽鍺這類多原子原胞材料，標準做法就是 Debye（聲學支）+ Einstein（光學支）的混合模型。

（三）Thermal Conductivity（p.121-128） ——熱容在這裡從「主角」變成「因子」。式(42) K = ⅓Cvℓ 裡的 C 就是本段的體積熱容（p.122 明寫 heat capacity per unit volume，Table 2 也用 J cm⁻³K⁻¹）。Kittel 在 p.127 直接把本段的結論接上去：「The only temperature-dependent term on the right is C, the heat capacity, which varies as T³ at low temperatures. We expect the thermal conductivity to vary as T³ at low temperatures.」——低溫熱導率的 T³ 完全來自熱容的 T³。Fig. 19 的鍺同位素實驗在 5 K 以下量到 K = 0.06T³，正是這條邏輯的實驗確認。

（四）第 5 章習題（p.128） ——四題裡有三題直接用本段的機器：習題 1（Singularity in density of states）要你算真實色散下的 D(ω)，直接打臉 Debye 的拋物線；習題 3（Zero point lattice displacement and strain）用 ∫₀^{ω_D}D(ω)(…)dω 算 ⟨R²⟩ = 3ℏω_D²/8π²ρv³，與前面推導的零點能同構；習題 4（Heat capacity of layer lattice）要你證明二維受限的晶格熱容 ∝ T²，考的正是前面那張「D(ω) ∝ ω^{d−1}」維度表。

（五）第 6 章 自由電子費米氣（p.138、p.145） ——兩個接點，而且是全書最漂亮的一組平行。

第一，p.138 把式(18) 的計數整套搬去數電子。原文：「From (10) we see that there is one allowed wavevector—that is, one distinct triplet of quantum numbers k_x, k_y, k_z—for the volume element (2π/L)³ of k space.」——與本段式(18) 前那句「there is one allowed value of K per volume (2π/L)³ in K space」幾乎逐字相同。接著式(15) 寫成 2·(4πk_F³/3)/(2π/L)³ = Vk_F³/3π² = N，「where the factor 2 on the left comes from the two allowed values of the spin quantum number for each allowed value of k」，於是式(16) 給出 k_F = (3π²N/V)^{1/3}。把它和本段的式(24) K_D = (6π²N/V)^{1/3} 並排：

兩者只差一個因子 2，而那個 2 就是自旋。同一套 V/8π³ 幾何，換一個簡併度就從聲子跨到電子。

第二，p.145 的 Experimental Heat Capacity of Metals：「At temperatures much below both the Debye temperature θ and the Fermi temperature T_F, the heat capacity of metals may be written as the sum of electron and phonon contributions: C = γT + AT³, where γ and A are constants characteristic of the material.」——那個 AT³ 就是本段的產物。而標準的實驗手法（p.145 式(37) C/T = γ + AT²，畫 C/T 對 T² 的直線，截距給 γ、斜率給 A；Fig. 9 用鉀的數據示範，擬合線是 C/T = 2.08 + 2.57T²）之所以可行，正是因為聲子項的冪次被本段釘死在 3、沒有任何自由度。若 T³ 是可調的冪次，這條直線法就垮了。

（六）第 10 章 超導（p.269-270、p.278） ——Debye 溫度直接進入 T_c。p.278 的 BCS 結果是 T_c = 1.14θ exp[−1/UD(ε_F)]，原文緊接著說明「where θ is the Debye temperature and U is an attractive interaction」。而 p.269-270 的 Isotope Effect 小節把這件事與同位素連起來：「The original BCS model gave the result T_c ∝ θ_Debye ∝ M^{−1/2}」——θ 與原子質量的這個 M^{−1/2} 關係，正是本段式(28) 的推論（θ ∝ v，而彈性連續體的聲速 v = (彈性常數/ρ)^{1/2} ∝ M^{−1/2}）。Kittel 同時給了誠實的但書：「but the inclusion of coulomb interactions between the electrons changes the relation. Nothing is sacred about α = ½.」

（七）第 18 章 奈米結構（p.557-562）與第 19 章 非晶固體（p.578-580） ——兩個本段會失效的場景。其一，當樣品尺寸接近 2π/K_D 時，本段假設的「K 空間連續」失效，p.561-562 給了一維通道的彈道熱導量子。其二，p.580 的 Fig. 8 畫了熔融石英（vitreous silica）的低溫熱容，圖說指出它在 1 K 以下大致線性於 T，並用虛線畫出「calculated Debye heat capacity of vitreous silica」作對照——非晶固體不遵守 T³ 定律，因為它沒有本段賴以成立的長程週期性，低頻端多出 Debye 模型無法產生的雙能階隧穿態。p.580 同時把第 5 章的熱導率式子照搬過來：「As in Chapter 5, the expression for the thermal conductivity K has the form (29) where c is the heat capacity per unit volume, v is an average phonon velocity, and ℓ is the phonon mean free path.」

一句話總結這一段的地位：p.111-113 只做了一件事——把「有多少個振動模態、它們的頻率怎麼分布」算清楚。但這一件事是後面所有熱性質（熱容、熱導、熱膨脹）、以及第 6 章電子熱容、第 10 章超導的共同地基。