06體模量與壓縮率；立方晶中的彈性波；本章摘要與習題

這一段在解決什麼問題我的判讀

第 3 章走到 p.80 的時候，彈性理論的骨架已經搭完了：應變六個分量 (31)(32)、應力六個分量、Hooke 定律 (37)(38) 把兩者用 36 個係數連起來，能量密度是應變的二次式 (39)，對稱性把 36 個降到 21 個 (42)，立方對稱再降到 3 個，寫成矩陣就是 (50)。

問題是：到此為止 C₁₁、C₁₂、C₄₄ 只是三個抽象的比例常數。它們既不是你在實驗室裡直接量的東西，也還看不出跟後面幾章要講的聲子有什麼關係。這一段做的正是把它們接回現實的兩件事。

第一件是靜態的：拿一個均勻的等向膨脹去打 (43)，得到體模量 B = ⅓(C₁₁ + 2C₁₂)。體模量是壓力與體積變化的關係，可以用靜水壓實驗或狀態方程直接量，Kittel 在 p.80 明說 B 與壓縮率 K 的數值列在 Table 3。

第二件是動態的：寫出彈性波的運動方程，解出三個主要方向的相速度。這一件在物理上重要得多，因為它是整本書往下走的橋——第 4 章的聲子色散關係在長波極限必須回到這裡算出的聲速，第 5 章的 Debye 模型則直接把「每個偏振的聲速取成常數」當作出發點。也就是說，這一段算出的 v = (C/ρ)^(1/2) 就是 Debye 溫度、晶格熱容、熱導率裡那個 v。

還有第三個副產品：能量密度 (43) 必須對所有應變組合為正，這個要求會把三個常數限制在一塊區域裡（習題 14）。而那塊區域的邊界，恰恰就是某個波速變成零的地方——晶格軟化、結構相變的數學前身。這條線索在這一段末尾才會收攏，但它是整段最漂亮的地方。

為什麼要在意這兩個本徵值：三個常數其實是「一個體積 + 兩個剪切」我的判讀

上面那段代數不只是技巧。它告訴你一件觀念上的事：立方晶的三個獨立彈性常數，在物理上其實是三個獨立的「彈簧」，各自對應一種形變模式。

第一種是等向膨脹（本徵值 C₁₁ + 2C₁₂，一重）。整個晶體均勻脹大或縮小，形狀不變、體積變。這條彈簧的硬度就是 3B，因為 B = ⅓(C₁₁ + 2C₁₂)。

第二種是四方剪切（本徵值 C₁₁ − C₁₂，二重）。晶體沿 x 拉長、沿 y 壓扁，體積不變、形狀變。對應的有效剪切模量是 ½(C₁₁ − C₁₂)，這就是習題 11 要你證的東西，後面也會看到它正是 [110] 方向第二支橫波的有效彈性常數。之所以是二重，是因為在保體積的三維空間裡，你有兩個獨立方向可以「拉長一軸壓扁另一軸」。

第三種是立方軸剪切（本徵值 C₄₄，三重）。把立方體的一個面沿著另一個立方軸推歪，也是保體積。

第二種與第三種都是剪切、都保體積，但它們的硬度不見得相同。它們的比值定義了 Zener 各向異性比 A ≡ 2C₄₄/(C₁₁ − C₁₂)（這個名字與定義是我補的，Kittel 這幾頁沒有提）。A = 1 時晶體在彈性上完全各向同性，這時三個常數退化成兩個，也就是連續介質彈性力學裡的 Lamé 常數。A 離 1 越遠，晶體的彈性越各向異性——這是後面看 Table 11、Table 12 時最該注意的一欄「隱藏欄位」。

體模量：為什麼要選「均勻膨脹」這個特定的變形

Kittel 的處理極為簡潔：考慮均勻膨脹 e_xx = e_yy = e_zz = δ/3，把它代進立方晶的能量密度 (43)，算出能量 (52)，再跟體模量的定義式 (53) 對照，就讀出 B。

為什麼選這個變形？因為 (35) 已經定義了膨脹（dilation）δ ≡ (V′ − V)/V = e_xx + e_yy + e_zz。要讓總膨脹等於 δ 而且不引入任何形狀改變，最自然的選法就是三個法向應變平分，各取 δ/3，剪切應變全部設成零。這樣一來 (43) 裡的 C₄₄ 項自動消失——體模量與 C₄₄ 無關，這是立方對稱的直接結果，也是為什麼 (54) 只含 C₁₁ 與 C₁₂。

