01聲子熱容：普朗克分布、簡正模計數、一維態密度

這一段在整章裡的位置：一條三步走的計算流水線我的判讀

第 4 章解完了一件事：給定晶格與力常數，晶體能振動的頻率有哪些，也就是色散關係 ω(K)。那是純古典的力學結果，跟溫度無關。但實驗室裡量得到的是熱容——把樣品加熱 1 K 要餵進多少能量。要從 ω(K) 走到 C_V，中間缺的東西很具體：每一個模式在溫度 T 下平均裝了多少能量，以及總共有多少個模式。

Kittel 在 p.107–110 這四頁就是在補這兩塊。整條流水線可以寫成三步：

第一步（占據數）：一個頻率為 ω 的簡諧振子，在溫度 T 的熱平衡下平均被激發到第幾階？答案是 Planck 分布 ⟨n⟩ = 1/(exp(ℏω/τ) − 1)，這是 p.107–108「Planck Distribution」小節的全部內容。

第二步（模態計數）：晶體裡到底有幾個模式、它們的 K 值怎麼排？答案是每個可動原子（每種偏振）分到一個模式，而 K 在 K 空間裡等間距排列——固定端邊界條件間距 π/L，週期性邊界條件間距 2π/L。這是「Normal Mode Enumeration」與「Density of States in One Dimension」兩小節。

第三步（加總）：把第一步乘上 ℏω 再對第二步的所有模式加總，得到總能量 U；對 T 微分得到 C_lat。這一步的形式結果就是式 (8)–(10)，但真正要動手算，必須先知道態密度 D(ω)——這正是 Kittel 在 p.108 用一句話點名的「中心問題」，也是後面 Debye 模型與 Einstein 模型要處理的東西。

把這三步分開看，有一個非常重要的觀念收穫：第一步跟材料完全無關（任何頻率為 ω 的簡諧振子都服從同一條 Planck 分布），第二步跟溫度完全無關（模式有幾個、K 怎麼排，是幾何與邊界條件決定的，跟熱不熱沒關係）。材料的個性全部集中在 ω(K) 這條曲線上，也就是集中在 D(ω) 裡面。這個「分工」是整個固態熱學的骨架，後面談電子比熱、談熱導率，都是同一個骨架換一批零件。

為什麼談 C_V 而不是 C_p：Kittel 開宗明義的一個選擇

這一節的第一句話就在處理一個看似瑣碎、其實常被忽略的分歧。原文寫得很直白：「By heat capacity we shall usually mean the heat capacity at constant volume, which is more fundamental than the heat capacity at constant pressure, which is what the experiments determine.」——理論上比較根本的是定容熱容 C_V，但實驗量到的是定壓熱容 C_p。

接著定義 C_V ≡ (∂U/∂T)_V，其中 U 是能量、T 是溫度。然後 Kittel 指出：聲子對晶體熱容的貢獻叫做晶格熱容（lattice heat capacity），記作 C_lat。

這個「理論給 C_V、實驗給 C_p」的落差，正文只用第一句話帶過（就是上面那句逐字引文），量化的大小則塞進 p.107 的腳註 1 裡處理掉。腳註給了熱力學關係 C_p − C_V = 9α²BVT（α 是線膨脹係數的溫度係數、V 是體積、B 是體積模數），並且說這個相對差（原文 the fractional difference）在固體中通常很小、常常可以忽略；而且當 T → 0 時 C_p → C_V（前提是 α 與 B 為常數）。

（以下為本站補充的解讀，非原文論證。）注意這個腳註其實已經預告了本章後半段的主題：α ≠ 0（熱膨脹）本身就是非諧效應。也就是說「C_p 與 C_V 有差」和「晶體會熱膨脹」和「聲子會互相散射而產生有限熱導率」是同一件事的三個面向——都源自位能中高於二次的項。第 4 章的諧近似世界裡 α = 0，所以 C_p = C_V，晶體不膨脹，熱導率無窮大。本章從 p.119 的 ANHARMONIC CRYSTAL INTERACTIONS 起會回頭收拾這件事。

Kittel 專屬記號 τ ≡ k_BT：不是筆誤，是一個刻意的簡化我的判讀

p.107 第二段出現了一個第一次讀 Kittel 很容易卡住的東西：「The total energy of the phonons at a temperature τ(≡ k_BT) in a crystal...」。τ 在這本書裡不是弛豫時間、不是應力，而是以能量為單位的溫度，定義就是 τ ≡ k_BT。式 (3) 又把這個定義再寫一次，可見 Kittel 很清楚讀者會被絆到。

這個記號來自 Kittel 自己的另一本書《Thermal Physics》（與 Kroemer 合著），在那裡 τ 被稱為 fundamental temperature。用它的理由很實際：波茲曼因子寫成 exp(−E/τ) 比 exp(−E/k_BT) 乾淨，而且 ℏω/τ 一眼就看得出是「一個無因次的比值：聲子能量比熱能」。這個無因次組合在本章會反覆出現，到式 (10) 乾脆被正式命名為 x ≡ ℏω/τ = ℏω/k_BT。

讀這一段時要建立的直覺是：τ 是熱浴一次能遞出來的能量份額，ℏω 是這個模式一次能收下的能量份額。 兩者的比值 ℏω/τ 決定一切：

• ℏω ≪ τ（比值遠小於 1）：熱浴給得起，模式被大量激發，行為古典。
• ℏω ≫ τ（比值遠大於 1）：熱浴一次湊不出一顆聲子的能量，模式被「凍結」，對熱容幾乎沒有貢獻。

