03聲子熱容：Debye T³ 定律、Einstein 模型、D(ω) 的一般結果

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

到 p.113 為止，Debye 模型已經完整建立：態密度 (22) D(ω) = Vω²/2π²v³、截止頻率 (23) ω_D³ = 6π²v³N/V、Debye 溫度 (28) θ = (ℏv/k_B)(6π²N/V)^{1/3}，熱能 (29) 與熱容 (30) 都寫成了一個無量綱積分乘上前因子。問題是：(29) 與 (30) 的積分在一般溫度下沒有初等閉形式。Fig 7 那條曲線是數值算出來的，不是解析式。

這對教科書是個尷尬的處境——你有一個漂亮的模型，卻沒有一條可以背、可以代、可以拿去和實驗比的公式。這一段做的三件事，正好把這個缺口補起來：

（一）在低溫極限把積分做掉。 T ≪ θ 時上限 x_D = θ/T 趨近無窮，而 ∫₀^∞ x³/(eˣ−1)dx 是一個標準的解析積分，值為 π⁴/15。於是 C_V ≅ 234 Nk_B(T/θ)³ 這條可以背的公式出現了。這就是 Debye T³ 定律，也是全書被引用次數最多的結果之一：任何一份低溫比熱的實驗數據，第一件事就是畫 C/T 對 T² 圖或 C 對 T³ 圖，從斜率反推 θ。

（二）用一個刻意做壞的模型當對照組。 Einstein (1907) 模型假設所有的模都是同一個頻率 ω₀。Kittel 的原始設定是一維 N 個同頻振子，給 D(ω) = Nδ(ω−ω₀)；三維時（p.117）N 換成 3N，也就是 D(ω) = 3Nδ(ω−ω₀)。它在高溫給出正確的 Dulong–Petit 值 3Nk_B，在低溫卻給出指數衰減而不是 T³。Kittel 把它放在 Debye 之後講，不是按歷史順序（Einstein 1907 早於 Debye 1912），而是因為要先知道正確答案是 T³，才看得出 Einstein 模型錯在哪裡、為什麼錯。而且這個「錯」是有用的：光學支的譜本來就窄而平，Einstein 模型今天仍然是處理光學聲子那一段最省事的工具。

（三）把態密度從 Debye 的特例解放出來。 (22) 是在 ω = vK 這個強假設下算的，等頻面剛好是球面、群速度剛好是常數，所以積分平凡。真實晶體的等頻面不是球面（見第 4 章的色散曲線），群速度在區邊界甚至會歸零。式 (35)–(37) 給出不依賴任何色散假設的一般式 D(ω) = (V/(2π)³)∫dS_ω/v_g，代價是要做一個等頻面上的面積分。這個式子的價值不只在聲子——Kittel 在 p.119 明說：「We can use this result also in electron band theory.」後面幾章講電子態密度時用的就是它，只是把 ω 換成 ε、把 v_g 換成 |∇_kε|。

順帶一提，(37) 立刻告訴我們一件 Debye 模型永遠看不到的事：只要某處 v_g = 0，被積函數就發散，D(ω) 在那個頻率附近會出現結構。這就是 Van Hove 奇異點（Fig 14），是真實聲子譜與 Debye 的平滑 ω² 之間最明顯的差別。

T³ 律成立的物理條件：只有彈性連續體那一段才算數

Kittel 在式 (32) 之後用一整段解釋 T³ 為什麼只在很低的溫度成立。原文逐字是：

「At sufficiently low temperature the T 3 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.」

這段話有三個層次，值得分開讀：

第一層：什麼模被激發。 低溫下只有長波聲學模（long wavelength acoustic modes）被熱激發。這是 Planck 分布的直接後果——ℏω ≫ k_BT 的模佔據數是指數小的。

第二層：為什麼那些模剛好符合 Debye 的假設。 長波模的波長遠大於晶格常數，所以它們「看不見」原子的離散性，晶體對它們而言就是一塊均勻的彈性連續體（elastic continuum），色散關係嚴格是 ω = vK，而且 v 就是用巨觀彈性常數（macroscopic elastic constants）算出來的聲速。（以下為我的判讀）Debye 模型對這些模不是近似，是精確的。這一點非常關鍵：Debye 模型不是「一個粗糙的內插」，它在低溫是漸近精確的，因為在低溫唯一有貢獻的那群模剛好完全滿足它的假設。

第三層：為什麼模型的錯誤部分不傷害結果。 Debye 模型對短波模的描述是錯的（真實色散彎曲、有截止、有光學支），但這些模「too high for them to be populated significantly at low temperatures」——錯誤被 Boltzmann 因子指數壓制掉了。模型錯在一個不影響答案的地方。

（以下為我的判讀）這是整個 Debye 理論的正當性論證，而不只是一段補充說明。同樣的邏輯結構會在後面反覆出現：一個模型只要在「唯一有貢獻的那個區域」是對的，它在別處錯得多離譜都無所謂。第 6 章的自由電子氣之所以在金屬裡好用，用的是完全一樣的論證——只有費米面附近 k_BT 寬度內的電子有貢獻。

Fig 10 的手算論證：不做任何積分，怎麼看出 T³