從上一段的本徵值語言看，這個選擇一點也不隨意：e = (δ/3)(1, 1, 1) 正好是法向區塊那個一重本徵向量的方向。沿本徵方向做形變，能量密度就是「½ × 本徵值 × 振幅²」，不會有交叉項——這才是 (52) 能一行寫完的真正原因。

書上接著說壓縮率 K 定義為 K = 1/B，並指出 B 與 K 的數值列在 Table 3。這裡要提醒一個抽取上的事實：Table 3 在書本 p.52，是一張旋轉 90° 的週期表形式的表，該頁沒有文字層（我用 PyMuPDF 算繪後目視判讀）。表頭寫的是等溫體模量與壓縮率，體模量的單位是 10¹² dyn/cm² 或 10¹¹ N/m²，壓縮率的單位是 10⁻¹² cm²/dyn 或 10⁻¹¹ m²/N，跟 Table 11、Table 12 的剛度常數同單位，可以直接對照。

Fig. 18 在說什麼：從「力」到「波動方程」的那一步

網頁上不放圖，所以把 Fig. 18 用文字說清楚。

畫面是一個邊長 Δx、Δy、Δz 的小立方體，切在晶體內部任意位置。看它在 x 方向受到的力。左邊那個面（位置 x）上有應力 −X_x(x)（負號是因為外法線指向 −x）；右邊那個平行面（位置 x + Δx）上有應力 X_x(x + Δx)。關鍵的一步是把右面的應力做泰勒展開：

X_x(x + Δx) ≅ X_x(x) + (∂X_x/∂x)Δx

兩個面的貢獻相減，X_x(x) 這一項精確抵消，只剩下梯度項乘上面積：淨力 = (∂X_x/∂x)Δx × ΔyΔz。這就是整個彈性波理論的核心洞見——推動物質的不是應力本身，而是應力的梯度。

除了這一對面之外，x 方向的力還有另外兩個來源：垂直於 y 的那一對面上的剪應力 X_y（x 方向的力作用在法線沿 y 的面上），與垂直於 z 的那一對面上的 X_z。圖上沒有畫這兩對，書裡用文字補充說明。三者加起來，作用在立方體上的 x 方向淨力是

F_x = (∂X_x/∂x + ∂X_y/∂y + ∂X_z/∂z) ΔxΔyΔz

最後套牛頓第二定律：力等於質量乘上 x 方向的加速度。立方體的質量是 ρΔxΔyΔz，加速度是 ∂²u/∂t²。兩邊的 ΔxΔyΔz 對消，就得到 (55)。ΔxΔyΔz 會對消這件事很重要：它代表結果與你切多大的體積元無關，也就是這是一條逐點成立的場方程，而不是某個特定小盒子的力平衡。

Fig. 19 的兩根彈簧：為什麼是梯度而不是應力

Fig. 19 是一張極簡的圖，但它承載的觀念比 Fig. 18 還基本，值得原樣講。

畫面是一個方塊，左右各接一根彈簧，左邊標 A、右邊標 B，兩根沿著 x 方向。書上的敘述分兩種情況。

情況一：A 與 B 被拉伸的量相同。這時方塊受到左右兩個大小相等、方向相反的拉力，淨力為零，方塊不動。書上用這個來說明一件事：固體裡均勻的應力 X_x 不會對體積元產生淨力。

情況二：B 比 A 拉得更多。這時方塊會被 X_x(B) − X_x(A) 這個力加速。

這張圖解釋了為什麼 (55) 右邊是應力的導數而不是應力本身。換個角度想：如果應力本身就能產生加速度，那麼一塊靜置在桌上被均勻大氣壓包圍的晶體就會自己動起來，顯然荒謬。只有應力在空間上不均勻時，兩側才會失衡，物質才會被推動。

（以下這兩個類比、以及最後偏振數目的對照，都是我補的；Kittel 這幾頁沒有寫，書上 Fig. 19 的圖說就只有上面那兩句。）這也是一個好的類比起點：電荷在均勻電位裡不受力，受力的是電位的梯度；流體不會因為壓力大而流動，會流動是因為壓力有梯度。彈性波在這個意義上就是「應力不均勻度」在晶體裡的傳播。但類比到此為止——彈性波與電磁波的本質差別在於，彈性波有三個偏振（兩橫一縱），而電磁波在真空中只有兩個橫偏振。多出來的縱波正是因為固體能抵抗壓縮。

