Kittel·固態物理

全書逐章 / 第 2 章 波繞射與倒晶格 / 小節 4

04布里淵區——sc/bcc/fcc 的倒晶格

BRILLOUIN ZONES — Reciprocal Lattice to sc / bcc / fcc Lattice

書本 p.33-38 · 小節 4 / 5 · 概念 12 推導 9 公式 14

把繞射條件 2k·G = G² 除以 4,它就從一條代數等式變成一句幾何陳述——k 的端點落在某條倒晶格向量的垂直平分面上;於是整個 Fourier 空間被這些平面切成布里淵區,而第一布里淵區就是倒晶格的 Wigner–Seitz 原胞。這一段接著把這個作圖法實際套到 sc、bcc、fcc 三種立方晶格上,算出立方體、菱形十二面體與截角八面體,並建立起「正空間原胞越小、倒空間布里淵區越大」這條後面每一章都要用的換算。

Brillouin 的重述:把一條代數條件換成一張幾何圖

第 2 章走到這裡,繞射條件其實已經講完三遍了。第一遍是 p.25 的 Bragg 定律 2d sinθ = nλ,第二遍是 p.31 的 Laue 條件 Δk = G(式 21),第三遍是 p.32 的 2k·G = G²(式 23)。三者互相等價,為什麼還要有第四遍?

Kittel 在 p.33 開宗明義給了理由:Brillouin 的版本是固態物理裡最被廣泛使用的那一個。原文是 “Brillouin gave the statement of the diffraction condition that is most widely used in solid state physics, which means in the description of electron energy band theory and of the elementary excitations of other kinds.” 注意這句話裡的用詞——「最廣泛使用」的理由不是它在 X 光繞射上比較方便(事實上並不),而是它在電子能帶理論與各類元激發(聲子、磁振子、電漿子……)的描述上不可取代。換句話說,這一段的價值不在第 2 章本身,而在它為第 4 章之後的整本書鋪了一個座標系統。

定義本身只有一句話:“A Brillouin zone is defined as a Wigner-Seitz primitive cell in the reciprocal lattice.” 布里淵區就是倒晶格的 Wigner–Seitz 原胞。書中並在括號裡提醒 “(The construction in the direct lattice was shown in Fig. 1.4.)”——同樣的作圖法你在第 1 章已經對正晶格做過一次了,現在只是把它原封不動搬到 k 空間。

Wigner–Seitz 作圖法的內容是:選一個格點當原點,畫出通往鄰近格點的連線,在每條連線的中點作垂直平分面,這些平面圍出的、包含原點的最小區域就是 WS 原胞。它的三個性質值得先記住:(i)它是原胞,體積等於一個格點所佔的體積;(ii)它有晶格的完整點群對稱性,這是平行六面體原胞做不到的;(iii)它裡面的每一點都「離原點比離任何其他格點近」。第三條在倒空間裡會直接翻譯成物理:布里淵區內部的 k 不滿足繞射條件,區邊界上的 k 才滿足。

這一段接下來做的事情非常樸素:先證明「繞射條件 = 落在垂直平分面上」(p.33),再把這個作圖法實際套到 sc、bcc、fcc 三種立方晶格上(p.34–p.38),算出立方體、菱形十二面體、截角八面體三個形狀。看似只是幾何練習,但這三個形狀之後會在第 7、8、9 章反覆出現,而算出它們體積的那個公式 (2π)³/V_c 會在第 5 章變成 Debye 模型的模數計數基礎。

tier 說明:引號內的句子、布里淵區的定義與 Fig. 1.4 的交叉引用逐字出自書本 p.33;Wigner-Seitz 作圖法三條性質的整理、以及最後一段關於後續章節的展望,是我補的教學說明。

Brillouin 形式的繞射條件

k · (½G) = (½G)²

(26)  p.33

k:晶體內入射波的波向量,由倒空間原點畫出。G:任一倒晶格向量,由原點指向某個倒晶格點。½G:該倒晶格向量的中點位置向量。左式是點積(純量),右式是 ½G 的長度平方。

出處推導
由式 (23) 2k·G = G² 兩邊同除以 4 得到。而 (23) 本身來自 (22) 2k·G + G² = 0(即彈性散射 k′ = k + G 配上 k′² = k²)在 G → −G 代換後的形式。整條鏈是:Bragg (1) ≡ Laue Δk = G (21) ≡ (22) ≡ (23) ≡ (26)。
極限行為
幾何意義:k 的端點必須落在「由原點到 G 這條線段的垂直平分面」上。取最短的 G 就給出離原點最近的一對平面,也就是第一布里淵區的邊界。由此可讀出一個重要極限——若 |k| < ½|G_min|,k 的端點連最近的那個平面都搆不到,晶體完全不可能產生繞射。這等價於 λ > 2d_max 時無繞射(因為 d = 2π/|G|,條件 2π/λ ≥ ½|G| 即 λ ≤ 4π/|G| = 2d)。可見光在晶體中不繞射、X 光才行,物理根源就在這一行。反過來,|k| 遠大於所有 ½|G| 時(高能電子繞射,k ~ 10² Å⁻¹)Ewald 球幾乎變成平面,同時掃到大量倒晶格點,這就是 TEM 電子繞射圖能一次出現整片斑點的原因。
典型數值
Si(diamond 結構,a = 5.430 Å,數值目視自書本 p.20 Table 3 影像):最短 G 為 ⟨111⟩ 型,|G| = √3·(2π/a) = 2.004 Å⁻¹,故 ½|G| = 1.002 Å⁻¹,對應 d(111) = 2π/|G| = 3.135 Å、λ_max = 2d = 6.27 Å。Cu Kα 的 λ = 1.5418 Å 遠小於此,可繞射;可見光 λ = 500 nm 對應 k = 1.26×10⁻³ Å⁻¹,比 ½|G_min| 小了約 800 倍,完全不可能繞射。

為什麼「除以 4」之後就變成垂直平分面:把 (23) → (26) 的幾何補完補推導

設定:書本 p.33 在 (23) 與 (26) 之間只有一句 “We divide both sides by 4 to obtain”,然後直接跳到「作一個垂直於 G 且過其中點的平面」。中間「為什麼除以 4 就會冒出垂直平分面」這一步書上沒有寫,但它是整個布里淵區概念的樞紐,值得完整攤開。

解讀:上面第 8、9 步得到的 |k| = |k − G| 是最值得記住的形式。它把一條物理條件(繞射)翻譯成一句純距離的話:k 的端點與原點、與倒晶格點 G 等距。而「與原點等距或更近」正是 Wigner–Seitz 原胞的定義。於是:

  • k 在第一布里淵區內部 ⟺ |k| < |k − G| 對所有 G ≠ 0 成立 ⟺ 不滿足任何繞射條件 ⟺ 波可以自由傳播。
  • k 落在第一布里淵區邊界 ⟺ 對某個最短的 G 有 |k| = |k − G| ⟺ 恰好滿足繞射條件 ⟺ 入射波與 k − G 這個散射波簡併,兩者必須混合。

第二條就是第 7 章「Origin of the Energy Gap」(書本 p.165)整節的出發點:兩個能量相同的平面波 e^{ikr} 與 e^{i(k−G)r} 在區邊界上簡併,晶體位能的 U_G 分量把簡併打開,於是在區邊界上開出能隙。所以「布里淵區邊界」與「能隙位置」是同一件事的兩個名字——這也是為什麼能帶圖一律畫在布里淵區上,而不是畫在任意一個平行六面體原胞上。

還有一個容易忽略的推論:(26) 對每一條 G 都定義一個平面,不只最短的那些。第二短的 G 給出第二組平面,圍出「第二布里淵區」,依此類推。第一區之所以特別,只是因為它是唯一一個完全被最短那批平面包住、且不含任何內部邊界的區域。書本 p.34 的正式定義用的正是這個措辭。

  1. 起點是書本 p.32 的式 (23):2k·G = G²。書上 p.33 只說 “We divide both sides by 4 to obtain (26)”,一步帶過,我們把中間展開。
  2. 左邊:2k·G ÷ 4 = k·G/2 = k·(½G)。這裡用到點積對純量的線性:k·(G/2) = (k·G)/2。
  3. 右邊:G² ÷ 4 = (G/2)·(G/2) = (½G)²。注意 G² 在 Kittel 的寫法裡就是 G·G = |G|²。
  4. 所以 (26) k·(½G) = (½G)² 與 (23) 是同一條式子,沒有加入任何新物理——新的只是「½G 被單獨拉出來當成一個幾何物件」。
  5. 現在讀幾何。令 n̂ = G/|G| 為 G 方向的單位向量。把 (26) 寫成 k·n̂·(|G|/2) = |G|²/4,兩邊除以 |G|/2 得 k·n̂ = |G|/2。
  6. k·n̂ 就是 k 在 G 方向上的投影長度。所以 (26) 說的是:k 沿 G 方向的投影恰好等於 |G| 的一半,與 k 垂直於 G 的那些分量完全無關。
  7. 「所有沿 n̂ 投影都等於 |G|/2 的點」構成一個平面,這個平面垂直於 G、並且通過點 ½G。而 ½G 正是線段 0→G 的中點。所以這個平面就是線段 0→G 的垂直平分面。
  8. 第二種等價讀法(更好用):計算 |k − G|² = k² − 2k·G + G²。把 (23) 的 2k·G = G² 代進去,後兩項相消,得 |k − G|² = k²,也就是 |k| = |k − G|。
  9. 這句話的意思是:k 的端點到原點的距離,等於它到倒晶格點 G 的距離。這正是 Wigner–Seitz 作圖法的定義性條件。

符號簿記:為什麼繞射後的波是 k − G,而不是 k + G補推導

設定:這是第 2 章最容易讓人卡住的一個符號問題。(21) 明明寫 Δk = G,p.33 卻寫 Δk = −G,而且兩處都沒有錯。Kittel 把符號代換的動作藏在 p.32 一句半行的話裡,讀者若跳過那句就會覺得 p.33 憑空翻了個號。

解讀:物理上根本不存在矛盾,因為倒晶格對原點中心對稱——G 是倒晶格向量時 −G 必然也是。真正的內容是:同一個繞射事件可以用 G 或 −G 兩套帳記,但 (23)/(26) 這個寫法選的是「讓 k·(½G) 為正」的那一套,也就是 k 與 G 夾銳角的那一套。這正好是 Fig. 9a 畫的樣子:k₁ 與 ½G_C 都指向右上方,夾角小於 90°。

實用上的記法:在 (26) 的約定下,永遠是「k 打在 G 的垂直平分面上,散射成 k − G」。幾何上 k 與 k − G 關於那個平面鏡射對稱(因為 |k| = |k − G|,而 k − (k−G) = G 垂直於該平面),這就是 Bragg「鏡面反射」圖像在倒空間裡的對應物。第 7 章寫簡併微擾時用的兩個基底 e^{ikr} 與 e^{i(k−G)r},取的也是這個約定,兩邊要對得起來。

  1. 書本 p.31 式 (19) 定義散射向量:k + Δk = k′,也就是 Δk ≡ k′ − k。
  2. 書本 p.31 式 (21) 給出繞射條件 Δk = G,代入 (19) 得 k′ = k + G。看起來繞射束方向應該是 k + G。
  3. 但書本 p.33 卻寫 “The diffracted beam will then be in the direction k − G, as we see from (19) with Δk = −G.” 出現了 Δk = −G。這不是筆誤,關鍵在 (22) → (23) 那一步做過一次符號代換。
  4. 回到 p.31 末。彈性散射要求 k′² = k²。把 k′ = k + G 代入:(k + G)² = k² → k² + 2k·G + G² = k² → 2k·G + G² = 0,這是式 (22)。
  5. Kittel 接著寫 “If G is a reciprocal lattice vector, so is −G, and with this substitution we can write (22) as (23)”。把 (22) 裡的 G 換成 −G:2k·(−G) + (−G)² = 0 → −2k·G + G² = 0 → 2k·G = G²,這就是 (23)。
  6. 所以 (23) 裡出現的那個 G,已經是原本 (22) 裡那個 G 的相反向量。既然原本滿足的是 Δk = G_舊 = −G_新,用 (23)/(26) 的語言講就是 Δk = −G。
  7. 代回 (19):k′ = k + Δk = k − G。與 p.33 的敘述一致。