Kittel 在 p.114 給了一個不用積分的論證，配圖是 Fig 10（書本 p.115）。因為網頁不放圖，先把圖描述清楚。

Fig 10 在畫什麼。 這是一張 K 空間的二維截面圖，畫出 K_x 與 K_y 兩條座標軸。圖上有兩個同心圓：外圈半徑標為 K_D（Debye 截止波向量），內圈半徑標為 K_T，內圈是塗滿灰色的實心圓，外圈與內圈之間是白的。右上角寫著兩個定義式：

ℏω_D = ℏvK_D = k_Bθ

ℏvK_T = k_BT

灰色＝被激發的模，白色＝完全沒被激發的模。二維畫面代表的其實是三維的兩個同心球。

論證本身。 Kittel 的原文推理是：只有 ℏω < k_BT 的晶格模會被明顯激發（「Only those lattice modes having ℏω < k_BT will be excited to any appreciable extent at a low temperature T.」）；而這些模的激發是近似古典的，每個模帶有接近 k_BT 的能量（「The excitation of these modes will be approximately classical, each with an energy close to k_BT, according to Fig. 1.」）。

接著是幾何：K 空間中被激發的模所佔的體積比例，量級上是 (ω_T/ω_D)³ 或 (K_T/K_D)³。由於 ℏvK_T = k_BT 且 ℏvK_D = k_Bθ，兩者相除給出 K_T/K_D = T/θ，所以體積比是 (T/θ)³。全部的模有 3N 個，於是被激發的模數是 3N(T/θ)³ 這個量級，每個帶 k_BT 的能量。

Fig 10 的圖說把結果寫死了，逐字是：「Of the 3N possible modes, the fraction excited is (K_T/K_D)³ = (T/θ)³, because this is the ratio of the volume of the inner sphere to the outer sphere. The energy is U ≈ k_BT · 3N(T/θ)³, and the heat capacity is C_V = ∂U/∂T = 12Nk_B(T/θ)³.」

12 是怎麼來的。 U ≈ 3Nk_BT⁴/θ³，對 T 微分得 12Nk_BT³/θ³ = 12Nk_B(T/θ)³。冪次 T³ 完全正確；係數 12 則不然——下一段專門處理這件事。

（以下為我的判讀）這個論證的價值在於它把 T³ 的來源攤成三個可以獨立檢驗的因子：（1）維度 3 來自 K 空間是三維的；（2）比值 T/θ 來自 ℏvK_T = k_BT 與 ℏvK_D = k_Bθ 這兩條線性色散關係——兩式相除即得 K_T/K_D = T/θ，把能量尺度換算成波向量尺度；（3）每模 k_BT 來自能量均分。任何一個因子改變，冪次就跟著變——這比記住 234 有用得多。

12 對 234：這是數量級論證，係數絕對不可引用我的判讀

這是本段最常見、也最致命的誤用，必須講到讀者不可能記錯。

事實。 精確的 Debye T³ 定律 (32) 給 C_V ≅ 234 Nk_B(T/θ)³；Fig 10 的手算論證給 C_V = 12 Nk_B(T/θ)³。兩者相差 234/12 = 19.5 倍。19.5 不是「小小的數值差異」，它是將近 20 倍——如果你拿 12 去擬合實驗，反推出來的 Debye 溫度會錯掉 19.5^{1/3} = 2.69 倍。矽的 645 K 會被算成 240 K。這不是誤差，是災難。

19.5 倍是從哪裡跑出來的。 兩個來源，而且都是結構性的：

（一）尖銳截斷 vs 平滑分布。 手算假設 K < K_T 的模「全滿」（每模 k_BT）、K > K_T 的模「全空」。真實的 Planck 分布 ⟨n⟩ℏω = ℏω/(e^{ℏω/k_BT}−1) 是一條平滑曲線：ℏω 略小於 k_BT 的模拿到的不到 k_BT，ℏω 略大於 k_BT 的模也不是完全為零。更重要的是，熱容量取的是導數，而導數對「分布在邊界附近怎麼變化」極度敏感——真正對 C_V 有貢獻的恰恰是那些 ℏω ~ k_BT、正處於過渡帶的模。手算論證把整個過渡帶抹平成一道階梯，等於把 C_V 最主要的貢獻區處理得最粗糙。

（二）K_T 的定義本身就是任意的。 「熱波向量」K_T 由 ℏvK_T = k_BT 定義——但為什麼是 1 倍的 k_BT？若改定義成 ℏvK_T = a·k_BT（a 是任意 O(1) 的數），整個論證一字不改地走下去，結果會乘上 a³。要讓答案從 12 變成 234，只要取 a = (234/12)^{1/3} = 2.69 就行了。換句話說，這個論證的係數可以透過一個沒有任何原理可以固定的選擇任意調整。這正是「數量級論證」的定義：它能定出冪次與參數組合，定不出前面的純數字。

