Kittel·固態物理

全書逐章 / 第 4 章 聲子 I:晶格振動 / 小節 2

02第一布里淵區/群速/長波極限/由實驗導出力常數

First Brillouin Zone / Group Velocity / Long Wavelength Limit / Derivation of Force Constants from Experiment

書本 p.93-95 · 小節 2 / 5 · 概念 12 推導 10 公式 11

上一段解出了色散關係 ω(K),但 K 是連續變數、原則上可以無限大——這一段回答「哪些 K 其實是同一件事」:原子面只在離散點上取樣波,使得相位 Ka 只有寬度 2π 的一段是獨立的,於是所有彈性波都塞得進第一布里淵區;接著用區邊界的駐波與 Bragg 對應、群速歸零、長波極限回到第 3 章的連續體彈性理論這三件事把這個框架釘死,最後把整條實驗量到的 ω(K) 做傅立葉反轉,直接讀出第 p 近鄰面的力常數——順帶發現金屬裡的有效作用力可以跨到 20 個原子面。

本段的問題:ω(K) 裡的 K 到底可以取多大?(書 p.93)

上一段已經把單原子基底的晶格振動解完了:p.92 的 (7) 給出 ω² = (2C/M)(1 − cos Ka),(9) 把它整理成 ω = (4C/M)^{1/2} |sin ½Ka|。這是一個定義在整條實數軸上的函數——把 K 代成 10π/a、100π/a,它照樣吐得出一個頻率。那麼,是不是每一個 K 都對應一個真實存在、彼此不同的晶格波?

Kittel 用一問一答開場,答案給得極短:「What range of K is physically significant for elastic waves? Only those in the first Brillouin zone.」(哪一段 K 對彈性波才有物理意義?只有第一布里淵區裡的那些。)整個 First Brillouin Zone 小節接下來做的事,就是把這句斷言證明給你看。

為什麼會有這個問題? 根源不在解出來的 ω(K),而在我們一開始設問題的方式。p.91 那裡就已經決定了:整面原子同相運動,所以「We can describe with a single coordinate u_s the displacement of the plane s from its equilibrium position」——一個面只用一個座標 u_s,s 是整數。換句話說,這個系統的自由度是可數的、離散的:每個面一個數。但我們拿來描述它的工具 exp(iKa) 卻是一個連續參數的函數。自由度比參數少,參數就一定有冗餘。第一布里淵區正是把這個冗餘除掉之後剩下的東西。

這件事的重要性遠超出第 4 章。「所有物理上不同的波都可以塞進第一布里淵區」這個結論,在第 5 章會變成「允許的 K 值填滿第一區、模態數等於自由度數」的計數基礎;在第 7 章會變成能帶的約化區表示;在第 9 章畫費米面時,整張圖就是畫在第一區裡面。這一段只有三頁,但它是後面十幾章共用的座標系。

相鄰兩原子面的位移比

u_{s+1} / u_s = u exp[i(s + 1)Ka] / u exp(isKa) = exp(iKa)

(10)  p.93

u_s 為第 s 個原子面離開平衡位置的位移;u 為振幅(與 s 無關的常數);K 為波向量;a 為相鄰原子面的間距,其數值依 K 的方向而定(p.92 在 (4) 底下寫明「The value to use for a will depend on the direction of K」);s 為標記原子面的整數。

出處推導
直接把 p.92 的 (4) 拿來作比。(4) 印的是 u_{s±1} = u exp(isKa) exp(±iKa),也就是行進波解 u_s = u exp(isKa) 的相鄰面寫法(u_s = u exp(isKa) 這個寫法書中在 p.93 才明文出現);Kittel 自己在 (10) 之前就寫「From (4) the ratio of the displacements of two successive planes is given by」,所以「由 (4) 作比」正是原書的做法。注意 (4) 的完整解還帶一個時間因子 exp(−iωt)(p.91 說「We look for solutions with all displacements having the time dependence exp(−iωt)」),但這裡比較的是同一時刻的兩個面,時間因子上下相消,所以不影響結論。
極限行為

這一式是整段的樞紐,因為它說明了:兩個相鄰面之間唯一的物理內容就是相位差 Ka,振幅比恆等於 1(每個面振幅都一樣,這是行進波不衰減的表現)。既然 exp(iθ) 對 θ 是 2π 週期的,Ka 與 Ka + 2πn 描述的是一模一樣的面間關係。

三個特例值得記住:Ka = 0 時比值為 +1,所有面同相位移,那是整塊晶體的剛體平移,對應 (7) 給的 ω = 0(平移不需要能量,是聲學支在 K → 0 必然歸零的原因);Ka = π 時比值為 −1,相鄰面恰好反相,這就是 (12) 的駐波;Ka 介於兩者之間時比值是單位圓上的複數,相位隨面號線性累積,那才是一般的行進波。

根因:原子面是波的離散取樣器(補充說明)我的判讀

Kittel 把理由壓縮成一句「The range −π to +π for the phase Ka covers all independent values of the exponential.」,讀起來像是純數學的陳述。把它翻成物理語言會清楚很多,而且這個翻譯值得記一輩子。

在連續彈性介質裡,位移場 u(x) 對每一個 x 都有定義,所以兩個不同波數的波 exp(iK₁x) 與 exp(iK₂x) 只要 K₁ ≠ K₂ 就是兩個不同的函數,永遠可以被區分。晶體裡不是這樣:位移只在 x = sa 這些離散點上有意義,s 是整數。把 exp(iKx) 在這些點上取樣,得到的是數列 exp(isKa)。現在把 K 換成 K + 2πn/a:

exp[i(K + 2πn/a)·sa] = exp(isKa)·exp(i2πns) = exp(isKa)

因為 s 與 n 都是整數,exp(i2πns) 嚴格等於 1。這不是近似、不是「差不多相等」,是恆等式。兩個相差 2πn/a 的波向量,在所有可觀測的點上給出完全相同的數值。

訊號處理的讀者會立刻認出這是取樣定理(Nyquist)在晶體裡的版本:取樣間距 a,可分辨的最大波數是 π/a,對應的最短波長是 2a;比這更短的波會「折疊」(alias)回低波數。這個類比很好用,但用完要回到精確陳述:在訊號處理裡,混疊是取樣造成的資訊損失,是我們選擇取樣率的後果;在晶體裡,離散性是物理本身的性質——原子之間真的沒有東西可以位移,所以那些「被折疊掉」的自由度從一開始就不存在,沒有任何資訊被損失。這是為什麼 Kittel 說區外的 K「merely reproduce lattice motions described by values within the limits ±π/a」(只不過是重複了 ±π/a 以內的值所描述的晶格運動),用的是 reproduce 而不是 approximate。

(取樣定理這個類比是本文補充的教學說明,Kittel 在這幾頁並未使用這個詞。)

Figure 5 的文字重建(書 p.93)

網頁不放圖,先把 p.93 下方那張圖描述出來,因為它是整段最有效的一句話。

圖是一條灰底的橫幅。上面畫了兩條曲線。一條是白色實線,短波長的高頻正弦波,在整個圖寬裡走了約九個完整週期(可以數到九個波峰),看起來像一排密集的鋸齒。另一條是黑色虛線,長波長的正弦波,整個圖寬只走了一個多(約 1.3 個)週期,是一條和緩的起伏。兩條曲線的振幅相同,但頻率差了約六倍——直覺上它們八竿子打不著。

關鍵在第三個元素:圖上有一排黑色實心圓點,標出原子面的位置。這些點同時落在白色實線與黑色虛線上。也就是說,兩條曲線在原子位置處的數值完全相同,只有在原子與原子之間才分道揚鑣。圖的左下角有一個雙箭頭線段,兩端各有一條短豎線,標注字母 a,指出相鄰兩個原子點的水平間距。

圖說逐字是:「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.」(實線所代表的波,沒有傳達任何虛線未給出的資訊。只需要波長長於 2a 的波就足以表示這個運動。)

看這張圖有一個必須先接受的前提:曲線在原子之間的值沒有物理意義。白色實線畫得再細緻,我們也永遠量不到它在兩個原子之間的值——那裡沒有東西可以位移,沒有任何實驗能測到那個數字。兩條曲線在唯一可觀測的地方(黑點)完全重合,所以它們描述的是同一個晶格運動。實線的 K 落在第一區之外,虛線的 K 落在第一區之內,圖說的第二句就是結論:只需要 λ > 2a 的波。

反過來說,這張圖也解釋了為什麼區邊界是 λ = 2a 而不是別的值:虛線再長一點沒關係(區內),但如果虛線縮短到比 2a 還短,它自己就會變成「實線」,被另一條更長的波取代。λ = 2a 是「不會被別人取代」的臨界波長。

獨立波向量的範圍(第一布里淵區)

−π < Ka ≤ π , 或 −π/a < K ≤ π/a

p.93

K 為波向量,a 為面間距。端點記為 K_max = ±π/a(原文:「The extreme values are K_max = ±π/a.」)。

出處推導
由 (10):兩面位移比只透過 exp(iKa) 依賴 K,而複指數對相位是 2π 週期的,因此任何寬度為 2π 的相位區間都涵蓋了全部獨立值。Kittel 取的是對稱於原點的那一段,原文寫「The range −π to +π for the phase Ka covers all independent values of the exponential.」接著一句「This range is the first Brillouin zone of the linear lattice, as defined in Chapter 2.」把定義的責任丟回第 2 章。此式在原書沒有編號。
極限行為

兩端各有一個極限值得分開想。K → 0(λ → ∞):波長遠大於 a,波「看不見」原子,系統退化成連續彈性體,這就是 Long Wavelength Limit 小節與第 3 章的接點。K → ±π/a(λ → 2a):這是這個晶格能承載的最短波長,行進波退化成駐波 (12),群速歸零 (14),對應 Bragg 反射。整個第一區就是這兩個極限之間的全部物理。

另外注意這是一個半開區間:左端 −π 被排除(嚴格小於),右端 +π 被保留(小於等於)。這不是排版疏忽,理由見下一個推導。