這個「湊不出一份就買不到」的量子化門檻，就是低溫熱容偏離 Dulong–Petit 的唯一原因。整章後面的 Debye 溫度 θ、Einstein 溫度 θ_E，本質上都只是把「哪些模式的 ℏω 大於 τ」這個界線換一種方式說出來。

另外提醒一個閱讀上的陷阱：因為 τ 有能量的量綱，所以像 exp(−ℏω/τ) 這種式子裡不會再出現 k_B。若在演算過程中發現 k_B 憑空冒出來或不見了，先檢查是不是把 τ 誤當成溫度（單位 K）在用。

Planck Distribution 小節的設定：一整組相同的簡諧振子

p.107 最下方的「Planck Distribution」小節，開場只有一句設定：「Consider a set of identical harmonic oscillators in thermal equilibrium.」——考慮一組相同的簡諧振子處於熱平衡。

接下來 Kittel 不是從配分函數出發，而是從一個更原始的東西出發：相鄰能階的布居數比值。原文寫的是「The ratio of the number of oscillators in their (n + 1)th quantum state of excitation to the number in the nth quantum state is」，接著就是式 (3)。

以下兩點是本站補充的解讀，原文沒有寫出來。

第一，為什麼比值只跟 n 無關？ 因為簡諧振子的能階是等距的：E_n = (n + ½)ℏω，所以 E_{n+1} − E_n = ℏω 對所有 n 都一樣。波茲曼因子只看能量差，於是 N_{n+1}/N_n = exp(−ΔE/τ) = exp(−ℏω/τ) 是一個與 n 無關的常數。等距能階 ⇒ 等比數列布居，這就是為什麼後面會出現幾何級數。如果能階不等距（例如非諧振子或原子的電子態），這一整套推導就垮了。

第二，零點能在這裡自動消失。 能階寫成 (n + ½)ℏω，那個 ½ℏω 對每一階都一樣，取能量差時被減掉了。所以式 (3) 不管你有沒有把零點能寫進能階，答案都一樣——這是本段一連串「零點能到底要不要算」的困惑的第一個埋伏點。

原文在寫下式 (3) 之後，下一頁（p.108）開頭才補上式 (3) 的依據：「by use of the Boltzmann factor.」——這是排版跨頁造成的斷句，讀 PDF 時容易漏掉。

Figure 1（書 p.106）的文字描述：曲線長什麼樣，虛線又是什麼

本小節（Planck 分布）唯一的一張圖在書 p.106（印在小節開始之前的整頁圖版上），圖說很短但埋了一個關鍵訊息。

圖的內容：橫軸標示為 x = k_BT/ℏω，範圍 0 到 4，有 1、2、3、4 的刻度；縱軸標示 ⟨n⟩，範圍 0 到 4，同樣有整數刻度與網格線。圖上畫了兩條線。

實線是 Planck 分布 ⟨n⟩ 本身。它從原點出發，一開始幾乎貼著橫軸（x < 0.5 的區域看起來就是零），然後在 x ≈ 1 附近開始翹起來，到 x ≈ 1.5 以後轉為一條近乎筆直的斜線，在 x = 4 時到達 ⟨n⟩ ≈ 3.5。

虛線是一條從原點出發、斜率為 1 的直線（通過 (1,1)、(2,2)、(3,3)、(4,4)）。實線始終位在虛線的下方，兩者的垂直間距在大 x 時看起來固定不變。

圖說原文：「Plot of Planck distribution function. At high temperatures the occupancy of a state is approximately linear in the temperature. The function ⟨n⟩ + 1/2, which is not plotted, approaches the dashed line as asymptote at high temperatures.」

翻成中文：普朗克分布函數的圖。高溫時一個態的占據數大致與溫度成線性。函數 ⟨n⟩ + ½（圖上沒有畫出來）在高溫時以那條虛線為漸近線。

這個圖說在暗示什麼（以下為本站解讀）：它明白告訴你，高溫下真正筆直的不是 ⟨n⟩，而是 ⟨n⟩ + ½。那個 ½ 不是隨手加的修正，它就是零點能的占據數對應物——把每個模式的能量寫成 (⟨n⟩ + ½)ℏω，高溫下才會乾淨地等於 τ。Kittel 在這張圖的圖說裡偷偷把零點能放進來了，但正文的式 (1) 卻沒有它。這個不一致值得單獨一節來談。

本章最容易打架的地方：式 (1) 沒有零點能，習題 5 卻說必須保留我的判讀

這是讀第 5 章時最容易產生內在矛盾的一個點，值得正面攤開。

事實一（書 p.99，第 4 章）：Kittel 在 QUANTIZATION OF ELASTIC WAVES 一節寫下式 (4.27) ε = (n + ½)ℏω，並且明講「The term ½ℏω is the zero point energy of the mode.」——一個模式的能量階梯是從 ½ℏω 起跳的，不是從 0 起跳。

事實二（書 p.107，本段）：式 (1) 寫的是 U_lat = Σ_K Σ_p ⟨n_{K,p}⟩ℏω_{K,p}。這裡面沒有 ½。照字面讀，絕對零度時 ⟨n⟩ = 0，所以 U_lat = 0。

事實三（書 p.106 圖說）：Fig 1 的圖說特別提到「The function ⟨n⟩ + 1/2」——那個 ½ 又冒出來了。