怎麼記住不要犯錯。 Kittel 自己在 p.114 的措辭已經給了信號——他說這是「a simple argument」，用的是 ~ 符號（原文寫 U 是 ~3Nk_BT(T/θ)³、C_V 是 ~12Nk_B(T/θ)³，那個 ~ 讀作「量級為」），而在 Fig 10 圖說裡能量寫成 U ≈ k_BT·3N(T/θ)³ 也是約等於。相對地，式 (31)(32) 用的是 ≅ 並且附上完整積分。符號本身就是分級：有 ~ 或 ≈ 的地方係數不可引用，有完整推導的地方才可以。不過這條規則有一個例外必須記住：Fig 10 圖說裡的 C_V = ∂U/∂T = 12Nk_B(T/θ)³ 用的是等號，而那正是本條目要你拒絕引用的那個係數。它仍然只是量級結果，因為它不過是把同一句圖說裡用 ≈ 寫的 U ≈ k_BT·3N(T/θ)³ 對 T 微分而來——微分不會把精度變高。所以真正的判準不是「看到什麼符號」，而是「這條式子的上游有沒有完整的推導」。

更廣的教訓。 這種「相空間比例 × 每單元能量」的手算在固態物理裡到處都是——第 6 章估算電子比熱時（只有費米面附近 k_BT/ε_F 比例的電子有貢獻，故 C_el ~ Nk_B(T/T_F)）用的是一模一樣的招式，而那裡的精確係數是 π²/2 而不是手算給的 1。學會用它估冪次、同時拒絕相信它的係數，是這一段真正要訓練的能力。

T³ 律的實際適用範圍：書上兩處說法與 Fig 9 的實例

把書裡關於「T³ 到底什麼時候成立」的三處陳述並排，這件事才會清楚。

第一處，Fig 7 圖說（書本 p.113），逐字：「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 3 law is below 0.1θ. The asymptotic value at high values of T/θ is 24.943 J mol⁻¹ deg⁻¹.」

（附帶一提，24.943 J mol⁻¹deg⁻¹ 就是 3R = 3 × 8.314 = 24.94，也就是 Dulong–Petit 值，換算成卡是 5.96 cal/mol-K。這個數字待會講 Einstein 模型時會再用到。）

第二處，正文（書本 p.114），逐字：「For actual crystals the temperatures at which the T 3 approximation holds are quite low. It may be necessary to be below T = θ/50 to get reasonably pure T 3 behavior.」

第三處，Fig 9 圖說（書本 p.115），逐字：「Low temperature heat capacity of solid argon, plotted against T 3. In this temperature region the experimental results are in excellent agreement with the Debye T 3 law with θ = 92.0 K.」

Fig 9 的圖面描述（網頁不放圖）：橫軸是 T³，單位 K³，刻度 0、1.33、2.66、3.99、5.32、6.65、7.98；縱軸是熱容，單位 mJ mol⁻¹K⁻¹，刻度 0、4.44、8.89、13.33、17.78、22.23。圖中有數十個空心方塊（實驗點）從原點附近沿著一條穿過原點的直線排列到右上角，最高點約在 (7.98, 20)。實驗點與直線幾乎完全重合，看不出系統性偏離。圖面被格線分成 6×5 的方格。

把三處串起來：Fig 7 的「0.1θ」是 Debye 模型自己算出來的曲線在哪裡變成直線；正文的「θ/50」是真實晶體的實驗數據在哪裡才會躺成直線；而 Fig 9 的氬（θ = 92.0 K，測到 T = 2.0 K = θ/46）正好就是後者的示範。前一個條目已經用數值積分證明了：這兩個數字之間的差距不可能來自 Debye 積分的截斷誤差，只能來自 Debye 態密度假設本身。

Table 1（書本 p.116）：無文字層的整頁影像，以下為算繪後目視判讀

校註：書本 p.116 是整頁的 Table 1，這一頁在 PDF 裡沒有文字層（用抽取器只能得到頁碼與印刷記號共 40 個字元）。以下所有數值都是我用 PyMuPDF 把該頁以 3–4 倍放大算繪成影像後目視判讀的結果，不是從文字層讀出來的。引用這些數字時請註明此點；若研究上要用精確值，請回查原書或原始資料來源。

表格標題逐字：Table 1　Debye temperature and thermal conductivity。表格內橫跨中段的兩行說明逐字：Low temperature limit of θ, in Kelvin 與 Thermal conductivity at 300 K, in W cm⁻¹K⁻¹。表格本身排成一張週期表的形狀，每個元素格子裡上列（灰底）是 θ、下列（白底）是 300 K 的熱導率 κ；資料缺漏的元素該欄留空。

下面是與本段討論直接相關的選錄（目視判讀）：

鹼金屬（第 IA 族，由上而下）

第 IV 族（由上而下）

惰性氣體固體

其他常用值

注意 Ar 的 θ = 92 K 與 Fig 9 圖說的 92.0 K 一致——書內自洽。另外惰性氣體固體的熱導率欄全部空白（書上未解釋原因；我的判讀是它們在 300 K 早已不是固體）。

讀 Table 1：鹼金屬規則與 θ 的高低是什麼在決定我的判讀

Kittel 在 p.114 對 Table 1 只給了一句評註，逐字是：「Note, for example, in the alkali metals that the heavier atoms have the lowest θ's, because the velocity of sound decreases as the density increases.」