補推導:為什麼區間是半開的(−π < Ka ≤ π 而不是 −π ≤ Ka ≤ π)補推導

這一步書上完全沒寫,但不補的話會在兩個地方出事。

第一個地方是計數。第 5 章要數一維鏈上有多少個獨立模態,答案必須剛好等於可以動的原子數,多一個少一個都代表模型錯了。如果把 ±π/a 都算進去,一維的計數立刻會多出一個,整套「模態數 = 自由度數」的論證就垮了。書 p.111 那裡在做的「one allowed value of K per volume in K space」這類計數,隱含的就是這個半開約定。

第二個地方是畫圖。ω(K) 對 K 是偶函數((7) 裡 K 只透過 cos Ka 出現),所以 Fig. 4 那張圖左右對稱。很多人因此以為「負 K 和正 K 是兩組不同的波,只是頻率相同」。在區內部(|K| < π/a)這個理解是對的:+K 與 −K 是往右跑與往左跑的兩個不同的行進波。但在端點不對:±π/a 這兩個「方向相反的波」其實是同一個不跑的駐波。「往右」與「往左」在這裡失去意義,這正是 (12) 要講的事。

  1. 假設兩端都算進來,那麼 K = −π/a 與 K = +π/a 會被當成兩個不同的波向量。
  2. 但這兩者相差 (+π/a) − (−π/a) = 2π/a,恰好是一維晶格的一個倒晶格向量(第 2 章 Fig. 11,書 p.35:倒晶格基向量 b 的長度就是 2π/a)。按照 (11) 的邏輯,它們必須是同一個態。
  3. 直接驗算。代入 (10):K = +π/a 給 u_{s+1}/u_s = exp(iπ) = −1;K = −π/a 給 u_{s+1}/u_s = exp(−iπ) = −1。兩者完全相同。
  4. 代入 (12):K = +π/a 給 u_s = u exp(+isπ) = u(−1)^s;K = −π/a 給 u_s = u exp(−isπ) = u(−1)^s。又是完全相同。這正是 Kittel 在 (12) 把兩個符號寫在一起成 u exp(±isπ) 的原因——他不是在列兩個解,是在說這兩個寫法給同一個答案。
  5. 頻率也一樣:由 (9),ω(±π/a) = (4C/M)^{1/2}|sin(±π/2)| = (4C/M)^{1/2}。同一個頻率。
  6. 結論:±π/a 是同一個物理態的兩個名字,只能算一次。半開區間 −π < Ka ≤ π 就是「算一次」最簡單的約定(保留右端、丟掉左端;反過來約定也完全可以,只是 Kittel 選了這一種)。

回扣第 2 章:布里淵區到底是什麼(書 p.34–35、p.37)

Kittel 在這裡只寫「as defined in Chapter 2」就過去了。把第 2 章的定義搬回來對照是值得的,因為兩章的動機完全不同,而它們在 p.94 才被接起來。

第 2 章 p.34 的正式定義是:「The first Brillouin zone is the smallest volume entirely enclosed by planes that are the perpendicular bisectors of the reciprocal lattice vectors drawn from the origin.」(第一布里淵區是被那些「由原點畫出的倒晶格向量之垂直平分面」完全圍住的最小體積。)同一頁還給了另一個等價說法:中央那一格「is a Wigner-Seitz cell of the reciprocal lattice」;p.37 再講一次「Another primitive cell is the central (Wigner-Seitz) cell of the reciprocal lattice which is the first Brillouin zone.」

一維的情形第 2 章 Fig. 11(p.35)已經畫好了。圖上是一條晶格常數 a 的線性晶格,下面畫出它的倒晶格:基向量 b 的長度等於 2π/a,由原點出發的最短倒晶格向量是 b 與 −b;它們的垂直平分點落在 k = ±π/a。圖說最後一句逐字是「The boundaries are at k = ±π/a.」——和第 4 章這裡算出來的 K_max = ±π/a 是同一件事。

兩章的動機差在哪裡? 第 2 章是從 X 光繞射出發的:p.34 說這些垂直平分面之所以重要,是因為「A wave whose wavevector drawn from the origin terminates on any of these planes will satisfy the condition for diffraction.」——終點落在這些面上的波向量滿足繞射條件。第 4 章是從晶格振動出發的:區邊界是「相鄰原子面反相、行進波變駐波」的 K。這兩件事聽起來毫無關係,但 Kittel 在 p.94 把它們縫在一起:區邊界的駐波「is equivalent to Bragg reflection of x-rays」。同一個幾何條件,一次作用在 X 光上,一次作用在聲子上,第 7 章還會再作用在電子上一次。

還有一句容易滑過去但很關鍵的話,在 p.93 (11) 的解釋裡:「We note that 2πn/a is a reciprocal lattice vector because 2π/a is a reciprocal lattice vector.」折回第一區時減掉的那個量,不是隨便挑的數字,它必須是倒晶格向量。這是為什麼「第一布里淵區」這個概念在三維、在非立方晶系裡依然成立:只要把 2πn/a 換成一般的倒晶格向量 G,整套論證一字不改。

把區外的 K 折回第一區

u_{s+1}/u_s = exp(iKa) ≡ exp(i2πn) exp[i(Ka − 2πn)] ≡ exp(iK′a) , 其中 K′ = K − 2πn/a

(11)  p.93

K 為落在第一區之外的原始波向量;n 為整數;K′ ≡ K − 2πn/a 為減去 n 個倒晶格向量之後落回第一區的等價波向量。原文:「Suppose K lies outside the first zone, but a related wavevector K′ defined K′ = K − 2πn/a lies within the first zone, where n is an integer.」

出處推導
技巧只有一個:把 exp(iKa) 乘上一個恆等於 1 的因子再拆開。中間那一步 exp(i2πn)·exp[i(Ka − 2πn)] 把原本的相位拆成「整數圈」加「餘下的相位」,而整數圈 exp(i2πn) = 1(原文:「because exp(i2πn) = 1」),剩下的 Ka − 2πn = (K − 2πn/a)a = K′a。
極限行為

n = 0 時 K′ = K,本來就在區內,這一式是恆等的空操作;|n| 愈大代表原本的 K 落在愈外面的區(n = ±1 對應第二區,以此類推)。這個操作是可逆且唯一的:給定任何 K,恰好只有一個整數 n 能把它送進半開區間 −π/a < K′ ≤ π/a。因此第一區裡的 K 不多不少,剛好把所有物理上不同的晶格波參數化一次。

另外請留意 Kittel 在這一式用的是恆等號 ≡ 而不是等號 =。(11) 的這兩步對任意整數 s、任意整數 n 都成立,本來就是恆等式,而不是某個條件下的相等;書上這裡剛好印了 ≡。不過 Kittel 的 ≡/≅ 用法並不嚴格——同一頁那個純粹是二階截斷的 cos Ka 展開也印成 ≡(見後面「長波展開所用的餘弦近似」那一式的【符號校註】)——所以不宜由符號反推作者的用意。

補推導:(11) 只證了一半——還得確認頻率也相同補推導

Kittel 省略這一步,因為對他來說「(7) 裡 K 只出現在 cos Ka 裡」一望即知。但初學者常常只記住「位移一樣」,然後在畫擴展區(extended zone)與約化區(reduced zone)圖時搞混:約化區表示法之所以合法,靠的就是 ω(K) 的週期性,不只是位移比的週期性。

這個折回操作在後面會被系統化。第 7 章 p.177 起正式引入 reduced zone scheme(約化區表示法),第 9 章 p.223–228 在畫費米面時反覆使用 reduced / periodic / extended 三種表示法互相換算。那裡做的事情,在數學上與這裡的 (11) 一模一樣,只是把 2πn/a 換成一般的倒晶格向量 G,把 ω(K) 換成 ε(k)。這一段的三行代數,是後面那幾十頁圖的全部理論基礎。

  1. (11) 證明的是「位移型態相同」。但一個模態不是只有位移型態,它還有頻率。如果 K 與 K′ 給同一個位移型態卻給不同的頻率,那就自相矛盾了——同一個運動不可能有兩個振動頻率。所以嚴格說,還欠一步。
  2. 由 p.92 的 (7):ω²(K) = (2C/M)(1 − cos Ka)。K 只透過 cos Ka 進入這個式子。
  3. 代入 K = K′ + 2πn/a:cos(Ka) = cos[(K′ + 2πn/a)·a] = cos(K′a + 2πn) = cos(K′a),因為餘弦函數的週期是 2π 而 n 是整數。
  4. 所以 ω²(K) = ω²(K′),也就是 ω(K) 本身就是以 2π/a 為週期的函數:ω(K + 2πn/a) = ω(K) 對任意整數 n 成立。
  5. 群速跟著週期:v_g = dω/dK 是週期函數的導函數,所以 v_g(K + 2πn/a) = v_g(K)。用 (14) 直接驗也行:cos[½(K + 2π/a)a] = cos(½Ka + π) = −cos ½Ka,看似差一個負號——但別忘了 (9) 裡還有絕對值 |sin ½Ka|,跨過區邊界時 sin ½Ka 的符號翻轉會把這個負號吃掉,淨結果仍是週期的。
  6. 結論:K 與 K′ 不只是「位移看起來一樣」,而是位移、頻率、群速、能量——所有可觀測量全部相同。它們是同一個模態的兩個名字。

區邊界的駐波

u_s = u exp(± isπ) = u (−1)^s

(12)  p.93

u_s 為第 s 面的位移;u 為振幅;s 為整數面號。±(正負兩個符號)對應區邊界的兩端 K = ±π/a,兩者給出同一個結果。

出處推導
在區邊界 K_max = ±π/a,於是每一面的相位是 sK_max·a = ±sπ(原文:「At the zone boundaries sK_max a = ±sπ, whence」),代回行進波解 u_s = u exp(isKa)(此寫法見 p.93;(4) 本身印的是等價的相鄰面形式 u_{s±1} = u exp(isKa) exp(±iKa))。再用 exp(±iπ) = −1 得到 (−1)^s。
極限行為