事實四（書 p.129，習題 5(a)）：要證明一個頻率為 ω 的聲子模式的自由能是 k_BT ln[2 sinh(ℏω/2k_BT)]，原文特別加了一句「It is necessary to retain the zero-point energy to obtain this result.」——必須保留零點能才能得到這個結果。

這四件事看起來互相打架，但其實不矛盾，關鍵在於問的是什麼物理量。

要算熱容 → 零點能可以丟。 零點能 Σ_K Σ_p ½ℏω_{K,p} 是一個與溫度無關的常數（只要 ω 不隨 T 變）。C_V = (∂U/∂T)_V 對它微分得零。所以式 (1) 少寫那一項，對本章從頭到尾要算的 C_lat 完全沒有影響。Kittel 省掉它是為了算熱容這個特定目的所做的合法簡化，不是疏忽。

要算自由能、熱膨脹、Grüneisen 常數 → 零點能絕對不能丟。 因為此時 ω 本身要隨體積變（δω/ω = −γΔ，習題 5(b) 的設定）。零點能 Σ ½ℏω_K 因此是體積的函數，對體積微分不為零，會產生一個真實的壓力（零點壓力）。習題 5(b) 要證明的 BΔ = γΣ ½ℏω coth(ℏω/2k_BT) 裡，那個 ½ℏω coth(…) 在 T → 0 時退化成 ½ℏω——絕對零度的晶體仍然被零點振動撐開，這就是零點能非丟不可的物理後果。

所以正確的讀法是：式 (1) 是「熱能」而不是「總能」。它算的是 U(T) − U(0) 這個差。而習題 5 問的是自由能的絕對值與它對體積的導數，差一個與體積有關的常數就會得到錯的答案。

一個實務上的自我檢查：每當你看到一個式子裡的聲子能量寫成 ⟨n⟩ℏω 而不是 (⟨n⟩ + ½)ℏω，就問自己「我接下來要對什麼微分？」對 T 微分（熱容、熵）——安全；對 V 或對晶格參數微分（壓力、體積模數、熱膨脹、聲子軟化）——不安全，必須把 ½ 補回去。這個判準在做第一原理聲子計算時同樣適用：計算聲子自由能時 ZPE（zero-point energy）項是否納入，直接決定了平衡晶格常數算不算得準。

校註：以全書搜尋確認，「zero point / zero-point」在第 5 章只出現在章首目錄頁（p.105，作為習題 3 的標題）、習題 3（p.128）與習題 5（p.129）三處，正文 p.107–126 完全沒有提到零點能。上面關於「為什麼可以省略」的解釋是本站補上的，不是原文的論證。

Normal Mode Enumeration：把式 (1) 與式 (2) 併起來

「Normal Mode Enumeration」小節的開頭只做一件很機械的事：把式 (2) 的 ⟨n⟩ 代進式 (1) 的 U_lat，得到式 (8)。原文的說法是「The energy of a collection of oscillators of frequencies ω_{K,p} in thermal equilibrium is found from (1) and (2)」。

（以下為本站補充。）這一步沒有新物理，但它標誌著推導的重心轉移：Planck 分布的部分結束了，接下來的問題全部是「模式怎麼數」。

這也是小節標題 Normal Mode Enumeration（簡正模計數）的意思。注意 enumeration 這個字用得很準確——它不是「近似」也不是「估計」，而是字面上的點數：一個一個把模式數出來。接下來兩頁 Kittel 真的就是在 N = 10 與 N = 8 的小系統上一個一個點。

從對 K 求和換成對 ω 積分：D_p(ω) 的定義就在這一句話裡

式 (8) 之後 Kittel 寫：「It is usually convenient to replace the summation over K by an integral. Suppose that the crystal has D_p(ω)dω modes of a given polarization p in the frequency range ω to ω + dω. Then the energy is」——接著就是式 (9)。

（以下的量綱討論與適用條件為本站補充，原文只寫了那一句 Suppose。）D_p(ω) 的定義完全包在那個 Suppose 裡：D_p(ω)dω 是偏振為 p、頻率落在 ω 到 ω + dω 之間的模式個數。注意它是先定義 D(ω)dω（一個純數），再由此反推 D(ω) 的量綱是「每單位角頻率的模式數」，也就是時間的量綱（秒）。這一點在檢查任何態密度公式時都很好用：D(ω) 的量綱必須是 [時間]，若算出來不是，中間一定掉了東西。

換成積分之後，式 (8) 的 Σ_K 被 ∫dω D_p(ω) 取代，而原本依賴 K 的 ω_{K,p} 變成積分變數 ω 本身。剩下的 Σ_p 保留——偏振沒有被積掉，因為不同偏振有不同的色散關係，因而有不同的 D_p(ω)。這是一個很多人抄公式時會漏掉的細節。

這個代換在什麼條件下合法？Kittel 沒有明說，但條件很清楚：相鄰 K 值的間距要遠小於 ω(K) 變化的尺度。下一小節算出的 ΔK = π/L 對宏觀樣品（L ~ 1 cm）是 10² m⁻¹ 等級，而 K 的整個範圍是 10¹⁰ m⁻¹ 等級，兩者差了七到八個數量級，所以這個代換在宏觀晶體裡好到不能再好。反過來，在奈米結構裡（第 18 章）L 小到 nm，這個代換就會失效，模式的離散性會直接顯現在光譜上。