先確認事實。 由 Table 1（目視判讀）：Li 344 → Na 158 → K 91 → Rb 56 → Cs 38。從 Li 到 Cs，θ 掉了 344/38 = 9.1 倍。趨勢單調，沒有例外，確實是「越重的原子 θ 越低」。

Kittel 的理由與 (28) 怎麼對上。 由 (28)，θ ∝ v·(N/V)^{1/3}。Kittel 點名的是 v 這個因子：聲速隨密度上升而下降。彈性波的聲速量級是 v ~ √(彈性常數/質量密度)——沿鹼金屬族往下走，原子質量從 Li 的 6.9 增加到 Cs 的 132.9（約 19 倍），而鹼金屬的鍵結（單一 s 電子的金屬鍵）反而越往下越弱，分子分母兩頭都在把 v 壓低。

但 (28) 裡還有第二個因子，Kittel 沒提。 (N/V)^{1/3} 沿族往下也在下降——鹼金屬原子越往下越大，原胞體積越大，數密度越小。所以兩個因子同向作用，一起把 θ 壓下去。這解釋了為什麼 9.1 倍這麼大：它不是單一因子造成的。

推廣到整張表：θ 高的元素長什麼樣。 掃過 Table 1 可以歸納出兩個條件——

（一）原子輕。 鑽石 2230、Be 1440 是全表最高的兩個；兩者都是週期表最前段的元素。質量小 ⇒ 同樣的力常數下振動頻率高 ⇒ ω_D 高。

（二）鍵強、原子排得密。 鑽石的 sp³ 共價鍵是已知最強的鍵之一，加上短鍵長給出大的 (N/V)^{1/3}，兩個因子疊加，θ 衝到 2230 K。另一端，惰性氣體固體（Ne 75、Ar 92、Kr 72、Xe 64）全部落在 64–92 K 的窄帶：它們靠的是微弱的 van der Waals 吸引（第 3 章），鍵極軟，v 極低，所以 θ 全部很低且彼此差距不大。

一個容易被誤讀的地方。 有人會以為「θ 高 ⇒ 熱導率高」，因為兩者都跟「鍵硬、聲子快」有關。Table 1 直接否定了這個直覺：Cu 的 θ 只有 343（比 Si 的 645 低），κ 卻是 4.01 W cm⁻¹K⁻¹（Si 只有 1.48）——因為銅的熱導主要由電子而不是聲子攜帶。而鑽石 θ = 2230 是全表最高，κ = 1.29 卻低於銅、銀、金、鋁。Table 1 把 θ 與 κ 並排放在同一格，正是要讓讀者自己發現「這兩個量沒有簡單的比例關係」。熱導率的機制要到 p.121 以後（Thermal Conductivity 一節）才處理，那裡會看到 κ = ⅓Cvℓ，其中平均自由程 ℓ 由非諧交互作用與缺陷決定，跟 θ 沒有直接關係。

Einstein 模型：為什麼在 Debye 之後才講一個更早、更差的模型我的判讀

Einstein 的 1907 年模型在歷史上早於 Debye 的 1912 年模型五年，是第一個用量子化解釋固體熱容為什麼在低溫掉下來的理論。Kittel 卻把它放在 Debye 之後，而且放在「Einstein Model of the Density of States」這個標題底下——注意標題強調的是態密度，不是熱容。這個編排選擇本身就是教學設計，值得先說清楚。

理由一：這一節真正的主角是 D(ω)。 整個 PHONON HEAT CAPACITY 大節的骨架是 p.108 的式 (9)：U = Σ_p ∫dω D_p(ω)·ℏω/(e^{ℏω/τ}−1)。這條式子裡 Planck 因子是普適的、對所有固體都一樣；唯一因材料而異的是 D(ω)。Kittel 在 p.108 講完 (9) 之後說得很直接：「The central problem is to find D(ω), the number of modes per unit frequency range.」所以整節其實是在展示三種 D(ω)：Debye 的 D(ω) ∝ ω²（連續彈性體）、Einstein 的 D(ω) = Nδ(ω−ω₀)（單一頻率）、以及 (37) 的一般式。Einstein 模型在這裡的身分是「態密度的第二個範例」，不是「歷史上的舊理論」。

理由二：它是一個乾淨的對照組。 Debye 的 D(ω) 與 Einstein 的 D(ω) 是兩個極端：一個把全部的模攤在 0 到 ω_D 之間、低頻端連續地趨近零；另一個把全部的模堆在同一個頻率、低頻端完全是空的。兩者的高溫極限一模一樣（都是 Dulong–Petit 3Nk_B，因為高溫時每個模不論頻率都拿 k_BT，跟 D(ω) 長什麼樣無關），低溫行為卻定性地不同（T³ vs 指數）。這個「高溫相同、低溫不同」的對比，把「低溫熱容只由 D(ω) 的低頻端決定」這件事示範得再清楚不過。