s 為偶數時 u_s = +u,s 為奇數時 u_s = −u:相鄰原子面反相。位移的空間型態退化成一個固定的、不隨時間移動的圖樣,只有整體振幅在隨 exp(−iωt) 振盪。Kittel 的原話是「This is a standing wave: alternate atoms oscillate in opposite phases, because u_s = ±1 according to whether s is an even or an odd integer. The wave moves neither to the right nor to the left.」(此為駐波:相間的原子以相反的相位振盪……這個波既不向右也不向左移動。)

務必分清楚兩件事:波形不動,但原子在動。每一個原子都在以 ω_max = (4C/M)^{1/2} 全幅振盪,振幅是 u,一點也沒有減少。停下來的是「圖樣的移動」,也就是能量的淨輸運——這與 (14) 在 K = π/a 給出 v_g = 0 是同一件事的兩種說法。

把 K 稍微拉離邊界一點(K = π/a − δ),u_s 就恢復成緩慢旋轉的相位,波又開始跑,但跑得很慢(v_g ≈ (Ca²/M)^{1/2}·½δa)。駐波不是一個孤立的異常點,而是群速連續趨近於零的終點。

典型數值
區邊界頻率 ω_max = (4C/M)^{1/2} = 2(C/M)^{1/2}。以銅沿 [100] 的縱波估算(詳細數字見下面「接回第 3 章」那一段):ω_max ≈ 4.9 × 10¹³ rad/s,對應 ν_max ≈ 7.8 THz,ħω_max ≈ 32 meV。這就是為什麼聲子能量在文獻裡總是以 meV 為單位、而中子(熱中子動能約 25 meV)剛好是量測聲子色散最合適的探針——兩者能量同一量級,散射時的能量轉移才量得準(第 4 章 p.100–102 的中子非彈性散射)。

區邊界 = Bragg 反射:Kittel 給的數字與一個抽取陷阱(書 p.94)

(12) 只說了「是駐波」,但沒說為什麼偏偏在 ±π/a 變成駐波。p.94 開頭補上這個為什麼,而且是用一個乍看毫不相干的領域來補的——X 光繞射。

原文逐字:「This situation is equivalent to Bragg reflection of x-rays: when the Bragg condition is satisfied a traveling wave cannot propagate in a lattice, but through successive reflections back and forth, a standing wave is set up.」(這個情形等價於 X 光的 Bragg 反射:當 Bragg 條件被滿足時,行進波無法在晶格中傳播,而是透過來回的連續反射建立起一個駐波。)

這句話的機制值得慢慢讀。一個往右跑的波,在每一個原子面上都被散射一點點。一般情況下這些散射波的相位雜亂,互相抵消,不影響主波前進。但當 Bragg 條件成立時,所有面的反向散射波同相疊加,於是往右的波被完整地轉換成往左的波,往左的波又被轉回往右——兩個等幅反向的行進波疊加起來,就是駐波。這是一個共振現象,不是吸收:能量沒有損失,只是被鎖在原地。

Kittel 接著把數字代進去。原文那一行是:「The critical value K_max = ±π/a found here satisfies the Bragg condition 2d sin θ = nλ: we have θ = ½π, d = a, K = 2π/λ, n = 1, so that λ = 2a.」

【抽取校註】 這份 PDF 的文字層會把分數拆成浮出的片段,½ 會被抽成分離的「1」與「2」並漂到段落末端,直接照抄會得到「θ = 1 2π」這種無意義的字串。我們把 p.94 這一行以 4 倍解析度算繪成影像放大目視,確認原書印的是 θ = ½π(分子 1、分母 2 的堆疊分數,後接 π),也就是 90°。本文採用目視判讀的結果。

補推導:把 Kittel 給的四個數字代進 Bragg 條件補推導

書上只給一串等號分隔的數值,沒有寫代入過程。補出來之後可以看到,「區邊界」「λ = 2a」「相鄰面反相」「Bragg 條件成立」這四句話是同一件事的四種說法,彼此可以自由換算。這種「一個物理事實、四個等價表述」的結構在固態物理裡到處都是,建議每遇到一次就把換算走一遍。

還有一個第 2 章給的限制,放在這裡看特別漂亮。p.25 在 (1) 之後緊接著寫:「This is the Bragg law, which can be satisfied only for wavelength λ ≤ 2d.」——Bragg 定律只有在 λ ≤ 2d 時才可能被滿足(因為 sin θ ≤ 1)。現在把 d = a 代入:Bragg 反射要求 λ ≤ 2a。而 Fig. 5 的結論剛好相反:獨立的晶格波要求 λ ≥ 2a。兩個條件只在唯一的一點 λ = 2a 交會——這就是區邊界。換句話說,在整個第一布里淵區裡,只有邊界那一個點同時滿足「是獨立的晶格波」與「滿足 Bragg 條件」。n ≥ 2 的高階 Bragg 反射對應 λ = 2a/n < 2a,全部落在「不獨立」的那一側,這就是下一段要講的「彈性波通常只能取 n = 1」的幾何版說法。

  1. Bragg 條件是第 2 章 p.25 的 (1):2d sin θ = nλ。
  2. d = a:發生散射的是原子面,相鄰面的間距就是我們一路在用的 a。
  3. θ = ½π,即 90°,sin θ = 1。這裡最容易出錯,而第 2 章早就把定義寫死了——p.25 逐字:「The path difference for rays reflected from adjacent planes is 2d sin θ, where θ is measured from the plane.」θ 是從晶面本身量起的,不是從面法線量起(後者是幾何光學的習慣)。所以 θ = 90° 表示入射線與晶面垂直,也就是波沿著面法線傳播——正是我們從 p.91 開始討論的情形:波向量沿 [100],整面原子同相運動。θ = ½π 不是隨手挑的角度,它是這個一維問題唯一可能的角度。
  4. n = 1(理由見下一段)。
  5. 代入:2 × a × 1 = 1 × λ ⇒ λ = 2a。
  6. 換成波向量:K = 2π/λ = 2π/(2a) = π/a = K_max。與 (12) 的區邊界完全一致。
  7. 反向檢查也成立:把 K = π/a 代回 λ = 2π/K 得 λ = 2a,代回 (10) 得相鄰面相位差 Ka = π,也就是反相——與 (12) 一致,與 Fig. 5 圖說的「Only wavelengths longer than 2a are needed」一致。三條路殊途同歸。

為什麼彈性波通常只取 n = 1,X 光卻可以有高階(書 p.94)

上面代入 Bragg 條件時直接寫了 n = 1,但 Bragg 條件本來允許 n = 1, 2, 3, …。n = 2 給 λ = a,n = 3 給 λ = 2a/3,依此類推。這些高階解到哪去了?

Kittel 的回答只有一句,但這一句是整段最有教學價值的一句:「With x-rays it is possible to have n equal to other integers besides unity 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.」(用 X 光時 n 可以取 1 以外的整數,因為電磁波的振幅在原子之間的空間裡也有意義;但彈性波的位移振幅通常只在原子本身的位置上才有意義。)

把這句話和 Fig. 5 對起來看就通了。λ = 2a/n(n ≥ 2)的波,在連續介質裡是合法的、與 λ = 2a 的波完全不同的波。但在晶格裡,它在原子位置上的取樣值會與某個 λ ≥ 2a 的波完全重合——它就是 Fig. 5 裡那條白色實線,而 λ = 2a 那條是虛線。所以 n ≥ 2 的解不是被禁止,是根本不是新的解,它們已經包含在 n = 1 裡了。

X 光不一樣。電磁場 E(x) 是一個定義在所有 x 上的連續場;原子之間的電場是真實的、原則上可測的物理量(它會對放在那裡的任何電荷施力)。所以 λ = 2a/2 = a 的 X 光與 λ = 2a 的 X 光,在原子之間的行為不同,是兩個真正不同的波,高階繞射峰因此才存在、才被觀測到。

這是一組很乾淨的對照:同一個晶格、同一個 Bragg 條件,因為場的「定義域」不同,而得出不同的階數規則。差別不在近似的好壞,在自由度的本質——連續場 vs. 離散自由度。

最後留意 Kittel 用了 usually(通常)這個保留字,沒有把話說死。他在這幾頁沒有解釋例外的情形,本文也不替他補。

Group Velocity:為什麼要引入第二個速度(書 p.94)

到這裡為止,我們手上只有單一平面波解 u_s = u exp[i(sKa − ωt)](空間部分 u_s = u exp(isKa) 見 p.93,(4) 印的是它的相鄰面寫法 u_{s±1} = u exp(isKa) exp(±iKa);時間因子 exp(−iωt) 來自 p.91)。這個解有一個物理上很不舒服的性質:它充滿整個晶體,從 s = −∞ 延伸到 s = +∞,而且從 t = −∞ 就一直在那裡振盪。這樣的波不「傳遞」任何東西——每個原子都在做穩態的簡諧振動,沒有任何訊號從一處送到另一處。

真正能傳遞能量與訊號的是波包:把 K₀ 附近一小段的平面波疊加起來,彼此的相位在某個位置相長、在別處相消,得到一個局域化的擾動。這個擾動的包絡(envelope)會移動,它的移動速度就是群速。Kittel 直接給結論,跳過了推導:「The transmission velocity of a wave packet is the group velocity」(波包的傳遞速度就是群速)。

他先寫純量式 v_g = dω/dK(未編號),再寫向量式 (13) v_g = grad_K ω(K)(v_g 與 K 都是粗體),並補上「the gradient of the frequency with respect to K. This is the velocity of energy propagation in the medium.」(頻率對 K 的梯度。這是介質中能量傳播的速度。)

一維時兩式當然相同,那為什麼要多寫一個向量式?因為 (13) 帶出一個一維完全看不見的事實:群速的方向是等頻面 ω(K) = 常數 的法線方向,一般並不平行於 K。在各向異性的晶體裡,你沿著某個方向送進一個波包,它的能量可以歪著跑出去(這個現象叫聲子聚焦,Kittel 在這幾頁沒有討論)。這個梯度形式在第 5 章計算聲子態密度與熱導率時是主角(書 p.118 的態密度積分式裡出現的就是「the magnitude of the group velocity of a phonon」),在第 7 章之後描述電子的群速時也是同一個公式換個符號而已。

