Kittel·固態物理

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

03每個原胞兩個原子(分支計數與久期方程)

TWO ATOMS PER PRIMITIVE BASIS (branch counting and the secular equation)

書本 p.96-98 · 小節 3 / 5 · 概念 19 推導 9 公式 10

單原子鏈只有一條色散曲線;一旦原胞裡放進第二個原子,色散關係就分裂成聲學支與光學支——這一段先用自由度把「到底該有幾條支」算死,再把雙原子鏈的久期方程從頭解到底,並解釋「光學」這個名字是怎麼來的。

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

第 4 章到 p.95 為止,Kittel 處理的一直是單原子基底(monatomic basis):每個原胞一個原子,所有平面等價,運動方程只有一條 (2),色散關係只有一條 (7)/(9):ω² = (4C/M)sin²(½Ka)。這條曲線在 K → 0 是線性的(聲速),在區邊界 K = π/a 斜率為零(駐波)。整件事乾淨得有點可疑——因為真實晶體幾乎沒有一個是這樣。

問題出在哪裡。NaCl、KBr、金剛石、矽、鍺、GaAs——這些晶體的原胞裡都有兩個原子。一旦如此,「平面」就不再只有一種:沿著某些方向看過去,你會看到 A 面、B 面、A 面、B 面交替排列。這時候 A 面與 B 面可以同相動,也可以反相動,而這兩種運動的頻率完全不同。單一條色散曲線裝不下這兩種行為。

這一段要建立三件事:

(一)數目。原胞有 p 個原子時,色散關係有幾條分支?答案是 3p 條,其中 3 條聲學、3p − 3 條光學。這個結論不需要解任何方程,純粹是自由度計數——而且每一步都可以驗算。這是本段第一個主題(書本 p.95 末的宣告 + p.96 的 numerology 段 + 腳註 1)。

(二)機制。為什麼會分裂成兩支?這需要真的把雙原子鏈解出來。Kittel 從 (18) 的兩條耦合運動方程出發,用試解 (19) 化成 (20) 的 2×2 線性系統,要有非零解就必須讓係數行列式為零 (21),展開得到久期方程 (22)。這是整章代數最密集的一段,而書上跳了兩大步沒寫。

(三)名字。「光學支」這個名字不是隨便取的。把 K = 0 的光學支解回代,會得到 u/v = −M₂/M₁ (26):兩個子晶格反相振動、質心不動。若兩個原子帶相反電荷,這個運動就是一個振盪的電偶極矩,可以被光波的電場直接驅動——這就是「光學」的由來,也是第 14 章 polariton 與 LST 關係的起點(書本 p.414)。

為什麼這一段值得花時間。後面所有關於熱容、熱傳導、紅外光學性質的討論,都建立在「一個晶體有 3pN 個振動模、分成聲學與光學兩類」這個骨架上。第 5 章的 Debye 模型只處理聲學支(線性色散),Einstein 模型則是光學支的極端理想化(頻率與 K 無關)——這兩個模型為什麼長那個樣子,答案全在這一段的 (23) 與 (24) 裡。

書上的宣告:p 個原子 ⇒ 3p 條分支

這一段的主結論,Kittel 在 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. Thus germanium (Fig. 8a) and KBr (Fig. 8b), each with two atoms in a primitive cell, have six branches: one LA, one LO, two TA, and two TO.」

注意這段話的結構:先給普遍公式(3p = 3 + (3p − 3)),再立刻用兩個真實材料驗算。p = 2 時,3p = 6,3 條聲學 + 3 條光學。Kittel 把這 6 條再按偏振拆開:一條 LA、一條 LO、兩條 TA、兩條 TO。

這個拆法背後的規則在 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).」——每個波向量有三個偏振:一個縱、兩個橫。這是三維空間裡位移向量的三個分量。所以 3p 條分支 = 3 個偏振 × p 條「每個偏振各自的分支」。

把這兩層合起來(以下這張一般 p 的偏振拆解表是我的推廣,不是書中內容:Kittel 只對 p = 2 寫出「one LA, one LO, two TA, and two TO」;而且「一縱兩橫」的分類本身只在高對稱方向上才嚴格成立):

偏振聲學光學小計
縱(L)1 條 LAp − 1 條 LOp
橫(T,兩個)2 條 TA2(p − 1) 條 TO2p
合計33p − 33p

p = 2 就給出 1 LA + 1 LO + 2 TA + 2 TO = 6 ✓。

「3 條聲學」這件事與 p 無關,這一點值得停下來想一秒:不管原胞塞了 2 個原子還是 20 個原子,聲學支永遠只有 3 條。理由在下一段的自由度計數裡——聲學支對應的是「整個原胞一起平移」的三個方向,原胞裡有幾個原子並不會增加平移的方向數。(這個剛體平移的論證是我的補充,Kittel 在此只給出計數結果,沒有解釋「為什麼恰好 3 條」;下一個 block 會再展開一次。)

自由度計數:3pN 從哪裡來(p.96 的 numerology 段)

Kittel 在 p.96 用一段話把分支數的來源講清楚。原文逐字是:

「The numerology of the branches follows from the number of degrees of freedom of the atoms. With p atoms in the primitive cell and N primitive cells, there are pN atoms. Each atom has three degrees of freedom, one for each of the x, y, z directions, making a total of 3pN degrees of freedom for the crystal. The number of allowed K values in a single branch is just N for one Brillouin zone.」

然後接到 p.97 開頭:

「LA and the two TA branches have a total of 3N modes, thereby accounting for 3N of the total degrees of freedom. The remaining (3p – 3)N degrees of freedom are accommodated by the optical branches.」

把這段拆成可以逐步驗算的四行:

第 1 行(原子數)。 原胞裡 p 個原子 × 晶體裡 N 個原胞 = pN 個原子。這裡的 N 是「原胞數」,不是「原子數」——這是最容易在期末考混掉的地方。

第 2 行(自由度)。 每個原子在三維空間裡有 3 個位移自由度(x、y、z)。所以晶體總共 3pN 個自由度。這是一個純粹的力學計數,跟晶體結構、鍵結、對稱性都無關。

第 3 行(每支的模數)。 單一條分支在一個布里淵區內,允許的 K 值恰好有 N 個。這是 Kittel 加腳註 1 的地方,下一個 block 完整補上。

第 4 行(配平)。 分支數 × 每支的 K 數 = 總自由度:

(分支數) × N = 3pN ⇒ 分支數 = 3p ✓

再把 3p 拆開:3 條聲學支承擔 3N 個自由度,剩下的 3pN − 3N = (3p − 3)N 由光學支承擔,也就是 3p − 3 條光學支各 N 個 K 值。

這個論證的力量在於它不需要解任何方程。你不需要知道力常數、不需要知道質量比、不需要知道晶體結構——只要知道原胞裡有幾個原子,分支數就定死了。這是「計數論證」在固態物理裡的第一個範例,同樣的邏輯後面會一再出現(第 5 章數模數算 Debye 頻率、第 7 章數能帶裡的軌道數)。

為什麼「聲學支恰好 3 條」。 Kittel 這裡沒有直接解釋,但邏輯藏在第 4 行裡:聲學支的定義特徵是 K → 0 時 ω → 0。K = 0 且 ω = 0 意味著整個晶體剛體平移而不產生回復力——這樣的運動在三維空間裡恰好有 3 個(x、y、z 三個平移方向),與 p 無關。所以不管原胞多複雜,ω(K → 0) → 0 的分支永遠是 3 條,其餘 3p − 3 條在 K = 0 都有有限的頻率。

腳註 1 補全:為什麼「一條分支在一個布里淵區內恰好 N 個 K」補推導

上一段第 3 行的「N」看起來像是變魔術變出來的。Kittel 把理由放在 p.96 的腳註 1 裡,而且是前向引用——它用的是第 5 章才會證明的結果。原文逐字是:

「We show in Chapter 5 by application of periodic boundary conditions to the modes of the crystal of volume V that there is one K value in the volume (2π)³/V in Fourier space. The volume of a Brillouin zone is (2π)³/Vc, where Vc is the volume of a crystal primitive cell. Thus the number of allowed K values in a Brillouin zone is V/Vc, which is just N, the number of primitive cells in the crystal.」

這段話其實是三行算術,但它壓縮得太緊,值得攤開。

  1. 第 1 步:K 空間的「配額」。 對體積為 V 的晶體施加週期邊界條件(Born–von Karman),允許的 K 形成一個離散格點。Kittel 說這個格點的密度是「one K value in the volume (2π)³/V」——也就是說,每一個 K 值在傅立葉空間裡佔據 (2π)³/V 的體積。注意這個體積與晶體多大成反比:晶體越大,K 點越密。這正是熱力學極限下 K 可以當連續變數的原因。
  2. 第 2 步:布里淵區的體積。 第一布里淵區是倒晶格的 Wigner–Seitz 原胞,其體積等於倒晶格的原胞體積。第 2 章的結果是:正晶格原胞體積 V_c 的倒晶格原胞體積為 (2π)³/V_c。Kittel 這裡直接引用:「The volume of a Brillouin zone is (2π)³/V_c」。
  3. 第 3 步:相除。 一個布里淵區裡塞得下幾個 K 值?

    (區的體積) ÷ (每個 K 佔的體積) = [(2π)³/V_c] ÷ [(2π)³/V] = V / V_c

  4. 第 4 步:辨認這個數。 V 是晶體總體積,V_c 是一個原胞的體積,所以 V/V_c 就是晶體裡的原胞數 N。(2π)³ 在第 3 步完全消掉了——這是這個論證漂亮的地方:結果與維度常數無關,也與你怎麼定義傅立葉變換的 2π 慣例無關。
  5. 第 5 步:量綱與極限檢查。 V/V_c 是無量綱的純數 ✓。當晶體變大(V → ∞)時 N → ∞,K 點變密但區的大小不變(因為 V_c 固定)——這正確地反映了「布里淵區的大小由晶體結構決定,與樣品大小無關」。

這一步為什麼重要。「每支 N 個模」是整個固態物理最常被使用、卻最少被重新推導的事實。它把「分支」(連續的 ω(K) 曲線)與「模」(離散的量子力學自由度)接起來:一條分支不是一個模,而是 N 個模的集合。第 5 章算熱容時要對模求和,用的就是這個計數;第 5 章的態密度 D(ω) 本質上就是「把這 N 個 K 點按 ω 分箱」。

一個常見的混淆。「N 個 K 值」是每一條分支各自有 N 個,不是所有分支合起來 N 個。3p 條分支各 N 個 K,總共 3pN 個模,與 3pN 個自由度一一對應——這個一一對應正是正規模(normal mode)理論的核心結論:一個 3pN 自由度的耦合諧振子系統,恰好可以對角化成 3pN 個獨立的簡諧振子。

一個布里淵區內的允許 K 值數(腳註 1)

(允許 K 數) = [(2π)³/V_c] ÷ [(2π)³/V] = V/V_c = N

p.96

V:晶體總體積;V_c:一個晶體原胞(primitive cell)的體積;(2π)³/V:傅立葉空間中每一個允許 K 值所佔的體積(由體積 V 的週期邊界條件給出,第 5 章證明);(2π)³/V_c:第一布里淵區的體積;N:晶體中的原胞數。

出處推導
書本 p.96 腳註 1。Kittel 只給了三個事實與最後的結論「which is just N」,中間的相除這一步沒有寫成式子;上式把它補齊。所用假設只有一個:週期(Born–von Karman)邊界條件,這是第 5 章的內容,本章屬前向引用。
極限行為
(2π)³ 在相除時完全消掉——這代表結果與傅立葉變換的 2π 慣例無關,也與空間維度的具體數值常數無關(在 d 維會變成 (2π)^d 相消,結論同樣是 N)。熱力學極限 V → ∞ 時 N → ∞,K 點密度 V/(2π)³ 發散,K 從離散變成連續,求和 Σ_K 可以換成積分 (V/(2π)³)∫d³K——這個替換是第 5 章態密度與 Debye 積分的全部基礎。反過來,對奈米尺度樣品(V 小、N 小),K 點稀疏,這個連續化就會失效,聲子譜出現可分辨的離散化——這是奈米結構熱傳中「聲子侷限(phonon confinement)」效應的源頭。
典型數值
以本段 Fig. 8b 的 KBr 為例:NaCl 結構,傳統立方晶胞邊長 a = 6.59 Å(書本 p.14 的表),每個立方晶胞含 4 個原胞,故 V_c = a³/4 = 6.59³/4 = 71.6 ų。一顆 1 mm³ 的 KBr 單晶:V = 10⁻⁹ m³ = 10²¹ ų,N = 10²¹/71.6 ≈ 1.4×10¹⁹ 個原胞,也就是每一條分支有約 1.4×10¹⁹ 個允許 K 值,六條分支合計約 8.4×10¹⁹ 個模。(V_c 與 N 為我依 p.14 的 a 值計算,非原書數字。)

分支計數驗算表:從 p = 1 到 p = 2,以及它怎麼對回 Fig. 8補推導

把 3p = 3 + (3p − 3) 這條規則對幾個實例走一遍,可以立刻看出它的內容有多硬。

材料結構原胞內原子數 p總分支 3p聲學 3光學 3p − 3按偏振拆
本章前半的單原子鏈每原胞 1 原子13301 LA + 2 TA
鍺 Ge(Fig. 8a)金剛石26331 LA + 2 TA + 1 LO + 2 TO
KBr(Fig. 8b)NaCl26331 LA + 2 TA + 1 LO + 2 TO

第一列的意義。 p = 1 時光學支數 = 3(1) − 3 = 0。這正確地回頭說明了為什麼第 4 章前半(p.91–95)的單原子鏈只有一條色散曲線:它是一維模型,只保留一個偏振,而 p = 1 使光學支全部消失。所以 (9) 的 ω² = (4C/M)sin²(½Ka) 是本段 (22) 在 p = 1 的退化情形——這件事會在後面的「等質量極限」block 被嚴格證明。

Fig. 8 只畫了四條曲線,六條跑到哪去了。 這是讀圖時第一個會卡住的地方。答案是橫波簡併:[111] 是立方晶體的三重旋轉軸,兩個垂直於 [111] 的橫向偏振方向可以被 120° 旋轉互換,因此必須同頻。所以兩條 TA 疊成一條、兩條 TO 疊成一條,圖上看到的四條曲線實際承載六條分支。Kittel 的 Fig. 8a 圖說也是這樣措辭的:「The two TA phonon branches are horizontal at the zone boundary position」——他說「兩條」,但圖上只有一條線。(橫波簡併的對稱性理由是我的補充說明,書中此處沒有明說。)