Kittel 把整章的目標壓縮成一句話

式 (10) 之後，Kittel 寫下本章最重要的一句路標：「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.」

中文：中心問題是找出 D(ω)，也就是每單位頻率範圍的模式數。這個函數叫做模態密度，或者更常見的說法，態密度。

以下三個細節是本站補充的解讀。

第一，名稱的雙軌。Kittel 明說這個函數的本名是 density of modes（模態密度），但大家更常叫它 density of states（態密度）。在聲子的脈絡下「態」就是「模式」，兩者同義。但要注意這個詞後來會被借去講電子（第 6 章的 D(ε)），而電子的「態」還要再乘自旋簡併度 2，計數規則不同——同一個詞在兩章裡的細節不一樣，是常見的混淆源。

第二，D(ω) 是「每單位頻率」而不是「每單位波向量」。前面數出來的 L/π 與 L/2π 是每單位 K 範圍的模式數，那是幾何量、跟材料無關；D(ω) 才是把色散關係吃進去之後的結果。式 (15) 做的正是這個轉換。

第三，這句話定義了本章剩下部分的結構。接下來 Kittel 會依序做：一維的 D(ω)（本段結尾）、三維的 D(ω)（p.111）、Debye 模型的 D(ω)（p.112）、Einstein 模型的 D(ω)（p.114）、以及一般化的 D(ω) 結果（p.117）。五種做法在回答同一個問題。

一維固定端問題的設定，與 Figure 2（書 p.109）的文字描述

「Density of States in One Dimension」小節從一個非常具體的邊界值問題開始：「Consider the boundary value problem for vibrations of a one-dimensional line (Fig. 2) of length L carrying N + 1 particles at separation a. We suppose that the particles s = 0 and s = N at the ends of the line are held fixed.」

注意數目：線上總共有 N + 1 個粒子（編號 s = 0, 1, 2, …, N），間距 a，所以線長 L = Na。兩端的 s = 0 與 s = N 被固定住，因此可以動的只有 N − 1 個。這個「N+1 個粒子、N−1 個可動」的計數是後面所有結論的根，一開始就要記清楚。

Figure 2 的內容（書 p.109 上方）：一條水平線，上面畫著 11 個實心圓點代表原子，最左端標 s = 0、最右端標 s = 10，左右兩側各有一個箭頭指向端點原子並標註「Fixed」。線的上方用雙箭頭標出全長 L，中間某處用一小段雙箭頭標出相鄰原子間距 a。圖中線段中段上方畫了一支向上的短箭頭，標註 u_s，箭頭指到一條短橫線，短橫線左側有一個小空心圓，標出原子位移後的位置；線上的實心黑點才是原子的平衡位置。（校註：箭頭底端落在兩個相鄰實心點的中間，不在任何一個原子正上方；空心圓代表的是位移後而非平衡位置，這一點與 p.110 Fig 4 的畫法一致。以上經算繪 p.109 原頁目視確認。）

Figure 2 的圖說原文：「Elastic line of N + 1 atoms, with N = 10, for boundary conditions that the end atoms s = 0 and s = 10 are fixed. The particle displacements in the normal modes for either longitudinal or transverse displacements are of the form u_s ∝ sin sKa. This form is automatically zero at the atom at the end s = 0, and we choose K to make the displacement zero at the end s = 10.」

這段圖說把整個邏輯講完了：取 sin 形式，左端自動滿足；右端則靠挑選 K 來滿足。而且特別點明 u_s 可以是縱向或橫向位移，兩者形式相同——這就是後面說「每個 K 有三個偏振」的伏筆。

正文接著說每個偏振為 p 的簡正振動模都具有駐波的形式，其中 u_s 是第 s 個粒子的位移，即式 (11)；並補一句 ω_{K,p} 由適當的色散關係與 K 相聯繫。這句「由適當的色散關係」是在把第 4 章的成果掛進來：模式的空間形狀由邊界條件決定，模式的頻率由晶格動力學決定，兩件事各自獨立。

式 (13) 與 K = Nπ/L 的退化解：Kittel 為什麼特地點出這一個

式 (12) 之後，原文做了兩件檢查。

檢查一（式 13）：「The solution for K = π/L has u_s ∝ sin(sπa/L) and vanishes for s = 0 and s = N as required.」——最長波長那個模式的位移形式是 u_s ∝ sin(sπa/L)，在 s = 0 與 s = N 兩端都為零，符合要求。把 L = Na 代進去會看得更清楚：u_s ∝ sin(sπ/N)，在 s = 0 與 s = N 處確實為零，而在 s = N/2（線的正中央）達到最大值。這就是基頻駐波：半個波長剛好塞滿整條線。

檢查二（退化解）：「The solution for K = Nπ/L = π/a = K_max has u_s ∝ sin sπ; this permits no motion of any atom, because sin sπ vanishes at each atom.」——K 取到布里元區邊界時，位移形式變成 sin(sπ)，而 s 是整數，所以每個原子的位移都是零，這個「解」不允許任何原子運動。

（以下的跨章對照為本站補充；所引第 4 章內容已回原書核對。）這第二個檢查看似技術性，其實非常重要，因為它連接到第 4 章的一個結論。在第 4 章 p.93 式 (4.12)，Kittel 已經指出在區邊界 K_max = ±π/a 處，行波解變成 u_s = u(−1)^s 的駐波，相鄰原子反相振動。同樣一個 K = π/a，在週期性邊界下（cos 形式）是一個合法的、相鄰原子反相的駐波；在固定端邊界下（sin 形式）卻是一個恆為零的偽解。 差別完全來自 sin 與 cos 的選擇，而 sin 的選擇是被左端固定條件逼出來的。這件事在 p.110 的 Fig 4 圖說裡會被再說一次（8π/8a 只對 cos 形式有意義）。