群速的定義(向量形式)

𝐯_g = grad_K ω(𝐊) ;一維純量形式 v_g = dω/dK

(13)  p.94

𝐯_g 為群速向量(原書以粗體印);𝐊 為波向量(粗體);ω(𝐊) 為色散關係;grad_K 表示對 K 的三個分量取梯度。純量式 v_g = dω/dK 印在 (13) 之前,原書未編號。

出處推導
不是從本章任何一式推出來的,是波包分析的標準結果,Kittel 在此直接引用。物理內容:把 ω(K) 在 K₀ 附近展開成 ω ≈ ω₀ + (dω/dK)|₀ (K − K₀),疊加出來的波包其包絡的移動速度就是那個一階係數。
極限行為

與相速的區別是這一式存在的全部理由。 相速 v_p = ω/K 描述的是等相位面移動得多快,那是一個記帳用的數學量;群速描述的是能量與訊號跑得多快,那才是可測的物理。只有在 ω 嚴格正比於 K(無色散)時兩者才相等——長波極限 (15) 正是這種情形。

在三維各向異性晶體中,grad_K ω 垂直於等頻面,因此 𝐯_g 一般不平行於 𝐊:能量流的方向與波的「傳播方向」可以不同。一維模型完全掩蓋了這件事,這是 Kittel 特地補寫向量形式的理由。

補推導:從 (9) 微分得到 (14)補推導

書上只寫「With the particular dispersion relation (9), the group velocity (Fig. 6) is (14)」,中間三行微分完全省略。補出來的價值不只是驗算:第二條路顯示 p.92 的 (8) 其實就是群速在區邊界歸零的預告片——Kittel 在還沒定義群速之前,就先用 (8) 把「零斜率」這個事實放在那裡,等到 p.94 才回頭告訴你那個斜率的名字叫群速、它為零的意思是能量不再傳播。讀第 4 章時把 (8)、(12)、(14) 三式放在一起看,會發現它們講的是同一句話。

  1. 出發點是 p.92 的 (9):ω = (4C/M)^{1/2} |sin ½Ka|。
  2. 先限定在 0 < K < π/a 這一段。此時 0 < ½Ka < π/2,sin ½Ka > 0,絕對值可以直接拿掉:ω = (4C/M)^{1/2} sin ½Ka。
  3. 對 K 微分,用鏈鎖律(內層 ½Ka 對 K 的導數是 ½a):dω/dK = (4C/M)^{1/2} · cos(½Ka) · (a/2)。
  4. 把 a/2 收進根號裡:(4C/M)^{1/2} · (a/2) = [(4C/M) · (a²/4)]^{1/2} = (Ca²/M)^{1/2}。這一步是 (14) 那個漂亮前置因子的來源——4 與 1/4 消掉了。
  5. 得 v_g = (Ca²/M)^{1/2} cos ½Ka,即 (14)。
  6. 量綱檢查:[C] = N/m = kg·s⁻²,[a²] = m²,[M] = kg,所以 Ca²/M 的量綱是 (kg·s⁻²·m²)/kg = m²·s⁻²,開根號得 m/s。✓
  7. 另一條路(順便解釋 p.92 的 (8)):直接對 (7) 兩邊微分。由 ω² = (2C/M)(1 − cos Ka) 得 dω²/dK = (2Ca/M) sin Ka,這就是書上的 (8)。再用 dω²/dK = 2ω·(dω/dK) 得 v_g = (Ca/M) sin Ka / ω。代入 ω = (4C/M)^{1/2} sin ½Ka 與倍角公式 sin Ka = 2 sin ½Ka cos ½Ka:v_g = (Ca/M)·2 sin ½Ka cos ½Ka / [(4C/M)^{1/2} sin ½Ka] = (Ca²/M)^{1/2} cos ½Ka。同一個答案。
  8. 這條路還順便說明了 (8) 為什麼重要:dω²/dK = (2Ca/M) sin Ka 在 K = ±π/a 時 sin(±π) = 0,所以 ω² 對 K 的斜率為零,ω 對 K 的斜率也為零——書 p.92 說「We show from (7) that the slope of ω versus K is zero at the zone boundary」,講的就是這件事,而且那裡就已經預告「The special significance of phonon wavevectors that lie on the zone boundary is developed in (12) below.」

單原子鏈的群速

v_g = (Ca²/M)^{1/2} cos ½Ka

(14)  p.94

v_g 群速;C 相鄰面之間、以一個原子為單位定義的力常數(p.91:「It is convenient hereafter to regard C as defined for one atom of the plane」);a 面間距;M 原子質量;K 波向量。

出處推導
把 p.92 的 (9) ω = (4C/M)^{1/2}|sin ½Ka| 對 K 微分(詳見上一段補推導)。前置因子 (Ca²/M)^{1/2} 由 (4C/M)^{1/2} × (a/2) 化簡而來。
極限行為

K → 0:cos → 1,v_g → v₀ ≡ (Ca²/M)^{1/2}。這是長波聲速,與 (15) 給的相速 ω/K 相同,也就是第 3 章那個「與頻率無關的聲速」。

K = π/2a(區的一半處):cos(π/4) = 0.707,群速已經掉到長波值的 71%。這個數字值得記住——色散在區的中段就已經很嚴重了,不是到邊界才突然出事。德拜近似(第 5 章)假設 ω = vK 一路線性到截止頻率,誤差的主要來源就是這裡:它把一個已經慢了 29% 的群速當成滿速。

K = π/a(區邊界):cos(π/2) = 0,v_g = 0,能量完全不傳播,與 (12) 的駐波互相印證。原文:「This is zero at the edge of the zone where K = π/a. Here the wave is a standing wave, as in (12), and we expect zero net transmission velocity for a standing wave.」

v_g = 0 的後果:第 5 章 p.119 指出「There is a special interest in the contribution to [態密度] from points at which the group velocity is zero. Such critical points produce singularities (known as Van Hove singularities) in the distribution function」。群速歸零 ⇒ 態密度 D(ω) 在該頻率發散。所以區邊界不只是一個幾何邊界,它在聲子頻譜上會留下一個可以被實驗看見的尖峰。

典型數值
以銅沿 [100] 的縱波為例,v₀ ≈ 4.4 km/s(算法見下面「接回第 3 章」)。在 K = π/2a 處 v_g ≈ 3.1 km/s,在 K = π/a 處 v_g = 0。整個從 4.4 km/s 掉到 0 的過程,發生在 K 從 0 走到 π/a ≈ 1.7 × 10¹⁰ m⁻¹ 的範圍內。

Figure 6 的文字重建(書 p.95 左上)

圖印在 p.95 的左上角,是一張單曲線的方框圖。

縱軸標的是無因次比值 v_g/(Ca²/M)^{1/2},也就是「群速除以長波聲速」,刻度標了 0、0.5、1.0 三格。橫軸是 K,刻度標了 0、π/2a、π/a 三處,軸標籤 K 寫在正下方。

曲線從左上角的 (0, 1.0) 出發。起點處曲線是水平的——斜率為零,這一點非常關鍵。往右走曲線逐漸下彎,在 K = π/2a 附近通過大約 0.71 的高度,接著愈來愈陡,最後在右邊界 K = π/a 處以整條曲線最陡的斜率交到 0——注意斜率是有限的,不是垂直:歸一化後 d(v_g/v₀)/d(K/(π/a)) = −π/2 ≈ −1.57,在原圖的長寬比下目視大約只有 45°–50°。

圖說逐字:「Group velocity v_g versus K, for model of Fig. 4. At the zone boundary K = π/a the group velocity is zero.」(群速 v_g 對 K 作圖,用的是 Fig. 4 的模型。在區邊界 K = π/a 處群速為零。)

這條曲線就是 cos ½Ka 在 0 ≤ ½Ka ≤ π/2 的四分之一週期,但「認出它是餘弦」之後還要再問兩個為什麼:

為什麼起點是水平的? 因為 d(cos x)/dx = −sin x 在 x = 0 為零。物理意義是:長波區有一整段 K,群速幾乎不隨 K 改變。這正是「聲速」這個概念能夠成立的原因——如果群速在 K → 0 時就已經敏感地隨 K 變化,「材料的聲速」就會是一個沒有意義的量,因為它會取決於你用什麼頻率去量。第 3 章的彈性常數理論,就靠這一小段水平。

為什麼終點最陡? 因為 −sin x 在 x = π/2 取到最大值 −1,那是整條曲線最陡的地方。物理意義是:在區邊界附近,群速對 K 極度敏感。一個實驗上的小誤差、一點溫度造成的晶格常數改變,都會讓區邊界附近的群速大幅變動。這也是為什麼區邊界附近的聲子輸運行為特別難算。

最後把 Fig. 4(p.92 的 ω–K 圖)與 Fig. 6 疊在腦子裡看,因為它們是同一個資訊的兩種畫法。

Fig. 4 的縱軸是 ω/(4C/M)^{1/2}(刻度到 1.2,曲線最高到 1.0),橫軸從 −π/a 一路畫到 2π/a,圖中並標示了曲線的解析式 ω/(4C/M)^{1/2} = |sin ½Ka|。圖的下方有一個雙箭頭括弧,從 −π/a 標到 π/a,注記「First Brillouin zone」。曲線在 K = 0 處為零且形成一個尖點(V 形,因為有絕對值),往兩側各以斜率 ±v₀ 近乎直線上升,然後逐漸彎平,在 K = ±π/a 達到水平的頂點 1.0,接著在 π/a 到 2π/a 之間又對稱地落回 0。

Fig. 4 特地畫到 2π/a,是 (11) 的圖像證據。 區外 π/a < K < 2π/a 那一段不是新資訊:它就是區內 −π/a < K < 0 那一段整體平移 2π/a 的結果,兩段完全重合。書上沒有明說這一點,但只要注意到圖被畫到 2π/a 而 caption 只講連續體近似,就會發現這張圖同時在做兩件事:左半講長波極限,右半講週期性。