理由三：它今天還在用。 Kittel 在 p.117 明說：「The Einstein model, however, is often used to approximate the optical phonon part of the phonon spectrum.」這不是客套話。真實材料的熱容常常被擬合成「一個 Debye 項（聲學支）＋ 一個或多個 Einstein 項（光學支）」，因為光學支的 ω(K) 在整個布里淵區內相對平坦（見第 4 章雙原子鏈的光學支），用單一頻率代表它並不是壞近似，而且光學支確實有能隙、確實在低溫指數凍結——Einstein 模型描述光學支時，那個「錯誤」的指數行為反而是正確的物理。

換句話說：Einstein 模型不是被淘汰的理論，是被降級使用的理論。它只是不能單獨拿來描述整個固體。

Fig 11（書本 p.117）：鑽石的實測與 Einstein 曲線

Kittel 用鑽石當 Einstein 模型的示範材料。圖說逐字是：「Comparison of experimental values of the heat capacity of diamond with values calculated on the earliest quantum (Einstein) model, using the characteristic temperature θ_E = ℏω/k_B = 1320 K. To convert to J/mol-deg, multiply by 4.186.」

圖面描述（網頁不放圖，以下為算繪後目視判讀）：縱軸標「C_p, in cal/mol-K」，刻度 0 到 6，每格 1；橫軸標「T/θ_E」，刻度 0 到 1.0，每格 0.1。圖中有一條虛線（Einstein 模型的計算曲線）與一組空心圓點（鑽石的實驗值）。虛線從左下角幾乎貼著橫軸出發（T/θ_E < 0.1 時實質為零），在 T/θ_E ≈ 0.2–0.4 之間急速上升，之後彎折趨緩，到 T/θ_E = 1.0 時逼近約 5.4 cal/mol-K 並趨於水平。實驗點從 T/θ_E ≈ 0.15（約 200 K）開始出現，在 0.15–0.4 這段與虛線幾乎完全重合，在 0.5 以上略高於虛線（最右邊三點在 0.65、0.8、0.95 附近，都落在虛線稍上方）。

為什麼選鑽石（以下為我的判讀，書上未說明選材理由）。 這是 Einstein 1907 年原論文就用的材料，理由很實際：鑽石的 θ_E = 1320 K 極高，所以室溫（300 K）就已經是 T/θ_E = 0.23 的「低溫」，熱容明顯低於 Dulong–Petit 值——量子效應在一般實驗室溫度下就看得見。換成鉛（Table 1 的 θ = 105 K），室溫是 T/θ ≈ 2.9，熱容早已飽和到 3Nk_B，什麼量子效應都看不到。鑽石是當時最容易觀察到「Dulong–Petit 定律失效」的材料之一（θ_E = 1320 K，室溫已是 T/θ_E = 0.23；Table 1 給的 Debye θ 更高達 2230 K，室溫是 T/θ = 0.13），這也是為什麼固體熱容的量子之謎最初是從鑽石身上被注意到的。它並不是唯一的——θ 高的材料在室溫都明顯低於 3R，例如本檔後面算出矽（θ = 645 K）在 300 K 也只有古典值的 0.80。

圖上看得到什麼、看不到什麼。 看得到的是：Einstein 曲線在 0.15 < T/θ_E < 0.4 這段與實驗吻合得相當好——這正是 Einstein 1907 年的重大成功，第一次用量子化解釋了熱容為什麼會掉。看不到的是這個模型在更低溫的失敗：圖的橫軸最低只到 0.15 左右就沒有資料點了，而 Einstein 與 Debye 的分歧要到 T/θ < 0.05 才變得災難性。Kittel 把這件事留給正文的文字陳述（「At low temperatures (34) decreases as exp(−ℏω/τ), whereas the experimental form … is known to be T 3」），圖上並沒有呈現。讀圖時不要以為「圖看起來吻合就代表模型對」——這張圖的溫度範圍剛好避開了模型最糟的區域。

順帶檢查一致性。 3R = 24.94 J mol⁻¹K⁻¹（Fig 7 圖說給的 24.943），除以 4.186 得 5.96 cal/mol-K，這就是圖上虛線最終要趨近的水平漸近線。曲線在 T/θ_E = 1.0 只到約 5.4（算繪量測 5.43），即 5.43/5.96 = 0.911；Einstein 函數 f(1) = 0.921 對應的是 5.49，兩者差約 1%，落在圖面繪製精度之內（這條虛線整體就比精確值低 1–2%）。方向與量值都對得上，可視為書內自洽。

General Result for D(ω)：問題是什麼

Kittel 開宗明義（書本 p.117），原文逐字：「We want to find a general expression for D(ω), the number of states per unit frequency range, given the phonon dispersion relation ω(K).」

這句話裡「given the phonon dispersion relation ω(K)」是重點：色散關係是輸入，態密度是輸出，而且不對 ω(K) 的形式做任何假設。這與前面兩個模型形成對比——Debye 假設 ω = vK（球對稱、線性），Einstein 假設 ω = ω₀（與 K 完全無關）。現在要處理的是任意形狀的 ω(K)，包括第 4 章那些在區邊界趨平、在不同方向上不同、分成聲學支與光學支的真實色散。

