Kittel·固態物理

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

05基底的傅立葉分析——結構因子、原子形狀因子;本章摘要與習題

FOURIER ANALYSIS OF THE BASIS — Structure Factors, Atomic Form Factor; SUMMARY; Problems

書本 p.39-45 · 小節 5 / 5 · 概念 15 推導 10 公式 9

前面幾節只用了晶格、只回答了「繞射峰出現在哪裡」;這一段把基底放回去,回答「峰有多亮、哪些峰根本不出現」——結構因子管原子之間的干涉,原子形狀因子管原子內部的干涉。

這一段在補哪一塊拼圖我的判讀

第 2 章走到這裡,前面幾節已經把一句話講完了:晶體的週期性把散射振幅壓縮到倒晶格點上。從 (18) 的 F = ∫ dV n(r) exp(−iΔk·r) 出發,代入電子密度的傅立葉展開 (9),得到 (20) F = Σ_G ∫ dV n_G exp[i(G − Δk)·r];只有當 Δk 正好等於某個倒晶格向量時指數才不振盪,於是有 (21) Δk = G。這條鏈的終點告訴你繞射峰出現在哪些方向——但它一個字也沒說峰有多亮。

問題出在哪?出在前面那幾節從頭到尾只用了「晶格」這個資訊。晶格是一組平移向量,它決定了倒晶格、決定了哪些 G 存在、決定了 Brillouin zone 的形狀。可是第 1 章早就說過,晶體結構 = 晶格 + 基底(basis)。基底的資訊還完全沒被用上。兩個晶體可以有一模一樣的晶格(因此有一模一樣的倒晶格、一模一樣的允許 G 集合),卻因為基底不同而有完全不同的繞射圖樣——甚至有些本來「允許」的 G 會整條消失。

這一段(書本 p.39–42)做的就是把基底放回去。做法是把 (21) 成立時的振幅再算一次、算精確:不只問「非零嗎」,而是問「等於多少」。答案分兩層:外層是結構因子 S_G,管一個晶胞裡 s 個原子之間的相位干涉;內層是原子形狀因子 f_j,管單一個原子內部電子雲的干涉。這兩層拆分之所以漂亮,是因為 f_j 只是原子的性質(跟晶體無關),而 S_G 只是幾何的性質(原子擺在哪裡)——把兩者乘起來相加,就得到峰的強度。

為什麼這件事值得花四頁?因為消光規則(extinction rule)就是從這裡長出來的。「bcc 看不到 (100)」「fcc 只看得到全奇或全偶」「鑽石結構連 (222) 都沒有」,這些是實驗上判別晶體結構的第一件武器:你不需要知道晶格常數,光看哪些峰缺席就能把 sc / bcc / fcc / diamond 區分開。這一段是整章唯一真正接上「拿到一張繞射圖後你要怎麼辦」的地方。

N 個晶胞的散射振幅與結構因子

F_G = N ∫_cell dV n(r) exp(−iG·r) = N S_G

(39)  p.39

F_G:滿足繞射條件 Δk = G 時的散射振幅;N:晶體中的晶胞數;n(r):電子濃度(單位體積的電子數);G:倒晶格向量;S_G:結構因子,定義成對單一晶胞的積分,積分原點 r = 0 取在晶胞的一個角落。

出處推導
由 (18) 的散射振幅 F = ∫ dV n(r) exp(−iΔk·r) 出發,代入繞射條件 (21) Δk = G。書上 p.39 開頭逐字寫的是「When the diffraction condition Δk = G of Eq. (21) is satisfied, the scattering amplitude (18) for a crystal of N cells may be written as」。關鍵步驟是:因為 exp(−iG·r) 在平移 T 之下不變(exp(−iG·T) = 1,這正是 (16)(17) 證過的),整個晶體的積分就等於一個晶胞的積分乘上晶胞數 N。
極限行為
N 的角色純粹是「有多少個晶胞在同調地一起散射」。強度 ∝ |F_G|² ∝ N²,這個 N² 不是筆誤:在峰的正中心,所有晶胞的貢獻同相位相加,振幅 ∝ N。但峰在倒空間每個方向上的寬度 ∝ 1/M,其中 M = N^(1/3) 是該方向上的晶胞數(Problem 4 證的就是這個);峰在三維倒空間所佔的體積因此 ∝ 1/M³ = 1/N,所以積分強度(峰下的體積)只 ∝ N² × (1/N) = N,也就是 ∝ 樣品體積——這才是能量守恆容許的結果。若 Δk 偏離 G,(20) 的指數開始振盪,F 迅速趨近零;書上 p.31 說這是「a simple exercise (Problem 4)」。
典型數值
(代入計算,書中未列)一顆 1 μm 的立方微晶、晶格常數 a ≈ 3 Å,每邊約 3.3×10³ 個晶胞,N ≈ 3.7×10¹⁰。峰中心強度 ∝ N² ≈ 1.4×10²¹,峰寬 ∝ 1/3.3×10³ ≈ 3×10⁻⁴ rad ≈ 0.017°。這就是繞射峰「又高又窄」的來源。

S_G 的定義為什麼要「積分原點取在角落」我的判讀

(39) 之後書上補了一句很容易被跳過的話:S_G「is defined as an integral over a single cell, with r = 0 at one corner」。書上就只寫了這一句,以下的理由是本站補的。這句話不是廢話,因為 S_G 的相位完全取決於你把原點放哪裡。

理由很直接:如果把所有原子位置平移一個常向量 d(等於把原點移動 −d),(43) 的每一項都多出同一個因子 exp(−iG·d),於是 S_G → S_G · exp(−iG·d)。S 的模長不變,相位整體旋轉。這帶來兩個後果:

第一,|S|² 是物理的、S 的相位不是。後面 p.40 會明講可觀測量是 S*S。換句話說,同一個晶體,你選不同的晶胞原點會算出不同的 S,但算出同樣的強度。這是自洽的,因為強度才是量得到的東西。

第二,做消光規則時務必用同一套原點。bcc 用 (000, ½½½)、fcc 用 (000, 0½½, ½0½, ½½0),這些都是書上採用的「角落為原點」慣例。如果你把 bcc 原點挪到體心,基底變成 (½½½, 000)——看起來一樣,但如果挪到別的位置(例如 ¼¼¼ 與 ¾¾¾),你會算出一個帶整體相位的 S,模長仍然對,可是逐項對照書上的式子就會對不起來。學生在這裡最常見的困惑,其實是「我算出來差一個相位因子,是不是錯了」——不是錯,是原點不同。

把晶胞電子密度拆成各原子的貢獻

n(r) = Σ_{j=1}^{s} nⱼ(r − rⱼ)

(40)  p.39

s:基底中的原子數;rⱼ:第 j 個原子中心的位置向量;nⱼ(r − rⱼ):第 j 個原子對位置 r 處電子濃度的貢獻,書上逐字說它「defines the contribution of that atom to the electron concentration at r」。

出處推導
這不是推導出來的,是一個建模選擇。書上 p.39 的說法是「Often it is useful to write the electron concentration n(r) as the superposition of electron concentration functions nⱼ associated with each atom j of the cell」——用的是 useful(有用),不是 exact(正確)。下一句立刻承認這個拆法不唯一。校註:書上這一式的求和下限印成 Σ_{j−1}^{s}(j 減 1),正確應為 j = 1。已用 PyMuPDF 14 倍放大目視確認,不是文字層抽取造成的。
極限行為
極限情形一:若各原子的電子雲完全不重疊(重原子、內層電子主導),這個拆法幾乎是精確的,nⱼ 就是自由原子的電子密度。極限情形二:共價鍵晶體(Si、鑽石)的鍵結電荷堆積在兩原子中間,「這團電荷屬於誰」沒有唯一答案——但這不影響 n(r) 的總和,所以也不影響 (39) 算出的 S_G。不唯一性被吸收進 fⱼ 的定義裡,互相抵消。

「分解不唯一」——Kittel 這句誠實話的分量

(40) 下面書上寫了兩句話,值得逐字讀:「The decomposition of n(r) is not unique, for we cannot always say how much charge density is associated with each atom. This is not an important difficulty.」(p.39)

這兩句其實包含一個完整的論證,只是壓縮得太狠,可以拆成三個層次來讀:

第一層(第一句的前半):n(r) 是一個連續分布,把它切成幾塊「屬於某個原子」的部分,切法有無限多種。你可以用 Voronoi 多面體硬切、可以用自由原子密度的疊加去擬合、可以用量子化學的 Bader 分析——每種切法給出不同的 nⱼ。

第二層(第一句的後半)是理由:電荷密度沒有「名牌」。這在共價晶體特別尖銳:Si—Si 鍵中間那團電荷,你說屬於左邊那顆還是右邊那顆?兩邊各半是一個約定,不是一個事實。

第三層(第二句)是關鍵,而書上沒解釋為什麼。原因是:可觀測量 S_G 完全不依賴這個選擇。(39) 定義的 S_G 只用到 n(r) 這個總和;(40)(41)(42)(43) 只是把同一個積分重新分組成 Σⱼ fⱼ exp(−iG·rⱼ)。你換一種分解,fⱼ 會變、rⱼ 甚至可以跟著變,但 Σⱼ fⱼ exp(−iG·rⱼ) 這個和不變。不唯一性是「表述」的不唯一,不是「結果」的不唯一。

這一點在實務上有具體後果:X 光結構精算(refinement)時,如果你用自由原子的 fⱼ 去擬合,擬合殘差裡就會留下鍵結電荷的痕跡(所謂 deformation density)。這不是模型錯,是模型的分解方式不含那塊電荷。p.42 末段對這件事有一段拿捏得很好的說明,我們稍後會逐字讀。

(39) → (41) → (43):把積分變數換掉的三行代數補推導

書上從 (39) 到 (43) 只印了兩個等號,中間跳過了「變數變換」與「為什麼可以把積分域從晶胞放大到全空間」這兩個步驟。這裡逐步補上。這段代數本身不難,但它是整條結構因子鏈的樞紐:把「晶胞的積分」重寫成「原子的積分 × 相位因子」之後,晶體幾何與原子性質就徹底分離了。

  1. 起點是 (39) 的定義:S_G = ∫_cell dV n(r) exp(−iG·r)。積分只跑一個晶胞。
  2. 代入 (40) 的分解 n(r) = Σⱼ nⱼ(r − rⱼ),並把求和拉到積分外(有限項求和與積分可交換):S_G = Σⱼ ∫ dV nⱼ(r − rⱼ) exp(−iG·r)。這就是書上 (41) 的第一行——注意書上這一行印的就是不帶 cell 下標的 ∫ dV(整段只有 (39) 印了 ∫_cell),這個下標的消失,正是下一步要講的「積分域悄悄放大」。
  3. 這裡有一個書上沒說的小手腳:積分域從 cell 悄悄放大成全空間。理由是 nⱼ 是局域的——第 j 顆原子的電子雲集中在 rⱼ 附近、在晶胞外迅速衰減;而落在鄰胞裡的那部分尾巴,等價於鄰胞裡對應原子的尾巴落進本胞。因為 exp(−iG·r) 對平移 T 不變(exp(−iG·T) = 1),這兩塊的相位相同,加起來剛好補齊成「一顆完整原子在全空間的積分」。這就是為什麼 (42) 下面書上寫「integrated over all space」。
  4. 作變數變換 ρ ≡ r − rⱼ(書上在 (41) 之後逐字定義 ρ = r − rⱼ)。這是平移,Jacobian = 1,所以 dV 不變,只是積分變數換名字。
  5. 指數拆開:exp(−iG·r) = exp[−iG·(ρ + rⱼ)] = exp(−iG·rⱼ) · exp(−iG·ρ)。關鍵在於 exp(−iG·rⱼ) 與積分變數 ρ 無關,可以提到積分外面。
  6. 於是 S_G = Σⱼ exp(−iG·rⱼ) ∫ dV nⱼ(ρ) exp(−iG·ρ)。這正是書上 (41) 的第二行。
  7. 把括號裡那個積分命名為 fⱼ,就是 (42):fⱼ = ∫ dV nⱼ(ρ) exp(−iG·ρ)。書上緊接著說「If nⱼ(ρ) is an atomic property, fⱼ is an atomic property」——意思是這個積分裡完全沒有晶體的資訊,只有原子自己的電子雲與 G 的大小。
  8. 合併得 (43):S_G = Σⱼ fⱼ exp(−iG·rⱼ)。乘積形式:原子性質 × 幾何相位,對基底求和。