Fig. 9a 逐筆讀圖(網頁不放圖,這裡用文字重建)

Fig. 9a 是一張倒空間的平面圖,背景是一片規則排列的黑點,代表原點 O 附近的倒晶格點(圖中排成方陣,五列五行)。圖上做了兩組完全平行的構造,用來示範同一件事做兩次:

第一組(垂直的那組)。 從原點 O 往右,同一列上最靠近的倒晶格點標為 C。連線 OC 就是倒晶格向量 G_C。從 O 沿 OC 畫一支短箭頭到中點,標記為 ½G_C。過這個中點作一條垂直於 OC 的白線,圖中在上端標上「1」——這就是平面 1。另有一支較粗的箭頭 k₁ 由 O 指向右上方,它的箭頭端恰好停在白線 1 上。圖說寫:“Any vector from the origin to the plane 1, such as k₁, will satisfy the diffraction condition k₁ · (½G_C) = (½G_C)².”

第二組(斜的那組)。 從 O 往左一行、往上兩列的那個倒晶格點標為 D。OD 就是 G_D,方向斜向左上。同樣畫出 ½G_D,並過中點作垂直於 OD 的白線,圖中在右上端標「2」——這是平面 2,它在圖上是一條由左下往右上傾斜的線。粗箭頭 k₂ 由 O 指向左上方,端點停在白線 2 上,對應條件 k₂ · (½G_D) = (½G_D)²。

這張圖要傳達的只有一句話:條件 (26) 裡的 G 每換一個,就換一個平面;k 的端點落在哪個平面上,就被哪個 G 繞射。圖中刻意讓 k₁ 與 k₂ 長度不同、方向不同,正是要強調滿足繞射的 k 不是孤立的幾個點,而是整整一張平面(三維中是平面,二維圖上畫成線)上的連續無窮多個。

一個讀圖細節值得指出:k₁ 與 ½G_C 的夾角是銳角,k₂ 與 ½G_D 的夾角也是銳角。這不是巧合,而是上一段符號簿記的結果——在 (26) 的約定下 k·(½G) = (½G)² > 0,所以 k 與 G 必然夾銳角。如果你畫出一支與 G 夾鈍角的 k,它滿足的是 −G 那一版的條件,對應的是另一片平面(原點另一側的那片)。

tier 說明:圖形內容與引號內的圖說逐字出自書本 p.34 的 Fig. 9a;最後一段關於銳角的讀圖推論是我補的。

Fig. 9b 與第一布里淵區的正式定義:從「一張平面」到「切碎整個 Fourier 空間」

Fig. 9a 一次只畫兩個 G。Fig. 9b 則把這件事做到底:對一個方形倒晶格,把原點連到附近所有倒晶格點的向量都畫出來(圖中是由中心放射出去的細黑線,帶箭頭,指向周圍一圈又一圈的黑點),再對每一條都作垂直平分線(圖中畫成白色粗線)。結果整張圖被白線切成密密麻麻的多邊形碎片——中央是一個正方形,外面一圈是四個三角形與四個更小的碎塊,再外面更碎。

Kittel 對這件事的描述是:“These planes divide the Fourier space of the crystal into fragments, as shown in Fig. 9b for a square lattice. The central square is a primitive cell of the reciprocal lattice. It is a Wigner-Seitz cell of the reciprocal lattice.”