結論：「Thus there are N − 1 allowed independent values of K in (12). This number is equal to the number of particles allowed to move. Each allowed value of K is associated with a standing wave.」——式 (12) 裡有 N − 1 個允許的獨立 K 值，這個數目等於可以運動的粒子數；每個允許的 K 對應一個駐波。

自由度守恆是這裡真正的物理原理：N − 1 個可動粒子在一維上有 N − 1 個自由度，簡正模變換是一個保維度的線性變換（把 N−1 個位移座標換成 N−1 個簡正座標），所以模式數必然等於 N − 1，一個不多一個不少。Kittel 沒有用「自由度守恆」這四個字，但整段計數的骨架就是它。

Figure 3（書 p.109）的文字描述，以及那句必須逐字記住的圖說

Figure 3 的內容：這是一張 K 空間的圖，不是實空間的圖。畫面上是一條水平的直線段（K 軸），下方有一個標註「K」的箭頭指向右方表示 K 增加的方向。線上散布著 9 個實心圓點，等間距排列。最左端標刻度 0（此處沒有點）；第一個點下方標 π/10a，第二個點下方標 2π/10a，第三個點之後用「....」省略；線的最右端標著 10π/10a，但這個位置上沒有點。

（校註：9 個點與「10π/10a 處無點」這兩項是把書 p.109 原頁以 500 dpi 算繪後目視點算確認的，文字層無法提供這個資訊。9 = N − 1，與 N = 10 的情形完全吻合。）

Figure 3 的圖說原文（逐字）：「The boundary condition sin sKa = 0 for s = 10 can be satisfied by choosing K = π/10a, 2π/10a, . . ., 9π/10a, where 10a is the length L of the line. The present figure is in K space. The dots are not atoms but are the allowed values of K. Of the N + 1 particles on the line, only N − 1 are allowed to move, and their most general motion can be expressed in terms of the N − 1 allowed values of K. This quantization of K has nothing to do with quantum mechanics but follows classically from the boundary conditions that the end atoms be fixed.」

這段圖說有三個各自獨立的重點：

(a) 點不是原子。「The dots are not atoms but are the allowed values of K.」Fig 2 與 Fig 3 長得很像（都是一條線上排著點），但一張在實空間、一張在 K 空間。這是初學者最容易看混的兩張圖，Kittel 顯然知道，所以在圖說裡明講。

(b) N − 1 個可動粒子的最一般運動，可以用 N − 1 個允許的 K 值展開。這是簡正模完備性的敘述：任何初始位移與速度組態，都能唯一分解成這 N − 1 個駐波的線性組合。這是「模式數必須等於自由度數」的另一種說法。

(c) 那句必須逐字記住的話：「This quantization of K has nothing to do with quantum mechanics but follows classically from the boundary conditions that the end atoms be fixed.」

正面打掉一個標準誤解：K 的量子化跟量子力學一點關係也沒有我的判讀

上一個區塊引的那句圖說值得單獨拿出來講，因為它打掉的是一個非常普遍、而且非常頑固的誤解。

誤解長這樣：「固態物理裡到處都是 2π/L 這種離散的 K 值，這一定是量子力學的效果，因為量子力學就是把連續的東西變成離散的。」

為什麼是錯的：Kittel 的原話是「This quantization of K has nothing to do with quantum mechanics but follows classically from the boundary conditions that the end atoms be fixed.」——K 的量子化完全來自古典的邊界條件。

證據很直接：整個式 (11) 到 (14) 的推導裡，ℏ 一次都沒有出現。用到的只有三件事：(a) 解是駐波 sin sKa；(b) 兩端必須為零；(c) s 是整數。這三件事在牛頓力學裡就成立。一根兩端綁死的吉他弦，它的允許波長同樣是離散的（λ = 2L/n），沒有人會說吉他弦是量子的。

那量子力學在哪裡？ 在別的地方，而且要分清楚：

這張表其實就是本段的全部內容按「古典／量子」重新分類一次。Planck 分布那半段是量子的，模態計數那半段是古典的，兩者在式 (9) 會合。

為什麼這個區分在實務上重要：因為它告訴你哪些結果可以放心地從連續彈性理論（第 3 章）借過來，哪些不行。模態計數 L/2π、態密度的幾何因子、布里元區的大小——這些在古典彈性波理論裡就已經成立，做有限元素法的聲學模擬也會算出同樣的東西。只有當你要問「這個模式裡裝了多少能量」時，才非得動用 ℏ 不可。

順帶一提，K 的離散性在宏觀樣品裡實際上觀察不到（間距 10² m⁻¹ 對上區大小 10¹⁰ m⁻¹），它在推導中的角色純粹是一種計數的腳手架。真正有物理後果的離散性是 (n+½)ℏω 那個能量階梯。這也是為什麼固定端與週期性邊界條件可以隨便挑一個用——腳手架換一副，蓋出來的房子一樣。

固定端算法的兩個結論：模態密度 L/π，以及三個偏振

p.109 在點完模式之後，給出兩個結論。