(57b)(57c)：為什麼可以「由對稱性直接寫出」

書上 p.81 開頭說 v 與 w 的運動方程可以由 (57a) 用對稱性直接得到，然後就把 (57b)(57c) 印出來。這句話值得拆開講，因為它用的是一個非常具體的對稱操作，不是含糊的「由對稱性」。

用的是 (45) 裡那個三重旋轉：x → y → z → x。這是立方晶必然具備的對稱操作（書上 p.78 說立方結構的最低對稱要求就是四條三重旋轉軸，沿 [111] 及等價方向）。在這個代換下，位移分量跟著輪換 u → v → w → u，而 C₁₁、C₁₂、C₄₄ 因為是立方不變量所以不變。

把 (57a) 整條做一次輪換：

ρ ∂²v/∂t² = C₁₁ ∂²v/∂y² + C₄₄(∂²v/∂z² + ∂²v/∂x²) + (C₁₂ + C₄₄)(∂²w/∂y∂z + ∂²u/∂y∂x)

把它跟書上的 (57b) 比對：ρ ∂²v/∂t² = C₁₁ ∂²v/∂y² + C₄₄(∂²v/∂x² + ∂²v/∂z²) + (C₁₂ + C₄₄)(∂²u/∂x∂y + ∂²w/∂y∂z) —— 完全一致（項的順序不同而已，且混合偏導可交換）。再輪換一次就得到 (57c)。

這裡有一個實用的檢查法：(57b) 的第一項一定是 C₁₁ ∂²v/∂y²（位移分量與微分方向同名配 C₁₁），C₄₄ 配的是另外兩個方向的二階導數，而耦合項一定是「(C₁₂ + C₄₄) × 另外兩個位移分量對『自己方向 × 對方方向』的混合導數」。用這個規則寫錯的機會很小。

[100] 方向：最乾淨的一個方向

從這裡開始 Kittel 尋找 (57) 的簡單特解。第一個方向是 [100]，也就是波向量沿著立方體的稜。

縱波：取 u = u₀ exp[i(Kx − ωt)]（式 58），v = w = 0。書上強調「波向量與粒子運動都沿著 x 立方稜」。代進 (57a)：因為 v = w = 0，耦合項整條消失；因為 u 只依賴 x，∂²u/∂y² = ∂²u/∂z² = 0，C₄₄ 項也消失。只剩 C₁₁ 那一項：

−ω²ρ u = C₁₁ (−K²) u ⟹ ω²ρ = C₁₁K²（式 59）

於是相速度 v_s = νλ = ω/K = (C₁₁/ρ)^(1/2)（式 60）。書上在 (58) 之後補了 K = 2π/λ 與 ω = 2πν 兩個定義，所以 νλ = (ω/2π)(2π/K) = ω/K 這個等式是自明的。

橫波：取波向量仍沿 x 立方稜，但粒子位移 v 在 y 方向：v = v₀ exp[i(Kx − ωt)]（式 61），u = w = 0。這次要代進 (57b)。C₁₁ ∂²v/∂y² = 0（v 不依賴 y）；耦合項因 u = w = 0 而消失；C₄₄(∂²v/∂x² + ∂²v/∂z²) 只剩 ∂²v/∂x² 那一半。得到 ω²ρ = C₄₄K²（式 62），故 v_s = (C₄₄/ρ)^(1/2)（式 63）。

第三支：位移改成沿 z（即 w 非零，用 (57c)），同樣的消項過程給出完全相同的速度 (C₄₄/ρ)^(1/2)。書上明確指出：K 平行於 [100] 時兩支獨立的剪切波速度相等，但在晶體中一般方向上這件事不成立。

這一句「一般方向不成立」是整段最該記住的話。兩支橫波在 [100] 簡併，不是因為橫波天生簡併，而是因為 [100] 軸上有四重旋轉對稱把 y 與 z 換過去——對稱性強迫它們相等。[110] 軸只有二重對稱，於是簡併立刻解除。

[110] 方向為什麼是重頭戲