原子形狀因子(atomic form factor)的定義

fⱼ = ∫ dV nⱼ(ρ) exp(−iG·ρ), ρ = r − rⱼ

(42)  p.39

fⱼ:第 j 個原子的原子形狀因子(也常譯作原子散射因子);nⱼ(ρ):以該原子中心為原點量的電子濃度;ρ:從原子中心量起的位置向量;積分跑遍全空間。

出處推導
由 (41) 第二行的括號內積分直接命名而來。數學上它就是原子電子密度 nⱼ 的三維傅立葉變換,取在波向量 G 上。
極限行為
G → 0(前向散射)時 exp(−iG·ρ) → 1,積分退化成 ∫ dV nⱼ = Z,即原子的電子總數;這正是稍後 (51) 的內容。另一個極限:G → ∞ 時被積函數劇烈振盪,fⱼ → 0——物理意義是高角度散射看到的是電子雲的精細結構,而電子雲在小尺度上是平滑的,所以沒有高頻成分。這就是為什麼繞射圖樣的峰強度隨 2θ 增大而整體下滑。順帶一提,這也解釋了為什麼中子繞射(散射自點狀的原子核)沒有這個衰減——Fig. 18 的鑽石中子圖樣裡 (400) 還看得見,就是這個道理的體現之一。
典型數值
(代入計算,書中未列)fⱼ 是無因次量,單位就是「等效電子數」。典型尺度:對 Cu Kα(λ = 1.5418 Å),最大可及的 G = 4π sinθ/λ ≈ 8.1 Å⁻¹;中等重量的原子在這個 G 之下 f 大約掉到 Z 的 20–30%。

基底的結構因子(向量形式)

S_G = Σⱼ fⱼ exp(−iG·rⱼ)

(43)  p.39

求和跑遍基底中的 s 個原子;rⱼ 是第 j 個原子在晶胞中的位置。

出處推導
由 (41) 與 (42) 合併。書上逐字說「We combine (41) and (42) to obtain the structure factor of the basis in the form」。
極限行為
s = 1(單原子基底,如 sc 晶格)時 r₁ = 0,S_G = f,對所有 G 都不為零——簡單晶格沒有消光,這是 Fig. 17 裡 KCl 「看起來像 sc」時只出現偶指標的邏輯前提。所有 fⱼ 相等(同種原子)時 S_G = f · Σⱼ exp(−iG·rⱼ),消光條件退化成純幾何條件,與原子是什麼完全無關——bcc 與 fcc 的規則就屬於這一類。反之若基底原子不同(KBr、閃鋅礦),幾何上「該抵消」的項抵消不乾淨,會留下 |f_A − f_B| 這種弱峰。這正是 Problem 7 的一維版本。 熱擾動的極限:本書附錄 A(p.641 起,Temperature Dependence of the Reflection Lines)就是拿這一式做熱平均——讓 rⱼ 帶上隨時間漲落的位移 u(t),取 ⟨exp(−iG·u)⟩ 後得到 Debye–Waller 因子 I = I₀ exp(−⅓⟨u²⟩G²)。也就是說,(43) 不只是第 2 章的結論,它是全書處理繞射強度的共同起點。

(44)(45)(46):倒晶格的正交關係在這裡收割補推導

從 (43) 到 (46) 是這一段最該講透、也最容易被當成「換個寫法」而滑過去的一步。它的實質是:把 (43) 裡那個難算的內積 G·rⱼ,用第 2 段辛苦建立的 bᵢ·aⱼ = 2πδᵢⱼ 化成純數字。如果你覺得第 2 段引入倒晶格基底 bᵢ 時那個 2π 的因子很突兀,這裡就是它的回報。

  1. 把原子位置用晶格基底展開((44),書上註明「as in (1.2)」,即第 1 章 (2) 式的形式):rⱼ = xⱼa₁ + yⱼa₂ + zⱼa₃。這裡 xⱼ, yⱼ, zⱼ 是分數座標,落在 0 到 1 之間,例如 bcc 的體心是 (½, ½, ½)。
  2. 把倒晶格向量用 (15) 展開:G = v₁b₁ + v₂b₂ + v₃b₃,v₁v₂v₃ 是整數(即反射的指標)。
  3. 作內積並全部展開,共九項:G·rⱼ = (v₁b₁ + v₂b₂ + v₃b₃)·(xⱼa₁ + yⱼa₂ + zⱼa₃) = Σ_{i,k} vᵢ xⱼ^{(k)} (bᵢ·a_k)。
  4. 現在用第 2 段的 (14):bᵢ·aⱼ = 2πδᵢⱼ(書本 p.29)。δᵢⱼ 把九項殺掉六項,只留下 i = k 的三項對角項,每一項帶一個 2π:
  5. G·rⱼ = 2π(v₁xⱼ + v₂yⱼ + v₃zⱼ)。這就是書上的 (45)。這一步是整條鏈唯一用到倒晶格定義的地方——沒有 (13)(14),(46) 寫不出來。
  6. 代回 (43),得到 (46):S_G(v₁v₂v₃) = Σⱼ fⱼ exp[−i2π(v₁xⱼ + v₂yⱼ + v₃zⱼ)]。
  7. 結果的意義:S 的相位只剩下「整數指標 × 分數座標」的組合,完全不含 a₁a₂a₃ 的長度與夾角。也就是說,消光規則只取決於原子的分數座標與反射指標,與晶格常數無關——這就是為什麼「bcc 看不到奇數和的反射」對所有 bcc 金屬都成立,不管是 Na、Fe 還是 W。

結構因子的指標形式(最常用的形式)

S_G(v₁v₂v₃) = Σⱼ fⱼ exp[−i2π(v₁xⱼ + v₂yⱼ + v₃zⱼ)]

(46)  p.40

v₁v₂v₃:反射的指標(整數,可正可負可為零);xⱼ, yⱼ, zⱼ:第 j 個原子在晶胞中的分數座標(以 a₁, a₂, a₃ 為單位);fⱼ:該原子的形狀因子。

出處推導
由 (43) 代入 (45) 的 G·rⱼ = 2π(v₁xⱼ + v₂yⱼ + v₃zⱼ),而 (45) 本身靠的是第 2 段 (14) 的正交關係 bᵢ·aⱼ = 2πδᵢⱼ。
極限行為
全部 vᵢ = 0(G = 0,前向方向)時所有指數為 1,S = Σⱼ fⱼ = 晶胞中的總電子數。這是 S 的上界——任何反射的 |S| 都不可能超過它,因為那是「全部電子完全同相」的情形。另一方面,若基底原子的分數座標都是有理數(實際晶體都是),指數只取有限多個離散相位值,這就是為什麼消光規則是「週期性」的整數規則(bcc 看 v₁+v₂+v₃ 的奇偶、鑽石看 mod 4)而不是連續的衰減。
典型數值
(代入計算,書中未列)以 bcc Fe 的 (110) 為例:兩顆 Fe,f(Fe) 在該 G 約 20 左右,x,y,z 分別是 (0,0,0) 與 (½,½,½),指數 = exp[−i2π(½+½+0)] = exp(−i2π) = 1,S = 2f ≈ 40。而 (100):exp[−i2π(½)] = −1,S = f(1 − 1) = 0。

S 一般是複數:相位在哪裡遺失了

(46) 之後書上緊接著一句話:「The structure factor S need not be real because the scattered intensity will involve S*S, where S* is the complex conjugate of S so that S*S is real.」(p.40)

這句話短,但它藏著晶體學一百年來的核心困難,值得展開。

為什麼 S 會是複數。 看 (46):每一項是 fⱼ(實數,至少在忽略異常散射時是)乘上一個單位複數 exp[−i2πφⱼ]。這些單位複數在複平面上指向不同方向,加起來的結果一般既有實部也有虛部。只有在特殊情況下 S 才保證為實:若晶體結構有反演對稱中心,而且你把原點選在反演中心上,則每個位於 rⱼ 的原子都有一個位於 −rⱼ 的同伴,exp(−iG·rⱼ) 與 exp(+iG·rⱼ) 成對出現,加起來是 2cos(G·rⱼ),是實數。bcc 與 fcc 都屬於這一類,這就是為什麼 (47)(48) 算出來都是實數。

為什麼這件事很嚴重。 實驗量到的是強度 I ∝ S*S = |S|²,也就是說你量到 |S|,但 S 的輻角完全丟失。而 (46) 告訴你,原子位置 (xⱼ, yⱼ, zⱼ) 的資訊主要就編碼在輻角裡。把電子密度從繞射資料反推回來,需要對 n_G 做逆傅立葉(見 (9):n(r) = Σ_G n_G exp(iG·r)),而 n_G 正比於 S_G——缺了相位,逆變換做不了。這就是晶體學裡的相位問題;本書在這一段沒有用這個名詞,但 Kittel 寫下「S need not be real」這句話時,講的就是它的來源。

實務上怎麼繞過去。 三條路:靠重原子(重原子的 f 大,它的相位主導整體)、靠同構取代、靠直接法(利用電子密度必須非負這個約束做統計推斷)。這些都超出本章,但研究生應該知道「量到 |S| 而不是 S」這件事有代價,而不是把 S*S 當成一個純粹的數學步驟。

bcc 晶格的結構因子

S(v₁v₂v₃) = f {1 + exp[−iπ(v₁ + v₂ + v₃)]}

(47)  p.40

f:該原子的形狀因子(bcc 兩顆原子相同,所以只有一個 f);v₁v₂v₃:以立方晶胞(conventional cubic cell)為參考的反射指標。

