這一段在解決什麼問題我的判讀
看一眼第 5 章的目錄（書本 p.105）就知道這一段站在哪裡：PHONON HEAT CAPACITY 從 p.107 開始，一路做 Planck 分布、簡正模式計數、一維與三維的態密度、Debye 模型、Debye T³ 定律、Einstein 模型，到 p.117 的 General result for D(ω) 為止。然後 p.119 才出現 ANHARMONIC CRYSTAL INTERACTIONS，接著是 Thermal expansion (120)、THERMAL CONDUCTIVITY (121)、Thermal resistivity of phonon gas (123)、Umklapp processes (125)、Imperfections (126)。
前半章（p.107–119 上半）處理的全部是熱容，而熱容這個量有一個很方便的性質：它只需要知道「有哪些模式、每個模式的頻率是多少、每個模式上有幾個聲子」。模式之間有沒有交互作用，對熱容的計算完全不影響——因為熱容只問總能量對溫度的導數，不問能量怎麼在模式之間流動。所以整個前半章可以躲在諧近似後面，一路躲到 p.119。
問題是，熱膨脹與熱傳導這兩件事，躲不過去。熱膨脹問的是「平均原子間距是多少」，熱傳導問的是「能量怎麼從熱端流到冷端」。前者需要位能不對稱，後者需要模式之間會互相散射——這兩件事在嚴格的諧晶體裡恰好都是零（不是很小，是嚴格為零、無窮大）。因此 Kittel 在這裡停下來，先做一件在教科書裡很少見的事：把自己前面用的近似列成清單，一條一條說它許諾了什麼，然後宣布這些許諾在真實晶體裡都不成立。
這份清單就是這一段的骨架。第 5 章後半所有的論證——熱膨脹的 (38)–(40)、熱阻的 (41)–(46)、Umklapp 的 (47)——全部都是在修補這五條的其中一條。讀這一段的正確方式不是「看過就算」，而是把五條記牢，之後每讀到一個新結果就回來問：這是在修哪一條？
（順帶一提，p.119 的上半頁還屬於前一節：那裡在講 D(ω) 的一般公式與 Van Hove 奇異點，以及 Fig. 13、Fig. 14。本段從同一頁的 ANHARMONIC CRYSTAL INTERACTIONS 標題開始。）
諧近似到底是什麼：原文在哪裡下刀
原文的定義句是：「The theory of lattice vibrations discussed thus far has been limited in the potential energy to terms quadratic in the interatomic displacements. This is the harmonic theory; among its consequences are:」（p.119）
請注意這句話把刀下在位能上，不是動能上。動能永遠是 ½Mv² 的二次式，沒有近似可言；所有的近似都發生在位能的 Taylor 展開。把一對原子（或整個晶體）的位能對平衡構型展開：
U = U₀ + Σᵢ (∂U/∂xᵢ)₀ xᵢ + ½ Σᵢⱼ (∂²U/∂xᵢ∂xⱼ)₀ xᵢxⱼ + (1/6) Σᵢⱼₖ (∂³U/∂xᵢ∂xⱼ∂xₖ)₀ xᵢxⱼxₖ + …
常數項 U₀ 只是能量零點，隨便丟。一次項恆為零，因為展開點就是平衡構型，一階導數按定義是零——這一點很重要，它保證了二次項是「最低的非平庸階」。二次項保留下來，就是第 4 章一路在用的力常數 C；把運動方程寫出來就是 (2)、(18) 那些線性耦合方程（單原子基是第 4 章 p.91 的 (2)，雙原子基是 p.97 的 (18)；p.95 的 (16a) 不是運動方程，而是把色散關係 (7) 推廣到 p 個近鄰面的結果），解出來就是 ω(K)。三次項以上全部丟掉，這就是諧近似（harmonic approximation）。
關鍵在於，「丟掉三次以上」不是一個小小的技術簡化，它把問題的性質整個換掉了。只保留二次項時，運動方程是線性的；線性方程可以做正規座標變換，把 3N 個耦合的自由度變成 3N 個完全獨立的簡諧振子。這個「完全獨立」就是後面五條後果的共同源頭：聲子氣體是一個理想氣體，粒子之間不碰撞、不交換能量、彼此看不見。
非諧項（anharmonic terms）做的事，用一句話說就是：把聲子從理想氣體變成有交互作用的氣體。三次項給聲子—聲子的三體交互作用，四次項給四體交互作用。第 5 章後半全部的物理，都是這個交互作用的後果。
（校註／分界：本節屬於原文的只有開頭那句定義。位能的 Taylor 展開、「一次項在平衡構型為零、三次以上可略」這一層意思，出現在第 4 章 p.91；而正規座標變換、「聲子理想氣體」的說法，以及「非諧＝把理想氣體變成有交互作用的氣體」這個框架，都是本站為了串起後面五條後果而補的教學鋪陳，不在 p.119–120 的原文裡。）
原文列的五個後果（逐字）
p.119 的條列共五條，原文逐字如下：
「• Two lattice waves do not interact; a single wave does not decay or change form with time.
- There is no thermal expansion.
- Adiabatic and isothermal elastic constants are equal.
- The elastic constants are independent of pressure and temperature.
- The heat capacity becomes constant at high temperatures T > θ.」
（校註：最後一條的 θ 在原書這裡沒有下標 D，是目視原頁確認的；第 5 章別處提到 Debye 溫度時多寫 θ。這裡的 θ 就是 Debye 溫度。）
這五條的排列其實是有邏輯的：第一條講動力學（波會不會互相作用），第二條講結構（平均間距會不會變），第三、四條講彈性響應（彈性常數的行為），第五條講熱力學（熱容）。也就是說，Kittel 是刻意挑了晶體物理四個不同面向各一個代表，來說明諧近似的影響有多廣。
下面五個小節逐條講透：每條先講「為什麼諧近似必然導致它」，再講「真實晶體怎麼破功、破得多嚴重」。
| 後果 | 諧近似下為什麼成立 | 真實晶體如何破功 | 本書哪裡處理 |
|---|---|---|---|
| 波不互相作用、不衰減 | 運動方程線性，簡正模式彼此正交 | 三聲子過程 ω₃ = ω₁ + ω₂ | p.120 下半、p.121–126 熱傳導 |
| 沒有熱膨脹 | 位能對 x = 0 對稱，⟨x⟩ ≡ 0 | 位能的 −gx³ 不對稱項 | p.120 Thermal expansion |
| 絕熱＝等溫彈性常數 | ω 不隨應變改變，聲子熵與應變無關 | ω 隨體積改變（Grüneisen γ ≠ 0） | Problem 5 (p.129) |
| 彈性常數與 P、T 無關 | 二階導數是常數，三階導數為零 | 三階以上非諧係數不為零 | Problem 5 (p.129) |
| 高溫熱容為常數 | 均分律：每個二次自由度 ½k_BT | 熱膨脹與非諧修正使 C 續升 | 本章未展開 |
後果一：兩個晶格波互不作用，單一波不衰減不變形我的判讀
為什麼諧近似必然給出這條。 諧位能讓運動方程變成線性齊次方程組。線性方程的解空間滿足疊加原理：如果 u⁽¹⁾(t) 與 u⁽²⁾(t) 各自是解，那麼 u⁽¹⁾ + u⁽²⁾ 也是解，而且是原封不動的疊加——兩個波穿過彼此之後，各自的振幅、頻率、波形完全不變，就像兩束光在真空中交錯而過。更正式地說，做完正規座標變換後，Hamiltonian 變成 H = Σ_K [P_K²/2M + ½Mω_K²Q_K²]，每一項只含一個 K，沒有任何 Q_K Q_K' 的交叉項。每個模式的能量因此是各自的守恆量。
「a single wave does not decay or change form with time」這半句其實包含兩件不同的事。不衰減：因為能量無處可去——這個模式的能量不能流到別的模式。不變形：一個單一 K 的行波本來就是 Hamiltonian 的本徵模式，時間演化只是乘一個相位因子 e^{−iω t}，空間波形完全不動。如果是一個波包（多個 K 的疊加），諧晶體裡它還是會變形，但那是色散（不同 K 走不同速度）造成的，不是非諧造成的——這兩件事不要混淆。
真實晶體怎麼破功。 就在同一段的下一頁，原文立刻給出反例：兩個聲子交互作用產生第三個聲子。一旦模式之間可以交換能量，一個被激發的模式就有有限壽命 τ，波在傳播中會衰減。
為什麼這一條最要緊。 因為它是整個熱傳導理論的存在理由。如果聲子真的互不碰撞，熱能從熱端出發會彈道式直線衝到冷端，完全不受阻礙——熱導率是無限大。原文在 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）也就是說，熱導率這個量能被定義出來，本身就預設了非諧性存在。把 K = ⅓Cvℓ 裡的平均自由程 ℓ 寫下來的那一刻，你已經承認了諧近似失效。
後果二：沒有熱膨脹我的判讀
為什麼諧近似必然給出這條，而且是嚴格為零。 取純諧位能 U(x) = cx²。Boltzmann 權重是 exp(−βcx²)，這是一個偶函數。要算的平均位移是
⟨x⟩ = ∫x e^{−βcx²}dx / ∫e^{−βcx²}dx。
分子的被積函數是「奇函數 × 偶函數 = 奇函數」，在對稱區間 (−∞, ∞) 上積分嚴格等於零。注意這個結論裡完全沒有出現 T：不管溫度多高、振幅多大，⟨x⟩ 都是零。這一點常被誤解，值得強調：諧晶體不是「膨脹得很少」，是「完全不膨脹」。
用圖像講：諧位能是一條左右對稱的拋物線。原子在裡面來回振盪，溫度升高只是讓它盪得更遠——但因為左右兩側一樣陡，它往右跑多遠就往左跑多遠，停留時間的加權中心永遠在正中央。溫度改變的是振幅，不是中心。
真實晶體怎麼破功。 真實的原子間位能（第 3 章的 Lennard-Jones 那一類）在排斥側陡、吸引側緩：把兩個原子壓近很費力，拉開相對容易。拋物線的對稱性被打破，加權中心就會往「緩」的那一側（也就是拉開的方向）偏移，而且溫度越高偏得越多。這就是熱膨脹。式 (38) 的 −gx³ 正是把這個不對稱最低階地寫下來。
這一條與後面幾條的關係。 「沒有熱膨脹」看起來像一個孤立的結論，其實它和第三條（絕熱＝等溫）是同一件事的兩種說法，這在下一節會講清楚。它也直接決定了一個實驗上很煩人的事實：C_p 與 C_V 的差 C_p − C_V = 9α²BVT 在諧晶體中恆為零（因為 α = 0），所以前半章計算 C_V 而實驗量 C_p 的那個落差，本質上也是非諧效應。
後果三：絕熱與等溫彈性常數相等我的判讀
先把兩個彈性常數的定義分清楚。 等溫彈性常數 C^T 是在「應變時讓晶體與熱庫保持同溫」的條件下量的：C^T = (∂²F/∂ε²)_T，F 是 Helmholtz 自由能。絕熱彈性常數 C^S 是在「應變快到來不及交換熱量、熵不變」的條件下量的：C^S = (∂²U/∂ε²)_S。超音波實驗量到的是絕熱值（聲波振盪太快，壓縮區來不及把熱傳給稀疏區），靜態壓縮實驗量到的是等溫值——所以這兩個值不相等會是一個實際會踩到的問題。
為什麼諧近似讓它們相等。 兩者的差別來自「應變會不會改變熵」。用 Maxwell 關係：
(∂S/∂ε)_T = −(∂σ/∂T)_ε
右邊是熱應力：固定應變時，升溫會不會產生額外的應力。在諧晶體中，力常數是一個與構型無關的常數，所以每個模式的頻率 ω_K 不依賴於應變。聲子的熵 S = Σ_K s(ℏω_K/k_BT) 只透過 ω_K 依賴於應變——既然 ω_K 不動，(∂S/∂ε)_T = 0，熱應力為零，於是等溫與絕熱兩條路徑給出同一個二階導數。
注意這和後果二是同一件事。 上式右邊為零就是「固定應變下升溫不產生應力」，等價於「固定應力下升溫不產生應變」，也就是 α = 0。所以「絕熱＝等溫」與「沒有熱膨脹」在熱力學上是互為表裡的兩句話，只要其中一個非零，另一個必然非零。標準的定量關係是
B_S / B_T = C_p / C_V = 1 + β_V γT，
其中 γ 是 Grüneisen 常數（本書 Problem 5，p.129），β_V 是體膨脹係數。（符號約定：本檔一律以 α 表示線膨脹係數、以 β_V = 3α 表示體膨脹係數；前面 C_p − C_V = 9α²BVT 裡的係數 9，正是因為那裡的 α 是線膨脹係數。）這條式子漂亮地說明了：兩個彈性常數的相對差正比於 β_VγT，而 β_V 與 γ 都是非諧的度量。
破功有多嚴重。 對室溫下的一般金屬，β_VγT 大約在百分之幾的量級——不大，但超音波量測的精度遠遠好過百分之幾，所以這個差是每天都要處理的實務問題。在接近熔點時它會長到 10% 以上。
後果四：彈性常數與壓力、溫度無關我的判讀
為什麼諧近似必然給出這條。 在諧近似中，位能被截斷成 U = cx²，於是有效彈簧常數就是 d²U/dx² = 2c，是一個寫死的常數。它不可能依賴溫度，因為位能裡根本沒有溫度這個變數；它也不可能依賴壓力，因為要讓彈性常數隨壓力改變，你需要 d³U/dx³ ≠ 0——壓力把平衡位置推到新的 x₀，新位置的二階導數才會不同；但諧位能的三階導數恆為零，所以不管把 x₀ 推到哪裡，二階導數都還是 2c。
這句話值得再說一遍，因為它是理解「為什麼三次項這麼關鍵」的最短路徑：彈性常數對壓力的導數，本質上就是位能的三階導數。同理，彈性常數對溫度的導數需要三階或四階導數（前者透過熱膨脹改變平衡位置，後者透過振幅改變有效力常數）。諧近似把這兩個導數都設成零，所以這兩個依賴性一起消失。
真實晶體破得很厲害。 實驗上體模數對壓力的無因次導數 dB/dP 對絕大多數固體落在 4 到 6 之間——這不是一個小修正，這表示壓力壓到 B 的量級時彈性常數會翻幾倍。而 dB/dP 與 Grüneisen 常數 γ 之間有直接的關聯（兩者都由三階非諧項主導），這也是為什麼 Problem 5 要你證明 γ = −∂ln θ/∂ln V：θ 隨體積改變，就是彈性常數隨體積改變的另一種寫法（因為 θ ∝ v_sound ∝ (C/ρ)^{1/2}）。
溫度方面，典型固體的彈性常數從 0 K 升到室溫會下降幾個百分比到十幾個百分比，而且這個下降在低溫時趨於平坦、在高溫時趨於線性——這個形狀本身就是聲子佔據數 ⟨n⟩ 的形狀，因為軟化的來源正是「被熱激發的聲子透過非諧項把有效力常數壓低」。
這一條在前面幾章被默默用了很多次。 第 3 章的彈性常數 C₁₁、C₁₂、C₄₄ 全部當成常數處理；第 4 章的聲速 v = (C/ρ)^{1/2} 因而也是常數；第 5 章 Debye 模型的整個架構（p.112–114）建立在「有一個與溫度無關的聲速 v」上。一旦承認這一條破功，Debye 溫度 θ 本身就是溫度的函數——這正是實驗上 θ(T) 曲線不是水平線的原因。
後果五：高溫熱容為常數（T > θ）我的判讀
為什麼諧近似必然給出這條。 這是能量均分律（equipartition）的直接結果，而均分律本身的成立條件就是「Hamiltonian 對該座標是二次式」。古典統計力學的定理說：Hamiltonian 中每一個以平方形式出現的自由度，在熱平衡下平均帶有 ½k_BT 的能量。一個簡諧振子有兩個這樣的自由度（動能的 p²、位能的 q²），所以平均能量是 k_BT。3N 個模式給出 U = 3Nk_BT，熱容 C = dU/dT = 3Nk_B，與溫度無關——這就是 Dulong–Petit 值。
請特別注意條件：均分律要求平方形式。如果位能是 U = cx² − gx³ − fx⁴，x 的貢獻就不再是 ½k_BT，均分律立刻失效。也就是說，第五條不是「諧近似恰好推出的一個結論」，而是「諧近似就是均分律的適用前提」。
「T > θ」這個條件在管什麼。 低溫時熱容偏離常數的原因是量子的（高頻模式凍結，p.107–114 的整個 Planck 分布與 Debye T³ 定律都在講這件事），與非諧無關。所以 Kittel 特別加上 T > θ，把量子效應排除掉，剩下的偏離才能歸咎於非諧。這是一個很值得學的寫法：條件是用來隔離變因的。
真實晶體怎麼破功。 實驗上在 T > θ 之後 C_p 並不會停在 3Nk_B，而是繼續緩慢上升，到熔點附近可以超出 Dulong–Petit 值 10–20%。這個超出有兩個來源，必須分開：
(i) C_p 與 C_V 的差。 C_p − C_V = 9α²BVT，這一項在高溫隨 T 線性成長。它本身就完全是非諧的產物（α = 0 時它為零），但它不是「熱容的非諧修正」，而是「量錯了對象」——實驗方便量 C_p，理論方便算 C_V。
(ii) C_V 本身的非諧修正。 把 (38) 那樣的位能做古典統計，能量的最低階非諧修正正比於 T²，於是 C_V 得到一個正比於 T 的修正項。這才是真正意義上的非諧熱容。
實務上把這兩者分開，需要獨立量到 α 與 B，這也是為什麼「高溫熱容的非諧修正」在實驗上一直是個難題。
五條其實是同一件事：諧位能的三個數學特徵我的判讀
把五條後果背下來不難，但那樣會覺得它們是五個不相干的事實。其實它們全部來自諧位能 U = cx² 的三個數學特徵，而且每一條都可以追到其中一個（或兩個）：
特徵一：可嚴格對角化，模式完全解耦。 線性運動方程 ⟹ 正規座標 ⟹ H = Σ_K H_K。這是後果一的全部來源，也是「聲子理想氣體」這個說法的精確意思。
特徵二：Boltzmann 權重是偶函數（宇稱對稱）。 e^{−βcx²} 在 x → −x 下不變，所以所有奇階矩都是零：⟨x⟩ = ⟨x³⟩ = ⟨x⁵⟩ = … = 0。這是後果二的全部來源。順帶一提，偶階矩不為零：⟨x²⟩ = k_BT/2c，這就是為什麼諧晶體雖然不膨脹，卻仍然有隨溫度成長的均方位移（Debye–Waller 因子、X 光繞射強度隨溫度衰減，都由 ⟨x²⟩ 決定，見本章 Problem 2 與 Problem 3）。「⟨x²⟩ ≠ 0 但 ⟨x⟩ = 0」這一句話，是諧近似最精準的一句自畫像。
特徵三：頻率與振幅無關、與體積無關。 古典地說，ẍ = −(2c/M)x 的解的頻率 ω = (2c/M)^{1/2} 完全不含振幅；量子地說，能階 (n + ½)ℏω 是等間距的，間距與 n 無關。既然 ω 與振幅無關，它就與溫度無關（溫度只透過振幅影響系統）；既然力常數是常數，ω 也與體積無關。這一條同時是後果三、後果四、後果五的來源：
- ω 與應變無關 ⟹ 聲子熵與應變無關 ⟹ 絕熱＝等溫（後果三）
- 力常數是常數 ⟹ 彈性常數與 P、T 無關（後果四）
- 能階等間距 ⟹ 高溫時配分函數變成古典的 ⟹ 均分律 ⟹ C = 3Nk_B（後果五）
所以真正該記的不是五條，而是三句話：模式解耦、宇稱對稱、頻率與振幅／體積無關。任何一個非諧項一進來，這三句話全部同時失效——這就是為什麼原文可以一口氣說「In real crystals none of these consequences is satisfied accurately」。
真實晶體：原文的轉折句，以及它用字的分寸
p.120 開頭的轉折是：「In real crystals none of these consequences is satisfied accurately. The deviations may be attributed to the neglect of anharmonic (higher than quadratic) terms in the interatomic displacements. We discuss some of the simpler aspects of anharmonic effects.」
這段話有三個用字的分寸值得注意。
第一，「is satisfied accurately」不是「is satisfied」。 Kittel 說的是「沒有一條精確成立」，不是「沒有一條成立」。這個分寸很重要，因為諧近似其實是非常好的近似——Debye 模型能把比熱算到幾個百分比以內，靠的就是它。非諧項是修正項，不是主角。但這一段要講的物理（熱膨脹、熱傳導）恰好是那些「諧項給零、修正項才是全部」的量，所以在這些量上，修正項就是主角。這是物理學裡很常見的情況：微擾理論的一階項為零時，二階項就不是修正而是答案。
第二，「anharmonic (higher than quadratic) terms」把定義寫死了。 非諧 = 高於二次。這排除了一個常見的混淆：「非線性」在不同領域有不同意思，這裡專指位能 Taylor 展開的三次及以上項。
第三，「some of the simpler aspects」是誠實的預告。 這一章對非諧的處理確實是最簡化的版本：熱膨脹用一個一維的古典振子，熱傳導用動力學理論的 K = ⅓Cvℓ。真正的非諧理論（微擾展開、聲子自能、Green 函數）完全不在這本書的範圍。讀者應該知道自己拿到的是什麼等級的工具。
三聲子過程的頻率條件
ω₃ = ω₁ + ω₂
p.120ω₁、ω₂ 為參與碰撞的兩個入射聲子的角頻率；ω₃ 為碰撞後產生的第三個聲子的角頻率。乘上 ℏ 之後，這條式子就是 ℏω₃ = ℏω₁ + ℏω₂，也就是能量守恆。
三聲子過程的物理圖像：移動中的三維光柵
原文用一段話把三聲子交互作用的機制講完，這段話值得逐字讀：
「The physics of the phonon interaction can be stated simply: the presence of one phonon causes a periodic elastic strain which (through the anharmonic interaction) modulates in space and time the elastic constant of the crystal. A second phonon perceives the modulation of the elastic constant and thereupon is scattered to produce a third phonon, just as from a moving three-dimensional grating.」（p.120）
把它拆成四步來讀：
第一步：一個聲子造成週期性彈性應變。 這一步是定義，不需要非諧性。聲子就是一個行波位移場，位移場的空間導數就是應變場，它在空間上以波長 2π/K₁ 週期變化、在時間上以 ω₁ 振盪。
第二步：透過非諧交互作用，這個應變調變了晶體的彈性常數。 這一步必須有非諧性。在嚴格的諧晶體裡，彈性常數是常數——不管你把晶體拉伸多少，它的彈性常數都不變（這正是後果四）。一旦位能有三次項，彈性常數（位能的二階導數）就變成位移的函數，於是第一個聲子造成的位移場會讓彈性常數在空間與時間上跟著波動。原文的「modulates in space and time」六個字，說的就是彈性常數變成 C(r, t)。
第三步：第二個聲子看到這個調變。 第二個聲子在傳播時，看到的不是一個均勻介質，而是一個「彈性常數（因而聲速）呈週期起伏、而且這個起伏還在移動」的介質。對波動來說，折射率（這裡是聲速）的週期起伏就是光柵。
第四步：被散射，產生第三個聲子。 波被光柵散射會產生新的波向量；因為光柵在移動（相速度 ω₁/K₁），散射還會伴隨頻率位移——這正是聲光調變（acousto-optic modulation）與 Brillouin 散射的機制。原文的類比「just as from a moving three-dimensional grating」用得非常精準：三維（因為晶體是三維的）、移動（所以頻率會變）、光柵（所以波向量會變）。
回到精確陳述。 類比講完之後要記住真正的內容是：三個場的乘積。三階位能項形如 x₁x₂x₃，它同時牽涉三個模式，因此描述的是三體過程。這一點下兩節會用代數與量子語言各講一次。
（校註／分界：屬於原文的只有上面那段逐字引用；「四步拆解」與「三個場的乘積」這個精確陳述，是本站為了解讀那段話而補的，p.120 並未作此拆解。）
把「光柵」算出來：三階項如何做三波混頻補推導
Kittel 只給了文字圖像，沒有寫任何式子。但這個圖像的代數只有五行，而且不寫出來就看不出「為什麼是三次項給三聲子、四次項不給」。這裡補上。
設定：取式 (38) 的前兩項 U(x) = cx² − gx³（四次項暫時放一邊），x 是某個晶格位移場。
- 第一步：把「被調變的彈性常數」精確化。 彈性常數（有效力常數）的定義是位能對位移的二階導數：k_eff(x) ≡ d²U/dx² = 2c − 6gx。注意在諧近似下 k_eff = 2c 是常數；一旦有 −gx³，k_eff 就變成位移的線性函數。這就是「非諧交互作用讓位移可以調變彈性常數」的精確意思。
- 第二步：把第一個聲子代進去。 設第一個聲子造成位移場 x₁(r, t) = A cos(K₁·r − ω₁t)。代入得 k_eff(r, t) = 2c − 6gA cos(K₁·r − ω₁t)。這就是原文「modulates in space and time the elastic constant」的式子版本：一個在空間週期為 2π/K₁、時間週期為 2π/ω₁ 的行進中的彈性常數起伏——這正是「移動的三維光柵」。
- 第三步：寫第二個聲子的運動方程。 第二個聲子 x₂ 在這個被調變的介質中運動：M ẍ₂ = −k_eff x₂ = −2c x₂ + 6g x₁ x₂。左邊加上 2c x₂ 移項後，右邊剩下的 6g x₁x₂ 就是一個外加驅動力，而它是兩個波的乘積。
第四步：把乘積展開成和頻與差頻。 設 x₂ = B cos(K₂·r − ω₂t)。用積化和差公式 cos α cos β = ½[cos(α+β) + cos(α−β)]：
x₁x₂ = ½AB { cos[(K₁+K₂)·r − (ω₁+ω₂)t] + cos[(K₁−K₂)·r − (ω₁−ω₂)t] }。
驅動力因此包含兩個新的時空頻率：(K₁+K₂, ω₁+ω₂) 與 (K₁−K₂, ω₁−ω₂)。
- 第五步：讀出結論。 如果晶體的色散關係 ω(K) 剛好在 K₃ = K₁ + K₂ 處有 ω(K₃) = ω₁ + ω₂，這個驅動力就會共振地把能量灌進第三個模式——這正是原文寫的 ω₃ = ω₁ + ω₂（兩個聲子合併成一個）。差頻那一項則對應 ω₃ = ω₁ − ω₂，改寫成 ω₁ = ω₂ + ω₃ 就看得出來它是一個聲子分裂成兩個。兩種過程都是三聲子過程，原文只寫了合併那一支。能量守恆與波向量守恆（含 G）同時滿足時，過程才真的發生——這個「同時滿足兩個守恆律」的幾何條件，就是 p.125 Fig. 17 在畫的東西。
- 第六步：為什麼四次項不給三聲子。 把 −fx⁴ 加回去：k_eff = 2c − 6gx − 12fx²。四次項對彈性常數的調變正比於 x²，也就是兩個波的乘積；驅動力因而是 三個 波的乘積 x₁x₂x₃，對應四聲子過程（三進一出或二進二出）。所以「三階項 → 三聲子、四階項 → 四聲子」是一個嚴格的對應：位能的第 n 階項對應 n 聲子過程。四聲子過程在微擾展開中高一階，通常次要，但在高溫、或在本身非諧性很強的材料（熱電化合物、鹵化鈣鈦礦）中並不能忽略。
量子版：為什麼「三階項」對應「三聲子」（跨章連結）我的判讀
上一節的古典三波混頻已經把機制講清楚了，但「聲子」這個字在古典語言裡其實是多餘的——古典版講的是波，不是粒子。要看到「三個聲子」這個數字從哪裡冒出來，最乾淨的是用產生湮滅算符，而這套工具本書在 Appendix C（書本 p.651，Creation and Annihilation Operators）就給了。
Appendix C 的式 (32) 把原子位置算符寫成
q_s = Σ_k (ℏ/2NMω_k)^{1/2} [a_k exp(iks) + a_k⁺ exp(−iks)]，
也就是說，位移算符是「一個湮滅算符加一個產生算符」的線性組合：q ~ (a + a⁺)。這一句是整個論證的樞紐。
現在看位能的各階項：
二次項 q² ~ (a + a⁺)²。 展開有四項：aa、a⁺a、aa⁺、a⁺a⁺。其中 a⁺a 與 aa⁺ 就是「數算符」，給出能階 (n + ½)ℏω；aa 與 a⁺a⁺ 在正規座標下被對角化吸收掉。兩個算符 = 兩體，而「兩體」在同一個模式上就是自己跟自己，所以沒有真正的交互作用。
三次項 q³ ~ (a + a⁺)³。 展開有八項，依「產生幾個、湮滅幾個」分成四類：
- aaa：三個聲子同時消失。能量守恆要求 ℏω₁ + ℏω₂ + ℏω₃ = 0，而所有 ω > 0，不可能，禁止。
- a⁺aa：兩個聲子消失、一個產生。這就是 ω₃ = ω₁ + ω₂，原文寫的那一支。
- a⁺a⁺a：一個消失、兩個產生。這就是分裂過程 ω₁ = ω₂ + ω₃。
- a⁺a⁺a⁺：無中生有三個聲子。能量守恆同樣禁止。
所以在量子語言下，「三階項造成三聲子過程」不是一個需要論證的物理事實，而是算符計數：三個場算符相乘，就一定改變三個模式的佔據數。能量守恆再把八項砍到只剩兩類（合併與分裂），這兩類正好就是古典版算出來的和頻與差頻。
四次項 q⁴ 同理給四個算符，對應四聲子過程；一般地，第 n 階非諧項給 n 聲子過程。
這個對應在做文獻閱讀時非常好用：看到「three-phonon scattering rate」就知道在算三階力常數，看到「four-phonon correction」就知道在算四階力常數，兩者在計算成本上差好幾個量級（三階力常數張量已經是 3N × 3N × 3N）。
Thermal expansion 小節的設定：原文刻意做的三個簡化
小節開頭的方法宣告是：「We may understand thermal expansion by considering for a classical oscillator the effect of anharmonic terms in the potential energy on the mean separation of a pair of atoms at a temperature T.」（p.120）
這句話裡藏了三個刻意的簡化，每一個都決定了結果的適用範圍，讀的時候要先標記出來，否則後面會不知道結果為什麼在低溫失效。
簡化一：classical oscillator。 整段計算裡不會出現 ℏ。用的是 Boltzmann 分布 e^{−βU}，不是量子的能階求和；沒有零點運動，沒有量子佔據數 ⟨n⟩。這保證了結果在高溫（T > θ）可信，在低溫必然錯——而且我們馬上就會看到它錯在哪裡（膨脹係數在 T → 0 不趨於零）。
簡化二：a pair of atoms。 這是一個一維、單一座標的模型，不是一個晶體。真實晶體有 3N 個模式，每個模式有自己的非諧係數，而且橫向模式與縱向模式對膨脹的貢獻符號可以相反（這是負熱膨脹材料的機制）。這個模型算出的是「一個典型鍵在溫度 T 下平均伸長多少」，把它當成晶格常數的相對變化是一個額外的、未言明的假設。
簡化三：mean separation。 這是三個簡化裡最重要、也最容易被讀者誤解的一個。要算的是平均位移 ⟨x⟩，不是「位能最小值的位置」。原文隨後的措辭是「The minimum at x = 0 is not an absolute minimum, but for small oscillations the form is an adequate representation of an interatomic potential.」——請注意，這句話把「極小點就在 x = 0」當成不言而喻的前提（U(x) 與溫度無關，極小永遠釘在 x = 0，一步都不動），句子本身真正的重點則是提醒這個極小不是絕對極小。熱膨脹不是因為位能谷底移動了，而是因為在一個不對稱的谷裡，停留位置的加權平均值不等於谷底。這句話值得抄下來反覆讀：熱膨脹是統計平均的偏移，不是平衡點的偏移。
含三次與四次非諧項的雙原子位能
U(x) = cx² − gx³ − fx⁴
(38) p.120x：相對於「絕對零度平衡間距」的位移（原文：at a displacement x from their equilibrium separation at absolute zero）。c、g、f：原文明言「with c, g, and f all positive」，三者皆為正的常數。注意兩件事：(一) 二次項寫的是 cx² 而不是 ½cx²，所以有效彈簧常數是 2c，古典頻率是 ω = (2c/M)^{1/2}；(二) 三次與四次項前面帶的是負號，正號被吸收進「g、f 皆正」的約定裡。
gx³ 與 fx⁴ 各在管什麼，以及「x = 0 不是絕對極小」
原文對 (38) 的三句註解是：「The term in x³ represents the asymmetry of the mutual repulsion of the atoms and the term in x⁴ represents the softening of the vibration at large amplitudes. The minimum at x = 0 is not an absolute minimum, but for small oscillations the form is an adequate representation of an interatomic potential.」（p.120）三句話對應三件不同的事，拆開看。
「asymmetry of the mutual repulsion」——為什麼是不對稱，而且是「互斥」的不對稱。 真實的原子間位能在 x < 0（互相壓近）那一側是被互斥力支配的，而互斥力來自電子雲重疊與 Pauli 不相容原理，隨距離變化極陡（Lennard-Jones 用 R⁻¹²、Born–Mayer 用指數）。x > 0（拉開）那一側是被 van der Waals 或庫侖吸引支配的，隨距離變化緩得多。因此位能谷是左陡右緩的。−gx³ 就是這個左右不對稱的最低階描述：它在 x > 0 把位能往下拉、在 x < 0 把位能往上推，剛好複製了「壓難拉易」。
「softening of the vibration at large amplitudes」——為什麼 x⁴ 管軟化。 由 (38) 算有效力常數：k_eff(x) = d²U/dx² = 2c − 6gx − 12fx²。最後一項 −12fx² 不管 x 是正是負都是負的，也就是說振幅一大，有效力常數就被壓低，頻率 ω = (k_eff/M)^{1/2} 跟著下降。這就是「大振幅下振動軟化」。中間那一項 −6gx 則是有正有負，在一階（週期平均）上不移動頻率——但這不表示三次項不造成軟化：它在二階微擾仍給出一個永遠為負的頻率移位，與 −fx⁴ 的一階貢獻同為振幅平方階。Kittel 之所以把「軟化」掛在 x⁴ 上，是因為 x⁴ 在一階就直接壓低有效力常數。至於平均位置偏移，則只能由 g 負責（宇稱論證）——這個分工在下面算 ⟨x⟩ 時會以「宇稱」的形式再出現一次。
（校註／分界：原文對這三句的說明只到「represents the softening of the vibration at large amplitudes」為止；上面 k_eff(x) = d²U/dx² 的計算，以及下面這個定量比較，都是本站補的，不是書上的內容。用非諧振子振幅相依頻率的標準結果（Landau–Lifshitz《Mechanics》§28）把 (38) 代進去，振幅為 a 時的頻率移位是 δω = −[3f/(2Mω₀) + 15g²/(4M²ω₀³)]a²，其中 ω₀ = (2c/M)^{1/2}。兩項都是負的——三次項與四次項都造成軟化，且同為振幅平方階；四次項要到 4fc ≳ 5g² 才佔優。所以「軟化的主力是誰」取決於係數大小，不能一概歸給 f。）
順帶把因果鏈補完：ω 依賴振幅 ⟹ ω 依賴溫度 ⟹ 後果四（彈性常數與溫度無關）與後果五（高溫熱容為常數）同時破功。所以 fx⁴ 這一項雖然不管（最低階的）熱膨脹，卻是原文點名的「大振幅軟化」來源，也牽涉彈性常數軟化與熱容的非諧修正。
「not an absolute minimum」——這句話不是客套話。 U(x) = cx² − gx³ − fx⁴ 在 |x| → ∞ 時趨於 −∞。這代表數學上 x = 0 只是一個亞穩的局部極小，粒子有機會「滾出去」。這件事在下一節做 ∫₋∞^∞ 的積分時會變成一個真正的技術問題，值得單獨講。
一個技術要點：這個積分其實是發散的我的判讀
這一節書上完全沒提，但它決定了你該怎麼理解 (39) 那一步展開的合法性，而且是研究生階段該看穿的。
把 (38) 代進 Boltzmann 積分：∫₋∞^∞ dx exp[−βU(x)] = ∫₋∞^∞ dx exp(−βcx² + βgx³ + βfx⁴)。當 |x| → ∞ 時，指數裡 +βfx⁴ 主導（因為四次贏過二次），被積函數 → e^{+∞}，這個積分嚴格來說是發散的。三次項在 x → +∞ 也發散。也就是說，原文寫下的 ⟨x⟩ 定義式，字面上並不存在。
這不是 Kittel 的疏漏，而是統計物理裡一個標準且常見的情況：這是一個漸近展開（asymptotic expansion），不是收斂展開。 正確的讀法是這樣的：
第一，高斯因子設定了物理的尺度。 e^{−βcx²} 把 x 有效地限制在 |x| ≲ x_th ≡ (k_BT/c)^{1/2} 這個範圍內（這就是熱振幅）。積分的「真正貢獻區」在這個範圍內，發散來自 |x| ≫ x_th 的區域。
第二，原文的條件正是在排除那個區域。 原文在展開前先寫下條件：「For displacements such that the anharmonic terms in the energy are small in comparison with k_BT, we may expand the integrands as …」（p.120）。把 x ~ x_th 代進去，這個條件就是 g x_th³ ≪ k_BT 與 f x_th⁴ ≪ k_BT。滿足時，發散區域的貢獻被高斯因子壓成 e^{−(大數)}，在漸近展開的任何有限階都看不到——它是所謂「超越小」（beyond all orders）的項。
第三，物理上它對應什麼。 那個「跑掉」的分支對應原子真的解離。在遠低於熔點的溫度下，解離機率確實是 e^{−E_b/k_BT} 這種指數小的量，所以把它丟掉在物理上完全正當。但如果你想算的就是熔化或缺陷生成，這一項就是全部的答案——那時候就不能用漸近展開了。
第四，實務上的判準。 展開到第幾階才停？漸近級數的規則是：加到項開始變大之前停。對這個問題，只算到 g 的一階就是最佳截斷，這也是為什麼原文只留 (40) 的線性項。
記住這一段，下一節看到「為什麼 exp(−βcx²) 不能展開、exp(βgx³ + βfx⁴) 卻可以展開」時，你就不會覺得那只是一個任意的技巧。
平均位移的 Boltzmann 定義
⟨x⟩ = ∫₋∞^∞ dx · x · exp[−βU(x)] ⁄ ∫₋∞^∞ dx · exp[−βU(x)]， β ≡ 1/k_BT
p.120⟨x⟩：位移的熱平衡平均值，也就是「平均原子間距相對於絕對零度平衡間距的偏移」。β ≡ 1/k_BT 是原文自己下的定義（原書此處用的是恆等號 ≡，表示這是定義而非近似）。分母是配分函數（一維、只含位形部分；動能部分的高斯積分在分子分母中對消，所以可以完全略去，這一點原文沒說明）。
展開後的分子與分母積分
∫ dx · x · exp(−βU) ≅ ∫ dx [exp(−βcx²)](x + βgx⁴ + βfx⁵) = (3π^{1/2}/4)(g/c^{5/2})β^{−3/2} ； ∫ dx · exp(−βU) ≅ ∫ dx · exp(−βcx²) = (π/βc)^{1/2}
(39) p.120兩條式子分別是 ⟨x⟩ 定義式的分子與分母。第一條括號裡的三項 (x + βgx⁴ + βfx⁵) 是「x 乘上 1 + βgx³ + βfx⁴」的結果——每一項的冪次都比原來的展開項高一階。積分範圍是 ±∞。 校註（PDF 字符問題，已算繪原頁目視確認）：本書 PDF 在這兩條式子中使用的關係符號為 ≅（約等於），但文字層抽取工具將該自訂字符解碼為 ≡，會讓人誤以為這是定義式。同一頁 β ≡ 1/k_BT 用的是另一個字符，該處確實是 ≡。另外積分上下限的 ∞ 在文字層中以未解碼代碼 /L50525 出現。
(39) 的展開逐步：宇稱是全部的關鍵補推導
原文從 ⟨x⟩ 的定義式直接跳到 (39) 的結果，中間的四個步驟一個都沒寫。這四步裡有一步（宇稱篩選）是整個熱膨脹結論的心臟，不補出來讀者會完全看不到「為什麼是 g 而不是 f」。
設 β ≡ 1/k_BT，U(x) = cx² − gx³ − fx⁴。
第一步：把指數拆成兩個因子。 −βU(x) = −βcx² + βgx³ + βfx⁴，所以
exp(−βU) = exp(−βcx²) · exp(βgx³ + βfx⁴)。
注意第一個因子帶負號（收斂），第二個因子帶正號（發散）。這個符號差異正是下一步要小心處理的原因。
第二步：只展開第二個因子，第一個因子原封不動。
exp(βgx³ + βfx⁴) ≅ 1 + βgx³ + βfx⁴ + ⋯
原文在 p.121 專門提醒這一點：「Note that in (39) we have left cx² in the exponential, but we have expanded exp(βgx³ + βfx⁴) ≅ 1 + βgx³ + βfx⁴ + ⋯.」
為什麼不能連 exp(−βcx²) 一起展開？ 因為那個因子是唯一讓積分收斂的東西。把它展成 1 − βcx² + ½β²c²x⁴ − ⋯ 之後，每一項 ∫x^n dx 在 ±∞ 都發散，整個方法瞬間垮掉。用漸近展開的語言說：exp(−βcx²) 是這個問題的「零階解」（高斯測度），非諧項才是微擾。 把測度本身展開，等於把微擾論的基準拆了。
第三步：寫出分子的被積函數。 分子要乘一個額外的 x：
x · exp(−βcx²) · (1 + βgx³ + βfx⁴) = exp(−βcx²) · (x + βgx⁴ + βfx⁵)。
這就是原文 (39) 第一行括號裡那三項的來源——每一項的冪次都比展開項高一階，因為多乘了一個 x。這一步沒有任何物理，純粹是代數，但如果沒看出「x 被乘進括號裡」，(39) 的 x⁴ 與 x⁵ 會顯得莫名其妙。
第四步（核心）：宇稱篩選。 積分區間 (−∞, ∞) 是對稱的，而 exp(−βcx²) 是偶函數。因此「偶函數 × 偶函數」才活得下來，「偶函數 × 奇函數」一律積分為零。逐項檢查分子的三項：
項 冪次 宇稱 積分結果 x 1，奇 奇 0 βgx⁴ 4，偶 偶 存活 βfx⁵ 5，奇 奇 0 只有 g 那一項活下來。 這就是「最低階的熱膨脹只牽涉反對稱項 gx³、不牽涉對稱項 fx⁴」的完整證明，而且只需要一行宇稱論證。注意這裡的邏輯順序：gx³ 是奇函數（打破 x → −x 對稱），乘上額外的 x 之後變成偶函數而存活；fx⁴ 是偶函數（不打破對稱），乘上 x 之後變成奇函數而死掉。能製造 ⟨x⟩ ≠ 0 的，必須是打破 x → −x 對稱的項——這是一個對稱性論證，不是一個計算結果。
- 第五步：分母。 分母的被積函數是 exp(−βcx²)(1 + βgx³ + βfx⁴)。同樣做宇稱篩選：βgx³ 是奇的 → 0，剩下 1 與 βfx⁴ 兩項。原文把 βfx⁴ 也丟掉，只寫 ∫dx exp(−βcx²)。這一步原文沒有說明理由，補上：分母的 f 修正會給 ⟨x⟩ 一個相對修正因子 [1 + O(f/βc²)]⁻¹，而分子本身已經正比於 g，兩者相乘是 O(fg) 的二階小量；既然我們只要 g 的一階結果，這一項就該丟。
- 第六步：確認下一階含什麼（為什麼「lowest order」這四個字必須加）。 把展開做到二階：exp(βgx³ + βfx⁴) = 1 + βgx³ + βfx⁴ + ½β²g²x⁶ + β²gf x⁷ + ½β²f²x⁸ + ⋯。乘上 x 之後做宇稱篩選，下一個存活的是 β²gf x⁸（偶）；三階展開的 (1/6)β³g³x⁹ 乘 x 後成為 x¹⁰ 也存活。量綱清點後這兩者都給出正比於 T² 的修正。結論：⟨x⟩ = (3g/4c²)k_BT + O(T²)，而 T² 項同時含 gf 與 g³ 兩種組合。 所以 f 並不是永遠不參與熱膨脹，它只是在最低階不參與——p.121 那句話裡的「In lowest order」四個字是不能省的。
兩個高斯積分：書上完全跳過的那一步補推導
原文 (39) 直接寫出兩個積分的答案 (3π^{1/2}/4)(g/c^{5/2})β^{−3/2} 與 (π/βc)^{1/2}，一個字都沒解釋怎麼算的。這一步用的是標準的高斯矩技巧（對參數微分），值得完整寫一次，因為固態物理後面還會用很多次（Debye–Waller 因子、擴散、漲落理論全都要）。
目標：算出 ∫₋∞^∞ e^{−ax²}dx 與 ∫₋∞^∞ x⁴e^{−ax²}dx，然後代 a = βc。
第一步：基本高斯積分。 I₀(a) ≡ ∫₋∞^∞ e^{−ax²} dx = (π/a)^{1/2} = π^{1/2} a^{−1/2}。
代 a = βc 立刻得到 原文 (39) 的分母 (π/βc)^{1/2} ✓。
第二步：用「對參數微分」生出高階矩。 關鍵觀察：對 a 微分一次，就從被積函數裡拉下一個 (−x²)：
dI₀/da = ∫₋∞^∞ (−x²) e^{−ax²} dx ⟹ ∫ x² e^{−ax²} dx = −dI₀/da。
這比硬做分部積分快得多，也不用查表。
第三步：算到四次矩。 對 I₀(a) = π^{1/2}a^{−1/2} 連續微分兩次：
dI₀/da = π^{1/2} · (−½) a^{−3/2} = −½ π^{1/2} a^{−3/2}；
d²I₀/da² = −½π^{1/2} · (−3/2) a^{−5/2} = (3/4) π^{1/2} a^{−5/2}。
而 d²I₀/da² = ∫ x⁴ e^{−ax²} dx（微分兩次拉下 (−x²)² = x⁴，兩個負號相乘為正）。所以
∫₋∞^∞ x⁴ e^{−ax²} dx = (3/4) π^{1/2} a^{−5/2}。
（順帶記下 ∫x²e^{−ax²}dx = ½π^{1/2}a^{−3/2}，這一條在算 ⟨x²⟩ = k_BT/2c 時就要用。）
第四步：組裝分子。 分子存活的只有 βg∫x⁴e^{−βcx²}dx：
βg · (3/4)π^{1/2}(βc)^{−5/2} = (3π^{1/2}/4) · g · β^{1} · β^{−5/2} · c^{−5/2} = (3π^{1/2}/4)(g/c^{5/2}) β^{−3/2}。
與原文 (39) 第一行的答案逐字符合 ✓。注意 β 的冪次：1 − 5/2 = −3/2，這是最容易算錯的一步。
第五步：相除得到 (40)。
⟨x⟩ = (3π^{1/2}/4)(g/c^{5/2})β^{−3/2} ÷ (π/βc)^{1/2}
= (3π^{1/2}/4)(g/c^{5/2})β^{−3/2} × π^{−1/2}(βc)^{1/2}
= (3/4) · g · c^{−5/2} · c^{1/2} · β^{−3/2+1/2}
= (3/4)(g/c²) β^{−1}
= (3g/4c²) k_BT ✓
π 為什麼消失了：分子帶 π^{1/2}、分母帶 π^{1/2}，正好對消——這就是 (40) 裡看不到任何 π 的原因。c 的冪次：c^{−5/2} × c^{1/2} = c^{−2}，所以是 g/c²，不是 g/c^{5/2}。β 的冪次：−3/2 + 1/2 = −1，所以 β^{−1} = k_BT，⟨x⟩ 線性正比於 T。
第六步：量綱檢查（強烈建議每次都做）。 由 U = cx² − gx³ − fx⁴ 且 U 是能量 E、x 是長度 L：
[c] = E/L²，[g] = E/L³，[f] = E/L⁴。
於是 [g/c²] = (E·L⁻³)/(E²·L⁻⁴) = L/E。乘上 [k_BT] = E 得到 L ✓。⟨x⟩ 確實是一個長度。如果你哪一步把 c^{5/2} 抄錯成 c²，這個檢查會立刻抓到。
古典區的熱膨脹
⟨x⟩ = (3g / 4c²) k_B T
(40) p.120⟨x⟩：溫度 T 下的平均位移，也就是一個鍵相對於絕對零度的平均伸長量。g：(38) 的三次非諧係數（不對稱性的度量）。c：(38) 的二次諧係數（有效彈簧常數為 2c）。k_B：Boltzmann 常數。原文在 (40) 之後緊接著補上限定語：「in the classical region.」（p.121）
p.121 的兩句提醒（逐字），以及它們為什麼放在這裡
熱膨脹小節的正文在 p.120 結束於 (40)，但 Kittel 把兩句關鍵的註解放到了 p.121 的開頭。做題或複習時很容易翻頁翻過去，這兩句其實是本節最該畫線的地方。逐字如下（經算繪原頁目視確認）：
第一句（關於展開的分寸）：「Note that in (39) we have left cx² in the exponential, but we have expanded exp(βgx³ + βfx⁴) ≅ 1 + βgx³ + βfx⁴ + ⋯.」
這一句在守的是方法論：哪一部分留在指數裡、哪一部分可以展開，不是隨便選的。留在指數裡的必須是保證收斂的高斯測度；可以展開的必須是在該測度支撐範圍內確實是小量的部分。這個「保留高斯核、展開微擾」的模式，是統計場論裡整套微擾展開的雛形——之後讀到 Feynman 圖展開時，「高斯部分＝自由傳播子、展開部分＝頂點」就是同一件事的專業版。
第二句（關於哪一項負責）：「In lowest order the thermal expansion does not involve the symmetric term fx⁴ in U(x), but only the antisymmetric term gx³.」
這一句是本節最容易被記反的結論，值得完整咀嚼。注意 Kittel 用的兩個形容詞：symmetric（對稱，指 fx⁴ 在 x → −x 下不變）與 antisymmetric（反對稱，指 gx³ 在 x → −x 下變號）。他不是在說「四次項不重要」，而是在說「四次項不破壞對稱性，所以它無法製造一個非零的平均位移」。這是一個對稱性論證，不是一個數值大小的比較——就算 f 比 g 大一百倍，最低階的熱膨脹還是完全由 g 決定。
另外請留意「In lowest order」這個限定詞。如前面第六步所算，展開到下一階時 f 會透過 gf 的交叉項進來，給出正比於 T² 的修正。所以精確的陳述是：f 不能單獨造成熱膨脹，但它可以修飾由 g 造成的熱膨脹。
「不牽涉 fx⁴」為什麼常被記反我的判讀
這是這一段最常見的記憶錯誤，值得專門拆解為什麼大腦會記反。
錯誤記憶的來源。 學生在別的地方（非線性振盪、Duffing 振子、彈性常數的溫度依賴、熱容的非諧修正）最常遇到的「非諧項」是四次項——因為四次項是最低階的、保持對稱性的非諧項，所以在對稱位能（例如雙井、單井的軟化）的討論裡，四次項是主角。於是「非諧 = x⁴」變成一個下意識的聯想，遇到熱膨脹時就直接套上去。
為什麼這裡反過來。 熱膨脹問的問題是 ⟨x⟩ 等不等於零，這是一個關於對稱破缺的問題。要讓一個原本對稱分布的平均值偏離中心，你必須有一個打破 x → −x 對稱的東西。四次項是偶的，它可以讓分布變胖變瘦、讓頻率軟化、讓 ⟨x²⟩ 改變，但它在原點左右做的事一模一樣，所以平均值仍然在原點。只有奇次項（x³、x⁵、…）能把分布往一側推。
一句話記法。
- 偶次非諧項（fx⁴）改變的是「寬度」與「頻率」 → 管彈性常數軟化、熱容的非諧修正、四聲子散射（高階修正）。
- 奇次非諧項（gx³）改變的是「中心」 → 管熱膨脹，也管三聲子過程，因而是聲子壽命與熱阻的主要來源（p.120 逐字：「Three-phonon processes are caused by third-order terms in the lattice potential energy.」）。
同一個 gx³ 為什麼同時管兩件事。 注意本段前面已經出現過兩次 gx³：一次是三波混頻（k_eff = 2c − 6gx 的線性調變），一次是熱膨脹（宇稱篩選存活的唯一項）。這不是巧合：三次項是唯一「既能單獨破壞對稱性、又只牽涉三個場」的最低階項，所以熱膨脹與三聲子過程是同一個係數的兩個後果。這正是 Grüneisen 常數 γ 能同時出現在熱膨脹公式與熱導率公式裡的深層原因——兩者共用同一個 g。
一個反向的檢查題。 如果某材料的位能剛好是 U = cx² − fx⁴（g = 0），它會不會熱膨脹？答：最低階不會（⟨x⟩ = 0），而且因為 (38) 的展開裡再高階的存活項都至少含一個 g，在這個模型內它完全不膨脹，但它的振動頻率仍然會隨溫度下降。這個「不膨脹卻軟化」的組合，正好把兩個後果的獨立性展示得很清楚。
Fig. 15 在說什麼（網頁無圖，用文字還原）
原文對這張圖的文字只有三句：「Measurements of the lattice constant of solid argon are shown in Fig. 15. The slope of the curve is proportional to the thermal expansion coefficient. The expansion coefficient vanishes as T → 0, as we expect from Problem 5.」（p.121）圖說是「Figure 15 Lattice constant of solid argon as a function of temperature.」
以下是把該圖算繪放大後的目視判讀（本站不放圖，故用文字還原；數值為讀圖估計，非原書印出的數字）：
軸的設定。 橫軸是溫度，0 到約 85 K，刻度每 20 K 一格。左縱軸是晶格常數，單位 Å，由下而上 5.30、5.34、5.38、5.42、5.46。右縱軸是密度，單位 g cm⁻³，由下而上 1.78、1.74、1.70、1.66、1.62——注意右軸是反向的（密度越往上越小），因為晶格常數變大就是密度變小，兩軸是同一組數據的兩種讀法。這個雙軸設計本身就是在提醒：熱膨脹等價於熱致密度下降。
曲線的形狀（這是全圖的重點）。 資料點從約 4 K 開始，一路到約 84 K 的三相點（圖右上角有標註「Triple point」並用箭頭指到曲線末端）。曲線在低溫端幾乎是水平的：4 K 與 15 K 兩個點都落在 5.300 Å 附近，肉眼幾乎看不出差別。從約 20 K 開始才明顯抬頭，之後斜率持續變大，呈現向上彎曲的形狀：20 K 約 5.306、30 K 約 5.318、40 K 約 5.335、50 K 約 5.360、60 K 約 5.383、70 K 約 5.415、80 K 約 5.447，到三相點約 5.46。整段從 5.30 漲到 5.46，相對變化約 3%。
這張圖在論證什麼。 因為斜率正比於膨脹係數（原文明說），所以「曲線低溫端平坦」就等於「膨脹係數在 T → 0 時趨於零」。而 (40) 預測膨脹係數是一個與溫度無關的常數——曲線應該是一條直線。它不是。這是本段刻意留下的一個正面衝突，也是原文立刻把讀者導向 Problem 5 的理由。
為什麼挑固態氬。 原文沒說明，但理由很清楚：氬是最接近「純 Lennard-Jones 固體」的真實材料（第 3 章就是用惰性氣體晶體來示範 LJ 位能），沒有共價鍵、沒有離子電荷、沒有導電電子的貢獻，非諧效應完全來自兩體中心力位能。用它來檢驗 (38) 這種兩體位能模型，是最乾淨的選擇。另外它的 θ 只有約 90 K，三相點 84 K 也在同一量級，所以整張圖橫跨了「量子區」到「接近古典區」，正好把 (40) 成立與失效的兩段一次畫出來。
古典理論在 T → 0 破功：(40) 與 Fig. 15 的正面衝突我的判讀
把前面兩節擺在一起，就得到一個乾淨的矛盾，這個矛盾是這一段最有教學價值的地方，值得寫清楚。
(40) 的預測：⟨x⟩ = (3g/4c²)k_BT，所以 d⟨x⟩/dT = 3gk_B/4c²，與 T 無關。畫成「晶格常數 vs 溫度」就是一條斜率固定的直線，一路延伸到 T = 0。
Fig. 15 的事實：曲線在低溫端平坦，斜率隨溫度上升而增大，明顯不是直線。原文直說「The expansion coefficient vanishes as T → 0」。
錯在哪裡？不在非諧，在古典。 這一點務必分清楚，因為很多人會誤以為是 (38) 那個位能形式不夠好。實際上問題出在前面標記的「簡化一」：整段計算用的是古典 Boltzmann 分布，其中沒有 ℏ。古典統計裡，每個自由度不管頻率多高，在任何非零溫度都帶有 k_BT 的能量——沒有任何模式會被凍結。因此膨脹的「驅動力」與溫度成正比，膨脹係數是常數。
真實晶體在低溫時，ℏω ≫ k_BT 的模式佔據數 ⟨n⟩ = 1/(e^{ℏω/k_BT} − 1) 指數地趨於零，這些模式完全不參與熱膨脹。能參與的只有 ℏω ≲ k_BT 的低頻模式，而它們的數目在 Debye 模型下隨 T³ 成長。這正是熱容在低溫遵守 T³ 定律（本書 p.114）的同一個機制。
正確的低溫結果。 Problem 5（p.129）的 Grüneisen 處理會給出
體膨脹係數 β_V = γ C_V / B，
也就是膨脹係數正比於熱容。既然 C_V ∝ T³（Debye T³ 定律），膨脹係數也 ∝ T³ → 0，Fig. 15 的平坦低溫端就被解釋了。而在高溫 C_V → 3Nk_B 為常數，膨脹係數也回到常數——與 (40) 的結論一致。所以 (40) 不是錯的，它是 Grüneisen 結果的高溫極限。
這個「膨脹係數 ∝ 熱容」的關係非常有用，因為它是一個跨材料的經驗律（Grüneisen 定律）：把實測的 α(T) 與 C_V(T) 相除，得到的 γ 在很寬的溫度範圍內近乎常數。實務上這是一個很好的資料檢查：如果你量到的 α(T)/C_V(T) 隨溫度劇烈變化，通常表示樣品有相變、缺陷或量測系統誤差。
Problem 5（Grüneisen 常數）：低溫的修法在哪裡
原文把低溫的正確處理推給 Problem 5（書本 p.129，標了星號表示「somewhat difficult」）。既然正文明白地說「as we expect from Problem 5」，這一題就不是選做題，而是本節論證的一部分。以下整理該題的結構（內容依 p.129 原文，符號按原書）。
(a) 要證明頻率為 ω 的一個聲子模式，其自由能為 k_BT ln[2 sinh(ℏω/2k_BT)]。原文特別註明：必須保留零點能 ½ℏω 才能得到這個結果。這一點很關鍵——它正是「量子」進入計算的入口，而 (40) 的古典算法裡完全沒有這一項。
(b) 設 Δ 為體積的相對變化量（fractional volume change），晶體的自由能寫成
F(Δ, T) = ½BΔ² + k_BT Σ ln[2 sinh(ℏω_K/2k_BT)]，
B 是體模數。第一項是彈性能（把晶體撐開要付的代價），第二項是所有聲子模式的自由能總和。接著假設 ω_K 的體積依賴性為
δω/ω = −γΔ，
γ 即 Grüneisen 常數。若把 γ 視為與模式 K 無關，要證明 F 對 Δ 取極小時滿足
BΔ = γ Σ ½ℏω coth(ℏω/2k_BT)，
並進一步用熱能量密度改寫成 Δ = γU(T)/B。
(c) 要證明在 Debye 模型下 γ = −∂ln θ / ∂ln V。
原文結尾的警語也要讀：「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.」
這一題怎麼修好 (40) 的毛病，用一句話說。 整個修正就藏在 coth 裡。Σ ½ℏω coth(ℏω/2k_BT) 恰好就是量子諧振子的平均能量總和（因為 ½ℏω coth(ℏω/2k_BT) = ℏω(⟨n⟩ + ½)）。在高溫它趨於 Σ k_BT = 3Nk_BT，於是 Δ ∝ T，膨脹係數是常數——回到 (40)。在低溫它趨於零點能加上按 Planck 分布加權的熱能，熱能部分 ∝ T⁴，於是 Δ 的溫度依賴部分 ∝ T⁴，膨脹係數 ∝ T³ → 0——解釋 Fig. 15。
為什麼 γ 是「非諧性」的度量。 由定義 δω/ω = −γΔ：在嚴格的諧晶體中，力常數與體積無關，因此 ω 與體積無關，γ ≡ 0，於是 Δ ≡ 0，沒有熱膨脹。所以 γ ≠ 0 與「位能有三次以上的項」是同一件事的兩種說法。這也把 p.119 的第二條與第四條後果統一了起來：彈性常數（因而 ω、因而 θ）隨體積改變，就是熱膨脹存在的充要條件。 (c) 小題要證的 γ = −∂ln θ/∂ln V 正是這句話的精確版本。
Grüneisen 關係（Problem 5 的結果）
Δ = γ U(T) / B ， 其中 BΔ = γ Σ ½ℏω coth(ℏω / 2k_BT) ， 且（Debye 模型）γ = −∂ ln θ / ∂ ln V
p.129Δ：體積的相對變化（fractional volume change）。γ：Grüneisen 常數，由 δω/ω = −γΔ 定義，衡量聲子頻率對體積的敏感度。B：體模數。U(T)：熱能量密度（由 Σ ½ℏω coth(ℏω/2k_BT) 給出，該式等於 Σ ℏω(⟨n⟩ + ½)，含零點能）。校註：原書此處未交代體積歸一化——BΔ 的量綱是能量密度，而 Σ 字面上是對所有模式求和（總量），兩者差一個體積因子；正確的讀法是把這裡的 Σ 理解為每單位體積的求和。本檔一律以 U(T)、C_V 表示每單位體積的熱能與熱容。θ：Debye 溫度。V：體積。
量級自洽檢查：用 Fig. 15 的斜率反推固態氬的 γ補推導
把 (40)、Fig. 15 與 Grüneisen 關係串起來做一次數字檢查，可以確認自己沒有把哪個因子記錯。注意：以下用到的體模數與莫耳體積不是本書這幾頁給的，屬於外部數據，這裡只作量級自洽性檢查，不是書上的結果。
- 第一步：從 Fig. 15 讀線膨脹係數。 由前面的讀圖估計，在 50–80 K 區間 da/dT ≈ 2.9 × 10⁻³ Å/K，a ≈ 5.40 Å，故 α = (1/a)(da/dT) ≈ 5.4 × 10⁻⁴ K⁻¹。體膨脹係數 β_V = 3α ≈ 1.6 × 10⁻³ K⁻¹。（這個值大得驚人——比典型金屬大一個數量級——因為氬是靠 van der Waals 鍵結的分子晶體，鍵很軟。）
第二步：算單位體積熱容。 fcc 氬，一個立方晶胞邊長 a ≈ 5.40 Å 含 4 個原子，故每原子體積 = a³/4 = 157.5/4 ≈ 39.4 ų = 3.94 × 10⁻²⁹ m³。在 80 K（接近 θ ≈ 90 K，已相當接近古典區）取 Dulong–Petit 值 C_V = 3k_B/原子：
C_V/V = 3k_B / 3.94×10⁻²⁹ m³ = 3 × 1.38×10⁻²³ / 3.94×10⁻²⁹ ≈ 1.05 × 10⁶ J m⁻³ K⁻¹。
第三步：用 Grüneisen 關係反解 γ/B。 由 β_V = γ (C_V/V)/B（C_V/V 是每單位體積的熱容）得
γ/B = β_V /(C_V/V) = 1.6×10⁻³ / 1.05×10⁶ ≈ 1.5 × 10⁻⁹ Pa⁻¹。
第四步：代入固態氬的體模數（外部數據）。 固態氬在低溫的體模數約 2 GPa 的量級。取 B ≈ 2 × 10⁹ Pa：
γ ≈ 1.5 × 10⁻⁹ × 2 × 10⁹ ≈ 3.0。
- 第五步：判讀。 大多數固體的 Grüneisen 常數落在 1.5–3 之間，惰性氣體晶體偏高端（因為 LJ 位能的不對稱性強），所以 γ ≈ 3.0 是一個完全合理的值。這代表整條鏈（Fig. 15 的讀圖 → 體膨脹係數 → Grüneisen 關係）在數量級上自洽，也順帶確認了 β_V = 3α 這個因子 3 沒有記錯——如果誤用 β_V = α，算出的 γ 會變成 1.0，落在惰性氣體晶體不該有的位置。
- 第六步：回頭看 (38) 的係數。 由 (40)，d⟨x⟩/dT = 3gk_B/4c²。取最近鄰間距 d = a/√2 ≈ 3.82 Å，d⟨x⟩/dT = αd ≈ 2.1 × 10⁻³ Å/K，於是 g/c² ≈ (4/3)(1.9×10⁻³ Å/K)/k_B。要把這個組合拆成 g 與 c 各自的值，還需要另一個獨立的資訊（例如由 θ 或聲速定出 c），本書這幾頁沒有提供——這也說明 (38) 這個模型的實用性在於提供標度關係（⟨x⟩ ∝ T、⟨x⟩ ∝ g/c²），而不是提供可直接查表的參數。
這一段的常見誤解清單我的判讀
誤解一：「熱膨脹是因為溫度高時原子振動振幅變大。」 振幅變大是必要但不充分。純諧振子的振幅也隨溫度變大（⟨x²⟩ = k_BT/2c），但它一點都不膨脹。真正的原因是位能不對稱，振幅變大只是把這個不對稱「曝光」出來。正確的一句話版本是：振幅變大 + 位能左陡右緩 = 平均位置右移。
誤解二：「熱膨脹來自四次項。」 最低階完全來自三次項。四次項是對稱的，不能製造非零平均位移。p.121 原文明說。四次項管的是頻率軟化（彈性常數的溫度依賴、熱容的非諧修正），不是膨脹。
誤解三：「溫度升高使位能極小值往外移。」 位能 U(x) 是一個與溫度無關的函數，它的極小值永遠在 x = 0（原文「The minimum at x = 0 is not an absolute minimum…」這句就是這樣預設的）。移動的是統計平均值 ⟨x⟩，不是極小點。這個區別在概念上是本節的核心。
誤解四：「(40) 說膨脹係數正比於 T。」 不對。(40) 說的是 ⟨x⟩ 正比於 T，因此膨脹係數 d⟨x⟩/dT 是常數。把兩者搞混，會在 Fig. 15 的討論上完全讀反。
誤解五：「Fig. 15 顯示古典理論不對，所以 (38) 的模型是錯的。」 不對。Fig. 15 顯示的是古典統計在低溫不適用，不是非諧位能模型錯。把同一個非諧假設（γ ≠ 0）配上量子統計（Problem 5 的 coth），就能同時解釋低溫的 T³ 與高溫的常數。
誤解六：「所有材料都熱膨脹。」 (38) 假設 g > 0。實際材料中，某些模式（特別是低頻的橫向彎曲模式，例如兩端被拉住的原子鏈橫向擺動）對體積的貢獻使有效 γ 為負，於是出現負熱膨脹。本書這幾頁（p.119–121）沒有討論負熱膨脹（我只確認了這幾頁，沒有對全書下結論）。
誤解七：「ω₃ = ω₁ + ω₂ 是唯一的三聲子過程。」 原文只寫了合併這一支，但分裂過程 ω₁ = ω₂ + ω₃ 同樣由三次項給出（見前面三波混頻推導的差頻項，與量子版的 a⁺a⁺a 類）。兩者在熱傳導的計算裡都要算進去。
誤解八：「三聲子過程一定造成熱阻。」 不對，這是第 5 章後面 Umklapp 一節的重點。滿足 K₁ + K₂ = K₃（N 過程）的碰撞保持總聲子動量，不會讓熱流衰減；只有 K₁ + K₂ = K₃ + G（U 過程，式 (47)，p.125）才會翻轉聲子流的方向而造成熱阻。「非諧 ⟹ 有散射」與「有散射 ⟹ 有熱阻」之間還差一個 G。
誤解九：「聲子是理想氣體。」 只在諧近似下成立。「聲子氣體」這個詞在第 5 章後半（Thermal resistivity of phonon gas, p.123）會繼續使用，但那時已經是一個有交互作用的氣體——用動力學理論 K = ⅓Cvℓ 的前提，正是粒子會碰撞。
往下讀：這一段接到哪裡我的判讀
直接往下（同章）。 p.121 的 THERMAL CONDUCTIVITY 一節接的是後果一。式 (41) j_U = −K dT/dx 一寫下來，就等於承認聲子會被散射；p.122 的式 (42) K = ⅓Cvℓ 用動力學理論把問題化約成「平均自由程 ℓ 是多少」；p.123 的 Thermal resistivity of phonon gas 接著說明 ℓ 由幾何散射與聲子—聲子散射決定、高溫下 ℓ ∝ 1/T，並用 (45)(46) 指出 N 過程守住總聲子動量因而不造成熱阻；p.125–126 的 Umklapp processes 指出只有 G ≠ 0 的碰撞（式 (47)）才貢獻熱阻，並在 p.126 給出「可參與 U 過程的高能聲子數目 ∝ exp(−θ/2T)」，因而 ℓ ∝ exp(θ/2T)；p.126 的 Imperfections 則處理當 ℓ 大到與樣品尺寸相當時的 size effect，此時 K ≅ CvD（式 (48)）。這一整串的起點，就是本段第一條後果的破功。
回頭連（前面幾章）。 後果四與後果五在第 3 章（彈性常數）與第 5 章前半（Debye 模型）被默默使用。特別是 Debye 模型的整個架構建立在「有一個固定的聲速 v、因而有一個固定的 θ」上；一旦承認 γ = −∂ln θ/∂ln V ≠ 0，θ 本身就隨體積（因而隨溫度）改變，這正是實驗上擬合出的 θ(T) 不是水平線的原因。
往後連（後面幾章與附錄）。 三聲子過程的算符語言在 Appendix C（p.651，Creation and Annihilation Operators）有完整的工具；⟨x²⟩ 這個「諧近似仍然給非零」的量在本章 Problem 2（Rms thermal dilation of crystal cell）與 Problem 3（Zero point lattice displacement and strain）處理，並在 X 光繞射的 Debye–Waller 因子中出現。
與熱容／熱阻萃取類研究的接點（本站補充）。 如果你的研究是元件層級的熱行為萃取（例如 FinFET 的熱容 C_th 與熱阻 R_th），這一段提供的是物理量的出處：
- C_th 的上限來自後果五。 高溫 Dulong–Petit 值 3Nk_B 給出體積比熱的天花板，任何萃取出「超過 3Nk_B 很多」的體積比熱都該先懷疑是量測把界面熱容或寄生熱質量算進去了。反過來，在薄鰭／奈米線中低頻模式被幾何截斷，低溫比熱會偏離體材料的 T³。
- R_th 的物理完全落在後果一。 元件自熱的熱阻由聲子平均自由程決定，而 p.125–126 的兩個機制在奈米尺度同時起作用：Umklapp 決定體材料的本徵 ℓ（矽室溫約數百 nm 量級），size effect（式 (48) K ≅ CvD）則在鰭寬只有幾 nm 時主宰一切——這正是 FinFET 熱阻遠高於體矽預期值的第一性原理來源。
- 熱膨脹決定應力與溫度係數。 (40) 與 Grüneisen 關係是 TCR、通道應變隨溫度漂移、以及不同材料層之間熱應力的物理起點。
讀完這一段之後，第 5 章剩下的部分（p.121–127）在邏輯上就只是「把後果一的破功量化」而已。