出發點是模的計數。 由 p.111 式 (18)，週期邊界條件下 K 空間中每單位體積有 (L/2π)³ = V/8π³ 個允許的 K 值（每個極化、每個分支各自獨立計數）。所以「頻率落在 ω 與 ω+dω 之間的模數」就等於「這個 K 密度」乘上「滿足此條件的 K 空間區域體積」。那個區域是兩個等頻面夾出來的殼層（shell）。這就是式 (35)。

Kittel 隨即點出困難所在，原文逐字：「The real problem is to evaluate the volume of this shell.」

為什麼這是個真問題。 在 Debye 模型裡殼層是兩個同心球面之間的薄殼，體積是 4πK²dK，小學幾何。但一般情形下等頻面是任意形狀的曲面（想像一顆凹凸不平的馬鈴薯，見 Fig 12），而且殼層的厚度沿著曲面是變化的——在色散陡的地方（頻率隨 K 變化快）殼層薄，在色散平緩的地方殼層厚。所以不能簡單地用「面積 × 厚度」，必須逐點積分。式 (36)、(37) 處理的就是這件事。

先預告結論的形狀，這樣讀推導時會知道每一步要往哪裡去：最終答案 (37) 是 D(ω) = (V/8π³)∫dS_ω/v_g——等頻面上的面積分，被積函數是群速度的倒數。「厚度」這個幾何量最後化成了「群速度倒數」這個動力學量，這個轉換正是整個推導的核心。

Van Hove 奇異點：v_g = 0 的地方會發生什麼

式 (37) 的推導在 p.118 結束，但這一小節的結論還有一段在 p.119。Kittel 在該處寫下（逐字）：

「The integral is taken over the area of the surface ω constant, in K space. The result refers to a single branch of the dispersion relation. We can use this result also in electron band theory. 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（書本 p.119）的圖面描述（網頁不放圖）：兩張並排的小圖，共用「D(ω) 對 ω」的座標。(a) 標題為 Debye solid：一條從原點出發、光滑向上彎的拋物線（ω² 形狀），在某處戛然而止（Debye 截止）——整條曲線完全平滑，沒有任何結構。(b) 標題為 an actual crystal structure：曲線在低頻端同樣以 ω² 起步，但往上走之後出現數個明顯的尖峰、轉折與台階，整體形狀不規則，最後在最高頻處落回零。圖說逐字：「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 ω 2 for small ω, but discontinuities develop at singular points.」

這張圖要說的三件事：

（一）低頻端 Debye 是對的。 真實晶體的 D(ω) 在小 ω 確實以 ω² 起步（「starts as ω² for small ω」）。這是前面 T³ 律正當性的圖像版——低溫只看得到低頻端，而低頻端兩張圖是一樣的。

（二）中高頻 Debye 完全錯了。 真實的 D(ω) 有豐富的結構，Debye 的光滑拋物線一個都沒抓到。這解釋了為什麼 Fig 7 的 Debye 曲線與真實材料的熱容在中間溫區會偏離，也解釋了為什麼從不同溫區擬合出來的「等效 θ」會不一樣。

（三）結構的來源是 v_g = 0。 由 (37)，被積函數是 1/v_g；群速度為零的臨界點讓被積函數發散，D(ω) 在對應的頻率上出現不光滑的特徵。Kittel 稱這些特徵為 Van Hove 奇異點。

與一維的對照（書本 p.111）。 Kittel 早在講一維態密度時就說過同一件事，逐字：「There is a singularity in D₁(ω) whenever the dispersion relation is horizontal; that is, whenever the group velocity is zero.」一維的 D₁(ω) = (L/π)/(dω/dK) 在 dω/dK → 0 時直接發散成無窮大。三維的情形要溫和一些——見下一條目。

為什麼 v_g = 0 一定會發生，以及三維為什麼只是「不連續」而不發散我的判讀

Kittel 指出 v_g = 0 產生奇異點，但沒有說明兩件事：這種點為什麼一定存在，以及為什麼 Fig 14(b) 畫的是「discontinuities」而不是發散到無窮。這兩點是我補的說明。

第一件事：v_g = 0 不是例外，是必然。

ω(K) 是倒晶格的週期函數——把 K 平移一個倒晶格向量 G，頻率不變（第 2 章的結果）。一個定義在有界週期域上的連續函數必定有極大值與極小值，而在光滑函數的極值處梯度為零。所以 ∇_Kω = 0 的點必然存在，而且至少有一個極小（K = 0 的聲學支底部）與一個極大。

更具體地，這些點通常落在布里淵區邊界與高對稱點上。第 4 章的單原子鏈給了最乾淨的例子：ω² = (4C/M)sin²(½Ka)，在 K = ±π/a 處 dω/dK = 0——該處的解是駐波，相鄰原子反相振動，波包不傳播。三維晶體在區邊界的每一個面、以及 Γ、X、L、W 等高對稱點附近都會出現類似的情形。

再加上鞍點（saddle point）：在三維中，若某個方向是極大、另一個方向是極小，梯度同樣為零。這種點不需要任何對稱性就會出現——只要函數在區內同時有極大與極小，拓樸上就強迫鞍點存在（這是 Morse 理論的內容）。所以每一條聲子分支都必定帶著一整組 Van Hove 奇異點。

第二件事：三維的奇異性長什麼樣。