出處推導
把 bcc 的基底代入 (46)。書上 p.40 逐字給的基底是「identical atoms at x₁ = y₁ = z₁ = 0 and at x₂ = y₂ = z₂ = ½」。第一顆給出 exp(0) = 1;第二顆給出 exp[−i2π(½v₁ + ½v₂ + ½v₃)] = exp[−iπ(v₁+v₂+v₃)]。
極限行為
v₁+v₂+v₃ 為奇數 → 指數 = exp(−iπ×奇數) = −1 → S = 0(消光);v₁+v₂+v₃ 為偶數 → 指數 = +1 → S = 2f。書上就是這樣逐字列出兩行的。注意 S = 2f 而不是 f,因為立方晶胞裡有兩顆原子——強度 ∝ 4f²,是相同 f 的 sc 晶格的四倍。重要前提:這是以立方晶胞為參考的結果。若改用 bcc 的原始晶胞(primitive cell,只含一顆原子),基底只有一顆,S = f 對所有反射都不為零——沒有任何消光。消光規則是「用了非原始晶胞」的產物,它並沒有消滅任何物理上的繞射峰,只是把「立方晶胞的倒晶格」裡那些不屬於真正倒晶格的點標記出來。
典型數值
(代入計算,書中未列)金屬鈉 bcc,a ≈ 4.23 Å。允許的最低角反射是 (110):d = a/√2 = 2.99 Å;用 Cu Kα(λ = 1.5418 Å)由 (24) 2d sinθ = λ 得 sinθ = 0.258,2θ ≈ 29.9°。而 (100) 若存在會落在 d = 4.23 Å、2θ ≈ 21.0°——實驗上那裡什麼都沒有。

bcc 消光規則的逐步驗算,與鈉的譜線清單補推導

書上 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; here the indices (v₁v₂v₃) are referred to a cubic cell.」

  1. 先把 (47) 的判準寫乾淨:令 σ ≡ v₁ + v₂ + v₃。exp(−iπσ) = (−1)^σ,所以 S = f[1 + (−1)^σ]。σ 奇 → S = 0;σ 偶 → S = 2f。整個 bcc 消光規則就是一個奇偶判斷,連 f 都不用知道。
  2. (100):σ = 1 + 0 + 0 = 1,奇 → S = 0。缺席,與書上一致。
  3. (300):σ = 3,奇 → S = 0。缺席。注意 (300) 與 (100) 同方向、d 值是 (100) 的三分之一,但因為 σ 同為奇數,一起消失。
  4. (111):σ = 3,奇 → S = 0。缺席。這一條特別值得記住,因為 (111) 是 fcc 最強的低角反射之一——「有沒有 (111)」是 bcc 與 fcc 的第一個分水嶺。
  5. (221):σ = 2 + 2 + 1 = 5,奇 → S = 0。缺席。
  6. (110):σ = 2,偶 → S = 2f。出現,而且這是 bcc 的最低角反射(d = a/√2 最大)。
  7. (200):σ = 2,偶 → S = 2f。出現。
  8. (222):σ = 6,偶 → S = 2f。出現。
  9. 把允許的反射按 v₁²+v₂²+v₃² 排序,得到 bcc 的特徵序列 2, 4, 6, 8, 10, 12, …:(110), (200), (211), (220), (310), (222), …。相鄰峰的 sin²θ 比值為 1 : 2 : 3 : 4 : 5 : 6,這是粉末繞射圖上辨認 bcc 的標準指紋。
  10. 對照組:sc 的序列是 1, 2, 3, 4, 5, 6, 8,…(缺 7,因為 7 無法寫成三個平方和);fcc 的序列是 3, 4, 8, 11, 12, 16, …。三者一眼可分。

Fig. 16:消光在實空間長什麼樣(本頁無圖,以下為文字還原)

算出 S = 0 是一回事,理解為什麼會是零是另一回事。書上 p.40 用一整段回到實空間解釋,並配上 Fig. 16。網頁不放圖,這裡把圖的內容描述清楚。

Fig. 16 畫的是什麼(書本 p.41,目視判讀):畫面是 bcc 晶格的側視剖面,畫成一列三個並排的方胞(每格四個角上有黑點=角原子,每格正中央有一個黑點=體心原子)。橫貫這一列的三條水平線由上而下標為「1st plane」「2nd plane」「3rd plane」:第 1 與第 3 條穿過角原子,第 2 條穿過體心原子,正好在前兩者的中間。右側有一個雙箭頭標著 a,表示第 1 與第 3 平面的間距就是立方晶格常數 a。左上方有三條平行的白色入射/反射光線,分別打在三個平面上;線之間用大括號標出相位差:第 1 與第 2 平面之間是 π,第 2 與第 3 之間也是 π,而第 1 與第 3 之間標的是「Phase difference 2π」。

圖說(書本 p.41 逐字):「Explanation of the absence of a (100) reflection from a body-centered cubic lattice. The phase difference between successive planes is π, so that the reflected amplitude from two adjacent planes is 1 + e^{−iπ} = 1 − 1 = 0.」

這張圖在說的邏輯鏈,書上 p.40 的正文寫得更完整:

(1)「The (100) reflection normally occurs when reflections from the planes that bound the cubic cell differ in phase by 2π.」——(100) 反射的定義就是「相鄰的立方晶胞邊界面相位差 2π」,這是 Bragg 條件在 d = a 時的內容。

(2)「In the bcc lattice there is an intervening plane (Fig. 16) of atoms, labeled the second plane in the figure, which is equal in scattering power to the other planes.」——中間插了一層,而且散射能力完全相同。

(3)「Situated midway between them, it gives a reflection retarded in phase by π with respect to the first plane, thereby canceling the contribution from that plane.」——位置在正中間 ⇒ 相位落後恰好 π ⇒ 逐一抵消。

關鍵在第 (2) 句,而不是第 (3) 句。 相位差 π 只保證兩者反相;要完全抵消還需要振幅相等。書上把這件事提出來當結論:「The cancellation of the (100) reflection occurs in the bcc lattice because the planes are identical in composition.」(p.40)

這句話的價值在於它馬上給出「什麼時候規則會失效」(以下的 CsCl 例子與「體心點陣消光」這個說法是本站補的,書上沒有):如果中間那層換成不同的原子(例如 CsCl 結構:角落是 Cs、體心是 Cl),幾何上仍是 bcc 排列,但兩層散射能力不同,S(100) = f_Cs − f_Cl ≠ 0,(100) 就回來了,只是很弱。所以「bcc 消光」嚴格說應該叫「體心點陣(I 型點陣)消光」,它是點陣類型的性質,不是「原子排成 bcc 的形狀」的性質。書上最後補了一句留給讀者:「A similar cancellation can easily be found in the hcp structure.」(hcp 的 (001) 也因為中間那層 ⅓⅔½ 的原子而消光。)

回頭對照第 3 段的指標陷阱:(200) 與 (100) 不是同一件事我的判讀

這裡要把兩段扣起來。第 3 段(書本 p.32)在推導 Bragg 形式時說了一段話,當時看起來只是技術性註記,現在才看得出它的重量:

「The integers hkl that define G are not necessarily identical with the indices of an actual crystal plane, because the hkl may contain a common factor n, whereas in the definition of the indices in Chapter 1 the common factor has been eliminated.」(p.32)

也就是說,第 1 章的 Miller 指標是約分過的:(200) 這組平面在第 1 章的語言裡就是 (100)。但在繞射的語言裡,(200) 與 (100) 是兩個不同的倒晶格向量,|G(200)| 是 |G(100)| 的兩倍,對應的 d 值差一半,落在完全不同的 2θ。書上 p.32 用 n 這個公因數把兩者接起來:(24) 的 2d sinθ = nλ 裡,「d is the spacing between adjacent parallel planes with indices h/n, k/n, l/n」。

現在把這件事和 bcc 消光放在一起看,就得到本章最容易被搞混的一組事實:

問題答案理由
bcc 的 (100) 反射存在嗎?不存在(47) 給 σ = 1 為奇,S = 0
bcc 的 (200) 反射存在嗎?存在σ = 2 為偶,S = 2f
(100) 與 (200) 是同一組晶面嗎?是(Miller 指標都是 {100})差一個公因數 n = 2
那為什麼一個消光一個不消光?因為它們是不同的 G、不同的繞射階數消光看的是 G,不是晶面

換句話說,「bcc 沒有 (100)」的正確讀法不是「bcc 沒有 {100} 這組晶面」——當然有,立方晶胞的六個面就是。正確讀法是「從 {100} 這組晶面出發的一階反射消失了,但二階反射還在」。實空間的說法是:因為中間插了一層等效的原子,實際的散射面間距是 a/2 而不是 a,所以第一個出現的反射就是 d = a/2 的那個,我們把它標成 (200)。

這個陷阱在考試與讀論文時都會咬人。看到文獻寫「the (200) peak of bcc Fe」,不要以為它跟 sc 的 (200) 是同一件事:bcc 的 (200) 是它的第二個峰(第一個是 (110)),而 sc 的 (200) 是它的第四個峰(前面依序是 (100)、(110)、(111))——同一組指標在兩種點陣裡排序完全不同,指紋不同。

fcc 晶格的結構因子

S(v₁v₂v₃) = f {1 + exp[−iπ(v₂ + v₃)] + exp[−iπ(v₁ + v₃)] + exp[−iπ(v₁ + v₂)]}

(48)  p.40

f:原子形狀因子(四顆原子相同);v₁v₂v₃:以立方晶胞為參考的反射指標。

出處推導
把 fcc 的基底代入 (46)。書上 p.40 逐字給的基底是「identical atoms at 000; 0½½; ½0½; ½½0」(原文中的分數在 PDF 文字層會被拆成「1 2」漂走,此處以 PyMuPDF 放大目視還原)。四顆原子依序給出 1、exp[−i2π(½v₂+½v₃)]、exp[−i2π(½v₁+½v₃)]、exp[−i2π(½v₁+½v₂)]。
極限行為
指標全為偶數或全為奇數 → 三個指數都是 +1 → S = 4f;指標有奇有偶 → 恰好兩個指數為 −1、一個為 +1 → 1 − 1 + 1 − 1 = 0 → S = 0。與 bcc 同理,這裡的「消光」也只是「用了四倍大的立方晶胞」造成的。極限檢查:(000) 給 S = 4f = 晶胞總電子數,符合 (46) 的一般結論。
典型數值
(代入計算,書中未列)fcc 的允許反射序列 v₁²+v₂²+v₃² = 3, 4, 8, 11, 12, 16, 19, 20, …,即 (111), (200), (220), (311), (222), (400), (331), (420)。Fig. 17 的 KBr 圖樣標出的正是 (111), (200), (220), (311), (222), (400), (331), (420),與這個序列完全吻合。

fcc:為什麼一定是「全奇或全偶」補推導