Fig. 6 就是 Fig. 4 在 0 ≤ K ≤ π/a 這一段的斜率。 Fig. 4 的「近直線段」對應 Fig. 6 的「水平段(v_g ≈ v₀)」,Fig. 4 的「水平頂點」對應 Fig. 6 的「落到零」。一張看頻率、一張看斜率,缺一不可:只看 Fig. 4 很容易把區邊界附近的平坦誤以為「沒什麼在發生」,Fig. 6 才讓你看見那裡其實是能量輸運完全停止的地方。

補推導:區邊界上相速不為零,群速卻為零補推導

兩個速度一個是零一個不是零,不衝突,但一定要想清楚它們各自在量什麼。

相速量的是「等相位面移動得多快」。在駐波裡,我們仍然可以用行進波的寫法 exp[i(sKa − ωt)] 去描述它,於是「相位」這個記帳量還是在跑。但這是描述方式的產物,不是物理:(12) 已經明白告訴我們波形本身不動(「The wave moves neither to the right nor to the left.」)。相速甚至可以超過光速而不違反相對論,因為它不輸運任何東西。

群速量的是「能量與訊號跑得多快」,Kittel 在 (13) 底下講得很明白:「This is the velocity of energy propagation in the medium.」在區邊界它為零,意思是這個模態完全不參與能量輸運——把它激發起來,能量就停在原地。

這個區分在第 7 章談電子會原封不動再用一次:能隙邊緣的電子態也是駐波(書 p.165 講的就是 k = ±π/a 處入射波與 Bragg 反射波疊加成駐波),群速為零,所以那些態不導電。聲子在區邊界不傳熱,電子在區邊界不導電,是同一個數學事實的兩個物理後果。

  1. 相速的定義是 v_p ≡ ω/K(等相位面的移動速度)。用 (9):v_p = (4C/M)^{1/2} |sin ½Ka| / K。
  2. 在區邊界 K = π/a:sin ½Ka = sin(π/2) = 1,所以 v_p = (4C/M)^{1/2} / (π/a) = (4C/M)^{1/2} · a/π。
  3. 用長波聲速 v₀ ≡ (Ca²/M)^{1/2} 來表示。注意 (4C/M)^{1/2} · a = (4Ca²/M)^{1/2} = 2v₀,所以 v_p(π/a) = 2v₀/π ≈ 0.637 v₀。
  4. 但同一點的群速由 (14) 給出:v_g(π/a) = v₀ cos(π/2) = 0。
  5. 所以在區邊界,相速還有長波聲速的 64%,群速卻是零。兩者差距最大的地方就是這裡。
  6. 順帶一提,在 K → 0 的另一端:v_p = (4C/M)^{1/2}·(½Ka)/K = (Ca²/M)^{1/2} = v₀ = v_g。長波極限下兩者相等——這就是「無色散」的精確定義。

Long Wavelength Limit:兩套語言必須在這裡對得上(書 p.94)

前面三個小節都在看區邊界那一端(K 最大、λ 最短)。這一小節轉過頭看另一端:K → 0。

物理動機非常具體,而且是本書結構上的需要。 第 3 章後半(ANALYSIS OF ELASTIC STRAINS 起,p.73–84;其中 Elastic Waves in Cubic Crystals 在 p.80–84)把晶體當成連續彈性介質處理:定義應力 X_x、應變 e_xx,寫出彈性勁度常數 C₁₁、C₁₂、C₄₄,導出連續體的波動方程,得到 p.82 的 (59) ω²ρ = C₁₁K²、(60) v_s = νλ = ω/K = (C₁₁/ρ)^{1/2}。這幾節裡完全沒有原子——p.73 開宗明義就寫「We consider the elastic properties of a crystal viewed as a homogeneous continuous medium rather than as a periodic array of atoms.」(我們把晶體看成均勻的連續介質、而不是原子的週期性排列,來討論它的彈性性質。)(第 3 章前半 p.49–72 談的是晶體鍵結,與這裡的對帳無關。)

第 4 章反過來,從一開始就把晶體當成一串離散的原子面,用彈簧 C 把它們連起來。裡面完全沒有連續介質。

兩套描述講的是同一塊晶體。如果它們在某處對不上,其中一套就是錯的。Long Wavelength Limit 這一小節就是對帳的地方,而對帳的邏輯很清楚:當波長遠大於原子間距時,波「看不見」個別原子(它只感受到一大群原子的平均行為),晶格模型就必須退化成連續體模型。這不是巧合,是一致性要求。如果 (7) 在 Ka ≪ 1 時給不出 ω ∝ K,那就代表第 4 章的模型從根本上錯了。

反過來說,這個對帳也告訴我們第 3 章的適用範圍在哪裡結束:彈性常數理論只在第一布里淵區中央那一小塊有效,而 Fig. 6 顯示到 K = π/2a 時群速已經偏離 29%。所以「連續體近似」不是在區邊界才失效,它在區的中途就已經明顯失準了。

「Ka ≪ 1」到底是多小? 整段的展開 cos Ka ≡ 1 − ½(Ka)² 截斷誤差是 +(Ka)⁴/24,把幾個 Ka 值代進去就知道這個「長波」的門檻在哪裡:

Ka對應波長 λ/a = 2π/(Ka)cos Ka 真值1 − ½(Ka)²相對誤差
0.2310.980070.980000.007%
0.512.60.877580.875000.29%
1.06.30.540300.500007.5%
π/240−0.2337完全失效

(表中的數值為本文計算,書上沒有給。)結論是:「Ka ≪ 1」在實務上大約是指 λ ≳ 10a,也就是 |K| ≲ 0.2 × (π/a)——第一布里淵區中央約 ±20% 的一小塊。第 3 章的整套彈性常數理論、以及第 5 章德拜模型的低溫 T³ 定律,都是在這一小塊裡建立起來的。這也解釋了為什麼低溫比熱的 T³ 定律只在遠低於德拜溫度時成立:溫度一高,被熱激發的聲子的 K 就跑出這一小塊了。

長波展開所用的餘弦近似

cos Ka ≡ 1 − ½(Ka)²

p.94

K 波向量,a 面間距,Ka 為相鄰兩面的相位差。原文脈絡:「When Ka ≪ 1 we expand cos Ka ≡ 1 − ½(Ka)², so that the dispersion relation (7) becomes」。此式在原書未編號。

出處推導
cos x 的泰勒展開 cos x = 1 − x²/2! + x⁴/4! − … 取到二階項,令 x = Ka。一階項不存在,因為 cos 是偶函數——這一點後面很重要:它保證了 ω² 的領導項是 K² 而不是 K。
極限行為

截斷誤差是 +x⁴/24。量化一下就知道「長波」在實務上是多長:Ka = 0.2(λ ≈ 31a)時相對誤差約 0.007%;Ka = 0.5(λ ≈ 12.6a)時約 0.29%;Ka = 1.0(λ ≈ 6.3a)時已經到 7.5%;Ka = π/2(λ = 4a)時真值是 0 而近似值是 −0.234,完全失效(詳細數字見「Long Wavelength Limit」那一段的表)。所以「Ka ≪ 1」在實務上大約是指 λ ≳ 10a,也就是 |K| ≲ 0.2 × (π/a),第一布里淵區中央大約 ±20% 的一小塊。第 3 章的彈性理論,以及第 5 章德拜模型的低溫 T³ 定律,都是在這一小塊裡建立起來的。

【符號校註】 書上此處印的是恆等號 ≡,但這顯然只是取到二階的近似,並非恆等式。同一個展開在 p.98(雙原子鏈的小 Ka 分析)印的卻是 ≅。我們把兩處分別以 6 倍解析度算繪成影像放大目視,確認 p.94 印的是三橫線的 ≡,p.98 印的是上方帶波浪的 ≅——是原書自身的排版不一致,兩處在本專案中都照原樣保留,不作統一「訂正」。(附帶一提,這也暴露了抽取器的一個缺陷:p.94 用的字符是 /H11013,p.98 用的是 /H33361,現行字符表把兩者都映射成 ≡,其中後者應為 ≅。)

補推導:從 (7) 得到 (15)補推導

這個推導只有三行,但它的結論是整本書最常被引用的結論之一。

Kittel 對 (15) 的解讀寫得很好,值得逐字讀:「The result that the frequency is directly proportional to the wavevector in the long wavelength limit is equivalent to the statement that the velocity of sound is independent of frequency in this limit. Thus v = ω/K, exactly as in the continuum theory of elastic waves—in the continuum limit Ka ≪ 1.」(「頻率在長波極限下直接正比於波向量」這個結果,等價於「聲速在此極限下與頻率無關」這個陳述。因此 v = ω/K,與彈性波的連續體理論完全相同——在連續體極限 Ka ≪ 1 之下。)

請注意他的用字是 exactly as in the continuum theory,不是「近似於」。在 K → 0 的極限下,晶格模型與連續體模型不是相似,是同一個理論。第 3 章與第 4 章在這裡無縫接合。

  1. 出發點 p.92 的 (7):ω² = (2C/M)(1 − cos Ka)。
  2. 代入 cos Ka ≈ 1 − ½(Ka)²,括號裡變成 1 − [1 − ½(Ka)²] = ½(Ka)² = ½K²a²。注意常數 1 被完全消掉了——這就是為什麼聲學支在 K → 0 時 ω → 0,而不是趨向某個有限值(那是光學支的行為,下一段的雙原子鏈才會出現)。
  3. 所以 ω² ≈ (2C/M) × ½K²a² = (C/M)K²a²,即 (15)。注意 2 與 ½ 相消,係數乾淨。
  4. 開根號:ω = a(C/M)^{1/2} · K。ω 正比於 K,比例常數 a(C/M)^{1/2} 完全不含 K。
  5. 從 (9) 走一遍檢查:ω = (4C/M)^{1/2}|sin ½Ka| → (4C/M)^{1/2} · ½Ka = (4C/M · a²/4)^{1/2} K = (Ca²/M)^{1/2} K。同一個答案。✓
  6. 與 (14) 對照:v_g(K → 0) = (Ca²/M)^{1/2} cos 0 = (Ca²/M)^{1/2},而相速 v_p = ω/K = a(C/M)^{1/2} = (Ca²/M)^{1/2}。兩者相等——這就是「無色散」的精確定義,也是「聲速」這個單一數字能夠代表整個低頻行為的理由。