結論一（K 空間的模態密度）：「For the one-dimensional line there is one mode for each interval ΔK = π/L, so that the number of modes per unit range of K is L/π for K ≤ π/a, and 0 for K > π/a.」

中文：一維線上每一個 ΔK = π/L 的區間有一個模式，因此每單位 K 範圍的模式數是 L/π，適用範圍是 K ≤ π/a；超出 π/a 則為 0。

這個「L/π」是本段最重要的數字之一，式 (15) 會直接用它。要注意兩件事：(a) 它是密度（每單位 K），不是總數；(b) 它有適用範圍，K > π/a 之外模態密度是零，因為那裡的波向量只是重複第一布里元區內的運動（第 4 章 p.93 的結論）。這個截斷是晶格離散性的直接後果，也是晶體與連續介質最根本的差別——連續介質的 K 可以無限大，晶格不行。

結論二（偏振）：「There are three polarizations p for each value of K: in one dimension two of these are transverse and one longitudinal. In three dimensions the polarizations are this simple only for wavevectors in certain special crystal directions.」

中文：每一個 K 值對應三個偏振 p；在一維情況下其中兩個是橫波、一個是縱波。在三維情況下，偏振只有在某些特定的晶體方向上才這麼單純。

以下為本站補充。第一，「一維」指的是波傳播方向是一維，但原子仍然可以往三個方向位移，所以偏振有三個。這是為什麼一條 N−1 個可動原子的鏈總共有 3(N−1) 個模式，而不是 N−1 個。第二，Kittel 特別加了一句免責聲明：三維時偏振只在特殊方向上才能乾淨地分成純縱與純橫，一般方向上的簡正模是混合的。這句話在第 4 章 Fig 8a（鍺沿 [111] 這個特殊方向，四支分別標成 LA/TA/LO/TO，縱橫分得很乾淨，書 p.96）已經先看到一個實例，本章後面的三維態密度討論會再回到它。

另一種數法：週期性邊界條件

p.109 最後一段開啟了第二套算法：「Another device for enumerating modes is equally valid. We consider the medium as unbounded, but require that the solutions be periodic over a large distance L, so that u(sa) = u(sa + L).」

中文：另一種計數模式的手法（device）同樣有效。我們把介質視為無界的，但要求解在一個大距離 L 上是週期的，也就是 u(sa) = u(sa + L)。

「device」這個用字很關鍵——Kittel 明白地把它定位為一種計算手法，而不是對真實樣品的物理描述。真實晶體的表面當然不是週期性接回去的。接著他加了一句保證：「The method of periodic boundary conditions (Figs. 4 and 5) does not change the physics of the problem in any essential respect for a large system.」——週期性邊界條件的方法，對一個大系統而言，在任何本質的層面上都不改變問題的物理。

為什麼要換一套？ 原文沒有明說理由，但從式 (14) 的形式可以看出來：週期性邊界允許的解是行波 u_s = u(0)exp[i(sKa − ω_K t)]，而不是駐波。行波形式在做群速度、熱傳導、聲子散射（本章後半的 THERMAL CONDUCTIVITY 與 Umklapp processes）時遠比駐波好用，因為「聲子往哪個方向跑」這件事在駐波裡根本問不出來。所以這不只是換個記號，而是換一組更適合後續計算的基底。

代價：K 現在有正有負（行波可以往左或往右跑），K 的間距也從 π/L 變成 2π/L。這兩個變化恰好互相補償——見後面的對照區塊。

Figure 4（書 p.110）：N = 8 的圓環，把八個模式一個一個點出來

這是本段最值得慢讀的一張圖，因為它把「為什麼每個可動原子恰好分到一個模式」用最小的例子完整驗證了一次。

Figure 4 的內容：一個大圓（代表環），圓周上均勻分布 8 個實心圓點代表粒子。圖上標了 s = 1（圓頂偏右）、s = 2（右上）、s = 7（左側）、s = 8（左上）幾個編號，其餘粒子未標號。圓周下方偏右處畫了一個小空心圓與一個小方框，用箭頭標註 u_s，表示某個粒子沿環的切線方向偏離平衡位置的位移；另有一段弧長用箭頭標註 a，表示相鄰粒子的間距。

Figure 4 的圖說（逐字）：「Consider N particles constrained to slide on a circular ring. The particles can oscillate if connected by elastic springs. In a normal mode the displacement u_s of atom s will be of the form sin sKa or cos sKa: these are independent modes. By the geometrical periodicity of the ring the boundary condition is that u_{N+s} = u_s for all s, so that NKa must be an integral multiple of 2π. For N = 8 the allowed independent values of K are 0, 2π/8a, 4π/8a, 6π/8a, and 8π/8a. The value K = 0 is meaningless for the sine form, because sin s0a = 0. The value 8π/8a has a meaning only for the cosine form, because sin (s8πa/8a) = sin sπ = 0. The three other values of K are allowed for both the sine and cosine modes, giving a total of eight allowed modes for the eight particles. Thus the periodic boundary condition leads to one allowed mode per particle, exactly as for the fixed-end boundary condition of Fig. 3. If we had taken the modes in the complex form exp(isKa), the periodic boundary condition would lead to the eight modes with K = 0, ±2π/Na, ±4π/Na, ±6π/Na, and 8π/Na, as in Eq. (14).」

把它逐一點算成一張表（這是圖說文字的重新整理，內容全部來自原文）：