書上 p.41 開頭用三句話帶過 (48) 的討論:「If all indices are even integers, S = 4f; similarly if all indices are odd integers. But if only one of the integers is even, two of the exponents will be odd multiples of −iπ and S will vanish. If only one of the integers is odd, the same argument applies and S will also vanish.」這段話是對的,但它用「只有一個是偶/只有一個是奇」來分類,第一次讀會卡住:為什麼這兩種情況就窮盡了「有奇有偶」的所有可能?這裡改用一個更乾淨的方式重算,並補上書上省略的窮舉。

  1. 先化簡符號。令 A = (−1)^(v₂+v₃),B = (−1)^(v₁+v₃),C = (−1)^(v₁+v₂),則 (48) 就是 S = f(1 + A + B + C)。注意每個都只能是 ±1。
  2. 關鍵恆等式:A·B·C = (−1)^{2(v₁+v₂+v₃)} = +1。也就是說 A、B、C 三個 ±1 的乘積恆為 +1,所以負號只能出現偶數個:要嘛三個都是 +1,要嘛恰好兩個是 −1。第三種可能(恰好一個 −1 或三個 −1)在數學上被排除。
  3. 情況一:A = B = C = +1 → S = f(1+1+1+1) = 4f。
  4. 情況二:三者中恰有兩個為 −1 → S = f(1 + 1 − 1 − 1) = 0。這就是為什麼 fcc 只有『全滅』與『全亮』兩種結果,沒有中間值。
  5. 現在把 A = B = C = +1 翻譯回指標條件:需要 v₂+v₃、v₁+v₃、v₁+v₂ 同時為偶數。兩數之和為偶 ⇔ 兩數同奇偶。三個條件合起來就是「v₁、v₂、v₃ 三者兩兩同奇偶」,即三者全奇或三者全偶。
  6. 反過來,只要有一個指標的奇偶性與另外兩個不同(無論它是唯一的偶數還是唯一的奇數),就會有兩個和變成奇數 → 兩個 −1 → S = 0。這正是書上那兩句『only one is even』『only one is odd』的內容,而用 A·B·C = 1 的論證可以一次把它們涵蓋,不必分兩種情況討論。
  7. 書上的總結逐字是:「Thus in the fcc lattice no reflections can occur for which the indices are partly even and partly odd.」
  8. 檢查一致性:全奇或全偶的三元組在所有整數三元組中恰佔 1/4(三個指標各有奇偶兩種,共 8 種奇偶組合,其中「全奇」「全偶」佔 2 種)。也就是說,用立方晶胞算出來的倒晶格點有 3/4 被消光規則刪掉,剩下的 1/4 才是 fcc 真正的倒晶格——書本 p.38 已經算過,「the bcc lattice is reciprocal to the fcc lattice」,它是 bcc 排列的點陣。消光規則沒有消滅任何物理峰,它只是在修正「我們用了一個太大的晶胞」這件事。同一個結論若直接用 fcc 的原始晶胞(一顆原子)計算,S = f 恆不為零,而允許的 G 一開始就只有那 1/4 —— 兩種算法給出完全相同的繞射圖樣。

Fig. 17:KCl 為什麼會假裝成 sc(本頁無圖,以下為文字還原)

書上 p.41 說「The point is beautifully illustrated by Fig. 17」,然後只給了一句解釋:「both KCl and KBr have an fcc lattice, but n(r) for KCl simulates an sc lattice because the K⁺ and Cl⁻ ions have equal numbers of electrons.」這裡把圖與邏輯一起補完。

Fig. 17 畫的是什麼(書本 p.42,目視判讀):上下兩格粉末 X 光繞射圖,共用一條橫軸,標「2θ」並有一個向左的箭頭——也就是 2θ 從右往左遞增,右端 20°、左端 80°,中間刻度 30°、40°、50°、60°、70°。縱軸是散射強度(未標刻度)。

上格標「KCl」,標出的峰由右(低角)到左(高角)是:(200)、(220)、(222)、(400)、(420)。(200) 遠高於其他,(220) 次之,其餘都很矮。指標全是偶數,一個奇指標的峰也沒有。

下格標「KBr」,標出的峰是:(111)、(200)、(220)、(311)、(222)、(400)、(331)、(420)。(200) 最高,(220) 次之,(111) 明顯可見但比 (200) 矮很多,(311)、(331) 是可辨認的小峰。全奇與全偶的反射都在。

圖說(書本 p.42 逐字):「Comparison of x-ray reflections from KCl and KBr powders. In KCl the numbers of electrons of K⁺ and Cl⁻ ions are equal. The scattering amplitudes f(K⁺) and f(Cl⁻) are almost exactly equal, so that the crystal looks to x-rays as if it were a monatomic simple cubic lattice of lattice constant a/2. Only even integers occur in the reflection indices when these are based on a cubic lattice of lattice constant a. In KBr the form factor of Br⁻ is quite different to that of K⁺, and all reflections of the fcc lattice are present. (Courtesy of R. van Nordstrand.)」

書上正文的說法更簡短:「both KCl and KBr have an fcc lattice, but n(r) for KCl simulates an sc lattice because the K⁺ and Cl⁻ ions have equal numbers of electrons.」(p.41)這句話的完整推導書上沒寫,下一個單元補上。

補上 KCl/KBr 的結構因子代數補推導

書上只說「K⁺ 與 Cl⁻ 電子數相同,所以 KCl 看起來像 a/2 的 sc」,但沒有把 NaCl 型結構的 S 算出來。這一步其實很短,而且算完以後「為什麼是 a/2」「為什麼 KBr 的奇指標會回來」都一目了然。這也是本章第一次遇到「兩種原子的基底」,值得完整走一遍。

  1. NaCl 型結構 = fcc 點陣 + 兩原子基底。取陽離子在分數座標 (0,0,0),陰離子在 (½,0,0)(兩者相距 a/2)。以立方晶胞計,共 8 顆原子:陽離子在 000, 0½½, ½0½, ½½0;陰離子在 ½00, 0½0, 00½, ½½½。
  2. 把這 8 項代入 (46)。注意陰離子的四個位置,正好是陽離子四個位置各加上 (½,0,0)。所以陰離子那四項,是陽離子四項乘上共同因子 exp[−i2π(½v₁)] = exp(−iπv₁) = (−1)^{v₁}。
  3. 提出公因式:S(v₁v₂v₃) = [f₊ + f₋(−1)^{v₁}] × S_fcc(v₁v₂v₃),其中 S_fcc 就是 (48) 的大括號部分(全奇或全偶時為 4,否則為 0)。
  4. 先看 S_fcc:它把「有奇有偶」的反射全部殺掉,所以只剩全奇與全偶兩族——這部分與原子種類無關,純粹是 fcc 點陣的貢獻。
  5. 再看方括號。在允許的反射裡指標同奇偶,所以 (−1)^{v₁} 有兩個值:全偶 → (−1)^{v₁} = +1 → S = 4(f₊ + f₋);全奇 → (−1)^{v₁} = −1 → S = 4(f₊ − f₋)。
  6. KCl:K 的原子序 19,失去一個電子後 K⁺ 有 18 個電子;Cl 的原子序 17,得到一個電子後 Cl⁻ 也有 18 個電子。電子數相同 ⇒ f₊ ≈ f₋(書上逐字:「The scattering amplitudes f(K⁺) and f(Cl⁻) are almost exactly equal」)⇒ 全奇那一族 S = 4(f₊ − f₋) ≈ 0,整族消失,圖上只剩全偶指標。
  7. 為什麼剩下的等價於 a/2 的 sc:全偶指標可以寫成 (2h, 2k, 2l)。由 (15),|G| ∝ √(v₁²+v₂²+v₃²)/a,代入得 |G| ∝ 2√(h²+k²+l²)/a = √(h²+k²+l²)/(a/2)。這正是晶格常數為 a/2 的簡單立方晶格的倒晶格。而且既然 f₊ ≈ f₋,這些「等效 sc 格點」上的散射體還全都一樣——所以 X 光「看到」的確實是一個單原子的 sc 晶體,書上說 「the crystal looks to x-rays as if it were a monatomic simple cubic lattice of lattice constant a/2」。
  8. KBr:Br 原子序 35,Br⁻ 有 36 個電子,與 K⁺ 的 18 個差很多 ⇒ f₊ − f₋ 明顯不為零⇒ 全奇反射回來,(111)、(311)、(331) 出現。書上逐字:「In KBr the form factor of Br⁻ is quite different to that of K⁺, and all reflections of the fcc lattice are present.」
  9. 強度比的量級估計(代入計算,書中未列):以 G → 0 的極限值代替 f,KBr 的 (111) 強度 ∝ |36 − 18|² = 324,(200) 強度 ∝ |36 + 18|² = 2916,比值約 1 : 9。Fig. 17 下格的 (111) 確實明顯低於 (200),量級吻合。(真實比值還要乘上多重度、Lorentz-偏極因子與 f 隨角度的衰減,這裡只取數量級。)

從 KCl 學到的一般原則:硬消光與軟消光我的判讀

KCl/KBr 這個例子的價值遠超過「一個有趣的巧合」。它示範了消光規則的兩種等級,而研究生常常把兩者混為一談。

硬消光(幾何消光):bcc 的 σ 為奇、fcc 的指標有奇有偶。這兩條規則裡 f 整個被提到括號外,剩下的相位和恆為零,所以這些反射的 S 嚴格等於零——抵消是點陣中心化保證的,與原子是什麼、與溫度、與波長都無關,實驗上這些位置只有背景。(鑽石的 (222) 常被拿來當硬消光的例子,但它其實不是,理由見下一段。)它們對應的是「我們選了一個非原始的晶胞」這件事,所以硬消光的完整清單可以只從空間群讀出來——這就是國際結晶表(International Tables)裡每個空間群都附一份 reflection conditions 的原因。

軟消光(對比消光):KCl 的全奇反射。S = 4(f₊ − f₋) 並不是零,只是因為兩種離子的電子數恰好相同而幾乎為零。這種消光:

  • 會隨散射角變化(f₊ 與 f₋ 的角度依賴略有不同,高角度時差異會露出來);
  • 換一種輻射就可能完全消失(中子散射自原子核,K 與 Cl 的中子散射長度完全不同,用中子看 KCl 立刻就看得到全奇反射);
  • 在吸收邊附近用異常散射可以刻意放大(f 變成複數,f = f₀ + f′ + if″)。

球對稱近似造成的零:鑽石的 (222)。它為零靠的不是點陣中心化,而是 (50) 背後的球對稱原子假設(原文 p.41「If the electron distribution is spherically symmetric about the origin」)。在這個近似下 S 嚴格為零;但四面體鍵結讓價電子分布明顯不是球對稱的,一旦承認後面那段要講的鍵結電荷重分布,(222) 就不再嚴格為零——X 光實測確實會出現一個極弱的殘留峰,這正是後面說「不是看不到、只是不敏感」的最佳例子。書上 Fig. 18 用的是中子(散射自點狀的原子核,沒有形狀因子的方向性),那裡的 (222) 才真的是零。

判準:看到一族反射「不見了」,先問是幾何(中心化)保證的零、球對稱近似造成的零,還是數值巧合的零。第一種查空間群就有答案,後兩種都要算 f。這個區分在做薄膜 XRD、超晶格繞射、或判斷有序/無序合金(例如 Cu₃Au 的 superlattice 峰)時是第一步。Problem 7 的第 (c) 小題「Explain what happens if f_A = f_B」問的就是這件事的一維版本。

原子形狀因子 f:它是一個「比值」,不是一個「量」