長波(連續體)極限下的色散關係

ω² = (C/M)K²a²

(15)  p.94

ω 角頻率;C 每個原子的面間力常數;M 原子質量;K 波向量;a 面間距。右邊可以整理成 ω = v K,其中 v = a(C/M)^{1/2} = (Ca²/M)^{1/2}。

出處推導
把 cos Ka ≡ 1 − ½(Ka)² 代入 p.92 的 (7) ω² = (2C/M)(1 − cos Ka)。適用條件是原文明寫的 Ka ≪ 1。
極限行為

這一式的意義不在它有什麼,而在它沒有什麼:右邊除了 K² 之外,沒有任何其他對 K 的依賴。因此 v = ω/K = a(C/M)^{1/2} 是一個常數,與頻率無關、與波長無關。這就是「聲音在固體裡不會走音」的微觀理由——低頻的低音與高頻的高音以相同的速度抵達,所以隔著一堵牆聽到的聲音只是變悶(高頻被吸收),而不會變調。

把它推到另一端就知道極限在哪裡:K 增大時 (9) 的 sin ½Ka 開始飽和,ω 不再跟著 K 線性成長,v = ω/K 隨頻率下降,聲音開始色散——不同頻率成分以不同速度傳播,波包會散開。在區邊界,v_p 掉到長波值的 2/π ≈ 64%,而 v_g 掉到 0。

與第 3 章的接點:本式的 v 必須等於 p.82 (60) 的 v_s = (C₁₁/ρ)^{1/2}。這個等式可以反過來用,把唯象的 C 換算成可量測的 C₁₁(見下一段補推導)。

典型數值
銅沿 [100] 的縱波:由書本 Table 11(p.84)取 C₁₁ = 1.762 × 10¹¹ N/m²、密度 ρ = 9.018 g/cm³ = 9018 kg/m³(皆為 0 K 的數值),得 v = (C₁₁/ρ)^{1/2} = (1.762 × 10¹¹ / 9018)^{1/2} ≈ 4.4 × 10³ m/s。這就是 (15) 裡的 a(C/M)^{1/2}。作為對照,Table 11 同頁的鎢 C₁₁ = 5.326 × 10¹¹ N/m²、ρ = 19.317 g/cm³ 給 v ≈ 5.3 km/s;鉛 C₁₁ = 0.555 × 10¹¹ N/m²、ρ = 11.599 g/cm³ 給 v ≈ 2.2 km/s。鎢的勁度是鉛的 9.6 倍、密度只有 1.7 倍,勁度贏過密度,所以聲速高出約 2.4 倍——符合 v ∝ (勁度/密度)^{1/2}。(注意鎢其實比鉛重,它聲速高完全是 C₁₁ 撐起來的,密度反而是往反方向作用。)

補推導:把 (15) 裡的唯象常數 C 接回第 3 章可量測的 C₁₁補推導

書上沒有做這個換算,但它回答了一個讀第 4 章時一定會冒出來的問題:C 到底是什麼?它可以量嗎? 答案是可以,而且上面這條式子就是換算關係。

這個推導還順便解釋了 p.91 那句容易被跳過的叮嚀:「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 是每個原子的力常數,不是「整個面」的。如果把 C 理解成整面的力常數,上面第 5 步就會多出一個面積因子,換算立刻錯掉。

代入銅的數值(fcc,沿 [100] 縱波):

量數值來源
C₁₁1.762 × 10¹¹ N/m²書本 Table 11,p.84(0 K)
ρ9.018 g/cm³書本 Table 11,p.84
a₀3.61 Å外部標準值(本書的晶格常數表在 p.20–21,該處無文字層、抽不到,故不引用書上數字)
M63.5 u = 1.055 × 10⁻²⁵ kg外部標準值

於是:

v = (C₁₁/ρ)^{1/2} ≈ 4.4 × 10³ m/s

C = C₁₁ a₀ = 1.762 × 10¹¹ × 3.61 × 10⁻¹⁰ ≈ 64 N/m

ω_max = 2(C/M)^{1/2} = 2 × (64 / 1.055 × 10⁻²⁵)^{1/2} ≈ 4.9 × 10¹³ rad/s ⇒ ν_max ≈ 7.8 THz

ħω_max ≈ 1.055 × 10⁻³⁴ × 4.9 × 10¹³ ≈ 5.2 × 10⁻²¹ J ≈ 32 meV,換成溫度是 ħω_max/k_B ≈ 375 K。

合理性檢查:銅的德拜溫度約 340 K(外部標準值,本書的德拜溫度表無文字層,未引用),與上面用最粗糙的一維最近鄰模型算出的 375 K 只差一成。一個只有兩個參數、把三維晶體壓成一條原子鏈的模型能做到這個程度,說明「原子質量 + 面間彈簧常數」確實抓到了晶格振動的主要物理。反過來說,剩下那一成的差距,正是後面幾頁要處理的東西——不只最近鄰、不只一個原子基底、不只一個方向。

  1. 第 3 章 p.82 的 (59)(60):沿 [100] 的縱波在連續彈性體裡滿足 ω²ρ = C₁₁K²,因此 v_s = νλ = ω/K = (C₁₁/ρ)^{1/2}。這裡 ρ 是質量密度,C₁₁ 是彈性勁度常數,兩者都可以用超音波實驗直接量到。
  2. 第 4 章 (15) 給同一個波:v = ω/K = a(C/M)^{1/2}。
  3. 既然描述的是同一個物理量,必須相等:a²C/M = C₁₁/ρ。
  4. 把 ρ 用微觀量表示。一個原子佔的體積 = 面間距 a × 該原子在面內佔的面積 S,所以 ρ = M/(aS)。
  5. 代入:a²C/M = C₁₁ · aS/M。兩邊的 M 消掉,再消一個 a,得 C = C₁₁ · S/a。
  6. 量綱檢查:[C₁₁] = N/m²,[S/a] = m²/m = m,所以 [C] = N/m。✓ 與 p.91 的 (1) F_s = C(u_{s+1} − u_s) + C(u_{s−1} − u_s) 所要求的「力/長度」一致。
  7. 套到 fcc 沿 [100]:{100} 面的間距是 a = a₀/2(a₀ 為傳統立方晶胞邊長),每個 {100} 面上每個原子佔的面積是 S = a₀²/2,於是 S/a = (a₀²/2)/(a₀/2) = a₀。得 C = C₁₁ a₀。

