這一段要補完的那個洞我的判讀
第 5 章講熱傳導的主幹只有一條式子，就是 p.122 的 (42)：K = ⅓Cvℓ。C 在 p.107–115 已經算完（Debye 模型、T³ 定律），v 是聲速、在 Debye 近似裡當常數，唯一沒交代的就是 ℓ——聲子的平均自由程。整個 p.121–128 的後半段，做的就是回答「ℓ 由什麼決定」。
書把答案拆成兩半，正好對應本段的兩個小節：
Umklapp processes（p.125–126） 回答的是「聲子彼此相撞時，ℓ 由什麼決定」。這裡有一個非常違反直覺的結論：不是所有的聲子碰撞都能產生熱阻。三聲子碰撞 K₁+K₂=K₃ 撞得再頻繁，熱阻仍然是零。真正製造熱阻的，是那些必須外加一個倒晶格向量 G 才能成立的碰撞。
Imperfections（p.126–128） 回答的是「當聲子彼此幾乎不撞時，ℓ 由什麼決定」。低溫下 U 過程被指數壓制，接手的是幾何——晶體邊界、試片直徑、同位素質量亂序。
這兩半合起來，就是整條熱導率對溫度的曲線：低溫端 K ∝ T³ 由尺寸效應撐著、中間有一個尖峰、高溫端 K ∝ 1/T 由 U 過程壓下來。Fig 18 那條氟化鈉曲線就是這條「駝峰」的實驗版本。讀完這一段，你應該要能只看一眼 K(T) 曲線就說出每一段是誰在主宰。
（附帶說明：第 5 章的目次（書本 p.105）列到 Imperfections 126 之後就直接是 PROBLEMS 128，本章沒有 SUMMARY 節。這與第 4 章不同，不必去找。）
Umklapp 過程：開場的那一句反直覺
小節標題 Umklapp Processes 印在 p.125 下半，開場第一句就是整段的重點：
「The important three-phonon processes that cause thermal resistivity are not of the form K₁ + K₂ = K₃ in which K is conserved, but are of the form」（逐字，p.125）
接著就是編號 (47) 的式子 K₁ + K₂ = K₃ + G，然後：
「where G is a reciprocal lattice vector (Fig. 17). These processes, discovered by Peierls, are called umklapp processes.」（逐字，p.125）
請注意這句話的語氣有多強：它不是說「U 過程也有貢獻」，而是說造成熱阻的重要三聲子過程「不是」K 守恆的那一種。N 過程不是次要，是根本不產生熱阻。
Peierls 是 1929 年提出這件事的（書上只給名字，沒給年份）。緊接著 Kittel 補了一句提醒：
「We recall that G may occur in all momentum conservation laws in crystals. In all allowed processes of the form of (46) and (47), energy is conserved.」（逐字，p.125）
校註（頁碼與式號已回原文核對）：這句話裡的 (46) 在 p.123 印的是 J = Σ_K n_K ℏK（聲子氣的總動量定義），而 K₁+K₂=K₃ 在 p.123 印的編號是 (45)。因此「processes of the form of (46) and (47)」按字面讀是對不上的——依上下文，作者想指的應該是 (45) 與 (47) 這兩類過程。這是書上的排版／引用疏漏，不是你讀錯。本講解後續一律用 (45) 指 N 過程、(47) 指 U 過程。
「umklapp」是德文，意思是「翻折過去」。Kittel 自己在第 6 章 p.152 講電子的 umklapp 散射時才解釋了這個字：「(In an umklapp process the wavevector of one particle may be "flipped over.")」（逐字，p.152）——第 5 章這裡沒有解釋字義，所以不少人讀到最後都不知道這個字在講什麼。它講的就是「波向量被折回第一區」這個動作本身。
N 過程（normal process）的波向量守恆
K₁ + K₂ = K₃
(45) p.123K₁、K₂：碰撞前被吸收的兩個聲子的波向量；K₃：碰撞後被放出的聲子的波向量。三者都落在第一布里淵區內。p.126 明確定義：「Collisions in which G = 0 are called normal processes or N processes.」（逐字）
U 過程（umklapp process）的波向量守恆
K₁ + K₂ = K₃ + G
(47) p.125K₁、K₂：被吸收的兩個聲子；K₃：被放出的聲子，必須落在第一布里淵區內；G：任一非零的倒晶格向量。在 Fig 17(b) 的二維正方晶格裡，圖上畫的 G 的長度是 2π/a（a 為晶格常數）、方向平行於 K_x 軸——這是 Fig 17 圖說的逐字內容：「The reciprocal lattice vector G as shown is of length 2π/a, where a is the lattice constant of the crystal lattice, and is parallel to the K_x axis.」
為什麼 G 一定會冒出來：第一布里淵區論證
p.126 開頭那一段是全章最值得逐字讀的論證。它分四步：
第一步，先把結論一般化。「We have seen examples of wave interaction processes in crystals for which the total wavevector change need not be zero, but may be a reciprocal lattice vector. Such processes are always possible in periodic lattices.」（逐字，p.126）——這不是聲子的專利。第 2 章的 X 光繞射（Laue 條件 Δk = G）、第 4 章的中子非彈性散射（k + G = k′ ± K），都是同一件事的不同面貌：在週期晶格裡，嚴格守恆的從來不是波向量本身，而是「波向量模 G」。
第二步，指出聲子的情況更強。「The argument is particularly strong for phonons: the only meaningful phonon K's lie in the first Brillouin zone, so that any longer K produced in a collision must be brought back into the first zone by addition of a G.」（逐字，p.126）
這裡的 particularly strong 值得停一下。為什麼對聲子「特別強」？因為對光子或電子，區外的 k 至少還對應到高能帶、有獨立的物理內容；但對聲子，區外的 K 完全沒有新的物理——第 4 章 p.93 已經證明，K 與 K + G 描述的是一模一樣的原子運動（因為原子只存在於離散的格點上，取樣定理限制了你能分辨的最短波長）。所以 K₁+K₂ 落在區外時，你沒有選擇：那個向量所代表的模態，就是區內的 K₃。加 G 不是一個「可以發生也可以不發生」的選項，而是記帳方式的必然。
第三步，給出最關鍵的一句話。「A collision of two phonons both with a negative value of K_x can by an umklapp process (G ≠ 0), create a phonon with positive K_x.」（逐字，p.126）
這一句是整個 U 過程的靈魂。兩個往 −x 跑的聲子，撞完之後可以生出一個往 +x 跑的聲子。這在自由空間的氣體分子碰撞裡是絕對不可能的（動量守恆禁止），但在晶格裡完全合法，因為晶格整體吸收了 ℏG 這筆動量。也正因為如此，熱流的方向可以被反轉——這就是熱阻的微觀機制。
第四步，命名。「Umklapp processes are also called U processes.」「Collisions in which G = 0 are called normal processes or N processes.」（逐字，p.126）
校註（回 PDF 字體層核對）：這兩句命名 N/U 的話裡，「Collisions in which G = 0」的 G 用的是粗體向量，與 (47) 一致。順帶一個排版細節：同一段稍早的「by an umklapp process (G ≠ 0)」裡，那個 G 在字體層是斜體純量（NewCaledonia-Italic），與三行後粗體的「G = 0」不一致——這是書上的排版小疏漏，語意不變。
Fig 17 到底畫了什麼（網頁無圖，逐一描述）我的判讀
Fig 17 印在 p.125 的中段，是兩張並排的二維 K 空間圖，共用同一套圖說。以下是對原圖的逐項目視描述（已放大原頁確認）。
共同元素：每張圖的中央有一組 K_x（水平，向右為正）與 K_y（垂直，向上為正）座標軸。以原點為中心畫了一個淺灰色的正方形，那就是二維正方晶格的第一布里淵區。圖說逐字：「The square in each figure represents the first Brillouin zone in the phonon K space; this zone contains all the possible independent values of the phonon wavevector.」
箭頭的讀法（這是最容易看錯的地方，圖說特別交代了）：「Vectors K with arrowheads at the center of the zone represent phonons absorbed in the collision process; those with arrowheads away from the center of the zone represent phonons emitted in the collision.」（逐字）——也就是說，箭頭尖端指向原點的 = 被吸收的入射聲子；箭頭尖端背離原點的 = 被放出的聲子。所有聲子 K 向量都以原點為公共頂點，只是箭頭朝向不同；(b) 裡的 G 是例外，它畫在 K₃ 的箭尖與 K₁+K₂ 的箭尖之間（見下文對 (b) 的描述）。
(a) Normal 過程：K₁ 從左上方畫進來、箭頭停在原點（吸收）；K₂ 從左下方畫進來、箭頭停在原點（吸收）；K₃ 從原點畫出去、指向右下方（放出），整支箭頭都還在灰色正方形內。圖下標 (a)，圖說標示「(a) Normal K₁ + K₂ = K₃」。把 K₁ 與 K₂ 這兩個向量（方向都是朝右下／右上）頭尾相接，合向量就是 K₃，仍在區內，不需要 G。
(b) Umklapp 過程：K₁ 從更遠的左上方畫進來、箭頭停在原點；K₂ 從上方偏左畫進來、箭頭停在原點。接著有一條虛線箭頭從原點畫向右下方，標記為 K₁ + K₂，它的尖端已經穿出灰色正方形的右邊界——這就是「碰撞產生了過長的 K」。然後圖上畫了一條水平向右的實線箭頭標記為 G，從 K₃ 的尖端指到 K₁+K₂ 的尖端；而 K₃ 是從原點畫向左下方、落在正方形內的實線箭頭。
把 (b) 的幾何關係讀出來就是：K₃ + G = K₁ + K₂，也就是 (47)。而重點在於：K₁ 與 K₂ 的 x 分量都是正的（兩支入射箭頭都是從左往右畫到原點），但 K₃ 的 x 分量是負的。圖說逐字：「We see in (b) that in the umklapp process the direction of the x-component of the phonon flux has been reversed.」
注意 (b) 圖與 p.126 正文的方向剛好相反：正文舉的例子是「兩個 K_x 為負的聲子生出一個 K_x 為正的聲子」，圖裡畫的是「兩個 K_x 為正的聲子生出一個 K_x 為負的聲子」。這兩個敘述是同一件事的鏡像（把 G 反向即可），不是矛盾，但初讀時很容易以為自己看錯圖。就熱傳導而言，重要的是後者：熱流從熱端往冷端流（+x），U 過程把它折回 −x，這才叫熱阻。
最後一行圖說是全圖最容易被跳過、卻最重要的一句：「For all processes, N or U, energy must be conserved, so that ω₁ + ω₂ = ω₃.」（逐字）——U 過程折的是波向量，不是能量。能量守恆在 N 與 U 完全一樣、沒有任何鬆動。這一點下一段會再強調一次。
Fig 16d 與那條「熱端—冷端」的棒子我的判讀
p.125 的最上方是 Fig 16 系列的最後一張 Fig 16d（前三張 a、b、c 在 p.124，屬於前一段的範圍）。這張圖是一個水平長方形，左端標 Hot end、右端標 Cold end，中間標 U processes，框內畫滿了朝向四面八方、幾乎完全隨機的白色小箭頭。
對照 Fig 16c（p.124，只有 N 過程）就看得出設計意圖：16c 的箭頭是整齊地一致向右（標 Net flow of phonons，左端是 Phonon source、右端是 Phonon sink），16d 的箭頭則是方向全亂。
圖說逐字（p.125）：「In U processes there is a large net change in phonon momentum in each collision event. An initial net phonon flux will rapidly decay as we move to the right. The ends may act as sources and sinks. Net energy transport under a temperature gradient occurs as in (b).」
這段圖說的三句話各有各的意思：
第一句講 U 過程每次碰撞的動量變化「大」——大到 ℏ|G| 的量級。
第二句「An initial net phonon flux will rapidly decay」是熱阻的定義本身：你在左端打進一束有方向的聲子流，它往右走幾個 ℓ 就把方向記憶忘光了。這正是擴散（diffusive）而非彈道（ballistic）輸運，也正是 p.121 在定義 (41) j_U = −K dT/dx 時所強調的——「This form implies that the process of thermal energy transfer is a random process.」（逐字，p.121）。沒有 U 過程，就沒有隨機化，就沒有 dT/dx 這個形式。
第三句 把讀者導回 Fig 16b（p.124 的封閉管子類比）：真正的熱傳導實驗裡沒有淨質量（聲子數）流，只有淨能量流。
Fig 16 整個系列在說的是同一個類比：a 是無摩擦管內的氣體（動量守恆 → 熱阻為零）；b 是兩端封死的管子（正常的熱傳導定義）；c 是只有 N 過程的晶體（等同 a，熱阻為零）；d 是有 U 過程的晶體（動量被晶格吃掉，終於有熱阻）。Kittel 花四張圖講這件事，就是因為「碰撞頻繁 ≠ 有熱阻」這一點太反直覺。
把動量收支算一遍：為什麼 N 過程守得住、U 過程守不住補推導
書上在 p.123 只寫了一句「the change in J is K₃ − K₂ − K₁ = 0」，在 p.125–126 也只說 U 過程有「large net change in phonon momentum」，中間的帳沒有列出來。這裡把兩種過程的動量收支完整列一次——這是判斷「哪些散射機制會產生熱阻」的通用工具，之後遇到雜質散射、邊界散射都可以套。
推完之後會得到一個可以直接用的判準：一個散射過程能不能產生熱阻，看的不是它改不改變單一聲子的狀態，而是看它改不改變聲子氣的總動量 J。
- 定義（p.123 的 (46)）：聲子氣的總動量 J = Σ_K n_K ℏK，其中 n_K 是波向量 K 的模態上的聲子數。注意這是「晶體動量」（crystal momentum），不是真動量；p.100 對此有過一句提醒——真正嚴格守恆的是整個系統的動量。
- N 過程的收支：碰撞前後聲子佔據數的變化是 Δn_{K₁} = −1、Δn_{K₂} = −1、Δn_{K₃} = +1。於是 ΔJ = ℏ(K₃ − K₂ − K₁)。
- 代入 (45) K₁ + K₂ = K₃，得 ΔJ = ℏ(K₃ − K₂ − K₁) = 0。每一次 N 碰撞，J 分毫不差地守恆。
- U 過程的收支：同樣 ΔJ = ℏ(K₃ − K₂ − K₁)，但這次要代入 (47) K₁ + K₂ = K₃ + G，即 K₃ = K₁ + K₂ − G。
- 得 ΔJ = ℏ(K₁ + K₂ − G − K₂ − K₁) = −ℏG。這筆動量沒有消失，它被整個晶體（所有原子的質心）接走了。因為晶體的質量是宏觀的，它接走 ℏG 之後速度變化可以忽略、能量變化更可以忽略——所以能量守恆完全不受影響，這就是「U 過程折波向量但不折能量」的真正原因。
- 把碰撞率寫進動量的演化方程：dJ/dt = −(1/τ_U)J，其中 1/τ_U 是 U 碰撞的頻率（N 碰撞完全不出現在這條式子裡，因為它對 J 的貢獻是零）。
- 解得 J(t) = J(0)exp(−t/τ_U)：聲子流的淨動量只靠 U 過程衰減。對應到空間上，聲子流在走過 ℓ = v τ_U 的距離之後才失去方向記憶。
- 把這個 ℓ 代回 (42) K = ⅓Cvℓ，就得到 p.126 最後那一句結論：進入 (42) 的 ℓ 是 umklapp 碰撞的平均自由程。
高溫行為：所有模態都醒著，熱阻率正比於 T
p.126 第二段處理高溫。原文逐字：
「At high temperatures T > θ all phonon modes are excited because k_BT > ℏω_max. A substantial proportion of all phonon collisions will then be U processes, with the attendant high momentum change in the collision.」
這兩句話的邏輯鏈是：
k_BT > ℏω_max ⟹ 所有模態都被激發 ⟹ 大 K 的聲子（K ~ ½G）的數目不再受 Boltzmann 因子壓制 ⟹ 能參與 U 過程的聲子俯拾即是 ⟹ 相當比例的碰撞是 U 過程。
注意 ω_max 就是色散關係的最高頻率，在 Debye 模型裡即 ω_D，而 ℏω_D = k_Bθ（p.112 的 (27) 定義）。所以「T > θ」與「k_BT > ℏω_max」講的是同一件事。
接下來這一句是結論：
「In this regime we can estimate the thermal resistivity without particular distinction between N and U processes; by the earlier argument about nonlinear effects we expect to find a lattice thermal resistivity ∝ T at high temperatures.」（逐字，p.126）
「without particular distinction between N and U」值得畫線。它的意思是：高溫下 U 過程佔的比例夠高，於是「總碰撞率」與「U 碰撞率」只差一個 O(1) 的常數因子，你可以偷懶地用總碰撞率去估熱阻而不會估錯冪次。但這是高溫才有的便宜，低溫下絕對不能這樣做——這正是下一段要講的。
「the earlier argument about nonlinear effects」指的是 p.123 的那一段：「The theory of the effect of anharmonic coupling on thermal resistivity predicts that ℓ is proportional to 1/T at high temperatures, in agreement with many experiments.」以及它的物理解釋——「at high temperature the total number of excited phonons is proportional to T. The collision frequency of a given phonon should be proportional to the number of phonons with which it can collide, whence ℓ ∝ 1/T.」（皆逐字，p.123）
校註：本專案使用的 PDF 文字層把「∝」這個符號存成自訂字符 /L50524，抽取時不會被解碼器處理，會直接印出代碼。上面兩處的「∝」都是回原頁算繪目視確認後還原的。
補齊高溫 K ∝ 1/T 的完整推導鏈補推導
書上把這條鏈拆在兩處（p.123 給 ℓ ∝ 1/T，p.126 給熱阻率 ∝ T），中間的代數沒有連起來。這裡補完，順便把每一步用到的假設標出來——考試最常考的就是「高溫下 K 的溫度依賴」，而漏掉「C 在高溫是常數」這一步就會答錯。
- 高溫下（T > θ）每一個模態的平均聲子數由 Planck 分布給出：⟨n⟩ = 1/(e^{ℏω/k_BT} − 1)。展開分母：當 ℏω ≪ k_BT 時 e^{ℏω/k_BT} − 1 ≈ ℏω/k_BT，故 ⟨n⟩ ≈ k_BT/ℏω。
- 全晶體的總聲子數 N_ph = Σ_K ⟨n(ω_K)⟩ ≈ k_BT Σ_K 1/(ℏω_K)。右邊的求和與溫度無關，所以 N_ph ∝ T。這就是 p.123 所說的「at high temperature the total number of excited phonons is proportional to T」。
- 一個給定聲子的碰撞頻率正比於它能撞上的對象數目：1/τ ∝ N_ph ∝ T。（這一步的物理假設是「散射截面本身與溫度無關」，也就是非諧耦合常數不隨溫度變。）
- 平均自由程 ℓ = vτ ∝ 1/T。這就是 p.123 的 ℓ ∝ 1/T。
- 高溫下熱容 C 趨於 Dulong–Petit 值 3Nk_B（p.113：「At T ≫ θ the heat capacity approaches the classical value of 3Nk_B」），與溫度無關；Debye 模型裡 v 也取常數。
- 代入 (42)：K = ⅓Cvℓ ∝ (常數)(常數)(1/T)，故 K ∝ 1/T。
- 晶格熱阻率 W ≡ 1/K ∝ T。這就是 p.126 那句「a lattice thermal resistivity ∝ T at high temperatures」。
- 反過來檢查一致性：p.122 Table 2 的石英，0 °C 時 ℓ = 40 Å、−190 °C（83 K）時 ℓ = 540 Å。若 ℓ ∝ 1/T 成立，273/83 = 3.3 倍，但實測是 540/40 = 13.5 倍——差了 4 倍。這不是表錯，而是 83 K 已經低於石英的 θ，ℓ 已經開始走指數而非 1/T。這張表本身就在提醒你 1/T 只在 T > θ 成立。
低溫行為：U 過程需要「半個 G」與「半個 k_Bθ」
p.126 第三段是本節的技術核心，也是最常被略讀的一段。原文逐字：
「The energy of phonons K₁, K₂ suitable for umklapp to occur is of the order of ½k_Bθ, because each of the phonons 1 and 2 must have wavevectors of the order of ½G in order for the collision (47) to be possible.」
「If both phonons have low K, and therefore low energy, there is no way to get from their collision a phonon of wavevector outside the first zone. The umklapp process must conserve energy, just as for the normal process.」
「At low temperatures the number of suitable phonons of the high energy ½k_Bθ required may be expected to vary roughly as exp(−θ/2T), according to the Boltzmann factor. The exponential form is in good agreement with experiment.」
三段話，三件事：
第一，「½G」的來源是幾何。要讓 K₁+K₂ 走出第一區，兩個入射聲子必須各自貢獻夠大的波向量。第一區的「半徑」大約是 ½G，所以最經濟的組合是 K₁ ≈ K₂ ≈ ½G（下一個 block 會把這個估計做細）。
第二，「½k_Bθ」是把 ½G 翻成能量。在 Debye 圖像裡 ω = vK 而 ℏω_D = k_Bθ 對應 K_D，K ≈ ½K_D 就對應 ℏω ≈ ½k_Bθ。
第三，Boltzmann 因子給出 exp(−θ/2T)。要找到一個能量為 ½k_Bθ 的聲子，機率權重是 exp(−(½k_Bθ)/k_BT) = exp(−θ/2T)。
Kittel 用的措辭很謹慎：「may be expected to vary roughly as」。這是估計，不是嚴格結果——嚴格的低溫 U 過程理論（Peierls）給出的前置冪次因子與晶體結構、聲子譜細節有關。但 Kittel 直接說：「The exponential form is in good agreement with experiment.」（逐字）——指數形式這件事本身，實驗是站在理論這一邊的。
這一段也順手強調了一句常被忽略的話：「The umklapp process must conserve energy, just as for the normal process.」（逐字）。搭配 Fig 17 圖說的 ω₁ + ω₂ = ω₃，Kittel 在一頁之內把「U 過程不破壞能量守恆」說了兩次。會這樣重複，就是因為學生太常誤以為 U 過程是某種「不守恆的碰撞」。
把「½G」與「½k_Bθ」算出來（書上只給結論）補推導
p.126 只說「each of the phonons 1 and 2 must have wavevectors of the order of ½G」，沒有說明為什麼是 ½ 而不是別的數。這個 ½ 其實可以用一維（或沿 G 方向的一維投影）幾何幾行算出來，而且算的過程會讓「為什麼低 K 聲子撞不出 U 過程」變得完全顯然。
先講結論會用在哪：算出來的 ½k_Bθ 直接決定了 Boltzmann 指數 exp(−θ/2T)，而這個指數是低溫熱導率上升沿的全部物理。若把 ½ 換成 1，指數變成 exp(−θ/T)，與實驗差很遠。
- 設定：沿 G 方向做一維討論。第一布里淵區在這個方向上是 −½G ≤ K ≤ ½G（區邊界在 ±½G = ±π/a，這是第 4 章 p.93 的結果）。
- U 過程的最低要求：K₁ + K₂ 必須落在區外，也就是 |K₁ + K₂| > ½G。
- 兩個聲子都在區內，故 K₁ ≤ ½G 且 K₂ ≤ ½G，於是 K₁ + K₂ ≤ G。可能的區間是 ½G < K₁ + K₂ ≤ G——U 過程的窗口存在，但很窄。
- 在「兩個聲子能量相當」這個最可能的配置下（因為要同時找到兩個高能聲子，機率由兩者能量之和決定，能量之和固定時對稱分配並不改變總和，但要跨過 ½G 這個門檻，最省能量的方式是兩者各出一半），取 K₁ ≈ K₂ ≈ K，門檻條件變成 2K > ½G，即 K > ¼G。書上取的 ½G 是「量級」的說法（p.126 用的字是 of the order of），介於門檻值 ¼G 與上限 ½G 之間。
- 把波向量翻成能量：Debye 近似下 ω = vK；區邊界對應的 Debye 截止為 K_D，且由 (27) ℏω_D = k_Bθ。在這個一維幾何裡 ½G = π/a 與 K_D 同量級，故取 K ≈ ½K_D。
- 於是 ℏω = ℏvK ≈ ½ℏvK_D = ½ℏω_D = ½k_Bθ。這就是書上說的「½k_Bθ 量級」——原文完整的句子是「The energy of phonons K₁, K₂ suitable for umklapp to occur is of the order of ½k_Bθ」（逐字，p.126）。
- 低溫下這種聲子的數目：Planck 分布在 ℏω ≫ k_BT 時退化為 Boltzmann 因子，⟨n⟩ = 1/(e^{ℏω/k_BT} − 1) ≈ e^{−ℏω/k_BT}。代入 ℏω = ½k_Bθ 得 ⟨n⟩ ≈ exp(−θ/2T)。
- 代回輸運：U 碰撞率 1/τ_U ∝ ⟨n⟩ ∝ exp(−θ/2T)，故 ℓ_U ∝ exp(+θ/2T)。這是一個隨溫度下降而爆炸性增長的量——T = θ/10 時 exp(5) ≈ 150 倍，T = θ/20 時 exp(10) ≈ 2.2×10⁴ 倍。
- 代入 (42)：K = ⅓Cvℓ_U ∝ T³·exp(θ/2T)。T³ 往下壓、指數往上拉，指數贏——這就是為什麼 K(T) 在降溫時會先急遽上升。但 ℓ_U 不可能無限長：一旦它長到與試片直徑相當，就換幾何接手，K 回到純 T³。兩條機制的交會點就是 Fig 18 那個尖峰。
低溫下可參與 U 過程的聲子數（Boltzmann 因子）
n_U ∼ exp(−θ / 2T)
p.126θ：Debye 溫度（由 (28) θ = (ℏv/k_B)(6π²N/V)^{1/3} 定義）；T：絕對溫度。指數裡的 2 來自「參與 U 過程的聲子能量約為 ½k_Bθ」這個估計。原文逐字：「the number of suitable phonons of the high energy ½k_Bθ required may be expected to vary roughly as exp(−θ/2T), according to the Boltzmann factor.」
全章最常被誤述的一句：(42) 裡的 ℓ 是什麼
Umklapp 小節的最後一句話，是整個第 5 章最容易被誤記、也最值得背下來的一句。原文逐字（p.126）：
「In summary, the phonon mean free path which enters (42) is the mean free path for umklapp collisions between phonons and not for all collisions between phonons.」
拆開來說：K = ⅓Cvℓ 裡的那個 ℓ，是 umklapp 碰撞的平均自由程，不是所有聲子碰撞的平均自由程。
為什麼這一句這麼重要？因為它把前面所有的論證收束成一句可操作的規則：
- N 碰撞完全不計入 ℓ。不管它們發生得多頻繁，它們對 J 的改變是零，對熱阻的貢獻也是零。
- 於是 ℓ 這個量不是「聲子平均走多遠會撞到別的聲子」，而是「聲子平均走多遠會失去方向記憶」。這兩者在數值上可以差好幾個數量級。
- 高溫時兩者差一個 O(1) 的因子（因為相當比例的碰撞是 U），所以 p.126 允許你「without particular distinction」地估算；低溫時兩者差指數倍，絕對不能混用。
這一句話在 p.122 就已經先埋了伏筆——Table 2 的表頭註釋逐字寫著：
「[Calculated from (44), taking v = 5 × 10⁵ cm/sec as a representative sound velocity. The ℓ's obtained in this way refer to umklapp processes.]」
也就是說，Table 2 那些 40 Å、540 Å、23 Å、100 Å 的數字，全部都是 umklapp 自由程。學生常把它們讀成「聲子撞一次要走多遠」，那是錯的。
同一個結構在第 6 章會再出現一次。金屬的電阻率也有一模一樣的邏輯：小角度的正常電子—聲子散射幾乎不改變電子的漂移動量、對電阻貢獻很小；真正把電阻做出來的是 umklapp 散射。p.151 逐字：「Umklapp scattering of electrons by phonons (Chapter 5) accounts for most of the electrical resistivity of metals at low temperatures.」而 p.152 給出對應的指數 exp(−θ_U/T)，並用鉀當例子：「The experimental data (Fig. 12) for potassium have the expected exponential form with θ_U = 23 K compared with the Debye θ = 91 K.」（逐字）。同樣是指數、同樣來自 G，只是主角從聲子換成電子。
晶格熱導率（氣體動力論形式）
K = ⅓ C v ℓ
(42) p.122K：熱導率係數；C：單位體積的熱容（p.122 逐字：「C is the heat capacity per unit volume」）；v：平均粒子速度（對聲子即聲速）；ℓ：粒子在兩次碰撞之間的平均自由程——依 p.126 的結論，這裡專指 umklapp 碰撞的自由程。
Imperfections：當聲子彼此不再相撞
小節標題 Imperfections 印在 p.126 下半。開場兩句逐字：
「Geometrical effects may also be important in limiting the mean free path. We must consider scattering by crystal boundaries, the distribution of isotopic masses in natural chemical elements, chemical impurities, lattice imperfections, and amorphous structures.」
這份清單有五項，值得逐項對應到實際材料：
| 機制 | 散射中心 | 典型主宰的場景（本欄為本講解補充，非書中內容） |
|---|---|---|
| crystal boundaries（晶體邊界） | 試片表面、晶粒界面 | 極低溫的高純單晶；奈米結構 |
| distribution of isotopic masses（同位素質量分布） | 質量與平均值不同的原子核 | 天然元素的高純晶體，熱導率峰值附近 |
| chemical impurities（化學雜質） | 摻雜原子 | 摻雜半導體 |
| lattice imperfections（晶格缺陷） | 空位、差排、堆疊差錯 | 形變過的晶體、輻照後 |
| amorphous structures（非晶結構） | 沒有長程週期性 | 玻璃、SiO₂ |
這些機制與 U 過程有一個根本的差別，而 Kittel 在 p.123 早就說清楚了：
「Phonon collisions with a static imperfection or a crystal boundary will not by themselves establish thermal equilibrium, because such collisions do not change the energy of individual phonons: the frequency ω₂ of the scattered phonon is equal to the frequency ω₁ of the incident phonon.」（逐字，p.123）
靜態缺陷是彈性散射：它改變聲子的方向，但不改變聲子的能量。所以它能貢獻熱阻（因為 J 被打亂），卻不能建立局部熱平衡（因為不能在不同頻率之間重新分配能量）。要定義溫度、要寫下 dT/dx，你還是需要非彈性的聲子—聲子過程。這是一個微妙但重要的分工：U 過程負責「產生熱阻」＋「建立平衡」，缺陷散射只負責前者。
尺寸效應：熱導率變成試片幾何的函數
p.126 下半的核心段落，逐字：
「When at low temperatures the mean free path ℓ becomes comparable with the width of the test specimen, the value of ℓ is limited by the width, and the thermal conductivity becomes a function of the dimensions of the specimen. This effect was discovered by de Haas and Biermasz. The abrupt decrease in thermal conductivity of pure crystals at low temperatures is caused by the size effect.」
「At low temperatures the umklapp process becomes ineffective in limiting the thermal conductivity, and the size effect becomes dominant, as shown in Fig. 18. One would expect then that the phonon mean free path would be constant and of the order of the diameter D of the specimen, so that」——接著就是 (48)。
這裡有三件事要抓住：
第一，這是一個「材料性質不再是材料性質」的時刻。熱導率照理說是 intensive 的材料常數，與你切多大塊無關。但在低溫下 K 變成試片尺寸的函數——把同一塊晶體切細，熱導率就會下降。de Haas 與 Biermasz（1930 年代在萊頓）就是靠比較不同粗細的試片發現這件事的。書上只給人名不給年份。
第二，「abrupt decrease」這個字眼。Kittel 說的 abrupt decrease 指的是降溫方向——過了峰頂往低溫走，K 沿著 T³ 直線墜下去（對數—對數圖上斜率 3），而這個陡峭來自 C ∝ T³，不是來自 ℓ——因為 ℓ 在這裡已經凍成常數了。至於峰的右側，在 Fig 18 的溫度範圍（約 17–80 K）裡是 U 過程指數解凍造成的陡降（資料點數位化後量得的對數—對數斜率約 −4），比低溫側還陡；真正的 K ∝ 1/T 漸近段要到 T > θ 才出現，而 NaF 的 θ 在 500 K 量級（Kittel 這幾頁未給，外部值），根本不在這張圖上。
第三，這是「彈道 vs 擴散」的分界線。當 ℓ ≥ D，聲子從試片一邊飛到另一邊都不會撞到任何東西，輸運不再是擴散而是彈道（ballistic）。p.121 在定義 (41) 時特別警告過這件事：「If the energy were propagated directly through the specimen without deflection, then the expression for the thermal flux would not depend on the temperature gradient, but only on the difference in temperature between the ends of the specimen, regardless of the length of the specimen.」（逐字，p.121）。
換句話說，在尺寸效應主宰的區域，「熱導率」這個概念其實已經在崩解——(48) 之所以還能寫成 K 的形式，是因為我們硬把試片直徑 D 塞進 ℓ 的位置，讓公式在形式上維持成立。這一點在奈米尺度元件（包括 FinFET 的鰭片）裡不是學術細節，而是設計現實。
尺寸效應主宰時的熱導率
K ≅ C v D
(48) p.126K：熱導率；C：單位體積熱容；v：聲速；D：試片直徑。原文逐字：「the phonon mean free path would be constant and of the order of the diameter D of the specimen」。注意書上用的是「≅」（約等於）而非「=」，因為 ℓ ≈ D 本身就只是量級關係。
為什麼低溫是 T³ 而不是別的冪次：一步一步補推導
書上在 p.127 只用兩句話帶過「The only temperature-dependent term on the right is C … We expect the thermal conductivity to vary as T³」。這裡把 C ∝ T³ 的來源接回 p.114 的 Debye T³ 定律，並且把三個 T³ 分清楚——這是初學者最混淆的地方：低溫熱容是 T³、低溫熱導率也是 T³，但它們是同一個 T³（熱導率的 T³ 就是熱容的 T³ 直接傳過來的），而不是兩個獨立的巧合。
- 從 (32)（p.114）：低溫 Debye 熱容 C_V = (12π⁴/5)Nk_B(T/θ)³ ≡ 234 Nk_B(T/θ)³。這裡 N 是晶體中的原子數——p.113 在 (29) 下方明寫「where N is the number of atoms in the specimen」（逐字），p.114 的簡單論證也是數 3N(T/θ)³ 個被激發的模態（3N = 3×原子數）。單原子晶格時原胞數與原子數相同，但鍺這類每原胞 2 個原子的結構必須用原子數，否則 C 會差 2 倍。
- 換成單位體積的熱容（(42)(48) 要的是這個）：C = 234 (N/V) k_B (T/θ)³。注意 (42) 的 C 是 per unit volume，(32) 的 C_V 是整塊晶體的——這個換算漏掉就會差一個體積因子。
- T³ 的物理來源（p.114 的簡單論證）：低溫下只有 ℏω < k_BT 的模態被激發，在 K 空間裡佔據的體積分率是 (K_T/K_D)³ = (T/θ)³，其中 ℏvK_T = k_BT。被激發的模態數約 3N(T/θ)³，每個帶能量 k_BT，故 U ≈ 3Nk_BT(T/θ)³ ∝ T⁴，C = ∂U/∂T ∝ T³。
- 把 C 代入 (48)：K ≅ CvD = 234 (N/V) k_B (T/θ)³ · v · D。
- 整理成 K = A T³，其中 A = 234 (N/V) k_B v D / θ³。這個 A 完全由材料常數與試片幾何決定，沒有任何可調參數。
- 推論一：A ∝ D。把試片直徑加倍，低溫熱導率的係數就加倍。這是 de Haas–Biermasz 的實驗簽名，也是檢驗「真的是邊界散射嗎」的判準。
- 推論二：A ∝ v/θ³；而由 (28) θ = (ℏv/k_B)(6π²N/V)^{1/3} 可知，固定原子密度時 v ∝ θ，所以跨材料比較（固定 N/V 與 D）時實際上是 A ∝ 1/θ²。Debye 溫度越高（越硬、越輕的晶體），低溫 T³ 段的係數越小。這聽起來反直覺（鑽石不是熱導率之王嗎？），但要記得高 θ 同時把 U 過程的指數 exp(θ/2T) 往後推得很遠，使得 K 的峰值極高——T³ 段的斜率係數小，但峰值出現得晚、爬得久、爬得高。
- 推論三：實驗上只要在對數—對數座標上量低溫段的斜率，斜率 = 3 就證明尺寸效應主宰；斜率明顯小於 3 通常代表有別的散射（同位素、雜質）在插手。Fig 19 的虛線 K = 0.06T³ 就是拿來做這個判斷的參考線。
Fig 18 到底畫了什麼（NaF 的熱導率駝峰）我的判讀
Fig 18 印在 p.127 的上半，圖說逐字：
「Thermal conductivity of a highly purified crystal of sodium fluoride, after H. E. Jackson, C. T. Walker, and T. F. McNelly.」
以下是對原圖的目視描述（已放大原頁確認軸標與刻度）。
這是一張對數—對數圖。橫軸標「Temperature, K」，刻度標示為 1、2、5、10、20、50、100。縱軸標「Thermal conductivity, in W cm⁻¹ K⁻¹」，刻度標示為 1、2、5、10、20、50、100、200。
曲線形狀是一個很尖的駝峰：
- 從左下角起（T ≈ 2 K，K ≈ 1.4），資料點沿著一條近乎直線往右上方爬。在對數—對數圖上是直線，代表冪次律；目視斜率接近 3，對應 K ∝ T³ ——這就是尺寸效應段。
- 在 T ≈ 15–17 K 附近達到峰值，K 略高於 200 W cm⁻¹K⁻¹（峰頂稍微超過縱軸最上面那條 200 的刻度線）。
- 越過峰頂後急遽下墜，而且下墜段比上升段更陡；到 T ≈ 80 K 時已經掉回 K ≈ 1.4 左右。
這張圖要你看懂的三件事：
第一，峰值高得驚人。200 W cm⁻¹K⁻¹ 是一個離子晶體（氟化鈉，食鹽結構的近親），在 17 K 時的熱導率比銅的最大值（p.128 給的 100 W cm⁻¹K⁻¹）還高一倍。這就是下一段「Dielectric crystals may have thermal conductivities as high as metals」的具體證據。
第二，峰的兩側是兩套完全不同的物理。左側（低溫）是 K ≅ CvD，斜率 3 由熱容給出，ℓ 被 D 鎖死；右側（高溫）是 U 過程開始解凍，ℓ_U ∝ exp(θ/2T) 從天文數字急速塌陷，把 K 一路壓下去。峰頂就是 ℓ_U ≈ D 的那個溫度。
第三，「highly purified」這個形容詞是必要條件。若晶體不純，雜質散射會在峰的位置搶先把 ℓ 鎖住，峰會變矮變鈍。Fig 18 之所以能長成這麼尖的形狀，是因為 Jackson、Walker、McNelly 用的是高純度樣品——峰的高度，本質上是樣品純度的度量。這一點在下一張 Fig 19 會被推到極致：連同位素的差別都看得出來。
介電晶體可以比金屬還會導熱
p.127 末到 p.128 初，Kittel 列了一組數字來打破「金屬才會導熱」的直覺。逐字：
「Dielectric crystals may have thermal conductivities as high as metals. Synthetic sapphire (Al₂O₃) has one of the highest values of the conductivity: nearly 200 W cm⁻¹ K⁻¹ at 30 K. The maximum of the thermal conductivity in sapphire is greater than the maximum of 100 W cm⁻¹ K⁻¹ in copper. Metallic gallium, however, has a conductivity of 845 W cm⁻¹ K⁻¹ at 1.8 K. The electronic contribution to the thermal conductivity of metals is treated in Chapter 6.」
把這幾個數字整理成表（全部來自 p.127–128 原文）：
| 材料 | 熱導率 | 溫度 | 類別 |
|---|---|---|---|
| 合成藍寶石 Al₂O₃ | 近 200 W cm⁻¹K⁻¹ | 30 K | 介電晶體（聲子輸運） |
| 銅（最大值） | 100 W cm⁻¹K⁻¹ | — | 金屬（電子輸運） |
| 金屬鎵 | 845 W cm⁻¹K⁻¹ | 1.8 K | 金屬（電子輸運） |
幾點要注意：
第一，單位是 W cm⁻¹K⁻¹，不是 W m⁻¹K⁻¹。200 W cm⁻¹K⁻¹ = 20,000 W m⁻¹K⁻¹。室溫銅是 4.01 W cm⁻¹K⁻¹（p.116 Table 1，該表為影像、此數字為目視判讀）= 401 W m⁻¹K⁻¹。所以低溫峰值比室溫值高兩個數量級以上——這一點在跨教科書比對數字時極容易搞錯。
第二，這些都是低溫峰值，不是室溫值。藍寶石的 200 是在 30 K，鎵的 845 是在 1.8 K。它們能這麼高，正是因為在那個溫度下散射機制幾乎全部關閉（U 過程凍結、缺陷極少），ℓ 長到宏觀尺度。
第三，Kittel 特地用 however 把鎵拉出來。銅雖然輸給藍寶石，但這不代表「電子輸運不如聲子輸運」——鎵在 1.8 K 的 845 說明電子在夠低的溫度、夠純的金屬裡可以壓倒性地贏。這個伏筆在第 6 章 p.156「THERMAL CONDUCTIVITY OF METALS」與 Wiedemann–Franz 定律（p.156「Ratio of thermal to electrical conductivity」）裡才會收。
第四，句末那句「treated in Chapter 6」是第 5 章正文唯一一次明示往後交棒。（全章 p.105–129 搜尋「Chapter」只有兩處：這裡，以及習題 1(a) 回頭引第 4 章色散關係的「derived in Chapter 4」。）第 5 章從頭到尾只算聲子對熱導的貢獻，金屬裡的電子貢獻被完整地留給第 6 章。讀第 6 章時記得回來對照：同樣是 umklapp、同樣是指數壓制，只是載子換人。
同位素散射：連原子核的體重差都看得見
p.128 第二段，逐字：
「In an otherwise perfect crystal, the distribution of isotopes of the chemical elements often provides an important mechanism for phonon scattering. The random distribution of isotopic mass disturbs the periodicity of the density as seen by an elastic wave. In some substances scattering of phonons by isotopes is comparable in importance to scattering by other phonons. Results for germanium are shown in Fig. 19. Enhanced thermal conductivity has been observed also in isotopically pure silicon and diamond; the latter has device importance as a heat sink for laser sources.」
這段話的關鍵是第二句，值得逐字拆。
「The random distribution of isotopic mass disturbs the periodicity of the density as seen by an elastic wave.」
為什麼是「density（密度）」而不是「force constants（力常數）」？因為同位素的差別只在原子核的質量，化學鍵（電子結構）完全一樣。所以彈性波看到的力常數場仍然是完美週期的，被破壞的只有質量密度 ρ 的週期性。而彈性波方程 ρ ∂²u/∂t² = C ∂²u/∂x² 裡，ρ 與 C 同等重要——ρ 不週期，波就會被散射。
這是一個非常「乾淨」的散射機制：沒有化學擾動、沒有應變場、沒有電荷，純粹只是質量亂序。正因為乾淨，它成為檢驗聲子散射理論的理想對象。
「comparable in importance to scattering by other phonons」這句話的份量常被低估。它是說：在某些材料的某個溫區，單單是天然同位素的分布，就足以與所有非諧聲子—聲子散射打成平手。Fig 19 的鍺就是這個宣稱的證據。
最後一句的工程意涵：「isotopically pure silicon and diamond; the latter has device importance as a heat sink for laser sources.」——同位素純化的鑽石被拿來當雷射光源的散熱片。這是第 5 章少數幾處直接指向元件應用的地方。天然鑽石有約 1.1% 的 ¹³C，把它降到 ppm 級就能顯著提升熱導率（Kittel 此處未給，外部知識）。
Fig 19 到底畫了什麼（鍺的同位素效應）我的判讀
Fig 19 印在 p.127 的下半。圖說逐字：
「Isotope effect on thermal conduction in germanium, amounting to a factor of three at the conductivity maximum. The enriched specimen is 96 percent Ge⁷⁴, natural germanium is 20 percent Ge⁷⁰, 27 percent Ge⁷², 8 percent Ge⁷³, 37 percent Ge⁷⁴, and 8 percent Ge⁷⁶. Below 5 K the enriched specimen has K = 0.06 T³, which agrees well with Casimir's theory for thermal resistance caused by boundary scattering. (After T. H. Geballe and G. W. Hull.)」
兩個樣品的同位素組成（整理自上述圖說）：
| 樣品 | 組成 |
|---|---|
| 濃化樣品（Enriched Ge⁷⁴） | 96 % Ge⁷⁴ |
| 天然鍺（Normal Ge） | 20 % Ge⁷⁰、27 % Ge⁷²、8 % Ge⁷³、37 % Ge⁷⁴、8 % Ge⁷⁶ |
（天然鍺的五項相加 = 20+27+8+37+8 = 100 %，數字自洽。）
以下是對原圖的目視描述（已放大原頁確認）。
同樣是對數—對數圖。橫軸標「Temperature, in K」，刻度 1、2、5、10、20、50、100、200、500。縱軸標「Thermal conductivity, K in W cm⁻¹ K⁻¹」，刻度 0.1、0.2、0.5、1、2、5、10、20、50，最上面還有一格。圖上有方格狀的格線。
校註（回原頁座標核對）：縱軸最上面那一格書上印成「10」，但依刻度間距它實際應該是 100：以 PDF 文字層座標算，0.1→1 相隔 72.3 pt、1→10 相隔 71.5 pt，而該標籤又在 10 之上 70.0 pt 處，正好是再一個十年。這是書本這張圖的排版錯誤。Enriched Ge⁷⁴ 的峰值 ≈ 38 是對著下面的 20 與 50 讀出來的，不受影響。
圖上有兩條實心曲線：
- 上面那條標記 Enriched Ge⁷⁴，資料點為實心圓。峰值出現在 T ≈ 15–20 K，K ≈ 35–40 W cm⁻¹K⁻¹。
- 下面那條標記 Normal Ge，資料點為空心圓。峰值出現在大致相同的溫度，K ≈ 11 W cm⁻¹K⁻¹。
兩條曲線在峰值處相差約三倍——這正是圖說開頭說的「amounting to a factor of three at the conductivity maximum」。
另外還有一條虛線，標記 K = 0.06T³，從左下方（T ≈ 2 K，K ≈ 0.5）往右上延伸，穿過濃化樣品的低溫資料點，一直畫到超過峰值的高度才截斷。它是純 T³ 的參考線。
高溫端（T > 100 K）兩條曲線逐漸靠攏，到 200–300 K 時都掉到 K ≈ 0.8–1 W cm⁻¹K⁻¹ 附近，幾乎重合。
這張圖要你看懂的四件事：
第一，三倍差距出現在峰值，不是在高溫。高溫時 U 過程主宰，同位素散射相形之下微不足道，兩個樣品沒有差別。只有在峰值附近——U 過程剛好凍結、邊界效應還沒完全接手的那個窗口——同位素散射才是決定性的。這正是「comparable in importance to scattering by other phonons」的具體意思。
第二，兩條曲線在極低溫（T < 5 K）也會靠攏。因為那裡邊界散射主宰，同位素散射（其散射率隨頻率急遽下降）已經不重要了。所以同位素效應有一個中間溫區的窗口。
第三，虛線 K = 0.06T³ 是「純邊界散射」的理論線。濃化樣品在 5 K 以下貼著它走，代表那裡已經沒有其他散射機制了。圖說明說這「agrees well with Casimir's theory for thermal resistance caused by boundary scattering」。
第四，天然鍺並不是一種「不純的材料」。它化學上是 100 % 的鍺、晶體結構完美，但熱導率就是只有濃化樣品的三分之一。「純度」在熱輸運裡有兩層意義：化學純度與同位素純度。這是本段最反直覺、也最容易在口試被追問的一點。
數值檢核：從 K = 0.06 T³ 反推試片直徑補推導
Fig 19 圖說給了一個具體的數字：濃化鍺在 5 K 以下 K = 0.06 T³（單位為 W cm⁻¹K⁻¹）。既然書上說這符合「邊界散射」的理論，那我們就用 (48) K ≅ CvD 反推 D，看看得到的是不是一個合理的試片尺寸。這是驗證一條公式是否真的在說它宣稱的事情的標準做法，強烈建議自己動手算一次。
外部輸入的聲明：以下用到兩個不在本段原文裡的數字——鍺的 Debye 溫度 θ = 374 K（來自 p.116 Table 1，該表在 PDF 裡是無文字層的影像，此數字為算繪後目視判讀）；鍺的原子數密度 N/V = 4.42×10²² cm⁻³（標準值，Kittel 這幾頁未給）。全程用 CGS。
- 先算低溫熱容（單位體積）。由 (32)：C = 234 (N/V) k_B (T/θ)³ = 234 × 4.42×10²² × 1.381×10⁻¹⁶ / (374)³ × T³。
- 分子：234 × 4.42×10²² = 1.034×10²⁵；× 1.381×10⁻¹⁶ = 1.428×10⁹。分母：374³ = 5.23×10⁷。
- 得 C = 27.3 T³ erg cm⁻³ K⁻¹。（檢查量綱：erg cm⁻³ K⁻¹ 正是單位體積熱容，對。）
- 再算聲速。由 (28) θ = (ℏv/k_B)(6π²N/V)^{1/3} 反解 v。先算 (6π² × 4.42×10²²)^{1/3} = (2.617×10²⁴)^{1/3} = 1.378×10⁸ cm⁻¹（這就是 K_D）。
- v = θ k_B / (ℏ K_D) = 374 × 1.381×10⁻¹⁶ / (1.055×10⁻²⁷ × 1.378×10⁸) = 5.17×10⁻¹⁴ / 1.454×10⁻¹⁹ = 3.55×10⁵ cm/s。（合理：鍺的聲速在 3–5 ×10⁵ cm/s 量級。）
- 把實驗值換成 CGS：K = 0.06 T³ W cm⁻¹K⁻¹，而 1 W = 10⁷ erg/s，故 K = 6.0×10⁵ T³ erg s⁻¹cm⁻¹K⁻¹。
- 用 (48) K ≅ CvD 反解：D = K/(Cv) = 6.0×10⁵ T³ / (27.3 T³ × 3.55×10⁵) = 6.0×10⁵ / 9.70×10⁶ = 6.2×10⁻² cm ≈ 0.6 mm。
- 若堅持保留 (42) 的 ⅓ 係數（K = ⅓CvD），則 D = 3 × 0.62 = 1.9 mm。
- 判讀：0.6–1.9 mm 正是 1950 年代量低溫熱導率所用的棒狀單晶試片的典型橫向尺寸。反推出來的「有效自由程」落在毫米量級，也就是宏觀尺度——這正是「聲子在 5 K 以下的鍺裡可以一路飛到試片表面才被散射」這句話的量化版本。(48) 確實在說它宣稱的事。
- 順帶得到一個可證偽的預測：把同一根晶體縱向剖成一半粗，0.06 這個係數應該降到約 0.03。這就是 de Haas–Biermasz 效應的定量形式，也是判別「是否真的是邊界散射」最直接的實驗。
- 注意 T³ 在兩邊完全抵消，所以反推出的 D 與溫度無關——這是一致性檢查通過的標誌。如果算出來的 D 隨溫度跑，就代表 ℓ 其實不是常數，尺寸效應的假設不成立。
補充：質量亂序為什麼會散射聲子（Kittel 未寫的機制）延伸知識
本 block 全部是 Kittel 這幾頁沒有的內容，提供給想知道「為什麼同位素效應集中在峰值附近」的讀者。書上只寫了「The random distribution of isotopic mass disturbs the periodicity of the density as seen by an elastic wave.」一句，沒有給任何定量結構。
把一個質量為 M + δM 的同位素放進由質量 M 的原子構成的完美晶格，等效於在運動方程裡加一個局部的微擾項 δM ∂²u/∂t² = −δM ω² u。關鍵在於這個微擾正比於 ω²，也就是頻率越高、擾動越強。
用費米黃金律算散射率，矩陣元平方帶來 ω⁴，再乘上末態密度（三維 Debye 的 D(ω) ∝ ω²）之後，點缺陷散射率的標準結果是
1/τ_iso ∝ V₀ Γ ω⁴ / v³
其中 Γ = Σᵢ fᵢ (1 − Mᵢ/M̄)² 是同位素的質量變異數因子（fᵢ 為第 i 種同位素的豐度、M̄ 為平均質量），V₀ 為每個原子的體積。這種 ω⁴ 依賴就是Rayleigh 散射的形式——與天空為什麼是藍色的同一條律則。（「Rayleigh」一詞在 Kittel 第 8 版中只出現於表面波與光學章節，與此處無關；這條公式是外部知識。）
ω⁴ 這個冪次解釋了 Fig 19 的一切：
- 高溫：U 過程的散射率隨溫度快速增長，把同位素散射蓋過去 → 兩條曲線重合。
- 極低溫：只有低頻聲子被激發，ω⁴ 使同位素散射率趨近於零 → 邊界散射主宰，兩條曲線又重合。
- 峰值附近：高頻聲子已被激發（ω⁴ 大），但 U 過程還沒解凍 → 同位素散射獨佔鰲頭，兩條曲線分得最開，差三倍。
拿天然鍺的組成代進 Γ：以 20 % Ge⁷⁰、27 % Ge⁷²、8 % Ge⁷³、37 % Ge⁷⁴、8 % Ge⁷⁶ 算出 M̄ = 0.20×70 + 0.27×72 + 0.08×73 + 0.37×74 + 0.08×76 = 14.0 + 19.44 + 5.84 + 27.38 + 6.08 = 72.74（與鍺的標準原子量 72.6 相符，組成自洽）。Γ ≈ Σ fᵢ(1 − Mᵢ/72.74)² ≈ 5.8×10⁻⁴。96 % 濃化的 Ge⁷⁴ 對應的 Γ 約小一個數量級以上——這個數量級的差別，就換來峰值三倍的熱導率。
PROBLEMS 的定位：五題各自補回正文省略的哪一塊我的判讀
習題從 p.128 中段開始（標題 Problems），到 p.129 結束，共五題。第 5 章沒有 SUMMARY 節，所以這五題實質上就是本章的收束。
它們不是練習題，是正文因篇幅而省略的五塊物理：
| 題號 | 標題 | 補回的物理 | 難度 |
|---|---|---|---|
| 1 | Singularity in density of states | Van Hove 奇異點的兩個具體案例（p.119 只給了 Fig 14 的定性圖） | 中等，純微積分 |
| 2 | Rms thermal dilation of crystal cell | 熱漲落的絕對大小——把 Debye 理論接到熔化判據 | 容易但需查外部數值 |
| 3 | Zero point lattice displacement and strain | 零點運動的振幅；一維晶體為何沒有長程序 | 中等偏難 |
| 4 | Heat capacity of layer lattice | 維度如何決定熱容冪次（C ∝ T^d） | 容易，(a) 是 T³ 推導的直接類比 |
| 5 | Grüneisen constant | 把零點能、非諧性、熱膨脹三者綁成一條式子 | 書上加星號：「somewhat difficult」 |
書上明確標註難度的只有第 5 題。p.129 的頁尾註腳逐字寫著：「*This problem is somewhat difficult.」題號前的星號就是指這個。
第 5 題也是唯一被正文回頭引用的一題。p.121 講熱膨脹時逐字寫道：「The expansion coefficient vanishes as T → 0, as we expect from Problem 5.」——正文把一個結論的證明外包給了習題。這代表第 5 題不是選作，是本章論證鏈的一環。
以下逐題把書上省略的代數全部補出來。每一題都先講「這題到底在問什麼」，再列步驟。
Problem 1(a)：一維單原子鏈的態密度奇異點補推導
題目逐字（p.128）：「From the dispersion relation derived in Chapter 4 for a monatomic linear lattice of N atoms with nearest-neighbor interactions, show that the density of modes is D(ω) = (2N/π)·1/(ω_m² − ω²)^{1/2}, where ω_m is the maximum frequency.」
這題在問什麼：p.119 的 Fig 14 只用定性的方式畫出「真實晶體的 D(ω) 會出現不連續」，並說這些是 Van Hove 奇異點，出現在群速度為零的臨界點。這一題要你親手把第一個具體案例算出來——一維鏈在區邊界的群速度為零（第 4 章 (8) 已證），於是 D(ω) 在 ω → ω_m 時發散。
- 取第 4 章的 (9)：ω = (4C/M)^{1/2}|sin ½Ka|。定義 ω_m ≡ (4C/M)^{1/2}，即最大頻率（在 K = ±π/a 取得）。於是在第一區的右半 0 ≤ K ≤ π/a 上 ω = ω_m sin(½Ka)。
- 一維的 K 空間態密度：長度 L = Na 的鏈，允許的 K 值間距為 2π/L，故單位 K 區間有 L/2π 個模態。只計 K > 0 的半區間（負 K 由 |sin| 的對稱性給出相同的 ω，兩邊各半、總數仍為 N），對應的模態密度為 L/π = Na/π。
- 驗證模態總數：∫₀^{π/a}(Na/π)dK = (Na/π)(π/a) = N ✓。一維 N 個原子確實有 N 個模態。
- 求 dω/dK：dω/dK = ω_m·(a/2)cos(½Ka)。
- 把 cos 換成 ω：cos(½Ka) = [1 − sin²(½Ka)]^{1/2} = [1 − (ω/ω_m)²]^{1/2} = (ω_m² − ω²)^{1/2}/ω_m。
- 代回：dω/dK = ω_m·(a/2)·(ω_m² − ω²)^{1/2}/ω_m = (a/2)(ω_m² − ω²)^{1/2}。（注意 ω_m 剛好消掉，這是這題算得乾淨的關鍵。）
- 態密度 D(ω) = (模態在 K 的密度)×(dK/dω) = (Na/π)·2/[a(ω_m² − ω²)^{1/2}]。
- a 消去，得 D(ω) = (2N/π)·1/(ω_m² − ω²)^{1/2} ✓ 與題目要求一致。
- 歸一化檢查（強烈建議做）：∫₀^{ω_m}D(ω)dω = (2N/π)[arcsin(ω/ω_m)]₀^{ω_m} = (2N/π)(π/2) = N ✓。
- 物理解讀：當 ω → ω_m 時分母 → 0，D(ω) 發散（雖然積分仍收斂，是可積的平方根發散）。發散的來源正是 dω/dK → 0——區邊界的群速度為零，代表在一個很窄的頻率區間內塞進了很多個 K 值。這就是 Van Hove 奇異點的原型，與 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)」（逐字）完全對上。
- 為什麼 Debye 模型看不到這個：Debye 取 ω = vK（線性、群速度處處等於 v），永遠不會有 dω/dK = 0，所以 Debye 的 D(ω) ∝ ω² 一路平滑到截止頻率再被硬切掉。Fig 14(a) 與 (b) 的差別就在這裡。
Problem 1(b)：光學支的不連續態密度補推導
題目逐字（p.128）：「Suppose that an optical phonon branch has the form ω(K) = ω₀ − AK², near K = 0 in three dimensions. Show that D(ω) = (L/2π)³(2π/A^{3/2})(ω₀ − ω)^{1/2} for ω < ω₀ and D(ω) = 0 for ω > ω₀. Here the density of modes is discontinuous.」
這題在問什麼：(a) 處理的是聲學支在區邊界的奇異點；(b) 處理的是光學支在 K = 0 的極大值。光學支在 K = 0 有 dω/dK = 0（頂點），所以也是一個臨界點，但奇異性的型態完全不同——這次是跳躍不連續，不是發散。
- 從 ω = ω₀ − AK² 反解 K：K = [(ω₀ − ω)/A]^{1/2}。注意這要求 ω ≤ ω₀，因為 K² ≥ 0。
- 三維 K 空間的態密度：體積 V = L³ 的晶體，每個 K 態佔 (2π/L)³，故單位 K 體積有 (L/2π)³ 個態。半徑 K 到 K + dK 的球殼體積是 4πK²dK。
- 於是 D(ω)dω = (L/2π)³ · 4πK² · |dK|，即 D(ω) = (L/2π)³ · 4πK² · |dK/dω|。
- 求 |dK/dω|：由 ω = ω₀ − AK² 得 dω/dK = −2AK，故 |dK/dω| = 1/(2AK)。
- 代入：D(ω) = (L/2π)³ · 4πK²/(2AK) = (L/2π)³ · 2πK/A。
- 把 K 換回 ω：D(ω) = (L/2π)³ · (2π/A) · [(ω₀ − ω)/A]^{1/2} = (L/2π)³ · (2π/A^{3/2}) · (ω₀ − ω)^{1/2} ✓ 與題目一致。
- ω > ω₀ 時沒有任何模態（K 會是虛數），故 D(ω) = 0。
- 奇異性的型態：在 ω → ω₀⁻ 時 D → 0（像 √ 一樣趨近於零），在 ω > ω₀ 恆為 0。所以 D(ω) 本身是連續的（都趨於 0），不連續的是它的導數 dD/dω——(ω₀−ω)^{1/2} 的斜率在 ω₀ 是無限大。題目說的「the density of modes is discontinuous」指的是這個斜率的突變（相對於 (a) 的函數值發散，這是較弱的奇異性）。
- 與 (a) 的對照：(a) 是一維、臨界點在區邊界、D(ω) 以 (ω_m² − ω²)^{−1/2} 發散；(b) 是三維、臨界點在 K = 0 的極大、D(ω) 以 (ω₀ − ω)^{+1/2} 趨零。冪次的正負號由維度決定：三維球殼的 4πK² 因子把發散壓下去了。這兩題合起來就是 Van Hove 分類的縮影——同樣是 ∇_K ω = 0，維度與臨界點型態（極大、極小、鞍點）決定了奇異性長什麼樣。
- 量綱檢查：(L/2π)³ 的量綱是 [長度³]，A 的量綱是 [頻率·長度²]（因為 AK² 要是頻率），故 A^{3/2} 是 [頻率^{3/2}·長度³]，而 (ω₀−ω)^{1/2} 是 [頻率^{1/2}]。合起來：長度³ × 頻率^{1/2}/(頻率^{3/2}長度³) = 1/頻率 ✓ 正是 D(ω) 該有的量綱（每單位頻率的態數）。
Problem 2：原胞的均方熱膨脹（把 Debye 接到熔化）補推導
題目逐字（p.128）：「(a) Estimate for 300 K the root mean square thermal dilation ΔV/V for a primitive cell of sodium. Take the bulk modulus as 7 × 10¹⁰ erg cm⁻³. Note that the Debye temperature 158 K is less than 300 K, so that the thermal energy is of the order of k_BT. (b) Use this result to estimate the root mean square thermal fluctuation Δa/a of the lattice parameter.」
這題在問什麼：前面整章算的都是平均值（平均能量、熱容）。這一題問的是漲落的大小——一個原胞的體積在 300 K 時會抖動多少？答案大到會讓人吃驚，而且直接連到固體為什麼會熔化。
題目給的提示怎麼用：「the Debye temperature 158 K is less than 300 K, so that the thermal energy is of the order of k_BT」——這是在說 300 K 已經是古典區（T > θ），所以可以直接用能量均分定理，每個自由度分到 ½k_BT，不必碰 Planck 分布。
外部輸入聲明：計算需要鈉的原胞體積，這要用晶格常數。Kittel 第 1 章的晶體結構表（書本 p.20–21）在本 PDF 裡是無文字層的影像，抽不到數值，本段原文也沒給。以下採用標準值 a ≈ 4.23 Å（bcc），並明確標示為外部輸入。
- 建立能量表達式：均勻膨脹 Δ ≡ ΔV/V 所儲存的彈性能密度是 ½BΔ²（這是體積模數 B 的定義）。對一個體積為 V_c 的原胞，能量 U = ½BΔ²V_c。
- 套能量均分：Δ 是一個平方型自由度，在古典區其平均能量為 ½k_BT。故 ½B⟨Δ²⟩V_c = ½k_BT。
- 解出 ⟨Δ²⟩ = k_BT/(B V_c)。這條式子本身就值得記住：熱漲落的大小 = 熱能 ÷ (彈性剛度 × 體積)。體積越小、材料越軟、溫度越高，漲落越大。
- 算原胞體積：鈉是 bcc，慣用立方晶胞含 2 個原子，故原胞（primitive cell）體積 V_c = a³/2。以 a = 4.23 Å = 4.23×10⁻⁸ cm：a³ = 7.57×10⁻²³ cm³，V_c = 3.78×10⁻²³ cm³。
- 算熱能：k_BT = 1.381×10⁻¹⁶ × 300 = 4.14×10⁻¹⁴ erg。
- 算分母：B V_c = 7×10¹⁰ × 3.78×10⁻²³ = 2.65×10⁻¹² erg。（量綱檢查：erg cm⁻³ × cm³ = erg ✓）
- 代入：⟨Δ²⟩ = 4.14×10⁻¹⁴ / 2.65×10⁻¹² = 1.56×10⁻²。
- (a) 的答案：⟨(ΔV/V)²⟩^{1/2} = (1.56×10⁻²)^{1/2} = 0.125，約 12.5 %。
- (b)：對小形變 V = a³ ⟹ ΔV/V = 3Δa/a，故 Δa/a = 0.125/3 ≈ 0.042，約 4.2 %。
- 判讀一：4 % 的晶格常數漲落非常大。鈉的熔點是 371 K，而 300 K 已經是熔點的 81 %。Lindemann 熔化判據（Kittel 這幾頁未提，外部知識）說晶體在原子均方位移達到最近鄰距離的 10 % 左右時熔化——4 % 在 300 K、10 % 在熔點附近，量級完全吻合。這一題真正的意義是：它讓你看到「軟金屬在室溫下離熔化並不遠」這件事可以從 Debye 理論估出來。
- 判讀二：⟨Δ²⟩ ∝ 1/(B V_c) 解釋了為什麼輕且軟的元素（鹼金屬、固態氫、固態氦）漲落大、熔點低，而重且硬的元素（鎢、鑽石）漲落小、熔點高。鑽石的 B（約 4.4×10¹² erg cm⁻³，外部標準值）比題目給鈉的 7×10¹⁰ erg cm⁻³ 大約 60 倍（近兩個數量級），θ = 2230 K（p.116 Table 1，目視判讀），室溫下連古典區都還沒進入。
- 常見錯誤：把 V_c 取成慣用立方晶胞體積 a³（而不是 bcc 的 a³/2），會得到 ⟨Δ²⟩ 少一半、Δ_rms ≈ 8.8 %。題目寫的是 primitive cell，所以必須用 a³/2。
Problem 3(a)：絕對零度的零點位移補推導
題目逐字（p.128）：「In the Debye approximation, show that the mean square displacement of an atom at absolute zero is ⟨R²⟩ = 3ℏω_D²/8π²ρv³, where v is the velocity of sound. Start from the result (4.29) summed over the independent lattice modes: ⟨R²⟩ = (ℏ/2ρV)Σω⁻¹. We have included a factor of ½ to go from mean square amplitude to mean square displacement.」
這題在問什麼：即使在 T = 0，原子也不是靜止的——每個模態都還有 ½ℏω 的零點能。這一題要你算出這個量子抖動有多大，答案會是「幾個百分點的鍵長」。它同時是 Problem 2 的 T = 0 對照組：Problem 2 算的是熱漲落，這題算的是量子漲落。
先把題目給的起點釐清（書上這裡寫得含糊，很多人卡在這）：第 4 章的 (29) 是 u₀² = 4(n + ½)ℏ/ρVω，那是振幅平方。對駐波 u = u₀ cos Kx cos ωt，均方位移要同時對時間與空間取平均：⟨cos²ωt⟩ = ½、⟨cos²Kx⟩ = ½，合起來 ⟨u²⟩ = ¼u₀²。在 T = 0（n = 0）時 u₀² = 2ℏ/(ρVω)，故每個模態貢獻 ¼ × 2ℏ/(ρVω) = ℏ/(2ρVω)——這正好就是題目給的 (ℏ/2ρV)Σω⁻¹ 裡的每一項。書上只說「included a factor of ½」，沒說另一個 ½ 來自空間平均。
- 起點：⟨R²⟩ = (ℏ/2ρV) Σ_K ω_K⁻¹，其中 ρ 是質量密度、V 是晶體體積，求和跑遍所有獨立的晶格模態（含三個極化）。
- 把求和換成積分：Σ_K → ∫₀^{ω_D} D(ω)dω，其中 Debye 態密度（三個極化一起算）由 p.112 的 (22) 乘 3 得 D(ω) = 3Vω²/(2π²v³)。
- 代入：Σω⁻¹ = ∫₀^{ω_D} [3Vω²/(2π²v³)]·ω⁻¹ dω = (3V/2π²v³)∫₀^{ω_D} ω dω。
- 積分：∫₀^{ω_D} ω dω = ω_D²/2。故 Σω⁻¹ = 3Vω_D²/(4π²v³)。
- 回代：⟨R²⟩ = (ℏ/2ρV) × 3Vω_D²/(4π²v³) = 3ℏω_D²/(8π²ρv³) ✓ 與題目一致。V 完全消掉，說明這是一個 intensive 的量（與晶體大小無關），符合物理直覺。
- 注意積分為什麼收斂：被積函數是 ω²·ω⁻¹ = ω，在 ω → 0 完全沒問題。三維的 ω² 態密度救了這個積分——這一點在 (b) 的一維情形就不成立了，那正是 (b) 要考的。
- 數值示範（鑽石）：取 θ = 2230 K（p.116 Table 1，目視判讀）、ρ = 3.51 g cm⁻³、N/V = 1.76×10²³ cm⁻³（外部標準值）。先算 K_D = (6π²N/V)^{1/3} = 2.18×10⁸ cm⁻¹；由 (28) 得 v = θk_B/(ℏK_D) = 2230×1.381×10⁻¹⁶/(1.055×10⁻²⁷×2.18×10⁸) = 1.34×10⁶ cm/s；ω_D = vK_D = 2.92×10¹⁴ rad/s。
- 代入：分子 3×1.055×10⁻²⁷×(2.92×10¹⁴)² = 3×1.055×10⁻²⁷×8.53×10²⁸ = 2.70×10²。分母 8π²ρv³ = 78.96×3.51×(1.34×10⁶)³ = 277×2.41×10¹⁸ = 6.67×10²⁰。
- 得 ⟨R²⟩ = 4.05×10⁻¹⁹ cm²，⟨R²⟩^{1/2} = 6.4×10⁻¹⁰ cm = 0.064 Å。鑽石的最近鄰距離是 1.54 Å，所以零點振幅約為鍵長的 4 %。
- 判讀：即使在絕對零度、即使是已知最硬的晶體，原子仍在以鍵長 4 % 的幅度抖動。這不是熱運動，是量子力學要求的。對照 Lindemann 的 10 % 熔化判據：這說明像固態氦那種 θ 極低、質量極輕的系統，零點運動大到足以在 T = 0 就把晶體「熔掉」——那正是量子固體與零點壓的起點（Kittel 此處未涉及）。
- 與 Problem 2 的對照：Problem 2 的鈉在 300 K 得到 Δa/a ≈ 4 %（熱）；這題的鑽石在 0 K 得到 ⟨R²⟩^{1/2}/d ≈ 4 %（量子）。兩個 4 % 的來源完全不同，但都逼近同一個熔化尺度——這不是巧合，而是兩者都被同一個彈性剛度尺度所控制。
Problem 3(b)：一維晶體的發散，與有限的均方應變補推導
題目逐字（p.128）：「Show that Σω⁻¹ and ⟨R²⟩ diverge for a one-dimensional lattice, but that the mean square strain is finite. Consider ⟨(∂R/∂x)²⟩ = ½ΣK²u₀² as the mean square strain, and show that it is equal to ℏω_D²L/4MNv³ for a line of N atoms each of mass M, counting longitudinal modes only. The divergence of R² is not significant for any physical measurement.」
這題在問什麼：這是本章最有深度的一題，雖然它藏在 Problem 3 的 (b) 裡。它在問：一維晶體到底存不存在？答案是「原子的絕對位置沒有意義（發散），但原子之間的相對距離完全良好（收斂）」。這正是 Peierls–Landau 定理的雛形——一維與二維系統沒有真正的長程位置序。
書上最後那句「The divergence of R² is not significant for any physical measurement」是關鍵提示：可觀測的是應變（相對位移），不是絕對位移。
- 先看發散。一維 Debye：長度 L、N 個原子、只算縱模。態密度 D(ω) = (L/π)(dK/dω) = L/(πv)，是一個常數（與 ω 無關——這是一維的特徵，對照三維的 ω²、二維的 ω）。
- Σω⁻¹ = ∫ D(ω)ω⁻¹dω = (L/πv)∫_{ω_min}^{ω_D} dω/ω = (L/πv) ln(ω_D/ω_min)。
- 最低頻模態由最長波長 λ = 2L 決定，ω_min = πv/L；而 ω_D = πvN/L（由 K_D = πN/L）。故 ln(ω_D/ω_min) = ln N。
- 於是 Σω⁻¹ = (L/πv)ln N，隨 N 對數發散；⟨R²⟩ = (ℏ/2ρL)Σω⁻¹ = ℏ ln N/(2πρv)（ρ 為線密度 NM/L，固定），同樣對數發散 ✓。
- 物理意義：發散來自 ω → 0 的長波聲學模。這些模態是整條鏈的緩慢起伏，能量極低所以振幅極大。鏈越長，能容納的長波越長，原子相對於「絕對座標原點」的位置就越不確定。這就是為什麼一維晶體沒有長程序。
- 再看應變。應變只在乎鄰居之間的相對位移，對整體平移完全不敏感——數學上就是多乘一個 K（因為 ∂/∂x → iK），把 ω → 0 的貢獻壓掉。
- 取題目給的定義 ⟨(∂R/∂x)²⟩ = ½Σ K²u₀²，並用 T = 0 的 (4.29)：u₀² = 2ℏ/(ρLω)。一維時 ρL = NM（總質量）。
- 得 ⟨(∂R/∂x)²⟩ = ½ Σ K² × 2ℏ/(NMω) = (ℏ/NM) Σ K²/ω。以 ω = vK 代入：= (ℏ/NMv) Σ K。
- 把求和換成積分：Σ K → ∫₀^{K_D}(L/π)K dK = (L/π)(K_D²/2) = LK_D²/(2π) = Lω_D²/(2πv²)。
- 代回得 ⟨(∂R/∂x)²⟩ = ℏLω_D²/(2πNMv³)。積分在 K → 0 完全收斂（被積函數 ∝ K），均方應變是有限的 ✓——題目要的定性結論成立。
- 與書上印的答案比對（校註）：題目寫的是 ℏω_D²L/(4MNv³)，也就是係數 1/4；用上述標準一維模態計數（密度 L/π 於 0 < K < K_D）得到的是 1/(2π)。兩者結構完全相同（同樣的 ℏ、ω_D²、L、M、N、v³ 與冪次），差的只是一個 2/π ≈ 0.64 的純數因子，來源是模態計數慣例與題目所給的 ½（若按駐波的時間＋空間平均，該係數應為 ¼）之間的不一致。這種 O(1) 係數的出入在 Kittel 的「show that」題裡不罕見；本題的物理重點在「發散 vs 收斂」，不在這個因子。 本講解據標準計數給出 1/(2π)，並據實標註與書上不符，不代為湊數。
- 把結論一般化（這是這題真正要你帶走的東西）：均方位移 ⟨R²⟩ ∝ ∫ D(ω)/ω dω。三維 D ∝ ω² → 積分 ∝ ∫ω dω 收斂；二維 D ∝ ω → 積分 ∝ ∫dω 收斂（但有限溫度下改為 ∫dω·T/ω 會對數發散）；一維 D = const → 積分 ∝ ∫dω/ω 對數發散。維度越低，長波漲落越致命。這條規則會在第 16 章（相變）與凝態物理的 Mermin–Wagner 定理裡再度出現。
Problem 4：層狀晶格的熱容——維度決定冪次補推導
題目逐字（p.128–129）：「(a) Consider a dielectric crystal made up of layers of atoms, with rigid coupling between layers so that the motion of the atoms is restricted to the plane of the layer. Show that the phonon heat capacity in the Debye approximation in the low temperature limit is proportional to T². (b) Suppose instead, as in many layer structures, that adjacent layers are very weakly bound to each other. What form would you expect the phonon heat capacity to approach at extremely low temperatures?」
這題在問什麼：Debye T³ 定律的「3」是從哪來的？答案是空間維度。這一題用二維把這件事逼出來：C ∝ T^d。(b) 則是一個很漂亮的陷阱——它問的是「當維度本身隨溫度改變時會怎樣」。
- (a) 二維 Debye。面積 A 的層，K 空間中每個態佔 (2π)²/A，半徑 K 到 K+dK 的環帶面積是 2πK dK，故態數 = (A/4π²)·2πK dK = (A/2π)K dK。
- 以 ω = vK 換元：D(ω)dω = (A/2π)(ω/v²)dω。原子運動限制在層平面內，故每個 K 有 2 個極化（兩個面內方向），D(ω) = (A/π)(ω/v²)。
- 記住這條規律：d 維 Debye 的態密度 D(ω) ∝ ω^{d−1}。一維常數、二維 ∝ ω、三維 ∝ ω²。
- 熱能 U = ∫₀^{ω_D} D(ω)·ℏω/(e^{ℏω/k_BT} − 1) dω = (Aℏ/πv²)∫₀^{ω_D} ω²/(e^{ℏω/k_BT} − 1) dω。
- 低溫極限：令 x = ℏω/k_BT，上限 → ∞。dω = (k_BT/ℏ)dx，ω² = (k_BT/ℏ)²x²。
- U = (Aℏ/πv²)(k_BT/ℏ)³ ∫₀^∞ x²dx/(e^x − 1) = (A k_B³/πv²ℏ²) T³ × 2ζ(3) ≈ (A k_B³/πv²ℏ²) T³ × 2.404。
- （用到 ∫₀^∞ x²/(e^x−1)dx = 2!·ζ(3) = 2×1.202 = 2.404。）
- C = ∂U/∂T ∝ T² ✓。
- 一般結果：d 維、線性色散下 U ∝ T^{d+1}、C ∝ T^d。積分 ∫x^d/(e^x−1)dx = d!·ζ(d+1) 永遠是一個純數。所以「T³」這個 3 從頭到尾就是空間維度，與材料無關。
- (b) 弱耦合層：關鍵在於「極低溫」這四個字。層間有一個很小的力常數，它給出一支跨層的色散，頻寬很小；設其最大頻率為 ω_⊥（對應能量尺度 ℏω_⊥，或溫度尺度 θ_⊥ = ℏω_⊥/k_B，θ_⊥ ≪ θ_∥）。
- 中間溫區（θ_⊥ ≪ T ≪ θ_∥）：所有跨層模態都已被完全激發（k_BT ≫ ℏω_⊥），跨層方向對熱容不再有「可再激發的自由度」，系統表現為一疊互相獨立的二維層 → C ∝ T²。
- 極低溫（T ≪ θ_⊥）：現在 k_BT < ℏω_⊥，連跨層方向的長波模態都開始被凍結／被解析出來。此時被激發的聲子波長遠大於層間距，它們看不見層狀結構，看到的是一塊（各向異性的）三維彈性連續體 → C 回到 ∝ T³。
- 所以答案是 T³。這題的教訓是：維度不是材料的固有屬性，而是「熱聲子波長 vs 結構特徵長度」比較出來的結果。溫度是決定你看到幾維的那把尺。
- 實例（Kittel 此處未提，外部知識）：石墨正是這個行為的教科書案例——中溫區熱容接近 T²，到數 K 以下轉為 T³。這也是為什麼討論石墨烯／二維材料的熱性質時，永遠要先問「在什麼溫度」。
Problem 5 為什麼是全章最好的整合練習
第 5 題在題號前印了一個星號，頁尾註腳逐字寫著：「*This problem is somewhat difficult.」——這是第 5 章唯一一道被作者標記難度的題目。
題目逐字（p.129）：
「5. Grüneisen constant. (a) Show that the free energy of a phonon mode of frequency ω is k_BT ln [2 sinh (ℏω/2k_BT)]. It is necessary to retain the zero-point energy ½ℏω to obtain this result. (b) If Δ is the fractional volume change, then the free energy of the crystal may be written as
F(Δ, T) = ½BΔ² + k_BT Σ ln [2 sinh (ℏω_K/2k_BT)] ,
where B is the bulk modulus. Assume that the volume dependence of ω_K is δω/ω = −γΔ, where γ is known as the Grüneisen constant. If γ is taken as independent of the mode K, show that F is a minimum with respect to Δ when BΔ = γΣ½ℏω coth (ℏω/2k_BT), and show that this may be written in terms of the thermal energy density as Δ = γU(T)/B. (c) Show that on the Debye model γ = −∂ ln θ/∂ ln V. Note: Many approximations are involved in this theory: the result (a) is valid only if ω is independent of temperature; γ may be quite different for different modes.」
為什麼說它是整合練習？因為它同時用到本章四塊互不相干的內容：
| 用到的東西 | 出處 | 在這題的角色 |
|---|---|---|
| 零點能 ½ℏω | 第 4 章 (27) | (a) 少了它算不出 sinh |
| Planck 分布／熱能 | 本章 (9)(26) | (b) 把 coth 認出是熱能 |
| 非諧性 | 本章 p.119–120 | γ ≠ 0 的唯一來源 |
| Debye 溫度 θ 的定義 | 本章 (28) | (c) 把 γ 接到可量的量 |
而且正文欠了它一筆債。p.121 逐字：「The expansion coefficient vanishes as T → 0, as we expect from Problem 5.」——正文在講熱膨脹時，把「為什麼熱膨脹係數在 T → 0 會消失」這個證明直接外包給本題。所以這題不做，第 5 章的熱膨脹段落就有一個洞。
以下把 (a)(b)(c) 的每一步代數補完。
Problem 5(a)：單一聲子模態的自由能補推導
要證：F = k_BT ln[2 sinh(ℏω/2k_BT)]。
這一步的重點不在技巧，在於「為什麼一定要留著零點能」。題目自己就把答案講出來了：「It is necessary to retain the zero-point energy ½ℏω to obtain this result.」下面會看到，零點能提供的正是 sinh 需要的那個 e^{+x/2}。
- 設 β ≡ 1/k_BT。單一模態的能階為 ε_n = (n + ½)ℏω，n = 0, 1, 2, …（第 4 章 (27)）。
- 配分函數 Z = Σ_{n=0}^∞ e^{−βε_n} = e^{−βℏω/2} Σ_{n=0}^∞ (e^{−βℏω})^n。
- 右邊是公比 e^{−βℏω} < 1 的幾何級數，和為 1/(1 − e^{−βℏω})。故 Z = e^{−βℏω/2}/(1 − e^{−βℏω})。
- 自由能 F = −k_BT ln Z = −k_BT[−βℏω/2 − ln(1 − e^{−βℏω})] = k_BT[βℏω/2 + ln(1 − e^{−βℏω})]。
- 把第一項塞進對數：βℏω/2 = ln e^{βℏω/2}，故 F = k_BT ln[e^{βℏω/2}(1 − e^{−βℏω})]。
- 展開括號：e^{βℏω/2} − e^{βℏω/2}·e^{−βℏω} = e^{βℏω/2} − e^{−βℏω/2}。
- 用 sinh 的定義 2 sinh y = e^y − e^{−y}（取 y = βℏω/2）：F = k_BT ln[2 sinh(ℏω/2k_BT)] ✓
- 零點能扮演的角色：若把 ε_n = nℏω（不含 ½ℏω），則 Z = 1/(1 − e^{−βℏω})、F = k_BT ln(1 − e^{−βℏω})。這個表達式在 T → 0 時趨於 0，而正確的 F 在 T → 0 時趨於 ½ℏω（把 sinh 展開：T → 0 時 2sinh(y) → e^y，F → k_BT·(ℏω/2k_BT) = ½ℏω ✓）。沒有零點能，(b) 裡對 Δ 微分就得不到 coth，只會得到 Planck 分布，最後算出來的 T = 0 膨脹會是零——那是錯的。
- 高溫極限檢查：T ≫ ℏω/k_B 時 sinh(y) ≈ y，F ≈ k_BT ln(ℏω/k_BT)，這正是古典諧振子的自由能 ✓。
- 順帶：由 F 可以直接回推平均能量 U = ∂(βF)/∂β = (ℏω/2)coth(ℏω/2k_BT)，這條式子在 (b) 會直接用到。
Problem 5(b)：對 Δ 取極小，得到 Δ = γU(T)/B補推導
要證：F 對 Δ 極小時 BΔ = γ Σ ½ℏω coth(ℏω/2k_BT)，且可寫成 Δ = γU(T)/B。
物理設定：晶體的自由能有兩項。第一項 ½BΔ² 是彈性能——它永遠反對膨脹，把 Δ 拉回 0。第二項是聲子的自由能——膨脹會讓頻率下降（δω/ω = −γΔ），而頻率下降會讓聲子自由能下降，所以它支持膨脹。熱膨脹就是這兩項的拔河結果。（注意 F 的兩項都是單位體積的量，故稍後 U 是熱能密度。）
- F(Δ, T) = ½BΔ² + k_BT Σ_K ln[2 sinh(ℏω_K/2k_BT)]，其中 ω_K = ω_K(Δ)。
- 把 δω/ω = −γΔ 寫成顯式的體積依賴：ω_K(Δ) = ω_K(0)(1 − γΔ)，故 ∂ω_K/∂Δ = −γ ω_K（假設 γ 與模態 K 無關，這是題目允許的）。
- 對 Δ 微分。第一項給 BΔ。第二項用鏈鎖律：∂/∂Δ {k_BT ln[2 sinh(ℏω_K/2k_BT)]} = k_BT · [cosh/sinh](ℏω_K/2k_BT) · (ℏ/2k_BT) · ∂ω_K/∂Δ。
- 整理：= (ℏ/2)coth(ℏω_K/2k_BT) · (−γω_K) = −γ·½ℏω_K coth(ℏω_K/2k_BT)。（k_BT 在這一步完全消掉，這是題目設計得漂亮的地方。）
- 所以 ∂F/∂Δ = BΔ − γ Σ_K ½ℏω_K coth(ℏω_K/2k_BT)。
- 令 ∂F/∂Δ = 0：BΔ = γ Σ_K ½ℏω_K coth(ℏω_K/2k_BT) ✓ 這正是題目要的第一個結果。
- 辨認右邊的求和。單一諧振子的平均能量（含零點能）是 ⟨ε⟩ = (⟨n⟩ + ½)ℏω，其中 ⟨n⟩ = 1/(e^x − 1)、x ≡ ℏω/k_BT。
- 代數恆等式：1/(e^x − 1) + ½ = (2 + e^x − 1)/[2(e^x − 1)] = (e^x + 1)/[2(e^x − 1)]。分子分母同除 e^{x/2}：= (e^{x/2} + e^{−x/2})/[2(e^{x/2} − e^{−x/2})] = ½ coth(x/2)。
- 故 ⟨ε⟩ = ℏω·½coth(ℏω/2k_BT) = ½ℏω coth(ℏω/2k_BT)——右邊求和裡的每一項，就是一個模態的平均能量（含零點能）。
- 於是 Σ_K ½ℏω_K coth(ℏω_K/2k_BT) = U(T)，晶體的熱能密度（含零點能）。
- 得 BΔ = γU(T)，即 Δ = γU(T)/B ✓ 第二個結果完成。
- 推論一（書上沒寫，但這是這題的重點）：對溫度微分得體積熱膨脹係數 α_V ≡ dΔ/dT = (γ/B)(dU/dT) = γC/B，其中 C 是單位體積熱容。這就是 Grüneisen 關係式。它說：熱膨脹係數與熱容成正比，比例常數 γ/B。
- 推論二（p.121 欠的那筆債）：低溫下 C ∝ T³ → 0，故 α_V → 0。這就是 p.121 那句「The expansion coefficient vanishes as T → 0, as we expect from Problem 5.」的證明。更強的敘述是：α 與 C 有相同的溫度依賴，所以量熱膨脹與量熱容在低溫是等價的實驗。
- 推論三：T = 0 時 U(0) = Σ½ℏω ≠ 0（零點能），故 Δ(0) = γU(0)/B ≠ 0。晶體在絕對零度就已經比「力學平衡位置」漲了一點——這是純量子效應，稱為零點膨脹。輕元素（固態氦、氫、鋰）的零點膨脹可達百分之幾。
- 題目自己的但書（逐字）：「Many approximations are involved in this theory: the result (a) is valid only if ω is independent of temperature; γ may be quite different for different modes.」後半句很重要——真實晶體裡橫聲學模的 γ 常常是負的，這就是矽、鍺、以及鑽石在低溫下出現負熱膨脹的原因（此現象 Kittel 此處未提，屬外部知識）。整體的 γ 是各模態以熱容加權的平均，低溫時負 γ 的橫模權重上升，總 α 就翻負。
Problem 5(c)：γ = −∂lnθ/∂lnV，以及它為什麼等於零就沒有熱膨脹補推導
要證：在 Debye 模型下 γ = −∂ ln θ/∂ ln V。
這一小題是把 γ 從「一個抽象的假設」變成「一個可以查表、可以量的材料常數」。θ 可以從低溫熱容量出來（Problem 的 T³ 係數），而 ∂θ/∂V 可以從高壓實驗量出來，於是 γ 變成實驗可及的量。
- 從定義出發：δω/ω = −γΔ，而 Δ = δV/V，故對每一個模態 γ = −(δω/ω)/(δV/V) = −∂ ln ω/∂ ln V。
- Debye 模型的關鍵性質：整個聲子譜由單一參數刻畫。所有模態都是 ω = vK，而允許的 K 值由邊界條件固定（K 的集合隨體積等比縮放），截止頻率由 (23) ω_D³ = 6π²v³N/V 給出，並由 (27) 與 (28) 對應到 θ = ℏω_D/k_B。
- 若 γ 與模態無關（題目的假設），則所有 ω 以同一比例隨體積變化，整個頻譜等比縮放，於是 ω_D 也以同樣的比例縮放，θ = ℏω_D/k_B 亦然。
- 因此 −∂ ln ω/∂ ln V = −∂ ln ω_D/∂ ln V = −∂ ln θ/∂ ln V = γ ✓
- 把 (28) 代進去看清楚結構：θ = (ℏv/k_B)(6π²N/V)^{1/3}，取對數：ln θ = ln(ℏ/k_B) + ln v + ⅓ln(6π²N) − ⅓ln V。
- 微分：γ = −∂lnθ/∂lnV = ⅓ − ∂ln v/∂ln V。
- 這條分解式非常有用。它把 γ 拆成兩塊：⅓ 是純幾何（體積變大、K 值變小、頻率跟著變小）；−∂ln v/∂ln V 是物理（聲速隨體積怎麼變，也就是力常數隨體積怎麼變）。
- 嚴格諧振晶體的檢查（這一步值得自己算一次）：諧振近似下力常數 C 是固定的，聲速 v = a(C/M)^{1/2}（一維鏈的長波極限），故 ln v = ln a + 常數，而 a ∝ V^{1/3} ⟹ ∂ln v/∂ln V = ⅓。
- 代入：γ = ⅓ − ⅓ = 0。於是 Δ = γU/B = 0——嚴格諧振晶體沒有熱膨脹，完全符合 p.119 那條列點：「There is no thermal expansion.」（逐字，p.119 的諧振理論後果清單）。
- 結論：γ ≠ 0 的唯一來源，是力常數本身隨體積改變——而這就是非諧性的定義。Grüneisen 常數是「非諧性有多強」的無因次度量。
- 典型數值（Kittel 此處未給，外部參考）：多數固體 γ 落在 1–3；矽在室溫約 0.45、鍺約 0.75，且兩者在低溫轉為負值。γ 的量級是 1 而不是 100 或 0.01，這件事本身就說明：非諧項雖然是「微擾」，但它與諧振項的比值在原子尺度並不算小——只是在位移上小，在體積依賴上不小。
- 與本段開頭的呼應：整個 Umklapp 小節的存在理由，也是非諧性（沒有非諧項就沒有三聲子過程、沒有 U 過程、熱導率無限大）。Problem 5 量化的 γ，與 p.123 那句「With anharmonic lattice interactions, there is a coupling between different phonons which limits the value of the mean free path」（逐字）講的是同一個非諧性的兩個面向：一個管熱膨脹，一個管熱阻。
常見誤解與自我檢查清單我的判讀
把這一段最容易錯的地方集中列出來。每一條後面都附上原文出處，可以回頭對。
誤解一：「U 過程不守恆能量。」
錯。Fig 17 圖說逐字：「For all processes, N or U, energy must be conserved, so that ω₁ + ω₂ = ω₃.」p.126 又重複一次：「The umklapp process must conserve energy, just as for the normal process.」U 過程折的是波向量，能量分毫不差。被破壞的是聲子氣的晶體動量（少了 ℏG），而那筆帳由整個晶體接走。
誤解二：「(42) 裡的 ℓ 是聲子的平均自由程。」
不完整。p.126 的最後一句逐字：「the phonon mean free path which enters (42) is the mean free path for umklapp collisions between phonons and not for all collisions between phonons.」是 umklapp 碰撞的平均自由程。高溫下兩者只差 O(1)，低溫下差指數倍。p.122 Table 2 的註釋也再說一次：「The ℓ's obtained in this way refer to umklapp processes.」
誤解三：「低溫下 U 過程的聲子數隨 exp(−θ/T) 變化。」
錯了一個 2。p.126 逐字：「may be expected to vary roughly as exp(−θ/2T)」。分母的 2 來自「參與 U 過程的聲子能量約為 ½k_Bθ」，而 ½ 又來自「波向量需要約 ½G」。這條 ½ 的鏈子不能斷。
誤解四：「碰撞越多，熱阻越大。」
錯。Fig 16c（p.124）畫的就是反例：只有 N 過程的晶體，不管撞多少次，聲子流都不衰減，熱阻為零。熱阻的判準是「總動量 J 有沒有被改變」，不是「有沒有碰撞」。
誤解五：「(48) K ≅ CvD 裡的 T³ 來自 ℓ。」
錯，來自 C。p.127 逐字：「The only temperature-dependent term on the right is C, the heat capacity, which varies as T³ at low temperatures.」在尺寸效應區，ℓ 被鎖成常數 D，完全不隨溫度變。
誤解六：「天然鍺是不純的材料。」
它化學上 100% 是鍺、晶格完美，但熱導率峰值只有 96% Ge⁷⁴ 濃化樣品的三分之一。同位素質量亂序本身就是一種散射中心，即使晶體「完美」。p.128 逐字：「In an otherwise perfect crystal, the distribution of isotopes of the chemical elements often provides an important mechanism for phonon scattering.」
誤解七：「熱導率是材料常數。」
在低溫不是。p.126 逐字：「the thermal conductivity becomes a function of the dimensions of the specimen. This effect was discovered by de Haas and Biermasz.」把同一塊晶體切細，K 會下降。
誤解八：「金屬一定比絕緣體會導熱。」
p.127–128 用三個數字反駁：合成藍寶石 30 K 時近 200 W cm⁻¹K⁻¹，高於銅的最大值 100 W cm⁻¹K⁻¹（但金屬鎵在 1.8 K 是 845）。注意這些都是低溫峰值，且單位是 W cm⁻¹K⁻¹ 不是 W m⁻¹K⁻¹。
自我檢查題（能全部答出就算讀懂這一段）：
- 為什麼 N 過程做不到「把 K_x 為負的聲子流變成正的」，U 過程做得到？
- 在 Fig 17(b) 裡，哪一支箭頭是虛線，為什麼？
- exp(−θ/2T) 裡的 2 是哪來的？
- 為什麼 Fig 18 的曲線在低溫段的對數—對數斜率是 3？
- 為什麼鍺的同位素效應在峰值最明顯，而在高溫與極低溫都看不見？
- Problem 5 的 γ 如果等於 0，會發生什麼事？為什麼？
延伸：這一段與奈米元件熱阻抽取的接點延伸知識
本 block 全部是 Kittel 沒有的內容，寫給在做元件熱模型（例如 FinFET 的熱容／熱阻萃取）的讀者，說明這五頁的物理在什麼地方會變成工程約束。
第一，尺寸效應不是低溫的專利。p.126 說「When at low temperatures the mean free path ℓ becomes comparable with the width of the test specimen」——注意成立的條件是 ℓ ≳ D，溫度只是達成這個條件的手段之一。在室溫的矽裡，主導熱輸運的聲子的平均自由程分布相當寬（有相當比例的熱量由 μm 量級自由程的聲子攜帶），而現代 FinFET 的鰭片寬度是數奈米。於是 ℓ ≫ D 在室溫就成立，de Haas–Biermasz 當年在 4 K 才看到的現象，在今天的元件裡是常態。實務後果是：塊材矽的 K ≈ 1.5 W cm⁻¹K⁻¹ 不能直接拿來用在鰭片上，有效熱導率可能掉到十分之一以下。
第二，(48) K ≅ CvD 正是「邊界主宰」時的估算式。它告訴你有效熱導率線性正比於特徵尺寸。這條線性關係是元件熱模型裡「尺寸相依熱導率」修正的最原始版本。
第三，(42) 的分解 K = ⅓Cvℓ 對應到熱等效電路的兩個元件。C（單位體積熱容）決定熱容 C_th、K 決定熱阻 R_th。值得注意的是：C 幾乎不受尺寸效應影響（它是模態計數的結果，只要結構還是三維的就仍是 Debye 那一套），但 K 受尺寸效應嚴重壓制。所以在奈米元件裡，R_th 被強烈放大、C_th 幾乎不變，熱時間常數 τ = R_th C_th 隨之拉長——這是自熱（self-heating）分析的核心。把塊材的 K 套進奈米結構，是元件熱模型最常見的系統性錯誤。
第四，U 過程的 exp(θ/2T) 在低溫量測裡會咬人。若你在低溫（例如 77 K 或 4 K）做熱量測來抽取參數，聲子自由程可能已經長到與整個結構相當，「熱導率」這個局域量的概念本身就在崩解（見 p.121 對彈道輸運的警告）。此時應改用彈道／準彈道的輸運描述，而不是 Fourier 定律。
第五，同位素工程是真的工程。p.128 已經點名了：「Enhanced thermal conductivity has been observed also in isotopically pure silicon and diamond; the latter has device importance as a heat sink for laser sources.」同位素純化矽的熱導率高於天然矽，這在高功率密度散熱上是一條實際的路線。
回到本章：Kittel 這五頁講的是 1950 年代的塊材低溫實驗，但它建立的三個判準至今未變——(1) 只有改變總動量的散射才產生熱阻；(2) 溫度透過 Boltzmann 因子開關 U 過程；(3) 當 ℓ 撞上幾何尺度，熱導率就不再是材料常數。做元件熱分析時，這三條每一條都會用到。