(42) 已經把 fⱼ 定義出來了((43) 只是把它組進 S_G),但那時它只是「括號裡那個積分」。到了 p.41 的 Atomic Form Factor 小節,書上才回頭說清楚它是什麼。開頭逐字寫:「In the expression (46) for the structure factor, there occurs the quantity fⱼ, which is a measure of the scattering power of the jth atom in the unit cell. The value of f involves the number and distribution of atomic electrons, and the wavelength and angle of scattering of the radiation.」

注意這句話裡列的四個依賴:電子數目、電子分布、輻射波長、散射角度。後兩者其實只透過 G = 4π sinθ/λ 這一個組合進入(由 (23) 與 (24) 可得),所以 f 實質上只是 G 的函數:f = f(G)。這就是為什麼結構因子表(international tables)裡的 f 值都列成 sinθ/λ 的函數,而不是分別列波長與角度。(這一句與下面「f 是比值」的解讀為本站補充,書上沒有。)

接著書上宣告「We now give a classical calculation of the scattering factor」。classical 這個詞值得停一下:這裡把電子當成連續的電荷雲,每一小塊 n(r)dV 都獨立地散射入射波,再把振幅相加。這是 Thomson 散射的疊加,沒有用到任何量子力學——電子的量子性只進入 n(r) 的形狀(那要解薛丁格方程式才知道),但散射本身的處理是純粹的波干涉。書上緊接著說:「The scattered radiation from a single atom takes account of interference effects within the atom.」——干涉發生在原子內部,這正是 f 隨角度衰減的原因。

最重要的一句話在 p.42 (51) 之後:「Therefore f is the ratio of the radiation amplitude scattered by the actual electron distribution in an atom to that scattered by one electron localized at a point.」

把這句話拆開(以下為本站補充,非書上內容):分母是「一個集中在一點的電子」的散射振幅。點狀電子沒有內部結構,散射到任何角度的振幅都一樣(Thomson 散射的角分布另外處理)。分子是真實原子的散射振幅。所以 f 是一個無因次的比值,單位是「等效點電子的個數」。f = 20 的意思是「這顆原子在這個 G 之下,散射得像 20 個集中在一點的電子」。這個讀法讓 (51) 的 f = Z 變得幾乎是同義反覆:在 G = 0(前向)時原子內沒有相位差,Z 個電子全部同相,當然等效於 Z 個點電子。

(49) → (50):球對稱角積分的每一步補推導

書上 p.41–42 把這個積分做完了,但中間有兩步跳得比較快:一是體積元素怎麼寫成 2π dr r² d(cos α)(2π 是怎麼來的),二是最後那個指數差怎麼變成 sin。這兩步都很短,但不補出來會卡住。

先交代設定。書上 p.41 逐字寫:「We defined the form factor in (42)」,然後重列一次成 (49):fⱼ = ∫ dV nⱼ(r) exp(−iG·r)。接著「Let r make an angle α with G; then G·r = Gr cos α. If the electron distribution is spherically symmetric about the origin, then …」——球對稱是唯一的物理假設,它對閉殼層原子是好近似,對有方向性的價電子則不然。

  1. 取球座標,極軸沿 G 的方向。體積元素 dV = r² dr sinα dα dφ。因為 nⱼ(r) 只依賴 r、指數只依賴 r 與 α,被積函數與方位角 φ 無關,φ 的積分直接給出 2π。這就是 (49) 之後那個式子裡 2π 的來源。
  2. 把 sinα dα 換成 −d(cos α)。當 α 從 0 跑到 π 時,cos α 從 +1 跑到 −1;負號與積分上下限對調互相抵消,所以 ∫₀^π sinα dα → ∫₋₁^{+1} d(cos α)。書上把這件事寫在 p.42 開頭:「after integration over d(cos α) between −1 and 1」。
  3. 於是 fⱼ = 2π ∫₀^∞ dr r² nⱼ(r) ∫₋₁^{+1} d(cos α) exp(−iGr cos α)。這就是書上緊接 (49) 印出的那一行(書上用的是 ≡ 而非 =)。
  4. 算內層積分。令 u = cos α,這是一個初等積分:∫₋₁^{+1} du e^{−iGru} = [e^{−iGru}/(−iGr)]₋₁^{+1} = (e^{−iGr} − e^{+iGr})/(−iGr) = (e^{iGr} − e^{−iGr})/(iGr)。
  5. 代回去:fⱼ = 2π ∫ dr r² nⱼ(r) · (e^{iGr} − e^{−iGr})/(iGr)。這正是書上 p.41 最後印的那一行。
  6. 用 Euler 公式 e^{iθ} − e^{−iθ} = 2i sin θ:(e^{iGr} − e^{−iGr})/(iGr) = 2i sin(Gr)/(iGr) = 2 sin(Gr)/(Gr)。i 在這裡完全消掉了——這是球對稱的直接後果:球對稱的電子雲有反演對稱,所以它的傅立葉變換是實數。 這也回應了前面 (46) 之後那句「S need not be real」:S 的複數性完全來自原子的擺放位置,不來自原子本身。
  7. 把 2 乘進 2π 得 4π:fⱼ = 4π ∫ dr nⱼ(r) r² · sin(Gr)/(Gr),即 (50)。
  8. 校註:書本 p.42 的 (50) 式,在 r² 與 sin(Gr)/(Gr) 之間印成了一個等號(原文字面是 fⱼ = 4π ∫ dr nⱼ(r)r² = sin Gr / Gr),但右邊的 r 是左邊積分的啞變數,寫成等號等於把束縛變數和自由變數混在一起,式子在數學上無定義;而且由 (51),左邊這個積分本身就等於 Z,不可能等於 sin Gr/(Gr)。此處應為乘號。已用 PyMuPDF 放大目視確認不是抽取誤差,是書上的排版錯誤。對照 p.41 最後一行同樣位置印的是「·」,可以確認原意是相乘。

球對稱原子的形狀因子

fⱼ = 4π ∫ dr nⱼ(r) r² · sin(Gr)/(Gr)

(50)  p.42

nⱼ(r):球對稱的徑向電子濃度(單位體積電子數);r:離原子核的距離;G = |G| = 4π sinθ/λ;4π r² dr:半徑 r 處厚度 dr 的球殼體積。

出處推導
由 (49) 在球對稱假設下對角度積分而得,過程見上一個推導單元。書上原文在 r² 與 sin Gr/Gr 之間印成等號,但右邊的 r 是左邊積分的啞變數,寫成等號在數學上無定義(且由 (51),左邊的積分等於 Z),此處依 p.41 末行的寫法校正為乘號。
極限行為
Gr ≪ 1(低角度或緻密電子雲):sin(Gr)/(Gr) → 1 − (Gr)²/6,所以 f ≈ Z[1 − G²⟨r²⟩/6]。這個展開式告訴你 f 在小角度的下降速率由電子雲的均方半徑 ⟨r²⟩ 決定——f 曲線的初始斜率就是原子大小的量度。價電子分布得遠(⟨r²⟩ 大),所以它們的貢獻在很小的 G 就掉光了;內層電子緊縮在核附近,一路撐到高 G。這是理解下一段「X 光對價電子重分布不敏感」的關鍵。 Gr ≫ 1:sin(Gr)/(Gr) 劇烈振盪並以 1/Gr 衰減,積分趨近零。物理上:散射向量的波長 2π/G 比電子雲小得多時,原子不同部位的散射彼此完全相消。 G = 0:退化成 (51) 的 f = Z。
典型數值
(代入計算,書中未列)Cu Kα λ = 1.5418 Å、2θ = 60° ⇒ G = 4π sin30°/1.5418 = 4.07 Å⁻¹。對一個 ⟨r²⟩^{1/2} ≈ 0.5 Å 的內層電子雲,Gr ≈ 2,sin(Gr)/(Gr) ≈ 0.45——仍有可觀貢獻;對 ⟨r²⟩^{1/2} ≈ 2 Å 的價電子雲,Gr ≈ 8,sin(Gr)/(Gr) ≈ 0.12,而且在 r 由 0 積到 2 Å 的範圍內 sin(Gr) 已經跑完一個以上的完整週期、正負相消,淨貢獻幾乎歸零。

前向散射極限:f = Z

fⱼ = 4π ∫ dr nⱼ(r) r² = Z

(51)  p.42

Z:原子(或離子)的電子總數。注意是電子數不是原子序——對離子而言 Z 要加減電荷數,這正是 KCl 裡 K⁺ 與 Cl⁻ 都等於 18 的原因。

出處推導
在 (50) 中取 Gr → 0 的極限。書上 p.42 的說法是:「If the same total electron density were concentrated at r = 0, only Gr = 0 would contribute to the integrand. In this limit (sin Gr)/Gr = 1, and …」注意書上是用「把電子全部集中到原點」這個思想實驗來引入的,不是單純取數學極限——這樣讀者才看得出 f 是一個「與點電子相比」的比值。
極限行為
G → 0 對應 θ → 0,即前向散射。書上補一句「In the forward direction G = 0, and f reduces again to the value Z」。物理意義:前向沒有任何路徑差,原子裡每一個電子都同相散射,所以振幅就是 Z 個電子振幅的直接相加。注意這是一個上界:對所有 G,|f(G)| ≤ f(0) = Z,因為任何非零的 G 都會引入相位差而造成部分相消。實務上 f(0) = Z 這個關係是形狀因子表的正規化條件,也是快速檢查自己算出來的 f(G) 是否合理的第一個檢驗。
典型數值
K⁺:19 − 1 = 18;Cl⁻:17 + 1 = 18(兩者相等,故 KCl 的全奇反射消失);Br⁻:35 + 1 = 36(與 K⁺ 差一倍,故 KBr 的全奇反射明顯存在)。這三個數字就是 Fig. 17 兩張圖差異的全部原因。

「X 光看不到鍵結」——這句話要說到什麼程度才算對

p.42 的最後一段是整節最需要小心讀的三句話。先逐字抄下來:

「The overall electron distribution in a solid as seen in x-ray diffraction is fairly close to that of the appropriate free atoms. This statement does not mean that the outermost or valence electrons are not redistributed somewhat in forming the solid; it means only that the x-ray reflection intensities are represented well by the free atom values of the form factors and are not very sensitive to small redistributions of the electrons.」

Kittel 在這裡做了一件很克制的事:他先給出一個容易被誤讀的命題(固體裡的電子分布跟自由原子差不多),然後立刻用一整句話把誤讀擋掉(這不表示價電子沒有重新分布),最後才給出真正的命題(反射強度用自由原子的 f 值就描述得很好,而且對電子的小幅重分布不很敏感)。

坊間常聽到的「X 光看不到鍵結電子」是對這段話的過度簡化。原文說的是 not very sensitive(不很敏感),不是 blind(看不到)。這個差別很重要:

為什麼不敏感。(以下到本段末為本站補充,書上沒有這些內容。)三個層次的理由:

(1) 數量:價電子只佔總電子數的一小部分。Si 有 14 個電子,只有 4 個是價電子,佔 29%;對過渡金屬比例更低。