書上開宗明義說明了理由：在立方晶的面對角線方向傳播的波之所以特別值得注意，是因為三個彈性常數可以單單從這個方向的三個傳播速度簡單地定出來。

把這句話的實驗意義說白一點。你手上有一塊立方晶，想量它的 C₁₁、C₁₂、C₄₄。

• 只做靜水壓：只能得到 ⅓(C₁₁ + 2C₁₂)，一個方程三個未知數。
• 只沿 [100] 打超音波：得到 C₁₁（縱波）與 C₄₄（橫波），兩個方程，C₁₂ 仍然拿不到，因為 [100] 的三支波裡根本沒有 C₁₂ 出現。
• 沿 [110] 打超音波：三支波分別給 ½(C₁₁ + C₁₂ + 2C₄₄)、C₄₄、½(C₁₁ − C₁₂)。三個獨立的方程、三個未知數，一次解完。

所以 [110] 是最方便、也最標準的那個「切一塊樣品、貼一組換能器就能量完全部彈性常數」的方向——這裡的「最方便」有明確的意思：三支波的偏振恰好是純縱與純橫，反解的代數可以手算。（書上只說 [110] 有 special interest、三個彈性常數可以由這個方向的三個速度 simply 求出，並沒有說它是唯一可行的方向；原則上一般低對稱方向同樣有三支波、三個獨立速度，也能反解，只是偏振不再純粹、代數不漂亮。）這也是為什麼超音波量測彈性常數的標準樣品多半取 (110) 切面（此工程慣例為我補的，Kittel 這幾頁沒有寫）。

數學上為什麼 [110] 這麼好用？因為它的波向量有兩個相等的非零分量 K_x = K_y = K/√2，於是 (57a)(57b) 的耦合項活了起來（K_xK_y ≠ 0），把 C₁₂ 從隱藏狀態拉出來；同時 z 分量為零，讓 w 那一支完全解耦，問題降成一個 2×2 行列式而不是 3×3。這是「夠複雜到能看見 C₁₂、又夠簡單到能手算」的甜蜜點。

實務：怎麼從 [110] 的三個速度反解三個彈性常數我的判讀

書上說 [110] 的價值是三個常數可以從三個速度簡單定出，但沒有把反解寫出來。既然這是這個方向存在的理由，把它補完。

量到三個速度 v_L、v_T1、v_T2（同一塊 (110) 切面的樣品，縱波換能器一次、兩個正交的橫波換能器各一次），並獨立量出密度 ρ，則：

最後兩列相加除以 2 得 C₁₁ = ρ(v_L² + v_T2²) − ρv_T1²，相減得 C₁₂ = ρ(v_L² − v_T2²) − ρv_T1²。三個常數全部拿到。

兩個實務上的提醒。第一，這三支波的偏振方向是固定的，換能器貼的方向錯了就量到別的模式，或者同時激發兩支而看到兩個回波。第二，兩支橫波的速度在 [110] 上不簡併，這反而是好事——它們在示波器上分得開，所以一次量測能得到兩個獨立資訊。如果拿 [100] 切面來量，兩支橫波簡併成同一個回波，你永遠分不出第二個獨立的剪切資訊，也就拿不到 C₁₂。

順帶一提，這也解釋了為什麼超音波量彈性常數需要密度：(C/ρ)^(1/2) 裡的 ρ 要另外量（阿基米德法或由晶格常數與分子量算）。Table 11 特地附上密度那一欄，就是為了讓讀者能自己做這個換算。

Fig. 20：三個主方向的有效彈性常數一覽

Fig. 20 是一張把結果收攏的圖，畫了三組小圖，每組一個傳播方向，箭頭標出 K 與三支波的偏振（L 代表縱波，T 代表橫波）。圖說指出這是立方晶三個主要傳播方向上三個彈性波模式的有效彈性常數，並明白寫出：在 [100] 與 [111] 方向傳播時，兩支橫波模式是簡併的。

把圖上的標示整理成表（下表的內容全部來自 Fig. 20 的標示）：

速度一律是 v = (有效常數/ρ)^(1/2)。