接著是這一整段的核心定義,書中用斜體印出(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.

逐字拆這句話的三個限定詞,每一個都有作用:

  • smallest volume:不是「任一個被圍住的區域」,是最小的那個。第二區、第三區也都是被同一批平面圍住的,但它們更大、而且是環狀(不單連通)。
  • entirely enclosed(完全被圍住):這排除了那些只被部分平面切到、缺一面的開放碎塊。
  • perpendicular bisectors ... drawn from the origin:平分的是「由原點出發」的那些向量,所以整個構造以原點為中心,自動具備倒晶格的點群對稱性。

還有一個關鍵而容易被跳過的邏輯步驟:為什麼這個 WS 原胞是「原胞」(primitive cell)而不只是「某個區域」?因為 WS 作圖法的區域彼此全等、平移之後恰好鋪滿整個空間,而且每個區域內恰好含一個倒晶格點(就是它自己的中心)。一個含一格點又鋪滿空間的區域,體積必然等於任何一個平行六面體原胞的體積。這條性質稍後會被拿來證明 Problem 3 的 V_BZ = (2π)³/V_c——書中對此並未鋪陳,但它是那題「Hint: The volume of a Brillouin zone is equal to the volume of the primitive parallelepiped in Fourier space」這句提示背後的理由。

最後把 p.34 開頭那句在物理上最有用的話單獨抄出來:“A wave whose wavevector drawn from the origin terminates on any of these planes will satisfy the condition for diffraction.” 注意主詞是 A wave,不是「an x-ray」。這句話對聲子、對電子、對任何在週期位能中傳播的波都成立——這就是為什麼同樣一個布里淵區,第 4 章拿去畫聲子色散、第 7 章拿去畫能帶、第 9 章拿去畫費米面。

tier 說明:圖形內容、引號內句子與斜體定義逐字出自書本 p.34;三個限定詞的拆解、以及「為什麼 WS 原胞是原胞」那一段(含它與 Problem 3 提示的關聯)是我補的,書中未鋪陳。

p.34 的歷史註記:一個 X 光繞射不需要、能帶論不能沒有的概念

p.34 在 sc 小節開始前插了一段只有兩行的歷史註記,很多人讀過就忘,但它其實是在解釋「這一節為什麼放在第 2 章」:

Historically, Brillouin zones are not part of the language of x-ray diffraction analysis of crystal structures, but the zones are an essential part of the analysis of the electronic energy-band structure of crystals.

為什麼繞射學不需要它。做結構分析的人要回答的問題是「哪些反射會出現、強度多少」。這只需要一份倒晶格點的清單(hkl)加上結構因子 S_G(本章 p.39 之後),完全用不到「哪些 k 離哪個格點最近」這種分割空間的資訊。Ewald 球構造(Fig. 8,p.33)就是繞射學真正在用的幾何工具——它關心的是球面掃過哪些點,而不是空間被切成哪些區。

為什麼能帶論非要它不可。能帶論的出發點是 Bloch 定理:波函數可以標記為 ψ_k(r) = e^{ik·r}u_k(r),而 k 與 k + G 標記的是同一個狀態。於是 k 這個量子數有無窮多個等價選擇,必須挑一個代表元。第一布里淵區就是這個代表元的標準選法:它每個等價類恰好取一個 k(邊界上成對重複,但這是零測度集),而且它保有晶體的完整點群對稱性,所以能帶 ε_n(k) 在其中的對稱性一目了然。若改用平行六面體原胞當代表元區域,能帶圖會被切得支離破碎,看不出對稱點。

這就是第 7 章(ENERGY BANDS,書本 p.161)必須使用布里淵區的直接理由,而且第 7 章 p.180「Number of Orbitals in a Band」與 p.181「Metals and Insulators」兩節的全部結論——一條能帶恰好容納 2N 個電子、因此每原胞偶數個電子才可能是絕緣體——都建立在「第一布里淵區裡恰好有 N 個允許的 k」這件事上,而那件事又建立在本段算出的 V_BZ = (2π)³/V_c 上。第 9 章(FERMI SURFACES AND METALS,書本 p.221)開頭的 reduced zone scheme 與 periodic zone scheme 兩節,講的也完全是「怎麼把 k 折回第一區」。

tier 說明:上面「為什麼繞射學不需要/為什麼能帶論需要」的展開是我補的教學說明(Bloch 定理、代表元、2N 電子的論證出自後面章節),書本 p.34 只印了那兩行結論句。

Fig. 10 與 Fig. 11:二維斜晶格與一維鏈——作圖法的兩個最小例子

Fig. 10(p.35,二維斜晶格)。 圖上是一片斜的(oblique,非直角、非等長)二維倒晶格點陣,原點標 O,附近兩個格點標 A 與 B。由 O 畫出數支粗箭頭通往鄰近格點,每支箭頭的中點處畫一條細線垂直於它(圖中在垂足處畫了小方角記號標示垂直)。這些細線互相交割,中央圍出的最小封閉區域(圖中以淺灰底標示,形狀是一個不規則六邊形)就是第一布里淵區。圖說:“We first draw a number of vectors from O to nearby points in the reciprocal lattice. Next we construct lines perpendicular to these vectors at their midpoints. The smallest enclosed area is the first Brillouin zone.”

這張圖選斜晶格是有用意的:它說明布里淵區的形狀完全由晶格本身決定,不需要任何對稱性當前提。二維一般情形下 WS 原胞是六邊形(只有在特殊比例下退化成矩形),這一點在第 9 章畫二維費米面時會直接用到。

Fig. 11(p.35,一維)。 上半是一條「Linear crystal lattice」:一排等距的點,間距標為 a。下半是對應的「Reciprocal lattice」:另一排等距的點,基向量標為 b,圖說明確寫 “The basis vector in the reciprocal lattice is b, of length equal to 2π/a.” 由原點 O 出發的最短倒晶格向量只有兩條:b 與 −b。它們的垂直平分「面」在一維退化成兩個點,位置在 ±½·(2π/a) = ±π/a。圖中以深色區塊標出 −π/a 到 +π/a 這一段,並在下方標註 k = −π/a 與 k = π/a。圖說收尾:“The boundaries are at k = ±π/a.”

一維這個例子雖然簡單,卻是後面用得最兇的一個。第 4 章(PHONONS I)p.93 有一整節叫「First Brillouin Zone」,講的就是單原子線性鏈的聲子色散 ω(k) 在 k = ±π/a 處群速度歸零、以及 k 超出這個範圍後不帶新資訊(因為 k 與 k + 2π/a 描述完全相同的原子位移組態)。第 7 章的 Kronig–Penney 模型(p.168)與近自由電子模型的能隙也都開在 k = ±π/a。換句話說,Fig. 11 這條 −π/a 到 π/a 的線段,是整本書出現次數最多的一張圖。

tier 說明:兩張圖的圖形內容與引號內圖說逐字出自書本 p.35;「為什麼選斜晶格」的詮釋、以及與第 4、7、9 章的連結(節名與頁碼依書中目錄)是我補的。

sc 晶格的原始平移向量

a₁ = a x̂ ; a₂ = a ŷ ; a₃ = a ẑ

(27a)  p.34

a:簡單立方晶格的晶格常數(立方體邊長)。x̂, ŷ, ẑ:互相正交的單位向量,書中 p.35 明言 “Here x̂, ŷ, ẑ are orthogonal vectors of unit length.”

出處推導
簡單立方晶格的定義。sc 的慣用立方胞本身就是原胞(每胞一個格點),所以立方體的三條邊直接就是原始平移向量,不需要取半。
極限行為
這是三種立方晶格中唯一「慣用胞 = 原胞」的情形。bcc 與 fcc 的慣用立方胞分別含 2 與 4 個格點,因此它們的原始平移向量 (29)、(34) 都帶有係數 ½a,而且不再互相正交——這正是後面兩個小節的叉積算起來比較費工的原因。

sc 晶格的倒晶格原始平移向量

b₁ = (2π/a) x̂ ; b₂ = (2π/a) ŷ ; b₃ = (2π/a) ẑ

(27b)  p.35

b₁, b₂, b₃:倒晶格的原始平移向量,量綱為 [1/長度]。2π/a:倒晶格的晶格常數。

出處推導
由書本 p.29 式 (13) b₁ = 2π (a₂×a₃)/(a₁·a₂×a₃) 等三式,代入 (27a) 直接算出。書中 p.35 先寫 “The volume of the cell is a₁·a₂×a₃ = a³”,再說 “found from the standard prescription (13)”。
極限行為
結論是「sc 的倒晶格仍是 sc,只是晶格常數由 a 變成 2π/a」——書中原句 “Here the reciprocal lattice is itself a simple cubic lattice, now of lattice constant 2π/a.” 這是三種立方晶格裡唯一自倒(self-reciprocal,容許尺度改變)的一種。反比關係要記牢:a 變大 → 倒晶格變密 → 布里淵區變小。極限上,a → ∞(孤立原子、無週期性)時 2π/a → 0,布里淵區縮成一點,能帶塌成離散能階;a → 0(連續介質)時布里淵區擴張到無窮大,倒晶格點退到無窮遠,於是沒有任何 Bragg 反射、沒有能隙,回到自由電子/連續彈性介質。
典型數值
取 a = 4 Å(一般金屬量級),2π/a = 1.571 Å⁻¹;布里淵區邊界在 π/a = 0.785 Å⁻¹。相較之下 Cu Kα X 光的 k = 2π/1.5418 Å = 4.075 Å⁻¹,遠大於區邊界,所以能輕易觸及許多反射。

sc 的第一布里淵區邊界

±½b₁ = ±(π/a) x̂ ; ±½b₂ = ±(π/a) ŷ ; ±½b₃ = ±(π/a) ẑ

(28)  p.36

式子給的是六個中點位置向量;真正的邊界是「過這六點且分別垂直於 ±b₁, ±b₂, ±b₃ 的六個平面」,即 k_x = ±π/a、k_y = ±π/a、k_z = ±π/a。

出處推導
把 (27b) 的六條最短倒晶格向量 ±b₁, ±b₂, ±b₃ 代入 (26) 的垂直平分面構造。書中原句:“The boundaries of the first Brillouin zones are the planes normal to the six reciprocal lattice vectors ±b₁, ±b₂, ±b₃ at their midpoints”。
極限行為
書中緊接著寫 “The six planes bound a cube of edge 2π/a and of volume (2π/a)³; this cube is the first Brillouin zone of the sc crystal lattice.” 注意 sc 是三者中唯一「六個平面就把區域完全封閉、不需要更多平面」的情形——因為六條最短 G 的方向 ±x̂, ±ŷ, ±ẑ 已經張出一個封閉凸體。bcc 與 fcc 的最短 G 方向分別是 ⟨110⟩ 與 ⟨111⟩,它們單獨圍出的形狀在 fcc 情形下還會被次短的 G 截掉角(式 38),這是三者最關鍵的差別。
典型數值
a = 4 Å 時,BZ 為邊長 1.571 Å⁻¹ 的立方體,體積 (2π/a)³ = 3.876 Å⁻³。用 Problem 3 的公式檢查:(2π)³/V_c = 248.05/64 = 3.876 Å⁻³,完全吻合(V_c = a³ = 64 ų)。

sc 的 b 向量實算,以及第一個 V_BZ = (2π)³/V_c 的驗證補推導

設定:書本 p.35 只寫「由標準處方 (13) 求得」就直接給出 (27b),中間的叉積沒有寫。sc 這個例子算起來最輕鬆,正好拿來把整套流程走一遍、並第一次驗證 Problem 3 的公式,後面 bcc 與 fcc 就只是同一套流程配上比較麻煩的向量。

解讀:sc 的特殊性在於三個 a_i 互相正交,所以 b_i 與 a_i 方向相同、只是長度取倒數(再乘 2π)。一旦 a_i 不正交,b_i 就不再平行於任何 a_i——bcc 的 b₁ = (2π/a)(ŷ+ẑ) 與 a₁ = ½a(−x̂+ŷ+ẑ) 方向明顯不同,這是初學者最常誤解的一點。正確的敘述是書本 p.29 的 “Each vector defined by (13) is orthogonal to two axis vectors of the crystal lattice”:b₁ 垂直於 a₂ 與 a₃,而不是平行於 a₁。

  1. 由 (13):b₁ = 2π (a₂×a₃)/(a₁·a₂×a₃)。
  2. 分母(原胞體積):a₁·a₂×a₃ = (a x̂)·(a ŷ × a ẑ) = a³ (x̂·(ŷ×ẑ)) = a³·1 = a³。書中 p.35 直接寫了這個結果。
  3. 分子:a₂×a₃ = (a ŷ)×(a ẑ) = a² (ŷ×ẑ) = a² x̂。
  4. 所以 b₁ = 2π·a² x̂ / a³ = (2π/a) x̂。同理輪換得 b₂ = (2π/a) ŷ、b₃ = (2π/a) ẑ,即 (27b)。
  5. 檢查正交歸一關係 (14) b_i·a_j = 2πδ_ij:b₁·a₁ = (2π/a)(a)(x̂·x̂) = 2π ✓;b₁·a₂ = (2π/a)(a)(x̂·ŷ) = 0 ✓。
  6. 倒晶格原胞體積:b₁·b₂×b₃ = (2π/a)³ (x̂·(ŷ×ẑ)) = (2π/a)³ = 8π³/a³。
  7. Problem 3 的預言:(2π)³/V_c = 8π³/a³。兩者相同 ✓。
  8. 再用 WS 原胞驗一次:由 (28),BZ 是邊長 2π/a 的立方體,體積 (2π/a)³ = 8π³/a³,與平行六面體原胞體積相同 ✓。這具體地印證了「WS 原胞與平行六面體原胞體積相等」。

Problem 3 的一般證明:V_BZ = (2π)³/V_c補推導

設定:Problem 3(書本 p.44)要求證 V_BZ = (2π)³/V_c。這條式子在本段落裡會被 sc、bcc、fcc 三次驗證,也是第 5、6、7 章所有「態密度」「模數計數」的地基,因此值得把證明完整寫出來而不只是驗算特例。書中只在習題裡給提示,正文完全沒有證。

解讀:這條式子的物理內容比它的代數更重要。取一塊體積 V_total = N V_c 的晶體、施加週期性邊界條件,允許的 k 在倒空間中均勻分布,密度為 V_total/(2π)³ 個 k 每單位倒空間體積。於是第一布里淵區內的 k 數目為

N_k = [V_total/(2π)³] × V_BZ = [N V_c/(2π)³] × [(2π)³/V_c] = N

也就是說:第一布里淵區裡恰好有 N 個允許的 k,N 就是晶體中的原胞數。這個結果本身不出現在第 2 章(週期性邊界條件要到第 5 章才引入),但它是 (2π)³/V_c 這條式子真正的用途:

  • 第 4/5 章:每個原胞含 p 個原子 → 3p 條聲子分支 × N 個 k = 3pN 個簡正模,恰好等於自由度總數。Debye 模型就是用一個等體積的球(Debye 球)去取代布里淵區,好把模數算對。
  • 第 6/7 章:一條能帶有 N 個 k、每個 k 兩個自旋態 → 一條能帶容納 2N 個電子。第 7 章 p.180「Number of Orbitals in a Band」與 p.181「Metals and Insulators」的全部推論由此而來:每原胞奇數個價電子必為金屬。
  • 第 9 章:費米球體積佔布里淵區體積的比例 = 每原胞價電子數 / 2。這是判讀費米面圖的第一個檢查。
  1. 目標:證明倒晶格原胞體積 |b₁·b₂×b₃| = (2π)³/V_c,其中 V_c = |a₁·a₂×a₃|。書本 p.44 Problem 3 給的提示是 “The volume of a Brillouin zone is equal to the volume of the primitive parallelepiped in Fourier space. Recall the vector identity (c×a)×(a×b) = (c·a×b)a.”
  2. 第一步(提示的前半):為什麼可以用平行六面體體積代替 BZ 體積?因為布里淵區是倒晶格的 WS 原胞,而同一個晶格的所有原胞體積相同——每個原胞恰含一個格點,且平移複本鋪滿全空間,故體積 = 1/(格點數密度)。這一步書中沒有展開,但它是提示成立的理由。
  3. 第二步:把 (13) 代入。令 V = a₁·a₂×a₃。則 b₁ = 2π(a₂×a₃)/V、b₂ = 2π(a₃×a₁)/V、b₃ = 2π(a₁×a₂)/V。
  4. 第三步:先算 b₂×b₃ = (2π/V)² (a₃×a₁)×(a₁×a₂)。
  5. 第四步:套用提示的恆等式 (c×a)×(a×b) = (c·a×b)a,對應 c → a₃、a → a₁、b → a₂:(a₃×a₁)×(a₁×a₂) = (a₃·a₁×a₂) a₁。
  6. 第五步:純量三重積在輪換下不變,a₃·a₁×a₂ = a₁·a₂×a₃ = V。所以 b₂×b₃ = (2π)² V a₁ / V² = (2π)² a₁ / V。
  7. 第六步:b₁·(b₂×b₃) = [2π(a₂×a₃)/V] · [(2π)² a₁/V] = (2π)³ (a₁·a₂×a₃)/V² = (2π)³ V/V² = (2π)³/V。
  8. 得證:V_BZ = |b₁·b₂×b₃| = (2π)³/V_c。∎
  9. 順帶一個副產品:第五步的 b₂×b₃ = (2π)²a₁/V 說明倒晶格的倒晶格就是原晶格。事實上由 (13) 套用到 b 上:2π(b₂×b₃)/(b₁·b₂×b₃) = 2π·[(2π)²a₁/V] / [(2π)³/V] = a₁ ✓。

bcc 晶格的原始平移向量

a₁ = ½a(−x̂ + ŷ + ẑ) ; a₂ = ½a(x̂ − ŷ + ẑ) ; a₃ = ½a(x̂ + ŷ − ẑ)

(29)  p.36

a:慣用立方胞的邊長,書中原句 “where a is the side of the conventional cube and x̂, ŷ, ẑ are orthogonal unit vectors parallel to the cube edges”。三支向量都沿體對角線方向、長度為體對角線的一半(½a√3),由一個體心指向三個相鄰立方胞的體心。

出處推導
bcc 的定義(立方體八個角 + 一個體心)。取體心為原點,三支向量分別指向三個相鄰立方體的體心,即 Fig. 12 所畫。這三支向量的長度都是 ½a√3(體對角線的一半),彼此夾角的餘弦為 −⅓。
極限行為
這組向量不是唯一的原始平移向量選擇——任何行列式為 ±1 的整數線性組合都同樣有效,例如常見的另一組 ½a(x̂+ŷ+ẑ)、a x̂、a ŷ。Kittel 挑這一組是為了讓 (31) 算出來的 b 恰好長得像 (34) 的 fcc 形式,方便一眼看出「bcc 的倒晶格是 fcc」。換了別組 a_i,算出的 b_i 會是 (31) 的整數線性組合,張出同一個倒晶格,但式子不會那麼漂亮。這也提醒:原胞的選法不唯一,但晶格、倒晶格、布里淵區都唯一。

bcc 原胞體積

V = |a₁ · a₂ × a₃| = ½a³

(30)  p.36

V:bcc 原胞(含一個格點)的體積。a³:慣用立方胞的體積。

出處推導
由 (29) 直接計算純量三重積,等價於「慣用立方胞含 2 個格點,故每格點佔 a³/2」。
極限行為
這個 ½ 是 bcc 一切數字的來源:原胞體積 ½a³ → 倒晶格原胞體積 (2π)³/(½a³) = 2(2π/a)³ → 布里淵區(菱形十二面體)體積也是 2(2π/a)³。若寫成慣用立方胞的語言,倒空間慣用立方胞邊長 4π/a、體積 (4π/a)³ = 8(2π/a)³ 含 4 個倒晶格點(fcc),每點 2(2π/a)³ ✓ 自洽。
典型數值
Na(bcc,a = 4.225 Å @ 5 K,目視自書本 p.20 Table 3 影像):V_c = ½a³ = ½(75.42) = 37.71 ų;原子數密度 1/V_c = 2.652×10²² cm⁻³。W(bcc,a = 3.16 Å):V_c = 15.78 ų。

實算 bcc 的 ½a³ 與三個叉積:由 (29) 走到 (31)補推導

設定:書本 p.36 在 (30) 與 (31) 之間只寫了一句 “The primitive translations of the reciprocal lattice are defined by (13). We have, using (28),”,三個叉積一個都沒算。這裡把它們全部算出來,順便驗證正交關係與倒晶格原胞體積。

校註(p.36):書中那句 “We have, using (28),” 的式號有問題。(28) 是上一小節 sc 的布里淵區邊界 ±½b₁ = ±(π/a)x̂,與 bcc 的 b 向量計算毫無關係;真正用到的是 (13)(定義)與 (30)(體積 ½a³ 當分母)。依推導流程判斷,此處應為 (30) 的誤植。此點經 PyMuPDF 將該頁算繪為影像目視核對,印刷版確實印的是「(28)」,因此是原書排印問題,不是文字層抽取誤差。

解讀:算完之後把 (31) 與 (34) 並排看,關係就浮出來了。(34) 的 fcc 原始平移向量是 ½A(ŷ+ẑ)、½A(x̂+ẑ)、½A(x̂+ŷ),其中 A 是 fcc 的立方邊長。而 (31) 可以寫成 ½·(4π/a)(ŷ+ẑ) 等三式。兩者對照,bcc(立方邊 a)的倒晶格是一個立方邊 A = 4π/a 的 fcc 晶格。書中 p.36 的說法是 “Note by comparison with Fig. 14 (p. 37) that these are just the primitive vectors of an fcc lattice, so that an fcc lattice is the reciprocal lattice of the bcc lattice.”——它要你翻到下一頁去比對圖,這裡直接把 4π/a 這個尺度寫清楚,省去讀者自己對。

  1. 先算原胞體積(式 30)。把 (29) 寫成矩陣形式,係數 ½a 提出來:V = (a/2)³ · det[[−1,1,1],[1,−1,1],[1,1,−1]]。
  2. 展開行列式:det = (−1)·[(−1)(−1) − (1)(1)] − (1)·[(1)(−1) − (1)(1)] + (1)·[(1)(1) − (−1)(1)] = (−1)(1−1) − (1)(−1−1) + (1)(1+1) = 0 + 2 + 2 = 4。
  3. 所以 V = (a³/8)·4 = ½a³,即 (30)。(物理檢查:bcc 慣用立方胞含 8×⅛ + 1 = 2 個格點,每格點 a³/2 ✓)
  4. 接著算 b₁ = 2π(a₂×a₃)/V。先算方括號內的叉積,把 ½a 提出:a₂×a₃ = (a/2)² [(1,−1,1) × (1,1,−1)]。
  5. (1,−1,1)×(1,1,−1) 的三個分量:x 分量 = (−1)(−1) − (1)(1) = 0;y 分量 = (1)(1) − (1)(−1) = 2;z 分量 = (1)(1) − (−1)(1) = 2。得 (0, 2, 2) = 2(ŷ + ẑ)。
  6. 所以 a₂×a₃ = (a²/4)·2(ŷ+ẑ) = (a²/2)(ŷ+ẑ),於是 b₁ = 2π·(a²/2)(ŷ+ẑ) / (a³/2) = (2π/a)(ŷ + ẑ) ✓
  7. 同理 a₃×a₁ = (a²/4)[(1,1,−1)×(−1,1,1)]。分量:x = (1)(1) − (−1)(1) = 2;y = (−1)(−1) − (1)(1) = 0;z = (1)(1) − (1)(−1) = 2。得 2(x̂+ẑ),故 b₂ = (2π/a)(x̂ + ẑ) ✓
  8. 最後 a₁×a₂ = (a²/4)[(−1,1,1)×(1,−1,1)]。分量:x = (1)(1) − (1)(−1) = 2;y = (1)(1) − (−1)(1) = 2;z = (−1)(−1) − (1)(1) = 0。得 2(x̂+ŷ),故 b₃ = (2π/a)(x̂ + ŷ) ✓
  9. 驗證 (14):b₁·a₁ = (2π/a)(a/2)[(0,1,1)·(−1,1,1)] = π(0+1+1) = 2π ✓;b₁·a₂ = π[(0,1,1)·(1,−1,1)] = π(0−1+1) = 0 ✓;b₁·a₃ = π[(0,1,1)·(1,1,−1)] = π(0+1−1) = 0 ✓
  10. 倒晶格原胞體積:b₁·b₂×b₃ = (2π/a)³ det[[0,1,1],[1,0,1],[1,1,0]] = (2π/a)³ · 2 = 2(2π/a)³,與書中 p.37 的敘述一致。
  11. Problem 3 檢查:(2π)³/V_c = (2π)³/(½a³) = 2(2π)³/a³ = 2(2π/a)³ ✓

bcc 的倒晶格原始平移向量

b₁ = (2π/a)(ŷ + ẑ) ; b₂ = (2π/a)(x̂ + ẑ) ; b₃ = (2π/a)(x̂ + ŷ)

(31)  p.36

a:bcc 正晶格慣用立方胞的邊長。三支 b 的長度都是 √2(2π/a),彼此夾角 60°。

出處推導
由 (13) 代入 (29),分母用 (30) 的 V = ½a³。三個叉積的完整計算見上一個推導區塊。
極限行為
把 (31) 改寫成 ½(4π/a)(ŷ+ẑ) 等,立刻看出這是立方邊 4π/a 的 fcc 晶格的原始平移向量。所以:bcc ↔ fcc 互為倒晶格,且尺度換算是 a ↔ 4π/a(不是 2π/a!這個因子 2 常被記錯,根源是 bcc/fcc 的慣用胞不是原胞)。檢查:a → ∞ 時 4π/a → 0,倒晶格塌縮;a → 0 時倒晶格無限膨脹。
典型數值
Na(a = 4.225 Å):倒晶格為 fcc,立方邊 4π/a = 2.974 Å⁻¹;|b| = √2(2π/a) = 2.103 Å⁻¹。W(a = 3.16 Å):4π/a = 3.977 Å⁻¹。

bcc 的一般倒晶格向量

G = v₁b₁ + v₂b₂ + v₃b₃ = (2π/a)[(v₂+v₃) x̂ + (v₁+v₃) ŷ + (v₁+v₂) ẑ]

(32)  p.36

v₁, v₂, v₃:任意整數(書中原句 “for integral v₁, v₂, v₃”)。括號內三個係數 (v₂+v₃)、(v₁+v₃)、(v₁+v₂) 就是 G 在立方座標系下以 2π/a 為單位的分量。

出處推導
把 (31) 的三支 b 代入 (15) G = v₁b₁ + v₂b₂ + v₃b₃ 並合併同類項。例如 x̂ 分量只來自 b₂ 與 b₃,故係數為 v₂ + v₃。
極限行為
這條式子的價值不在它本身,而在它的奇偶結構:三個係數之和 (v₂+v₃)+(v₁+v₃)+(v₁+v₂) = 2(v₁+v₂+v₃),恆為偶數。也就是說 bcc 的倒晶格恰好是「以 2π/a 為單位、三個整數分量之和為偶數」的那些點——而這正是立方邊 4π/a 的 fcc 晶格的座標描述(見下一個推導)。同時這也與 p.40 的結構因子消光規則 “S = 0 when v₁+v₂+v₃ = odd integer” 是同一件事的兩種說法。
典型數值
取 (v₁,v₂,v₃) = (1,0,0) → G = (2π/a)(0,1,1),係數和 2(偶)✓;取 (1,1,0) → G = (2π/a)(1,1,2),和 4 ✓;取 (1,1,1) → G = (2π/a)(2,2,2),和 6 ✓。永遠湊不出 (2π/a)(1,0,0) 這種和為奇數的點。

嚴格證明「bcc 的倒晶格就是 fcc」:奇偶論證的雙向版本補推導

設定:書本 p.36 證明「bcc 的倒晶格是 fcc」的方式是 “Note by comparison with Fig. 14 (p. 37)”——請讀者翻頁看圖比對形狀。這對第一次讀的人其實不太夠:比對的是「三支向量長得一樣」,但嚴格說要證的是「兩個點集合相同」。這裡補上代數證明。

解讀:這個奇偶論證有兩個延伸用途,都很實用。

第一,它直接解釋了 p.40 的 bcc 消光規則。用慣用立方胞(含 2 原子,位於 000 與 ½½½)做結構因子,書中得到 “S = 0 when v₁+v₂+v₃ = odd integer ; S = 2f when v₁+v₂+v₃ = even integer”。所謂「消光」,在倒晶格語言裡就是:以簡單立方倒晶格 (2π/a)(h,k,l) 去標記所有可能反射時,那些 h+k+l 為奇的點根本不是 bcc Bravais 晶格的倒晶格點。兩種說法完全等價——結構因子消光是「用了過大的慣用胞」所付出的代價,而不是什麼額外的物理。書中 p.40 也用實例印證:“Metallic sodium has a bcc structure. The diffraction pattern does not contain lines such as (100), (300), (111), or (221), but lines such as (200), (110), and (222) will be present”——逐一檢查指標和:100→1 奇 ✗、300→3 奇 ✗、111→3 奇 ✗、221→5 奇 ✗;200→2 偶 ✓、110→2 偶 ✓、222→6 偶 ✓。完全吻合。

第二,同樣的論證反過來用在 fcc 上:fcc 的倒晶格是「三個整數座標同奇偶」的集合,也就是 bcc。這對應 p.41 的 “in the fcc lattice no reflections can occur for which the indices are partly even and partly odd”。

tier 說明:上面的代數證明與奇偶等價論證是我補的(書中只做目視比對);引號內的句子與 (100)/(200) 等實例則逐字取自書本 p.36、p.40、p.41。

  1. 把倒晶格點以 2π/a 為單位寫成整數三元組 (n₁, n₂, n₃),即 G = (2π/a)(n₁x̂ + n₂ŷ + n₃ẑ)。
  2. 方向一(由 (32) 推出偶和):由 (32),n₁ = v₂+v₃、n₂ = v₁+v₃、n₃ = v₁+v₂。三者相加 n₁+n₂+n₃ = 2(v₁+v₂+v₃) = 偶數。所以每個 bcc 倒晶格點的整數座標和必為偶。
  3. 方向二(由偶和反推得到):給定任意整數三元組 (n₁,n₂,n₃) 且 n₁+n₂+n₃ 為偶,解 v₁ = (−n₁+n₂+n₃)/2、v₂ = (n₁−n₂+n₃)/2、v₃ = (n₁+n₂−n₃)/2。
  4. 驗算:v₂+v₃ = [(n₁−n₂+n₃) + (n₁+n₂−n₃)]/2 = n₁ ✓,另兩式同理。
  5. 這三個 v 是整數嗎?−n₁+n₂+n₃ = (n₁+n₂+n₃) − 2n₁,偶數減偶數 = 偶數,除以 2 得整數 ✓。另兩式同理。
  6. 所以「bcc 的倒晶格」與「整數座標和為偶的立方格點集合」是完全相同的集合,不多不少。
  7. 最後確認這個集合就是 fcc。以 2π/a 為單位、立方邊長取 2(即實際邊長 4π/a)的 fcc 晶格,其格點為 (0,0,0)、(1,1,0)、(1,0,1)、(0,1,1) 加上 (2,0,0) 型平移。列舉可見這四個基座標的和分別是 0、2、2、2,皆偶;平移量 (2,0,0) 之和也是偶,故 fcc 全部格點座標和為偶。
  8. 反之,任一座標和為偶的整數三元組:若三個都偶,模 2 同餘 (0,0,0) → 是角點的平移;若恰兩個奇,模 2 同餘 (1,1,0) 之一 → 是面心點的平移。沒有第三種可能(三個都奇的和為奇,被排除)。故兩集合相等 ∎
  9. 因此 bcc(立方邊 a)的倒晶格 = fcc(立方邊 4π/a),與書中 p.36 由「比對 Fig. 14」得到的結論一致,但這裡是代數證明而非目視比對。

bcc 倒晶格中最短的 12 條 G

(2π/a)(±ŷ ± ẑ) ; (2π/a)(±x̂ ± ẑ) ; (2π/a)(±x̂ ± ŷ)

(33)  p.37

每一組給出 4 條向量(兩個正負號獨立),三組共 12 條。書中特別註明 “where all choices of sign are independent”。每條的長度都是 √2(2π/a)。

出處推導
由 (31) 張出的倒晶格(立方邊 4π/a 的 fcc)中,最近鄰的 12 個格點。在 fcc 中最近鄰距離是立方邊除以 √2,即 (4π/a)/√2 = √2(2π/a) ✓,且最近鄰數(配位數)恰為 12 ✓。
極限行為
12 這個數字直接決定布里淵區的面數:每條最短 G 貢獻一個邊界面,於是 bcc 的第一布里淵區有 12 個面,形狀是菱形十二面體(Fig. 13)。次短的 G 是 (4π/a)(±x̂) 型共 6 條,長度 2(2π/a) > √2(2π/a);它們的垂直平分面距原點 2π/a,恰好切在菱形十二面體的六個四價頂點上而不切入,所以不會截角——這是 bcc 與 fcc 最重要的差別(fcc 那邊就真的被截了,見式 38)。
典型數值
Na(a = 4.225 Å):|G_min| = √2(2π/a) = 2.103 Å⁻¹,對應 d = 2π/|G| = 2.988 Å = a/√2 ✓(即 (110) 面間距)。BZ 面心(習稱 N 點)距原點 1.052 Å⁻¹。

菱形十二面體是怎麼長出來的:從 12 條最短 G 到頂點座標與體積補推導

設定:書本 p.37 對 bcc 布里淵區只寫了結論:“This zone (for the bcc lattice) is bounded by the planes normal to the 12 vectors of Eq. (33) at their midpoints. The zone is a regular 12-faced solid, a rhombic dodecahedron, as shown in Fig. 13.” 為什麼是 12 個面、頂點在哪裡、體積多少、以及為什麼次短的 G 不會再把它截掉——全都沒寫。這裡補齊。

解讀:整段最值得帶走的觀念是——布里淵區的形狀完全由「最短的那批 G 有幾條、指向哪裡」決定。bcc 的倒晶格是 fcc,fcc 的配位數是 12,所以 BZ 有 12 個面;而次短殼層的 G 太長(2 vs √2,比值 1.414),它們的平分面退到 2π/a 之外,剛好與頂點相切,對形狀毫無貢獻。fcc 的情況(下一節)就不同:最短的 8 條與次短的 6 條長度比是 √3 : 2 = 1.732 : 2,差距只有 15%,於是次短那批真的切進來,把八面體的角削掉。

另外,兩個頂點類型(6 個四價、8 個三價)在能帶論裡有名字:四價頂點是 H 點(⟨100⟩ 方向,距 Γ 為 2π/a),三價頂點是 P 點(⟨111⟩ 方向,距 Γ 為 √3π/a),面心是 N 點(⟨110⟩ 方向,距 Γ 為 √2π/a)。這些字母標記不出現在第 2 章,但你在第 7、9 章看到的 bcc 金屬能帶圖(Γ–H–N–P–Γ)走的就是這幾個點之間的路徑。tier 說明:Γ/H/N/P 的字母命名是標準慣例,書本 p.33–38 並未使用,屬我補充。

  1. 由 (26),每條 G 貢獻一個平面。取 G = (2π/a)(x̂+ŷ),則 ½G = (π/a)(1,1,0),平面方程式為 k·(π/a)(1,1,0) = (π/a)²·2,化簡得 k_x + k_y = 2π/a。其餘 11 條輪換與變號,共得 12 個平面。
  2. 這 12 個平面圍出的凸體,頂點就是三個或以上平面的交點。先找四面交點:點 (2π/a, 0, 0) 同時滿足 k_x+k_y = 2π/a、k_x−k_y = 2π/a、k_x+k_z = 2π/a、k_x−k_z = 2π/a,共 4 個平面 → 這是一個四價頂點。這類頂點沿 ⟨100⟩ 有 6 個,距原點 2π/a。
  3. 再找三面交點:點 (π/a, π/a, π/a) 滿足 k_x+k_y = 2π/a、k_x+k_z = 2π/a、k_y+k_z = 2π/a,共 3 個平面 → 三價頂點。這類頂點沿 ⟨111⟩ 有 8 個,距原點 √3π/a ≈ 1.732π/a。
  4. 頂點總數 6 + 8 = 14,面數 12,由 Euler 公式 V − E + F = 2 得棱數 E = 12 + 14 − 2 = 24。這正是菱形十二面體的規格(12 面 / 24 棱 / 14 頂點)✓
  5. 檢查每個面確實是菱形:面 k_x+k_y = 2π/a 上的四個頂點是 (2π/a,0,0)、(π/a,π/a,π/a)、(0,2π/a,0)、(π/a,π/a,−π/a)。對角線長度分別為 |(2π/a,0,0)−(0,2π/a,0)| = 2√2π/a 與 |(π/a,π/a,π/a)−(π/a,π/a,−π/a)| = 2π/a,比值 √2 —— 這就是菱形十二面體的標準面形 ✓
  6. 體積算法一(切割法):以八個 ⟨111⟩ 頂點 (±π/a,±π/a,±π/a) 為角,內含一個邊長 2π/a 的立方體,體積 (2π/a)³ = 8π³/a³。再在立方體六個面上各加一個四角錐,錐底是立方體的面(面積 (2π/a)² = 4π²/a²)、錐頂是 ⟨100⟩ 頂點(距該面 π/a)。每錐體積 = ⅓·(4π²/a²)·(π/a) = (4/3)π³/a³,六錐共 8π³/a³。總體積 = 8π³/a³ + 8π³/a³ = 16π³/a³。
  7. 體積算法二(角錐分解法,交叉檢查):V = ⅓ Σ(面積 × 面到原點距離)。每個菱形面積 = ½·d₁·d₂ = ½(2√2π/a)(2π/a) = 2√2π²/a²,面到原點距離 = ½|G| = √2π/a。於是 V = ⅓ × 12 × (2√2π²/a²) × (√2π/a) = ⅓ × 12 × 4π³/a³ = 16π³/a³ ✓ 兩法一致。
  8. 與 Problem 3 對照:16π³/a³ = 2 × 8π³/a³ = 2(2π/a)³,正是書中 p.37 給的倒晶格平行六面體體積;而 (2π)³/V_c = 8π³/(½a³) = 16π³/a³ ✓ 三者完全吻合。
  9. 次短 G 不截角的驗證:次短 G = (4π/a)x̂,其垂直平分面為 k_x = 2π/a。菱形十二面體上 k_x 的最大值出現在頂點 (2π/a,0,0),恰為 2π/a。所以該平面與多面體只相切於一點,沒有切掉任何體積 ✓

p.37 的「八分之一」計數:為什麼平行六面體原胞只含一個倒晶格點

p.37 在介紹菱形十二面體之前,先花了四句話做一件看似瑣碎的事——確認由 b₁, b₂, b₃ 張出的平行六面體真的是原胞:

One primitive cell of the reciprocal lattice is the parallelepiped described by the b₁, b₂, b₃ defined by (31). The volume of this cell in reciprocal space is b₁·b₂×b₃ = 2(2π/a)³. The cell contains one reciprocal lattice point, because each of the eight corner points is shared among eight parallelepipeds. Each parallelepiped contains one-eighth of each of eight corner points (see Fig. 12).

計數邏輯是:平行六面體有 8 個角,每個角上的格點同時屬於 8 個相鄰的平行六面體,所以每個平行六面體分到 8 × ⅛ = 1 個格點。內部與面上、棱上都沒有其他格點(因為 b₁,b₂,b₃ 是原始平移向量,依定義如此)。

接著書中把兩種原胞並列:

Another primitive cell is the central (Wigner-Seitz) cell of the reciprocal lattice which is the first Brillouin zone. Each such cell contains one lattice point at the central point of the cell.

這兩句話並列的用意,是 Problem 3 那句提示的伏筆。平行六面體原胞容易算體積(純量三重積)卻沒有對稱性;WS 原胞(= 布里淵區)有完整對稱性卻不好直接積分。兩者都「每胞一個格點、平移後鋪滿空間」,因此體積必然相同。於是要算布里淵區體積時,可以偷懶去算平行六面體——這就是 p.44 Problem 3 提示 “The volume of a Brillouin zone is equal to the volume of the primitive parallelepiped in Fourier space” 的全部內容。書中沒有把這個邏輯明說,但它就藏在這兩段的並列裡。

一個讀圖上的小困擾要先講明:書中在「八分之一」那句後面寫 (see Fig. 12),而 Fig. 12(p.36)畫的是 bcc 正晶格的原始基矢,不是倒空間的平行六面體。它能當作「角點被八個胞共用」的視覺示意(圖中確實畫了多個相鄰立方胞、角點共用),但讀者若期待看到一個倒空間的平行六面體會撲空。這是交叉引用選得不太理想,不是內容錯誤。

tier 說明:兩段引號內的文字逐字出自書本 p.37;「兩種原胞並列是 Problem 3 提示的伏筆」這個解讀,以及關於 Fig. 12 交叉引用的提醒,是我補的。

Fig. 13 讀圖:菱形十二面體長什麼樣子我的判讀

Fig. 13(p.36 右上)畫的是 bcc 的第一布里淵區。圖說只有兩句:“First Brillouin zone of the body-centered cubic lattice. The figure is a regular rhombic dodecahedron.” 以下是把圖形內容轉成文字(本區塊標 mine,因為是我對圖的描述與補充,不是書中文字)。

圖中用細線暗示了一個邊長 4π/a 的立方框架,但立方體本身的 12 條稜並沒有畫出來;畫出來的是 12 段「把同一個面上兩條對邊的中點連起來、交於該面面心」的細線(6 個面各 2 段)。12 顆外圍小球就落在這個立方體 12 條稜的中點上,也就是 (2π/a)(±x̂±ŷ)、(2π/a)(±x̂±ẑ)、(2π/a)(±ŷ±ẑ) 這 12 個最近鄰倒晶格點(式 33),每一顆對應菱形十二面體的一個面。框架正中央還有一顆球,那是原點 Γ,即這個布里淵區所屬的格點。「中心一顆 + 12 條稜中點各一顆」正是把 fcc 倒晶格的慣用立方胞以格點為中心擺放時的樣子——立方體的 8 個角 (±2π/a,±2π/a,±2π/a) 與 6 個面心 (±2π/a,0,0) 三個座標之和為奇數,都不是 bcc 的倒晶格點。

中央的多面體是一個帶有明顯稜角的立體,可以辨認出幾個菱形面(不是正方形,也不是六邊形):每個面是四邊形,四個角分別是兩個「鈍角」與兩個「銳角」,相鄰面之間的稜線在頂點處三條或四條交會。整體輪廓像是把一個立方體的十二條稜「鼓起來」成十二個菱形面。

配合上一段算出的數據,這個形狀的規格是:

項目規格
面數12(菱形),法線方向 ⟨110⟩
面到中心距離√2π/a,即最短 G 長度的一半
頂點14 個:6 個在 ⟨100⟩ 上距中心 2π/a(四面交會)、8 個在 ⟨111⟩ 上距中心 √3π/a(三面交會)
稜數24
每面對角線比1 : √2
體積16π³/a³ = 2(2π/a)³

值得注意的是這個多面體剛好內接於那個立方框架:它的 6 個 ⟨100⟩ 頂點伸到 k = 2π/a,恰好等於框架的半邊長(框架邊長 4π/a,即倒晶格慣用立方胞的尺寸),所以那六個頂點恰好頂到框架的六個面心;另外 8 個 ⟨111⟩ 頂點距 Γ 只有 √3π/a ≈ 1.73π/a,落在立方體內部,沒有碰到任何一個角。這個細節在下一節 fcc 的情形會變成關鍵——那裡未截角的正八面體在 ⟨100⟩ 方向伸到 3π/a,超過同樣是 2π/a 的半邊長,於是真的被切掉六個角;這裡則剛好相切、不切入。

fcc 晶格的原始平移向量

a₁ = ½a(ŷ + ẑ) ; a₂ = ½a(x̂ + ẑ) ; a₃ = ½a(x̂ + ŷ)

(34)  p.37

a:fcc 慣用立方胞的邊長。三支向量由立方體某一角指向三個相鄰的面心,長度都是 a/√2,彼此夾角 60°。書中以 Fig. 14 標示。

出處推導
fcc 的定義(八個角 + 六個面心)。以一個角為原點,三支最短的格點連線就是通往三個共用該角的面心。
極限行為
把這組 a 與 bcc 的 (31) 的 b 對照即可看出互倒關係。另外注意 (34) 與 (31) 形式完全相同(都是 (ŷ+ẑ)、(x̂+ẑ)、(x̂+ŷ) 型),只是前面的係數從 ½a 換成 2π/a——這不是巧合,而是 bcc↔fcc 互倒的代數體現:½a ↔ ½(4π/a) = 2π/a。

fcc 原胞體積

V = |a₁ · a₂ × a₃| = ¼a³

(35)  p.37

V:fcc 原胞(含一個格點)的體積。

出處推導
由 (34) 計算純量三重積;等價於「慣用立方胞含 8×⅛ + 6×½ = 4 個格點,故每格點 a³/4」。
極限行為
這個 ¼ 決定了 fcc 一切倒空間數字:倒晶格原胞體積 = (2π)³/(¼a³) = 4(2π/a)³,布里淵區(截角八面體)體積也是 4(2π/a)³。與 bcc 的 2(2π/a)³ 相差正好 2 倍,因為正空間原胞體積差 2 倍。
典型數值
Cu(fcc,a = 3.61 Å,目視自書本 p.20 Table 3 影像):V_c = ¼a³ = 11.76 ų,原子密度 8.50×10²² cm⁻³。Si(diamond 結構,晶格為 fcc,a = 5.430 Å):V_c = ¼a³ = 40.03 ų,但每原胞含 2 個 Si 原子,故原子密度 = 2/40.03 ų = 5.00×10²² cm⁻³ —— 這個「每原胞 2 原子」在第 4 章會生出光學支聲子,在第 5 章會影響 Debye 模數計數,在 Cth(熱容)的估算裡是必須先弄對的第一個數字。

fcc 的倒晶格原始平移向量

b₁ = (2π/a)(−x̂ + ŷ + ẑ) ; b₂ = (2π/a)(x̂ − ŷ + ẑ) ; b₃ = (2π/a)(x̂ + ŷ − ẑ)

(36)  p.38

a:fcc 正晶格慣用立方胞邊長。三支 b 的長度都是 √3(2π/a)。

出處推導
由 (13) 代入 (34),分母用 (35) 的 V = ¼a³。完整叉積見下一個推導區塊。
極限行為
書中結論 “These are primitive translation vectors of a bcc lattice, so that the bcc lattice is reciprocal to the fcc lattice.” 把 (36) 寫成 ½(4π/a)(−x̂+ŷ+ẑ) 即可看出是立方邊 4π/a 的 bcc 晶格的原始平移向量,與 (29) 形式相同。連同上一節,兩個結論合起來構成一個完整的對偶:bcc ↔ fcc,且尺度 a ↔ 4π/a。這個對偶在「倒晶格的倒晶格是原晶格」下自洽:對 fcc(a) 取倒得 bcc(4π/a),再取倒得 fcc(4π/(4π/a)) = fcc(a) ✓
典型數值
Cu(a = 3.61 Å):倒晶格為 bcc,立方邊 4π/a = 3.481 Å⁻¹;|b| = √3(2π/a) = 3.015 Å⁻¹。Si(a = 5.430 Å):4π/a = 2.314 Å⁻¹,|b| = 2.004 Å⁻¹。

實算 fcc 的 ¼a³ 與三個叉積:由 (34) 走到 (36)補推導

設定:書本 p.38 同樣只給結果 (35)、(36) 與「體積是 4(2π/a)³」,中間的行列式與叉積一步也沒寫。這裡補齊,順便完成 Problem 3 在三種立方晶格上的第三次驗證。

解讀:把三次驗證放在一起看,V_BZ = (2π)³/V_c 這條式子已經不是「習題的答案」而是一條可以直接使用的工具:

晶格(立方邊 a)原胞體積 V_c倒晶格倒晶格立方邊V_BZ = (2π)³/V_c第一 BZ 形狀
sca³sc2π/a(2π/a)³立方體,邊長 2π/a
bcc½a³(式 30)fcc4π/a2(2π/a)³菱形十二面體(Fig. 13)
fcc¼a³(式 35)bcc4π/a4(2π/a)³截角八面體(Fig. 15)

注意最後一欄的體積比 1 : 2 : 4 恰好是正空間原胞體積比 1 : ½ : ¼ 的倒數——這是 V_BZ ∝ 1/V_c 的直接體現,也是檢查自己有沒有把 bcc 與 fcc 記反的最快方法:正空間原胞越小,倒空間布里淵區越大。fcc 的原胞最小(¼a³),所以 fcc 的布里淵區最大。

  1. 體積(式 35):V = (a/2)³ det[[0,1,1],[1,0,1],[1,1,0]]。
  2. 展開行列式:det = 0·(0·0 − 1·1) − 1·(1·0 − 1·1) + 1·(1·1 − 0·1) = 0 − 1·(−1) + 1·(1) = 2。
  3. 所以 V = (a³/8)·2 = ¼a³ ✓(物理檢查:慣用胞含 4 格點,a³/4 ✓)
  4. b₁ = 2π(a₂×a₃)/V。先算 a₂×a₃ = (a/2)²[(1,0,1)×(1,1,0)]。分量:x = (0)(0) − (1)(1) = −1;y = (1)(1) − (1)(0) = 1;z = (1)(1) − (0)(1) = 1。得 (−1,1,1)。
  5. 故 a₂×a₃ = (a²/4)(−x̂+ŷ+ẑ),b₁ = 2π(a²/4)(−x̂+ŷ+ẑ)/(a³/4) = (2π/a)(−x̂+ŷ+ẑ) ✓
  6. b₂:a₃×a₁ = (a²/4)[(1,1,0)×(0,1,1)]。分量:x = (1)(1) − (0)(1) = 1;y = (0)(0) − (1)(1) = −1;z = (1)(1) − (1)(0) = 1。得 (1,−1,1),故 b₂ = (2π/a)(x̂−ŷ+ẑ) ✓
  7. b₃:a₁×a₂ = (a²/4)[(0,1,1)×(1,0,1)]。分量:x = (1)(1) − (1)(0) = 1;y = (1)(1) − (0)(1) = 1;z = (0)(0) − (1)(1) = −1。得 (1,1,−1),故 b₃ = (2π/a)(x̂+ŷ−ẑ) ✓
  8. 驗證 (14):b₁·a₁ = (2π/a)(a/2)[(−1,1,1)·(0,1,1)] = π(0+1+1) = 2π ✓;b₁·a₂ = π[(−1,1,1)·(1,0,1)] = π(−1+0+1) = 0 ✓;b₁·a₃ = π[(−1,1,1)·(1,1,0)] = π(−1+1+0) = 0 ✓
  9. 倒晶格原胞體積:b₁·b₂×b₃ = (2π/a)³ det[[−1,1,1],[1,−1,1],[1,1,−1]] = (2π/a)³ × 4 = 4(2π/a)³,與書中 p.38 一致。
  10. Problem 3 檢查:(2π)³/V_c = (2π)³/(¼a³) = 4(2π/a)³ ✓ 第三次驗證通過。

fcc 倒晶格中最短的 8 條 G

(2π/a)(±x̂ ± ŷ ± ẑ)

(37)  p.38

三個正負號獨立選取,共 2³ = 8 條向量,方向為八條 ⟨111⟩ 體對角線。每條長度 √3(2π/a)。書中原句 “The shortest G's are the eight vectors”。

出處推導
由 (36) 張出的倒晶格是立方邊 4π/a 的 bcc。bcc 的最近鄰是體心到八個角,距離 = (√3/2)×立方邊 = (√3/2)(4π/a) = √3(2π/a) ✓,配位數 8 ✓。
極限行為
8 條最短 G 的平分面法線指向 ⟨111⟩,八個平面圍出的是一個正八面體(書中:“the corners of the octahedron thus formed”)。但這個八面體不是第一布里淵區,因為它的六個角伸得太遠、被次短的 G(式 38)截掉。這一點是 fcc 與 bcc 的分水嶺:bcc 的次短 G 恰好相切不切入,fcc 的次短 G 則真的切進去。判別方式是比較兩個殼層的「平分面距離」與「多面體頂點距離」,見下一個推導。
典型數值
Cu(a = 3.61 Å):|G| = √3(2π/a) = 3.015 Å⁻¹,平分面距原點 1.507 Å⁻¹(習稱 L 點)。Si(a = 5.430 Å):|G| = 2.004 Å⁻¹,L 點距 Γ 為 1.002 Å⁻¹。這組 G 對應的是 (111) 反射——diamond 結構的 (111) 是允許反射(指標全奇),也是 Si 晶圓最常用的解理面。

截掉八面體六角的次短 G

(2π/a)(±2x̂) ; (2π/a)(±2ŷ) ; (2π/a)(±2ẑ)

(38)  p.38

共 6 條,方向 ⟨100⟩,長度 2(2π/a) = 4π/a,即倒空間慣用立方胞的邊長。

出處推導
書中明確驗證了這些確實是倒晶格向量:“Note that (2π/a)(2x̂) is a reciprocal lattice vector because it is equal to b₂ + b₃.” 檢查:b₂ + b₃ = (2π/a)[(x̂−ŷ+ẑ) + (x̂+ŷ−ẑ)] = (2π/a)(2x̂) ✓ 這一步 Kittel 罕見地把驗算寫出來了,因為 (2π/a)(2x̂) 的形式看起來不像 (36) 的任何一支,容易讓人懷疑。
極限行為
這 6 個平分面位於 k_x = ±2π/a 等六處。八面體在 ⟨100⟩ 方向的頂點在 3π/a 處(見下一推導),2π/a < 3π/a,故確實切入 → 六個角被削平,形成 6 個正方形面。最終形狀是 8 個正六邊形 + 6 個正方形 = 14 面的截角八面體。極限判讀:純幾何上,若次短殼層的平分面退到 3π/a 之外就不會截角;但對真實晶格而言這不會發生,因為第一布里淵區是 WS 原胞、必須平移後無縫填滿空間,而正八面體做不到——這正是 fcc 一定會被 {200} 面截角的深層理由(可與後面「截角八面體是唯一能單獨填滿空間的 Archimedes 立體」一段對照)。bcc 那邊則是「差一點就截到」的臨界情形(相切)。
典型數值
Cu:截角平面在 k = 2π/a = 1.741 Å⁻¹(習稱 X 點),未截的八面體頂點原本在 3π/a = 2.611 Å⁻¹。Si:X 點在 1.157 Å⁻¹。

八面體被削角的完整幾何:為什麼會截、截掉多少、剩下多少補推導

設定:書本 p.38 用兩句話描述整個截角過程:“The boundaries of the central cell in the reciprocal lattice are determined for the most part by the eight planes normal to these vectors at their midpoints. But the corners of the octahedron thus formed are cut by the planes that are the perpendicular bisectors of six other reciprocal lattice vectors”。「for the most part」「are cut」這兩個措辭都是定性的——切了多少、為什麼一定會切、切完剩多少,全部沒有數字。這裡把它算完。

解讀:第四步的比較(2π/a vs 3π/a)是這整段最值得記住的判別法。一個殼層的 G 會不會對布里淵區的形狀有貢獻,取決於它的平分面(距離 ½|G|)有沒有落在前面幾個殼層已經圍出的多面體內部。用這個判別法你可以自己決定要考慮到第幾個殼層,而不必死記「fcc 要用兩層、bcc 只要一層」。

第五步的體積算法也值得注意:積分 ∫2(v−t)²dt 用的是「切片面積」的技巧,對任何凸多面體被平面截斷都適用,比硬算頂點座標再做四面體分解快得多。而第九步的 V = ⅓Σ(面積×距離) 則是針對「所有面都可從原點看見的凸多面體」的通用公式(把多面體分解成以原點為頂、各面為底的角錐)。兩種算法得到同一個 32π³/a³,加上 Problem 3 的公式,等於做了三重交叉驗證。

  1. 第一步:八個 {111} 平面圍出什麼。取 G = (2π/a)(x̂+ŷ+ẑ),½G = (π/a)(1,1,1),由 (26) 得平面 k·(π/a)(1,1,1) = (π/a)²·3,即 k_x + k_y + k_z = 3π/a。八個符號組合給出八個平面 ±k_x ±k_y ±k_z = 3π/a,圍出一個正八面體。
  2. 第二步:這個八面體的頂點在哪裡。令 k_y = k_z = 0,得 k_x = 3π/a。所以八面體的六個頂點在 (±3π/a,0,0)、(0,±3π/a,0)、(0,0,±3π/a),距原點 3π/a。
  3. (順帶:面到原點的距離是 ½|G| = (√3/2)(2π/a) = √3π/a ≈ 1.732π/a,比頂點距離 3π/a 小得多——八面體是很「尖」的形狀。)
  4. 第三步:次短 G 的平分面在哪裡。由 (38),G = (4π/a)x̂,½G = (2π/a)x̂,平面為 k_x = 2π/a。
  5. 第四步:判斷是否截入。2π/a < 3π/a,平面比頂點更靠近原點 → 確實切入,六個頂點被削掉。(對比 bcc:那裡的次短平分面在 2π/a,而多面體頂點也恰在 2π/a,相切、不截。)
  6. 第五步:算被削掉的體積。先算整顆八面體:頂點在 ±v(v = 3π/a)的正八面體體積 = (4/3)v³ = (4/3)(27π³/a³) = 36π³/a³。
  7. 每個被削掉的角是一個「小八面體帽」。在 x = t 處的截面是 |k_y| + |k_z| ≤ v − t 的正方形(旋轉 45°),面積 2(v−t)²。從 t = c(截面位置)到 t = v 積分:∫ 2(v−t)² dt = (2/3)(v−c)³。
  8. 代入 v = 3π/a、c = 2π/a:每個角錐體積 = (2/3)(π/a)³。六個角共削掉 6 × (2/3)(π/a)³ = 4π³/a³。
  9. 第六步:剩下的體積。36π³/a³ − 4π³/a³ = 32π³/a³。
  10. 第七步:與 Problem 3 對照。4(2π/a)³ = 4 × 8π³/a³ = 32π³/a³ ✓ 完全吻合;而 (2π)³/V_c = 8π³/(¼a³) = 32π³/a³ ✓ 三方一致。
  11. 第八步:確認形狀規格。截角八面體的頂點形如 (2π/a, π/a, 0) 及其排列與變號,共 24 個。檢查:它滿足 k_x = 2π/a(截面)與 k_x+k_y+k_z = 3π/a({111} 面)✓。相鄰頂點距離例如 (2π/a,π/a,0) 到 (2π/a,0,π/a) 為 √2π/a,到 (π/a,2π/a,0) 也是 √2π/a —— 所有棱等長,因此六邊形面是正六邊形、正方形面邊長也是 √2π/a。這是 Archimedes 截角八面體。
  12. 第九步:用面積法交叉驗算體積。V = ⅓Σ(面積×距離)。6 個正方形:邊長 √2π/a → 面積 2π²/a²,距離 2π/a;8 個正六邊形:面積 (3√3/2)(√2π/a)² = 3√3π²/a²,距離 √3π/a。V = ⅓[6·(2π²/a²)(2π/a) + 8·(3√3π²/a²)(√3π/a)] = ⅓[24π³/a³ + 72π³/a³] = 32π³/a³ ✓

p.38 最後一句的陷阱:「(4π/a)³」不是布里淵區的體積我的判讀

書本 p.38 的最後一句是:

The six planes bound a cube of edge 4π/a and (before truncation) of volume (4π/a)³.

這句話極容易被誤讀成「fcc 的第一布里淵區體積是 (4π/a)³」,理由有三:它出現在整節最後、它緊接在「the truncated octahedron shown in Fig. 15」之後、而且它是這一節唯一一個帶立方體積數字的句子。但它講的完全是另一回事。

把句子拆開看:

  • The six planes 指的是式 (38) 那六個平面 k_x = ±2π/a、k_y = ±2π/a、k_z = ±2π/a。
  • 這六個平面確實「圍出」一個立方體:邊長 = 2×(2π/a) = 4π/a ✓
  • (before truncation) 這個括號是關鍵限定詞。它的意思是「在截角這件事發生之前(也就是若只考慮這六個平面、不管那八個 {111} 平面)」,這六個平面圍出的立方體體積是 (4π/a)³。
  • 這個立方體其實就是倒空間的慣用立方胞(倒晶格是 bcc,立方邊 4π/a)。它含 2 個倒晶格點(bcc),因此體積是布里淵區的 2 倍。

做一個數字檢查,誤讀立刻現形:

量數值與 a 的關係
第一布里淵區(截角八面體)體積4(2π/a)³32π³/a³
六個平面圍出的立方體體積(4π/a)³64π³/a³
比值恰好 1/2

32π³/a³ 正好是 64π³/a³ 的一半。這個 ½ 不是巧合,它就是「倒空間 bcc 慣用胞含 2 個格點」的直接後果——每個格點分到胞體積的一半,而每個格點所佔的區域就是它的 WS 原胞,也就是布里淵區。

還有第二層誤讀要避免:那個立方體也不是「截角前的八面體」。截角前的形狀是八個 {111} 平面圍出的正八面體,體積 36π³/a³(見上一個推導),既不是 32 也不是 64。所以這一頁其實有三個不同的體積在打轉:

  • 32π³/a³ = 4(2π/a)³ —— 真正的第一布里淵區(截角八面體)
  • 36π³/a³ —— 未截角的正八面體(書中未給數字)
  • 64π³/a³ = (4π/a)³ —— 六個 {200} 平面圍出的立方體,即倒空間慣用胞

記住哪個是哪個,最可靠的辦法是永遠回到 Problem 3:V_BZ = (2π)³/V_c,V_c 用的是原胞體積(fcc 是 ¼a³),不是慣用立方胞體積。只要 V_c 沒有代錯,就不會拿到 64π³/a³。

Fig. 15 讀圖:截角八面體與「兩個胞」的畫法我的判讀

Fig. 15(p.38 上方)與 Fig. 13 的畫法不同,值得專門說明。圖說是:“Brillouin zones of the face-centered cubic lattice. The cells are in reciprocal space, and the reciprocal lattice is body centered.” 注意 zones 與 cells 都是複數。

圖中畫了一個立方框架(細線),右側用雙箭頭標註框架的邊長為 4π/a。框架的八個角與體心上有小球,構成 bcc 排列——這就是圖說最後一句在講的「倒晶格是體心的」。

框架裡畫了兩個大小相同的多面體,一個在左後方、一個在右前方,彼此部分重疊地擺著。每個多面體都有清楚可見的正六邊形面(圖中以淺灰與深灰交替上色)與較小的正方形面(在多面體的「頂部」與側面可以看到被截平的小面)。這就是截角八面體。

畫兩個而不是一個,用意是展示空間填充:這些胞是倒晶格的 WS 原胞,平移複本會無縫鋪滿整個倒空間、不留空隙也不重疊。兩個胞的位置分別對應框架的角點與體心——也就是 bcc 倒晶格上相鄰的兩個格點。這同時也視覺化地印證了上一段的體積關係:一個邊長 4π/a 的立方體裡裝得下 2 個這樣的胞,所以每個胞是 (4π/a)³/2 = 32π³/a³ ✓

規格整理(本區塊標 mine,數字來自前面的推導,書中只給了 4(2π/a)³ 與「truncated octahedron」):

項目規格
面數14 = 8 個正六邊形(法線 ⟨111⟩)+ 6 個正方形(法線 ⟨100⟩)
六邊形面到中心√3π/a,即式 37 那批 G 長度的一半
正方形面到中心2π/a,即式 38 那批 G 長度的一半
頂點24 個,形如 (2π/a, π/a, 0) 及其排列變號,距中心 (√5/2)(2π/a)
棱長全部相等,√2π/a
體積4(2π/a)³ = 32π³/a³

最後一個觀察:截角八面體是三維空間中唯一一個「單獨一種多面體就能填滿空間、且是某個 Bravais 晶格的 WS 原胞」的 Archimedes 立體。它在數學上又叫 permutohedron。這解釋了為什麼 fcc 材料(Cu、Al、Ag、Au、Si、Ge、GaAs……)的能帶圖看起來都很像——它們共用同一個布里淵區形狀,只差一個尺度 2π/a。

三組互倒關係的總整理與自洽檢查我的判讀

把 p.34–38 三個小節的結果收攏成一張表,並補上書中沒有並列的幾個交叉檢查(本區塊標 mine,表中「書中頁碼/式號」欄指出每一項的原文出處,其餘為推導所得)。

項目scbccfcc
正晶格原始平移向量(27a) a x̂ 等(29) ½a(−x̂+ŷ+ẑ) 等(34) ½a(ŷ+ẑ) 等
慣用立方胞含格點數124
原胞體積 V_ca³(p.35 正文)½a³(式 30)¼a³(式 35)
倒晶格scfccbcc
倒晶格慣用立方邊2π/a4π/a4π/a
倒晶格原始平移向量(27b)(31)(36)
倒晶格原胞體積(2π/a)³2(2π/a)³(p.37)4(2π/a)³(p.38)
最短 G:條數612(式 33)8(式 37)
最短 G:長度2π/a√2(2π/a)√3(2π/a)
次短 G 是否截角不適用否(相切)是(式 38,6 條)
第一 BZ 形狀立方體菱形十二面體(Fig. 13)截角八面體(Fig. 15)
第一 BZ 面數61214 (8+6)
V_BZ = (2π)³/V_c 驗證✓ 8π³/a³✓ 16π³/a³✓ 32π³/a³

三個自洽檢查,用來抓自己的錯:

檢查一:倒晶格的倒晶格要回到原晶格。fcc(a) 的倒晶格是 bcc(4π/a)。再取一次倒:bcc(A) 的倒晶格是 fcc(4π/A),代 A = 4π/a 得 fcc(4π/(4π/a)) = fcc(a) ✓ 若你算出來不是原晶格,一定是 4π/a 這個尺度記成了 2π/a。

檢查二:最短 G 的條數要等於倒晶格的配位數。bcc 的倒晶格是 fcc,fcc 配位數 12 → 12 條最短 G → 12 面 ✓;fcc 的倒晶格是 bcc,bcc 配位數 8 → 8 條最短 G → 八面體(再被截)✓;sc 配位數 6 → 6 面立方體 ✓。這是最快的一致性檢查。

檢查三:體積比要是原胞體積比的倒數。V_c 之比 sc : bcc : fcc = 1 : ½ : ¼,則 V_BZ 之比必為 1 : 2 : 4 ✓。

最容易記反的一件事,單獨提出來:bcc 的倒晶格是 fcc,fcc 的倒晶格是 bcc。記憶樁子是「BZ 的形狀」——bcc 的 BZ 是十二面體、fcc 的 BZ 是十四面體,而 fcc 的原胞比較小(¼a³)所以 BZ 比較大、面比較多。形狀較複雜(截角八面體)的那個,屬於原胞較小的 fcc。

這一段接到哪裡去:模數計數、能帶、費米面,以及熱容我的判讀

本區塊全部標 mine:書本 p.33–38 只給了幾何,以下是把這些幾何接到後面章節的說明,出處以書中目錄與各章標題頁為準。

(一)第 4 章 PHONONS I(p.89):k 的範圍。第 4 章 p.93 有一節就叫「First Brillouin Zone」。一維鏈的位移場 u_s = u exp(iska) 在 k 與 k + 2π/a 下給出完全相同的原子位移,所以 k 超出 (−π/a, π/a] 不帶新資訊——這正是 Fig. 11 那條線段的物理意義。三維情形下「不帶新資訊的重複」由整個倒晶格描述,可獨立的 k 就落在第一布里淵區內。

(二)第 5 章 PHONONS II(p.105):模數計數與 Debye 模型。這是 V_BZ = (2π)³/V_c 最重要的應用。晶體有 N 個原胞、每原胞 p 個原子 → 自由度 3pN → 必須有 3pN 個簡正模。而聲子分支有 3p 條、每條分支在第一布里淵區裡有 N 個允許的 k(見 Problem 3 那個推導區塊的算式),3p × N = 3pN ✓ 對得起來。Debye 模型做的近似,就是把形狀複雜的布里淵區換成一個等體積的球(Debye 球),好把積分做出來。等體積條件 (4/3)πk_D³ = (2π)³/V_c 給出 k_D = (6π²/V_c)^(1/3)。

(三)第 6/7 章:一條能帶 2N 個態。每個 k 一個軌域、自旋二重 → 2N。第 7 章 p.180「Number of Orbitals in a Band」與 p.181「Metals and Insulators」的結論(每原胞奇數個價電子必為金屬)完全建立在這個計數上。

(四)第 9 章 FERMI SURFACES AND METALS(p.221):費米球與布里淵區的相對大小。自由電子費米球體積 (4/3)πk_F³ 除以 V_BZ,恰好等於「每原胞價電子數 ÷ 2」。用 Na 檢查:bcc,a = 4.225 Å,每原子 1 個價電子、每原胞 1 個原子。n = 1/V_c = 1/37.71 ų = 0.02652 Å⁻³,k_F = (3π²n)^(1/3) = 0.923 Å⁻¹。最近的布里淵區邊界是 N 點面,距離 √2π/a = 1.052 Å⁻¹。k_F/k_N = 0.877 < 1,費米球完全裝在布里淵區裡而且不碰到邊界——這就是 Na 被當成近自由電子金屬教科書範例的定量理由。體積比檢查:(4/3)π(0.923)³/6.578 = 3.289/6.578 = 0.500 ✓ 恰好半滿。對比 Cu(fcc,a = 3.61 Å,1 價電子):k_F = 1.360 Å⁻¹,而 L 點在 1.507 Å⁻¹、X 點在 1.741 Å⁻¹;自由電子球仍未觸及邊界,但實際 Cu 的能帶在 L 點附近被 d 帶推高,費米面在 ⟨111⟩ 方向鼓出「頸部」接觸 L 面——書中第 9 章 p.249 有一整節「Fermi Surface of Copper」在講這件事。

(五)與 FinFET 熱容(Cth)研究的接點。若你的目標是從第一原理/聲子層次理解矽的熱容與熱時間常數,本段的幾何就是計數的起點,而且有兩個容易踩的坑:

  • Si 是 diamond 結構,晶格是 fcc、每原胞 2 個原子。原胞體積 V_c = ¼a³ = 40.03 ų(a = 5.430 Å),不是 a³ = 160.1 ų。用錯會讓模數差 4 倍。
  • 布里淵區體積 V_BZ = (2π)³/V_c = 6.197 Å⁻³。若只把「每原胞一個 k」的那 N 個點換成等體積球,半徑 k_D = (6π²/V_c)^(1/3) = 1.139 Å⁻¹,落在 L 點(1.002 Å⁻¹)與 X 點(1.157 Å⁻¹)之間——這說明 Debye 球是布里淵區的一個相當合理的「平均半徑」近似。
  • 但標準 Debye 模型把 3pN = 6N 個模全塞進一個球時,k_D = (6π²·n_atom)^(1/3)(n_atom = 8/a³ = 0.04997 Å⁻³)= 1.436 Å⁻¹,已經超出第一布里淵區的最遠邊界。這不是計算錯誤,而是 Debye 模型在多原子基底下必然付出的代價:它用一個超出 BZ 的球來裝原本應該分佈在光學支上的模。做 Cth 的定量工作時,這是「Debye 溫度擬合值會隨溫度區間漂移」的根本原因之一,而不是量測誤差。

一句話總結本段在全書中的位置:第 2 章把倒空間切成布里淵區,之後每一章需要對 k 積分或對模計數時,積分區域都是它。

常見誤解與本段落校註我的判讀

誤解一:「布里淵區是原胞,所以它是平行六面體。」錯。布里淵區是 WS 原胞,不是平行六面體原胞。兩者體積相同(都是 (2π)³/V_c),但形狀不同。WS 原胞的價值在於它保有晶格的完整點群對稱性,能帶圖才畫得出對稱點。

誤解二:「bcc 的倒晶格是 fcc,所以尺度換算是 a → 2π/a。」錯,是 4π/a。因子 2 的來源是:bcc/fcc 的慣用立方胞不是原胞。只有 sc 才是 a → 2π/a。

誤解三:「(26) 是 (23) 的近似。」錯,兩者嚴格等價,只差兩邊除以 4。書中原句是 “We divide both sides by 4 to obtain (26)”,沒有任何近似。

誤解四:「Δk = G 和 Δk = −G 有一個是錯的。」都對。倒晶格中心對稱,G 與 −G 同時存在。(21) 與 p.33 的差異來自 p.32 那次 G → −G 的代換,見前面的符號簿記推導。

誤解五:「(4π/a)³ 是 fcc 的布里淵區體積。」錯,那是六個 {200} 平面圍出的立方體(= 倒空間慣用胞)的體積,是 BZ 的 2 倍。詳見前面專門討論這一句的區塊。

誤解六:「fcc 的 BZ 是八面體。」只有在忽略式 (38) 時才是。書中明言八面體的角 “are cut”,最終形狀是 14 面的截角八面體。

誤解七:「布里淵區邊界上的 k 都會被繞射,所以晶體對這些波不透明。」要小心。(26) 說的是「滿足繞射的運動學條件」;強度還要乘上結構因子 S_G(p.39 之後)。若 S_G = 0(例如 bcc 的 v₁+v₂+v₃ 為奇),該倒晶格點的反射強度為零。不過要注意:對 Bravais 晶格本身而言,消光的那些點根本不在倒晶格裡(見前面的奇偶論證),所以這個顧慮只在「用過大的慣用胞去標記」時才出現。


校註一(書本 p.36)。原文 “The primitive translations of the reciprocal lattice are defined by (13). We have, using (28),” 中的式號 (28) 應為 (30)。理由:(28) 是 sc 的布里淵區邊界 ±½b₁ = ±(π/a)x̂,與 bcc 的 b 向量計算無關;計算 (31) 實際需要的是 (13)(定義式)與 (30)(V = ½a³,當分母)。依據:已用 PyMuPDF 將該頁算繪成影像目視核對,印刷版確實印「(28)」,故判定為原書排印誤植,非文字層抽取誤差。

校註二(書本 p.36)。原文寫 “The boundaries of the first Brillouin zones are the planes...”,zones 為複數;依上下文(講的是 sc 的單一第一布里淵區)應為單數 zone。不影響理解。

校註三(書本 p.37)。「八個角各佔 ⅛」那句後面標註 (see Fig. 12),但 Fig. 12 畫的是 bcc 正晶格的原始基矢,不是倒空間的平行六面體。圖中確有多個相鄰立方胞共用角點的畫面可作示意,但交叉引用容易讓讀者找錯圖。

校註四(Γ/X/L/W/N/H/P 等標記)。本講解在數個地方使用了 Γ、X、L、N、H、P 這些布里淵區對稱點的字母標記,方便與後面章節的能帶圖對照。這些標記不出現在書本 p.33–38,屬標準慣例的補充。

校註五(晶格常數數值)。本講解用到的晶格常數(Cu 3.61 Å、Si 5.430 Å、Na 4.225 Å(5 K)、W 3.16 Å)取自書本 p.20 Table 3。該頁在 PDF 中沒有文字層(純影像),數值係以 PyMuPDF 算繪放大後目視判讀,並非文字抽取。