（以下的標度分析是我依 (37) 做的局部展開；Kittel 在 p.119 只寫「discontinuities develop at singular points」，沒有給形式。）

在臨界點 K_c 附近把色散展開成二次型，最簡單的極小值情形是 ω ≈ ω_c + α|q|²（q ≡ K − K_c）。於是：

• 群速度 v_g = |∇_Kω| = 2α|q|，線性趨近零；
• 等頻面是半徑 |q| = √((ω−ω_c)/α) 的小球，面積 ∮dS = 4π|q|² = 4π(ω−ω_c)/α，趨近零。

代入 (37)：貢獻 ∝ 面積/v_g ∝ |q|²/|q| = |q| ∝ √(ω − ω_c)。

結果是有限的。面積趨零的速度（|q|²）比群速度趨零的速度（|q|）更快，兩者相除還剩一個 |q|，所以 D(ω) 本身不發散，只是以 √(ω−ω_c) 的方式從零長出來——函數連續，但導數在 ω_c 處無窮大。這就是 Fig 14(b) 裡那些尖銳的轉折。鞍點的情形類似，給出 D(ω) 中的 √ 型扭折或台階。Kittel 用「discontinuities」一詞指的正是這類斜率的不連續，不是函數值跳躍。

維度是關鍵。 同樣的分析在不同維度給出定性不同的結果：

這解釋了 p.111 與 p.119 兩處說法為什麼不同：一維（p.111）真的發散，三維（Fig 14）只是不連續。Kittel 兩處用詞的差異不是筆誤。

為什麼這件事重要。 低維系統的 Van Hove 發散會被實驗直接看見——石墨烯、碳奈米管、二維材料的光學吸收譜與 Raman 譜上那些又高又尖的峰，就是二維／一維 Van Hove 奇異點。在三維塊材裡它只是溫和的轉折，容易被忽略；一旦降維，它就變成主導的特徵。

(37) 在電子能帶論中的角色：一句原文，加上我補的跨章連結我的判讀

p.119 有一句話容易被讀過去，但它其實是整個第 5 章與後半本書的接點。原文逐字：「We can use this result also in electron band theory.」

為什麼可以直接搬用。 回頭看 (35)–(37) 的推導，從頭到尾只用了兩樣東西：

（一）週期邊界條件下 K 空間的態密度是均勻的 V/8π³（p.111 式 (18)）。這件事與波是什麼完全無關——它只取決於「樣品是體積 V 的盒子」與「用週期邊界條件」。聲子、電子、光子、磁子都一樣。

（二）梯度垂直於等值面。這是純數學。

沒有任何一步用到「這是晶格振動」。所以只要把 ω(K) 換成任何一個 K 空間中的能量（或頻率）函數，(37) 就成立。對電子，把 ω(K) 換成 ε(k)，群速度換成電子的群速度 v = (1/ℏ)∇_kε，並且每個 k 態可以容納兩個自旋方向（因子 2），得到

D(ε) = (2V/8π³) ∫ dS_ε/|∇_kε| = (V/4π³) ∫ dS_ε/|∇_kε|

以下全部是我補的跨章連結，Kittel 在 p.117–119 只寫了上面那一句話：

費米面附近的 D(ε_F) 幾乎決定了金屬的一切低溫性質——電子比熱係數 γ ∝ D(ε_F)、Pauli 順磁磁化率 ∝ D(ε_F)、超導轉變溫度透過 BCS 的 λ = ND(ε_F)V 進入指數。而 D(ε_F) 就是費米面上 ∮dS/|∇_kε| 這個面積分，其中 |∇_kε| = ℏv_F 是費米速度。「費米面的面積除以費米速度」這個組合之所以到處出現，根源就在 (37)。

半導體的能帶邊。 導帶底是一個 ε ≈ ε_c + ℏ²k²/2m* 的極小值，正是上一條目分析過的臨界點，給出 D(ε) ∝ √(ε−ε_c)——半導體物理最基本的態密度形式，直接來自 (37) 加拋物線近似。

光吸收譜的結構。 電子的 Van Hove 奇異點出現在聯合態密度（joint density of states）上，決定介電函數 ε₂(ω) 的峰位。實驗上量光吸收、反推能帶臨界點，用的就是這套邏輯。

一句話總結：(37) 不是聲子專用的公式，它是「把 K 空間的幾何翻譯成能量空間的態密度」這個操作的通用形式。Kittel 在這裡只用了一句話帶過，但這條式子在後面幾章會反覆出現。

本段最容易犯的五個錯誤我的判讀

把前面各條目的警告集中在一處，方便複習與自我檢查。

（一）引用 Fig 10 的係數 12。Fig 10 的手算論證給 C_V = 12Nk_B(T/θ)³，精確結果是 234Nk_B(T/θ)³，差 19.5 倍。手算只能定冪次與參數組合，不能定係數（因為 K_T 的定義 ℏvK_T = k_BT 裡那個「1 倍」是任意的，改成 a 倍結果就乘 a³）。只要看到 ~ 或 ≈ 就不要引用係數。