驗算: 圖上曲線數 4,其中 LA 與 LO 各 1 條不簡併、TA 與 TO 各 1 條代表 2 條 ⇒ 1 + 1 + 2 + 2 = 6 ✓ 與 3p = 6 相符。

若換一個方向會怎樣。 在一般的低對稱方向(例如 [110] 或任意方向),兩個橫向偏振不再等價,簡併解除,圖上就會看到完整的六條曲線。這也是為什麼中子散射實驗一定沿著高對稱方向做——曲線少、指認容易。(此段為我的補充,書中未討論。)

Fig. 8a 逐筆讀圖:鍺在 [111] 方向的聲子色散(80 K)

網頁不放圖,以下用文字重建。圖說逐字是:「Phonon dispersion relations in the [111] direction in germanium at 80 K. The two TA phonon branches are horizontal at the zone boundary position, K_max = (2π/a)(½ ½ ½). The LO and TO branches coincide at K = 0; this also is a consequence of the crystal symmetry of Ge. The results were obtained with neutron inelastic scattering by G. Nilsson and G. Nelin.」

校註:圖說中的 K_max = (2π/a)(½ ½ ½) 在 PDF 文字層被拆成「(2π/a)( 1 2 1 2 1 2)」漂到段末;以上為對該頁算繪影像放大目檢後還原的形式。

座標軸。橫軸是 K/K_max([111] 方向),刻度 0、0.2、0.4、0.6、0.8、1.0。縱軸是「Phonon frequency, in 10¹² Hz」,刻度 0 到 10,每 2 一格。注意縱軸標的是頻率 ν(單位 Hz),不是角頻率 ω;兩者差一個 2π。

四條曲線(以下數值為我對該頁算繪影像的目視讀值,誤差約 ±0.1×10¹² Hz):

分支K = 0區邊界 K/K_max = 1形狀
TO≈ 9.15≈ 8.7幾乎水平,全程只降約 0.5
LO≈ 9.15≈ 7.3單調下降,末端轉平
LA0≈ 6.55從原點線性上升,逐漸彎平
TA0≈ 1.9從原點線性上升,很早就趨平

三個必須讀出來的特徵。

(一)LO 與 TO 在 K = 0 併攏。兩條曲線在縱軸上從同一點 ≈ 9.15×10¹² Hz 出發。Kittel 把理由歸給對稱性:「this also is a consequence of the crystal symmetry of Ge」。這是與 Fig. 8b 的 KBr 最關鍵的差別,下一個 block 會對照。

(二)LO 與 TO 在中途分岔。兩條曲線在 K = 0 到 K/K_max ≈ 0.45 之間簡併,圖上畫成同一條線;到 K/K_max ≈ 0.45(ν ≈ 8.7)附近才分岔,之後 LO 一路下降、TO 幾乎維持水平。兩條線從頭到尾沒有交叉,TO ≥ LO 全程成立。這意味著「光學支」不等於「頻率較高的支」——在區邊界 TO(8.7) 明顯高於 LO(7.3)。分支的標籤來自偏振與 K = 0 的行為,不是來自頻率排序。

(三)區邊界斜率為零。TA 曲線在 K/K_max = 1 明顯水平(圖說特別點名),LA、LO 在末端也轉平。這與單原子鏈 (8) 的結論一致:dω/dK = 0 在區邊界,群速度為零,該處是駐波(p.93 的 (12))。

(四)聲學/光學之間的頻隙。在區邊界,最高的聲學支 LA 到 ≈ 6.55,最低的光學支 LO 在 ≈ 7.3——兩者之間有約 0.7×10¹² Hz 的間隙。這對應 Fig. 7 與 (25) 所描述的頻隙結構(但 Ge 的兩個原子質量相同,這件事後面會專門討論)。

Fig. 8b 逐筆讀圖:KBr 在 [111] 方向的色散(90 K),以及 ω_T、ω_L 的定義

網頁不放圖,以下用文字重建。圖說逐字是:「Dispersion curves in the [111] direction in KBr at 90 K, after A. D. B. Woods, B. N. Brockhouse, R. A. Cowley, and W. Cochran. The extrapolation to K = 0 of the TO, LO branches are called ω_T, ω_L.」

與 Fig. 8a 的呈現差異。Fig. 8a 是畫成連續曲線的(平滑過的結果),Fig. 8b 則是散佈的實驗資料點,沒有連線。點分兩種:空心圓與實心圓。從標籤位置判斷,空心圓是縱波(LO 與 LA),實心圓是橫波(TO 與 TA)。縱軸同樣是「Phonon frequency, in 10¹² Hz」,但範圍只到 5,刻度 0、1、2、3、4、5。橫軸是 K/K_max([111] 方向),只在原點標了 0,其餘只有刻度線,沒有數字。

資料點(以下為我對該頁算繪影像的判讀值。橫座標以圖框左右邊界與四條內部刻度線 0.2/0.4/0.6/0.8 校準後量得;縱座標誤差約 ±0.05×10¹² Hz):

分支符號讀到的資料點(K/K_max, ν in 10¹² Hz)
LO空心(0, 4.98)、(0.10, 4.92)、(0.17, 5.02)、(0.23, 4.87)、(0.50, 4.73)、(0.73, 4.43)、(1.0, 4.30)
TO實心(0, 3.57)、(0.24, 3.53)、(0.50, 3.37)、(0.75, 3.20)、(1.0, 3.07)
LA空心(0.36, 1.35)、(0.50, 1.82)、(0.72, 2.45)、(1.0, 2.76)
TA實心(0.35, 0.97)、(0.47, 1.29)、(0.54, 1.46)、(0.61, 1.60)、(0.72, 1.86)、(0.87, 2.08)、(1.0, 2.15)

ω_T 與 ω_L 是怎麼定義的。這是圖說唯一在做定義的一句話:把 TO 與 LO 這兩條光學支外插回 K = 0,得到的兩個頻率分別叫做 ω_T 與 ω_L。從上表讀:ω_T ↔ ν ≈ 3.57×10¹² Hz,ω_L ↔ ν ≈ 5.0×10¹² Hz(換成角頻率要乘 2π;Kittel 用 ω 符號但縱軸畫的是 ν,這是書中的符號鬆散處,讀的時候要自己補上 2π)。

為什麼要「外插」而不是直接量。中子非彈性散射在 K → 0 量不到光學支:散射截面正比於某種結構因子,且 K → 0 時能量轉移大而動量轉移小,運動學上難以達到。所以實驗只能在有限 K 取點再外推。(此為我的補充說明,書中未解釋為何需要外插。)

最重要的觀察:ω_L > ω_T。KBr 的 LO 與 TO 在 K = 0 沒有併攏,而是分裂成 5.0 與 3.57,比值約 1.40。這與 Ge 的「LO and TO branches coincide at K = 0」形成直接對比——兩張圖並排放,就是為了讓你看到這個差別。下一個 block 處理它。

聲學/光學頻隙。在區邊界,橫波這一組:TA 頂到 2.15,TO 底到 3.07,中間空著;縱波這一組:LA 頂到 2.76,LO 底到 4.30。兩組都有明顯頻隙,與 Fig. 7 的定性圖一致。

Ge 的 LO–TO 簡併 vs KBr 的 ω_L > ω_T:這是本段唯一沒被解釋的實驗事實我的判讀

把兩張圖的 K = 0 行為並排:

Ge(Fig. 8a)KBr(Fig. 8b)
K = 0 的光學支頻率LO = TO ≈ 9.15×10¹² Hz(簡併)ω_L ≈ 5.0、ω_T ≈ 3.57×10¹² Hz(分裂)
兩個基底原子都是 Ge,電中性K⁺ 與 Br⁻,帶相反電荷
Kittel 給的理由「a consequence of the crystal symmetry of Ge」(本段沒有解釋)

Kittel 在這一段只說了一半。他在 Fig. 8a 圖說裡把 Ge 的簡併歸給對稱性,但完全沒說 KBr 為什麼不簡併。這不是疏漏,而是安排:理由要到第 14 章才給。

理由(第 14 章 p.414,本段的前向連結)。第 14 章 Fig. 14 的圖說逐字說:「In the longitudinal optical phonon mode the particle displacements and hence the dielectric polarization P are parallel to the wavevector.」——LO 模的離子位移(因而極化 P)平行於 K,這使得 div E ≠ 0,晶體內部出現一個宏觀縱向電場;而 TO 模的位移垂直於 K,「It follows that div E = 0 for a TO phonon」,沒有宏觀電場。這個額外的靜電回復力只作用在 LO 上,於是把 LO 的頻率推高,ω_L > ω_T。

這個分裂有多重要。第 14 章 p.414 的 (62) 就是 Lyddane–Sachs–Teller(LST)關係:ω_L²/ω_T² = ε(0)/ε(∞),而 Kittel 特別註明「The derivation assumed a cubic crystal with two atoms per primitive cell」——正是本段的模型。所以本段 Fig. 8b 讀出來的兩個數字,可以直接換算成介電常數比:

(ω_L/ω_T)² = (5.0/3.57)² ≈ 1.96

也就是說,光憑一張中子散射的色散圖,就能預測 KBr 的靜態介電常數約是高頻介電常數的兩倍。(這個換算是我用本段讀到的數值與 p.414 的 (62) 做的,書中沒有對 KBr 做這個計算。)

回頭看 Ge。Ge 的兩個基底原子都是中性的 Ge 原子,光學模的位移不產生極化,沒有宏觀電場,LO 與 TO 就沒有理由分開——這與 Kittel 說的「crystal symmetry」是同一件事的兩種說法。第 14 章 p.414 的 LST 關係在此退化成 ω_L = ω_T ⇒ ε(0) = ε(∞),正確地反映了同極晶體沒有離子極化貢獻。

這也預告了本段最後 (26) 的物理意義:光學模之所以叫「光學」,就是因為在離子晶體裡它是一個振盪偶極;而 LO/TO 分裂正是這個偶極的自場效應。兩件事是同一個物理。

模型設定(p.97):交錯平面、對稱方向、最近鄰近似

從 p.97 開始,Kittel 把上面的計數論證換成一個真的可以解的模型。原文逐字:

「We consider a cubic crystal where atoms of mass M₁ lie on one set of planes and atoms of mass M₂ lie on planes interleaved between those of the first set (Fig. 9). It is not essential that the masses be different, but either the force constants or the masses will be different if the two atoms of the basis are in nonequivalent sites.」

第一句是幾何設定:質量 M₁ 的原子住在一組平面上,質量 M₂ 的原子住在穿插其間(interleaved)的另一組平面上。這是把三維問題再一次壓成一維的老招數(p.91 用過一次):只要波沿著某個方向傳播時整個平面同相運動,每個平面就只需要一個座標。

第二句是一個很重要的免責聲明:質量不同不是必要的。真正需要的是「兩個基底原子處在不等價位置時,力常數或質量其中之一會不同」。換句話說,分裂成兩支的根源是基底有兩個位置,不是「兩個原子重量不一樣」。這一點在討論 Ge 時會變成關鍵(Ge 的兩個原子質量完全相同,卻仍然有光學支)。

第二段是方向的選擇:

「Let a denote the repeat distance of the lattice in the direction normal to the lattice planes considered. We treat waves that propagate in a symmetry direction such that a single plane contains only a single type of ion; such directions are [111] in the NaCl structure and [100] in the CsCl structure.」

「一個平面只含一種離子」是模型成立的前提。如果平面裡混著兩種原子,就沒辦法用單一座標 u_s 描述整個平面。這個條件對方向有很強的限制,Kittel 給了兩個例子:NaCl 結構要走 [111],CsCl 結構要走 [100]。這也解釋了為什麼 Fig. 8a、8b 都是 [111] 方向——KBr 是 NaCl 結構,Ge 是金剛石結構(也是沿 [111] 分成純 Ge 的交錯平面)。

第三段是動力學近似:

「We write the equations of motion under the assumption that each plane interacts only with its nearest-neighbor planes and that the force constants are identical between all pairs of nearest-neighbor planes.」

兩個假設,兩個都會在後面付出代價:

  1. 只與最近鄰平面作用。p.94 明確警告過這在金屬裡不成立:「In metals the effective forces may be of quite long range... Interactions have been found between planes of atoms separated by as many as 20 planes.」
  2. 所有最近鄰平面對的力常數都相同(都是同一個 C)。這是本模型的單參數化,也是它最脆弱的地方——對 KBr 的 [111] 成立得相當好,對 Ge 的 [111] 則根本不成立。這件事後面有專門的 block。

a 的定義:「nearest identical planes」——雙原子鏈最大的陷阱

Kittel 為了這件事特地在 (19) 之後單獨留了一句話,這在他的寫作風格裡等於用紅筆圈起來。原文逐字:

「We define a in Fig. 9 as the distance between nearest identical planes, not nearest-neighbor planes.」

a 是「最近的同種平面」之間的距離,不是「最近鄰平面」之間的距離。在等間距的雙原子鏈裡,這兩個量差一倍:最近鄰平面間距是 a/2,而 a 是一個完整的重複週期(一個 M₁ 面 + 一個 M₂ 面)。

為什麼一定要這樣定義。因為 K 是晶格平移對稱性下的量子數,而晶格的平移週期是「一個原胞」,也就是 a。試解 (19) 寫的是 u_s = u·exp(isKa):這裡的 s 標的是原胞編號,不是平面編號。第 s 個原胞裡有一個 M₁(位移 u_s)和一個 M₂(位移 v_s),兩者共用同一個相位因子 exp(isKa),振幅差別由 u 與 v 承擔。如果你把 a 誤取成最近鄰平面間距,那麼 u_s 與 v_s 的相位就不該相同,整個試解結構就垮了。

它決定了布里淵區有多大。p.98 說得很明白:「The extent of the first Brillouin zone is −π/a ≤ K ≤ π/a, where a is the repeat distance of the lattice.」——區邊界在 π/a,不是 2π/a。雙原子鏈的布里淵區恰好是同樣原子密度的單原子鏈的一半寬。這件事會在「等質量極限」的 block 裡被用來做嚴格驗算。

四個最常見的算錯,全部源自這裡:

錯誤症狀正確做法
把 a 取成最近鄰間距區邊界算成 2π/a,聲速 (24) 差 2 倍,ω² 差 4 倍a = 同種平面間距 = 2 × 最近鄰間距(等間距時)
用傳統立方晶胞邊長當 a在 [111] 方向差一個 √3a = 該方向上的重複距離,NaCl [111] 是 a_cube/√3
把 a 當成最近鄰鍵長鍵長是三維距離,平面間距是它在 K 方向的投影兩者一般不等(KBr 的 K–Br 鍵長 3.298 Å,[111] 面間距只有 1.902 Å)
以為 exp(iKa) 裡的 a 和 (25) 裡的 a 是不同的 a前後不一致從 (18) 到 (26) 全部都是同一個 a

具體數字(我依書中他處數據計算,非本段原文)。KBr 是 NaCl 結構,書本 p.14 給 a_cube = 6.59 Å。沿 [111],K 面與 Br 面等間距交替,間距 a_cube/(2√3) = 1.902 Å,而同種平面的重複距離 a = a_cube/√3 = 3.805 Å。另可核對書本 p.66 Table 7 的 KBr 最近鄰距離 R₀ = 3.298 Å,這是三維的 K–Br 鍵長,與 1.902 Å 的面間距是完全不同的量——這正是上表第三列要防的錯。

鍺是金剛石結構,書本 p.17 給 a_cube = 5.658 Å,沿 [111] 的同種平面重複距離同樣是 a_cube/√3 = 3.266 Å。這與 Fig. 8a 圖說的 K_max = (2π/a_cube)(½ ½ ½) 一致:該向量的長度是 (2π/a_cube)(√3/2) = π√3/a_cube = π/(a_cube/√3) = π/3.266 Å⁻¹ ✓ 恰好是 π/a(本段意義下的 a)。這個一致性本身就是「a 必須是同種平面間距」的證明。

Fig. 9 逐筆讀圖:雙原子鏈的編號約定

網頁不放圖,以下用文字重建。圖說逐字是:「A diatomic crystal structure with masses M₁, M₂ connected by force constant C between adjacent planes. The displacements of atoms M₁ are denoted by u_{s−1}, u_s, u_{s+1},..., and of atoms M₂ by v_{s−1}, v_s, v_{s+1}. The repeat distance is a in the direction of the wavevector K. The atoms are shown in their undisplaced positions.」

圖面內容(以下為我對該頁算繪影像的目視描述):

圖分上下兩排,畫的是同一條鏈,上排帶位移箭頭、下排不帶。

  • 沿著 K 的方向(向右),原子交替排列,大的灰色圓是 M₁,小的黑點是 M₂。
  • 順序由左至右是:M₁(u_{s−1})、M₂(v_{s−1})、M₁(u_s)、M₂(v_s)、M₁(u_{s+1})、M₂(v_{s+1})。
  • 上排每個原子都有一個向右的水平箭頭,長度相同;符號 u_{s−1}、v_{s−1}、u_s、v_s、u_{s+1}、v_{s+1} 標在各自箭頭的上方。
  • 右側另有一個較粗的箭頭,標 K,指向右方——這是波向量的方向。
  • 下排的 M₁、M₂ 標在右邊兩個原子的上方,用來指認哪個符號對應哪種球(圖的下方是標 a 的尺寸線)。
  • 最關鍵的是圖底的距離標示:一條雙箭頭線段標著 a,它的兩端落在兩個相鄰的小黑點(M₂ 面)上,也就是 v_s 與 v_{s+1} 的位置。這就是圖說所說「repeat distance is a」的視覺定義——a 跨過一個 M₁ 和一個 M₂,是兩個同種平面之間的距離,與 p.97 那句 「nearest identical planes, not nearest-neighbor planes」完全一致。

這張圖唯一要記住的事:編號的偏心。第 s 個原胞包含 u_s(左)與 v_s(右)。因此:

  • M₁(u_s) 的兩個鄰居是 v_s(右邊,同一個原胞)與 v_{s−1}(左邊,前一個原胞)。
  • M₂(v_s) 的兩個鄰居是 u_s(左邊,同一個原胞)與 u_{s+1}(右邊,下一個原胞)。

這個不對稱純粹來自編號約定,不是物理上的不對稱。它直接決定了 (18) 兩條方程長得不一樣,也決定了 (20) 裡一條帶 exp(−iKa)、另一條帶 exp(+iKa)。如果你把編號改成「v_s 在 u_s 左邊」,兩個指數就會交換,但最終的久期方程 (22) 完全不變——因為 (22) 只含 cos Ka。

雙原子鏈的運動方程

M₁ d²u_s/dt² = C(v_s + v_{s−1} − 2u_s) ; M₂ d²v_s/dt² = C(u_{s+1} + u_s − 2v_s) .

(18)  p.97

u_s:第 s 個原胞中 M₁ 平面的位移;v_s:第 s 個原胞中 M₂ 平面的位移;M₁、M₂:兩種平面上原子的質量;C:相鄰平面之間的力常數(依 p.91 的約定,C 是對平面中「一個原子」定義的,故 M 也是單一原子質量,量綱一致)。

出處推導
由 p.91 的 (1)「F_s = C(u_{s+1} − u_s) + C(u_{s−1} − u_s)」推廣而來,把「所有平面等價」換成「兩種平面交替」,並依 Fig. 9 的鄰居關係指定每一項。使用的假設有三:位移線性(彈性能為二次式)、只與最近鄰平面作用、所有最近鄰平面對共用同一個 C。
極限行為
M₁ = M₂ = M 時,兩式形式相同,系統退化回等間距的單原子鏈 (2),但要注意此時的「原胞」仍含兩個平面,所以布里淵區仍是 π/a——這是 zone folding 的情形,後面有專門推導。C → 0 時兩式解耦,每個平面獨立、ω = 0(自由粒子,無回復力)。M₂ → ∞ 時第二式給 d²v_s/dt² → 0,重原子面不動,M₁ 面在固定的兩端之間振盪,ω² → 2C/M₁ 與 K 無關——這就是 Einstein 振子的極限,光學支變平。
典型數值
由 Fig. 8b 的 KBr 資料擬合(見後面的數值 block):橫波這一組 C ≈ 11 N/m、縱波這一組 C ≈ 21.5 N/m。作為對照,典型的化學鍵力常數在 10–500 N/m 量級(共價鍵在高端,離子與凡得瓦鍵在低端),KBr 的 10 N/m 量級與離子鍵的軟度相符。

(18) 的受力分析:為什麼兩式的下標不對稱補推導

Kittel 只寫「We refer to Fig. 9 to obtain (18)」,把受力分析留給讀者。這一步其實只有三行,但下標很容易寫錯,值得逐項核對。回憶 p.91 的 (1):一個平面受到的力,是它與每個鄰居的相對位移乘上力常數的總和。

  1. 第 1 步:先寫通式。任一平面受到的力 = C × Σ(鄰居位移 − 自己的位移),對兩個最近鄰平面求和。這是 (1) 的形式,也是 Hooke 定律。
  2. 第 2 步:M₁(u_s) 的鄰居。從 Fig. 9 讀:u_s 的右鄰是 v_s(同一原胞內的 M₂),左鄰是 v_{s−1}(前一原胞的 M₂)。所以

    F(u_s) = C(v_s − u_s) + C(v_{s−1} − u_s) = C(v_s + v_{s−1} − 2u_s)

    代入牛頓第二定律 ⇒ M₁ d²u_s/dt² = C(v_s + v_{s−1} − 2u_s) ✓ 與 (18) 第一式相符。

  3. 第 3 步:M₂(v_s) 的鄰居。同樣從 Fig. 9 讀:v_s 的左鄰是 u_s(同一原胞內的 M₁),右鄰是 u_{s+1}(下一原胞的 M₁)。所以

    F(v_s) = C(u_{s+1} − v_s) + C(u_s − v_s) = C(u_{s+1} + u_s − 2v_s)

    ⇒ M₂ d²v_s/dt² = C(u_{s+1} + u_s − 2v_s) ✓ 與 (18) 第二式相符。

  4. 第 4 步:不對稱的來源。第一式用 {s, s−1},第二式用 {s+1, s}。差別只來自「在第 s 個原胞裡,u 在左、v 在右」這個編號約定。若把約定反過來,兩式的下標就對調,(20) 裡的 exp(±iKa) 也對調——但 (22) 完全不受影響。
  5. 第 5 步:三個健全性檢查。 (a) 平移不變性:令所有 u_s = v_s = d(常數),則 (18) 右側 = C(d + d − 2d) = 0,兩式都給零加速度 ✓。整體平移不產生回復力,這保證了 K → 0 時聲學支的 ω → 0。 (b) 動量守恆:把 (18) 的兩式對所有 s 求和,每一項 C·u 或 C·v 都出現兩次、一正一負,總和為零 ⇒ d²/dt²(Σ M₁u_s + Σ M₂v_s) = 0,總動量守恆 ✓。 (c) 量綱:[C] = N/m,[C·u] = N ✓ 與 [M d²u/dt²] = N 一致。

為什麼要這麼小心。這兩條方程的下標若寫錯一個,(20) 的兩個指數就不再是共軛對,(21) 的行列式展開就不會產生實數的 cos Ka,最後 (22) 會變成複數方程,兩個根也就不再都是實的正數——整個物理圖像會崩掉。下標的一致性不是形式問題,它是 ω² 為實數的保證。

順帶一提:這個「共軛對」結構其實是 Hermiticity 的表現。把 (20) 寫成矩陣方程 ω²Mx = Dx,動力學矩陣 D 的非對角元是 −C[1 + exp(∓iKa)],兩者互為複共軛 ⇒ D 是 Hermitian ⇒ 本徵值 ω² 必為實數。這個結構在第 4 章之後的所有晶格動力學裡都成立,不限於這個玩具模型。(Hermitian 這個說法是我的補充,書中此處未提。)

行波試解(不同振幅的兩個子晶格)

u_s = u·exp(isKa)·exp(−iωt) ; v_s = v·exp(isKa)·exp(−iωt) .

(19)  p.97

u、v:兩個子晶格的複振幅(常數,與 s 無關);s:原胞編號(整數);a:同種平面之間的重複距離;K:波向量;ω:角頻率。exp(isKa) 是空間相位,exp(−iωt) 是時間相位。

出處推導
由 p.92 的 (4)「u_{s±1} = u·exp(isKa)·exp(±iKa)」推廣。原文引言逐字:「We look for a solution in the form of a traveling wave, now with different amplitudes u, v on alternate planes」——關鍵字是 now with different amplitudes:相對於單原子鏈的唯一新東西,就是兩個子晶格可以有不同的振幅(也可以有不同的相位,因為 u、v 是複數)。這個多出來的自由度正是第二條分支的來源。
極限行為
u = v 時(同相同幅)兩個子晶格一起動,這是 K → 0 的聲學支(p.99:「Another solution for the amplitude ratio at small K is u = v」)。u/v 為負實數時兩個子晶格反相,這是 K = 0 的光學支,比值由 (26) 給出。一般的 K,p.99 明說:「At a general K the ratio u/v will be complex」——u 與 v 之間有一個既非 0 也非 π 的相位差。這正是為什麼 u、v 必須允許是複數,而不能只寫成實振幅。

校註:(19) 印成 v₂,應為 v_s

這一條是讀原書時會踩到的印刷錯誤,必須指出來,否則後面的代數對不上。

書本 p.97 的 (19) 第二式印的是 v₂ = v·exp(isKa)·exp(−iωt)——下標是阿拉伯數字 2,不是斜體的 s。我把該頁放大到 10 倍算繪目檢確認過:左式的 u_s 下標是斜體的 s,右式的 v₂ 下標是直立的數字 2,兩者字形明顯不同。這是書上的排版錯誤,不是 PDF 文字層的抽取失誤。

正確的應該是 v_s。三個理由:

  1. (18) 的第二式左側是 M₂ d²v_s/dt²,方程要能代入就必須有 v_s 的表達式,v₂ 在這裡沒有任何意義。
  2. Fig. 9 的圖說明確定義了「the displacements... of atoms M₂ by v_{s−1}, v_s, v_{s+1}」,用的是 s 標號。
  3. 若下標真的是固定的 2,則 v 與 s 無關,(20) 就無法從 (18) 導出。

讀這一段時的實際影響。如果你照書抄成 v₂,代入 (18) 時會發現 v_{s−1} 與 v_s 都變成同一個 v₂,於是第一式右側變成 C(2v₂ − 2u_s),不會產生 exp(−iKa) 這一項,最後拿不到 (21) 的非對角元。所以只要你真的動手算,這個錯誤會自己暴露出來——但若只是看,很容易滑過去。

(19) 代入 (18):逐步得到 (20)補推導

Kittel 只寫「On substitution of (19) in (18) we have (20)」,中間三行沒有。這三行不難,但它示範了一個關鍵技巧:指數因子 exp(isKa)·exp(−iωt) 會在每一項出現,可以整體約掉——這正是 Bloch 形式試解的全部威力所在,它把無窮多條耦合的常微分方程化成兩條代數方程。

  1. 第 1 步:處理時間微分。由 (19),u_s 對 t 的依賴只有 exp(−iωt),所以

    d²u_s/dt² = (−iω)²·u_s = −ω²u_s,同理 d²v_s/dt² = −ω²v_s。

    (這一步與單原子鏈 p.91 的 (3) 完全相同:「We look for solutions with all displacements having the time dependence exp(−iωt). Then d²u_s/dt² = −ω²u_s」。)

  2. 第 2 步:處理空間移位。由 (19),

    v_{s−1} = v·exp[i(s−1)Ka]·exp(−iωt) = v·exp(isKa)·exp(−iωt)·exp(−iKa)

    u_{s+1} = u·exp[i(s+1)Ka]·exp(−iωt) = u·exp(isKa)·exp(−iωt)·exp(+iKa)

    移位一格 = 乘上一個相位因子 exp(±iKa)。

  3. 第 3 步:代入 (18) 第一式。

    M₁(−ω²)·u·exp(isKa)e^{−iωt} = C[ v·exp(isKa)e^{−iωt} + v·exp(isKa)e^{−iωt}·exp(−iKa) − 2u·exp(isKa)e^{−iωt} ]

  4. 第 4 步:約掉共同因子 exp(isKa)·exp(−iωt)。這個因子在每一項都出現,且永不為零,可以整體消去:

    −ω²M₁u = C[ v + v·exp(−iKa) − 2u ] = Cv[1 + exp(−iKa)] − 2Cu ✓ 這是 (20) 第一式。

  5. 第 5 步:對 (18) 第二式做同樣的事。

    M₂(−ω²)·v·exp(isKa)e^{−iωt} = C[ u·exp(isKa)e^{−iωt}·exp(iKa) + u·exp(isKa)e^{−iωt} − 2v·exp(isKa)e^{−iωt} ]

    約掉共同因子:

    −ω²M₂v = C[ u·exp(iKa) + u − 2v ] = Cu[exp(iKa) + 1] − 2Cv ✓ 這是 (20) 第二式。

  6. 第 6 步:檢查共軛結構。第一式的耦合係數是 C[1 + exp(−iKa)],第二式是 C[exp(iKa) + 1]——兩者互為複共軛。這不是巧合,它是第 4 步約掉的相位因子在兩個子晶格之間相差半個原胞所留下的印記,也是下一步 (21)→(22) 能產生實數 cos Ka 的唯一原因。