Derivation of Force Constants from Experiment:最近鄰模型在金屬裡不夠用(書 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. Interactions have been found between planes of atoms separated by as many as 20 planes.」(在金屬中有效作用力可以是相當長程的,並透過傳導電子海由一個離子傳到另一個離子。已經發現相隔多達 20 個原子面之間仍有交互作用。)

回想 p.91 的設定:「For brevity we consider only nearest-neighbor interactions, with p = ±1.」——只考慮最近鄰,而且他明說理由是 for brevity(為了簡潔)。現在 Kittel 把帳算回來:在金屬裡這個簡化與事實差得很遠。以銅的 {100} 面間距 1.8 Å 估算,20 個面就是 36 Å 的作用範圍,是最近鄰面間距的二十倍(若換算成最近鄰原子距離 a₀/√2 ≈ 2.55 Å,約為十四倍)。

為什麼金屬會這樣? 原文只給了機制的名字:透過傳導電子海。機制本身是:一個離子移動會擾動它周圍的電子密度,而電子氣會重新分布來屏蔽這個擾動;這個重新分布不是在離子附近就結束,它會一路延伸出去,於是遠處的離子也「感覺到」第一個離子動了。離子之間並沒有直接的長程力,長程性完全是電子海這個中介造成的。絕緣體裡沒有傳導電子,作用力主要靠短程的閉殼層排斥與離子靜電力,衰減快得多。

Kittel 在本章習題 4(p.103,題名 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.」注意這個形式只隨 1/p 衰減(外面還乘一個振盪的 sin),不是指數衰減。1/p 衰減正是「20 個面之後還沒消失」的數學來源。習題接著要你證明 ∂ω²/∂K 在 K = k₀ 處發散,也就是色散曲線在 k₀ 處有一個垂直切線的扭結(kink)——這就是 Kohn 異常。

(補充說明:習題把常數寫成 k₀ 而沒有指明它是什麼;在自由電子理論裡它對應 2k_F,k_F 是第 6 章的費米波向量,因此 Kohn 異常的位置可以反過來用來量金屬的費米面。Kittel 在這幾頁與該習題中都沒有講到這一步,這一句是本文補充的背景。)

所以問題變成:既然不能假設任何特定的力常數模型,有沒有辦法直接從實驗量到的色散曲線把力常數一個一個讀出來?接下來的 (16a)(16b)(17) 就是這個方法,而且它是整章唯一一條「從實驗反推微觀參數」的式子。

補推導:(16a) 的「easily found」到底怎麼 found補推導

Kittel 只寫了「The generalization of the dispersion relation (7) to p nearest planes is easily found to be」就跳到答案。整段代數只有六行,但值得自己寫一次,因為 (16b) 的積分完全建立在 (16a) 的形式上——如果 (16a) 的括號不是 (1 − cos pKa),後面正交性那一步會得到不同的結果(特別是那個常被誤以為是筆誤的負號)。

兩個容易搞錯的地方:

第一,求和只跑 p > 0。+p 與 −p 的貢獻已經在第 1 步合併了(u_{s+p} 與 u_{s−p} 共用同一個 C_p,因為晶格有反演對稱),所以 C_p 是「相隔 p 個面的那一對面之間、以每個原子為單位的力常數」,每一對只算一次。如果讓 p 跑遍所有非零整數又不除以 2,答案會差一倍。

第二,(16a) 保證了 ω² 是 K 的偶函數,且以 2π/a 為週期(每一項 cos pKa 都是偶的、週期 2π/(pa),其最小公週期為 2π/a)。這兩個性質正是下一步能做傅立葉反轉的前提:一個定義在 −π/a < K ≤ π/a 上的週期偶函數,可以唯一地展開成 cos pKa 的級數,而展開係數就是力常數。(16a) 不只是一條色散關係,它是一個傅立葉級數。

  1. 把 p.91 的 (1) 推廣。原本只有 s ± 1 兩個鄰居,現在讓所有 s ± p 的面都對 s 施力,每一組鄰距用自己的力常數 C_p,力仍然正比於位移差(線性彈性假設不變):F_s = Σ_{p>0} C_p [(u_{s+p} − u_s) + (u_{s−p} − u_s)] = Σ_{p>0} C_p (u_{s+p} + u_{s−p} − 2u_s)。
  2. 運動方程(p.91 的 (2) 的推廣):M d²u_s/dt² = Σ_{p>0} C_p (u_{s+p} + u_{s−p} − 2u_s)。
  3. 代入行進波解 u_s = u exp[i(sKa − ωt)]。位移的二次時間微分給 −ω²u_s;而 u_{s±p} = u exp[i((s±p)Ka − ωt)] = u_s · exp(±ipKa)。
  4. 兩邊同除 u_s(時間因子與 exp(isKa) 一起消掉):−Mω² = Σ_{p>0} C_p [exp(ipKa) + exp(−ipKa) − 2]。
  5. 用歐拉恆等式 exp(ix) + exp(−ix) = 2 cos x:−Mω² = Σ_{p>0} C_p (2 cos pKa − 2) = −2 Σ_{p>0} C_p (1 − cos pKa)。
  6. 兩邊除以 −M:ω² = (2/M) Σ_{p>0} C_p (1 − cos pKa),即 (16a)。
  7. 自洽檢查:只保留 p = 1 且記 C₁ ≡ C,得 ω² = (2C/M)(1 − cos Ka),正是 p.92 的 (7)。✓ 所以 (16a) 確實是 (7) 的推廣,沒有寫錯。

推廣到 p 個近鄰面的色散關係

ω² = (2/M) Σ_{p>0} C_p (1 − cos pKa)

(16a)  p.95

C_p 為相隔 p 個面的兩個原子面之間、以一個原子為單位定義的面間力常數;p 為正整數(鄰距,實際距離為 pa);M 原子質量;a 面間距;K 波向量。求和只跑 p > 0。

出處推導
把 p.91 的 (1)(2) 從最近鄰推廣到任意鄰距,再代入行進波解(詳見上一段補推導)。原文的說法是「The generalization of the dispersion relation (7) to p nearest planes is easily found to be」。
極限行為

K → 0(長波極限):每一項 1 − cos pKa → ½(pKa)²,所以 ω² → (a²K²/M) Σ_{p>0} p² C_p,聲速 v = a[Σ_{p>0} p² C_p / M]^{1/2}。注意遠鄰的貢獻被 p² 加權放大:第 20 鄰的力常數即使只有最近鄰的 1%,它對聲速的貢獻也被放大 400 倍而變成可比擬的量。這是為什麼在金屬裡拿實驗聲速去硬套最近鄰公式 v = a(C/M)^{1/2} 會得到嚴重錯誤的 C。

K = π/a(區邊界):cos pπ = (−1)^p,所以 1 − cos pπ 對奇數 p 等於 2、對偶數 p 等於 0,得 ω²_max = (4/M) Σ_{p 為奇數} C_p。偶數鄰居在區邊界完全沒有貢獻——理由很具體:由 (12) 駐波的位移是 u_s = u(−1)^s,相隔偶數個面的兩個面 u_s 與 u_{s+2m} 同號同值,相對位移為零,既然力正比於位移差,這對面之間就沒有力,自然不影響頻率。

把這兩個極限並排看很有啟發性:長波聲速偏重遠鄰(p² 加權),區邊界頻率只看奇數鄰。同一組 C_p,在色散曲線的兩端以完全不同的組合方式現身,所以量一整條 ω(K) 才能把它們分離出來——這正是 (17) 在做的事。

補推導:「The integral vanishes except for p = r」底下的兩個積分補推導

書上這一整套只有一句「The integral vanishes except for p = r.」,是全章省略得最兇的地方之一。補出來之後可以看到,這裡沒有任何深奧的東西,就是餘弦函數在一個完整週期上的正交性。但不補的話,讀者會卡在兩個點:為什麼常數 1 那一項不見了,以及為什麼會冒出一個負號。

把整件事看成傅立葉變換對。 (16a) 是正向變換:給定一組力常數 {C_p},合成出色散曲線 ω²(K)。(17) 是逆變換:給定實驗量到的色散曲線 ω²(K),分解出力常數 {C_p}。能這樣做的條件在上一段已經備齊——ω²(K) 是定義在第一區上、以 2π/a 為週期的偶函數,所以它的完整傅立葉基底就是 {1, cos Ka, cos 2Ka, cos 3Ka, …},而力常數就是展開係數(差一個 −2/M 的因子)。

這個觀點也解釋了為什麼必須量「整條」曲線。 傅立葉係數是全域資訊:要把 C_p 定出來,積分必須跑遍整個第一區。只量長波段(超音波實驗做得到的區域)只能定出 Σ p²C_p 一個組合數,無法把個別的 C_p 分開。要分開就必須量到區邊界,而那需要波長 ~2a 量級的探針——這正是為什麼聲子色散關係要用中子非彈性散射來量(熱中子的德布羅意波長恰好是 Å 量級,能量又恰好是 meV 量級,波長與能量同時匹配,見 p.100–102)。

  1. 整段的策略是傅立葉級數的標準伎倆:想取出第 r 個係數,就乘上第 r 個基底函數再積分,讓其他所有項自己消掉。Kittel 說得很白:「We solve for the interplanar force constants C_p by multiplying both sides by cos rKa, where r is an integer, and integrating over the range of independent values of K」。積分範圍就是第一布里淵區 −π/a 到 π/a,長度 2π/a。
  2. 換元讓算式乾淨:令 x = Ka,則 dK = dx/a,積分範圍變成 −π 到 π。
  3. 第一個積分(來自 (16a) 括號裡那個常數 1):∫_{−π/a}^{π/a} cos rKa dK = (1/a)∫_{−π}^{π} cos rx dx = (1/a)[sin rx / r]_{−π}^{π} = (1/a)·(sin rπ − sin(−rπ))/r = 0,因為 r 是非零整數、sin(整數 × π) = 0。物理上:常數項在區內對任何非零諧波的投影為零。
  4. 第二個積分(來自 cos pKa × cos rKa):用積化和差 cos A cos B = ½[cos(A − B) + cos(A + B)]:∫ cos pKa cos rKa dK = ½ ∫ [cos(p − r)Ka + cos(p + r)Ka] dK。
  5. 由第 3 步的結果,∫_{−π/a}^{π/a} cos mKa dK = 0 對任何非零整數 m 成立;而 m = 0 時被積函數是 1,積分等於區間長度 2π/a。
  6. 因為 p、r 都是正整數,p + r 永遠 ≥ 2 不為零,第二項恆為零。第一項只有在 p = r 時 m = 0 而不為零,此時值為 ½ × (2π/a) = π/a。
  7. 合起來:∫_{−π/a}^{π/a} cos pKa cos rKa dK = (π/a) δ_{pr}。這就是 Kittel 那一句話的全部內容。
  8. 現在把 (16a) 兩邊乘 M cos rKa 再對 K 積分。左邊 = M ∫_{−π/a}^{π/a} dK ω²_K cos rKa。右邊 = 2 Σ_{p>0} C_p ∫_{−π/a}^{π/a} dK (1 − cos pKa) cos rKa = 2 Σ_{p>0} C_p [0 − (π/a)δ_{pr}] = −2πC_r/a。這就是 (16b)。
  9. 負號的來源:(16a) 括號裡是 (1 − cos pKa)。常數 1 的貢獻被第 3 步殺成零,活下來的是 −cos pKa 這一項,負號就是從這裡帶出來的。(16b) 與 (17) 前面的負號完全是這麼來的,不是筆誤。
  10. 最後把 (16b) 解出 C_r(再把啞指標 r 改回 p):C_p = −(a/2π) · M ∫_{−π/a}^{π/a} dK ω²_K cos pKa = −(Ma/2π) ∫_{−π/a}^{π/a} dK ω²_K cos pKa,即 (17)。

乘上 cos rKa 後對第一區積分

M ∫_{−π/a}^{π/a} dK ω²_K cos rKa = 2 Σ_{p>0} C_p ∫_{−π/a}^{π/a} dK (1 − cos pKa) cos rKa = −2πC_r/a

(16b)  p.95

ω_K 表示波向量為 K 的模態頻率(原書以 ω 帶下標 K 印出);r 為整數(原文只說 where r is an integer;因 (16a) 的求和僅跑 p > 0,實際取用的是正整數),是我們想取出的那個鄰距;其餘符號同 (16a)。

出處推導
把 (16a) 兩邊乘 M cos rKa,再對第一布里淵區積分。第二個等號用到餘弦正交性(原文只寫「The integral vanishes except for p = r.」,完整計算見上一段補推導)。
極限行為

這一式本身不是最終答案,它是一個中間步驟,但它的形式值得看清楚:左邊完全由實驗量(M 與整條 ω(K) 曲線)構成,右邊只剩下單獨一個 C_r。無窮多個未知數 {C_1, C_2, C_3, …} 被一次抽出一個——這就是正交展開的威力,也是為什麼傅立葉方法在物理裡無所不在。

如果實驗的 ω(K) 只在有限個 K 點上量到(實際上永遠如此),這個積分就只能用數值方法近似。此時能可靠取出的最大 p 受限於 K 點的密度:要解析出鄰距 p 的成分,K 的取樣間隔必須細於 π/(pa)。所以「在金屬中發現相隔 20 個面仍有交互作用」這個結論,背後代表的是一套 K 點取樣夠細、統計夠好的中子散射實驗。

由實驗色散曲線反演出面間力常數

C_p = − (Ma/2π) ∫_{−π/a}^{π/a} dK ω²_K cos pKa

(17)  p.95

C_p 為相隔 p 個面的面間力常數(作用範圍 pa,原文:「gives the force constant at range pa, for a structure with a monatomic basis」);M 原子質量;a 面間距;ω²_K 為實驗量到的色散關係的平方;積分跑遍第一布里淵區。

出處推導
把 (16b) 解出 C_r 並把啞指標改寫成 p。前置因子 −(Ma/2π) 由 (16b) 右邊的 −2π/a 取倒數乘 M 得到。
極限行為

這是本段真正的成果,而且它的性質值得逐條列出:

(一)不預設任何力學模型。 式子裡沒有「最近鄰」「彈簧」「Lennard-Jones」這些假設,只要你量得到整條 ω(K),就能把力常數讀出來。這是一個實驗資料到微觀參數的直達管道。

(二)逐一取出。 改變 p 就換一個鄰距:p = 1 給最近鄰面、p = 2 給次近鄰面……實驗上發現 C_p 到 p = 20 還不可忽略,就知道最近鄰模型在該材料失效。反過來,如果量出來的 C_p 在 p ≥ 2 都幾乎為零,就證實了最近鄰近似在該材料成立。

(三)逆向檢查。 若 ω²(K) 恰好只含 cos Ka 這一個諧波成分,則除 C_1 外全部為零,回到 (7)。詳細驗算見下一段。

(四)適用限制。 原文最後一句寫明「for a structure with a monatomic basis」——單原子基底。基底裡有兩個以上原子時色散關係會分裂成聲學支與光學支,(17) 這個單一積分不再夠用(那是 p.95 下半 TWO ATOMS PER PRIMITIVE BASIS 開始的故事)。

補推導:用最近鄰模型自檢 (17)(強烈建議自己算一次)補推導

這個檢驗花不到五分鐘,但它能一次抓到三種常見錯誤:

負號記錯:如果把 (17) 前面的負號漏掉,上面會算出 C₁ = −C,力常數變成負的(意味著晶格不穩定),一眼就看得出不對。

係數 a 與 2π 的位置弄反:如果寫成 −(M/2πa) 或 −(2πM/a),C₁ 就不會乾乾淨淨地等於 C,而會多出 a² 之類的因子,量綱也會壞掉。(量綱檢查:[M] = kg,[a] = m,[dK] = m⁻¹,[ω²] = s⁻²,所以 [Ma · dK · ω²] = kg·m·m⁻¹·s⁻² = kg·s⁻² = N/m。✓ 正是力常數的量綱。)

積分範圍記錯:如果誤用 0 到 π/a(只積半個區),正交性積分的值會減半,C₁ 會變成 C/2。

養成「推完一條反演公式,立刻拿一個已知解代回去檢驗」的習慣,在後面幾章(特別是第 6 章之後大量出現的積分變換)會省下非常多的時間。

  1. 拿一個我們已經知道答案的例子去檢驗 (17):把 p.92 的 (7) ω²_K = (2C/M)(1 − cos Ka) 代進去,看能不能把 C 吐回來。
  2. 先算 p = 1:C₁ = −(Ma/2π) ∫_{−π/a}^{π/a} (2C/M)(1 − cos Ka) cos Ka dK。常數 M 消掉,提出 (2C):C₁ = −(aC/π) ∫_{−π/a}^{π/a} (1 − cos Ka) cos Ka dK。
  3. 展開括號成兩個積分:∫ cos Ka dK − ∫ cos²Ka dK。第一個由正交性為 0;第二個 ∫_{−π/a}^{π/a} cos²Ka dK = π/a(即 δ_{11} 的情形)。
  4. 所以 C₁ = −(aC/π) × (0 − π/a) = −(aC/π) × (−π/a) = C。✓ 正確吐回最近鄰力常數,連負號都對上了。
  5. 再算 p = 2:C₂ = −(aC/π) ∫_{−π/a}^{π/a} (1 − cos Ka) cos 2Ka dK = −(aC/π)[∫cos 2Ka dK − ∫cos Ka cos 2Ka dK] = −(aC/π)[0 − 0] = 0。✓ 最近鄰模型確實沒有次近鄰貢獻,符合預期。
  6. 同理所有 p ≥ 2 都給零,因為 (7) 的傅立葉展開裡只有 cos Ka 這一個諧波。

這三頁最常見的五個誤解我的判讀

誤解一:「第一布里淵區之外的 K 不存在。」

不對。區外的 K 在數學上完全合法,exp(iKx) 對任何 K 都是一個好好的函數。問題在於它不是新的物理狀態:它在原子位置上的取樣值與區內某個 K 完全相同,頻率、群速、能量也全部相同。Fig. 5 的白色實線是一條真實存在的曲線,只是它相對於虛線多出來的資訊,全部落在「原子與原子之間」,而那裡沒有自由度可以承載資訊。Kittel 的用字是 reproduce(重現),不是 forbidden(禁止)。

誤解二:「區邊界群速為零,表示原子不動。」

不對。(12) 的駐波裡每一個原子都在以 ω_max 全幅振盪,振幅仍然是 u,而且這是整個色散關係裡頻率最高的模態。為零的是能量的淨傳輸速度,不是振幅。打個比方:兩端固定的吉他弦上的駐波,弦在劇烈振動,但波形不往任何一邊跑。

誤解三:「長波極限 (15) 只是近似,所以由它推出的聲速也只是近似。」

要分清楚兩件事。(15) 本身確實是把 cos Ka 截斷到二階的近似,誤差 O((Ka)⁴)。但它推出的陳述「v = ω/K 與頻率無關」是一個極限陳述:在 K → 0 時嚴格成立。Kittel 的用字是「exactly as in the continuum theory of elastic waves」——exactly,不是 approximately。第 3 章整套彈性常數理論就建立在這個極限上,如果它只是近似,那麼所有用超音波量出來的彈性常數都會是可疑的,而事實並非如此。

誤解四:「(17) 前面那個負號是筆誤。」

不是。負號來自 (16a) 括號裡的 −cos pKa:常數項 1 在積分中被正交性殺掉,活下來的只有帶負號的那一項。可以用最近鄰模型自檢(上一段),把負號拿掉會得到 C₁ = −C,力常數變負、晶格不穩定,馬上知道錯了。

誤解五:「±π/a 是兩個不同的態,只是頻率碰巧相同。」

不對。它們相差一個倒晶格向量 2π/a,由 (11) 是同一個態。這就是為什麼獨立 K 的範圍要寫成半開區間 −π < Ka ≤ π,以及為什麼 (12) 要把兩個符號寫在一起成 u exp(±isπ)——他在說這兩個寫法給同一個答案。在區的內部,+K 與 −K 確實是兩個不同的態(往右跑與往左跑的行進波);只有在端點,這個「方向」的概念隨著行進波退化成駐波而一併消失。

這三頁在全書的位置:往前接什麼、往後被誰用我的判讀

往前接:

第 2 章(p.34–37)給了布里淵區的幾何定義——倒晶格向量的垂直平分面圍出的最小體積,也就是倒晶格的 Wigner-Seitz 原胞;一維的情形在 Fig. 11(p.35)畫出邊界在 k = ±π/a。第 2 章是從 X 光繞射的需求引入這個概念的,當時它看起來像一個純幾何的輔助工具。本段給了它第一個與繞射無關的物理用途,並且在 p.94 把兩者接上(區邊界駐波 ≡ Bragg 反射)。

第 3 章(p.80–84)的連續體彈性理論,在本段的 (15) 處與晶格模型接合。接合點的代數(C = C₁₁ S/a)書上沒寫,本文補在前面。接合之後可以看得很清楚:第 3 章不是一個「錯的但方便的」理論,它是第 4 章在 Ka ≪ 1 的精確極限。

往後被用:

第 4 章後半(p.95 起):TWO ATOMS PER PRIMITIVE BASIS。第一區、區邊界、群速這三個工具原封不動搬過去,只是色散關係分裂成聲學支與光學支。(17) 的反演公式則明文限定在 monatomic basis,不能直接搬。

第 5 章(Phonons II):本段的每一個結論都變成計算工具。允許的 K 值填滿第一區,成為模態計數的基礎(p.109–111);群速出現在態密度的積分式裡(p.118,「the magnitude of the group velocity of a phonon」);而區邊界群速為零這件事,直接造成態密度的發散——p.119 逐字:「Such critical points produce singularities (known as Van Hove singularities) in the distribution function」。德拜模型則是把本段的 (15)(ω = vK 線性色散)一路外推到截止頻率,Fig. 6 告訴我們這個外推在區的中段就已經高估群速 29%,這就是德拜模型在中溫區偏差的來源。

第 7 章(Energy Bands):同一個 K = ±π/a 的 Bragg 條件,這次作用在電子的德布羅意波上。書 p.165 逐字:「the Bragg reflection condition k = ±π/a is satisfied by the wavevector, a wave traveling to the right is Bragg-reflected to travel to the left, and vice versa.」——與 p.94 那段講聲子的話講的幾乎是同一句物理,只是把 K 換成 k、把彈性波換成電子波。結果是能隙的出現。p.177 起正式引入 reduced zone scheme,那就是把本段 (11) 的「折回操作」制度化。

第 9 章(Fermi Surfaces):p.223–228 反覆使用 reduced、periodic、extended 三種區表示法,全部建立在 (11) 這三行代數上。

一句話總結本段:晶體是離散的,所以描述晶體的波向量是週期性的;週期性的基本區間就是第一布里淵區;區邊界對應最短波長 2a,在那裡行進波退化成駐波、能量停止傳播;區中央對應長波,在那裡晶格退化成連續體、第 3 章接手。剩下的整本書,大半是在這個框架裡填不同的內容。