第 4 章的地圖：章首目錄頁逐條轉錄（書 p.89）
章首頁印的標題是 4 Phonons I. Crystal Vibrations（聲子 I：晶體振動）。照原文印的順序與頁碼，這一章的小節是：VIBRATIONS OF CRYSTALS WITH MONATOMIC BASIS（p.91，下轄 First Brillouin zone p.93、Group velocity p.94、Long wavelength limit p.94、Derivation of force constants from experiment p.94）；TWO ATOMS PER PRIMITIVE BASIS（p.95）；QUANTIZATION OF ELASTIC WAVES（p.99）；PHONON MOMENTUM（p.100）；INELASTIC SCATTERING BY PHONONS（p.100）；SUMMARY（p.102）；PROBLEMS（p.102）。
習題共七題，橫跨 p.102–103：1. Monatomic linear lattice（102）、2. Continuum wave equation（103）、3. Basis of two unlike atoms（103）、4. Kohn anomaly（103）、5. Diatomic chain（103）、6. Atomic vibrations in a metal（103）、7. Soft phonon modes（103）。章首頁最後單獨一行寫著：“Chapter 5 treats the thermal properties of phonons.”——第 5 章處理聲子的熱性質。
從這張目錄讀出的結構訊號：一章其實只做兩個模型，而量子化排在最後我的判讀
這張目錄裡藏了一個很重要的結構訊號，值得在讀正文之前先看出來：「量子化」排在第三順位，在兩個純古典的彈簧－質點模型之後。QUANTIZATION OF ELASTIC WAVES 要到 p.99 才出現；也就是說 Kittel 整整八頁（p.91–98）都在解古典的力學問題，而「phonon（聲子）」這個詞要到第九頁才被正式定義成量子。
這不是編排上的疏忽，而是這門學科真實的邏輯順序：色散關係 ω(K) 本身是純古典的結果。把彈性波量子化，改變的只是「每一個模式（K, 偏振）裡能裝多少能量」——能量變成 (n + ½)ℏω 的階梯——而完全不改變 ω(K) 這條曲線長什麼樣。這一點到第 5 章算熱容時會兌現成一句非常關鍵的話：模式由古典晶格動力學決定，佔據數由 Planck 分布決定，兩件事可以分開算。本段（p.89–92）做的正是前者當中最基本的那一步。
本段要解決的問題，以及為什麼從「一個原子的原胞」開始我的判讀
第 3 章的結尾把晶體當成連續彈性介質，得到的色散關係是一條直線：ω = v K，聲速 v 由彈性勁度與密度決定（例如 p.82 的式 (60)，[100] 方向縱波 v_s = (C₁₁/ρ)^{1/2}）。那套理論漂亮，但它從頭到尾不知道晶體是由離散的原子組成的——波長要多短都可以，頻率就一路線性上去，沒有上限。
這顯然是錯的。真實晶體裡波長不可能短於原子間距的兩倍（比原子還密的波形，原子看不到也表現不出來），頻率因此一定有上限。本段的任務就是把「離散性」放回去，看看色散關係會被改成什麼樣子。 Kittel 選的模型是最小可行模型：一串等間距的原子面，只有最近鄰之間有彈簧。
為什麼是「一個原子的原胞（monatomic basis）」？因為這是唯一能讓整個問題只剩下一個未知振幅的情形。原胞裡放兩個原子，方程式立刻變成 2×2 的矩陣問題，色散關係分裂成聲學支與光學支——那是 p.95 起 TWO ATOMS PER PRIMITIVE BASIS 一節的事。先做單原子鏈，是為了把「離散性造成什麼後果」這件事單獨隔離出來看清楚，不要被「基底內部自由度」的效應混進來。
本段涵蓋書 p.89–92，終點是式 (9)：ω = (4C/M)^{1/2}|sin ½Ka|。這條式子是整個聲子理論的地基。之後的第一布里淵區、群速、長波極限、Debye 模型、聲子熱導——全部站在它上面。
Fig. 1：固體中的基本激發譜系（書 p.90）——聲子在整本書裡的位置
第 4 章正文開始之前，Kittel 先在 p.90 放了一張只有兩欄、七列的表，圖說是：“Important elementary excitations in solids.”（固體中重要的基本激發。）兩欄的表頭分別是 Name 與 Field。內容如下：
| Name（基本激發的名字） | Field（它所對應的場） |
|---|---|
| Electron（電子） | —— （本欄在原書中為空白） |
| Photon（光子） | Electromagnetic wave（電磁波） |
| Phonon（聲子） | Elastic wave（彈性波） |
| Plasmon（電漿子） | Collective electron wave（集體電子波） |
| Magnon（磁子） | Magnetization wave（磁化波） |
| Polaron（極化子） | Electron + elastic deformation（電子＋彈性形變） |
| Exciton（激子） | Polarization wave（極化波） |
校註（配對依據）：這份 PDF 的文字層把兩欄的閱讀順序打亂了——直接抽出來會得到「Electron / Photon / Electromagnetic wave / Elastic wave / Polarization wave / Electron + elastic deformation / Magnetization wave / Collective electron wave / Phonon / Plasmon / Magnon / Polaron / Exciton」這種左右欄交錯又跳行的序列，照字面讀會把 Phonon 配到 Polarization wave 上去。上表的配對是用 PDF 文字區塊的 y 座標還原的：Name 欄的文字都落在 x ≈ 261–290，Field 欄落在 x ≈ 315–416，同一列的兩個字串 y 座標幾乎完全相同——Photon 162.87 ↔ Electromagnetic wave 162.87；Phonon 180.53 ↔ Elastic wave 180.87；Plasmon 197.53 ↔ Collective electron wave 197.53；Magnon 215.37 ↔ Magnetization wave 215.53；Polaron 232.70 ↔ Electron + elastic deformation 232.53；Exciton 250.04 ↔ Polarization wave 249.87。而 Electron 位於 y = 145.53，Field 欄在這個高度上沒有任何文字，所以它的對應欄是真的空白，不是抽取遺漏。
讀懂 Fig. 1：這張表其實在說「準粒子 ↔ 場」的對偶我的判讀
這張表的價值不在於背七個名詞，而在於它把整本書的後半段一次攤開，並且宣告了一個貫穿全書的思考模式：在固體裡，每一種可以攜帶能量與動量的集體運動，都可以用兩種語言描述——一種是「波」（場），一種是「粒子」（激發）。兩欄講的是同一件事。
順著這個角度重讀每一列：
- Phonon ↔ Elastic wave。聲子就是彈性波的量子。彈性波是第 3 章的主角（連續介質），聲子是第 4 章 p.99 才正式定義的量子。所以本章的整條主線可以一句話概括：先把「彈性波」修正成晶格上的離散波（p.91–98，純古典），再把它量子化成「聲子」（p.99 起）。表上那一列的左欄與右欄，就是這一章的前後兩半。
- Photon ↔ Electromagnetic wave。這是唯一一列讀者在進固態物理之前就已經熟悉的對偶。Kittel 把它放在聲子正上方，用意很明顯：聲子之於彈性波，就像光子之於電磁波。這個類比在計算上是嚴格的——兩者都是諧振子場的量子，都服從 Planck 分布，都用 (n + ½)ℏω。差別只有兩個，而且都很重要：(a) 電磁波的色散關係 ω = cK 在任意 K 都成立，聲子的 ω(K) 在布里淵區邊界會彎下來；(b) 電磁場有兩個橫向偏振，彈性波有一個縱向加兩個橫向共三個偏振。這兩個差別正是 Debye 模型與黑體輻射公式長得像卻不完全一樣的來源。
- Electron 的 Field 欄為什麼空白。這是全表最耐人尋味的一格。其他六列都是「多體系統的集體運動」——磁化、極化、電子密度的集體振盪、晶格形變——這些量在古典物理裡本來就是一個連續的場，量子化之後才變成粒子。電子不是：電子是基本粒子，它不是某個古典場的量子化產物。在單電子近似的層次上，電子波函數的「場」是 Schrödinger 場，但那不是一個古典可觀測的巨觀場，所以 Kittel 留白。（第二量子化的語言裡當然可以把電子場也寫成場算符，但那超出這本書的層次。）
- Polaron ↔ Electron + elastic deformation。這一列是唯一一個「複合」的激發：一個電子拖著它自己造成的晶格形變一起跑。它同時用到左欄的第一列（電子）與第三列（彈性波），是這張表裡唯一明示「兩種激發可以綁在一起變成第三種」的例子。這也暗示了電子－聲子耦合的存在——而電子－聲子耦合正是金屬電阻率、超導配對機制、以及半導體載子遷移率的溫度相依性的共同根源。
為什麼這張表值得花力氣讀：它是 Kittel 這本書的目錄的另一種寫法。Plasmon 在第 14 章（PLASMONS, POLARITONS, AND POLARONS：電漿子、極化激元、極化子），Magnon 在第 12 章（磁振子），Exciton 在第 15 章（光學過程與激子），Polaron 也在第 14 章。聲子被排在第 4 章，是這一族裡最早登場的——理由是它最容易算（彈簧和質點而已），而且它的數學結構（色散關係、布里淵區、群速、量子化、動量守恆選擇律）會原封不動地套用到後面每一個成員身上。把第 4 章學紮實，後面那幾章有一半的力氣可以省下來。
Fig. 2 與 Fig. 3：把「整個原子面一起動」畫出來（書 p.90）
網頁上不放圖，所以這兩張圖用文字還原。兩張圖都在書 p.90，是為了替 p.91 的「三維退化成一維」這個關鍵動作提供視覺依據。
Fig. 2（縱波），圖說逐字為：“(Dashed lines) Planes of atoms when in equilibrium. (Solid lines) Planes of atoms when displaced as for a longitudinal wave. The coordinate u measures the displacement of the planes.”
圖面上是六條大致垂直的實線，由左到右標記為 s − 1、s、s + 1、s + 2、s + 3、s + 4，每條線上畫了三個黑點代表原子。虛線是平衡位置，彼此間距相等，這個間距在圖中間用一個雙箭頭標成 a。實線是波通過時被推移後的位置。六條實線各自對應一條平衡虛線，但 s 面在該瞬間的位移恰好為零（u_s = 0），所以圖上沒有畫出它的虛線、也沒有畫箭頭，只留下一個孤立的標籤 u_s；實際畫出來的虛線因此只有五條（s−1、s+1、s+2、s+3、s+4）。其餘五面的位移用水平箭頭標出並命名為 u_{s−1}、u_{s+1}、u_{s+2}、u_{s+3}、u_{s+4}——u_{s−1} 的箭頭向左，u_{s+1} 到 u_{s+4} 的箭頭向右且逐面變長，正是波形在空間中的取樣。圖的右側畫一個粗體的水平箭頭 K，指向右方，表示波向量的方向。
這張圖要傳達的核心是：位移 u 的方向與 K 的方向平行（都是水平的），這就是「縱向（longitudinal）」的定義；而且每一條線是整面剛性地平移，面內的原子之間沒有相對運動。
Fig. 3（橫波），圖說逐字為：“Planes of atoms as displaced during passage of a transverse wave.”
圖面上同樣是一排垂直線代表原子面，K 仍然是水平向右的粗箭頭。不同的是：每條線整體被垂直地抬起一小段，抬升量用向上的箭頭標成 u_{s−2}、u_{s−1}、u_s、u_{s+1}、u_{s+2}，其中中間那一面 u_s 的箭頭最長。平衡位置用一條水平參考線上的小虛圓圈標出。
這張圖的核心是：位移 u 的方向與 K 垂直，這就是「橫向（transverse）」的定義。注意橫波有兩個獨立的偏振（可以往上下抬，也可以往紙面內外抬），而縱波只有一個——這正是 p.91 原文那句 “For each wavevector there are three modes as solutions for u_s, one of longitudinal polarization (Fig. 2) and two of transverse polarization (Fig. 3).” 的圖像依據：每個波向量對應三個模式，一縱兩橫。這個「3」會在第 5 章的態密度與 Debye 模型裡變成一個處處出現的因子。
問題的陳述：要求的是 ω 對 K 的函數（書 p.91 第一段）
正文的第一句把題目講得很精確：“Consider the elastic vibrations of a crystal with one atom in the primitive cell. We want to find the frequency of an elastic wave in terms of the wavevector that describes the wave and in terms of the elastic constants.”——考慮原胞中只有一個原子的晶體的彈性振動；我們要用描述該波的波向量、以及彈性常數，把彈性波的頻率表示出來。
讀這一句話：它同時寫死了「要什麼」與「不要什麼」我的判讀
這一句話其實已經把後面九條式子的目標寫死了：要得到的東西是一個函數 ω(K)，它的參數只能是彈性常數（後面具體化成力常數 C）與原子質量 M。不要的東西同樣清楚：不求振幅（線性系統裡振幅是自由的）、不求相位（可以任意選時間原點）、不求任何依賴激發強度的量（諧近似下所有東西都與振幅無關）。
「ω 是 K 的函數」這件事本身就是一個物理結論，不是廢話。 在連續介質裡 ω = vK，v 與 K 無關，所有波長的波都以同一速度傳播，波包不會散開——這叫「無色散」。在離散晶格裡 ω(K) 會變成非線性的，不同波長的波跑不一樣快，波包會散開——這才叫「色散關係（dispersion relation）」。Kittel 是在引出式 (7) 的那一句話裡稱它為 “the dispersion relation ω(K)”，也就是說這個詞印在式 (7) 之前，不是之後。要注意這不是全書第一次用到它：第 3 章 p.82 就已經出現過一次——“On substitution in (57b) this gives the dispersion relation (62)”，指的是 ω²ρ = C₄₄K²，而那條關係其實是無色散的直線。到了 (7)，這個詞才第一次名副其實。
建模假設一（原文）：沿三個高對稱方向，整面同相，三維塌成一維
Kittel 接著寫：“The mathematical solution is simplest in the [100], [110], and [111] propagation directions in cubic crystals. These are the directions of the cube edge, face diagonal, and body diagonal.”——在立方晶體中，數學解在 [100]、[110]、[111] 三個傳播方向最簡單；它們分別是立方體的稜、面對角線、體對角線的方向。
關鍵的那一句是下一句：“When a wave propagates along one of these directions, entire planes of atoms move in phase with displacements either parallel or perpendicular to the direction of the wavevector.”——當波沿著這些方向之一傳播時，整個原子面同相地運動，位移不是平行就是垂直於波向量的方向。
於是：“We can describe with a single coordinate u_s the displacement of the plane s from its equilibrium position. The problem is now one dimensional.”——我們可以用單一個座標 u_s 描述第 s 個面相對於其平衡位置的位移。問題現在是一維的了。
最後補上模式計數：“For each wavevector there are three modes as solutions for u_s, one of longitudinal polarization (Fig. 2) and two of transverse polarization (Fig. 3).”
假設一的深讀：為什麼「整面同相」不是近似，以及 a 為什麼會隨方向改變我的判讀
這個假設是整段推導中被讀者低估得最嚴重的一個。Kittel 用兩行帶過，但它做的事情非常激烈：把 3N 個自由度的問題砍成 N 個。值得把它拆成三個層次看。
第一層：為什麼在這三個方向上「整面同相」是嚴格的，不是近似。 波向量 K 決定了相位如何隨位置變化：相位因子是 exp(iK·r)。如果一個原子面上的所有原子 r 都給出相同的 K·r，那它們的相位就完全相同，整面必然同相。而「K·r 在整個面上為常數」的意思正是：這個面垂直於 K。所以只要選擇的原子面是垂直於 K 的晶面，整面同相就是幾何上的恆等式，不需要任何近似。
真正需要對稱性幫忙的不是「同相」，而是位移方向。一般方向的 K，三個本徵模式的位移方向既不平行也不垂直於 K，而是斜的（稱為 quasi-longitudinal 與 quasi-transverse），這時三個偏振會互相耦合，一維化就失效。[100]、[110]、[111] 是立方晶體的高對稱軸（分別具有四重、二重、三重旋轉對稱），對稱性強制本徵模式必須是純縱或純橫。這件事第 3 章其實已經預告過：p.82 講 [100] 方向時原文說 “Thus for K parallel to [100] the two independent shear waves have equal velocities. This is not true for K in a general direction in the crystal.”
順帶一個容易被忽略的細節：這三個方向的對稱性強度不一樣。[100]（四重軸）與 [111]（三重軸）的兩支橫波簡併（頻率相同）；[110] 只有二重軸，兩支橫波分裂成兩條不同的曲線。第 3 章 p.82 的 “Waves in the [110] Direction” 一節之所以特別有價值（原文說 “because the three elastic constants can be found simply from the three propagation velocities in this direction”），正是因為 [110] 這條線上有三個互不相同的速度，剛好可以解出 C₁₁、C₁₂、C₄₄ 三個常數。
第二層：a 是「面間距」，不是「晶格常數」。 這一點 Kittel 自己在式 (4) 下方就寫明了：“where a is the spacing between planes and K is the wavevector.”，並補一句 “The value to use for a will depend on the direction of K.”（a 的取值將取決於 K 的方向。）這兩句極容易滑過去。書上沒說的是：三個方向的面間距各是多少、fcc 與 bcc 又要怎麼算——那才是實際算數字時最常出錯的地方。a 的定義是垂直於 K 的相鄰原子面之間的距離，對同一個晶體，三個方向的 a 完全不同；而且對 fcc 與 bcc，有些方向不等於把立方晶格常數 a₀ 除以 √(h²+k²+l²)——要看該方向最短的允許倒晶格向量是哪一個（見下表）。面心／體心位置的原子會在「名目上的」兩個面之間再插入一層面，被插入的那幾個方向就對不上；沒有被插入的方向（例如 fcc 的 [111]、bcc 的 [110]）則仍然等於 a₀/√(h²+k²+l²)。
以立方晶格常數 a₀ 為單位，三種立方布拉菲晶格沿三個方向的面間距是：
| 晶格 | K ∥ [100] | K ∥ [110] | K ∥ [111] |
|---|---|---|---|
| 簡單立方 sc | a₀ | a₀/√2 ≈ 0.707 a₀ | a₀/√3 ≈ 0.577 a₀ |
| 面心立方 fcc | a₀/2 | a₀/(2√2) ≈ 0.354 a₀ | a₀/√3 ≈ 0.577 a₀ |
| 體心立方 bcc | a₀/2 | a₀/√2 ≈ 0.707 a₀ | a₀/(2√3) ≈ 0.289 a₀ |
（算法：面間距 d = 2π/|G_min|，G_min 是沿該方向最短的倒晶格向量；fcc 的倒晶格是 bcc，允許的指標為全奇或全偶，所以 [110] 方向最短的是 (220) 而非 (110)；bcc 的倒晶格是 fcc，允許 h+k+l 為偶，所以 [111] 方向最短的是 (222) 而非 (111)。這是第 2 章結構因子的直接後果。）
這個 a 會決定布里淵區邊界在哪裡（K = ±π/a），也會決定最高頻率對應的波長（λ = 2a）。 用錯 a，整條色散曲線的橫軸就整個縮放錯了。而且——這正是任務裡提示的那個伏筆——到了 p.95 的雙原子鏈，「a」的意義還會再被改寫一次：那時候 a 變成原胞的長度（含兩個原子面），而不是相鄰面的距離，布里淵區邊界因此落在半個位置上。讀到那裡時如果沒有在這裡先把 a 的定義釘死，很容易得到一條差兩倍的曲線。
第三層：這個假設什麼時候失效。 (a) K 不在高對稱方向：本徵模式混合，必須回到完整的 3×3 動力學矩陣求解（本書不做，但這是晶格動力學標準的下一步）。(b) 晶體不是立方的：連 [100] 這樣的方向都不保證是對稱軸。(c) 原胞裡不只一個原子：即使 K 在高對稱方向，整面同相依然成立，但「一個面一個座標」不夠用了——原胞裡有幾個原子就要幾個座標，這正是 p.95 起要處理的事。
建模假設二（原文）：彈性能是相對位移的二次式——諧近似
第二段假設在 p.91，原文逐字是：“We assume that the elastic response of the crystal is a linear function of the forces. That is equivalent to the assumption that the elastic energy is a quadratic function of the relative displacement of any two points in the crystal.”——我們假設晶體的彈性響應是力的線性函數；這等價於假設彈性能是晶體中任意兩點相對位移的二次函數。
接著是兩句交代「為什麼可以只留二次項」：“Terms in the energy that are linear in the displacements will vanish in equilibrium—see the minimum in Fig. 3.6. Cubic and higher-order terms may be neglected for sufficiently small elastic deformations.”——能量中對位移為線性的項在平衡時會消失（見 Fig. 3.6 的極小值）；對於足夠小的彈性形變，三次與更高階的項可以忽略。
（書中指回去的 Fig. 3.6 是第 3 章 p.58 的圖，圖說寫的是 Lennard-Jones 位能的形狀：“Form of the Lennard-Jones potential (10) which describes the interaction of two inert gas atoms. The minimum occurs at R/σ = 2^{1/6} ≈ 1.12. Notice how steep the curve is inside the minimum, and how flat it is outside the minimum. The value of U at the minimum is −ε; and U = 0 at R = σ.”——注意曲線在極小值內側多陡、外側多平；極小值處 U = −ε，而 R = σ 處 U = 0。）
補出諧近似的代數：把 U 在平衡點展開，看清楚 C 到底是什麼補推導
Kittel 只說「線性項在平衡時消失、三次以上可略」，沒有把泰勒展開寫出來。補出來只要五行，但補出來之後，力常數 C 的身分會從「一個被引入的符號」變成「位能曲線在極小值的曲率」，而這正是後面能用實驗量到的 ω 反推 C 的理由。
- 設定：把晶體的總位能 U 看成所有原子面位移 {u_s} 的函數。平衡構型定義為 u_s = 0（對所有 s），這是 U 的一個極小點。
- 在極小點做泰勒展開：U({u}) = U₀ + Σ_s (∂U/∂u_s)₀ u_s + ½ Σ_{s,s′} (∂²U/∂u_s∂u_{s′})₀ u_s u_{s′} + （三次以上）。
- 零階項 U₀ 是平衡時的內聚能，是一個常數，對力沒有貢獻（力是 U 的梯度）。這一項是第 3 章的主題，第 4 章完全用不到——這就是為什麼鍵結的強弱（井有多深）與振動的頻率（井有多彎）是兩個獨立的問題。
- 一階項為零：在極小點，所有一階偏導數 (∂U/∂u_s)₀ = 0，這是極小值的定義。Kittel 指回 Fig. 3.6 就是要你看這件事——Lennard-Jones 曲線在 R/σ = 2^{1/6} 的地方切線水平。物理意義：平衡構型下每個原子受力為零；如果一階項不為零，原子就會自己跑掉，那個構型根本不是平衡構型。
- 二階項是唯一留下來的：U ≈ U₀ + ½ Σ_{s,s′} Φ_{ss′} u_s u_{s′}，其中 Φ_{ss′} ≡ (∂²U/∂u_s∂u_{s′})₀ 就是力常數矩陣。「彈性能是二次式」這句話的精確意思就是這一行。
- 由能量得力：F_s = −∂U/∂u_s = −Σ_{s′} Φ_{ss′} u_{s′}。力是位移的線性函數——這就是原文「elastic response is a linear function of the forces」與「is of the form of Hooke's law」兩句話的共同來源。能量的二次式 ⇔ 力的線性式，這兩件事在數學上是同一件事，Kittel 用 “That is equivalent to” 連接它們，用詞是精確的。
- 平移不變性的約束：把整個晶體剛性平移（所有 u_s 都加同一個常數 c），能量不能改變，所以每一行的力常數必須加總為零：Σ_{s′} Φ_{ss′} = 0。在只保留最近鄰的情形下這給出 Φ_{ss} = 2C 而 Φ_{s,s±1} = −C，代回去正好得到式 (1) 的 F_s = C(u_{s+1} − u_s) + C(u_{s−1} − u_s)。這解釋了式 (1) 裡為什麼是「差」而不是「絕對位移」：力只能依賴相對位移，否則整體平移會產生無中生有的力。
這一步換到的東西：C 現在有了一個與模型無關的定義——C 是位能在平衡點對「相鄰兩面相對位移」的二階導數，也就是位能井底的曲率。這件事有三個直接後果。
第一，C 與內聚能沒有固定關係。井深（內聚能）與井底曲率（C）是同一條曲線的兩個獨立讀數。這解釋了為什麼有些材料鍵很強卻很軟，有些相反。
第二，C 可以用實驗量到的 ω 反推回去，這正是章首目錄裡 “Derivation of force constants from experiment”（p.94）那一節要做的事，也是中子非彈性散射（p.100）在實務上的最大用途。
第三，諧近似的失效條件是「三次項不可忽略」，不是「位移很小」。這兩件事不完全等價：位移小但位能曲線的三次導數很大（強非諧材料），諧近似一樣會垮。
諧近似失效會怎樣：三個「純諧晶體不存在」的證據我的判讀
諧近似是這整章的立身之本，但它有一個非常尷尬的性質：如果它完全成立，固體的許多基本性質會是錯的。把這三件事記起來，你就知道第 5 章為什麼要專門開一節講非諧效應。
其一：熱膨脹係數會是零。 在嚴格的諧位能裡，位能井對平衡點左右對稱，原子振幅再大，時間平均的位置仍然停在井底，平均鍵長不隨溫度改變。真實晶體會膨脹，是因為位能井的外側比內側平（就是 Fig. 3.6 圖說特別要你注意的那件事：“Notice how steep the curve is inside the minimum, and how flat it is outside the minimum.”），振幅一大，平均位置就往外偏。熱膨脹是三次項的直接證據。
其二：聲子的壽命會是無限長。 諧系統的簡正模式彼此完全獨立——這正是後面能把每個 (K, 偏振) 當成一個獨立諧振子來量子化的前提。獨立的意思是：一個聲子被激發之後，永遠不會衰變，也不會散射到別的模式去。
其三：熱導率會發散。 這是前一項的推論，也是三者中最有工程意義的一個。若聲子永遠不散射，它們就像真空中的光子一樣自由穿越晶體，熱阻為零，熱導率無窮大。真實晶體的本徵熱導率有限；第 5 章 p.122 的式 (42) 就是用氣體動力論的形式處理這件事——原文說明該式中的 C 是單位體積的熱容、v 是平均粒子速度、另一個因子是粒子在兩次碰撞之間的平均自由程，並指出 Debye 最先把它套用到介電固體上，此時 C 取聲子熱容、v 取聲子速度、平均自由程取聲子的平均自由程。而限制平均自由程的機制在高溫下主要是聲子－聲子的 Umklapp 散射，那完全是三次項（三聲子過程）的產物——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.”
這一點直接關係到 FinFET 的熱阻萃取：奈米尺度元件的自加熱之所以嚴重，一半原因是邊界散射把 ℓ 壓到與鰭寬同量級（尺寸效應，與非諧性無關），另一半原因是熱點的局部溫度把非諧散射拉強（Umklapp 過程隨溫度增強，κ 大致隨 1/T 下降）。第 4 章建立的諧近似框架給出的是聲子有哪些模式、速度多少；把這些模式的壽命算出來，才是實際熱導的另一半。兩半缺一不可，而本段負責的是前一半。
建模假設三（原文）：只取最近鄰，且 C 是「對平面中一個原子」定義的
第三段假設在 p.91，原文逐字是：“We assume that 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 of their displacements. For brevity we consider only nearest-neighbor interactions, with p = ±1. The total force on s from planes s ± 1:”——我們假設由 s + p 面的位移在 s 面上造成的力，正比於兩者位移之差 u_{s+p} − u_s。為簡潔起見，我們只考慮最近鄰交互作用，p = ±1。
給出式 (1) 之後，Kittel 接了三句對 C 的說明，這三句每一句都要看清楚：
“This expression is linear in the displacements and is of the form of Hooke's law.”——這個表示式對位移是線性的，且具有虎克定律的形式。
“The constant C is the force constant between nearest-neighbor planes and will differ for longitudinal and transverse waves.”——常數 C 是最近鄰面之間的力常數，且對縱波與橫波會不同。
“It is convenient hereafter to regard C as defined for one atom of the plane, so that F_s is the force on one atom in the plane s.”——此後方便起見，把 C 視為是對平面中的一個原子定義的，如此 F_s 就是作用在 s 面中一個原子上的力。
三個關於 C 的陷阱，以及最近鄰假設什麼時候垮掉我的判讀
陷阱一：C 的「每原子」約定。 「C is defined for one atom of the plane」這句話看起來只是換單位，實際上是在為式 (2) 鋪路。式 (2) 的左邊是 M d²u_s/dt²，M 是一個原子的質量（原文下一行明說 “where M is the mass of an atom”）。左邊既然是一個原子的質量乘加速度，右邊就必須是一個原子受的力，兩邊才對得起來。如果你把 C 當成「整個面對整個面」的力常數，那右邊是整面的力，左邊卻是單原子的質量，量綱雖然還對得上（因為兩邊都可以再除以面內原子數），但你反推出來的 C 會差一個「面內原子數」的因子——而這個因子在不同方向、不同晶格是不同的。這是實務上用實驗色散曲線反推力常數時最常見的錯誤來源。
C 的單位是 N/m（力／長度），與一根彈簧的彈性係數同單位。這與第 3 章的彈性勁度 C₁₁（單位 N/m²，是體積性質）是不同的量，雖然符號都叫 C。兩者的橋樑是長波極限下的聲速，本段最後會實際做這個換算。
陷阱二：縱波與橫波的 C 不同。 原文寫得很直白 “will differ for longitudinal and transverse waves”。物理上很好理解：把兩個相鄰的面壓近或拉遠（縱向），抵抗來自鍵的伸縮，那是很硬的；把兩個面互相滑動（橫向），抵抗來自鍵角的改變，那通常軟得多。所以典型地 C_L > C_T，反映到色散曲線上就是縱支的最高頻率高於橫支，聲速也是縱波比較快。這與第 3 章的結果一致：立方晶體 [100] 方向縱波 v = (C₁₁/ρ)^{1/2}、橫波 v = (C₄₄/ρ)^{1/2}（p.82 式 (60) 與 (63)），而 Table 12（p.84）裡幾乎所有材料的 C₁₁ 都明顯大於 C₄₄——例如 Si 的 C₁₁ = 1.66、C₄₄ = 0.796（單位 10¹¹ N/m²）。所以這一整段推導其實要做三次（一縱兩橫），只是三次的數學形式一模一樣，只有 C 的數值換掉。
陷阱三：最近鄰假設在金屬裡是嚴重錯誤的。 這不是我的推測，書上自己在 p.94 就打臉了。“Derivation of Force Constants from Experiment” 一節的開頭逐字寫著：“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.”——在金屬中有效力可以有相當長的作用程，經由傳導電子海從離子傳到離子；已經發現相隔多達 20 個原子面之間仍有交互作用。
所以 p = ±1 是一個為了「簡潔（for brevity）」而做的計算便利，不是物理事實。一般化很直接：保留任意 p，式 (7) 會變成 ω² = (2/M)Σ_{p>0} C_p (1 − cos pKa) 的形式，也就是把色散曲線寫成傅立葉級數；每多一組 C_p，曲線就多一個諧波成分，可以描述更複雜的形狀（例如中間出現凹陷、或是 Kohn 異常那種尖點）。這正是本章習題 4「Kohn anomaly」（p.103）在練的事，也是實驗上用中子散射資料擬合力常數的標準做法。
反過來說：只保留最近鄰時，色散曲線被迫是 |sin ½Ka| 這一個形狀，沒有任何自由度可以調。整條曲線只有一個參數 C，它同時決定了聲速與最高頻率，兩者的比值被鎖死。真實材料若量到的聲速與最高頻率比值不符，就是在告訴你需要更遠程的力常數。
最近鄰模型下第 s 面所受的總力（虎克定律形式）
F_s = C(u_{s+1} − u_s) + C(u_{s−1} − u_s)
(1) p.91F_s：作用在第 s 面中一個原子上的力（依原文 p.91 的約定）；u_s：第 s 面偏離其平衡位置的位移（縱波時沿 K 方向，橫波時垂直於 K）；C：最近鄰面之間的力常數，單位 N/m，對縱波與橫波取值不同；下標 s ± 1：緊鄰的左右兩個面。
第 s 面中一個原子的運動方程式
M d²u_s/dt² = C(u_{s+1} + u_{s−1} − 2u_s)
(2) p.91M：一個原子的質量（原文 p.91：“where M is the mass of an atom”）；其餘符號同 (1)。右邊是把 (1) 的兩項展開合併後的結果：C(u_{s+1} − u_s) + C(u_{s−1} − u_s) = C(u_{s+1} + u_{s−1} − 2u_s)。
補齊 (2) → (3)：為什麼可以直接把 d²/dt² 換成 −ω²補推導
原文只寫 “We look for solutions with all displacements having the time dependence exp(−iωt). Then d²u_s/dt² = −ω²u_s, and (2) becomes …”。這一步初學者常常照抄，但其中有兩個值得講清楚的決定：為什麼可以分離變數，以及為什麼用複數指數而不用 cos。
- 分離變數的合法性：式 (2) 是線性且係數不含時間的微分方程組。這類方程的通解可以展開成 exp(−iωt) 這種諧時間相依的解的疊加（這就是對時間做傅立葉變換）。所以「只找 exp(−iωt) 形式的解」不會漏掉任何東西——這是完備性論證，不是近似。注意這一步用掉的是諧近似（線性）與系統不隨時間改變（係數為常數）這兩個前提。
- 寫成 u_s(t) = u_s exp(−iωt)（等號右邊的 u_s 現在只是與時間無關的振幅）。
- 求導兩次：d/dt[exp(−iωt)] = −iω exp(−iωt)；再一次得 (−iω)² = −ω²。所以 d²u_s/dt² = −ω² u_s exp(−iωt) = −ω² u_s(t)。這就是原文那句 d²u_s/dt² = −ω²u_s 的由來。
- 代回 (2)：左邊 M(−ω²u_s) exp(−iωt)，右邊 C(u_{s+1} + u_{s−1} − 2u_s) exp(−iωt)。時間因子 exp(−iωt) 在兩邊是公因子，直接消掉，因為它對任何 t 都不為零。
- 得到 (3)：−Mω²u_s = C(u_{s+1} + u_{s−1} − 2u_s)。方程式中的時間不見了——微分方程組變成了純代數的差分方程組。這是整段推導的第一個關鍵簡化。
為什麼用 exp(−iωt) 而不用 cos ωt：因為後面要用 exp(isKa) 這個空間因子，兩個指數函數相乘只是指數相加，代數會非常乾淨；如果用三角函數，(5) 那一步會變成一堆和角公式。物理上真實的位移當然是實數，正確的理解是：取複數解的實部，或者說把 ω 與 −ω 的兩個複數解配成共軛對相加。這是物理計算裡的標準操作，但要記得在最後計算能量、強度這類非線性的量之前先取實部，否則會出錯。
符號約定的一致性：Kittel 選 exp(−iωt) 搭配 exp(+isKa)，於是完整的波是 exp[i(sKa − ωt)]，這是沿 +x 方向傳播的波（相位固定要求 sKa − ωt = const，即 x 隨 t 增加）。若寫成 exp(+iωt)，同樣的 K > 0 會變成向 −x 傳播。這個約定貫穿全書，第 3 章 p.81 的式 (58) 與 p.82 的式 (61)、(64) 用的也是 exp(−iωt) 這個時間因子，一致。
消去時間後的差分方程
−Mω²u_s = C(u_{s+1} + u_{s−1} − 2u_s)
(3) p.91ω：角頻率（待求）；u_s：現在代表第 s 面的振幅（與時間無關的複數，含相位）；M、C 同前。
行進波試解（原文只寫出鄰面的形式）
u_{s±1} = u exp(isKa) exp(±iKa)
(4) p.92u：與 s 無關的共同振幅（一個複數常數，含整體相位）；K：波向量；a：垂直於 K 的相鄰原子面之間的距離，原文明說 “The value to use for a will depend on the direction of K.”；s：整數，面的編號。隱含的本體是 u_s = u exp(isKa)，而完整的解是 u_s(t) = u exp(isKa) exp(−iωt)。
為什麼指數函數一定是對的試解：離散平移對稱性（1D 版的布洛赫定理）補推導
Kittel 只說這個差分方程「有行進波解」，讀起來像是猜的。它不是猜的——指數函數是唯一可能的形式，理由完全來自對稱性。這一步補起來很值得，因為同一個論證在第 7 章講能帶時會原封不動地再用一次。
- 觀察對稱性：式 (3) 對每一個 s 都長得一模一樣——係數 M、C 與 s 無關。也就是說，把整條鏈的編號整體平移一格（s → s + 1），方程組回到自己。這叫離散平移不變性。
- 引入平移算符 T：定義 (Tu)_s ≡ u_{s+1}。上一步的意思就是 T 與運動方程的算符對易。線性代數的標準結論：兩個對易的算符可以同時對角化，所以我們可以只在 T 的本徵向量中尋找解。
- 解 T 的本徵方程：Tu = λu 即 u_{s+1} = λ u_s，遞迴得 u_s = λ^s u_0。
- 用有界性挑出 λ 的模：無限長的鏈上，物理解不能在 s → +∞ 或 s → −∞ 發散。|λ| > 1 則往右發散，|λ| < 1 則往左發散。所以必須 |λ| = 1，於是 λ 可以寫成 λ = exp(iKa)，其中 K 為實數，a 只是把無因次的相位換算成波向量時墊進去的長度尺度。
- 得到 (4)：u_s = u_0 λ^s = u exp(isKa)，而 u_{s±1} = λ^{±1}u_s = u exp(isKa) exp(±iKa)。這不是一個聰明的猜測，而是對稱性的必然結果。
- （有限鏈的版本）：若鏈有 N 個面並施加週期性邊界條件 u_{s+N} = u_s，則 λ^N = 1，於是 exp(iKNa) = 1，K = 2πn/(Na)，n 為整數。K 變成離散的，且在第一布里淵區內恰好有 N 個允許值。這個計數在第 5 章算態密度與熱容時是核心：N 個原子的鏈有 N 個模式，自由度守恆。
跨章連結：這整套論證（對稱性 → 平移算符 → 本徵值是模為 1 的相位因子 → 解具有 exp(iKx) × 週期函數 的形式）就是布洛赫定理。在本章因為每個原胞只有一個自由度，週期函數退化成常數 u，所以看不出「plane wave × periodic function」的結構；到了 p.95 的雙原子鏈，每個原胞有兩個振幅，週期部分就變成一個二分量向量，結構才顯形。完整的陳述要到書 p.167 的 “BLOCH FUNCTIONS” 一節才出現，Kittel 在那裡寫的是波函數形式 ψ_k(r) = u_k(r) exp(ik·r)，其中 u_k(r) 具有晶格的週期性。
值得記住的一句話：晶格振動與電子能帶之所以有那麼多平行的概念（布里淵區、色散曲線、區邊界的零斜率、能隙／頻隙），根源只有一個——兩者服從同一個離散平移對稱性。差別只在於波動方程一個是牛頓方程、一個是薛丁格方程。先在聲子這裡把幾何學弄熟，第 7 章會省下大量力氣。
把行進波試解代入差分方程
−ω²Mu exp(isKa) = Cu{exp[i(s + 1)Ka] + exp[i(s − 1)Ka] − 2 exp(isKa)}
(5) p.92各符號同 (3)、(4)。原文的引導語是 “With (4), we have from (3):”。
補齊 (5) → (6)：消去 u exp(isKa) 的完整代數補推導
原文對這一步只給了一句話：“We cancel u exp(isKa) from both sides, to leave”。這是本段最該補的一步，因為 (6) 右邊突然冒出一個負號，很多讀者在這裡對不起來。整個過程只有四行，但每一行都有值得說明的地方。
- 拆開右邊的指數：利用 exp(A + B) = exp(A)exp(B)，把 s ± 1 的指數拆成 s 的部分乘以 ±1 的部分。 exp[i(s + 1)Ka] = exp(isKa)·exp(iKa) exp[i(s − 1)Ka] = exp(isKa)·exp(−iKa)
- 在右邊把 exp(isKa) 提出來： 右邊 = Cu·exp(isKa)·{exp(iKa) + exp(−iKa) − 2} （第三項 −2exp(isKa) 本身就含 exp(isKa)，提出來之後留下 −2。）
- 兩邊對照： −ω²M·[u exp(isKa)] = C·[u exp(isKa)]·{exp(iKa) + exp(−iKa) − 2} 兩邊的中括號完全相同。
- 確認可以消：u ≠ 0（否則是振幅為零的平庸解，沒有振動）；exp(isKa) 對任何實數 sKa 的模都等於 1，永遠不為零。所以公因子可以安全地除掉： −ω²M = C[exp(iKa) + exp(−iKa) − 2]
- 移項成書上的寫法：兩邊同乘 −1，把負號從左邊搬到右邊： ω²M = −C[exp(iKa) + exp(−iKa) − 2] ……這就是式 (6)。 那個讓人困惑的負號，來源單純就是「左邊原本是 −ω²M」。Kittel 選擇把 ω²M 留在左邊乾淨的位置，代價是右邊掛一個負號。
- 檢查符號是否合理：方括號裡是 2cos Ka − 2，這個量恆為負（因為 cos Ka ≤ 1）；前面的負號把它翻正，所以 ω²M ≥ 0，ω² ≥ 0。若這裡算出負的 ω²，就是代數出錯了——這是最好的自我檢查點。
這一步真正的份量：無限多條耦合方程（每個 s 一條）在這裡壓成一條純代數式，而且未知數只剩 ω 與 K 兩個。振幅 u 完全消失了——這是線性系統的招牌特徵：頻率與振幅無關。諧振子不管盪得多大，週期都一樣；諧近似下的晶格也是如此。（一旦保留非諧項，頻率就會開始依賴振幅，這是非諧效應最直接的實驗簽名之一。）
面編號 s 也完全消失了，意思是：這個解對鏈上的每一個面都同時成立，它是一個真正的集體模式，不是某個局部原子的運動。這正是「簡正模式（normal mode）」的定義。
消去公因子後的代數關係
ω²M = −C[exp(iKa) + exp(−iKa) − 2]
(6) p.92只剩 ω、K 兩個變數與 M、C 兩個參數。振幅 u 與面編號 s 都已消去。
補齊 (6) → (7)：歐拉恆等式那一步，以及為什麼結果必然是實數補推導
原文寫 “With the identity 2 cos Ka = exp(iKa) + exp(−iKa), we have the dispersion relation ω(K).” 代數只有兩行，但「複數指數怎麼變回實數」這件事本身值得說明一句。
- 歐拉公式：exp(iθ) = cos θ + i sin θ，exp(−iθ) = cos θ − i sin θ。兩式相加，虛部相消：exp(iθ) + exp(−iθ) = 2cos θ。取 θ = Ka 就是原文引用的恆等式。虛部之所以必然相消，是因為 (6) 右邊的 +Ka 與 −Ka 兩項是成對出現的——而它們成對出現，又是因為每個面的左右兩邊各有一個最近鄰。不過真正保證 ω² 是實數的，並不是晶體有沒有反演中心，而是力常數矩陣本身的對稱性：二階混合偏導數可交換給出 Φ_{ss′} = Φ_{s′s}，再加上平移不變性就得到 Φ_p = Φ_{−p}（本模型即 Φ_{s,s+1} = Φ_{s,s−1} = −C），於是動力學矩陣是厄米的，ω² 必為實數；Φ_p = Φ_{−p} 同時讓它是實矩陣，這才是 +Ka 與 −Ka 成對出現的來源。
- 代入 (6)：ω²M = −C[2cos Ka − 2] = −2C cos Ka + 2C = 2C(1 − cos Ka)。
- 兩邊除以 M：ω² = (2C/M)(1 − cos Ka) ……這就是式 (7)。
現在這條式子的物理內容第一次完全顯露：ω² 正比於 (1 − cos Ka)，也就是正比於「相鄰兩面位移的相位差所造成的相對位移大小」。Ka = 0 時相鄰面同相，沒有相對位移，回復力為零，ω = 0。Ka = π 時相鄰面反相，相對位移最大，回復力最強，ω 最大。整條色散曲線在說的就只有這一件事：頻率高低完全由「相鄰面被拉扯得多厲害」決定。
Kittel 就是在引出這一步的那句話裡寫下 “we have the dispersion relation ω(K)”——該詞印在式 (7) 之前；而且它並不是全書首次登場，第 3 章 p.82 已經先用在無色散的 (62) ω²ρ = C₄₄K² 上。不過到了 (7) 這個詞才名副其實：ω² 與 K² 不再成正比，波速開始依賴波長，色散正式登場。
單原子鏈的色散關係（餘弦形式）
ω² = (2C/M)(1 − cos Ka)
(7) p.92ω：角頻率；K：波向量；a：面間距；C：最近鄰面間力常數（每原子）；M：原子質量。這是第一次出現完整的 ω(K) 關係，原文在引出它時明白稱之為 “the dispersion relation ω(K)”。
ω² 對 K 的斜率，及其在區邊界為零
dω²/dK = (2Ca/M) sin Ka = 0 （在 K = ±π/a）
(8) p.92dω²/dK：ω 平方對 K 的導數（注意不是 dω/dK）；其餘同 (7)。原文的前導句是 “The boundary of the first Brillouin zone lies at K = ±π/a. We show from (7) that the slope of ω versus K is zero at the zone boundary:”，後續句是 “at K = ±π/a, for here sin Ka = sin (±π) = 0.”
式 (8) 的一個陷阱：dω²/dK = 0 不一定代表 dω/dK = 0補推導
Kittel 用 (8) 論證「ω 對 K 的斜率在區邊界為零」，但他微分的是 ω² 而不是 ω。這一步在區邊界是對的，在 K = 0 卻是錯的——而書上沒有提醒這件事。這個細節值得講，因為它同時解釋了為什麼色散曲線在原點是有斜率的尖角、在區邊界卻是平滑的水平切線。
- 連鎖律：dω²/dK = 2ω·(dω/dK)。因此 dω²/dK = 0 只能推出「ω = 0 或 dω/dK = 0」，兩者取其一。
- 在區邊界 K = ±π/a：此時 ω = (4C/M)^{1/2} ≠ 0，所以必然是 dω/dK = 0。Kittel 的論證成立。直接驗算：由 (9) ω = (4C/M)^{1/2}|sin ½Ka|，dω/dK = (4C/M)^{1/2}·(a/2)cos(½Ka)·sgn(...)，在 Ka = π 時 cos(π/2) = 0，確實為零。✓
- 在 K = 0：sin(0) = 0，所以 dω²/dK 一樣是零——但這裡 ω = 0，所以連鎖律給不出任何關於 dω/dK 的結論。實際上此處 dω/dK = ±a(C/M)^{1/2} = ±v_s，等於聲速，明顯不為零。
- 幾何圖像：ω(K) = (4C/M)^{1/2}|sin ½Ka| 在 K = 0 是一個尖角（V 字形），左右斜率為 ∓v_s 與 ±v_s；在 K = ±π/a 則是一個光滑的極大值，切線水平。Fig. 4 把這兩個特徵都畫得很清楚：原點是尖的，π/a 處是圓的。
為什麼要在意：因為這兩個點對應兩種完全不同的物理。
K = 0 的尖角 是聲學支的標誌：ω ∝ |K|，群速等於相速等於聲速，波包不散開，這是聲音之所以能保真傳播的原因。尖角來自絕對值——嚴格說 ω(K) 在 K = 0 是不可微的，但兩側的單邊導數各自有定義，分別對應向左與向右傳播的波。
K = π/a 的水平切線 是駐波的標誌：群速為零，能量原地打轉。在熱傳導上這件事有直接後果——接近區邊界的聲子幾乎不輸送熱量，雖然它們在高溫下佔據數最多。這解釋了一個常被誤解的現象：高溫下熱容飽和（每個模式都被激發滿了），熱導率卻還在下降，因為主要被激發的那些高 K 模式群速小、又被非諧散射得厲害。做熱阻萃取時，如果用 Debye 模型的單一常數群速去估，在高溫段會系統性地高估熱導。
補齊 (7) → (9)：半角恆等式那一步（書上只寫 “By a trigonometric identity”）補推導
原文在這裡只留下一句 “By a trigonometric identity, (7) may be written as”，完全沒說是哪一個恆等式。這一步是標準結果，但書上沒逐步寫，所以在這裡補齊——連同「為什麼最後會冒出一個絕對值」一起。
- 要用的恆等式：1 − cos θ = 2 sin²(θ/2)。證明：二倍角公式 cos(2φ) = cos²φ − sin²φ = (1 − sin²φ) − sin²φ = 1 − 2sin²φ。令 φ = θ/2，即得 cos θ = 1 − 2sin²(θ/2)，移項就是 1 − cos θ = 2sin²(θ/2)。
- 代入 (7)，取 θ = Ka： ω² = (2C/M)(1 − cos Ka) = (2C/M)·2sin²(½Ka) = (4C/M) sin²(½Ka) 這是 (9) 的前半。
- 開平方：ω = ±(4C/M)^{1/2}·|sin(½Ka)|。注意開平方時 √(sin²x) = |sin x| 而不是 sin x——當 −2π/a < K < 0 時 sin(½Ka) 是負的。
- 挑正根：物理上頻率取非負（ω ≥ 0）。負頻率解並沒有消失，它只是與正頻率解描述同一個實數運動：exp[i(sKa − ωt)] 與 exp[−i(sKa − ωt)] 互為複共軛，相加才是實位移。所以把 ω ≥ 0 定死，方向資訊全部交給 K 的正負去攜帶。於是 ω = (4C/M)^{1/2}|sin ½Ka| ……這就是 (9) 的後半，書上確實印了絕對值符號。
為什麼要特地換成半角形式，明明 (7) 已經是完整的色散關係？有三個理由。
第一，最大值一眼可見。|sin| 的最大值是 1，所以 ω_max = (4C/M)^{1/2}，不需要再做微分。這個量是整章最重要的一個數字：它是這個晶格能支撐的最高振動頻率，第 5 章的 Debye 截止頻率就是它的量級。
第二，週期性與對稱性一眼可見。|sin ½Ka| 對 K 是偶函數（ω(−K) = ω(K)；這是時間反演對稱的普遍結果，對所有晶體都成立，連非中心對稱晶體的聲子色散也是偶函數），週期是 2π/a（因為半角使 sin 的週期變成 4π/a，但絕對值把它折半回 2π/a）。這兩個性質決定了為什麼只畫 0 到 π/a 就夠了。
第三，Fig. 4 的縱軸就是這條式子的無因次化。圖上標的是 ω/(4C/M)^{1/2} = |sin ½Ka|，一條與所有材料參數無關的普適曲線。所有單原子最近鄰鏈的色散曲線長得一模一樣，材料只決定縱軸的尺度 (4C/M)^{1/2} 與橫軸的尺度 π/a。這是這個模型最漂亮也最侷限的地方——漂亮是因為普適，侷限是因為真實材料的色散曲線形狀各不相同，要描述那些差異就必須放進更遠程的力常數。
單原子鏈的色散關係（半角形式，本段的終點）
ω² = (4C/M) sin²(½Ka) ; ω = (4C/M)^{1/2} |sin ½Ka|
(9) p.92(4C/M)^{1/2}：最高角頻率 ω_max，發生在 Ka = π；|sin ½Ka|：無因次的形狀函數，只依賴相位 Ka，與材料無關。原文在這條式子之後只補了一句 “A plot of ω versus K is given in Fig. 4.”
Fig. 4 的文字還原（書 p.92）：一條普適曲線
圖說逐字為：“Plot of ω versus K. The region of K ≪ 1/a or λ ≫ a corresponds to the continuum approximation; here ω is directly proportional to K.”
圖面的內容（依原書算繪後目視判讀）：
縱軸標的不是 ω 本身，而是無因次比值 ω/(4C/M)^{1/2}，刻度從 0 到 1.2，每 0.2 一格。這樣畫的用意是把材料參數全部吸收進去，讓曲線變成普適的。
橫軸是 K，畫出的範圍從 −π/a 經過 0、π/a，一直延伸到 2π/a，總共涵蓋一個半週期。
曲線：在 K = −π/a 處值為 1.0（頂點），往右一路下降，在 K = 0 處觸底為 0，而且是一個尖角（V 字形），不是圓滑的谷底；接著往右上升，在 K = π/a 處達到最大值 1.0，這裡是圓滑的極大值，切線水平；再往右又下降，在 K = 2π/a 回到 0。也就是說，π/a 右邊那一段與左邊那一段完全是鏡像——曲線是週期為 2π/a 的函數，π/a 到 2π/a 這一段沒有攜帶任何新資訊，這正是下一段（p.93）要論證的「只有第一布里淵區有物理意義」。
圖中央偏右處寫著曲線的方程式：ω/(4C/M)^{1/2} = |sin ½Ka|，絕對值符號在圖上是明確畫出來的。
圖的底部有一個左右帶箭頭的水平線段，從 −π/a 拉到 π/a，中間標註 “First Brillouin zone”。在 K = 0 與 K = π/a 的位置各有一條細的垂直參考線。
這張圖要你看見三件事：(1) 原點附近是直線（連續近似成立），(2) 越靠近區邊界曲線越平（群速趨近零），(3) 超過區邊界的部分只是重複。三件事都會在接下來兩頁被正式處理。
數值落地：用書上 Table 11 的鋁資料反推力常數 C 與最高頻率補推導
這一整段推導的參數只有 C 與 M 兩個，但書上到這裡一個數字也沒給。把數字代進去，色散曲線才會從一條符號式變成你對材料的直覺。下面用第 3 章 Table 11（書 p.84）真的印出來的鋁資料，走完整條換算鏈：彈性勁度 → 聲速 → 力常數 → 最高頻率 → 等效溫度。
- 輸入（來自書本 p.84 Table 11）：鋁在 0 K 的 C₁₁ = 1.143、300 K 的 C₁₁ = 1.068，單位是 10¹² dyne/cm²（= 10¹¹ N/m²）；密度 ρ = 2.733 g/cm³（表中只在 0 K 那一列印出密度，300 K 列的密度欄是「—」）。取室溫值 C₁₁ = 1.068×10¹¹ N/m²、ρ = 2733 kg/m³。
- 輸入（不在本段頁面上，屬於一般常數）：鋁的原子量 26.98 → M = 26.98 × 1.6605×10⁻²⁷ = 4.48×10⁻²⁶ kg；鋁是 fcc，立方晶格常數 a₀ = 4.05 Å。
- 第一步：a 取多少？ K ∥ [100]，fcc 沿 [100] 的原子面間距是 a = a₀/2 = 2.025 Å = 2.025×10⁻¹⁰ m，不是 a₀。（理由見前面的面間距表：fcc 的面心原子在名目上的兩層之間再插了一層。這正是原文 “The value to use for a will depend on the direction of K.” 在實務上的意思。）
- 第二步：聲速。由第 3 章 p.82 的式 (60)，[100] 縱波 v = (C₁₁/ρ)^{1/2} = (1.068×10¹¹ / 2733)^{1/2} = (3.908×10⁷)^{1/2} = 6.25×10³ m/s。
- 第三步：由長波極限反推 C。本段的 (7) 在 Ka ≪ 1 給出 v = a(C/M)^{1/2}，反解得 C = M v²/a²。 C = 4.48×10⁻²⁶ × (6.25×10³)² / (2.025×10⁻¹⁰)² = 4.48×10⁻²⁶ × 3.906×10⁷ / 4.101×10⁻²⁰ = 1.750×10⁻¹⁸ / 4.101×10⁻²⁰ = 42.7 N/m
- 第四步：最高頻率。由 (9)，ω_max = (4C/M)^{1/2} = (4 × 42.7 / 4.48×10⁻²⁶)^{1/2} = (3.81×10²⁷)^{1/2} = 6.17×10¹³ rad/s，即 ν_max = ω_max/2π = 9.8×10¹² Hz = 9.8 THz。
- 第五步：等效溫度。ℏω_max/k_B = (1.0546×10⁻³⁴ × 6.17×10¹³)/1.3807×10⁻²³ = 6.51×10⁻²¹/1.3807×10⁻²³ = 471 K。
- 第六步：區邊界波長。λ_min = 2a = 4.05 Å，正好等於一個 fcc 立方晶格常數。
這組數字告訴你什麼：
C ≈ 43 N/m 是原子間力常數的尺度（再提醒一次：依 p.91 的約定，它是一個原子對相鄰整面的等效勁度，不是單一化學鍵）。典型固體的原子間力常數落在 10–100 N/m，這個數字落在中間——鋁是軟金屬，所以偏下緣。作為對比：Table 12（p.84）給鑽石 C₁₁ = 10.76，而 Table 11（p.84）給鋁 300 K 的 C₁₁ = 1.068，兩者同為 10¹¹ N/m² 的單位，鑽石硬了十倍；再加上碳的原子量只有鋁的 44%（12.01 對 26.98），由 ω ∝ (C/M)^{1/2} 可知鑽石的聲子頻率會高出好幾倍。這與 p.116 Table 1 的目視判讀完全一致：鑽石 θ = 2230 K，是全表最高；矽 θ = 645 K；鋁 θ = 428 K。
與書上 Debye 溫度的比對：本模型算出鋁的 ℏω_max/k_B = 471 K，書上 p.116 給 428 K，差 10%。用一個只有一個可調參數、而且那個參數還是由聲速定死的模型，能把 Debye 溫度估到 10% 以內，這是這個玩具模型最有說服力的一次驗證。
9.8 THz 與 471 K 的量級要記住。它們說明：聲子的能量尺度是數十 meV（ℏω_max ≈ 41 meV），而電子在半導體裡的能隙是 eV 等級（Si 約 1.1 eV）。兩者差了 25 倍以上——這個比例正是為什麼電子與晶格可以近似地分開處理（Born–Oppenheimer 近似），也是為什麼電子跨過間接能隙時必須同時借助聲子提供動量（所以間接能隙材料的吸收係數遠小於直接能隙材料）。
一個誠實的檢查：這個模型只用「聲速對得上」這一個條件去定 C，剩下的全是預測。但要注意它的自由度有多低——只保留最近鄰時，聲速與最高頻率的比值被鎖死為 ω_max/(v·K_max) = 2/π ≈ 0.64，任何材料都一樣，沒有任何調整空間。對某些材料（特別是有 Kohn 異常的金屬、或共價鍵方向性很強的材料）這個比值明顯不對，那就是在告訴你必須放進更遠程的力常數——正如 p.94 所說，金屬中的交互作用可以延伸到 20 個原子面之遠。
五個容易踩的誤解我的判讀
誤解一：「u_s 是第 s 個原子的位移」。 不是。u_s 是第 s 個原子面的位移，整個面上所有原子共用這一個數。這個等同只在 K 沿高對稱方向、且整面同相時成立。忘了這件事，在算「一個模式對應幾個自由度」時就會數錯——N 個面給出 N 個 K 值，而不是 N 個原子。
誤解二：「ℏK 是聲子的動量」。 本章 p.100 有一節就叫 “PHONON MOMENTUM”（見章首目錄，書 p.89），本段還沒碰到，但這裡先埋一個提醒：K 只在模 2π/a 的意義下有定義（因為 exp(i2πns) = 1），所以 ℏK 不可能是真正的動量——真正的動量是唯一確定的，不會「等價於另一個值」。ℏK 是晶體動量（crystal momentum），它在散射過程中守恆到差一個倒晶格向量為止（Umklapp 過程），而這個「差一個倒晶格向量」正是聲子熱導有限的關鍵機制。
誤解三：「K > π/a 的波不存在」。 準確的說法是：它存在，但它不是新的東西。p.93 的式 (11) 會證明，K 與 K − 2πn/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.” 這在數位訊號處理裡就是取樣定理（aliasing）：用間距 a 的格點取樣，任何比 2a 短的波形都會被混疊回長波。離散取樣 ⇒ 頻譜週期化 ⇒ 布里淵區，三者是同一件事。
誤解四：「色散是因為原子之間有交互作用」。 不是。連續彈性介質裡原子之間也有交互作用（那就是彈性常數的來源），但它沒有色散。色散的唯一來源是離散性——是「波只在間距 a 的點上取樣」這件事。證據很直接：把 (7) 取 a → 0 的極限（同時固定聲速 v = a(C/M)^{1/2}，也就是讓 a²C/M 保持不變），對任何固定的 K 都有 Ka → 0，於是 ω² → (C/M)a²K² = v²K²，色散完全消失，回到連續介質的直線。a 才是色散的唯一開關。
誤解五：「(9) 的絕對值只是為了好看」。 不是。開平方得到的是 √(sin²x) = |sin x|，而在 −2π/a < K < 0 的區間 sin(½Ka) 確實為負。若漏掉絕對值，會得到負頻率，那在物理上沒有意義（或者說，它與正頻率解描述的是同一個實運動）。Fig. 4 的圖上也明確畫出了絕對值符號。
接到 FinFET 熱阻萃取：這九條式子為什麼跟 Cth 有關我的判讀
本段建立的東西看起來離元件熱阻很遠，但它其實是整條鏈的第一環。把關係理一次：
元件熱容 Cth 的微觀來源就是聲子的模式數。 第 5 章會把熱容寫成對所有 (K, 偏振) 模式求和，每個模式貢獻一個 Einstein 函數。要求和，必須先知道有哪些模式、頻率各是多少——那正是 ω(K) 這條曲線在做的事。本段的 (9) 是這條曲線最簡單的版本。在 Debye 近似裡，這條真實的 |sin ½Ka| 曲線被粗暴地換成直線 ω = vK 再加上一個截止頻率，而截止頻率的量級正是本段算出來的 (4C/M)^{1/2}。用 Debye 模型估 Cth 時，你其實是在用本段這條曲線的兩個端點（斜率與最大值）去代表整條曲線。 知道它是怎麼被簡化的，就知道誤差會出現在哪個溫度區間（中溫段最糟，因為那裡被激發的模式正好落在曲線彎曲最厲害的地方）。
熱阻 Rth 需要的是群速與壽命，而群速直接來自本段。 v_g = dω/dK 在 K = 0 是聲速、在區邊界是零。這個變化在奈米元件裡特別重要：當鰭寬縮到 5–10 nm，邊界散射把平均自由程壓到與鰭寬同量級，此時不同 K 的聲子貢獻差異被放大——長波聲子（群速大但波長長，容易被邊界散射）與短波聲子（群速小但波長短，比較不受邊界影響）的競爭決定了有效熱導。用單一常數群速的 Debye 圖像會系統性地估錯這個競爭。
尺寸效應的判準也在本段。 本段給出的最短波長是 λ_min = 2a（矽約 3–5 Å），另一端是元件尺寸。（以下兩個數字不在 Kittel 書內，屬元件熱傳文獻的常見量級，僅供對照）：矽在室溫下主要熱載聲子的波長約 1–3 nm、平均自由程約 100–300 nm。鰭寬 5 nm 已經遠小於平均自由程、但仍大於主要聲子波長——這個區間通稱「準彈道輸運」，傅立葉定律失效，但波動效應（干涉、侷限）還沒完全主導。萃取 Cth 與 Rth 時如果沿用塊材熱導率，誤差可以到數倍。Kittel 自己在 p.126 也點到同一類機制：“Geometrical effects may also be important in limiting the mean free path. We must consider scattering by crystal boundaries”。本段的意義在於：它讓你知道 a、λ_min、K_max 這幾個長度尺度分別落在哪裡，才能判斷手上的元件屬於哪個輸運區間。
最後一個實務提醒：矽不是單原子基底——金剛石結構每個原胞有兩個原子，所以矽的色散關係要用 p.95 起的雙原子鏈模型，有聲學支也有光學支。光學支的群速幾乎為零，幾乎不輸送熱，但它承載可觀的熱容。這就是為什麼在矽裡「熱容大的模式」與「導熱的模式」不是同一批——這個分離對暫態熱分析（也就是 Cth 萃取真正在量的東西）影響很大。本段的單原子鏈是理解那件事的必要準備，但不是終點。