(2) 分布:由 (50) 的小 G 展開 f ≈ Z[1 − G²⟨r²⟩/6],價電子的 ⟨r²⟩ 大,所以它們對 f 的貢獻集中在極小的 G。而 Bragg 反射只發生在離散的 G 上,最小的那個 G 通常已經大到價電子的貢獻衰減得差不多了。價電子的資訊主要落在 G 很小的地方,但那裡沒有 Bragg 峰可以量。

(3) 強度:S 的平方使得小的修正更不顯眼——f 變 2% 只造成強度變 4%,而實驗誤差、吸收修正、Debye–Waller 修正都是這個量級。

但它確實看得到,只是要很努力。 現代的電荷密度研究(charge-density study)就是靠量測極高精度的低角度反射、扣掉自由原子模型、看殘差圖(deformation density map)上的鍵結電荷堆積。這種實驗需要低溫(壓抑 Debye–Waller 因子)、同步輻射、以及非常小心的吸收與消光修正。所以正確的說法是:「在一般的結構解析精度下,用自由原子 f 值就夠了;要看到鍵結電荷需要專門設計的高精度實驗。」

對做半導體與元件的人這一點有實際意義:如果你想從 XRD 判斷 Si 的鍵結狀態或缺陷的電子結構,這條路很難走;XRD 擅長的是幾何(晶格常數、應變、取向、晶粒大小、缺陷密度),電子結構要交給光譜(XPS、EELS、XANES)。

SUMMARY:四個等價的繞射條件,其實是同一句話的四種語言

書本 p.43 的 SUMMARY 把整章壓縮成七個項目符號。這一節不是背誦清單,它真正的價值在於讓你看清楚:前面推了十幾頁的東西,其實只是同一個條件的四種說法。以下先逐字列出書上的內容,再解釋每一種說法的適用場合。

書上前兩項逐字為:

「• Various statements of the Bragg condition: 2d sin θ = nλ ; Δk = G ; 2k·G = G².」

「• Laue conditions: a₁·Δk = 2πv₁ ; a₂·Δk = 2πv₂ ; a₃·Δk = 2πv₃.」

說法式號語言什麼時候用它
2d sin θ = nλ(24)實空間、晶面、角度手算 2θ、指標化粉末圖、判斷某個峰落在哪
Δk = G(21)倒空間、向量推導、證明、任何要談「動量轉移」的場合
2k·G = G²(23)倒空間、純量導出 Bragg 形式、Ewald 球作圖的代數版本
aᵢ·Δk = 2πvᵢ(25)三個錐面的交線單晶繞射的幾何直觀、Laue 法
k·(½G) = (½G)²(26)Brillouin zone 邊界能帶理論,第 7 章以後全部用這個

它們為什麼等價。 從 Δk = G 出發:

彈性散射要求 |k′| = |k|,而 k′ = k + Δk = k + G,所以 (k+G)² = k²,展開得 2k·G + G² = 0,即 (22)。把 G 換成 −G(倒晶格對原點對稱,−G 也是倒晶格向量)得到 (23) 的 2k·G = G²。再用 Problem 1 的結果 d(hkl) = 2π/|G| 與 k = 2π/λ 代入,就得到 2d sinθ = λ,考慮公因數 n 之後成為 (24) 的 2d sinθ = nλ。至於 Laue 三式,是把 Δk = G 分別與 a₁, a₂, a₃ 做內積,再用 bᵢ·aⱼ = 2πδᵢⱼ——又是 (14)。

注意這四種說法都只回答「峰在哪裡」。 這一段前半(結構因子)回答的是「峰有多亮」,SUMMARY 用另一個項目符號單獨列出:「• The scattered amplitude in the direction k′ = k + Δk = k + G is proportional to the geometrical structure factor: S_G ≡ Σ fⱼ exp(−irⱼ·G) = Σ fⱼ exp[−i2π(xⱼv₁ + yⱼv₂ + zⱼv₃)],where j runs over the s atoms of the basis, and fⱼ is the atomic form factor (49) of the jth atom of the basis.」

「位置」與「強度」這兩件事分屬晶格與基底,是整章最該帶走的結構性結論。

SUMMARY:三個立方晶格的第一布里淵區

SUMMARY 最後給了一張三行的小表,它是第 4、5 段那些幾何計算的結晶。書上逐字列的是:

Crystal latticeFirst Brillouin zone
Simple cubicCube
Body-centered cubicRhombic dodecahedron (Fig. 13)
Face-centered cubicTruncated octahedron (Fig. 15)

這張表要配著「誰是誰的倒晶格」一起記,否則很容易記反:書本 p.36 說「an fcc lattice is the reciprocal lattice of the bcc lattice」,p.38 說「the bcc lattice is reciprocal to the fcc lattice」。於是:

實空間晶格倒晶格第一 BZ = 倒晶格的 Wigner–Seitz 胞形狀
sc(常數 a)sc(常數 2π/a)sc 的 WS 胞立方體,邊長 2π/a
bcc(立方邊 a)fccfcc 的 WS 胞菱形十二面體(12 面)
fcc(立方邊 a)bccbcc 的 WS 胞截角八面體(8 個六邊形 + 6 個正方形)

常犯的錯:把「bcc 的 BZ」想成「bcc 的 WS 胞」。不是。BZ 是倒晶格的 WS 胞,而 bcc 的倒晶格是 fcc,所以 bcc 的 BZ 是 fcc 點陣的 WS 胞(菱形十二面體)。書本 p.37 把這件事講得很明確:bcc 的 BZ「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」。

SUMMARY 還附了一句對後面幾章極其重要的話:「Only waves whose wavevector k drawn from the origin terminates on a surface of the Brillouin zone can be diffracted by the crystal.」以及那條把整章的傅立葉觀點濃縮成一行的敘述:「Any function invariant under a lattice translation T may be expanded in a Fourier series of the form n(r) = Σ_G n_G exp(iG·r).」——請注意這句話的主詞是「any function」,不是「electron concentration」。電子密度只是第一個例子;晶體中的位勢 U(r)(第 7 章能帶)、力常數、磁化率都可以這樣展開。第 2 章真正教會你的是一個方法,不是一個現象。

Problems 1–3:三題純幾何,但第 1 題是本章正文的欠債

書本 p.43–44 的前三題都是倒晶格的幾何練習,其中第 1 題特別重要——正文在 p.32 推 Bragg 形式時已經先用了它的結果,卻把證明丟給習題。如果你在 p.32 覺得「d(hkl) = 2π/|G| 是哪來的」,答案就在這裡。

(說明:以下題目敘述為原文逐字,書上只給題目不給解答;每題後面的解題要點與評註是本站補的,不是書上的內容。)

Problem 1. Interplanar separation.(逐字)「Consider a plane hkl in a crystal lattice. (a) Prove that the reciprocal lattice vector G = hb₁ + kb₂ + lb₃ is perpendicular to this plane. (b) Prove that the distance between two adjacent parallel planes of the lattice is d(hkl) = 2π/|G|. (c) Show for a simple cubic lattice that d² = a²/(h²+k²+l²).」

這一題把第 1 章的 Miller 指標與第 2 章的倒晶格焊接在一起。

(a) hkl 平面與三個軸交於 a₁/h、a₂/k、a₃/l。取平面內的兩個向量 u = a₁/h − a₂/k 與 v = a₂/k − a₃/l,用 (14) 的 bᵢ·aⱼ = 2πδᵢⱼ 算:G·(a₁/h) = (hb₁+kb₂+lb₃)·(a₁/h) = (1/h)(h·2π) = 2π,同理 G·(a₂/k) = 2π、G·(a₃/l) = 2π。所以 G·u = 2π − 2π = 0、G·v = 0。G 垂直於平面內兩個不平行的向量,故垂直於該平面。

(b) 相鄰平面的間距等於「從原點到第一個平面」沿 G 方向的投影:d = (a₁/h)·Ĝ = G·(a₁/h)/|G| = 2π/|G|。

(c) 對 sc,bᵢ = (2π/a)x̂ᵢ,故 |G|² = (2π/a)²(h²+k²+l²),代入 d = 2π/|G| 得 d² = a²/(h²+k²+l²)。

(c) 的結果是粉末繞射指標化的核心公式:sin²θ ∝ (h²+k²+l²),前面 bcc/fcc 的「特徵序列」就是靠它把 2θ 翻譯成指標。

Problem 2. Hexagonal space lattice. 書上給的原始平移向量是a₁ = (3^{1/2}a/2)x̂ + (a/2)ŷ;a₂ = −(3^{1/2}a/2)x̂ + (a/2)ŷ;a₃ = cẑ。要證 (a) 原始胞體積 = (3^{1/2}/2)a²c;(b) 倒晶格的原始平移是b₁ = (2π/3^{1/2}a)x̂ + (2π/a)ŷ;b₂ = −(2π/3^{1/2}a)x̂ + (2π/a)ŷ;b₃ = (2π/c)ẑ,書上逐字說「so that the lattice is its own reciprocal, but with a rotation of axes」;(c) 描述並畫出六方晶格的第一 BZ(是一個正六角柱)。「六方的倒晶格還是六方,但轉了 30°」這件事在第 7 章畫石墨烯/六方半導體的能帶時天天用到。

Problem 3. Volume of Brillouin zone. 要證 BZ 體積 = (2π)³/V_c,V_c 是實空間原始胞體積。書上給的提示是「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」。(校註:原文這個恆等式中間的點積符號,在 PDF 文字層被解碼成 ∑,對照書中同一符號在 p.30、p.31 的其他出現位置可確認是 ·。)這一題的結論是本書最常被默默使用的定理之一:BZ 體積固定,所以晶體裡的 k 點密度(每單位 k 空間的容許態數)只跟晶體體積有關,與晶格細節無關。第 6 章數自由電子態、第 5 章數聲子模態,用的都是這條。

Problem 4:峰寬 ∝ 1/M,與「為什麼 Bragg 峰這麼窄」補推導

這是本章最有物理內容的一題,也是 p.31 那句「It is a simple exercise (Problem 4) to show that F is negligibly small when Δk differs significantly from any reciprocal lattice vector」的兌現。書上把題目設得很完整(給了級數和、給了要證的答案),但中間的代數還是留給你。以下逐步補齊。

