這一段在解決什麼問題我的判讀
第 4 章從 p.91 一路算到 p.100，做的都是同一件事：給定力常數，算出 ω(K)。單原子鏈給 (7)(9)，長程力給 (16a)(17)，雙原子鏈給 (22)(23)(24)(25)，量子化給 (27)(29)，聲子動量給 (30)(31)(32)。到目前為止，所有的 ω(K) 都是理論輸出。
但早在 p.94 的 Derivation of Force Constants from Experiment 小節就埋了一句話：「We can make a statement about the range of the forces from the observed experimental dispersion relation for ω.」（此句跨 p.94／p.95 分頁，「experimental」與「dispersion」之間即為頁界）——要用 (17) 反解力常數 C_p，你得先有實驗量到的 ω(K)。這一段（p.100 下半到 p.102 上半）就是來還這筆債的：ω(K) 從哪裡來。
答案只有一頁半，而且只有兩條式子 (33) 與 (34)。但這兩條式子合起來就是整個中子散射光譜學的全部運動學基礎，值得逐字拆開來看。拆完之後，這一段還附了兩張圖（Fig. 11 的鈉色散曲線、Fig. 12 的三軸譜儀照片），一份四條的 SUMMARY，以及七道習題。
這七道習題不是裝飾。正文為了篇幅省略了幾塊重要的物理，而習題把它們補了回來——特別是 Problem 3（區邊界的兩個子晶格解耦，補完 (25) 的物理圖像）、Problem 4（Kohn anomaly，長程力的直接後果，接費米面）、以及 Problem 7（soft phonon modes，ω² 變負、晶格失穩，接第 16 章的結構相變）。這三題在後面會逐題講透，包含書上完全沒寫的每一步代數。
引用體例（請先讀一次）。本文中凡標「逐字」的英文引句，都直接取自書本該頁，並只做三種不改變任何文字的還原：(i) 去斷行連字號——PDF 文字層把跨行的字拆成「experimen- tally」，還原為 experimentally；(ii) 把行間公式放回它在紙面上的位置——PDF 文字層會把 displayed equation 浮到段末，例如 SUMMARY 第二條在文字層讀作「…the wavevector selection rule that governs the process is where G is a reciprocal lattice vector.」，紙面上「is」與「where」之間印著 k = k′ + K + G，引用時放回原位；(iii) 還原上下標與分數——文字層的「p 2/2Mn」還原為 p²/2M_n，「1 2 ℏω」還原為 ½ℏω。除此之外一字不改，跨頁斷開的句子會另外註明。
為什麼是中子：原文怎麼說我的判讀
這一節的標題 INELASTIC SCATTERING BY PHONONS 印在 p.100 下半，內容延續到 p.102。開場兩句話是整段的定調：
「Phonon dispersion relations ω(K) are most often determined experimentally by the inelastic scattering of neutrons with the emission or absorption of a phonon.」
「A neutron sees the crystal lattice chiefly by interaction with the nuclei of the atoms.」（p.100 末至 p.101 首，這句被分頁切在兩頁之間）
第二句是關鍵，而且常被讀過去。它說的是：中子沒有電荷，所以它不跟電子雲玩，它直接撞原子核。這一點的重要性有三層：
第一，X 光（第 2 章）的散射振幅正比於電子密度的傅立葉分量，所以 X 光量到的是電荷分布；中子的散射振幅由核散射長度決定，所以中子量到的是原子核的位置。而晶格振動就是原子核在動——中子看的正好是你想量的那個東西。
第二，因為交互作用是短程的核力（作用範圍 ~fm，遠小於中子波長 ~Å），散射是弱的、可以用一階微擾（Born 近似）處理，也就是單次散射主導。這讓 (33) 這種「一進一出、外加一個聲子」的簡單運動學圖像真的成立，不必擔心多重散射把訊號攪糊。
第三，弱交互作用的代價是穿透力強、截面小——中子可以穿過整塊樣品（所以量的是塊材性質，不是表面），但也因此需要反應器或散裂源等級的中子通量，實驗才做得動。Fig. 12 那台在 Brookhaven 的機器就是這麼來的。
（第二、三層是本講解補充的說明；Kittel 在這裡只寫了「chiefly by interaction with the nuclei」這一句。）
中子散射的波向量選擇律
k + G = k′ ± K
(33) p.101k：入射中子的波向量；k′：散射後中子的波向量；K：散射過程中被產生（取 + 號）或吸收（取 − 號）的聲子波向量；G：任一倒晶格向量。原文原話：「Here K is the wavevector of the phonon created (+) or absorbed (−) in the scattering process, and G is any reciprocal lattice vector. For a phonon we choose G such that K lies in the first Brillouin zone.」
中子散射的能量守恆
ℏ²k²/2M_n = ℏ²k′²/2M_n ± ℏω
(34) p.102M_n：中子質量；ℏ²k²/2M_n：入射中子動能；ℏ²k′²/2M_n：散射中子動能；ℏω：被產生（+）或吸收（−）的聲子能量。原文的鋪陳：「The kinetic energy of the incident neutron is p²/2M_n, where M_n is the mass of the neutron. The momentum p is given by ℏk, where k is the wavevector of the neutron.」
把 (33) 與 (34) 聯立：為什麼「兩條式子」剛好夠用補推導
Kittel 只寫了一句操作指令：「To determine the dispersion relation using (33) and (34) it is necessary in the experiment to find the energy gain or loss of the scattered neutrons as a function of the scattering direction k – k′.」（p.102）為什麼是這兩條、為什麼剛好夠，書上沒有數。這裡把帳算清楚，因為這是理解整台譜儀為什麼長那個樣子的起點。
- 第一步：數未知數。 一次散射事件裡，我們想知道的是聲子的波向量 K（三維，3 個分量）與它的頻率 ω（1 個數）。未知數共 4 個。
- 第二步：數方程式。 (33) 是向量式，給 3 條純量方程；(34) 是純量式，給 1 條。方程式共 4 條。4 = 4，所以只要 k 與 k′ 都知道，(K, ω) 就被唯一決定（G 的選法只是把 K 平移一個倒晶格向量，不算新的自由度）。
- 第三步：數「可以量的量」。 k 是實驗者設定的（入射單色化之後，方向與大小都已知，3 個數）；k′ 要量（散射方向 2 個角度 + 能量 1 個數，共 3 個數）。所以每一次「設定 k、量 k′」都換來一組 (K, ω)。
- 第四步：把 (33) 整理成可用的形式。 定義散射向量 Q ≡ k − k′（這是實驗室裡真正在掃的東西）。由 (33)，±K = k − k′ + G = Q + G。因為 G 任意且 −G 也是倒晶格向量，這等價於標準寫法 Q = G ± K。讀法：你在倒空間裡選定一個布拉格點 G，聲子就住在它周圍 ±K 的地方。
- 第五步：把 (34) 整理成可用的形式。 ℏω = ±(ℏ²k²/2M_n − ℏ²k′²/2M_n) = ±(E − E′)。上號：ℏω = E − E′ > 0（中子失能，產生聲子）；下號：ℏω = E′ − E > 0（中子得能，吸收聲子）。兩種情況下 ℏω 都是正的，這與 p.100 的宣告一致：「the energy of a phonon must be positive, so it is conventional and suitable to view ω as positive.」
- 第六步：關鍵的一步——「決定」不等於「量到」。 上面的 4 = 4 只告訴你：給定 (k, k′)，運動學上對應唯一一組 (K, ω)。但散射不一定發生。只有當這組 (K, ω) 恰好滿足晶體真正的色散關係 ω = ω_j(K)（j 是支標號）時，才有可觀測的散射強度。所以實驗做的事是：讓 (k, k′) 掃過一條軌跡，這條軌跡在 (K, ω) 空間裡畫出一條線，線與色散曲面的交點就是你看到的峰。整個實驗就是在數據裡找這些交點，一個點一個點把 ω_j(K) 拼出來。Fig. 11 上的每一顆黑點或白圈，都是這樣一次掃描的一個峰。
- 第七步：為什麼非得同時有 (33) 與 (34)。 只有 (33)（像 X 光繞射那樣只量方向）你只能知道 K，不知道 ω；只有 (34)（只量能量、不管方向）你只能知道 ω，不知道它屬於哪個 K。色散關係是 ω 與 K 的對應，所以兩者缺一不可。這正是「非彈性」散射相對於第 2 章「彈性」散射的全部增量：多量一個能量，就從結構資訊升級成動力學資訊。
散射向量與聲子波向量的關係（由 (33) 整理，書上未寫成這個形式）
Q ≡ k − k′ , Q = G ± K , K = ±(Q + G)|_{摺回第一區}
p.101Q：散射向量（實驗室裡真正掃描的倒空間座標，也常記作動量轉移 ℏQ）；G：倒晶格向量；K：聲子波向量，依原文約定選 G 使其落在第一布里淵區內。
為什麼可以、也應該把 K 摺回第一區我的判讀
原文那句「For a phonon we choose G such that K lies in the first Brillouin zone.」讀起來像一句約定俗成的技術細節，其實它是把本章 p.93 的 First Brillouin Zone 一節整個引用進來。
回想 p.93 的論證：位移比只依賴 exp(iKa)（式 (10)），而把 K 換成 K′ = K − 2πn/a 之後，exp(iK′a) = exp(iKa)exp(−i2πn) = exp(iKa)（式 (11)）。原文的結論是：「Thus the displacement can always be described by a wavevector within the first zone.」也就是說，區外的 K 不描述任何新的晶格運動。既然如此，聲子的 K 就沒有理由寫在區外。
這件事在 (33) 裡的角色是這樣的：實驗量到的 Q = k − k′ 通常很大，往往跨好幾個布里淵區（用剛剛算的熱中子例子，k = 3.49 Å⁻¹，而鈉沿 [001] 的區邊界只到 2π/a ≈ 1.49 Å⁻¹，所以 |Q| 可以輕易到 2–3 個區邊界那麼遠）。你不能說「我量到了一個波向量 3 Å⁻¹ 的聲子」，因為那個聲子跟某個第一區內的聲子是同一個模式。(33) 裡那個「任意的 G」就是讓你做這件減法的工具。
這裡有一個學生常犯的錯：以為 G 是被實驗「決定」的。不是。G 是你選的記帳方式。物理只認 K 摺回第一區之後的那個值與對應的 ω。同一個聲子可以說成「Q = G₁ + K」也可以說成「Q = G₂ − K′」，只要摺回去是同一個模式，就是同一件事。這個「記帳自由度」跟第 9 章 reduced zone scheme / periodic zone scheme（書本 p.223, p.225）是完全同一個概念，只是那裡換成電子的 Bloch 波向量。第 5 章的 Umklapp 過程（書本 p.125）也是同一個 G 在起作用，只是那裡的 G ≠ 0 會改變聲子流的方向，造成熱阻。
量級估計：為什麼偏偏是「熱」中子（本節為原文未寫的補充）補推導
Kittel 在這一段沒有做任何數值估計，但不做這個估計就沒辦法回答一個很自然的問題：為什麼是中子，而不是用第 2 章已經很熟的 X 光？第 2 章其實給過線索——書本 p.24 的 Fig. 1 標題是「Wavelength versus particle energy, for photons, neutrons, and electrons」，橫軸同時印了三條刻度：「Photon energy, keV」「Neutron energy, 0.01 eV」「Electron energy, 100 eV」。同一個波長 ~1 Å，三種粒子的能量差了好幾個數量級。下面把這件事算成數字。
判準有兩條，而且必須同時滿足：
（一）波長要匹配布里淵區的尺度，否則 K 掃不到區邊界；
（二）能量要匹配聲子的能量尺度，否則 (34) 裡的 ℏω 相對於中子動能是個可以忽略的小量，量不出來。
- 第一步：中子的能量—波長關係。 由 (34) 的動能式與 k = 2π/λ：E = ℏ²k²/2M_n = h²/(2M_n λ²)。代 h = 6.626 × 10⁻³⁴ J·s、M_n = 1.675 × 10⁻²⁷ kg、λ = 1.8 Å = 1.8 × 10⁻¹⁰ m。
- 第二步：算出來。 h² = 4.390 × 10⁻⁶⁷；2M_nλ² = 2 × 1.675 × 10⁻²⁷ × (1.8 × 10⁻¹⁰)² = 1.085 × 10⁻⁴⁶。相除得 E = 4.05 × 10⁻²¹ J = 25.2 meV。換算成溫度 T = E/k_B = 4.05 × 10⁻²¹ / 1.381 × 10⁻²³ = 293 K。這就是「熱中子」名稱的由來：與室溫熱浴平衡的中子，波長恰好是 Å 的量級。反應器裡的慢化劑（水、重水、石墨）做的就是把分裂出來的 MeV 快中子慢化到這個能量。
- 第三步：檢查判準（一）。 k = 2π/1.8 Å = 3.49 Å⁻¹。鈉是體心立方，晶格常數約 4.23 Å（此值為外部資料，非本書抽取所得），沿 [001] 的區邊界在 2π/a ≈ 1.49 Å⁻¹。所以 k/K_邊界 ≈ 2.3——中子的波向量是區邊界的兩三倍。這正是需要的：太小則 Q = k − k′ 永遠掃不到區邊界（|Q| ≤ 2k），太大則所有聲子都擠在極小的散射角裡分不開。✔
- 第四步：檢查判準（二）。 由 Fig. 11 目視，鈉的最高聲子頻率 ν ≈ 3.8 × 10¹² Hz，ℏω = hν = 6.626 × 10⁻³⁴ × 3.8 × 10¹² = 2.52 × 10⁻²¹ J = 15.7 meV。與中子動能 25.2 meV 相比，ℏω/E ≈ 0.62。中子散射一次就掉了六成的能量，這是用最粗糙的飛行時間測量都看得到的變化。✔
- 第五步：同一個波長的 X 光呢。 E = hc/λ = (6.626 × 10⁻³⁴ × 3.00 × 10⁸)/(1.8 × 10⁻¹⁰) = 1.10 × 10⁻¹⁵ J = 6.89 keV。判準（一）當然滿足（波長一樣）。但判準（二）慘敗：ℏω/E = 15.7 meV / 6890 eV = 2.3 × 10⁻⁶。要看到聲子，你的能量解析度必須好到 ΔE/E ~ 10⁻⁶ 到 10⁻⁷，也就是要在 7 keV 的光子上分辨 1 meV 的變化。這比 X 光管或一般單光儀能做到的高出四到五個數量級。✘
- 第六步：結論的正確說法。 不是「X 光不能散射聲子」——它當然可以，非彈性 X 光散射每天都在發生。問題是你分不出來：散射前後的光子能量幾乎相同，能量損失被埋在譜儀的解析度裡。中子的優勢完全來自一個巧合：M_n 這個質量，使得「波長 ~ Å」與「能量 ~ 10 meV」這兩件事同時成立。光子因為 E = hc/λ（能量與波長成反比、且係數是 c）做不到；電子因為質量太小（同樣 1 Å 對應約 150 eV，見第 2 章 Fig. 1 的「Electron energy, 100 eV」刻度）也做不到，何況電子帶電、穿透深度只有奈米量級，量到的是表面而非塊材。
中子的能量—波長關係（由 (34) 與 de Broglie 關係整理）
E = ℏ²k²/2M_n = h²/(2M_n λ²) , λ[Å] ≈ 9.045 / √(E[meV])
p.102E：中子動能；λ：中子的 de Broglie 波長；k = 2π/λ；M_n = 1.675 × 10⁻²⁷ kg 為中子質量；h = 6.626 × 10⁻³⁴ J·s。
Fig. 11 是什麼：原文說了什麼、圖上真正有什麼
因為這份網頁不放圖，以下用文字把 Fig. 11 完整描述出來。圖形內容為對書本 p.101 原圖放大後的目視判讀；圖說文字則是原文逐字。
原文圖說（p.101，逐字）：「Figure 11 The dispersion curves of sodium for phonons propagating in the [001], [110], and [111] directions at 90 K, as determined by inelastic scattering of neutrons, by Woods, Brockhouse, March and Bowers.」
圖的版面（目視判讀）：三個並排的矩形面板，共用同一條縱軸與同一條橫軸標題。
- 縱軸：「Frequency, in 10¹² Hz」，刻度標 0、2、4（另有 1 與 3 的無標記刻度），所以整張圖的頻率範圍是 0 到 4 × 10¹² Hz。
- 橫軸標題（印在三個面板下方正中）：「Wavevector, in units 2π/a」。
- 三個面板由左到右分別標 [001]、[110]、[111]，各自的橫軸都從 0 起算，右端終點分別印為 100、½½0、½½½（皆以 2π/a 為單位）。
- 圖例（印在左邊 [001] 面板內）：「● Longitudinal」「○ Transverse」，即實心圓為縱支、空心圓為橫支。
注意抽取陷阱：這份 PDF 的文字層會把分數拆成浮出的碎片，Fig. 11 的三個橫軸終點在文字層裡變成一串「00 1 20100 1 2 1 2 1 2 1 2」，完全無法照抄。上面的 100 / ½½0 / ½½½ 是把該頁以 350–900 dpi 算繪成影像後目視還原的結果。
Fig. 11 讀圖訓練（一）：縱支為什麼比橫支高、橫支為什麼在 [110] 分裂成兩條我的判讀
這張圖是全章第一次出現「真實材料、真實數據」，值得慢慢讀。以下每一項都標明是目視判讀還是可由書中他處推得。
（1）三個面板各有幾條曲線。目視：[001] 面板上實心圓與空心圓各一串，而且兩串幾乎貼在一起；[110] 面板上實心圓一串、空心圓明顯分成上下兩串；[111] 面板上實心圓一串、空心圓一串。
這不是隨機的。第 3 章 p.83 的 Figure 20 圖說逐字寫著：「The two transverse modes are degenerate for propagation in the [100] and [111] directions.」——立方晶體沿 [100] 與 [111] 傳播時，兩個橫模簡併；沿 [110] 則不簡併。Fig. 20 還列出三個方向的有效彈性常數：[100] 的 L 是 C₁₁、T 是 C₄₄；[110] 的 L 是 ½(C₁₁+C₁₂+2C₄₄)、T₁ 是 C₄₄、T₂ 是 ½(C₁₁−C₁₂)；[111] 的 L 是 ⅓(C₁₁+2C₁₂+4C₄₄)、T 是 ⅓(C₁₁−C₁₂+C₄₄)。
所以 Fig. 11 的 [110] 面板出現兩條橫支，不是數據雜亂，而是對稱性的直接後果——T₁ 的位移沿 [001]，T₂ 的位移沿 [11̄0]（第 3 章 p.83 明寫第二根的位移「is along [11̄0] and perpendicular to the K vector」；校註：這個負指標在 PDF 文字層被抽成「[11 -0]」，上橫線退化成一個連字號，正確寫法是 [11̄0]）。
（2）縱支為什麼總是比橫支高。物理上，縱波壓縮鍵長，橫波只剪切鍵角；對中心力主導的固體，伸縮硬度遠大於剪切硬度，所以 L 的恢復力常數大、頻率高。書上的定量版本就在 Fig. 20 的那些組合裡。用書本 p.84 Table 12 給的鈉室溫值 C₁₁ = 0.073、C₁₂ = 0.062、C₄₄ = 0.042（單位 10¹² dyne/cm²）代進去：
| 方向 | 縱支有效常數 | 橫支有效常數 | 預測 v_L/v_T |
|---|---|---|---|
| [100] | C₁₁ = 0.073 | C₄₄ = 0.042 | √(0.073/0.042) = 1.32 |
| [110] | ½(C₁₁+C₁₂+2C₄₄) = 0.1095 | T₁: C₄₄ = 0.042；T₂: ½(C₁₁−C₁₂) = 0.0055 | 對 T₁ 為 1.61；對 T₂ 為 4.46 |
| [111] | ⅓(C₁₁+2C₁₂+4C₄₄) = 0.1217 | ⅓(C₁₁−C₁₂+C₄₄) = 0.0177 | 2.62 |
這張表把 Fig. 11 的形狀解釋得乾乾淨淨：[001] 的 L 與 T 只差 32%，所以兩串點幾乎重疊；[110] 的 T₂ 比 L 慢四倍以上，所以最下面那串空心圓明顯壓得很低。
（3）[110] 的下橫支為什麼特別軟。因為 ½(C₁₁ − C₁₂) 這個組合在鈉身上幾乎歸零：½(0.073 − 0.062) = 0.0055，只有 C₄₄ = 0.042 的八分之一。這是體心立方鹼金屬的著名特徵（彈性異向比 C₄₄/½(C₁₁−C₁₂) ≈ 7.6，可直接由 Table 12 算出）。它的物理意義是：鈉對「把立方體剪成長方體」這種變形幾乎沒有抵抗力，因為鹼金屬的鍵結來自幾乎均勻的傳導電子海（第 3 章 Metals 一節，p.69），沒有方向性的共價骨架。
（4）區邊界附近的頻率（目視讀值，僅供數量級參考）。[001] 面板右端（H 點）：L 與 T 都接近 3.6 × 10¹² Hz，幾乎重合。[110] 面板右端（N 點）：L 約 3.8（全圖最高）、上橫支約 2.5、下橫支約 1.0。[111] 面板：L 在橫軸約 0.3 處衝到約 3.5 的高點後反而向下彎，到右端（P 點）降到約 3.0；橫支在右端約 2.9，與 L 幾乎簡併。
（5）[111] 的 L 支為什麼會「先升後降」。單原子近鄰鏈的 (9) ω = (4C/M)^{1/2}|sin ½Ka| 是單調上升到區邊界的；Fig. 11 的 [111] 縱支卻在區內就出現極大再回落。這正是本章 p.94–95 所說的長程力的證據：原文寫「In metals the effective forces may be of quite long range and are carried from ion to ion through the conduction electron sea. Interactions have been found between planes of atoms separated by as many as 20 planes.」有了 p > 1 的 C_p，(16a) 的 ω² = (2/M)Σ_{p>0} C_p(1 − cos pKa) 就可以產生各種非單調的形狀。想反推 C_p，就用 (17) 的傅立葉反轉式去積分 Fig. 11 的曲線——這就是 p.94 那一節標題 Derivation of Force Constants from Experiment 的實際做法，而 Fig. 11 正是它的原料。
Fig. 11 讀圖訓練（二）：為什麼三個橫軸終點剛好都是區邊界補推導
圖上三個面板的右端分別是 (1,0,0)、(½,½,0)、(½,½,½)（單位 2π/a）。這三個數字看起來毫無規律，但它們其實是同一件事：鈉是體心立方，這三點正好落在體心立方第一布里淵區的表面上。書上沒有做這個檢查，這裡補出來——會讀這個，以後看任何一張色散圖的橫軸都不會迷路。
所需材料全在第 2 章 p.36–37：體心立方的倒晶格是面心立方，最短的 12 個倒晶格向量是 (2π/a)(±ŷ±ẑ)、(2π/a)(±x̂±ẑ)、(2π/a)(±x̂±ŷ)（式 (33) 於第 2 章），而第一區「is bounded by the planes normal to the 12 vectors of Eq. (33) at their midpoints」，形狀是「a regular 12-faced solid, a rhombic dodecahedron」（p.37，並見 p.43 的摘要表）。
區邊界條件因此是：對某個最短的 G，k·Ĝ = |G|/2，等價於 2k·G = G²。
- 第一步：把 G 與 k 都寫成 2π/a 的倍數。 令 G = (2π/a)g、k = (2π/a)q，則條件 2k·G = G² 化為 2q·g = g²。所有最短的 g 都是 {110} 型，g² = 2。所以條件簡化成 q·g = 1。
- 第二步：檢查 [001] 的終點 q = (1,0,0)。 取 g = (0,1,1)：q·g = 0 ✘。取 g = (1,1,0)：q·g = 1 ✔。所以 (1,0,0) 在垂直於 (1,1,0) 的那個面上——確實在區邊界。這就是體心立方的 H 點。（注意：如果晶格是簡單立方，區邊界會落在 (½,0,0)；(1,0,0) 之所以還在第一區的表面上，正是因為體心立方的區比簡單立方大一倍。）
- 第三步：檢查 [110] 的終點 q = (½,½,0)。 取 g = (1,1,0)：q·g = ½ + ½ = 1 ✔。在區邊界，而且這是垂直於 (110) 那面的面心，即體心立方的 N 點。N 點是第一區上離原點最近的一批點（|q| = 0.707），所以 [110] 方向的色散曲線最早到頂。
- 第四步：檢查 [111] 的終點 q = (½,½,½)。 取 g = (0,1,1)：q·g = ½ + ½ = 1 ✔。在區邊界。這是三個菱形面交會的頂角，即體心立方的 P 點（|q| = 0.866）。
- 第五步：讀出物理。 三個面板的橫軸長度不一樣——[001] 走到 |q| = 1，[110] 只走到 0.707，[111] 走到 0.866。所以三張圖不能直接把橫軸並排比較斜率；要比長波聲速，必須回到 ω/K 而不是 ω/（格子刻度）。這是初學者看多方向色散圖最常犯的錯。
- 第六步：驗一下群速度。 本章 p.94 的 (13) v_g = grad_K ω(K)。在區邊界上，若該點是由對稱性保證的極值（例如 H 點沿 [001]），則群速度為零、波是駐波（p.93–94 的 (12) 與 Bragg 反射論證）。Fig. 11 的 [001] 面板右端曲線確實明顯轉平 ✔；但 [111] 面板的 L 支在 P 點仍在下降、沒有轉平，這並不矛盾。判準不是「頂角還是面心」，而是 2K 本身是不是倒晶格向量。把本區塊開頭列的那 12 個最短倒晶格向量 (2π/a)(±ŷ±ẑ)、(2π/a)(±x̂±ẑ)、(2π/a)(±x̂±ŷ) 當生成元，以 2π/a 為單位時，倒晶格向量恰好是「三個分量皆為整數、且三者之和為偶數」的那些點（這就是體心立方的倒晶格＝面心立方）。H 點：2H = (2π/a)(2,0,0)，分量和 2 為偶數 ⇒ 是倒晶格向量 G。於是由 ω(−K) = ω(K)（力常數為實數，時間反轉不變）與 ω(K + G) = ω(K)（週期性）得 ω(H + q) = ω(−H − q) = ω(H − q)：ω 在 H 點左右對稱，斜率必須為零。N 點：2N = (2π/a)(1,1,0)，和為 2，同樣是倒晶格向量，同樣被保證零斜率。P 點：2P = (2π/a)(1,1,1)，分量和 3 為奇數 ⇒ 不是倒晶格向量，這個對稱性論證失效，所以沒有任何理由要求 ∂ω/∂K 在 P 點為零。（特別注意：H 點本身也是頂角——四個菱形面交會處，不是面心——所以「頂角 vs 面心」不能當判準，真正在起作用的是 2K = G。）這個細微差別書上沒講，但它是讀圖時真正會被考的東西。
Fig. 12：三軸中子譜儀的三個軸，各對應 (33)(34) 裡的哪個量延伸知識
先說清楚證據等級。Kittel 在這一整段裡對 Fig. 12 只寫了一行圖說，逐字是：「Figure 12 A triple axis neutron spectrometer at Brookhaven. (Courtesy of B. H. Grier.)」正文完全沒有解釋「三軸」是哪三個軸、機器怎麼運作。全書搜尋「triple axis」也只有這一處。所以以下內容是 Kittel 完全沒有的外部知識，標為 ext；但它是把 (33)(34) 從紙上變成儀器的唯一橋樑，不講的話這張照片對讀者毫無資訊量。
照片上看得到什麼（對書本 p.101 原圖算繪後目視判讀）：一張黑白照片。畫面右側是一座巨大的圓柱形屏蔽鼓，上方接著一圈環狀管線與閥件；中央偏左是一具被厚重屏蔽包覆的水平臂，末端是一個圓柱形的探測器屏蔽桶；左側站著一位穿毛衣的人在看手上的資料，旁邊是一座控制／電子儀器機櫃；地板上可以看到弧形的軌道，供臂繞著轉。照片本身沒有任何可抽取的文字或數據。
三個軸分別是什麼，以及它們在 (33)(34) 裡的角色：
| 軸 | 硬體 | 決定什麼 | 對應 (33)/(34) 的哪一項 |
|---|---|---|---|
| 第一軸 | 單色化晶體（照片右側屏蔽鼓出口處）繞自身轉動，並帶動下游整個臂繞它公轉 | 用布拉格反射從反應器的白光中子束中選出單一波長 λ_i，即固定 k 的大小 | (33) 與 (34) 裡的 k |
| 第二軸 | 樣品台（中央），可繞垂直軸轉（樣品定向）並帶動探測器臂公轉（散射角 2θ） | 決定樣品晶軸相對入射束的取向，以及散射方向，即固定 k′ 的方向 | (33) 裡 Q = k − k′ 的方向，也就是你在倒空間裡選哪個 G、沿哪個晶軸掃 K |
| 第三軸 | 分析晶體（探測器臂上），用另一次布拉格反射只讓單一波長 λ_f 進入探測器 | 量出散射中子的 k′ 大小，也就是它的能量 | (34) 裡的 ℏ²k′²/2M_n，也就是能量增益／損失 |
把這張表和前面的第七步對照，整件事就通了：第一軸釘住 k，第二軸釘住 k′ 的方向，第三軸量 k′ 的大小；三者合起來就把 (K, ω) 完全決定。Kittel 那句「find the energy gain or loss of the scattered neutrons as a function of the scattering direction k – k′」——「energy gain or loss」是第三軸的工作，「scattering direction k – k′」是第二軸的工作，而「入射單色化」是第一軸的前提。一句話裡的三個要素，正好就是三個軸。
為什麼要三個軸而不是兩個。如果沒有第三軸（分析晶體），你量到的是所有能量的散射中子總和，(34) 裡的 ℏω 就被積分掉了——那台機器就退化成第 2 章的中子繞射儀，只能量結構。第三軸是「非彈性」三個字的硬體化身。
實務上怎麼掃。兩種標準掃描：固定 Q 掃 ℏω（constant-Q scan，在色散曲線上沿垂直方向切，得到一個能量峰），或固定 ℏω 掃 Q（constant-E scan，水平方向切）。Fig. 11 上那六十來顆點，就是六十來次這樣的掃描各自找到的峰位。這解釋了為什麼中子色散曲線是離散的資料點而不是連續曲線——每一顆點都是一次數小時等級的計數實驗。
SUMMARY 逐條回扣（一）：聲子能量與選擇律
p.102 的 SUMMARY 只有四條，短得像順口溜，但每一條都精準對應前面一整節。逐條核對一次，順便把「同一條物理寫成不同式子」的混淆解掉。
第一條（逐字）：「The quantum unit of a crystal vibration is a phonon. If the angular frequency is ω, the energy of the phonon is ℏω.」
回扣 QUANTIZATION OF ELASTIC WAVES（p.99）。那一節的 (27) 寫的是 ε = (n + ½)ℏω，不是 ℏω。兩者不衝突：(27) 是整個模式被激發到量子數 n 時的能量，而摘要說的是「a phonon」——多一個聲子，模式能量增加 ℏω。摘要為了簡潔省掉了零點能 ½ℏω，這是對的，因為零點能與聲子數無關、不參與 (34) 的能量收支。原文 p.99 對此說得很清楚：「The term ½ℏω is the zero point energy of the mode.」（依前述體例，½ℏω 是從文字層浮出的「1 2 ℏω」還原回句中原位。）
第二條（逐字）：「When a phonon of wavevector K is created by the inelastic scattering of a photon or neutron from wavevector k to k′, the wavevector selection rule that governs the process is k = k′ + K + G, where G is a reciprocal lattice vector.」
這一條最容易讓人困惑，因為它與正文的 (31) 和 (33) 長得不一樣。詳見下一個區塊的對照。另外注意摘要說的是「a photon or neutron」——它把 (31)（光子）與 (33)（中子）合併成一條通則，因為選擇律的來源是晶格的離散平移對稱性，跟探針是什麼完全無關。
SUMMARY 的波向量選擇律
k = k′ + K + G
p.102k：入射粒子（光子或中子）的波向量；k′：散射後的波向量；K：被產生的聲子波向量；G：倒晶格向量。原文限定情境為「a phonon of wavevector K is created」，也就是只寫了產生的情形。
把 (30)(31)(32)(33) 與摘要式擺在一起：它們其實是同一條補推導
很多人第一次讀這一章，會以為 Kittel 給了四、五條不同的選擇律。其實只有一條，寫了五次。把它們全部移項到「k 在左、其餘在右」的標準型，一切就清楚了。書上沒有做這張對照表。
- 原文列出的五種寫法（逐字照抄式子本身）：(30) k′ = k + G（第 2 章的彈性散射，p.100）；(31) k′ + K = k + G（光子非彈性散射、產生聲子，p.100）；(32) k′ = k + K + G（吸收聲子，p.100）；(33) k + G = k′ ± K（中子散射，產生取 +、吸收取 −，p.101）；SUMMARY 式 k = k′ + K + G（產生，p.102）。
- 統一移項成 k = … 的形式。(30)：k = k′ − G。(31)：k = k′ + K − G。(32)：k = k′ − K − G。(33) 上號：k = k′ + K − G。(33) 下號：k = k′ − K − G。摘要式：k = k′ + K + G。
- 用 G → −G 消除表面差異。 因為 G 是「任一」倒晶格向量，−G 也是，所以 −G 與 +G 在陳述上完全等價。於是：(30) ⇔ k = k′ + G（純彈性）；(31) ⇔ (33) 上號 ⇔ 摘要式 ⇔ k = k′ + K + G（產生聲子）；(32) ⇔ (33) 下號 ⇔ k = k′ − K + G（吸收聲子）。
- 結論（一句話）：全章只有兩條選擇律——產生與吸收——再加上 K = 0 的彈性特例。(33) 用一個 ± 號把兩條寫成一條；SUMMARY 只寫產生那一支。任何看起來不同的地方，都只是 G 的正負號約定。
- 為什麼值得花這個力氣。 因為到了第 5 章談 Umklapp、第 14 章談電子—聲子交互作用、第 15 章談 Raman 散射（書本 p.444），同一條式子還會再換好幾種寫法出現。如果你把「移項」誤認成「不同的物理」，每換一次符號就要重學一次。
- 一個常見的符號陷阱。 有些教科書（與絕大多數中子散射文獻）用 Q = k_i − k_f 並寫成 Q = G ± K，有些則用 Q = k_f − k_i（差一個負號）。Kittel 這裡兩種都沒有明確定義 Q，他直接寫 (33)。跨書比對時一定要先確認對方的 Q 定義，否則產生與吸收會整個對調。
SUMMARY 逐條回扣（二）：第一區的充分性與支數計數
第三條（逐字）：「All elastic waves can be described by wavevectors that lie within the first Brillouin zone in reciprocal space.」
回扣 First Brillouin Zone（p.93）。那一節的論證鏈是：(10) 位移比 u_{s+1}/u_s = exp(iKa) → 只有相位 Ka ∈ (−π, π] 是獨立的 → (11) 減掉 2πn/a 之後位移比不變 → 「Thus the displacement can always be described by a wavevector within the first zone.」原文對區外的 K 說得很不客氣：「Values of K outside of the first Brillouin zone (Fig. 5) merely reproduce lattice motions described by values within the limits ±π/a.」Fig. 5 的圖說也強化了這一點：「The wave represented by the solid curve conveys no information not given by the dashed curve. Only wavelengths longer than 2a are needed to represent the motion.」
這一條為什麼放在摘要裡：因為它是 (33) 那句「we choose G such that K lies in the first Brillouin zone」的理論依據。摘要把「可以摺」與「一定摺得進去」講成一件事。
一個容易被忽略的前提：這條敘述對晶格成立，是因為位移只定義在原子上（離散取樣）。原文在 p.94 特別對比過：X 光可以有 n > 1 的高階布拉格反射，「because the amplitude of the electromagnetic wave has a meaning in the space between atoms, but the displacement amplitude of an elastic wave usually has a meaning only at the atoms themselves.」「只在原子上有意義」＝取樣定理＝第一區就夠用。這是整個布里淵區概念在聲子問題上的物理根源。
第四條（逐字）：「If there are p atoms in the primitive cell, the phonon dispersion relation will have 3 acoustical phonon branches and 3p − 3 optical phonon branches.」
回扣 TWO ATOMS PER PRIMITIVE BASIS（p.95–96）。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.」p.96 給的計數理由值得複述：p 個原子、N 個原胞 → pN 個原子 → 每個原子三個自由度 → 共 3pN 個自由度；而「The number of allowed K values in a single branch is just N for one Brillouin zone.」（並附註腳說明第 5 章會用週期性邊界條件證明這一點）。3pN 個自由度除以每支 N 個 K 值，得到 3p 支。其中 LA 加兩個 TA 共 3 支佔掉 3N，剩下 (3p − 3)N 由光學支承擔。
這條計數律的實用價值：拿到任何一張色散圖，先數有幾條曲線，就能反推原胞裡有幾個原子。Fig. 8a 的鍺與 Fig. 8b 的 KBr 都是 p = 2，所以「six branches: one LA, one LO, two TA, and two TO」（p.95）。而 Fig. 11 的鈉是 p = 1（體心立方的原胞只有一個原子），所以只有三支聲學支、沒有任何光學支——這正是 Fig. 11 每個面板上只有 L 與 T（[110] 是 L、T₁、T₂）而沒有更高的第四條曲線的原因。用第四條摘要去核對第一張真實數據圖，這就是摘要該有的用法。
習題導讀：七題各補回什麼，為什麼只講透三題我的判讀
第 4 章的七道習題印在書本 p.102–103（第 7 題的題幹延續到 p.104），而且它們不是複習題，是補課題。正文為了篇幅省略的幾塊物理，恰好各由一題補回：
| 題號 | 標題 | 補回正文省略的哪一塊 |
|---|---|---|
| 1 | Monatomic linear lattice | (9) 的能量學：驗證簡諧振子的動能＝位能時均分，並反過來驗證 (9) |
| 2 | Continuum wave equation | (2) → 連續介質波動方程的極限，也就是 (15) 的另一種說法 |
| 3 | Basis of two unlike atoms | (25) 兩個區邊界根的物理圖像：兩個子晶格解耦，一動一靜 |
| 4 | Kohn anomaly | (17) 長程力的直接後果：色散曲線出現垂直切線的扭結，接費米面 |
| 5 | Diatomic chain | 等質量、力常數交替（C 與 10C）的鏈：說明「雙原子基底」不必是兩種原子 |
| 6 | Atomic vibrations in a metal | 金屬中離子振動的 Coulomb 估算，得 ω = (e²/MR³)^{1/2} |
| 7 | Soft phonon modes | p.100「ω² 為負則 ω 為虛數」那句話的具體實現：庫侖長程力使晶格失穩 |
第 7 題還帶一個原文腳註（p.103 頁底，逐字）：「*This problem is rather difficult.」——全章唯一被加星號的題目。
下面挑 3、4、7 講透，理由是這三題各自補回一塊正文真的沒給、而且後面章節會用到的東西。第 1、2、5、6 題會在最後用一個區塊速記帶過。
Problem 3（Basis of two unlike atoms）：題目在問什麼
原文逐字（p.103）：「3. Basis of two unlike atoms. For the problem treated by (18) to (26), find the amplitude ratios u/v for the two branches at K_max = π/a. Show that at this value of K the two lattices act as if decoupled: one lattice remains at rest while the other lattice moves.」
為什麼這題重要。正文 p.98 在區邊界只給了兩個數字，式 (25)：ω² = 2C/M₁；ω² = 2C/M₂。給完就跳過去了。學生看到這兩個數只會記得「上下兩支在區邊界的頻率」，但完全不知道那兩個模式長什麼樣子。而 p.99 接著說的「There is a frequency gap at the boundary K_max = ±π/a of the first Brillouin zone.」——這個間隙為什麼存在、間隙上下兩端各是什麼運動，正文一個字都沒交代。
這題就是把 (25) 從「兩個數字」變成「兩張圖像」。而且它的答案會直接解釋兩件事：為什麼有間隙，以及為什麼 M₁ = M₂ 時間隙會消失（後者接回第一區折疊的概念）。
Problem 3 完整解：區邊界的兩個子晶格為什麼解耦補推導
所需材料是正文 p.97 的 (20)（把行波試解 (19) 代入運動方程 (18) 之後的聯立式）：
−ω²M₁u = Cv[1 + exp(−iKa)] − 2Cu ；
−ω²M₂v = Cu[exp(iKa) + 1] − 2Cv 。
以及 p.98 的行列式條件 (21) 與其展開 (22) M₁M₂ω⁴ − 2C(M₁+M₂)ω² + 2C²(1 − cos Ka) = 0。書上完全沒有做下面這幾步。
- 第一步：在 Ka = π 把耦合項算出來。 exp(iπ) = exp(−iπ) = −1，所以 1 + exp(−iKa) = 1 − 1 = 0，同樣 exp(iKa) + 1 = 0。(20) 的兩式右邊的第一項（唯一把 u 與 v 綁在一起的項）整個消失。
- 第二步：兩式因此完全脫鉤。 剩下 −ω²M₁u = −2Cu 與 −ω²M₂v = −2Cv，即 (2C − M₁ω²)u = 0 與 (2C − M₂ω²)v = 0。這是兩個各自獨立的單自由度方程。
- 第三步：確認這與 (25) 一致。 行列式 (21) 在 Ka = π 時非對角元為 −C[1 + exp(iπ)] = 0，所以行列式退化成 (2C − M₁ω²)(2C − M₂ω²) = 0，兩根正是 ω² = 2C/M₁ 與 ω² = 2C/M₂ ✔ 與 (25) 完全相符。
- 第四步：取第一個根 ω² = 2C/M₁，求 u/v。 代入第二式：(2C − M₂ · 2C/M₁)v = 2C(1 − M₂/M₁)v = 0。只要 M₁ ≠ M₂（題目說的 unlike atoms），括號不為零，故 v = 0。而第一式 (2C − M₁ · 2C/M₁)u = 0 恆成立，對 u 沒有任何限制。所以 u/v = ∞：M₁ 子晶格在動，M₂ 子晶格完全靜止。
- 第五步：取第二個根 ω² = 2C/M₂。 對稱地，第一式給 (2C − M₁ · 2C/M₂)u = 2C(1 − M₁/M₂)u = 0 ⇒ u = 0，第二式對 v 無限制。所以 u/v = 0：M₂ 子晶格在動，M₁ 子晶格完全靜止。
- 第六步：對上圖 7。 Fig. 7 的圖形內標註，以及 p.98 的正文（逐字：「The dependence of ω on K is shown in Fig. 7 for M₁ > M₂.」），都指明 M₁ > M₂（圖說本身沒有這句話）。既然 ω² = 2C/M₁ < 2C/M₂，較低的那個根（＝聲學支在區邊界的端點）是「重原子動、輕原子不動」；較高的那個根（＝光學支在區邊界的端點）是「輕原子動、重原子不動」。這個對應關係很值得記：光學支的高頻來自輕原子。
- 第七步：不用代數、用力學直接看懂為什麼。 在 Ka = π，位移是 (12) 式的駐波 u_s = u(−1)^s、v_s = v(−1)^s。算 M₁ 原子受的力（(18) 第一式右邊）：C(v_s + v_{s−1} − 2u_s) = C[v(−1)^s + v(−1)^{s−1} − 2u(−1)^s] = C[v(−1)^s − v(−1)^s − 2u(−1)^s] = −2C u_s。v 的貢獻精確抵消——因為一個 M₁ 原子左右兩側的 M₂ 鄰居位移方向正好相反，它們拉扯 M₁ 的力大小相等方向相反。於是每個 M₁ 原子看起來就像「被兩根彈簧綁在不動的牆上」的單一質點，頻率當然是 ω² = 2C/M₁。這就是「decoupled」三個字的全部內容。
- 第八步：解釋間隙。 既然區邊界的兩個模式一個只讓 M₁ 動、一個只讓 M₂ 動，而兩者的「等效彈簧常數」都是 2C、質量卻不同，頻率必然是兩個不同的數，中間那一段沒有任何模式——這就是 p.99 說的 frequency gap，原文寫「Wavelike solutions do not exist for certain frequencies, here between (2C/M₁)^{1/2} and (2C/M₂)^{1/2}.」間隙寬度純粹由質量比決定。
- 第九步：M₁ → M₂ 的極限（書上沒講，但這是本題最有價值的延伸）。 若 M₁ = M₂ = M，第四步的括號 (1 − M₂/M₁) = 0，論證失效；兩個根合併成 ω² = 2C/M，間隙關閉。這不是巧合：等質量、等力常數時，真正的原胞只含一個原子、真正的區邊界在 K = 2π/a 而不是 π/a。我們看到的「兩支」其實是同一條單原子鏈色散曲線在 K = π/a 處被摺回來的兩半。所以間隙的存在等價於「基底真的不對稱」。這一點直接預告第 7 章：電子在週期位能中的能隙，也是同一個機制——區邊界的簡併被微擾解除。Problem 5 用「等質量、力常數交替 C 與 10C」把這件事再演一次：那裡質量相同但力常數不同，間隙照樣打開，證明不對稱性來源可以是質量、也可以是鍵強。
區邊界的振幅比（Problem 3 的答案，書上只給頻率不給振幅）
Ka = π 時：ω² = 2C/M₁ ⇒ v = 0（u/v → ∞）；ω² = 2C/M₂ ⇒ u = 0（u/v = 0）
p.103u、v：M₁ 與 M₂ 兩個子晶格的複振幅（試解 (19)）；C：相鄰平面間力常數；a：重複距離（p.97 定義為「the distance between nearest identical planes, not nearest-neighbor planes」）。
Problem 4（Kohn anomaly）：題目在問什麼，以及它在補哪一塊
原文逐字（p.103）：「4. Kohn anomaly. We suppose that the interplanar force constant C_p between planes s and s + p is of the form C_p = A sin pk₀a / pa, where A and k₀ are constants and p runs over all integers. Such a form is expected in metals. Use this and Eq. (16a) to find an expression for ω² and also for ∂ω²/∂K. Prove that ∂ω²/∂K is infinite when K = k₀. Thus a plot of ω² versus K or of ω versus K has a vertical tangent at k₀: there is a kink at k₀ in the phonon dispersion relation ω(K).」
強宣稱查證：「Kohn anomaly」這個字串在全書出現三次，分別是第 4 章章前目錄（p.89）、這道習題（p.103）、以及索引（p.677 指向 p.103）。Kittel 的正文從頭到尾沒有展開討論過 Kohn anomaly——它整個被放進一道習題裡。正文唯一的相關提及在第 9 章 p.242（EXPERIMENTAL METHODS IN FERMI SURFACE STUDIES），逐字：「Further information on the momentum distribution is given by positron annihilation, Compton scattering, and the Kohn effect.」——只是把「the Kohn effect」列為量費米面動量分布的手段之一，一句話帶過，沒有任何解釋。這一句其實正好佐證後面 k₀ ↔ 2k_F 的論證：Kittel 自己把 Kohn 效應歸在「費米面的動量分布」那一類。同樣地，「thermal neutron」這個字串在全書一次都沒出現過。（以上為全書文字層搜尋結果。）
這題在補哪一塊。正文 p.94–95 的 Derivation of Force Constants from Experiment 說了兩件事：一是金屬裡的力可以很長程，「Interactions have been found between planes of atoms separated by as many as 20 planes.」；二是有了 (16a) ω² = (2/M)Σ_{p>0} C_p(1 − cos pKa)，就可以用 (17) C_p = −(Ma/2π)∫dK ω²_K cos pKa 把 C_p 反解出來。
但正文沒有回答一個很自然的問題：長程力會讓色散曲線長成什麼樣子？如果 C_p 只到 p = 1，ω(K) 就是平滑的 |sin ½Ka|；如果 C_p 一路延伸到 p = 20 甚至無窮，會發生什麼事？Problem 4 給的答案非常戲劇化：會出現一個奇異點——在某個特定的 K = k₀，導數發散、曲線出現垂直切線。而這個 k₀ 不是隨便的數，它由傳導電子的費米面尺度決定。
所以這題是本章與第 6、9 章之間唯一的明示橋樑。
Problem 4 假設的長程層間力常數
C_p = A sin(p k₀ a) / (p a)
p.103C_p：相隔 p 層的兩個原子面之間的力常數；A、k₀：常數；a：層間距；p：整數（原文明說「p runs over all integers」）。原文對這個形式的評註只有一句：「Such a form is expected in metals.」
Problem 4 完整解：證明 ∂ω²/∂K 在 K = k₀ 發散，並求出 ω²(K) 的封閉形式補推導
書上只說「Prove that ∂ω²/∂K is infinite when K = k₀」，一步也沒給。下面先做最短的證明（這是 Kittel 想要的），再補上完整的封閉解（這個結果相當驚人，而且書上完全沒有提）。
出發點是正文 p.95 的 (16a)：ω² = (2/M) Σ_{p>0} C_p (1 − cos pKa)。
- 第一步：代入 C_p。 ω² = (2/M) Σ_{p>0} [A sin(pk₀a)/(pa)] (1 − cos pKa) = (2A/Ma) Σ_{p≥1} sin(pk₀a)(1 − cos pKa)/p。
- 第二步：直接對 K 微分（逐項）。 ∂/∂K [1 − cos pKa] = pa sin(pKa)。注意 C_p 前面的 1/(pa) 與這個 pa 精確抵消——這正是題目把 C_p 設成 sin(pk₀a)/(pa) 而不是別的形式的用意。於是 ∂ω²/∂K = (2A/M) Σ_{p≥1} sin(pk₀a) sin(pKa)。
- 第三步：化成餘弦。 用 sin X sin Y = ½[cos(X − Y) − cos(X + Y)]：∂ω²/∂K = (A/M) Σ_{p≥1} [cos(p(k₀ − K)a) − cos(p(k₀ + K)a)]。
- 第四步：看 K = k₀ 發生什麼事。 第一個級數的相位 p(k₀ − K)a 變成 0，於是每一項都等於 cos 0 = 1，級數成為 Σ_{p≥1} 1 = ∞。第二個級數的相位是 2pk₀a，只要 2k₀a 不是 2π 的整數倍，它就是一個有界（Abel 可和，值為 −½）的振盪級數。兩者相減仍然發散。故 ∂ω²/∂K → ∞ 於 K = k₀，證畢。這正是 Kittel 要的「vertical tangent」。
- 第五步（書上沒有）：求 ω² 本身的封閉形式。 用標準鋸齒波傅立葉級數 S(θ) ≡ Σ_{p≥1} sin(pθ)/p = (π − θ)/2（對 0 < θ < 2π，且 S 為奇函數、週期 2π）。先把第一步的和式拆開：Σ sin(pk₀a)(1 − cos pKa)/p = S(k₀a) − Σ sin(pk₀a)cos(pKa)/p，再用 sin X cos Y = ½[sin(X+Y) + sin(X−Y)]，得 Σ sin(pk₀a)cos(pKa)/p = ½[S((k₀+K)a) + S((k₀−K)a)]。
- 第六步：代入鋸齒波公式（取 0 < k₀a < π、0 < Ka < π）。 情形 A：K < k₀，此時 (k₀−K)a > 0，S((k₀−K)a) = (π − k₀a + Ka)/2，而 S((k₀+K)a) = (π − k₀a − Ka)/2。兩者相加除以 2 得 (π − k₀a)/2，與 S(k₀a) 相減恰好為零。所以 ω² = 0。
- 第七步：情形 B：K > k₀。 此時 (k₀−K)a < 0，必須用奇函數性質：S(負值 θ) = −(π + θ)/2。於是 ½[S((k₀+K)a) + S((k₀−K)a)] = ½[(π − k₀a − Ka)/2 − (π + k₀a − Ka)/2] = −k₀a/2。與 S(k₀a) = (π − k₀a)/2 相減得 (π − k₀a)/2 + k₀a/2 = π/2。所以 ω² = (2A/Ma)(π/2) = πA/(Ma)。
- 第八步：結論與數值驗算。 ω²(K) = 0（K < k₀）、πA/(Ma)（k₀ < K < π/a）——是一個完美的階梯函數。以 Abel 求和（乘 r^p、r → 1⁻，取 r = 0.9999995、截到 4 × 10⁶ 項）數值驗算：(k₀a, Ka) = (1.0, 0.5) 得和 0.00000；(1.0, 2.0) 得 1.57080；(2.0, 1.0) 得 0.00000；(2.0, 3.0) 得 1.57080；(0.8, 0.3) 得 0.00000；(0.8, 3.0) 得 1.57080。理論值 0 與 π/2 = 1.57080 ✔ 完全吻合。
- 第九步：怎麼讀這個結果。 「垂直切線」在這個理想模型裡其實是不連續跳躍（比 Kittel 所說的 kink 更劇烈）。原因在第五步那句「無窮長程」：C_p 只以 1/p 衰減，這在物理上是一個極端理想化。真實金屬的有效力常數帶有 Friedel 振盪的 cos(2k_F pa)/(pa)³ 型式，衰減快得多，所以真實的 Kohn anomaly 是一個對數型的微弱扭結，需要高解析度中子散射才看得到。Kittel 這題把衰減放慢到 1/p，就是為了讓奇異性變得乾淨、容易證明。
- 第十步：量綱檢查。 [A] 由 C_p = A sin(·)/(pa) 得 [A] = [力常數][長度] = N。於是 πA/(Ma) 的量綱是 N/(kg·m) = s⁻² ✔ 確實是 ω² 的量綱。
Problem 4 的 ∂ω²/∂K（書上要求「find an expression」，未給答案）
∂ω²/∂K = (2A/M) Σ_{p≥1} sin(p k₀ a) sin(p K a) = (A/M) Σ_{p≥1} [cos(p(k₀−K)a) − cos(p(k₀+K)a)]
p.103A、k₀：題目給定的常數；a：層間距；M：原子質量；K：波向量。
Problem 4 的 ω²(K) 封閉形式（書上完全未提及）
ω²(K) = 0 （0 < K < k₀） ； ω²(K) = π A /(M a) （k₀ < K < π/a）
p.103A、k₀：題目給定常數；M：原子質量；a：層間距。值域限定在第一布里淵區 0 < Ka < π，並假設 0 < k₀a < π。
k₀ 是什麼：Kohn anomaly 為什麼是費米面的指紋（跨章連結）我的判讀
Kittel 對 k₀ 的唯一提示是那句「Such a form is expected in metals.」——沒有說 k₀ 是什麼。但把書裡散落的幾處拼起來，答案很明確。
線索一（本章 p.94）：「In metals the effective forces may be of quite long range and are carried from ion to ion through the conduction electron sea.」離子之間的有效交互作用是傳導電子傳遞的。
線索二（第 14 章 p.407）：Kittel 在討論屏蔽庫侖位能時寫道，Thomas–Fermi 近似只在 q ≪ k_F 成立，「More complete calculations with all wavevectors included exhibit spatial oscillations, called Friedel oscillations, in 2k_F r」。電子氣對點電荷的屏蔽響應會以 2k_F 為週期在實空間振盪。
線索三（第 14 章 p.409，Screening and Phonons in Metals）：Kittel 明確把電子的介電函數接到縱聲學聲子上——「For longitudinal modes the total dielectric function, ions plus electrons, must be zero」，並由此導出 (50) ω = vK、v = (m/3M)^{1/2} v_F，還拿鉀去對數據。所以「電子屏蔽決定聲子頻率」這件事，Kittel 是在第 14 章才正式處理的。
線索四（第 14 章 p.422 與第 18 章 p.532–533）：Peierls 指出一維金屬對波向量 G = 2k_F 的靜態晶格畸變不穩定，因為那會在費米面上開能隙、降低電子能量。
把四條線索接起來：k₀ 扮演的就是 2k_F 的角色。物理機制是這樣的——電子氣對外加擾動的響應（極化率 χ(K)）在 K = 2k_F 處是非解析的，因為那是「把電子從費米面一側散射到另一側」這個過程剛好開始／停止可能的臨界波向量。費米面越銳利（溫度越低、越乾淨），非解析性越尖。聲子的力常數被這個響應重整，於是 ω(K) 在 K = 2k_F 繼承了同樣的非解析性。Kohn anomaly 因此是一張「用聲子量電子」的照片：從色散曲線上扭結的位置，可以直接讀出費米面在該方向的直徑。
數字感（以鈉為例）：書本 p.139 的 Table 1 給鈉的費米波向量 k_F = 0.92 × 10⁸ cm⁻¹ = 0.92 Å⁻¹，所以 2k_F = 1.84 Å⁻¹。鈉沿 [001] 的區邊界在 2π/a ≈ 1.49 Å⁻¹（a ≈ 4.23 Å，此晶格常數為外部資料）。2k_F 比區邊界還大，所以扭結不會直接出現在第一區裡沿 [001] 的位置上，而是要摺回來——出現在 |G − 2k_F| 對應的 K 上。這正好又是前面「選 G 把 K 摺回第一區」那條規則的又一次應用。
後面哪裡會用到。第 6 章（Free Electron Fermi Gas，書本 p.131）建立 k_F 與費米面的概念；第 9 章（Fermi Surfaces and Metals，書本 p.221）教你怎麼量真實金屬的費米面（de Haas–van Alphen 等）；第 14 章 p.409 把屏蔽與聲子正式接上。Kohn anomaly 是這三章與第 4 章之間的接點，而 Kittel 只把它放在一道習題裡——這就是為什麼這題非講不可。
Problem 7（Soft phonon modes）：題目全文，含續頁與兩處校註我的判讀
這題的題幹跨兩頁：p.103 的後半到 p.104 的上半。要完整讀懂必須一起看這兩頁。
原文逐字（p.103）：「*7. Soft phonon modes. Consider a line of ions of equal mass but alternating in charge, with e_p = e(−1)^p as the charge on the pth ion. The interatomic potential is」——句子在此被分頁切斷，接續在 p.104：「the sum of two contributions: (1) a short-range interaction of force constant C_1R = γ that acts between nearest neighbors only, and (2) a coulomb interaction between all ions. (a) Show that the contribution of the coulomb interaction to the atomic force constants is C_pC = 2(−1)^p e²/p³a³, where a is the equilibrium nearest-neighbor distance. (b) From (16a) show that the dispersion relation may be written as ... where ω₀² ≡ 4γ/M and σ = e²/γa³. (c) Show that ω² is negative (unstable mode) at the zone boundary Ka = π if σ > 0.475 or 4/7ζ(3), where ζ is a Riemann zeta function. Show further that the speed of sound at small Ka is imaginary if σ > (2 ln 2)⁻¹ = 0.721. Thus ω² goes to zero and the lattice is unstable for some value of Ka in the interval (0, π) if 0.475 < σ < 0.721. Notice that the phonon spectrum is not that of a diatomic lattice because the interaction of any ion with its neighbors is the same as that of any other ion.」
p.103 頁底的腳註逐字：「*This problem is rather difficult.」
校註一（印刷錯誤，經獨立推導與 400–900 dpi 目視雙重確認）。p.104 那條色散關係在書上印的是：
ω²/ω₀² − sin² ½Ka + σ Σ_{p−1}^{∞} (−1)^p (1 − cos pKa) p⁻³ ,
其中第一個關係符號印成減號、求和下標印成 p−1。這兩處都是排版錯誤：正確的式子是
ω²/ω₀² = sin² ½Ka + σ Σ_{p=1}^{∞} (−1)^p (1 − cos pKa) p⁻³ 。
判斷依據有三：(i) 本講解後面從 (16a) 獨立推出的結果就是等號版本；(ii) 若取減號，整條式子沒有等號、根本不是一個方程；(iii) 求和變數當然從 p = 1 開始（p = 0 沒有意義）。把該頁以 900 dpi 算繪放大目視，確認紙面上印的確實是「−」與「p−1」，因此這是原書的錯字，不是抽取錯誤。
校註二（記號歧義）。「σ > 0.475 or 4/7ζ(3)」中的「4/7ζ(3)」要讀成 4/(7ζ(3)) 而非 (4/7)ζ(3)。驗算：ζ(3) = 1.2020569，4/(7 × 1.2020569) = 0.47538 ✔ 與印出的 0.475 相符；而 (4/7)ζ(3) = 0.687 ✘。
校註三（符號）。「ω₀² ≡ 4γ/M」用的是恆等定義符號 ≡（PDF 中的 /H11013，已正確解碼），而同一句的「σ = e²/γa³」用的是普通等號。這個區別是刻意的：ω₀² 是一個新定義的縮寫，σ 則是由既有量組成的組合量。不過這個區別只在本頁成立，不能推廣成全書體例：Kittel 用 ≡ 並不嚴格，例如 p.94 把小 Ka 展開印成「cos Ka ≡ 1 − ½(Ka)²」、p.98 的 (23)(24) 兩個近似根也印成「ω² ≡ …」，那幾處語意都是約等於。
Problem 7 (a)：庫侖交互作用的力常數為什麼是 2(−1)^p e²/p³a³補推導
書上只說 Show that，一步沒給。這一步其實只有三行，但需要先把「力常數」的定義釘清楚，否則因子 2 會算錯。
- 第一步：釘住力常數的定義。 本章 p.91 定義力常數時說「the force on the plane s caused by the displacement of the plane s + p is proportional to the difference u_{s+p} – u_s」，寫成 F = C_p(u_{s+p} − u_s)。對一個只依賴距離的兩體位能 φ(r)，沿鏈方向（縱模）的展開給 C_p = φ″(r)|_{r = pa}。（這是標準結果：把 φ(pa + Δu) 展到二階，二階係數即 φ″，而恢復力 = −∂U/∂Δu = −φ″Δu。）
- 第二步：寫下第 0 個離子與第 p 個離子之間的庫侖位能。 電荷是 e_p = e(−1)^p，所以 e₀ = +e、e_p = e(−1)^p。CGS 制下 φ(r) = e₀e_p/r = e²(−1)^p / r。
- 第三步：微分兩次。 dφ/dr = −e²(−1)^p/r²；d²φ/dr² = +2e²(−1)^p/r³。
- 第四步：代 r = pa。 C_p^C = 2(−1)^p e²/(p³a³) ✔ 與題目給的答案完全一致。因子 2 來自 1/r 的二階微分（r⁻¹ → 2r⁻³），不是來自左右兩個鄰居——這是最容易搞錯的一點：(16a) 的 Σ_{p>0} 已經把左右對稱性處理掉了。
- 第五步：讀懂符號的交替。 (−1)^p 表示偶數距離的離子同號（互斥，C_p^C > 0，等效於正的彈簧常數）、奇數距離的離子異號（互吸，C_p^C < 0）。最近鄰 p = 1 是異號，庫侖項給的是負的力常數 −2e²/a³——庫侖吸引在縱向上是一個不穩定化的貢獻，因為把兩個相吸的離子拉近，吸引力還會變大。全靠短程項 γ 撐住穩定性。這就是整題的物理核心。
- 第六步：量綱檢查。 e²/a³ 在 CGS 下量綱為 (靜電單位)²/cm³ = erg/cm² = dyne/cm ✔ 正是力常數的量綱。因此 σ = e²/γa³ 是無量綱的，是「庫侖硬度／短程硬度」的純比值。
庫侖交互作用對原子力常數的貢獻（Problem 7a 的答案，原文給出結果）
C_p^C = 2(−1)^p e² / (p³ a³)
p.104C_p^C：相隔 p 個間距的兩個離子之間、源自庫侖交互作用的力常數；e：離子電荷量（第 p 個離子帶 e_p = e(−1)^p）；a：平衡最近鄰距離（原文：「where a is the equilibrium nearest-neighbor distance」）；採 CGS 制。
Problem 7 (b)：把兩種力常數代進 (16a)，導出色散關係補推導
這一步是全題最機械但也最能看出 ω₀² ≡ 4γ/M 為什麼要那樣定義的地方。書上只寫「From (16a) show that」。
總力常數是兩部分相加：C_p = C_p^R + C_p^C，其中 C_1^R = γ、C_p^R = 0（p ≥ 2），而 C_p^C = 2(−1)^p e²/(p³a³)。
- 第一步：搬出 (16a)。 ω² = (2/M) Σ_{p>0} C_p (1 − cos pKa)。
- 第二步：算短程部分。 只有 p = 1 有貢獻：(2/M)·γ·(1 − cos Ka)。用半角恆等式 1 − cos Ka = 2 sin²(½Ka)，得 (2γ/M)·2sin²(½Ka) = (4γ/M) sin²(½Ka)。
- 第三步：這就是 ω₀² 的來源。 定義 ω₀² ≡ 4γ/M，短程部分正好寫成 ω₀² sin²(½Ka)。對照本章 p.92 的 (9)：ω² = (4C/M) sin²(½Ka)——ω₀ 就是「只有最近鄰力常數 γ 的單原子鏈」的區邊界最高頻率。題目把它取作單位，是為了讓最後的答案只剩一個參數 σ。這也解釋了 ≡ 這個符號：它是定義，不是推導結果。
- 第四步：算庫侖部分。 (2/M) Σ_{p≥1} [2(−1)^p e²/(p³a³)] (1 − cos pKa) = (4e²/(Ma³)) Σ_{p≥1} (−1)^p (1 − cos pKa) p⁻³。
- 第五步：兩部分相加並同除以 ω₀² = 4γ/M。 庫侖部分的前因子變成 (4e²/(Ma³)) × (M/4γ) = e²/(γa³) = σ ✔ ——σ 的定義也是為了讓這個因子剛好化簡成一個字母。
- 第六步：得到答案。 ω²/ω₀² = sin²(½Ka) + σ Σ_{p=1}^{∞} (−1)^p (1 − cos pKa) p⁻³ 。這與 p.104 印出的式子逐項相符，唯二的差別是書上把關係符號印成「−」、把求和下標印成「p−1」——見前面的校註一。既然本推導是獨立從 (16a) 出發的，這兩處確定是原書排版錯誤。
- 第七步：讀懂式子的結構。 第一項 sin²(½Ka) 恆為正，是短程力提供的穩定化；第二項的每一項在 K ≠ 0 時 (1 − cos pKa) ≥ 0，乘上交替的 (−1)^p 之後奇數 p 貢獻負、偶數 p 貢獻正，而 1/p³ 讓 p = 1 的負貢獻壓倒性地大。所以庫侖項整體是負的、是不穩定化的，σ 則是它相對於短程穩定化的權重。整條式子就是一場拔河：σ 小則短程勝、晶格穩；σ 大則庫侖勝、ω² 變負。
- 第八步：為什麼「不是雙原子晶格」。 原文最後一句提醒：「Notice that the phonon spectrum is not that of a diatomic lattice because the interaction of any ion with its neighbors is the same as that of any other ion.」意思是——雖然電荷交替（看起來像 NaCl 型的雙原子基底），但質量相同、力常數也對每個離子一模一樣（C_p 只依賴 p，不依賴 s 的奇偶）。所以原胞其實只含一個離子、布里淵區是 −π/a < K ≤ π/a、只有一支聲學支，沒有光學支、沒有頻率間隙。這與 Problem 3 的雙原子鏈形成完美對照：那裡是質量不同造成基底不對稱，這裡電荷雖然交替但動力學是均勻的。
Problem 7(b) 的色散關係（含原書印刷校註）
ω²/ω₀² = sin²(½Ka) + σ Σ_{p=1}^{∞} (−1)^p (1 − cos pKa) p⁻³
p.104ω₀² ≡ 4γ/M（原文用恆等定義符號 ≡）：只含最近鄰短程力 γ 的單原子鏈之區邊界最高頻率平方；σ = e²/γa³：無量綱參數，庫侖硬度與短程硬度之比；M：離子質量（題目設兩種離子等質量）；a：最近鄰距離；K：波向量，第一區為 −π/a < K ≤ π/a。
Problem 7 (c)：兩個臨界 σ 怎麼算出來（0.475 與 0.721）補推導
書上只給了兩個數字與它們的封閉形式，一步推導都沒有。這兩個數各自來自一個經典的級數求和，值得完整寫出來。
- 第一步：區邊界 Ka = π。 sin²(½ · π) = sin²(π/2) = 1。至於和式：cos(pπ) = (−1)^p，所以 1 − cos pKa = 1 − (−1)^p，偶數 p 給 0、奇數 p 給 2。
- 第二步：把和式縮到奇數項。 Σ_{p≥1} (−1)^p (1 − (−1)^p) p⁻³ = Σ_{p 奇} (−1)^p · 2 · p⁻³ = −2 Σ_{p 奇} p⁻³（奇數 p 的 (−1)^p = −1）。
- 第三步：求奇數項的 ζ。 Σ_{p 奇} p⁻³ = Σ_{所有 p} p⁻³ − Σ_{p 偶} p⁻³ = ζ(3) − Σ_{m≥1}(2m)⁻³ = ζ(3) − (1/8)ζ(3) = (7/8)ζ(3)。
- 第四步：合起來。 和式 = −2 × (7/8)ζ(3) = −(7/4)ζ(3)。所以 ω²/ω₀²|_{Ka=π} = 1 − (7/4)ζ(3) σ。
- 第五步：求臨界值。 ω² < 0 ⟺ 1 − (7/4)ζ(3)σ < 0 ⟺ σ > 4/(7ζ(3))。數值：ζ(3) = 1.2020569，4/(7 × 1.2020569) = 0.47538 ✔ 與原文的 0.475 相符，也確認了原文「4/7ζ(3)」該讀成 4/(7ζ(3))。
- 第六步：小 Ka 展開。 sin²(½Ka) ≈ (Ka)²/4；1 − cos pKa ≈ (pKa)²/2。代入和式：Σ_{p≥1} (−1)^p · (pKa)²/2 · p⁻³ = ((Ka)²/2) Σ_{p≥1} (−1)^p p⁻¹。
- 第七步：那個交錯調和級數。 Σ_{p≥1} (−1)^p/p = −(1 − ½ + ⅓ − ¼ + …) = −ln 2。所以和式 ≈ −(Ka)²(ln 2)/2。
- 第八步：合起來看聲速。 ω²/ω₀² ≈ (Ka)²[¼ − σ(ln 2)/2]。既然 ω² = v²K²，v² ∝ ¼ − σ(ln 2)/2。v 為虛數 ⟺ ¼ < σ(ln 2)/2 ⟺ σ > 1/(2 ln 2)。數值：1/(2 × 0.693147) = 0.72135 ✔ 與原文的 0.721 相符。
- 第九步：把兩個臨界值擺在一起讀。 σ < 0.475：兩端都穩，晶格穩定。0.475 < σ < 0.721：長波端穩（v 實數）、區邊界端不穩（ω² < 0），由連續性，ω² 必在 (0, π) 內某處過零。σ > 0.721：連長波都不穩，整條聲學支塌掉。這就完整證明了原文的結論：「Thus ω² goes to zero and the lattice is unstable for some value of Ka in the interval (0, π) if 0.475 < σ < 0.721.」
- 第十步：數值核對整條曲線。 把 ω²/ω₀²(Ka) 用前 4000 項數值求和：σ = 0.475 時 Ka = π 處為 +0.00079 ≈ 0（臨界點 ✔）；σ = 0.55 時零點落在 Ka ≈ 2.32；σ = 0.60 落在 1.88；σ = 0.65 落在 1.41；σ = 0.70 落在 0.75；σ = 0.72 落在 0.19；σ = 0.721 時已幾乎處處 ≤ 0。σ 從 0.475 增大到 0.721 的過程中，軟化的波向量從區邊界 π 連續地掃到 0。
Problem 7(c) 的兩個臨界值（原文給出數值，推導與級數為本講解補出）
區邊界失穩：σ_c1 = 4/(7ζ(3)) = 0.47538 ； 聲速虛數：σ_c2 = 1/(2 ln 2) = 0.72135
p.104σ = e²/γa³：無量綱的庫侖／短程硬度比；ζ(3) = 1.2020569…（Apéry 常數）；ln 2 = 0.693147…。
Problem 7 的物理：ω² 為負到底是什麼意思，以及它接到哪一章我的判讀
這題之所以值得花力氣，是因為它把本章 p.100 一句看起來像註腳的話，變成了一個可以算的具體現象。
那句話（p.99–100 交界，逐字）：「What is the sign of ω? The equations of motion such as (2) are equations for ω², and if this is positive then ω can have either sign, + or −. But the energy of a phonon must be positive, so it is conventional and suitable to view ω as positive. If the crystal structure is unstable, then ω² will be negative and ω will be imaginary.」（此段跨 p.99 與 p.100 分頁，「But the」與「energy of a phonon」之間即為頁界；最後一句是本題的直接依據。）
「ω 為虛數」在數學上是什麼。試解是 u_s ∝ exp(−iωt)。若 ω = i|ω|，則 exp(−iωt) = exp(|ω|t)——位移不再振盪，而是指數成長。簡諧近似下這個成長沒有上限；真實晶體裡它會一直長到非簡諧項（第 5 章 Anharmonic Crystal Interactions，書本 p.119）把它剎住，晶格於是停在一個新的、位移過的平衡結構上。這就是結構相變。
「軟模」為什麼叫軟模。在相變溫度附近，某個模式的 ω² 隨溫度線性趨近於零（而不是突然跳到負值）。ω² → 0 表示該模式的恢復力消失，晶格對這個特定的位移圖樣「軟」了。跨過臨界點後 ω² 變負，晶格就沿著那個圖樣塌陷下去。Problem 7 用 σ 取代溫度當作控制參數，把這整個機制濃縮成一個可以手算的模型。
軟化波向量的三種情形，對應三種相變：
| 軟化發生在 | 新結構 | 例子／章節 |
|---|---|---|
| K = 0（Γ 點） | 整個原胞內的原子相對位移，產生自發極化 | 鐵電相變，第 16 章 Soft Optical Phonons，書本 p.473 |
| K = 區邊界（π/a） | 原胞加倍（cell doubling），週期變兩倍 | 反鐵電／二聚化；一維情形即 Peierls 不穩定（第 14 章 p.422、第 18 章 p.532–533） |
| K = 區內某個 K*（本題 0.475 < σ < 0.721 的情形） | 週期 2π/K* 的調製結構，一般與晶格不可公度 | 電荷密度波、不可公度相 |
第三種正是 Problem 7 最漂亮的地方：那個窗口不是「穩定」與「不穩定」之間的模糊地帶，而是一個有自己物理內容的相。上一個區塊算出的 K*a 表顯示，σ 從 0.475 走到 0.721 時 K* 從 π 連續掃到 0——換句話說，調製週期可以是任何值，不受晶格週期約束。
Kittel 在書中哪裡正式處理這件事。第 16 章（Dielectrics and Ferroelectrics，書本 p.453 起）的 Structural Phase Transitions（p.467）與 Soft Optical Phonons（p.473）。那裡 p.469 有一句很關鍵的話（逐字）：「One may define the character of the transition in terms of the dynamics of the lowest frequency (“soft”) optical phonon modes. If a soft mode can propagate in the crystal at the transition, then the transition is displacive. If the soft mode is only diffusive (non-propagating) there is really not a phonon at all, but is only a large amplitude hopping motion between the wells of the order-disorder system.」（最後這句跨 p.469／p.470 分頁，「not a phonon at all,」與「but is only」之間即為頁界。）位移型（displacive）與序—無序型（order-disorder）的分野，就建立在這個 ω² → 0 的判準上。第 16 章的習題 9 標題也叫 Soft modes and lattice transformations（書本 p.485）。
所以這道「rather difficult」的習題，實質上是第 16 章的先修課。Kittel 把它放在第 4 章，是為了讓讀者在只有 (16a) 這一個工具的時候，就親手算出一次晶格失穩。
其餘四題速記（Problems 1、2、5、6）我的判讀
這四題沒有補回正文缺口，但各有一個值得記住的重點。題目文字為原文逐字，解題方向為本講解補充。
Problem 1. Monatomic linear lattice（p.102–103）。原文：「Consider a longitudinal wave u_s = u cos(ωt − sKa) which propagates in a monatomic linear lattice of atoms of mass M, spacing a, and nearest-neighbor interaction C. (a) Show that the total energy of the wave is E = ½M Σ_s (du_s/dt)² + ½C Σ_s (u_s − u_{s+1})², where s runs over all atoms. (b) By substitution of u_s in this expression, show that the time-average total energy per atom is ¼Mω²u² + ½C(1 − cos Ka)u² = ½Mω²u², where in the last step we have used the dispersion relation (9) for this problem.」
重點：(a) 只是把動能與位能寫出來；(b) 的關鍵是 ⟨cos²⟩ = ⟨sin²⟩ = ½ 的時間平均，以及「每個原子」意味著要除以原子數。最後一步把 ½C(1 − cos Ka) 換成 ¼Mω²，書上說用的是 (9)（原文：「where in the last step we have used the dispersion relation (9) for this problem.」）；由於 (9) ω² = (4C/M)sin²(½Ka) 與 p.92 的 (7) ω² = (2C/M)(1 − cos Ka) 只差一個半角恆等式，兩者完全等價，實際代數時用 (7) 最直接。得到的結論 E = ½Mω²u² 意味著動能時均 ¼Mω²u² 恰等於位能時均——簡諧振子的能量均分定理在晶格上的驗證。這個結果直接餵給 p.99 的量子化 (28)(29)，也是第 5 章熱容計算的起點。
Problem 2. Continuum wave equation（p.103）。原文：「Show that for long wavelengths the equation of motion (2) reduces to the continuum elastic wave equation ∂²u/∂t² = v² ∂²u/∂x², where v is the velocity of sound.」
重點：把 (2) 的 u_{s+1} + u_{s−1} − 2u_s 視為二階中心差分，泰勒展開到 a² 階得 a²∂²u/∂x²，於是 M∂²u/∂t² = Ca²∂²u/∂x²，即 v² = Ca²/M。這與 p.94 的 (15) ω² = (C/M)K²a² 是同一件事的兩種寫法（一個在實空間、一個在倒空間）。這題的意義是把離散晶格與第 3 章的連續彈性理論正式接起來。
Problem 5. Diatomic chain（p.103）。原文：「Consider the normal modes of a linear chain in which the force constants between nearest-neighbor atoms are alternately C and 10C. Let the masses be equal, and let the nearest-neighbor separation be a/2. Find ω(K) at K = 0 and K = π/a. Sketch in the dispersion relation by eye. This problem simulates a crystal of diatomic molecules such as H₂.」
重點：這是 Problem 3 的鏡像。那裡是質量不同、力常數相同，這裡是質量相同、力常數不同。兩者都打開頻率間隙，證明間隙的來源是原胞內的不等價性本身，不管這個不等價來自哪個物理量。最後一句「simulates a crystal of diatomic molecules such as H₂」點出了物理圖像：10C 是分子內的強鍵，C 是分子間的弱鍵；光學支就是分子的內部伸縮振動，聲學支是分子整體的平移振動。注意最近鄰距離是 a/2，所以重複距離仍是 a、區邊界仍在 π/a。
Problem 6. Atomic vibrations in a metal（p.103）。原文：「Consider point ions of mass M and charge e immersed in a uniform sea of conduction electrons. ... If one ion is displaced a small distance r from its equilibrium position, the restoring force is largely due to the electric charge within the sphere of radius r centered at the equilibrium position. Take the number density of ions (or of conduction electrons) as 3/4πR³, which defines R. (a) Show that the frequency of a single ion set into oscillation is ω = (e²/MR³)^{1/2}. (b) Estimate the value of this frequency for sodium, roughly. (c) From (a), (b), and some common sense, estimate the order of magnitude of the velocity of sound in the metal.」
重點：(a) 用高斯定律，半徑 r 球內的電子電荷是 −(4/3)πr³ · n·e = −e(r/R)³，恢復力 = −e²r/R³，故 ω² = e²/MR³。注意這不是離子的等離子頻率，而是「單一離子在剛性電子背景中振盪」的愛因斯坦（單粒子）頻率：把題目給的數密度 n = 3/(4πR³)（原文：「Take the number density of ions (or of conduction electrons) as 3/4πR³, which defines R.」）代進離子等離子頻率的定義，得 ω_pi² = 4πne²/M = 3e²/(MR³)，是本題答案的 3 倍，也就是 ω = ω_pi/√3。（第 14 章 p.410 的 (52) 則是另一件事：ω² = 4πne²/μ、1/μ = 1/M + 1/m，原文明說那是「the electron plasma frequency (20), but with the reduced mass correction for the motion of the positive ions」，其中 n 是電子密度，與本題的 e²/MR³ 不是同一個量。）(c) 的「common sense」指的是：ω 是區邊界的頻率量級、區邊界波向量約 π/a ~ 1/R，所以 v ~ ωR = (e²/MR)^{1/2}。這題與第 14 章 p.409 的 (50) v = (m/3M)^{1/2}v_F 是同一個物理的兩種近似，而且它示範了本章 Fig. 11 那些 10¹² Hz 的頻率為什麼是這個量級。
常見誤解與自我檢核我的判讀
讀完這一段，用下面幾條檢查自己有沒有真的懂。每一條都是實際上很常見的誤解。
（1）「(33) 裡的 G 是實驗決定的。」 錯。G 是記帳方式。物理只認摺回第一區的 K 與對應的 ω。同一個聲子可以掛在任何一個布拉格點旁邊去量；選哪個 G 只影響訊號強度，不影響 (K, ω)。
（2）「(31)、(33)、SUMMARY 的選擇律是三條不同的式子。」 錯。移項並利用「G 與 −G 都是倒晶格向量」之後，它們是同一條。全章其實只有兩條：產生（k = k′ + K + G）與吸收（k = k′ − K + G）。
（3）「聲子帶動量 ℏK。」 錯，而且 Kittel 在 p.100 專門開一節（PHONON MOMENTUM）擋這句話。聲子表現得像帶 ℏK，但不帶真實線性動量，因為聲子座標是相對座標（原文用 H₂ 分子的 r₁ − r₂ 做類比）。只有 K = 0 的均勻模式對應質心座標，才帶真動量。
（4）「中子散射時，聲子創生對應中子得到能量。」 錯，恰好相反。(34) 上號（創生）⇒ ℏ²k²/2M_n − ℏ²k′²/2M_n = +ℏω ⇒ k′ < k ⇒ 中子失去能量。而且創生過程在 T = 0 也存在（因為是從真空創生），吸收過程則需要既有聲子，強度正比於 n(ω)，T → 0 時消失。
（5）「X 光波長對、所以也可以量聲子。」 半對。波長確實對（判準一過關），但 7 keV 的光子要分辨 15 meV 的變化，需要 ΔE/E ~ 10⁻⁶ 的解析度，一般 X 光源做不到（判準二不過關）。中子的優勢在於 M_n 讓「Å 波長」與「10 meV 能量」同時成立。
（6）「Fig. 11 的三個面板橫軸長度一樣，可以直接比較斜率。」 錯。[001] 走到 |q| = 1、[110] 走到 0.707、[111] 走到 0.866（皆以 2π/a 為單位），因為體心立方第一區是菱形十二面體，H、N、P 三點到原點的距離不同。比較聲速必須換算成 ω/K。
（7）「Fig. 11 的 [110] 面板有兩條橫支，是數據品質問題。」 錯。這是對稱性的必然：第 3 章 p.83 Fig. 20 圖說明講「The two transverse modes are degenerate for propagation in the [100] and [111] directions.」——唯獨 [110] 不簡併。兩支的有效彈性常數分別是 C₄₄ 與 ½(C₁₁ − C₁₂)，而鈉的後者只有前者的八分之一（Table 12，p.84），所以下面那支特別軟。
（8）「Problem 7 的電荷交替，所以它是雙原子晶格。」 錯，而且原文特別提醒了：「the interaction of any ion with its neighbors is the same as that of any other ion」。質量相同、力常數對每個離子相同 ⇒ 原胞只含一個離子 ⇒ 只有一支聲學支、沒有光學支、沒有間隙。電荷交替只影響 C_p 的符號，不影響動力學的均勻性。
（9）「Kohn anomaly 是 Kittel 正文講過的東西。」 錯。全書文字層搜尋顯示「Kohn anomaly」這個字串只出現在第 4 章章前目錄（p.89）、Problem 4（p.103）、索引（p.677）三處，正文從未展開；唯一的相關提及是第 9 章 p.242 把「the Kohn effect」列為量費米面動量分布的方法之一，一句話帶過。同樣地「thermal neutron」全書零次出現。相關的物理要到第 14 章 p.407（Friedel 振盪）與 p.409（Screening and Phonons in Metals）才正式登場。
（10）與熱物性研究的接點（寫給做熱量測的人）。 這一段量到的 ω_j(K) 不只是「一張漂亮的圖」——它是所有聲子熱物性的唯一輸入。第 5 章的熱容 C_V 由態密度 D(ω) 決定，而 D(ω) 完全由 ω_j(K) 產生（第 5 章 p.108 起）；熱導率 κ = ⅓Cvℓ 裡的群速度 v 就是 (13) 的 grad_K ω。所以任何從第一原理算熱容或熱導的流程，第一步都是算出這張色散圖，第二步才是積分。Fig. 11 的形狀直接決定了鈉的低溫熱容偏離 Debye T³ 的程度——[110] 那支特別軟的橫支會在低頻堆出額外的態密度，讓實際的 Debye 溫度隨溫度漂移。這是「量色散」與「量熱」之間的直接橋樑。