[111] 那一列書上正文沒有推導，而是留給習題 9 與 10（下面會做）。至於簡併與否的規律，可以用旋轉軸的階數一眼看出：[100] 是四重軸、[111] 是三重軸，階數 ≥ 3 的旋轉軸會強迫垂直於它的兩個偏振簡併；[110] 只有二重軸，沒有這個保護，於是分裂。這條規則（旋轉軸階數 ≥ 3 ⟹ 橫波簡併）是我補的群論說明，Kittel 這幾頁沒有寫，但它讓表格的三列不必死背。

一般方向：三個簡正模，但偏振不再是純橫或純縱

書上 p.85 在進入 SUMMARY 之前有一段總結性的陳述，值得逐句看清楚。

第一句：給定波向量 K 的大小與方向，晶體中有三個簡正模（normal modes）。三這個數字不是巧合——它等於每個原子的自由度數（x、y、z 三個位移分量），也等於 (57a)(57b)(57c) 這三條耦合方程。在第 4 章你會看到，若原胞裡有 p 個原子，簡正模會變成 3p 支（3 支聲學、3p − 3 支光學），而這一段處理的連續介質對應的正是那 3 支聲學支在 K → 0 的極限。

第二句也是最重要的一句：一般而言，這些模式的偏振（也就是粒子位移的方向）並不恰好平行或垂直於 K。這推翻了很多人從高中物理帶來的直覺——「波不是分成縱波與橫波兩種嗎？」在各向異性介質裡，嚴格的縱／橫分類只在特殊方向成立。一般方向上三支波是「準縱波」（quasi-longitudinal，偏振偏離 K 幾度到幾十度）與兩支「準橫波」。

第三句給出例外：在立方晶的 [100]、[111]、[110] 這三個特殊傳播方向上，三個模式中有兩個的粒子運動恰好橫向於 K，第三個恰好縱向（平行於 K）。書上並直接說，在這些特殊方向上分析比一般方向簡單得多。

為什麼恰好是這三個方向？因為它們是立方晶的三種高對稱軸（四重、三重、二重）。當 K 沿著一條旋轉軸時，位移沿該軸（縱向）與垂直該軸（橫向）分別屬於不同的對稱表示，動力學矩陣在這兩個子空間之間的矩陣元被對稱性強制為零，所以不混。K 一旦偏離這些軸，混合就恢復。這個群論解釋是我補的，但它正是書上「analysis is much simpler in these special directions」這句話背後的機制。

這件事在實驗上有後果：用超音波量彈性常數時，樣品切面的取向誤差會讓你激發到準縱波而非純縱波，量到的速度不再對應表上那個乾淨的有效常數。這是超音波彈性常數量測中最主要的系統誤差來源之一。

Table 11 與 Table 12：讀表的三個角度補推導

（先分清來源：緊接的第一小段純粹是書上 p.83 末尾正文與 Table 11 表註的轉述，不是我的推論；下面三個「角度」裡的表格、百分比與 Zener 比則都是我代入書上兩張表自己算出來的，這也是整塊標成「補推導」的原因。）

書上 p.83 末尾說明：立方晶在低溫與室溫的絕熱彈性剛度常數選錄在 Table 11，並提醒讀者注意彈性常數隨溫度升高而下降的一般趨勢；室溫的更多數值則在 Table 12。兩張表都有完整文字層，單位是 10¹² dyne/cm²（= 10¹¹ N/m²）。Table 11 的表註說 0 K 的數值是由降到 4 K 的量測外插得到的，並致謝 Charles S. Smith 教授協助編表。

角度一：溫度趨勢。拿 Table 11 的 Cu 來看，0 K 是 C₁₁ = 1.762、C₁₂ = 1.249、C₄₄ = 0.818，300 K 是 1.684、1.214、0.754，密度 9.018 g/cm³。三個常數全部下降，C₄₄ 降得最多（−7.8%），C₁₂ 降得最少（−2.8%）。這個趨勢在表中幾乎每一列都成立（Pd 是 Table 11 裡唯一的例外列：C₄₄ 由 0.712 升到 0.717，而 C₁₂ 在 0 K 與 300 K 都印 1.761、完全沒有下降）。物理原因是非簡諧：溫度升高、原子振幅變大，晶體平均而言坐到勢能曲線較平坦的位置。這也意味著聲速隨溫度下降，第 5 章的 Debye 溫度嚴格說來也是溫度的函數。

角度二：C₄₄ 與 ½(C₁₁ − C₁₂) 的比值（Zener 各向異性比 A = 2C₄₄/(C₁₁ − C₁₂)，此定義為我所補）。以下由書上兩張表的數字算出：