題目設定(書本 p.44 逐字):「We suppose that in a linear crystal there are identical point scattering centers at every lattice point ρ_m = ma, where m is an integer. By analogy with (20), the total scattered radiation amplitude will be proportional to F = Σ exp[−ima·Δk].」

  1. 把和寫成幾何級數。令 x ≡ exp(−i a·Δk),則 F = Σ_{m=0}^{M−1} x^m。用書上給的公式 Σ_{m=0}^{M−1} x^m = (1 − x^M)/(1 − x) 得F = (1 − exp[−iM(a·Δk)]) / (1 − exp[−i(a·Δk)])。(校註:書上這一式的分子少印一個右括號,印成 1 − exp[−iM(a·Δk];已放大目視確認是書上的排版錯誤,不是抽取誤差。)
  2. (a) 要算 |F|² = F*F。關鍵的小恆等式(書上沒寫):|1 − e^{iφ}|² = (1 − cos φ)² + sin²φ = 2 − 2cos φ = 4 sin²(φ/2)。(最後一步用倍角公式 1 − cos φ = 2sin²(φ/2)。)
  3. 把這個恆等式分別用在分子(φ = M a·Δk)與分母(φ = a·Δk):|F|² = 4sin²(½M a·Δk) / 4sin²(½ a·Δk) = sin²(½M(a·Δk)) / sin²(½(a·Δk))。這就是書上 (a) 要證的結果。這個函數叫 Fejér 核或「M 縫干涉函數」。
  4. 先看峰高。當 a·Δk → 2πh 時,分母 sin²(½(2πh+ε)) = sin²(ε/2) → (ε/2)²,分子 sin²(½M(2πh+ε)) = sin²(πMh + Mε/2) = sin²(Mε/2) → (Mε/2)²。(注意 Mh 是整數,πMh 可以整個丟掉。)比值 → M²。峰高正比於 M²,這對應 (39) 裡那個 F_G = N S_G 的 N。
  5. (b) 求峰寬。書上定義(p.44 逐字,ε 在 PDF 文字層被解碼成 ώ,此處還原):「We change Δk slightly and define ε in a·Δk = 2πh + ε such that ε gives the position of the first zero in sin ½M(a·Δk).」
  6. 第一個零點:sin(Mε/2) = 0 且 ε ≠ 0 ⇒ Mε/2 = π ⇒ ε = 2π/M。得證。
  7. 物理翻譯。ε 是 a·Δk 的偏離量,所以沿 a 方向的波向量偏離是Δk 的分量改變 ε/a = 2π/(Ma) = 2π/L,其中 L = Ma 是晶體在該方向的長度。繞射峰在倒空間的寬度 = 2π / 晶體尺寸。 這是傅立葉變換的普遍性質(有限長度 ↔ 有限頻寬),在不同脈絡下叫不確定性原理、叫 Scherrer 公式、叫 finite-size broadening。
  8. 數量級(代入計算,書中未列):a = 3 Å 的晶體、M = 10⁵(L = 30 μm),相對峰寬 ε/2π = 1/M = 10⁻⁵。換算成角度要對 Bragg 定律微分,Δθ = (Δd/d)·tan θ,以 θ = 15° 計得 Δ(2θ) ≈ 5×10⁻⁶ rad ≈ 3×10⁻⁴ 度。這就是書上說「can be extremely narrow for macroscopic values of M」的意思,也解釋了為什麼實驗上量到的峰寬幾乎全部來自儀器解析度與晶粒/缺陷,而不是這個本質下限。書上最後補一句「The same result holds true for a three-dimensional crystal.」
  9. 與 Fig. 18 對照(本項為依圖目視讀數所做的驗算,書中沒有這段):鑽石中子繞射的峰寬標在圖上是 0.45°、0.65°、0.70°、0.80°,比上面估的 3×10⁻⁴ 度大了三個數量級以上。所以 Fig. 18 的峰寬完全不是有限尺寸效應,而是中子繞射儀的角度解析度(以及粉末晶粒的貢獻)。這是一個很好的提醒:Problem 4 給的是下限,不是預測值。

Problem 5:鑽石結構的 4n 規則,與 Fig. 18 的逐峰驗算補推導

這一題是整章消光規則的集大成,而且結論對做半導體的人直接有用:Si 與 Ge 都是鑽石結構,它們的 XRD 圖樣缺席清單就是這一題的答案。

題目(書本 p.44 逐字):「The crystal structure of diamond is described in Chapter 1. The basis consists of eight atoms if the cell is taken as the conventional cube. (a) Find the structure factor S of this basis. (b) Find the zeros of S and show that the allowed reflections of the diamond structure satisfy v₁ + v₂ + v₃ = 4n, where all indices are even and n is any integer, or else all indices are odd (Fig. 18). (Notice that h, k, l may be written for v₁, v₂, v₃ and this is often done.)」

  1. 鑽石結構 = fcc 點陣 + 兩原子基底,第二顆位在 (¼, ¼, ¼)。以立方晶胞計就是 8 顆原子:fcc 的四個位置 000, 0½½, ½0½, ½½0,以及它們各加 (¼,¼,¼) 後的四個位置 ¼¼¼, ¼¾¾, ¾¼¾, ¾¾¼。(PDF 文字層會把 ¼¼¼ 抽成「1 4 1 4 1 4」並漂到段末,需依上下文重組。)
  2. 如同 KCl 的算法,後四顆是前四顆乘上共同相位因子 exp[−i2π(¼)(v₁+v₂+v₃)]。(a) 的答案:S = S_fcc × [1 + exp(−iπσ/2)],其中 σ ≡ v₁+v₂+v₃,S_fcc 是 (48) 的大括號(全奇或全偶為 4,否則為 0)。
  3. 注意 exp(−iπ/2) = −i,所以 exp(−iπσ/2) = (−i)^σ。整個結構因子是S = f · S_fcc · [1 + (−i)^σ](八顆原子相同,f 可提出)。
  4. (b) 逐種情形檢查。第一層篩選:S_fcc ≠ 0 需要指標全奇或全偶。有奇有偶者一律消光,這一層與 fcc 完全相同。
  5. 全奇的情形:σ = 奇數。(−i)^σ = ±(−i),於是 1 + (−i)^σ = 1 ∓ i,模長 √2 ≠ 0。全奇指標一律允許,|S| = 4f√2,|S|² = 32f²。
  6. 全偶的情形:σ = 偶數,(−i)^σ = (−1)^{σ/2}。若 σ ≡ 0 (mod 4) → (−1)^{σ/2} = +1 → 1 + 1 = 2 → 允許,|S| = 8f,|S|² = 64f²。若 σ ≡ 2 (mod 4) → (−1)^{σ/2} = −1 → 1 − 1 = 0 → 消光。
  7. 把兩層合起來就是書上要的結論:允許的反射是「全奇指標」或「全偶且 σ = 4n」。最有名的犧牲者是 (222)(全偶、σ = 6 ≡ 2 mod 4)與 (200)(全偶、σ = 2)。「有 (111) 但沒有 (222)」是鑽石結構在繞射圖上的簽名。
  8. 與 Fig. 18 對照(書本 p.45,以下為依圖目視判讀):圖是粉末鑽石的中子繞射圖樣,橫軸「Counter position 2θ」從 20° 到 75°,縱軸「Neutron intensity, counts per minute」0 到 800,背景約 50 cpm。標出的峰只有四個:(111) 在 2θ ≈ 30°,峰高約 620;(220) 在 ≈ 48°,約 520;(311) 在 ≈ 58°,約 370;(400) 在 ≈ 72°,約 120。每個峰旁用一對相向的箭頭標出寬度:0.45°、0.65°、0.70°、0.80°。圖說逐字:「Neutron diffraction pattern for powdered diamond. (After G. Bacon.)」
  9. 逐峰檢查:(111) σ=3 全奇 → 允許 ✓;(220) σ=4 全偶且 4n → 允許 ✓;(311) σ=5 全奇 → 允許 ✓;(400) σ=4 全偶且 4n → 允許 ✓。而 (200)(σ=2)與 (222)(σ=6)在圖上確實找不到,與規則一致。
  10. 進一步的強度驗算(本項為本人依 Fig. 18 目視讀數所做的驗算,書中沒有這段):粉末圖的積分強度 I ∝ m · |S|² · L,m 是多重度({111}: 8、{220}: 12、{311}: 24、{400}: 6),L = 1/(sin²θ cosθ) 是 Lorentz 因子。中子散射長度與角度無關,所以不必管 f 的衰減。算出 m|S|²L 的相對值為 (111) : (220) : (311) : (400) ≈ 2.9 : 3.7 : 2.7 : 1——預測 (220) 最強、(111) 與 (311) 居中且接近、(400) 明顯最弱。圖上讀數(峰高 × 峰寬)給出的排序與相對關係確實是 (220) ≳ (111) > (311) ≫ (400),排序完全相符;至於數值比例,未扣背景時約 2.9 : 3.5 : 2.7 : 1、扣掉約 50 cpm 背景後約 4.6 : 5.5 : 4.0 : 1,可見從縮圖目視讀數的不確定度足以讓比值變動一倍以上,只能說量級與排序吻合,不能當成定量驗證。(另外也未計 Debye–Waller 與吸收修正。)即使如此,這個對照已經足以確認 (b) 算出的「全偶 |S|² = 64f²、全奇 |S|² = 32f²」這個因子 2 不是算錯:若把它拿掉,(111) 與 (311) 的預測值會各翻一倍,(111) 就會遠超過 (220),與圖不符。

Problem 6:氫原子的形狀因子——唯一能手算到底的例子補推導

(50) 看起來簡單,但對真實原子這個積分沒有解析解(要數值解 Hartree–Fock 的 n(r))。氫原子基態是唯一可以一路手算到封閉形式的情形,所以這一題的價值在於讓你親眼看到 f(G) 的長相。

題目(書本 p.44 逐字):「For the hydrogen atom in its ground state, the number density is n(r) = (πa₀³)^{−1} exp(−2r/a₀), where a₀ is the Bohr radius. Show that the form factor is f_G = 16/(4 + G²a₀²)².」

  1. 先檢查正規化:∫ dV n = 4π ∫₀^∞ dr r² (πa₀³)^{−1} e^{−2r/a₀} = (4/a₀³) ∫₀^∞ r² e^{−2r/a₀} dr = (4/a₀³)·(2!/(2/a₀)³) = (4/a₀³)·(a₀³/4) = 1。Z = 1,符合 (51)。 這一步值得做,因為它確認了題目給的係數,也順便練了 ∫ r^n e^{−pr} dr = n!/p^{n+1}。
  2. 代入 (50):f_G = 4π ∫₀^∞ dr r² · (πa₀³)^{−1} e^{−2r/a₀} · sin(Gr)/(Gr)。把常數整理出來,並注意 r²/(Gr) = r/G:f_G = (4/(a₀³G)) ∫₀^∞ dr r e^{−2r/a₀} sin(Gr)。
  3. 用標準積分(書上沒給,這是本題唯一需要查表或自己推的一步):∫₀^∞ r e^{−pr} sin(qr) dr = 2pq/(p² + q²)²。(推法:寫成 Im ∫₀^∞ r e^{−(p−iq)r} dr = Im [1/(p−iq)²],再取虛部即得。)
  4. 代入 p = 2/a₀、q = G:∫ = 2·(2/a₀)·G / [(4/a₀²) + G²]² = (4G/a₀) · a₀⁴/(4 + G²a₀²)² = 4Ga₀³/(4 + G²a₀²)²。
  5. 乘回前面的係數:f_G = (4/(a₀³G)) · 4Ga₀³/(4 + G²a₀²)² = 16/(4 + G²a₀²)²。得證。
  6. 檢查極限(這才是本題的重點)。G = 0:f = 16/16 = 1 = Z ✓,與 (51) 一致。小 G 展開:f ≈ 1 − G²a₀²/2 + …,對照 (50) 的通用展開 f ≈ Z[1 − G²⟨r²⟩/6] 得 ⟨r²⟩ = 3a₀²,而氫原子基態的 ⟨r²⟩ 確實是 3a₀²——兩條路算出同一個數,交叉驗證通過。
  7. 半高點(代入計算,書中未列):f = ½ 時 (4 + G²a₀²)² = 32 ⇒ G²a₀² = 1.657 ⇒ Ga₀ = 1.29。以 a₀ = 0.529 Å 計,G = 2.43 Å⁻¹,對應 sinθ/λ = G/4π = 0.194 Å⁻¹。用 Cu Kα(λ = 1.5418 Å)就是 2θ ≈ 35°。也就是說,氫的散射能力在中等角度就掉掉一半了——這是 X 光繞射幾乎定不出氫原子位置的根本原因(既只有 1 個電子,又散得特別快),也是有機/生物晶體要用中子繞射定氫位置的理由。
  8. 與真實原子對照:f_G = 16/(4+G²a₀²)² 在大 G 時以 1/G⁴ 衰減。這個 1/G⁴ 尾巴來自 n(r) 在原點的尖點(cusp)——指數函數 e^{−2r/a₀} 在 r = 0 的導數不為零。所有真實原子都有這個尖點(核庫侖奇異性),所以所有原子的 f 在大 G 都是 1/G⁴ 衰減。這是一個相當普遍的結論。

Problem 7:雙原子鏈——把 KCl 的故事壓成一維補推導

最後一題是整段的縮影:它用一維的 ABAB… 鏈,把「幾何決定哪些階次、對比決定它們多亮」這件事講得乾乾淨淨。做完這題,Fig. 17 的 KCl/KBr 對比就不再需要背了。

題目(書本 p.45 逐字):「Consider a line of atoms ABAB . . . AB, with an A—B bond length of ½a. The form factors are f_A, f_B for atoms A, B, respectively. The incident beam of x-rays is perpendicular to the line of atoms. (a) Show that the interference condition is nλ = a cos θ, where θ is the angle between the diffracted beam and the line of atoms. (b) Show that the intensity of the diffracted beam is proportional to |f_A − f_B|² for n odd, and to |f_A + f_B|² for n even. (c) Explain what happens if f_A = f_B.」

  1. 先把幾何弄清楚。A—B 鍵長是 ½a,所以同種原子的週期是 a——這條鏈的晶格常數是 a,基底是 {A 在 0, B 在 ½a}。(PDF 文字層會把 ½ 拆成「1 2」漂到段末,很容易誤讀成鍵長 a,務必小心。)
  2. (a) 入射束垂直於原子鏈,所以入射方向對相鄰原子沒有路徑差(入射波前同時到達)。只有出射方向有:相距 a 的兩顆同種原子,沿出射方向(與鏈夾角 θ)的路徑差是 a cos θ。建設性干涉要求它等於波長的整數倍:nλ = a cos θ。得證。
  3. 順帶一提,這一維的條件就是 Laue 三式 (25) 中的一式:a·Δk = 2πn。入射垂直 ⇒ a·k = 0,所以 a·Δk = a·k′ = ka cos θ = (2π/λ)a cos θ = 2πn,整理即得 nλ = a cos θ。Laue 條件在一維就長這樣。
  4. (b) 用 (46) 的一維版本。分數座標:A 在 x = 0,B 在 x = ½。S(n) = f_A exp(0) + f_B exp[−i2π(n)(½)] = f_A + f_B(−1)^n。
  5. n 為偶數 → (−1)^n = +1 → S = f_A + f_B → I ∝ |f_A + f_B|²。n 為奇數 → (−1)^n = −1 → S = f_A − f_B → I ∝ |f_A − f_B|²。得證。(寫成模平方是因為 f 在吸收邊附近可以是複數;一般情形下 f 為實數,取平方即可。)
  6. (c) 若 f_A = f_B:所有奇數階完全消失,只剩偶數階。物理上,這時兩種原子在 X 光眼中無法區分,這條鏈變成一條週期為 ½a 的單原子鏈——而週期 ½a 的鏈,其干涉條件是 mλ = (½a) cos θ,也就是 2mλ = a cos θ,正好只給出偶數 n = 2m。兩種算法完全一致。
  7. 這就是 Fig. 17 的一維版本:KCl 的 f(K⁺) ≈ f(Cl⁻) 讓全奇反射消失,晶體看起來像常數 a/2 的 sc;這裡 f_A = f_B 讓奇數階消失,鏈看起來像週期 a/2 的單原子鏈。「消光」的通用機制只有一個:當兩套散射體變得無法區分時,真實的週期比你以為的小,於是你以為的一部分反射根本不存在。

校註:本段原文抽取與印刷的六處異常我的判讀

本專案以 PDF 文字層抽取原文,這一段(書本 p.39–45)共發現六處需要說明的地方。前三處是書上本身的排版錯誤(已用 PyMuPDF 放大算繪目視確認),後三處是文字層/解碼器的失真(原書印刷是正確的)。

書上的排版錯誤

位置書上印的應為判定依據
p.39 式 (40)Σ 下標印成 j−1j = 114 倍放大目視;求和從第 1 顆原子跑到第 s 顆
p.42 式 (50)r² 與 sin Gr/Gr 之間印成 =· (相乘)右邊的 r 是左邊積分的啞變數,等號無定義(該積分依 (51) 等於 Z);p.41 同一推導的前一行印的是「·」
p.44 Problem 4分子印成 1 − exp[−iM(a·Δk]…(a·Δk)]少一個右括號;與同式分母的寫法不一致

這三處都不影響物理,但 (50) 那一處如果照抄會讓人完全看不懂,值得特別標出。

文字層的失真(原書無誤)

整理本章時發現了三處抽取失真,其中兩處來自本站解碼器的缺口,當場回頭修好了解碼器:

位置原本抽成實際內容現況
p.40 式 (48)未知字符跨行的大括號 { }解碼器已補上(自訂字符 20902/20903)
p.44 Problem 4(b)ώε(epsilon)解碼器已修正。這個字形全書用了 1074 次(ε_F、ε₀、ε_k 都是它),偏移規則算錯了一格

第三處原先記為「p.44 Problem 3 的點積被解成 ∑」,經回原文查證後確認這條判斷有誤:

p.30 的 Σ n_G exp(iG · r) 抽取結果完全正確,求和符號與點積各自對應不同的字符,

沒有互相誤映。Problem 3 的算式落在未被抽取的顯示區塊裡,無法從文字層判定,

因此不作宣稱。

另外兩類已知陷阱在本段確實出現,已在對應處重組:

(1) 分數被拆散:p.40 的 bcc 位置「x₂ = y₂ = z₂ = ½」抽成「1 2」漂到段末;fcc 的「000; 0½½; ½0½; ½½0」抽成「000; 0 1 2 1 2; 1 20 1 2; 1 2 1 20.」;Problem 7 的鍵長「½a」抽成孤立的「1 2」。以上全部依放大目視還原。

(2) 圖是無文字層的影像或標籤錯位:Fig. 16、17、18 的內容都以 PyMuPDF 算繪後目視判讀,本文中的峰位、峰高、峰寬數值均為目視讀數,不是書上印出的數字,引用時請注意這一點。

本段沒有發現上橫線掉字的問題:p.39–45 出現的指標 (100)、(300)、(111)、(221)、(200)、(110)、(222)、(220)、(311)、(331)、(400)、(420) 全部是正指標,對稱性與數量也都對得上(bcc 奇偶規則、fcc 全奇全偶規則、鑽石 4n 規則逐條驗算皆一致),沒有「同一組指標重複」或「數量對不上」的徵兆。

三個常見誤解,與往後幾章的接點我的判讀

誤解一:「bcc 沒有 (100) 峰,因為 bcc 沒有 {100} 晶面。」

錯。bcc 當然有 {100} 晶面(立方晶胞的六個面就是)。消失的是一階反射,因為體心那層原子把它抵消掉了;二階反射 (200) 好端端地在那裡。正確的說法是「從 {100} 出發的繞射,第一個看得到的是 d = a/2 的那一階」。這件事的根源是書本 p.32 那段關於公因數 n 的警告——繞射指標 (hkl) 與 Miller 指標 (hkl) 是兩套不同的東西,只是寫法一樣。

誤解二:「消光規則是原子排列造成的。」

半對。消光規則是點陣中心化(centering)造成的,不是原子種類造成的。證據:CsCl 的原子排成 bcc 的形狀,但因為角落與體心是不同的離子,(100) 並不消光(只是很弱,S = f_Cs − f_Cl)。同理,閃鋅礦(ZnS)與鑽石的原子排法幾何上相同,但 ZnS 的 (222) 不完全消光。幾何決定「哪些可能為零」,成分決定「是否真的為零」。這就是前面說的硬消光與軟消光之分。

誤解三:「S*S 只是取強度的數學手續。」

不是。(46) 之後書上明白說 S 一般是複數,而實驗只給你 |S|。丟掉的那個相位裡裝的正是原子位置——這就是晶體學一百年來的核心困難。「我們量到的是電子密度傅立葉係數的模長,不是係數本身」這句話應該跟 Δk = G 一樣被記住。


往後幾章會怎麼用這一段

去處用到什麼具體怎麼用
附錄 A(p.641 起)(43) 的 S_G = Σ fⱼ exp(−iG·rⱼ)讓 rⱼ 帶上熱擾動 u(t) 做熱平均,得到 Debye–Waller 因子 I = I₀ exp(−⅓⟨u²⟩G²):升溫時峰變弱但不變寬
第 4–5 章(聲子)BZ 的形狀與 (2π)³/V_c 的體積數聲子模態、算 Debye 溫度與熱容,都要先知道 BZ 裡裝得下多少 k 點
第 7 章(能帶)(26) 的 k·(½G) = (½G)² 與週期位勢的傅立葉展開近自由電子模型的能隙開在 BZ 邊界,位勢的傅立葉係數 U_G 扮演的角色與這裡的 n_G 完全平行
第 17 章(表面與界面)、第 20 章(點缺陷)、第 22 章(合金)消光規則與峰寬 ∝ 1/M超晶格峰與有序/無序轉變(Cu₃Au 型)、晶粒尺寸與缺陷造成的峰展寬

與 Cth(熱容 / 熱阻抗)研究的接點。 對做 FinFET 熱參數萃取的人,這一段有兩條線值得記著:

(1) 附錄 A 的 Debye–Waller 因子 exp(−⅓⟨u²⟩G²) 是一把溫度計。⟨u²⟩ 是原子均方位移,它與晶格振動的振幅直接相關;量不同溫度下 Bragg 峰的強度衰減,就能反推 ⟨u²⟩(T),進而連到 Einstein/Debye 模型的參數——這是一條非電性的、直接量測局部晶格溫度的路徑(實務上稱為 X 光或電子繞射測溫)。對奈米尺度元件而言,電性方法(如閘極二極體溫度計)量的是某個平均,而繞射方法量的是真正的晶格自由度。

(2) BZ 的體積 (2π)³/V_c(Problem 3)是所有熱容計算的入場券。聲子熱容、Debye 溫度、熱導率的 kinetic 公式,全部從「BZ 裡有 N 個 k 點」這個計數開始。第 2 章這幾頁看起來只在談繞射,但它建立的倒空間框架,才是第 5 章能把 C = ∂U/∂T 算出來的前提。