八個粒子，八個模式。原文的結論句是：「Thus the periodic boundary condition leads to one allowed mode per particle, exactly as for the fixed-end boundary condition of Fig. 3.」

兩端各掉一個，中間各補一個——這個平衡不是巧合（以下為本站解讀）。 K = 0 與 K = 8π/8a = π/a 這兩個特殊點各只給一個模式，而中間三個 K 值各給兩個（sin 與 cos 獨立）。總數 1 + 2 + 2 + 2 + 1 = 8。若換成複數形式 exp(isKa)，記帳方式改變但總數不變：K = 0 給一個，±2π/Na、±4π/Na、±6π/Na 各給兩個（共六個），8π/Na 給一個（因為 +8π/Na 與 −8π/Na 是同一個解），合計 1 + 6 + 1 = 8。同一組物理模式的兩種等價記法。

注意圖說裡 K 用 Na 當分母（2π/Na 等），而式 (14) 用 L 當分母（2π/L 等），因為 L = Na，兩者相同。

Figure 5（書 p.110）：同一件事畫在 K 空間裡

Figure 5 的內容：一條水平的 K 軸，下方標註「K」與向右的箭頭。軸上有 9 個實心圓點，由左至右依序標註：−Nπ/L、−6π/L、−4π/L、−2π/L、0、2π/L、4π/L、6π/L、Nπ/L。點的間距均勻。

Figure 5 的圖說（逐字）：「Allowed values of wavevector K for periodic boundary conditions applied to a linear lattice of periodicity N = 8 atoms on a line of length L. The K = 0 solution is the uniform mode. The special points ±Nπ/L represent only a single solution because exp(iπs) is identical to exp(−iπs); thus there are eight allowed modes, with displacements of the sth atom proportional to 1, exp(±iπs/4), exp(±iπs/2), exp(±i3πs/4), exp(iπs).」

以下三個要點中，引號內為原文，其餘為本站解讀。

要點一：畫了 9 個點，卻只有 8 個模式。 圖說把原因說得很清楚——最左與最右那兩個點 ±Nπ/L 其實是同一個解，因為 exp(iπs) 與 exp(−iπs) 對整數 s 完全相同（兩者都等於 (−1)^s）。所以點數 9，模式數 8。這是 K 空間計數最常見的重複計算陷阱：布里元區的兩個邊界面在物理上是同一個點（差一個倒晶格向量），只能算一次。

要點二：位移的完整清單。圖說直接把八個模式的位移形式列出來：1、exp(±iπs/4)、exp(±iπs/2)、exp(±i3πs/4)、exp(iπs)。數一數：1 個 + 2 個 + 2 個 + 2 個 + 1 個 = 8。這與 Fig 4 圖說最後給的複數形式清單是同一組東西：K = 0 對應位移 1（所有原子同相，均勻模式）；K = ±2π/L 對應 exp(±iπs/4)（因為 Ka = 2πa/L = 2π/8 = π/4）；K = Nπ/L = π/a 對應 exp(iπs) = (−1)^s，相鄰原子反相——正是第 4 章 p.93 式 (4.12) 那個區邊界駐波。

要點三：K = 0 是均勻模式。圖說明講「The K = 0 solution is the uniform mode」。這個模式是整個晶體的剛體平移，頻率為零，不儲存彈性能。它在熱容計算中貢獻可忽略（一個模式對上 10²³ 個），但在形式上它是週期性邊界條件比固定端多出來的那一個模式。

群速度為零 ⇒ 態密度發散：一維的 Van Hove 奇異點我的判讀

校註：「群速度為零 ⇒ D₁(ω) 發散」這個結論本身是原文寫的，不是本站的推論。 Kittel 在式 (15) 之後的下一句（p.111 第一句）就明文寫出：「There is a singularity in D₁(ω) whenever the dispersion relation is horizontal; that is, whenever the group velocity is zero.」以下只是把這句話的機制與跨章後果展開。真正要到 p.119 才出現的，是 Van Hove 這個名字（以及 critical point 的說法），不是發散這件事本身。

式 (15) 把 D₁(ω) 寫成 (L/π) 除以 dω/dK。這個寫法立刻暴露一個結構性的後果：只要色散曲線在某處變平（dω/dK → 0），態密度就在該頻率發散。

這種情形一定會發生，不是特例。 第 4 章 p.94 式 (4.14) 給出單原子鏈的群速度 v_g = (Ca²/M)^{1/2} cos(½Ka)，並明說「This is zero at the edge of the zone where K = π/a」。第 4 章 Fig 6（書 p.95）畫的就是這件事：群速度從 K = 0 的最大值單調下降，在區邊界 K = π/a 處觸地為零。而區邊界的存在是晶格週期性的必然後果。所以：只要是晶格，色散關係在布里元區邊界就一定變平，態密度就一定在最高頻率處發散。

（第 4 章 p.92 也用另一條路說過同一件事：由式 (4.7) 可得 dω²/dK = (2Ca/M)sin Ka，在 K = ±π/a 處因 sin(±π) = 0 而為零。）

書上哪裡驗證了這個發散？ 本章習題 1（書 p.128）「Singularity in density of states」第 (a) 小題要求：由第 4 章導出的單原子線性晶格色散關係，證明態密度為

D(ω) = (2N/π) · 1/(ω_m² − ω²)^{1/2}

