①主要內容與研究問題
第 4 章要解決的問題，Kittel 在 p.91 第一句就寫死了：We want to find the frequency of an elastic wave in terms of the wavevector that describes the wave and in terms of the elastic constants.——要的是 ω 對 K 的函數，輸入是彈性常數。前三章剛好把材料準備齊：第 1 章說原子站在哪裡，第 2 章給了倒晶格與布里淵區這個定義域，第 3 章把「原子之間有多硬」化成可量測的 C₁₁、C₄₄。第 4 章做的事，就是把這三樣東西送進牛頓第二定律，搖出一條色散關係。
有一件結構上的事值得先講明：這一章幾乎全是古典力學。十五頁裡只有 p.99 那半頁（QUANTIZATION OF ELASTIC WAVES）是量子的，而且只有一行 ε = (n + ½)ℏω。對本章逐頁做文字搜尋，heat、temperature、capacity、Debye、conductivity 五個字出現次數皆為 0（本站搜尋），唯一的 thermal 是 p.89 章首那句 Chapter 5 treats the thermal properties of phonons. 這不是疏漏，是分工：簡諧晶格的簡正模式頻率，古典與量子算出來完全相同——量子力學改變的是「每個模式裡裝多少能量」，不改變「有哪些模式、頻率各是多少」。所以第 4 章只負責交出模式清單，第 5 章才把 Planck 分佈接上去。認清這個分工，你就不會因為「這章是古典的」而懷疑它的結論在熱物性裡還能不能用。
章內其實只有兩個模型。單原子鏈（p.91–95）靠三個假設——高對稱方向整面同相、諧近似、只取最近鄰——把三維壓成一維，九步得到 ω = (4C/M)^½|sin ½Ka|。雙原子鏈（p.95–99）只多放一個原子進原胞，色散就分裂成聲學支與光學支，並在區邊界開出頻率間隙。兩個都是玩具模型，但它們夾著一條完全不依賴任何近似的結論，p.95 逐字：If there are p atoms in the primitive cell, there are 3p branches to the dispersion relation: 3 acoustical branches and 3p – 3 optical branches. 這條只來自自由度計數，對任何材料、任何力常數都成立。這一章如果只能帶走一件事，帶走這條。
它撐起哪幾章？第 5 章是直接的繼承者：ω(K) 決定態密度，p.118 式 (37) 明白寫成 D(ω) = V/(2π)³∮dS_ω/v_g——分母就是本章 (13)(14) 的群速；D(ω) 再決定熱容與熱導。第 5 章從頭到尾沒有一個量不是從本章這條曲線積出來的。連「為什麼需要兩個模型」的答案也在本章：聲學支近乎線性、光學支近乎平坦，一個函數裝不下兩種形狀，所以 p.112 用 Debye 處理聲學支、p.117 用 Einstein 近似光學支（該頁逐字：The Einstein model, however, is often used to approximate the optical phonon part of the phonon spectrum.）。往後：第 7 章的電子能隙與本章的頻率間隙是同一套區邊界簡併微擾；第 14 章的 LST 關係明說建立在「每原胞兩原子的立方晶體」，也就是本章的模型；第 16 章的軟模鐵電是本章 Problem 7 的 ω² < 0；第 12 章的磁振子量子化與 p.99 的聲子量子化數學同構；第 18 章把本章的模式計數搬進奈米結構。
最後一節（非彈性中子散射，p.100–102）是整章的轉向：前面都在「算」ω(K)，這一節說明它怎麼被「量」出來，而書中三張色散圖（Fig. 8a 的 Ge、Fig. 8b 的 KBr、Fig. 11 的 Na）全部是中子資料。這一章真正的產出不是那兩條公式，而是理論與實驗共用的那張圖。往後任何一篇聲子熱傳的論文，第一張圖仍然是它。
②③核心概念與公式 · 逐節深講
依書中實際小節切分，一節一頁。未標層級者為書中內容，頁碼與式號逐條核對；書中省略而由本站補上的代數步驟標為「補推導」。
⊘常見誤解與易錯點 （我的判讀）
- 「ω(K) 是晶體的材料性質，量一次就定了。溫度只決定有多少聲子被激發，不會改變色散關係本身——所以用 0 K 算出來的聲子譜，拿去做任何溫度的熱分析都可以。」在本章的世界裡這句話完全正確，這正是它危險的地方。整章建立在諧近似上（p.91：彈性能是相對位移的二次式），而諧近似的後果第 5 章 p.119 逐條列出來，其中一條逐字就是 The elastic constants are independent of pressure and temperature. 問題是真實晶體不是這樣。最直接的反證就在書裡：本章三張色散圖全都標了溫度——Fig. 8a 的 Ge 是 80 K、Fig. 8b 的 KBr 是 90 K、Fig. 11 的 Na 是 90 K。如果 ω(K) 真與溫度無關，標溫度就沒有意義。第 3 章 p.84 Table 11 更是白紙黑字：鋁的 C₁₁ 從 0 K 的 1.143 掉到 300 K 的 1.068（10¹¹ N/m²），−6.6%；鉛掉 10.8%、銅掉 4.4%。彈性常數隨溫度變，聲速就隨溫度變，ω(K) 當然跟著移。把它拆成「對 Cth 影響多大」與「對 Rth 影響多大」兩問，答案完全不同。對 Cth：可以忽略。由 ω ∝ √C，鋁的 C₁₁ 掉 6.6% 只讓頻率尺度掉 √0.934 = 3.3%（本站計算）；而 300 K 的矽 dlnC_v/dlnθ = −0.42（本站由書 p.116 的 θ = 645 K 算 Debye 函數的對數導數），所以 θ 移動 3.3% 只改變 C_v 約 1.4%——遠小於「有效受熱體積是多少」這個幾何不確定度。對 Rth：不可忽略，而且是質變。非諧項不只是讓 ω 漂移，它是聲子平均自由程 ℓ 存在的唯一原因；p.119 同一張清單的第一條逐字寫著 Two lattice waves do not interact; a single wave does not decay or change form with time.——嚴格諧晶體的 ℓ = ∞、κ = ∞、Rth = 0。操作判準：聽到「用 0 K 的 DFT 聲子譜算 500 K 的熱容」，那大致沒問題（誤差百分之一二）；聽到「用同一份聲子譜算 500 K 的熱導」，就要追問三階力常數怎麼處理的——沒有三階項，熱導根本是無限大。
- 「矽的光學支頻率很高，室溫下 ℏω ≫ k_BT，還沒被激發，所以室溫熱容基本上是三條聲學支在貢獻。光學支要到高溫才需要考慮。」「頻率高所以沒被激發」這個推論漏掉了一個定量檢查：要凍結，需要的是 T ≪ θ_E = ℏω/k_B，不是「頻率看起來很大」。用書中真的印出來的數字算一次就穿幫了。Fig. 8a 的鍺（本章唯一有色散資料的半導體）光學支落在 8.7–9.15 × 10¹² Hz（該圖為影像，此為目視讀值，見本章第 3 段），換算 θ_E = hν/k_B = 418–439 K，ℏω 只有 36–38 meV（本站計算）。300 K 時 T/θ_E ≈ 0.7，代進 p.117 的 Einstein 式 C = Nk_B x²e^x/(e^x−1)²，每個光學模已經帶到古典值的 0.84–0.85（本站計算）。這離「凍結」差得很遠。矽的光學支高一些（Γ 點約 15.3 THz，θ_E ≈ 734 K，外部值、書中未給），300 K 的 Einstein 因子仍有 0.62。
做一次明確的分支拆解，不要靠直覺。把 3 條聲學支當 Debye、3 條光學支當 Einstein（這正是 p.112 與 p.117 的分工）；注意聲學支專屬的 Debye 溫度不是書上表裡那個 θ：書中 θ 是把全部 3N_atom 個模式一起當 Debye 算的，只留 3 條分支時 θ_acoustic = θ·(N_cell/N_atom)^⅓ = θ/2^⅓（本站推導）。矽 → θ_a = 512 K，鍺 → θ_a = 297 K。結果（本站計算，單位 J·cm⁻³·K⁻¹）：
300 K 聲學 300 K 光學 光學佔比 Si（θ_a = 512 K、θ_E = 734 K） 0.90 0.64 42% Ge（θ_a = 297 K、θ_E = 418 K） 0.87 0.78 47% 誠實說明：這個兩段式拆解把矽的總熱容算成 1.54，比單一 Debye(645) 的 1.66 與文獻塊材值 1.64（外部值）低約 6%，所以「42%」的不確定度大約 ±5 個百分點。但方向不會錯：室溫下光學支貢獻接近一半的熱容。判準只有一句——先算 θ_E = hν/k_B，再看 T/θ_E。T/θ_E > 0.3 就不能說「凍結」；真正凍結要到 T < θ_E/5，矽大約是 150 K 以下。
- 「Fig. 6 顯示群速在區邊界降到零，代表區邊界的聲子既不搬能量、也幾乎沒有存在感。所以態密度 D(ω) 在 ω 接近 ω_max 時應該趨近零，高頻端可以忽略。」把「單一模式搬多少能量」與「有多少模式」混為一談了，而這兩件事在 D(ω) 裡是互為倒數的關係。p.118 式 (37) 寫得毫不含糊：D(ω) = V/(2π)³∮ dS_ω/v_g——群速在分母。v_g → 0 讓態密度發散，不是消失。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. 4 那條 |sin ½Ka| 在區邊界的「平頂」，到了第 5 章就是 D(ω) 的一根尖峰；第 5 章 Problem 1 要你親手把一維的 D(ω) 算出來，答案在 ω → ω_max 處發散。
換一條讀圖規則：在 ω(K) 圖上，曲線越平的地方，態密度越大。這條規則一口氣解釋三件事。(a) 區邊界的 Van Hove 奇點。(b) 為什麼近乎水平的光學支能貢獻巨大的態密度尖峰，也就是為什麼 p.117 說 Einstein 模型「常用來近似光學聲子那一段」——平坦 = 所有模式同頻 = Einstein 的前提。(c) 對熱分析最要緊的推論：熱容權重最大的地方，正好是群速最小、最不導熱的地方。
把這條規則拿去檢查 Debye 模型，還可以量化它錯在哪。用本章 (9) 的真實色散對比線性外推（本站計算）：
K/K_max 真實 ω/ω_max Debye 線性外推 v_g 真值 / 聲速 0.25 0.383 0.393（+3%） 0.92 0.50 0.707 0.785（+11%） 0.71（Debye 高估 41%） 1.00 1.000 1.571（+57%） 0（Debye 高估無限倍） 所以 Debye 模型的系統性偏差是：把頻率往上推、把群速往上推、把區邊界那根態密度尖峰完全抹平。低溫時只用到 K ≈ 0 的部分所以準（T³ 律），高溫時被 Dulong–Petit 上限救了也準，中溫段（T ≈ θ/3 到 θ）最糟——對矽正好是 200–650 K，也就是你整個工作溫度區間。
- 「矽的立方晶胞有 8 個原子，依 p.95 的 3p 規則，矽應該有 3 × 8 = 24 條聲子分支。」
3p 裡的 p 是原胞（primitive cell）內的原子數，不是慣用立方晶胞（conventional cell）內的原子數。p.96 逐字：With p atoms in the primitive cell and N primitive cells, there are pN atoms. 而腳註 1 定義的 N = V/V_c 用的也是原胞體積。矽是金剛石結構——書 p.16 逐字：The space lattice of diamond is face-centered cubic. The primitive basis of the diamond structure has two identical atoms at coordinates 000 and ¼¼¼ ... the conventional unit cube of the diamond structure contains 2 × 4 = 8 atoms. There is no way to choose a primitive cell such that the basis of diamond contains only one atom. 所以 p = 2，分支數是 6：3 條聲學（1 LA + 2 TA）+ 3 條光學（1 LO + 2 TO）——正是 p.95 拿鍺與 KBr 當例子時寫的 six branches: one LA, one LO, two TA, and two TO。
這個錯誤有兩層。總自由度會被自由度守恆救回來（改用慣用晶胞則 p 大 4 倍、N 小 4 倍，3pN 不變），所以你可能一直沒發現；但分支數與布里淵區的形狀全錯。fcc 的倒晶格是 bcc，第一區是截角八面體；若誤把慣用立方晶胞當原胞，你會得到一個立方體的區，沿 [111] 的區邊界位置差 √3 倍。
三連問，每次遇到新材料都跑一次：(i) Bravais 格子是什麼？（矽：fcc，書 p.16）(ii) 基底有幾個原子？（矽：2 個，在 000 與 ¼¼¼，書 p.16）(iii) 3p = ?（矽：6）。
第二道獨立檢查在圖上：Fig. 8a 的圖說逐字寫 K_max = (2π/a)(½ ½ ½)，那是 fcc 倒晶格（bcc）的 L 點，其中 a 是慣用立方晶格常數。若你腦中是「簡單立方原胞」的圖像，區邊界會落在 (2π/a)(½ 0 0)，跟圖對不起來。第三道檢查用本章自己的 (25)：等質量時區邊界兩個根重合、頻隙關閉，所以任何雙原子鏈的結果代入 M₁ = M₂、a → 2a′ 都必須退化回單原子鏈的 (9)——這條同時抓得到 a 的定義錯誤（p.97 逐字：the distance between nearest identical planes, not nearest-neighbor planes）。
- 「這一章通篇是一維鏈，真實晶體是三維的、要對角化 3×3 的動力學矩陣。所以本章算出來的數字（力常數、ω_max、區邊界頻率）只有教學意義，沒有定量價值。」這句話高估了「一維化」的代價，同時低估了「最近鄰」的代價——兩個近似被搞反了。一維化在這三個方向上不是近似，是對稱性的結果。p.91 逐字：When a wave propagates along one of these directions, entire planes of atoms move in phase with displacements either parallel or perpendicular to the direction of the wavevector. 沿 [100]、[110]、[111]，整面同相是晶體對稱性強制的，於是每個偏振各自嚴格滿足一條一維差分方程。真正的近似是「只取最近鄰」與「所有最近鄰平面對共用同一個 C」。而這兩個近似的成敗有實測可查：KBr 沿 [111] 用一個力常數就把區邊界四個頻率預測到 5% 以內（見本章第 3 段的數值練習）；鍺卻失敗，因為 [111] 的平面間距以 3:1 交錯，「等力常數」假設垮掉。至於最近鄰夠不夠遠，p.94 直接給了警告：Interactions have been found between planes of atoms separated by as many as 20 planes.
先問「一維化合法嗎」，再問「近似夠準嗎」，順序不能反。第一問的判準：在這個傳播方向上，一個平面內是否只含一種原子、且面內所有原子等價？是 → 一維化嚴格（這就是為什麼書中三張色散圖只畫 [001]、[110]、[111]，以及為什麼 p.97 指定 NaCl 結構用 [111]、CsCl 結構用 [100]）；否 → 連建模都不對，此時任何數值都沒有意義。第二問的判準：有效交互作用延伸幾層？方法書上就有——p.94–95 的 (16a)(17)：把實驗量到的整條 ω²(K) 做傅立葉反演，一次把 C₁、C₂、C₃… 全部解出來，看第幾項才衰減掉。（注意 (17) 明文限定 monatomic basis，雙原子鏈不能直接套。）
定量價值的實證：本章第 1 段用第 3 章 p.84 的鋁資料，只靠聲速定一個參數，算出 ℏω_max/k_B = 471 K，對上書 p.116 的 Debye 溫度 428 K，差 10%。一個單參數玩具模型能做到這樣，不叫「沒有定量價值」。它的極限在別處：只保留最近鄰時，ω_max/(v·K_max) 被鎖死為 2/π ≈ 0.64，任何材料都一樣，沒有調整空間——這個比值對不上實驗，就是在告訴你必須加遠程力常數。
④本章最重要的結論
晶格的離散性使彈性波的色散關係在布里淵區邊界飽和（群速度為零、解成為駐波），因此所有彈性波都可以用第一布里淵區內的波向量描述；長波極限下 ω 正比於 K，回歸連續介質彈性理論。每原胞 p 個原子時色散關係分裂成 3 個聲學支與 3p − 3 個光學支，兩原子基底在區邊界出現頻率間隙。晶格振動的能量量子化為聲子，能量為 (n + ½)ℏω；聲子帶有 ℏK 的準動量並遵守含倒晶格向量 G 的波向量選擇律，這使得中子非彈性散射成為量測 ω(K) 的標準實驗手段。
書中逐字
"The quantum unit of a crystal vibration is a phonon." (p.102)
"However, a phonon does not carry physical momentum." (p.100)
"The transmission velocity of a wave packet is the group velocity" (p.94)
⑤與前後章的關係 （我的判讀）
用 Ch2 的 Brillouin 區當定義域、用 Ch3 的彈性常數當長波極限的檢核。它產出的色散關係 ω(K)、以及每原胞兩原子時分裂出的聲學支與光學支，是 Ch5 態密度 D(ω) 的唯一輸入。Ch12 的磁振子量子化與本章的聲子量子化在數學上同構。