W 的 A = 0.995 幾乎精確等於 1，代表鎢在彈性上幾乎完全各向同性——所有方向的縱波都是 5.25 km/s、所有橫波都是 2.91 km/s（用 ρ = 19.317 g/cm³ 算）。這是 Fig. 20 那張表的一個極限情形：當 A = 1，三列的六個有效常數塌縮成兩個。另一個極端是 Na，A = 7.6，[110] 的兩支橫波速度差了快 3 倍。NaCl 的 A < 1，是少見的反向各向異性。

角度三：與週期表位置的關聯。鹼金屬（Na、Li、K）的彈性常數比過渡金屬小一到兩個數量級，而且 C₁₁ 與 C₁₂ 幾乎相等——這正是金屬鍵（第 3 章前面講過，結合來自價電子動能的降低）軟而無方向性的表現。共價晶體（Diamond、Si、Ge）則有大的 C₄₄，因為共價鍵有方向性，抵抗鍵角改變。離子晶體（NaCl、MgO）介於中間。彈性常數的模式，某種程度上就是結合類型的指紋——這是第 3 章前半（van der Waals、離子鍵、金屬鍵、共價鍵）與後半（彈性常數）之間的暗線。

SUMMARY 的五條，逐條對回前文

第 3 章的摘要只有五條，每一條都對應章內的一個小節。這裡逐條說明它在講什麼、對應哪一段推導。

第一條：惰性氣體原子的晶體由 van der Waals 交互作用（誘發偶極—偶極交互作用）束縛，其強度隨距離以 1/R⁶ 變化。這對應 p.53–59 的討論。1/R⁶ 的來源是兩個瞬時偶極之間的二階微擾能，不是靜電偶極—偶極（那會是 1/R³ 且平均為零）。這條也是 Lennard-Jones 位能 (10) 裡吸引項的來源，而 (12) 的晶格和 Σ′p⁻⁶ = 14.45392（fcc）就是把它對全晶體加總的結果。

第二條：原子間的排斥交互作用一般來自重疊電荷分布的靜電排斥，以及 Pauli 原理——它迫使重疊的平行自旋電子進入能量較高的軌域。這條有兩個機制而不是一個，Kittel 特地並列。純靜電排斥不足以解釋短程的陡峭排斥；Pauli 不相容原理才是主因。這條與第五條（共價鍵）互為表裡：Pauli 排斥對反平行自旋是減弱的，所以反平行自旋才允許大幅重疊。

第三條：離子晶體由異號離子的靜電吸引束縛。電荷為 ±q 的 2N 個離子構成的結構，其靜電能（CGS）為 U = −Nαq²/R = −NΣ(±)q²/r_ij，其中 α 是 Madelung 常數，R 是最近鄰距離。這對應 p.60–66。注意這條裡出現的是 2N 個離子、N 出現在係數上——與習題 5 的一維鏈用的是同一套記號約定（2N 個離子 = N 個「分子」）。

第四條：金屬的結合來自價電子動能相對於自由原子時的降低。這是四種結合裡唯一一個動能主導的機制。物理圖像是：電子從束縛在單一原子上（波函數侷限、動量不確定度大、動能高）變成在整塊金屬裡離域（波函數延展、動能低）。這條會在第 6 章自由電子費米氣裡被量化。

第五條：共價鍵的特徵是反平行電子自旋的電荷分布重疊。Pauli 對排斥的貢獻對反平行自旋是減弱的，這使得更大程度的重疊成為可能；重疊的電子以靜電吸引把相關的離子實束縛在一起。這條把第二條反過來用：既然 Pauli 排斥只針對平行自旋，那麼把兩個電子配成自旋相反，就能讓它們的電荷雲大幅重疊而不付出 Pauli 代價；重疊區的負電荷同時吸引兩邊的離子實，鍵就形成了。

摘要裡沒有提彈性常數。這是這一章結構上一個值得注意的地方：第 3 章的後半（p.73–85，應力應變、彈性常數、彈性波）在 SUMMARY 五條裡完全沒有出現。原因大概是這部分在體例上比較像是往第 4 章的過渡，而不是「結合」這個主題的一部分。但從習題的分佈（14 題裡有 7 題是彈性）看得出來，Kittel 並不認為它次要。