這一步在做什麼(物理意義)。(18) 是無窮多條耦合方程(每個 s 一條,共 2N 條)。Bloch 形式的試解 (19) 假設了「相鄰原胞之間只差一個固定相位 exp(iKa)」,這把 2N 條方程一次壓成 2 條(每個 K 值各一組 2×2 系統)。付出的代價是:原本一個解要指定 2N 個振幅,現在要對 N 個不同的 K 各解一次 2×2 系統——總數仍是 2N,但每一組只有 2 維,可以手算。

這就是為什麼晶體比玻璃好算。平移對稱性把一個 2N 維的本徵值問題分塊對角化成 N 個 2 維問題。玻璃沒有平移對稱性,K 不是好量子數,這個分塊就不存在——這是非晶固體聲子計算困難的根本原因。(此註為我的補充。)

代入試解後的 2×2 代數方程

−ω²M₁u = Cv[1 + exp(−iKa)] − 2Cu ; −ω²M₂v = Cu[exp(iKa) + 1] − 2Cv .

(20)  p.97

u、v:兩個子晶格的複振幅;其餘符號同 (18)、(19)。exp(−iKa) 與 exp(+iKa) 分別來自 v_{s−1} 與 u_{s+1} 的相位移位。

出處推導
把 (19) 代入 (18),利用 d²/dt² → −ω² 與移位 → exp(±iKa),再把每一項共有的因子 exp(isKa)exp(−iωt) 約掉。完整五步見前一個補推導 block。
極限行為
K = 0:兩個指數都變成 1,方程化為 −ω²M₁u = 2Cv − 2Cu 與 −ω²M₂v = 2Cu − 2Cv。此時要嘛 u = v 且 ω = 0(聲學支),要嘛 u/v = −M₂/M₁ 且 ω² = 2C(1/M₁+1/M₂)(光學支,即 (26))。Ka = π(區邊界):exp(±iπ) = −1,兩個耦合係數都變成 C(1 − 1) = 0——兩條方程完全解耦!這是區邊界最重要的結構特徵:M₁ 面與 M₂ 面不再互相驅動,各自獨立振盪,直接給出 (25) 的兩個根 2C/M₁ 與 2C/M₂。理解了這一點,(25) 根本不需要計算就能寫出來。

久期行列式(secular determinant)

| 2C − M₁ω² −C[1 + exp(−iKa)] | | −C[1 + exp(iKa)] 2C − M₂ω² | = 0

(21)  p.98

矩陣的兩列分別對應 (20) 的兩條方程,兩行分別對應未知數 u 與 v。對角元 2C − M_iω² 是「自身」項,非對角元 −C[1 + exp(∓iKa)] 是兩個子晶格之間的耦合。

出處推導
把 (20) 移項成齊次形式:(2C − M₁ω²)u − C[1 + exp(−iKa)]v = 0 與 −C[1 + exp(iKa)]u + (2C − M₂ω²)v = 0。引言逐字是「The homogeneous linear equations have a solution only if the determinant of the coefficients of the unknowns u, v vanishes」——齊次線性方程組要有非平庸解(u、v 不全為零),係數矩陣必須奇異,行列式為零。這是 Cramer 法則的直接推論。(校註:書上此處兩個非對角元皆印為 −C[1 + exp(+iKa)];本欄的右上角已依 (20) 修正為 −C[1 + exp(−iKa)],理由詳見下一個 block。)
極限行為
這是一個以 ω² 為未知數的本徵值問題。矩陣可以寫成 D(K) − ω²M 的形式(M = diag(M₁, M₂)),ω² 是廣義本徵值。由於兩個非對角元互為複共軛,D(K) 是 Hermitian 矩陣,保證兩個 ω² 根都是實數;又因為 D 是半正定(它來自一個非負的彈性能二次式),兩個根都非負,ω 才有物理意義。Ka = π 時非對角元為零,矩陣對角化,兩個根直接讀出 2C/M₁ 與 2C/M₂。C = 0 時矩陣變成 −diag(M₁ω², M₂ω²),只有 ω = 0 的解。

校註:(21) 兩個非對角元在書上都印成 exp(+iKa)

這是本段第二個、也是更嚴重的一個印刷錯誤,我把 p.98 該處放大算繪目檢確認過。

書本 p.98 的 (21) 印出來是:

2C − M₁ω² −C[1 + exp(iKa)]
−C[1 + exp(iKa)] 2C − M₂ω²= 0

兩個非對角元都印成 exp(+iKa)。這與前一頁(p.97)的 (20) 不一致:(20) 第一式的耦合係數是 C[1 + exp(−iKa)],第二式才是 C[exp(+iKa) + 1]。所以 (21) 的右上角應該是 −C[1 + exp(−iKa)]。

為什麼可以確定是印錯,而不是我理解錯。因為書上緊接著給的 (22) 只有在「共軛對」的情況下才成立:

  • 若照書上印的(兩個都是 +iKa),行列式展開會得到 −C²[1 + exp(iKa)]² = −C²[1 + 2exp(iKa) + exp(2iKa)],這是一個複數,不可能化成 (22) 裡的實數項 2C²(1 − cos Ka)。
  • 若照 (20) 應有的共軛對,非對角元乘積是 C²[1 + exp(iKa)][1 + exp(−iKa)] = C²|1 + exp(iKa)|² = 2C²(1 + cos Ka),實數,展開後恰好給出 (22)。

所以 (22) 本身就是 (21) 有印刷錯誤的證據。下一個 block 會把這個展開完整做一遍。

對讀者的實際影響。如果你照著書上的 (21) 硬算,會得到一個複係數的四次方程,兩個「ω²」根會是複數,物理上無法解釋,而且怎麼樣都湊不出 (22)。這是一個典型的「書上的錯誤讓你以為自己算錯」的情境——這也正是要把每一步補出來的理由。

(21) → (22):展開 2×2 行列式(書上直接跳過的一步)補推導