（二）把「T < 0.1θ」當成實驗判準。Fig 7 的「0.1θ」是 Debye 模型自己的曲線何時變成 T³ 直線（數學問題，我算出誤差 2.8%）；正文的「θ/50」是真實晶體的數據何時躺成直線（物理問題，取決於真實色散何時開始偏離線性）。兩者差五倍，而且差距不是來自積分截斷。做實驗時用後者。

（三）把 Einstein 模型的低溫失敗歸因於「近似太粗」。它不是精度問題，是譜的拓樸問題：δ 函數態密度人為製造了一個能隙，有能隙 ⇒ 指數凍結；真實聲學支無能隙 ⇒ 冪次律。再怎麼調 ω₀ 都補不回來（我算過：T = θ/50 時兩者差十五個數量級）。反過來，正因為光學支確實有能隙，Einstein 模型用在光學支上是恰當的——Kittel 逐字：「The Einstein model, however, is often used to approximate the optical phonon part of the phonon spectrum.」

（四）忘記 (37) 與 (22) 都是「單一分支」。Kittel 在 p.119 明說「The result refers to a single branch of the dispersion relation.」，p.112 的 (22) 也註明是 for each polarization type。Debye 模型可以簡單乘 3（因為假設三支聲速相同），真實晶體不行——三支聲學支的等頻面形狀與群速度都不同，必須分別積分。原胞有 p 個原子時還有 3p−3 條光學支要算。

（五）把 Table 1 的 θ 當成唯一的 Debye 溫度。表格標頭逐字是「Low temperature limit of θ, in Kelvin」——這是低溫極限值，正是配 (32) 用的那一個。從高溫熱容擬合、從彈性常數計算、從熔點經驗式估計得到的 θ 都不一樣（因為真實 D(ω) 不是 ω²，不同溫區的熱容對譜的不同部分加權）。論文裡引用 θ 一定要說明它是怎麼定出來的。

（附加）本頁資料的引用紀律。Table 1（書本 p.116）在 PDF 中是無文字層的整頁影像，本條目所有 θ 與熱導率數值都是算繪後目視判讀的結果。要在研究中使用請回查原書。

延伸：這一段對元件熱容（C_th）萃取的意義延伸知識

本條目是延伸內容（tier: ext）：元件熱容萃取不在 Kittel 第 5 章 p.114–119 的討論範圍內，以下是我為了連結讀者的研究主題（FinFET 熱容萃取）所補，請與書本內容分開看待。書中的數值（θ 值）來自 Table 1 的目視判讀。

第一件事：室溫的矽既不在 T³ 區、也還沒到 Dulong–Petit。由 Table 1，Si 的 θ = 645 K，所以室溫 300 K 對應 T/θ = 0.465。我用式 (30) 數值積分得到此處 C_V/3Nk_B = 0.80——也就是說室溫矽的晶格熱容只有古典飽和值的八成，而且還在隨溫度上升。這有兩個直接後果：

（a）不能用 T³ 律。T³ 只在 T ≲ θ/50 = 13 K 以下乾淨（而且 (32) 在 T = 0.465θ 會給出 234 × 0.465³ = 23.5 Nk_B，是實際值 2.4 Nk_B 的十倍，完全荒謬）。

（b）也不能用常數 3Nk_B。若在瞬態熱模擬中把 C_th 設成常數，在自熱使溫度從 300 K 升到 400 K 的過程中會系統性低估熱容（400 K 時 T/θ = 0.62，C_V/3Nk_B 上升到約 0.88）。溫度相依的 C_th(T) 必須保留。

第二件事：不同材料在同一個溫度處在完全不同的區域。這對 FinFET 的多層材料堆疊特別重要——同一個元件裡的不同材料，在同一個溫度下可能一個早已飽和、另一個還深在量子區。用 Table 1 的 θ 值（目視判讀）與我的 Debye 積分計算：

鑽石在室溫只有 17% 的古典熱容——這是它常被誤判的原因之一。

第三件事：(35) 的連續化假設在奈米尺度會鬆動。式 (35) 把對 K 的求和換成積分，前提是 K 點密到可以當連續分布。FinFET 的鰭寬已經到個位數奈米，在受限方向上允許的 K 值只剩下十幾個，聲子譜在該方向離散化、態密度出現子帶結構——這就是聲子侷限。此時 (37) 的三維面積分要換成「受限方向求和、自由方向積分」的混合形式，D(ω) 的低頻行為也不再是 ω²（二維薄膜是 ω¹、一維線是 ω⁰）。低頻 D(ω) 的冪次改變會直接改變低溫熱容的冪次——前面「U ∝ T^{d+1}」那條規則在這裡就用得上。

第四件事：θ 與熱導率無關，不要混用。Table 1 把 θ 與 300 K 熱導率放在同一格，很容易讓人以為兩者連動。它們不是：Cu（θ = 343，κ = 4.01）與 Si（θ = 645，κ = 1.48）的排序完全相反，因為銅的熱導由電子主導。熱容（本段主題）由 D(ω) 與 Planck 分布決定；熱導率還要乘上平均自由程，由非諧散射與邊界散射決定（書本 p.121 之後）。萃取 C_th 與萃取 R_th 是兩個獨立的問題，用的是兩套完全不同的物理。