習題 1–14 導讀：哪些是推導鏈的直接延伸我的判讀

第 3 章有 14 題，橫跨整章。依主題可以分成三組：

惰性氣體與量子固體（1、2、3）：延伸 Lennard-Jones 與晶格和。第 1 題是一維盒中粒子的零點能估計，第 2、3 題用 (12) 的晶格和技術。

離子晶體（4、5、6、7）：延伸 Madelung 能。第 5 題是全章最值得手算的一題，因為一維鏈的 Madelung 常數可以用 ln 2 的級數精確求出，不需要 Ewald。

彈性（8、9、10、11、12、13、14）：七題，全部是這一段的直接延伸。其中 9、10 補上 Fig. 20 的 [111] 那一列，11 定義有效剪切常數，12、13 是行列式方法的一般化，14 是穩定性判據。8 則把 (51) 的矩陣反演接到工程上最常用的 Young's modulus 與 Poisson's ratio。

下面把其中最值得解到底的幾題做完。這些不是「額外練習」，而是正文推導鏈上被省略的環節——Kittel 把它們放到習題，是因為技巧與正文相同、只是需要動手。

這一段往後怎麼用：接到第 4、5 章與熱物性我的判讀

最後把這一段放回整本書的脈絡，說明它接下來會在哪裡被用到。

接第 4 章（聲子 I：晶格振動）。第 4 章會把原子的離散性放回去，得到單原子鏈的色散關係 ω(K)。那條曲線在 K → 0 的斜率必須等於這一段算出的聲速 (C/ρ)^(1/2)——這是連續介質極限與離散晶格模型之間的一致性檢查，也是驗算第 4 章結果的標準手法。反過來說，這一段的三支波（一縱兩橫）就是第 4 章三支聲學支在長波極限的樣子。

接第 5 章（聲子 II：熱性質）。第 5 章的 Debye 模型把每個偏振的聲速取成常數，而那個常數就是這裡算出來的 v。具體地說，Debye 溫度 θ_D 與聲速的關係是 θ_D ∝ v·n^(1/3)（n 是原子數密度），所以彈性常數大、密度小的材料（鑽石）有高 Debye 溫度，低溫熱容小、熱導率高。Table 11 顯示彈性常數隨溫度上升而下降，這也意味著 θ_D 嚴格說並非常數。

把 Si 的數字串起來看：C₁₁ = 1.66 × 10¹² dyn/cm² ⟹ 縱波 8.4 km/s、橫波 5.8 km/s ⟹ 平均聲速 ~6.5 km/s ⟹ Debye 溫度約 640–650 K。這條「彈性常數 → 聲速 → Debye 溫度 → 熱容」的鏈條，起點就在這一段。

接熱傳導與熱電容研究（與 FinFET Cth 萃取的接點）。晶格熱導率的動力學理論形式是 κ = ⅓ C v ℓ，其中 v 正是這裡的聲速、C 是晶格熱容、ℓ 是聲子平均自由程。在奈米尺度元件裡（FinFET 的 fin 寬度已經進到十奈米量級，遠小於室溫 Si 的聲子平均自由程約 100–300 nm），邊界散射讓 ℓ 被幾何尺寸截斷，κ 大幅下降，而 v 與 C 基本上不變——因為它們是材料的體性質，由這一段的彈性常數與第 5 章的熱容決定。

這對熱容 Cth 的萃取有一個直接的含意：元件層級的熱容 Cth ≈ ρ·c_p·V 幾乎不受奈米尺度效應影響，受影響的主要是熱阻 Rth。也就是說，在做 Cth 萃取時，「材料熱容是可信的體值」這個假設有相當紮實的物理基礎，而熱阻那一側才是需要小心尺寸效應的地方。這一段提供的彈性常數與聲速，正是判斷這條界線落在哪裡的第一手資料。

一個容易誤用的地方。Fig. 20 的表只適用於立方晶、且座標軸沿立方稜。如果你在做應變矽（strained Si）或 SiGe，晶體受到雙軸應變之後立方對稱被破壞（變成四方對稱），獨立彈性常數從 3 個增加到 6 個，(57) 這組方程與 Fig. 20 的表都不再適用。這在 FinFET 的模擬裡是真實存在的情境——通道區的 Si 幾乎總是處於應變狀態。