書上在 (21) 與 (22) 之間只寫了一個「or」。這一步是本段代數的核心,而且它用到一個小恆等式 |1 + exp(iKa)|² = 2(1 + cos Ka),值得單獨記下來。以下用修正過的 (21)(右上角為 exp(−iKa),理由見前一個校註 block)。

  1. 第 1 步:寫下 2×2 行列式的定義。對矩陣 [[A, B],[C', D]],行列式為 AD − BC'。這裡

    A = 2C − M₁ω², D = 2C − M₂ω²

    B = −C[1 + exp(−iKa)], C' = −C[1 + exp(iKa)]

    所以條件是 (2C − M₁ω²)(2C − M₂ω²) − C²[1 + exp(−iKa)][1 + exp(iKa)] = 0(兩個負號相乘變正)。

  2. 第 2 步:展開對角乘積。

    (2C − M₁ω²)(2C − M₂ω²) = 4C² − 2C·M₂ω² − 2C·M₁ω² + M₁M₂ω⁴

    = M₁M₂ω⁴ − 2C(M₁ + M₂)ω² + 4C²

    第一項與第二項已經與 (22) 的前兩項一模一樣了。

  3. 第 3 步:展開非對角乘積——關鍵的小恆等式。

    [1 + exp(iKa)][1 + exp(−iKa)]

    = 1 + exp(−iKa) + exp(iKa) + exp(iKa)exp(−iKa)

    = 1 + [exp(iKa) + exp(−iKa)] + 1

    = 2 + 2cos Ka

    = 2(1 + cos Ka)

    用到 p.92 (6)→(7) 那裡同一條恆等式:2cos Ka = exp(iKa) + exp(−iKa)。

    另一種看法:這兩個因子互為複共軛,所以乘積就是模方 |1 + exp(iKa)|²。用歐拉式 1 + e^{iθ} = 2cos(θ/2)·e^{iθ/2} 可得模方 = 4cos²(Ka/2) = 2(1 + cos Ka) ✓ 兩條路同結果。

  4. 第 4 步:合併。把第 2、3 步代回第 1 步:

    M₁M₂ω⁴ − 2C(M₁ + M₂)ω² + 4C² − C²·2(1 + cos Ka) = 0

    展開常數項:4C² − 2C² − 2C²cos Ka = 2C² − 2C²cos Ka = 2C²(1 − cos Ka)

  5. 第 5 步:得到 (22)。

    M₁M₂ω⁴ − 2C(M₁ + M₂)ω² + 2C²(1 − cos Ka) = 0 ✓ 與書上的 (22) 完全一致。

  6. 第 6 步:健全性檢查。 (a) K = 0:1 − cos 0 = 0,常數項消失 ⇒ ω²[M₁M₂ω² − 2C(M₁+M₂)] = 0,兩根為 ω² = 0(聲學)與 ω² = 2C(M₁+M₂)/(M₁M₂) = 2C(1/M₁ + 1/M₂)(光學,即 (23))。 (b) 量綱:[M²ω⁴] = kg²·s⁻⁴;[C·M·ω²] = (N/m)(kg)(s⁻²) = (kg s⁻²)(kg)(s⁻²) = kg²s⁻⁴ ✓;[C²] = kg²s⁻⁴ ✓ 三項同量綱。 (c) 週期性:(22) 只透過 cos Ka 依賴 K ⇒ ω(K) 以 2π/a 為週期,且為 K 的偶函數 ⇒ 第一布里淵區 −π/a < K ≤ π/a 已包含全部資訊 ✓ 與 p.93 的論證一致。

這一步的教訓。(22) 只含 cos Ka,所有複數在第 3 步都消掉了。這件事不是運氣:它是 (20) 那兩個耦合係數互為共軛的直接後果,而共軛又來自 Fig. 9 的編號約定(一個鄰居在 s−1、一個在 s+1)。整條因果鏈是:

晶格平移對稱 ⇒ Bloch 試解 ⇒ 動力學矩陣 Hermitian ⇒ 行列式為實 ⇒ ω² 為實。

如果你在某一步算出複數的 ω²,回頭查 (20) 的兩個指數是不是共軛對——這是最常見的出錯點,也是書上 (21) 印錯的地方。

久期方程(雙原子鏈的色散關係,隱式)

M₁M₂ω⁴ − 2C(M₁ + M₂)ω² + 2C²(1 − cos Ka) = 0

(22)  p.98

這是一個以 ω² 為未知數的二次方程(對 ω 而言是四次,但只含偶次冪)。M₁M₂ 是 ω⁴ 的係數,2C(M₁+M₂) 是 ω² 的係數,2C²(1 − cos Ka) 是常數項,也是唯一含 K 的地方。

出處推導
把 (21) 的 2×2 行列式展開,並用 [1 + exp(iKa)][1 + exp(−iKa)] = 2(1 + cos Ka)。書上從 (21) 直接跳到 (22),中間六步見前一個補推導 block。(注意:展開必須使用共軛對形式的 (21);書上 (21) 兩個非對角元都印成 exp(+iKa),是排版錯誤。)
極限行為
Vieta 關係給了兩個免費的檢查工具:兩根之和 ω₊² + ω₋² = 2C(M₁+M₂)/(M₁M₂) = 2C(1/M₁ + 1/M₂),與 K 完全無關——這意味著光學支往下掉多少,聲學支就往上升多少,兩條曲線的高度和是常數。兩根之積 ω₊²ω₋² = 2C²(1 − cos Ka)/(M₁M₂),在 K → 0 時趨於零,強迫其中一根(聲學支)趨於零。K → 0:常數項 → 0,退化成 ω²[M₁M₂ω² − 2C(M₁+M₂)] = 0。Ka = π:1 − cos π = 2,常數項為 4C²,方程可因式分解為 (M₁ω² − 2C)(M₂ω² − 2C) = 0 ⇒ 直接得到 (25) 的兩根。C → 0:只剩 M₁M₂ω⁴ = 0 ⇒ ω = 0(雙重根),無回復力則無振動。
典型數值
以 KBr 橫波為例(我依 Fig. 8b 擬合,見後面的數值 block):C = 11.0 N/m、M₁ = M(Br) = 79.90 u、M₂ = M(K) = 39.10 u。Vieta 和:ω₊² + ω₋² = 2(11.0)(1/1.327×10⁻²⁵ + 1/6.492×10⁻²⁶) = 22.0 × (7.536×10²⁴ + 1.5404×10²⁵) = 22.0 × 2.294×10²⁵ = 5.047×10²⁶ s⁻²,對應「兩支頻率平方之和」恆為 5.05×10²⁶ s⁻²——可以拿 Fig. 8b 的任一個 K 值去核對這條守恆律。

(22) 的精確解:書上說「可以精確解」但沒寫出來補推導

Kittel 在 (22) 之後寫:「We can solve this equation exactly for ω², but it is simpler to examine the limiting cases Ka ≪ 1 and Ka = ±π at the zone boundary.」——他明說可以精確解,但選擇不寫。精確解其實只有兩行,而且寫出來之後,(23)、(24)、(25) 全部可以當成它的極限一次驗完,比分開算還省事。

  1. 第 1 步:套二次公式。把 (22) 看成 Aω⁴ + Bω² + D = 0,其中 A = M₁M₂,B = −2C(M₁+M₂),D = 2C²(1 − cos Ka)。

    ω² = [2C(M₁+M₂) ± √( 4C²(M₁+M₂)² − 8C²M₁M₂(1 − cos Ka) )] / (2M₁M₂)

  2. 第 2 步:把 2C(M₁+M₂) 提出根號。根號內 = 4C²(M₁+M₂)²·[1 − 2M₁M₂(1 − cos Ka)/(M₁+M₂)²],開根號得 2C(M₁+M₂)·√[...]。於是

    ω²_± = [C(M₁+M₂)/(M₁M₂)] · { 1 ± √[ 1 − 2M₁M₂(1 − cos Ka)/(M₁+M₂)² ] }

  3. 第 3 步:引入約化質量把式子變乾淨。定義 μ = M₁M₂/(M₁+M₂)(約化質量)。則 C(M₁+M₂)/(M₁M₂) = C/μ,且 2M₁M₂/(M₁+M₂)² = 2μ/(M₁+M₂)。精確解變成

    ω²_± = (C/μ)·{ 1 ± √[ 1 − (2μ/(M₁+M₂))(1 − cos Ka) ] }

    上號(+)是光學支,下號(−)是聲學支。

  4. 第 4 步:驗 K = 0。1 − cos 0 = 0 ⇒ 根號 = 1 ⇒ ω²₊ = 2C/μ = 2C(1/M₁ + 1/M₂) ✓ 就是 (23);ω²₋ = 0 ✓ 聲學支從原點出發。
  5. 第 5 步:驗 Ka = π。1 − cos π = 2 ⇒ 根號內 = 1 − 4μ/(M₁+M₂) = 1 − 4M₁M₂/(M₁+M₂)² = [(M₁+M₂)² − 4M₁M₂]/(M₁+M₂)² = (M₁−M₂)²/(M₁+M₂)²

    根號 = |M₁ − M₂|/(M₁+M₂)。設 M₁ > M₂(與 Fig. 7 一致):

    ω²₊ = (C/μ)·[(M₁+M₂) + (M₁−M₂)]/(M₁+M₂) = C·(M₁+M₂)/(M₁M₂) · 2M₁/(M₁+M₂) = 2C/M₂

    ω²₋ = (C/μ)·[(M₁+M₂) − (M₁−M₂)]/(M₁+M₂) = C·(M₁+M₂)/(M₁M₂) · 2M₂/(M₁+M₂) = 2C/M₁

    兩者正是 (25) 的兩個值 ✓,而且順帶回答了書上沒說的問題:哪一根屬於哪一支。

  6. 第 6 步:判別式恆為正。根號內 = 1 − (2μ/(M₁+M₂))(1 − cos Ka)。因為 1 − cos Ka ≤ 2,最壞情況下根號內 ≥ 1 − 4μ/(M₁+M₂) = (M₁−M₂)²/(M₁+M₂)² ≥ 0。

    ⇒ 兩根永遠是實數且非負,ω 永遠有物理意義 ✓。等號只在 M₁ = M₂ 且 Ka = π 時成立——此時兩根相等,頻隙關閉。

精確解揭露了三件在極限式裡看不見的事。

(一)頻隙的寬度由質量差決定。由第 5 步,區邊界的兩根是 2C/M₁ 與 2C/M₂,頻隙的存在與否完全取決於 M₁ ≠ M₂(在 C 相同的前提下)。質量越接近,頻隙越窄;質量相等,頻隙關閉。

(二)哪個根配哪一支,是由質量大小決定的。聲學支的頂端是 2C/M₁(重的那個),光學支的底端是 2C/M₂(輕的那個)。這與 Fig. 7 的標註完全吻合:圖上寫著 M₁ > M₂,聲學支頂端標 (2C/M₁)^{1/2}、光學支底端標 (2C/M₂)^{1/2},而且圖上聲學支的頂端確實比光學支的底端低。書上的 (25) 只並列寫出「ω² = 2C/M₁ ; ω² = 2C/M₂」,沒有說哪個屬於哪一支——要看 Fig. 7 才知道,或者像這裡一樣從精確解推出來。

(三)區邊界模的物理圖像。由 (20) 的極限分析:Ka = π 時耦合係數 C[1 + exp(±iπ)] = 0,兩個子晶格完全解耦。ω² = 2C/M₁ 的模是「M₁ 面振盪、M₂ 面完全靜止」,ω² = 2C/M₂ 則相反。重的那個振盪得慢(頻率低 ⇒ 聲學支頂端),輕的振盪得快(頻率高 ⇒ 光學支底端)——這就是為什麼 2C/M₁ 配聲學支,不需要記,用「重的慢」就能推出來。

小 Ka 極限:用主導平衡一行拿到 (23) 與 (24)補推導

書上寫:「For small Ka we have cos Ka ≅ 1 − ½K²a² + . . . , and the two roots are」——一句話帶過,緊接著就是兩條展開式,式號 (23)、(24) 分別排在它們的右緣(不是句子的一部分)。但這一步有一個真正的陷阱:兩個根不能用同一種方式展開。若你把 cos 的展開代進二次公式然後對根號做泰勒,聲學支的領導項會發生相消,粗心的話會得到 ω² = 0。以下給兩條路:一條是 Kittel 暗示的、一條是更快也更不容易錯的「主導平衡」(dominant balance)。

校註:書本原文此處印的是「cos Ka ≅ 1 − ½K²a² + . . .」(約等於符號),(23)、(24) 左側也都是 ≅。PDF 文字層把這些 ≅ 全部輸出成 ≡(恆等於),是字符解碼的失真;我已對該頁算繪影像目檢確認為 ≅。另外 (24) 分子的 ½ 在文字層被拆成浮出的「1 2」碎片,目檢確認分子是 ½C,整式為 ω² ≅ [½C/(M₁+M₂)]K²a²。

  1. 第 0 步:把展開式寫清楚。cos Ka ≅ 1 − ½K²a² + …,所以 1 − cos Ka ≅ ½K²a²。代入 (22):

    M₁M₂ω⁴ − 2C(M₁+M₂)ω² + 2C²·(½K²a²) = 0,即

    M₁M₂ω⁴ − 2C(M₁+M₂)ω² + C²K²a² = 0

  2. 第 1 步(光學支,主導平衡)。假設 ω² 在 K → 0 時趨於一個有限的非零值。那麼三項的大小是:ω⁴ 項 O(1)、ω² 項 O(1)、常數項 O(K²)。丟掉 O(K²) 的常數項:

    M₁M₂ω⁴ − 2C(M₁+M₂)ω² ≅ 0 ⇒ ω²[M₁M₂ω² − 2C(M₁+M₂)] ≅ 0

    非零解:ω² ≅ 2C(M₁+M₂)/(M₁M₂) = 2C(1/M₁ + 1/M₂) ✓ 這是 (23)。

  3. 第 2 步(聲學支,主導平衡)。假設 ω² 在 K → 0 時趨於零,且 ω² = O(K²)。那麼 ω⁴ 項是 O(K⁴),ω² 項是 O(K²),常數項是 O(K²)。丟掉 O(K⁴) 的 ω⁴ 項(這是關鍵:丟的是最高次項,不是最低次項):

    −2C(M₁+M₂)ω² + C²K²a² ≅ 0

    ⇒ ω² ≅ C²K²a² / [2C(M₁+M₂)] = [½C/(M₁+M₂)]·K²a² ✓ 這是 (24)。

  4. 第 3 步:驗證自洽。第 1 步假設 ω² = O(1),結果 ω² = 2C/μ,確實 O(1) ✓。第 2 步假設 ω² = O(K²),結果正比於 K² ✓。兩個假設都自洽,而且它們窮盡了二次方程的兩個根(因為 Vieta 的積 = C²K²a²/(M₁M₂) = O(K²),兩根的積是 O(K²),一根 O(1) 則另一根必為 O(K²))。
  5. 第 4 步(另一條路:從精確解展開)。由前一個 block 的精確解,令 x ≡ (2μ/(M₁+M₂))(1 − cos Ka) ≅ (2μ/(M₁+M₂))·½K²a² = μK²a²/(M₁+M₂),這是一個 O(K²) 的小量。

    光學支:ω²₊ = (C/μ)(1 + √(1−x)) ≅ (C/μ)(1 + 1 − x/2) ≅ 2C/μ = 2C(1/M₁+1/M₂) ✓

    聲學支:ω²₋ = (C/μ)(1 − √(1−x))。這裡 1 − √(1−x) 的領導項相消為零,必須展到一階:√(1−x) ≅ 1 − x/2 ⇒ 1 − √(1−x) ≅ x/2。所以

    ω²₋ ≅ (C/μ)·(x/2) = (C/μ)·μK²a²/(2(M₁+M₂)) = CK²a²/(2(M₁+M₂)) ✓ 與 (24) 相同。

    這一步就是陷阱所在:如果只把 √(1−x) 近似成 1,會得到 ω²₋ ≅ 0,整個聲學支消失。

  6. 第 5 步:Vieta 交叉驗算。(23) × (24) = [2C(M₁+M₂)/(M₁M₂)] × [CK²a²/(2(M₁+M₂))] = C²K²a²/(M₁M₂)。

    而 (22) 的兩根之積 = 常數項/首項 = 2C²(1 − cos Ka)/(M₁M₂) ≅ 2C²(½K²a²)/(M₁M₂) = C²K²a²/(M₁M₂) ✓ 完全相符。

    再驗和:(23) + (24) ≅ 2C(1/M₁+1/M₂) + O(K²),而 Vieta 的和 = 2C(M₁+M₂)/(M₁M₂) = 2C(1/M₁+1/M₂) ✓((24) 就是這個和的 O(K²) 修正被分給了聲學支——這也再次說明「兩支的頻率平方和守恆」:光學支在小 K 掉下去的量,正好等於聲學支升上來的量)。

為什麼要學「主導平衡」這個技巧。第 1、2 步各只有一行,而且不需要用到二次公式。它的邏輯是:先猜結果的量級,用這個猜測判斷哪一項可以丟,解出來之後再回頭驗證猜測自洽。這個方法在整本書後面一再出現(第 6 章自由電子氣的低溫展開、第 8 章半導體的載子統計、第 14 章介電函數的極點附近),值得在這裡練熟。

兩個結果的物理落差有多大。(23) 完全不含 K——光學支在 K → 0 附近是平的,群速度 dω/dK → 0,這種聲子不輸運能量。(24) 給 ω ∝ K,群速度是常數 = 聲速——聲學支才是輸運熱的主角。這個對比是第 5 章 Debye 模型只保留聲學支、以及熱傳導中「光學聲子幾乎不貢獻 κ」的全部根據。

光學支在長波極限(K → 0)的頻率

ω² ≅ 2C(1/M₁ + 1/M₂) = 2C/μ , μ ≡ M₁M₂/(M₁+M₂)

(23)  p.98

ω:光學支在 K → 0 的角頻率;C:最近鄰平面之間的力常數;M₁、M₂:兩種平面的原子質量;μ:約化質量(1/μ = 1/M₁ + 1/M₂,此符號為我引入,書中未用)。

出處推導
由 (22) 在 Ka ≪ 1 的極限取大根。書上只寫「For small Ka we have cos Ka ≅ 1 − ½K²a² + . . . , and the two roots are」就直接列出。補推導見前一個 block:對 (22) 丟掉 O(K²) 的常數項,得 ω²[M₁M₂ω² − 2C(M₁+M₂)] ≅ 0 的非零根。(書上此處印的是 ≅,PDF 文字層誤輸出為 ≡。)
極限行為
完全不含 K,這是光學支最重要的性質:在 K = 0 附近色散是平的,群速度 v_g = dω/dK → 0。物理後果是光學聲子幾乎不輸運熱——它們可以吸收能量(貢獻熱容),但不能把能量搬到別處。這是第 5 章 Debye 模型只用聲學支算熱傳導的正當性,也是為什麼 Einstein 模型(所有模同一個頻率、與 K 無關)其實是「光學支」的理想化。M₂ → ∞:1/M₂ → 0,ω² → 2C/M₁,重原子固定不動,輕原子在兩個固定牆之間振盪——標準的 Einstein 振子。M₁ = M₂ = M:ω² → 4C/M,這恰好是單原子鏈 (9) 的最大頻率 (4C/M)^{1/2}(發生在單原子鏈的區邊界)——這不是巧合,而是 zone folding,後面有專門的推導。注意 ω² 正比於 1/μ 而不是 1/(M₁+M₂):光學模是兩體相對運動,質心不動(見 (26)),所以出現的是約化質量,正如雙原子分子的振動頻率 √(k/μ)。
典型數值
KBr(我依 Fig. 8b 的 ω_T 反推):μ = M(K)M(Br)/(M(K)+M(Br)) = 39.10×79.90/119.00 = 26.25 u = 4.360×10⁻²⁶ kg。取 Fig. 8b 讀到的 ν_T(K=0) = 3.57×10¹² Hz ⇒ ω_T = 2π×3.57×10¹² = 2.243×10¹³ rad/s。由 (23):C = μω²/2 = 4.360×10⁻²⁶ × (2.243×10¹³)²/2 = 11.0 N/m。同法用縱波的 ν_L = 5.0×10¹² Hz 得 C_L = 21.5 N/m。兩者相差約兩倍,正確反映了 p.91 的提醒「The constant C is the force constant between nearest-neighbor planes and will differ for longitudinal and transverse waves」。

聲學支在長波極限(K → 0)的色散

ω² ≅ [½C/(M₁ + M₂)]·K²a² ⇒ ω ≅ a√[C/(2(M₁+M₂))]·K

(24)  p.98

ω:聲學支的角頻率;K:波向量;a:同種平面之間的重複距離(不是最近鄰平面間距!);C:力常數;M₁+M₂:一個原胞的總質量。分子的係數是 ½C。

出處推導
由 (22) 在 Ka ≪ 1 的極限取小根。補推導見前一個 block:因為 ω² = O(K²),可丟掉 O(K⁴) 的 ω⁴ 項,得 −2C(M₁+M₂)ω² + C²K²a² ≅ 0。(校註:PDF 文字層把 (24) 分子的 ½ 拆成浮出的「1 2」碎片,且把 ≅ 誤輸出成 ≡;我已對該頁算繪影像放大目檢,確認分子為 ½C、關係符為 ≅。)
極限行為

ω ∝ K,這就是聲波。斜率 v = ω/K = a√[C/(2(M₁+M₂))] 是常數,即聲速與頻率無關(無色散)——p.94 對單原子鏈的 (15) 講過同一句話:「the velocity of sound is independent of frequency in this limit」。出現的是原胞總質量 M₁+M₂ 而不是約化質量,因為聲學模是整個原胞一起平移(見 p.99:「The atoms (and their center of mass) move together」)——這與 (23) 的 μ 形成完美對照:光學支看相對運動(μ),聲學支看整體運動(M₁+M₂)。

M₁ = M₂ = M 時:ω² ≅ [½C/(2M)]K²a² = (C/4M)K²a²。乍看與單原子鏈 (15) 的 ω² = (C/M)K²a² 差 4 倍——但兩式裡的 a 不是同一個!雙原子鏈的 a 是同種平面間距 = 2×(單原子鏈的平面間距 a')。代入 a = 2a':ω² = (C/4M)K²(2a')² = (C/M)K²a'² ✓ 完全相同。這個檢查同時驗證了 (24) 與「a = nearest identical planes」的定義。

這條式子是第 5 章的入口。Debye 模型假設所有模都遵守 ω = vK,就是把 (24) 外推到整個布里淵區;Debye 溫度 θ = ℏv K_D/k_B 裡的 v 正是這裡的 a√[C/(2(M₁+M₂))]。

典型數值
KBr(我的計算,用前式擬合的 C 與書本 p.14 的 a_cube = 6.59 Å):沿 [111] 的同種平面重複距離 a = a_cube/√3 = 3.805 Å。M₁+M₂ = 119.00 u = 1.976×10⁻²⁵ kg。 橫波:v_TA = a√[C_T/(2(M₁+M₂))] = 3.805×10⁻¹⁰ × √[11.0/(3.952×10⁻²⁵)] = 3.805×10⁻¹⁰ × 5.27×10¹² = 2.0×10³ m/s。 縱波:v_LA = 3.805×10⁻¹⁰ × √[21.5/(3.952×10⁻²⁵)] = 3.805×10⁻¹⁰ × 7.38×10¹² = 2.8×10³ m/s。 兩者都落在離子晶體聲速的典型範圍(2–4 km/s)。誠實的對照:把 Fig. 8b 最低的 LA 資料點 (K/K_max ≈ 0.36, ν = 1.35×10¹² Hz) 拿來算弦斜率,v = ω/K = 2νa/(K/K_max) = 2×1.354×10¹²×3.805×10⁻¹⁰/0.357 ≈ 2.9×10³ m/s,只比模型預測的 2.8×10³ 高約 3%。也就是說,單參數最近鄰模型連中段 K 的弦斜率都抓得住;它在各處的偏差見後面的數值 block。

區邊界 Ka = π:(25) 的兩根,以及哪一根屬於哪一支補推導

書上關於區邊界只寫了兩句:「The extent of the first Brillouin zone is −π/a ≤ K ≤ π/a, where a is the repeat distance of the lattice. At K_max = ±π/a the roots are (25)」,然後列出 ω² = 2C/M₁ 與 ω² = 2C/M₂。書上沒有說哪一根屬於哪一支——這個資訊只出現在 Fig. 7 的標註裡。以下把計算與歸屬都補齊,而且給出一條完全不用計算的捷徑。

  1. 第 1 步(捷徑:從 (20) 直接看出來)。在 Ka = π 時,exp(iKa) = exp(−iKa) = exp(±iπ) = −1。所以 (20) 的兩個耦合係數都變成

    C[1 + exp(∓iKa)] = C[1 + (−1)] = 0

    兩條方程完全解耦:

    −ω²M₁u = −2Cu ⇒ ω² = 2C/M₁(此時 v 任意,取 v = 0)

    −ω²M₂v = −2Cv ⇒ ω² = 2C/M₂(此時 u 任意,取 u = 0)

    這就是 (25),一行算完,不需要碰行列式。

  2. 第 2 步(從 (22) 驗算)。Ka = π ⇒ 1 − cos π = 2,(22) 變成

    M₁M₂ω⁴ − 2C(M₁+M₂)ω² + 4C² = 0

    這個四項式可以直接因式分解:(M₁ω² − 2C)(M₂ω² − 2C) = 0

    (展開驗證:M₁M₂ω⁴ − 2CM₁ω² − 2CM₂ω² + 4C² = M₁M₂ω⁴ − 2C(M₁+M₂)ω² + 4C² ✓)

    兩根 ω² = 2C/M₁、2C/M₂ ✓ 與第 1 步一致。

  3. 第 3 步:本徵向量說明歸屬。第 1 步順帶給出了本徵向量:

    • ω² = 2C/M₁ 的模對應 v = 0:M₂ 面完全靜止,只有 M₁ 面在振盪。
    • ω² = 2C/M₂ 的模對應 u = 0:M₁ 面完全靜止,只有 M₂ 面在振盪。

    這是區邊界駐波的圖像(p.93 (12) 說明了區邊界的解是駐波 u_s = u(−1)^s):兩個子晶格交替「輪流」承擔全部的運動。

  4. 第 4 步:歸屬判定。設 M₁ > M₂(Fig. 7 的設定)。則 2C/M₁ < 2C/M₂——重原子振盪得慢。而聲學支永遠在光學支下面(由 (22) 的兩根排序,或由精確解的 ± 號),所以:

    區邊界的根哪個子晶格在動屬於Fig. 7 的標註
    ω² = 2C/M₁(小)只有重的 M₁ 面動聲學支頂端(2C/M₁)^{1/2} 標在聲學支末端
    ω² = 2C/M₂(大)只有輕的 M₂ 面動光學支底端(2C/M₂)^{1/2} 標在光學支末端

    (我已對 p.96 的 Fig. 7 算繪目檢,確認標註位置如上。)

  5. 第 5 步:頻隙寬度。禁止傳播的頻率區間是

    (2C/M₁)^{1/2} < ω < (2C/M₂)^{1/2}

    寬度由質量比決定:ω_上/ω_下 = √(M₁/M₂)。M₁ = M₂ 時比值為 1,頻隙關閉。

這個結果給了一個實驗上的量測手段。如果你在中子散射資料上讀到聲學支的頂端與光學支的底端,兩者平方的比值就直接給出質量比 M₁/M₂(在單參數最近鄰模型的框架下)。反過來,質量已知時可以求 C,而且可以與 K = 0 的 (23) 求出的 C 交叉檢驗——同一個模型的兩端若給出一致的 C,表示最近鄰近似在這個材料上站得住。後面的 KBr 數值 block 就是做這件事。

與單原子鏈的對照。單原子鏈在區邊界只有一個頻率 (4C/M)^{1/2},沒有頻隙,因為只有一條支。雙原子鏈把布里淵區折半,把原本落在外半區的部分折回來變成光學支,兩支在區邊界處分開——分開的幅度就是頻隙。若 M₁ = M₂,折回來的兩段在區邊界剛好接上(都等於 2C/M),頻隙為零,這時你看到的其實還是同一條單原子鏈的曲線,只是被畫在半個區裡。下一個 block 嚴格證明這件事。

區邊界(K = ±π/a)的兩個頻率

ω² = 2C/M₁ (聲學支頂端,只有重原子面在動); ω² = 2C/M₂ (光學支底端,只有輕原子面在動)

(25)  p.98

M₁ > M₂(採 Fig. 7 的設定);C:力常數;a:同種平面重複距離;區邊界 K_max = ±π/a。(書上的 (25) 只並列寫出兩個 ω² 值,括號內的歸屬說明是我依 Fig. 7 的標註與 (20) 在 Ka = π 的解耦結構補上的。)

出處推導
由 (22) 在 Ka = π 取值(1 − cos π = 2),方程可因式分解為 (M₁ω² − 2C)(M₂ω² − 2C) = 0。更快的路是回到 (20):Ka = π 時 exp(±iKa) = −1,兩個耦合係數 C[1 + exp(∓iKa)] 同時歸零,兩條方程解耦,各自給出一個獨立振子。
極限行為
M₁ = M₂ = M 時兩根合而為一(都等於 2C/M),頻隙關閉——這正是 Ge 這種同種原子基底在本模型下的預測(但實驗上 Ge 仍有頻隙,見後面的討論 block)。M₂/M₁ → 0(極輕的第二種原子):光學支底端 2C/M₂ → ∞,頻隙無限寬。真實例子是含氫的化合物:H 比其他原子輕一個量級,光學支被推得極高,與聲學支之間出現巨大的頻隙。這個頻隙是「彈性波的禁帶」,p.99 直接這樣說:「Wavelike solutions do not exist for certain frequencies, here between (2C/M₁)^{1/2} and (2C/M₂)^{1/2}. This is a characteristic feature of elastic waves in polyatomic lattices.」——落在這個區間的頻率,K 會變成複數,波在空間中指數衰減。這是聲子晶體(phononic crystal)與聲學超材料的物理基礎。
典型數值
KBr,用 (23) 擬合出的 C(我的計算,Fig. 8b 讀值見前): 橫波組(C_T = 11.0 N/m):ω² = 2C/M(Br) ⇒ ν = 2.05×10¹² Hz(預測 TA 頂端),實測讀值 2.15 ⇒ 差 5%;ω² = 2C/M(K) ⇒ ν = 2.93×10¹² Hz(預測 TO 底端),實測讀值 3.07 ⇒ 差 5%。 縱波組(C_L = 21.5 N/m):ω² = 2C/M(Br) ⇒ ν = 2.87×10¹² Hz(預測 LA 頂端),實測 2.76 ⇒ 差 4%;ω² = 2C/M(K) ⇒ ν = 4.10×10¹² Hz(預測 LO 底端),實測 4.30 ⇒ 差 5%。 四個獨立預測全部落在 5% 以內,而 C 只由 K = 0 的一個數字定出——對一個只有一個可調參數的模型而言,這是相當強的驗證。

Fig. 7 逐筆讀圖:雙原子鏈色散的定性全貌與頻隙

網頁不放圖,以下用文字重建。圖說逐字是:「Optical and acoustical branches of the dispersion relation for a diatomic linear lattice, showing the limiting frequencies at K = 0 and K = K_max = π/a. The lattice constant is a.」

圖面內容(以下為我對書本 p.96 算繪影像的目視描述):

  • 橫軸標 K,只畫正半邊(0 到區邊界);區邊界處畫一條垂直虛線,底下標 π/a。縱軸沒有標籤,只用文字標出各個極限值。
  • 上面一條曲線標 Optical phonon branch。它從縱軸上的一個高點出發,這個高點被標成 [2C(1/M₁ + 1/M₂)]^{1/2}——正是 (23) 開根號。曲線一開始幾乎水平,中段轉折下彎,到區邊界處趨平,末端標 (2C/M₂)^{1/2}。
  • 下面一條曲線標 Acoustical phonon branch。它從原點出發,一開始是直線(斜率即聲速,對應 (24)),逐漸彎曲、飽和,在區邊界處水平地停在 (2C/M₁)^{1/2}。
  • 兩條曲線之間、圖的中段標著 M₁ > M₂ — 這是本圖的前提設定。
  • 因為 M₁ > M₂,(2C/M₁)^{1/2} < (2C/M₂)^{1/2},兩條曲線在區邊界沒有接上,中間留下一段空白——那就是頻隙。

這張圖要記住的四個座標:

位置聲學支光學支
K = 00[2C(1/M₁+1/M₂)]^{1/2},即 (23)
K = π/a(2C/M₁)^{1/2},即 (25) 較小的根(2C/M₂)^{1/2},即 (25) 較大的根

頻隙的物理(p.99,緊接本段之後):「Wavelike solutions do not exist for certain frequencies, here between (2C/M₁)^{1/2} and (2C/M₂)^{1/2}. This is a characteristic feature of elastic waves in polyatomic lattices. There is a frequency gap at the boundary K_max = ±π/a of the first Brillouin zone.」

落在頻隙裡的頻率,(22) 解出來的 K 是複數,exp(iKx) 變成指數衰減——波進不去晶體。這與第 7 章電子能隙的數學結構完全相同(都是「波在週期位勢中的禁帶」),也是為什麼 Kittel 在 p.94 把區邊界的駐波類比成 Bragg 反射:「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.」

Fig. 7 與 Fig. 8 的關係。Fig. 7 是模型的定性圖(一個偏振、兩條支),Fig. 8 是真實材料的實驗資料(三個偏振、六條支)。要把兩者對起來,需要把 Fig. 8 拆成「橫波組」(TA + TO)與「縱波組」(LA + LO),每一組各自長得像一張 Fig. 7,只是兩組的 C 不同。這正是後面數值 block 的做法。

(23) 代回 (20) 得 (26):兩個子晶格的振幅比補推導

書上寫:「For the optical branch at K = 0 we find, on substitution of (23) in (20), (26)」——「on substitution」四個字包了三行代數。這三行是本段的終點,因為它把抽象的 ω² 值翻譯成具體的運動圖像。

  1. 第 1 步:取 K = 0 的 (20) 第一式。exp(−iKa)|_{K=0} = 1,所以

    −ω²M₁u = Cv(1 + 1) − 2Cu = 2Cv − 2Cu

  2. 第 2 步:代入 (23) 的 ω²。ω² = 2C(1/M₁ + 1/M₂):

    −2C(1/M₁ + 1/M₂)·M₁·u = 2Cv − 2Cu

    左邊展開:−2C(1 + M₁/M₂)u

  3. 第 3 步:兩邊除以 2C(C ≠ 0)。

    −(1 + M₁/M₂)u = v − u

    即 −u − (M₁/M₂)u = v − u

  4. 第 4 步:兩邊的 −u 相消。

    −(M₁/M₂)u = v

    ⇒ u/v = −M₂/M₁ ✓ 這是 (26)。

  5. 第 5 步:用 (20) 第二式獨立驗算。K = 0 時:−ω²M₂v = 2Cu − 2Cv。代入同樣的 ω²:

    −2C(1/M₁ + 1/M₂)M₂v = 2Cu − 2Cv ⇒ −(M₂/M₁ + 1)v = u − v ⇒ −(M₂/M₁)v = u

    ⇒ u/v = −M₂/M₁ ✓ 兩式給出同一個比值——這正是行列式為零所保證的(兩條方程線性相依)。若兩式給出不同的比值,表示 ω² 代錯了。

  6. 第 6 步:順帶得到聲學支的振幅比。K = 0 的聲學支是 ω = 0。代入 (20) 第一式:0 = 2Cv − 2Cu ⇒ u = v。這與 p.99 的說法一致:「Another solution for the amplitude ratio at small K is u = v, obtained as the K = 0 limit of (24). The atoms (and their center of mass) move together, as in long wavelength acoustical vibrations, whence the term acoustical branch.」
  7. 第 7 步:質心檢查。原胞的質心位移正比於 M₁u + M₂v。

    光學支:代入 u = −(M₂/M₁)v ⇒ M₁u + M₂v = M₁(−M₂/M₁)v + M₂v = −M₂v + M₂v = 0 ✓ 質心不動。

    聲學支:u = v ⇒ M₁u + M₂v = (M₁+M₂)u ≠ 0,質心跟著一起移動 ✓。

    這正好解釋了為什麼 (23) 出現約化質量 μ(相對運動)而 (24) 出現總質量 M₁+M₂(質心運動)。

(26) 的三個讀法。

(一)代數。u/v = −M₂/M₁ 是負數,所以兩個子晶格反相。而且比值的大小是質量的反比:輕的原子振幅大、重的原子振幅小。若 M₂ ≪ M₁(例如 H 與重原子),則 |u/v| = M₂/M₁ ≪ 1——重原子幾乎不動,光學模幾乎全是輕原子在跳。

(二)力學。這就是一個雙原子分子的伸縮振動:兩個質點用彈簧連著,質心不動,頻率 √(k/μ)。(23) 的 ω² = 2C/μ 完全是這個形式(2C 是因為每個原子兩側各有一根彈簧)。換句話說,K = 0 的光學模就是整個晶體裡每一個原胞同步做同一個分子振動。

(三)電學。若 M₁、M₂ 帶相反電荷,質心不動但電荷中心會動 ⇒ 原胞出現一個振盪的電偶極矩。這是下一個 block 的主題,也是「光學」這個名字的由來。

光學支在 K = 0 的振幅比

u/v = −M₂/M₁

(26)  p.98

u:M₁ 子晶格的振幅;v:M₂ 子晶格的振幅;負號表示反相;比值大小為質量的反比。

出處推導
把 (23) 的 ω² = 2C(1/M₁ + 1/M₂) 代回 (20)(在 K = 0,兩個指數都等於 1)。書上只寫「on substitution of (23) in (20)」;完整七步(含用第二式獨立驗算、以及質心檢查)見前一個補推導 block。
極限行為
M₁ = M₂ 時 u/v = −1:兩個子晶格等幅反相,這是 Ge 這類同種原子基底的光學模。M₂ ≪ M₁ 時 |u/v| ≪ 1:重原子幾乎不動,光學模幾乎全由輕原子承擔——這是所有含氫化合物(以及含輕元素的半導體,如 GaN 中的 N)光學支頻率特別高、且與聲學支分得特別開的原因。質心固定:M₁u + M₂v = M₁(−M₂/M₁)v + M₂v = 0,書上的下一句就是「The atoms vibrate against each other, but their center of mass is fixed.」對照聲學支:K = 0 的聲學支給 u = v(p.99),質心跟著動,這正是「整個原胞一起平移」——兩者的對比就是「相對座標 vs 質心座標」的分解。這個比值只在 K = 0 成立;p.99 明說「At a general K the ratio u/v will be complex, as follows from either of the equations (20)」。
典型數值
KBr:M₁ = M(Br) = 79.90 u,M₂ = M(K) = 39.10 u ⇒ u/v = −39.10/79.90 = −0.489。也就是說在 K = 0 的光學模裡,K⁺ 的振幅大約是 Br⁻ 的兩倍,兩者方向相反。(此為我依標準原子量的計算,非書中數字。)

為什麼叫「光學支」:從 (26) 到光的電場

這是本段的收尾,也是整個名稱的由來。書本 p.98 最後一段接到 p.99 開頭,完整的句子是:

「The atoms vibrate against each other, but their center of mass is fixed. If the two atoms carry opposite charges, as in Fig. 10, we may excite a motion of this type with the electric field of a light wave, so that the branch is called the optical branch.」

(此句橫跨 p.98 與 p.99:p.98 止於「we may excite a motion of this」,其餘在 p.99。)

把這句話拆成四個邏輯步驟:

第 1 步:(26) 說兩個原子反相振動、質心固定。這一步純力學,不需要電荷。

第 2 步:若兩個原子帶相反電荷(離子晶體,如 KBr 的 K⁺ 與 Br⁻),反相的位移就把正負電荷中心拉開 ⇒ 每個原胞出現一個振盪的電偶極矩 p(t) ∝ (u − v)·q。注意:質心固定但電偶極矩不為零,兩者不衝突——質心看的是質量加權,偶極矩看的是電荷加權。

第 3 步:光波的電場可以直接抓住這個偶極矩。外加電場 E 對正離子施力 +qE、對負離子施力 −qE,這是一組反相的驅動力——與光學模的位移模式完全匹配。相對地,聲學模的 u = v(同相),電場的反相力無法驅動它。所以光只能激發光學支,不能激發聲學支。

第 4 步:於是這條分支被命名為「optical branch」。對照 p.99 給聲學支的命名理由:「The atoms (and their center of mass) move together, as in long wavelength acoustical vibrations, whence the term acoustical branch.」——因為長波下整個原胞一起動,就像空氣中的聲波,所以叫「聲學」。

必須澄清的一個誤解。Kittel 的條件句用的是「If the two atoms carry opposite charges」——帶電荷不是光學支存在的條件,只是它能被光激發的條件。Ge 的兩個原子都是中性的 Ge,光學支照樣存在(Fig. 8a 的 LO、TO),只是它不是紅外活性的(不能被光直接一階激發,但可以被拉曼散射看到)。名字是從離子晶體來的,套用到所有材料上。(「不是紅外活性但拉曼活性」這一句是我的補充,本段未提;Kittel 在此只說「If ... we may excite」。)

這條線往下接到哪裡。第 14 章 p.414 就是把這件事做嚴格化:光子與 TO 聲子的耦合產生 polariton,LO 模因為 div E ≠ 0 而被宏觀電場推高,兩者的頻率比給出 LST 關係 ω_L²/ω_T² = ε(0)/ε(∞)(該章的 (62)),而 p.415 的 Fig. 15、Fig. 16 展示了 NaCl 與 LiF 在這個頻率區間的高反射率——那就是本段的頻隙在光學上的直接後果。

Fig. 10 逐筆讀圖:光學模與聲學模的位移長什麼樣

網頁不放圖,以下用文字重建。圖說逐字是:「Transverse optical and transverse acoustical waves in a diatomic linear lattice, illustrated by the particle displacements for the two modes at the same wavelength.」

注意圖說的最後一句:兩個模畫的是同一個波長。這是刻意安排的——同樣的 λ(同樣的 K),兩種完全不同的運動。

圖面內容(以下為我對書本 p.98 算繪影像的目視描述):

圖分上下兩格,各有一個向右的粗箭頭標 K,表示波沿水平方向傳播,而位移是垂直的(橫波)。

上格:Optical mode。

  • 畫了兩條分開的曲線,彼此交叉,形成一連串像「眼睛」或「透鏡」的閉合形狀。
  • 一條曲線上串著實心黑圓,圓內是「+」號(正離子);另一條串著空心白圓,圓內是「−」號(負離子)。
  • 兩條曲線在若干點交會,交會點上是 + 離子(該處位移為零,兩個子晶格都通過平衡位置)。
  • 關鍵在於:在兩交會點之間,+ 曲線往上凸的時候 − 曲線往下凹。也就是說正負離子子晶格彼此反相 ⇒ 這正是 (26) 的 u/v < 0 的圖像。
  • 圖上 + 與 − 的振幅看起來相近(圖是示意,沒有按 M₂/M₁ 的比例畫)。

下格:Acoustical mode。

  • 只畫了一條平滑的正弦曲線,所有離子(+ 與 −)都串在同一條曲線上,交替排列 +、−、+、−。
  • 相鄰的 + 與 − 位移幾乎相同(只差曲線上一小段的差別)⇒ 同相 ⇒ 這是 u = v 的圖像。
  • 波長與上格相同(可以數波峰對照)。

這張圖一眼看懂的三件事:

  1. 光學 = 兩條曲線(兩個子晶格各走各的);聲學 = 一條曲線(兩個子晶格共用一條)。這是分辨兩種模最快的視覺特徵。
  2. 同樣的波長,光學模在「原胞內部」有結構,聲學模沒有。光學模的短程劇烈變化(相鄰原子反向)就是它頻率高的原因——彈簧被拉伸得更厲害。
  3. 電偶極矩。上格中每一對相鄰的 +/− 都因反相而分開了一段距離 ⇒ 產生淨偶極;下格中 +/− 一起移動 ⇒ 偶極矩不變。這就是上一個 block 第 2、3 步的圖解。

與 p.98 正文的接點。正文說「The particle displacements in the transverse acoustical (TA) and transverse optical (TO) branches are shown in Fig. 10」——Kittel 特地畫橫波,因為橫波的位移垂直於傳播方向,畫在紙面上最清楚。縱波的位移平行於 K,畫出來會變成疏密圖,不好看。物理是一樣的。

等質量極限:雙原子鏈嚴格退化回單原子鏈(zone folding)補推導

Kittel 說過「It is not essential that the masses be different」,但沒有把 M₁ = M₂ 的情形算出來。這個極限值得算,有三個理由:(一)它是驗證整套代數的最強檢查;(二)它嚴格證明了「a = 最近同種平面間距」這個定義是自洽的;(三)它給出 zone folding(能帶折疊) 的第一個例子——這個概念在第 7 章談能帶、以及在超晶格與二維材料的文獻裡到處都是。

設定:M₁ = M₂ = M。雙原子鏈的重複距離是 a(同種平面間距),所以最近鄰平面間距是 a' = a/2。若把它當成單原子鏈看,平面間距就是 a'。

  1. 第 1 步:把 M₁ = M₂ = M 代入精確解。由前面的精確解 ω²_± = (C/μ)·{1 ± √[1 − (2μ/(M₁+M₂))(1 − cos Ka)]},此時 μ = M/2、M₁+M₂ = 2M,所以 2μ/(M₁+M₂) = M/(2M) = 1/2,而 C/μ = 2C/M。

    ω²_± = (2C/M)·{1 ± √[1 − ½(1 − cos Ka)]}

  2. 第 2 步:化簡根號內。1 − ½(1 − cos Ka) = ½(1 + cos Ka) = cos²(Ka/2)(用半角公式 1 + cos θ = 2cos²(θ/2))。

    根號 = |cos(Ka/2)|。於是

    ω²_± = (2C/M)·[1 ± |cos(Ka/2)|]

  3. 第 3 步:聲學支(下號)化成單原子鏈形式。在第一布里淵區 0 ≤ Ka ≤ π 內 cos(Ka/2) ≥ 0,可去絕對值:

    ω²₋ = (2C/M)[1 − cos(Ka/2)] = (2C/M)·2sin²(Ka/4) = (4C/M)·sin²(Ka/4)

    現在把 a = 2a' 代入:Ka/4 = K(2a')/4 = ½Ka'。所以

    ω²₋ = (4C/M)sin²(½Ka')

    這與書本 p.92 的 (9)「ω² = (4C/M)sin²(½Ka)」逐字相同(該處的 a 是單原子鏈的平面間距,即我們的 a')✓

  4. 第 4 步:光學支(上號)是折回來的那一段。

    ω²₊ = (2C/M)[1 + cos(Ka/2)] = (2C/M)·2cos²(Ka/4) = (4C/M)·cos²(Ka/4)

    用 cos θ = sin(π/2 − θ):cos²(Ka/4) = sin²(π/2 − Ka/4) = sin²[(2π − Ka)/4]。

    定義 K' ≡ 2π/a − K(把 K 對單原子鏈的區邊界 π/a' = 2π/a 做鏡射),則 K'a/4 = (2π − Ka)/4,於是

    ω²₊ = (4C/M)sin²(K'a/4) = (4C/M)sin²(½K'a')

    這也是單原子鏈的 (9),只是取在 K' 而非 K。

  5. 第 5 步:解讀——這就是 zone folding。單原子鏈的第一布里淵區是 |K| ≤ π/a' = 2π/a。雙原子鏈的第一布里淵區是 |K| ≤ π/a,恰好是前者的一半。

    • 單原子鏈曲線落在內半區 0 ≤ K ≤ π/a 的那一段,就是雙原子鏈的聲學支(第 3 步)。
    • 單原子鏈曲線落在外半區 π/a ≤ K ≤ 2π/a 的那一段,以 K → 2π/a − K 折回內半區,就是雙原子鏈的光學支(第 4 步)。
  6. 第 6 步:三點驗算。 (a) K = 0 的光學支:ω²₊ = (4C/M)cos²(0) = 4C/M。而 (23) 在 M₁=M₂=M 給 ω² = 2C(1/M + 1/M) = 4C/M ✓。同時 (4C/M)^{1/2} 正是單原子鏈 (9) 的最大頻率(發生在 K = π/a' = 2π/a,即被折回 K = 0 的那一點)✓。 (b) 區邊界:Ka = π ⇒ cos(π/2) = 0 ⇒ 兩根都等於 2C/M,頻隙關閉 ✓ 與 (25) 在 M₁=M₂ 的結果一致。而單原子鏈在 K = π/a = ½(2π/a)(區中點)的值是 (4C/M)sin²(π/4) = 2C/M ✓ 完全吻合。 (c) 小 K 的聲速:ω²₋ ≅ (4C/M)(½Ka')² = (C/M)K²a'² ✓ 與 p.94 的 (15) 逐字相同。而由 (24) 在 M₁=M₂=M 得 ω² ≅ [½C/(2M)]K²a² = (C/4M)K²(2a')² = (C/M)K²a'² ✓ 同一個答案。

這個檢查同時驗了三件事。(一)從 (18) 到 (25) 的整條代數鏈沒有算錯,因為它在 M₁ = M₂ 時精準退化回本章前半已經獨立推導過的結果。(二)「a 是最近同種平面間距」這個定義是唯一自洽的選擇——如果 a 取成最近鄰間距,第 3 步就會多出一個 4 倍因子,與 (9) 對不上。(三)把布里淵區折半、多出來的分支從外半區折回來——這個「zone folding」機制是普遍的。

zone folding 為什麼重要。只要你把原胞人為放大(例如把兩個相同的原胞硬說成一個大原胞),分支數就會加倍,但物理沒有改變——多出來的「光學支」只是原來那條支被折進去的部分。真正的物理判準是:折回來的兩段在新的區邊界是否分開。若晶體真的有兩種不等價的位置(質量不同、或力常數不同),就會在區邊界打開頻隙,光學支才是真的;若兩個位置完全等價,頻隙為零,折疊只是記帳方式的改變。

這個判準在第 7 章會以完全相同的形式出現:週期位勢在布里淵區邊界打開能隙,能隙的大小正比於位勢的傅立葉分量 U_G;若 U_G = 0,能隙關閉,「能帶」只是自由電子拋物線的折疊。晶格振動與電子能帶,在這一點上是同一個數學。(與第 7 章的類比為我的補充。)

把 KBr 的數字真的代進去:一個參數、四個獨立預測補推導

這一段的模型只有一個可調參數 C(每個偏振各一個)。這意味著:只要用一個實驗數字把 C 定出來,其餘所有預測都是零參數的。以下用 Fig. 8b 的 KBr 資料做這件事。(所有數值為我的計算;Fig. 8b 的讀值為對算繪影像的目視判讀,原子量為標準值,a_cube = 6.59 Å 取自書本 p.14。)

步驟 1:定出 C。用 (23) 與 K = 0 的外插值。

μ = M(K)·M(Br)/[M(K)+M(Br)] = 39.10 × 79.90 / 119.00 = 26.25 u = 4.360×10⁻²⁶ kg

偏振Fig. 8b 讀到的 K=0 值ω = 2πνC = μω²/2
橫(ω_T)ν = 3.57×10¹² Hz2.243×10¹³ rad/s11.0 N/m
縱(ω_L)ν = 5.0×10¹² Hz3.142×10¹³ rad/s21.5 N/m

縱波的力常數約為橫波的兩倍——這與 p.91 的提醒一致:「The constant C ... will differ for longitudinal and transverse waves.」

步驟 2:用 (25) 預測區邊界,與圖上的資料點對照。這是零參數預測。M(Br) = 1.327×10⁻²⁵ kg(重,配聲學支),M(K) = 6.492×10⁻²⁶ kg(輕,配光學支)。

預測對象公式模型給的 νFig. 8b 實測讀值誤差
TA 頂端(2C_T/M_Br)^{1/2}/2π2.05×10¹² Hz2.15×10¹² Hz−5%
TO 底端(2C_T/M_K)^{1/2}/2π2.93×10¹² Hz3.07×10¹² Hz−5%
LA 頂端(2C_L/M_Br)^{1/2}/2π2.87×10¹² Hz2.76×10¹² Hz+4%
LO 底端(2C_L/M_K)^{1/2}/2π4.10×10¹² Hz4.30×10¹² Hz−5%

四個獨立預測全部落在 5% 以內。對一個只有最近鄰、單一力常數、把三維壓成一維的玩具模型來說,這個表現遠超預期。它說明 KBr 沿 [111] 的動力學確實被「交錯的 K 面與 Br 面 + 一個最近鄰彈簧」抓住了主要部分。

步驟 3:誠實報告中段 K 的偏差。把同樣的 C 代回精確解,算中段 K 的聲學支:

  • 橫波最低的 TA 點 K/K_max ≈ 0.35(Ka = 0.345π):模型給 ν = 0.89×10¹² Hz,實測 0.97 ⇒ 低估 8%。
  • 縱波最低的 LA 點 K/K_max ≈ 0.36(Ka = 0.357π):模型給 ν = 1.29×10¹² Hz,實測 1.35 ⇒ 低估 5%。

整段 [111] 色散在這個單參數模型下,偏差都落在 10% 以內——區邊界四個點在 5% 以內,中段 K 略微偏低(4–8%)。換句話說,這份資料上還看不出非加次近鄰不可的必要性:殘餘的幾個百分點與讀圖誤差同量級。若真要把這幾個百分點也做掉,做法就在 p.94–95 那一節(Derivation of Force Constants from Experiment,(16a)、(16b)、(17)):把實驗的 ω²(K) 做傅立葉分析,一口氣把 C₁、C₂、C₃… 全部反解出來。

這個練習示範的方法論。「用一個數定參數 → 用其餘的數檢驗 → 誠實報告偏差 → 指出下一階修正在哪裡」——這是物理建模的標準流程,比記住 (23)(24)(25) 三條式子重要得多。

Ge 為什麼在區邊界仍有頻隙:Kittel 這個模型漏掉的東西我的判讀

這是本段最值得深究的一個矛盾,而 Kittel 沒有明說。

矛盾在哪裡。Ge 的兩個基底原子都是鍺,質量完全相同:M₁ = M₂。依 (25),兩個根都等於 2C/M,頻隙應該關閉。但 Fig. 8a 的實驗資料清楚顯示:在區邊界,LA 停在 ≈ 6.55×10¹² Hz,LO 停在 ≈ 7.3×10¹² Hz,中間有約 0.7×10¹² Hz 的間隙(目視讀值)。模型預測與實驗矛盾。

線索就在 p.97 的那句免責聲明。Kittel 寫:「It is not essential that the masses be different, but either the force constants or the masses will be different if the two atoms of the basis are in nonequivalent sites.」——他明白指出還有另一個來源:力常數。但他隨後又假設「the force constants are identical between all pairs of nearest-neighbor planes」,把這個來源關掉了。所以 (25) 的頻隙只剩質量差這一條路。

沿 [111] 的金剛石結構,平面間距是不等的。用書本 p.17 的 a_cube(Ge) = 5.658 Å 算(以下為我的計算):

  • fcc 格子的 (111) 面間距 = a_cube/√3 = 3.266 Å,這就是本段意義下的 a(同種平面的重複距離),與 Fig. 8a 圖說的 K_max = (2π/a_cube)(½½½) 相符(其長度 = π√3/a_cube = π/3.266 Å⁻¹ = π/a ✓)。
  • 第二個基底原子在 (¼,¼,¼)a_cube,它沿 [111] 的投影 = (3a_cube/4)/√3 = a_cube√3/4 = 2.450 Å(這也正是 Ge 的最近鄰鍵長 √3a_cube/4)。
  • 所以在一個週期 3.266 Å 之內,兩種平面分別位於 0 與 2.450 Å ⇒ 相鄰平面的間距在 2.450 Å 與 0.816 Å 之間交替,比值恰為 3 : 1。

結論:Ge 的 [111] 方向根本不滿足「所有最近鄰平面對力常數相同」。兩種間距差三倍,對應的力常數必然差很多。這是一個「等質量、不等力常數」的雙原子鏈——它同樣有光學支,同樣在區邊界有頻隙,只是頻隙由 |C₁ − C₂| 而非 |M₁ − M₂| 決定。

對照 KBr 為什麼就準。KBr 是 NaCl 結構,沿 [111] 的 K 面與 Br 面是等間距的(都是 a_cube/(2√3) = 1.902 Å,我的計算)。所以「所有最近鄰平面對共用一個 C」在 KBr 上是合理的,而頻隙由質量差 79.90 : 39.10 提供——這正是前一個 block 的預測能準到 5% 的原因。

兩個材料放在一起,剛好示範了 Kittel 那句話的兩半:

質量差力常數差頻隙來源Kittel 的模型適用嗎
KBr [111]有(2:1)無(等間距)質量適用,預測準到 5%
Ge [111]無(同種原子)有(間距 3:1)力常數不適用,(25) 會預測頻隙為零

這是本段最容易被跳過、卻最有教學價值的一點。教科書給的模型不會標示自己的失效範圍;讀者要自己把模型的假設一條一條拿去對材料。(本 block 的幾何計算與結論為我的分析;書中只給了那句「either the force constants or the masses will be different」的提示,沒有對 Ge 做這個討論。)

常見誤解與自我檢查清單我的判讀

把這一段最容易出錯的地方集中列出來。左欄是常見說法,右欄是為什麼錯、以及書上哪裡可以核對。

常見誤解為什麼錯
「a 是最近鄰平面的間距」p.97 明說「the distance between nearest identical planes, not nearest-neighbor planes」。取錯會讓區邊界差 2 倍、(24) 的 ω² 差 4 倍。用「等質量極限必須退化回 (9)」來自我檢查。
「光學支就是頻率高的那一支」Fig. 8a 的 Ge:兩條光學支在 K/K_max ≲ 0.45 之間簡併(畫成同一條線),之後分岔而非交叉,區邊界 TO(8.7) 高於 LO(7.3)。分支的標籤來自偏振與 K→0 行為,不是頻率排序。
「原胞裡 p 個原子就有 p 條分支」是 3p 條。漏掉的因子 3 是空間的三個偏振方向(p.91)。
「聲學支有 3 條是因為有 3 個原子」聲學支永遠 3 條,與 p 無關。它對應三個剛體平移方向。
「每條分支上的 K 值數是 3N」是 N(腳註 1:N = V/V_c)。3pN 是總模數 = 3p 條分支 × N 個 K。
「M₁ = M₂ 就沒有光學支」錯。p.97:「It is not essential that the masses be different」。Ge 就是反例。等質量只讓區邊界的頻隙關閉(在等力常數的前提下),不讓光學支消失。
「(25) 的 2C/M₁ 是光學支」反了。M₁ > M₂ 時 2C/M₁ 是聲學支頂端(重的振盪慢),2C/M₂ 是光學支底端。核對 Fig. 7 的標註,或用 (20) 在 Ka = π 解耦的本徵向量判斷。
「光學支一定能被光激發」只有兩個原子帶相反電荷時才行。p.98 的句子是條件句「If the two atoms carry opposite charges」。Ge 的光學支不是紅外活性的。
「(26) 的 u/v = −M₂/M₁ 對所有 K 成立」只在 K = 0 成立。p.99:「At a general K the ratio u/v will be complex」。
「(22) 是對 ω 的四次方程,要解四個根」它只含 ω 的偶次冪,是對 ω² 的二次方程。兩個 ω² 根、四個 ω 根(±ω₊、±ω₋),負頻率不是新的物理解(p.99 末有討論「What is the sign of ω?」)。
「Fig. 8 只有四條曲線,所以只有四條分支」沿 [111] 兩個橫偏振簡併,TA 與 TO 各代表兩條。Fig. 8a 圖說自己就寫「The two TA phonon branches」。

兩個隨手可用的驗算工具:

(一)Vieta 和。任何 K 下,兩支的頻率平方和恆為 ω₊² + ω₋² = 2C(1/M₁ + 1/M₂),與 K 無關。拿實驗色散圖的任兩支相加,應該得到常數——若不是常數,表示最近鄰單參數模型在該材料上不成立。

(二)等質量退化。任何你推導出的雙原子鏈結果,代 M₁ = M₂ = M、a = 2a',都必須回到單原子鏈 (9) 的對應結果。這條檢查抓得到絕大多數的代數錯誤與 a 定義錯誤。

這一段往後接到哪裡我的判讀

把本段的結論一條一條接到後面的章節,可以看出它在全書結構裡的位置。

接第 4 章後半(p.99 起,QUANTIZATION OF ELASTIC WAVES)。本段建立了「3pN 個模」這個骨架;下一節立刻把每個模量子化成諧振子,ε = (n + ½)ℏω(該節的 (27))。模的數目是本段給的,能量是下一節給的——兩者相乘才是晶體的振動能。

接第 5 章(Phonons II,熱性質)。

  • 腳註 1 的「一支 N 個 K」在第 5 章被正式證明(週期邊界條件),並直接長成態密度 D(ω)。
  • (24) 的 ω = vK 是 Debye 模型的唯一輸入:Debye 假設所有 3 條聲學支都遵守這條線性色散一路到 K_D,由此得到低溫的 T³ 熱容。
  • (23) 的「ω 與 K 無關」是 Einstein 模型的理想化:所有模同一個頻率。真實的光學支確實近乎平坦(Fig. 8a 的 TO 全程只從 9.15 降到 8.7),所以 Einstein 模型對光學支的熱容描述得意外地好。
  • 本段解釋了為什麼需要兩個模型:一個晶體同時有線性的聲學支與平坦的光學支,單一模型裝不下。

接第 14 章(Plasmons, Polaritons, and Polarons,p.393 起)。

  • 本段 Fig. 8b 定義的 ω_T、ω_L,在 p.414 的 (62) 組成 LST 關係 ω_L²/ω_T² = ε(0)/ε(∞),而 Kittel 特別註明該推導「assumed a cubic crystal with two atoms per primitive cell」——就是本段的模型。
  • p.414 的 Fig. 14 解釋了本段沒說的 LO–TO 分裂機制(LO 有宏觀電場、TO 沒有)。
  • p.415 的 Fig. 15(NaCl 反射率)把本段的頻隙變成可以直接量測的光學現象。
  • p.414 還有一句與本段呼應:「For soft modes with ω_T → 0 we see that ε(0) → ∞, a characteristic of ferroelectricity.」——軟模鐵電的起點就是本段的 TO 支。

接第 7 章(能帶)。「等質量極限 = zone folding」這個結構與第 7 章「弱週期位勢在區邊界打開能隙」是同一個數學:把原胞放大 ⇒ 布里淵區縮小 ⇒ 分支折疊 ⇒ 若兩個位置真的不等價,折疊點就打開一個隙。晶格振動的頻隙與電子的能隙,來源相同。

接熱傳導與元件熱模型(與 Cth 萃取研究的接點)。這是我加的一段跨領域連結,書中沒有:

  • 只有聲學支輸運熱。(23) 說光學支在 K → 0 是平的 ⇒ v_g ≈ 0 ⇒ 幾乎不搬能量。熱導率 κ = ⅓Cvℓ 裡的 v 基本上由 (24) 的聲速決定。光學支的角色是「熱容的倉庫」與「聲學聲子散射的通道」(Umklapp 過程常需要一個光學聲子參與),不是輸運者。
  • 原胞越大,熱導越低。原胞含 p 個原子時,3 條聲學支不變,但光學支增加到 3p − 3 條。光學支多、聲學支的可用散射通道多、且平均群速度被稀釋 ⇒ 複雜原胞的材料熱導低。這是熱電材料刻意選用大原胞化合物的設計原則之一。
  • 對 FinFET 的熱阻(Cth/Rth)萃取的意義:矽的原胞有 2 個原子(金剛石結構,本段的 p = 2),光學支占了 6 條裡的 3 條;在室溫附近,光學模已被大量佔據,貢獻可觀的熱容但幾乎不貢獻熱導。所以「熱容」與「熱導」在聲子層次上是由不同的分支主導的——這在做集總參數熱模型(Cth–Rth 網路)的物理詮釋時是一個容易被忽略的區別。
  • 奈米尺度的額外效應:腳註 1 的 N = V/V_c 在樣品極小時 N 不再是大數,K 點離散化,聲子譜出現可分辨的間隔;同時聲學支的長波模(就是熱傳的主力)波長可能超過結構尺寸而被截斷。這是鰭寬進入 10 nm 量級後,聲子熱傳與塊材行為分道揚鑣的根本原因。