其中 ω_m 是最大頻率。這個式子在 ω → ω_m 時以 (ω_m − ω)^{−1/2} 發散——逆平方根發散，一維 Van Hove 奇異點的標準形式。習題 1(b) 接著處理三維的情形（光學支 ω(K) = ω₀ − AK² 在 K = 0 附近），答案是 D(ω) = (L/2π)³(2π/A^{3/2})(ω₀ − ω)^{1/2}，在 ω < ω₀ 有限、ω > ω₀ 為零——三維是階梯式不連續而非發散。同一個機制在不同維度給出不同強度的奇異性，這是很好的對照。

這個名字在書裡的出處：Kittel 要到本章 p.119「General Result for D(ω)」小節的末尾才點名。他在導出三維一般式 D(ω) = V/(2π)³ ∫dS_ω/v_g（式 (37)，書 p.118）之後指出：群速度為零的那些點特別值得注意，這些臨界點（critical points）會在 D(ω) 中產生奇異性（Van Hove singularities），並指向 Fig 14。所以本段的式 (15) 是那個結論的一維前身——同一個機制（分母上的群速度歸零），一維版本先在這裡露面，三維的完整處理留到 p.117–119。

為什麼這件事在整本書裡一直回來：因為分母放群速度這個結構，不只屬於聲子。書 p.119 明說這個結果也可以用在電子能帶理論上（「We can use this result also in electron band theory.」）。書中真正兌現的地方，以及一項本站補充，如下：

• （本站補充，非書中內容）電子態密度的同構結構：同樣的 |∇_k ε| 分母結構在電子能帶的態密度理論裡也成立，能帶極值與鞍點同樣產生 Van Hove 奇異點。不過 Kittel 在第 6、7 章並未寫出 D(ε) = ∫dS/|∇_k ε| 這個一般式，也沒有對電子提臨界點或 Van Hove 奇異點；全書只在 p.119（聲子）與 p.528（一維電子子能帶）兩處用到 Van Hove 這個名字（索引條目即寫「Van Hove singularities, 119, 528」）。
• 第 15 章（書 p.435）的光學吸收：兩個能帶平行（∇_k(ε_c − ε_v) = 0）處聯合態密度發散，吸收譜出現尖銳結構；原文明講這用的是與式 (5.37) 相同的論證。
• 第 18 章（書 p.528–529）的一維奈米結構：量子線每個子能帶的閾值處 D(ε) 以 (ε − ε_{i,j})^{−1/2} 發散，原文直接稱之為 van Hove singularities，並說這與三維的行為形成對比。碳奈米管的 STM 穿隧譜（書 p.530 Fig 10）量到的尖峰就是它。

注意第 18 章那個 (ε − ε_{i,j})^{−1/2} 的指數與習題 1(a) 的 (ω_m² − ω²)^{−1/2} 是同一個 −1/2。這不是巧合：兩者都是一維、都是在能量／頻率的極值處（群速度為零）、都是把 dK/dω 放在分母。「一維系統的態密度在能帶邊緣以 −1/2 冪次發散」是一條跨越聲子與電子的通則，本段的式 (15) 是它最早、最乾淨的出處。

本段的檢查清單：讀完應該能回答的問題，與最常見的四個錯誤我的判讀

讀完這一段應該能立刻回答的問題（每一題後面標出答案在書上的位置）：

1. 為什麼 Planck 分布的推導只需要「相鄰能階等距」這一個關於振子的性質？（p.107 式 (3) 的討論）
2. 式 (6) 的第二條恆等式為什麼可以用 x·d/dx 從第一條生出來？（p.108 式 (6)）
3. 固定端邊界條件下為什麼是 N − 1 個模式而不是 N 個或 N + 1 個？（p.109 式 (12) 後的討論）
4. 週期性邊界條件的 K 值間距為什麼是固定端的兩倍，總模式數卻一樣？（p.109–110）
5. Fig 4 的 N = 8 圓環上，K = 0 與 K = 8π/8a 為什麼各只貢獻一個模式？（p.110 Fig 4 圖說）
6. 式 (15) 在什麼情況下會發散，這在物理上代表什麼？（p.110 式 (15) 與 p.111 第一句，對照 p.119）

最常見的四個錯誤：

錯誤一：以為 K 的量子化是量子效應。 已在前面專節打掉。書 p.109 Fig 3 圖說的原句是「This quantization of K has nothing to do with quantum mechanics but follows classically from the boundary conditions that the end atoms be fixed.」

錯誤二：把 L/π 與 L/2π 混用。 這兩個數字對應兩套不相容的記帳系統。用 L/2π 就必須同時數 ±K；用 L/π 就只數正 K。混用會差兩倍。檢驗方法永遠是：算出來的模式總數必須等於（可動原子數 × 偏振數），對不上就是記帳錯了。

錯誤三：在區邊界重複計數。 ±π/a 是同一個點（差一個倒晶格向量 2π/a），只能算一次。Fig 5 畫了 9 個點卻只有 8 個模式，就是這件事的視覺提醒。在三維裡這個陷阱更隱蔽——布里元區的相對面要黏合起來算。

錯誤四：把式 (1) 當成晶格的總能量。 式 (1) 沒有零點能，它算的是 U(T) − U(0)。要算熱容（對 T 微分）沒問題；要算自由能、壓力、熱膨脹、Grüneisen 常數（對 V 微分），就必須把 Σ ½ℏω_K 補回去。書 p.129 習題 5(a) 的原話是「It is necessary to retain the zero-point energy to obtain this result.」

與後續章節的接點：