Kittel·固態物理

全書逐章 / 第 2 章

02

波繞射與倒晶格

Wave Diffraction and the Reciprocal Lattice

書本 p.23–455 個小節55 條公式36 段推導約 68,541 字熱容關聯 間接

①主要內容與研究問題

第 1 章教你怎麼描述一個晶體,第 2 章回答一個更難堪的問題:你怎麼知道它真的長那樣?

沒有顯微鏡看得到原子的週期排列——可見光的 5000 Å 對 Å 級的晶格而言太粗,書本 p.25 開宗明義說這種波長下各原子的彈性散射疊加起來「只是普通的光學折射」。唯一的門路是繞射。所以第 2 章表面上在講 X 光、中子、電子怎麼被晶體散射,實際上是在建立整本書唯一的結構證據鏈:後面二十章所有「矽是 diamond 結構、a = 5.430 Å」這類陳述,全部是這一章的方法量出來的。

這一章的敘事結構值得先看懂,因為它是全書少見的「先給一個會用但不成立的答案,再拆掉重蓋」。第一步是 Bragg 定律 2d sin θ = nλ,三行幾何就推完,而且 Kittel 自己在 p.25 就承認這個推導 is convincing only because it reproduces the correct result——它假設了原子面像鏡子一樣鏡面反射,這個假設沒有任何正當性。更要命的是它只算得出峰在哪個角度,算不出峰有多亮,而強度才是真正裝著原子位置的地方。

第二步是重蓋。把電子密度 n(r) 攤成傅立葉級數,逼問「週期性允許哪些傅立葉分量存在」,答案唯一:波向量必須落在一組由 bᵢ·aⱼ = 2πδᵢⱼ 定義的點陣上。這就是倒晶格。注意這個邏輯方向——倒晶格不是誰為了方便發明的工具,它是 n(r + T) = n(r) 這一條件的唯一解。接著散射振幅被算出來就是 n(r) 的傅立葉變換,於是 Δk = G 這個繞射條件掉出來,Bragg 定律反過來成為它的推論。這一步的哲學重量在於:繞射圖是倒晶格的地圖,不是實空間的照片。你量到的是 |n_G|,不是 n_G——相位在取強度時被丟掉了,這正是晶體學一百年來的核心困難。

第三步是幾何化。把 2k·G = G² 兩邊除以 4 得到 k·(½G) = (½G)²,它說 k 的端點落在某條倒晶格向量的垂直平分面上。這些平面把整個 Fourier 空間切碎,最靠近原點的那一塊就是第一布里淵區。書中 p.34 很誠實地指出:X 光繞射其實不需要布里淵區這個概念,是能帶論不能沒有它。

第四步是把基底放回去,得到結構因子 S_G 與原子形狀因子 f,回答「哪些幾何上允許的峰其實是零」——也就是消光。到這裡,「峰在哪裡」歸晶格管、「峰多亮」歸基底管,第 1 章那個 lattice + basis 的拆分在這一章收割。

這一章因此撐起後面三條主線:Ch4–5 的聲子色散與模數計數活在第一布里淵區裡;Ch7 的能隙就是把 2k·G = G² 這同一條式子用在電子波上的結果(能隙不是一個獨立假設);Ch9 的費米面畫在同一個倒空間。連 Ch4 p.100 的「聲子動量」與 p.101 的選擇律 k + G = k′ ± K,都只是本章 Δk = G 放寬彈性條件之後的推廣。一句話:第 2 章不是一章實驗方法,它是全書的座標系。

②③核心概念與公式 · 逐節深講

依書中實際小節切分,一節一頁。未標層級者為書中內容,頁碼與式號逐條核對;書中省略而由本站補上的代數步驟標為「補推導」。

01晶體對波的繞射——Bragg 定律DIFFRACTION OF WAVES BY CRYSTALS — The Bragg law(含 SCATTERED WAVE AMPLITUDE 標題頁)書本 p.23-26 · 概念 17 推導 5 公式 5用一條三行就能推完的式子 2d sin θ = nλ 打開晶體結構的大門,同時誠實交代這條式子為什麼不夠用——它只告訴你峰在哪裡,不告訴你峰有多強,而這正是後面整章要改用 Fourier 分析與倒晶格重做一次的理由。02散射波振幅——傅立葉分析與倒晶格向量SCATTERED WAVE AMPLITUDE — Fourier Analysis; Reciprocal Lattice Vectors書本 p.27-29 · 概念 13 推導 5 公式 17Bragg 定律只算得出繞射「往哪個角度去」,算不出「有多亮」。這一段把晶體的電子密度攤成傅立葉級數,逼問「哪些傅立葉分量被晶格週期性允許」——答案就是倒晶格;而 bᵢ·aⱼ = 2πδᵢⱼ 這條定義,存在的唯一理由是要讓 exp(iG·T) 剛好等於 1。03散射波振幅——繞射條件與 Laue 方程SCATTERED WAVE AMPLITUDE — Diffraction Conditions; Laue Equations書本 p.30-32 · 概念 12 推導 7 公式 10把 Bragg 的「鏡面反射」圖像換成真正的散射理論:散射振幅 F 其實就是電子密度 n(r) 的傅立葉變換,於是繞射峰只在 Δk = G 時出現、而峰的振幅正好等於 n(r) 的某一個傅立葉分量——繞射圖因此是倒晶格的地圖,不是實空間的照片。04布里淵區——sc/bcc/fcc 的倒晶格BRILLOUIN ZONES — Reciprocal Lattice to sc / bcc / fcc Lattice書本 p.33-38 · 概念 12 推導 9 公式 14把繞射條件 2k·G = G² 除以 4,它就從一條代數等式變成一句幾何陳述——k 的端點落在某條倒晶格向量的垂直平分面上;於是整個 Fourier 空間被這些平面切成布里淵區,而第一布里淵區就是倒晶格的 Wigner–Seitz 原胞。這一段接著把這個作圖法實際套到 sc、bcc、fcc 三種立方晶格上,算出立方體、菱形十二面體與截角八面體,並建立起「正空間原胞越小、倒空間布里淵區越大」這條後面每一章都要用的換算。05基底的傅立葉分析——結構因子、原子形狀因子;本章摘要與習題FOURIER ANALYSIS OF THE BASIS — Structure Factors, Atomic Form Factor; SUMMARY; Problems書本 p.39-45 · 概念 15 推導 10 公式 9前面幾節只用了晶格、只回答了「繞射峰出現在哪裡」;這一段把基底放回去,回答「峰有多亮、哪些峰根本不出現」——結構因子管原子之間的干涉,原子形狀因子管原子內部的干涉。

⊘常見誤解與易錯點 (我的判讀)

  • 「第一布里淵區是 k 的允許範圍,k 超出 BZ 就沒有物理意義」
    把「冗餘」誤讀成「禁止」。BZ 外的 k 完全合法,只是對定義在格點上的量而言不帶新資訊:一維鏈的位移 u_s = u·exp(iska) 只在整數 s 上有定義,k 與 k + 2π/a 給出逐點完全相同的位移,所以折回 BZ 不損失任何東西。但對定義在連續空間的量(電子波函數 ψ(r)、X 光的波向量 k)就完全不是這樣——k 與 k + G 是兩個真正不同的波,只是能量本徵值可以用能帶指標重新編號而已。
    判準只有一句:問這個波的振幅定義在哪裡。只在離散格點上取值 → 折回第一 BZ(聲子、緊束縛振幅);在連續 r 上取值 → 延伸區圖像才是物理的(光子、電子的真實動量)。立刻用一個數字檢查:Cu Kα(λ = 1.5418 Å)的 k = 2π/λ = 4.075 Å⁻¹,而矽第一布里淵區最遠的 X 點只在 2π/a = 1.157 Å⁻¹(a = 5.430 Å 出自書本 p.20 Table 3,兩個 k 值為本站計算)。Ewald 球的半徑是第一 BZ 尺度的 3.5 倍——如果 k 真的不能超出 BZ,X 光繞射根本做不成。這也是為什麼繞射畫在延伸的 Fourier 空間裡(書本 p.30:wavevectors are always drawn in Fourier space),而聲子色散卻只畫在 BZ 內。兩者不是矛盾,是兩種不同的量。
  • 「Δk = G 就是動量守恆,ℏG 是光子交給晶格的動量,所以繞射會讓晶格升溫」
    前半句對、後半句錯,錯在沒追問「交給誰」。Kittel 在 Ch4 p.100 把這件事講得很白:In the reflection process the crystal as a whole will recoil with momentum −ℏG,並強調 The true momentum of the whole system always is rigorously conserved。反衝的是整塊晶體(約 10²³ 個原子),不是單一原子。真正的動量嚴格守恆;被放寬成「模 G 守恆」的是晶體動量,因為晶格只有離散平移對稱而非連續平移對稱。
    用能量帳自己檢驗一次就記得住:若把 ℏG 的反衝丟給一個矽原子,反衝能 (ℏG)²/2M 對 (400) 反射(G = 8π/a = 4.628 Å⁻¹,M = 28.09 u)是 1.59 meV(本站計算),相對於 12.4 keV 的光子能量是 1.3 × 10⁻⁷。但「彈性散射」在實驗上準確得遠比 10⁻⁷ 還好——這只有在反衝質量是宏觀的時候才說得通。這與 Mössbauer 效應是同一個機制(書本 p.643 明確把 Debye–Waller 的理論同時適用到 the Mössbauer effect)。操作判準:彈性 Bragg 繞射不加熱晶體;會加熱的是非彈性的那一份,而它跑到瀰散背景裡去了(見下一條)。
  • 「Bragg 峰很窄代表晶體完美;樣品升溫或受損時峰會變寬」
    把峰寬和峰強兩個機制混成同一個。溫度進入的是 Debye–Waller 因子,書本附錄 A p.643 式 (6) 為 I(hkl) = I₀·exp(−k_B T G²/Mω²)——它只乘在強度上,完全不碰峰的角位置與寬度。原文也說得很清楚:升溫時 the intensity of the diffracted line decreases, but not catastrophically,而失去的強度 appears as a diffuse background。峰變寬的原因是另一組完全不同的東西:有限晶粒尺寸(本章 Problem 4 的峰寬 ∝ 1/M)、應變分布、儀器解析度、mosaic spread。
    三種效應各有指紋,記住就能當場分辨:(a) 溫度 → 所有峰強度下降,且高 |G| 的峰掉得比低 |G| 快(p.643:Reflections of low G are affected less),寬度不變;(b) 晶粒尺寸 → 展寬 β ∝ 1/cos θ;(c) 應變 → 展寬 β ∝ tan θ。因此把 β·cos θ 對 sin θ 作圖,截距給尺寸、斜率給應變。【延伸知識】這張圖叫 Williamson–Hall plot,尺寸那條關係叫 Scherrer 公式——這兩個名詞在 Kittel 第 8 版中零次出現(已全書逐頁搜尋確認),書上只有 Problem 4 的 1/M 定性結果。直接的實務後果:你量 fin 或應變層的 HRXRD 時看到峰變弱但沒變寬,那是溫度,不是缺陷,不要當成樣品壞了。
  • 「繞射斑點之間的瀰散背景是儀器雜訊,扣掉就好」
    扣掉的其實是聲子。本章「強度只集中在倒晶格點上」這個結論,前提是無限大、完美、靜止的晶格;只要原子在動,強度就會漏到 G 點以外。書本 p.643 給的數字很有說服力:即使在絕對零度,光是零點運動就已經讓 90 percent of the beam is elastically scattered and 10 percent is inelastically scattered。那 10% 不是雜訊,它是晶格振動的訊號。
    把本章的 Δk = G 看成彈性極限,真正的一般式在 Ch4 p.101 式 (33):k + G = k′ ± K,K 是被產生或吸收的聲子波向量。於是「背景」其實是一張 ω(K) 的地圖——非彈性中子散射(書本 p.101 起)與非彈性 X 光散射量的就是它。操作判準:你的分析若只用峰位與峰強(結構精修),扣背景是對的;你若想知道晶格動力學——而 Cth 正是晶格動力學對頻率的積分——背景才是你要的東西。Debye 模型用一條直線色散去近似的那個 ω(K),實驗上就是從這裡量出來的。
  • 「TEM 選區繞射(SAED)的斑點亮度就是 |S_G|²,所以可以從斑點強度反推原子位置」
    本章整套理論(散射振幅 = V·n_G、I ∝ |S_G|²)建立在一個被 p.25 明講出來的前提上:每個原子面只反射 10⁻³ 到 10⁻⁵ 的入射輻射,所以一個光子從進到出頂多散射一次。電子與物質的作用強了好幾個數量級(本章 Fig. 1 顯示電子只要約 150 eV 就有 1 Å 波長,代價是極強的交互作用、穿透深度只有幾個原子層),單次散射的前提直接破功。
    斑點的位置仍然可信,強度不可信。位置由 Δk = G 決定,那是純幾何,本章保證了它;強度則被多重散射重新分配。最好用的當場診斷是看「不該出現的峰」有沒有出現:矽是 diamond 結構,(200) 滿足 fcc 的全偶條件卻不滿足 4n 規則(本章 Problem 5),結構因子嚴格為零——粉末 XRD(書本 Fig. 4)裡它確實不見;但在 SAED 裡 (200) 幾乎一定會亮。那不是樣品有異常,那是電子散射了兩次(先 (111) 再 (1̄1̄1))。【延伸知識】「kinematic / dynamical theory」「multiple scattering」「double diffraction」這組詞在 Kittel 第 8 版中零次出現(已全書逐頁搜尋確認;全書唯一的 kinematic 出現在 p.101,指的是中子散射的運動學關係,不是這個意思)。要量化電子繞射強度必須跑 Bloch wave 或 multislice 模擬,不能套本章的公式。

④本章最重要的結論

繞射的完整條件是 Δk = G:散射向量必須等於某個倒晶格向量,因此繞射圖樣就是晶體倒晶格的一張地圖。這個條件有三種等價敘述——Bragg 的 2d sinθ = nλ、Laue 方程 aᵢ·Δk = 2πvᵢ、以及 Brillouin zone 的 k·(½G) = (½G)²,後者把「可被繞射的 k」幾何化為布里淵區邊界。繞射峰的位置由晶格決定,強度則由基決定:structure factor S_G 與 atomic form factor fⱼ 一起解釋了 bcc(v₁+v₂+v₃ 為奇數)與 fcc(指標奇偶混雜)的系統性消光。第一布里淵區作為倒晶格的 Wigner-Seitz cell,後續成為能帶理論與各種元激發分析的基本舞台。

書中逐字

"A diffraction pattern of a crystal is, as we shall show, a map of the reciprocal lattice of the crystal." (p.30)

"This is the central result of the theory of elastic scattering of waves in a periodic lattice." (p.32)

"A Brillouin zone is defined as a Wigner-Seitz primitive cell in the reciprocal lattice." (p.33)

"Only waves whose wavevector k drawn from the origin terminates on a surface of the Brillouin zone can be diffracted by the crystal." (p.43)

⑤與前後章的關係 (我的判讀)

直接建在 Ch1 的 a₁,a₂,a₃ 上。它產出的第一 Brillouin 區成為後續三條主線的共同定義域:Ch4–5 的聲子色散與模式計數、Ch7 的能帶、Ch9 的費米面。Ch12 的磁振子色散也畫在同一個倒空間裡。

⑥與熱容 Cth 的關聯程度 (我的判讀)

間接

Ch5 的 Debye 截止波向量 K_D、以及「模式總數 = 3N」這個論證,全部活在第一 Brillouin 區這個空間裡。不需精讀,知道 Brillouin 區是什麼即可。

結論先講:間接——本章不提供任何熱物理,它提供的是你之後算 Cth 時所在的那個空間

先把證據放上來,免得聽起來像客套話。對本章書本 pp.23–45(PDF pp.43–65)逐頁搜尋,「Debye」「heat capacity」「temperature」「phonon」四個詞出現次數皆為 0。這一章從頭到尾沒有出現過溫度。所以任何「第 2 章教你算熱容」的說法都是硬扯。

它真正的貢獻只有一條,但這一條繞不過去:V_BZ = (2π)³/V_c(本章 Problem 3)。所有熱容計算都從「第一布里淵區裡裝得下幾個 k 點」開始。


一、矽的倒空間數字:算一次,之後每一章都在用

量來源數值
晶格常數 a書本 p.20 Table 3(該頁無文字層,目視判讀)5.430 Å
原胞體積 V_c = ¼a³(diamond = fcc + 2 原子基底)本站計算40.03 ų
布里淵區體積 (2π)³/V_c本站計算6.197 Å⁻³
X 點(2π/a)本站計算1.157 Å⁻¹
L 點(√3π/a)本站計算1.002 Å⁻¹
Debye 球半徑(每原胞 1 個 k:(6π²/V_c)^⅓)本站計算1.139 Å⁻¹
Debye 球半徑(每原子 1 個 k:(6π²n_atom)^⅓,n_atom = 8/a³)本站計算1.436 Å⁻¹
Debye 溫度 θ(低溫極限)書本 p.116 Table 1(影像頁,200 dpi 算繪後目視判讀)645 K
熱導率 κ(300 K)書本 p.116 Table 1(同上)1.48 W·cm⁻¹·K⁻¹

倒數第三列的 1.436 Å⁻¹ 是本章唯一一個會直接咬到你 Cth 擬合的數字。它超出了矽第一布里淵區最遠的邊界(X 點 1.157 Å⁻¹)。這不是算錯:標準 Debye 模型要把 3pN = 6N 個模全塞進一個球,而第一布里淵區只裝得下 3N 個(3 支聲學支 × N 個 k 點)。多出來的 3N 個光學模被硬塞進一層超出 BZ 的假想球殼。

可操作的後果:你如果用 Debye 形式的 C(T) 去擬合量到的 Cth(T),擬出來的 θ 會隨溫度區間漂移——低溫段擬出的 θ 偏低(只有聲學支被激發,模型還算對),高溫段偏高(光學支開始貢獻,但 Debye 球把它們放在錯的頻率上)。這是模型的系統偏差,不是你的量測誤差。要避開就別用單一 θ,改用 Debye + Einstein 兩項(光學支近似為平頻)。


二、本章的實驗手段能不能拿來量 FinFET 的溫度?算過之後:不能

這是最容易被誤導的地方。繞射確實內建一支溫度計——附錄 A p.643 的 Debye–Waller 因子 I(hkl) = I₀·exp(−k_B T G²/Mω²)(書本式 (6),古典高溫近似)。量峰強隨溫度的衰減就能反推晶格溫度,而且量到的是真正的晶格自由度,不是電性方法給你的某個等效平均。聽起來很適合 self-heating 研究。把數字代進去就知道為什麼不行。

以矽計算,取 ω ≈ ω_D = k_Bθ/ℏ = 8.44 × 10¹³ s⁻¹(由書本 θ = 645 K 換算),M = 28.09 u:

反射G (Å⁻¹)300 K 的 I/I₀靈敏度 d(ln I)/dT
(111)2.0040.9950.0017 %/K
(400)4.6280.9740.0089 %/K
(800)9.2570.8990.036 %/K

(G、I/I₀ 與靈敏度皆為本站以書本式 (6) 計算。注意 300 K < θ = 645 K,古典式低估 ⟨u²⟩,真值會再小一點。)

要分辨 1 K,(400) 峰的強度變化只有 8.9 × 10⁻⁵。光是光子統計就要約 10⁸ 個計數,系統穩定度還必須優於 10⁻⁵——這在同步輻射上勉強做得到,在一般實驗室做不到。更致命的是時間尺度:FinFET 的熱時間常數在 ns 量級,而一般繞射的積分時間是秒,差了九個數量級。【延伸知識】要在 ns 解析熱響應必須做 pump–probe 時間解析繞射(XFEL 或同步輻射 slicing);「time-resolved diffraction」「pump-probe」這組概念在 Kittel 第 8 版中零次出現(已全書逐頁搜尋確認)。

所以:本章的方法不是你萃取 Cth 的路。你的路仍然是 Zth(f) 或脈寬掃描這類電熱量測。這一節值得算過一次,是為了讓你以後看到「用 XRD 量元件溫度」的提案時,能立刻說出為什麼不行。


三、本章真正會碰到你實驗的地方:量出 Cth 的那個分母

Cth = c_v × V。c_v 是材料性質(Ch5 的事),V 是幾何,而奈米級幾何在產線與實驗室都是靠繞射量的:HRXRD 的倒空間圖(reciprocal space map,就是本章 Δk = G 的直接視覺化)量 fin 與 SiGe 應力層的應變與週期,小角散射量 fin 的臨界尺寸。本章 Problem 4 的「峰寬 ∝ 1/M」正是這件事的理論根。【延伸知識】具體技術名稱(CD-SAXS、reciprocal space mapping)不在書中。

順帶把量級感補完,因為它會影響你怎麼解讀量到的 Cth。用書本 p.116 的 κ = 1.48 W·cm⁻¹·K⁻¹ 與矽 300 K 的體積熱容 ≈ 1.64 J·cm⁻³·K⁻¹(外部值,非書中),熱擴散率 α = κ/c_v = 0.90 cm²/s(本站計算)。1 ns 內的熱擴散距離 √(αt) ≈ 300 nm,遠大於 fin 本身。換算成被加熱的體積 (300 nm)³ ≈ 2.7 × 10⁻¹⁴ cm³,對應 Cth ≈ 4 × 10⁻¹⁴ J/K,比 fin 本體的約 10⁻¹⁷ J/K(見第 1 章)大了三個數量級以上。

這代表你從 Zth(f) 拆出來的「Cth」不是一個材料常數,而是一個與量測頻率相依的有效值:頻率越低、熱擴散得越遠、看到的有效體積越大。報告 Cth 時務必連同對應的時間尺度一起報。(這個估計用的是塊材 κ;fin 尺度下邊界散射會把 κ 壓低數倍——書本 Ch5 p.126–127 討論的正是這件事——所以 300 nm 是上限。)


四、讀本章時力氣該怎麼分配

小節對 Cth 研究建議
1. Bragg 定律(p.23–26)無略讀。知道 λ ≤ 2d、以及 θ 從晶面量起即可
2. 傅立葉分析與倒晶格(p.27–29)間接(它定義 b₁b₂b₃)要讀,但只需讀到會由 a 算出 b、會算 V_c
3. 繞射條件與 Laue 方程(p.30–32)無Laue 三式與 Ewald 作圖可以完全跳過
4. 布里淵區(p.33–38)這一節是重點精讀。V_BZ = (2π)³/V_c 與 fcc 的截角八面體要真的懂
5. 結構因子與形狀因子(p.39–45)無(除非你要親自做 XRD 量測)略讀;Problem 4(峰寬 ∝ 1/M)值得看

一句話:第 2 章對 Cth 研究而言讀第 4 小節就夠了,其餘是為做繞射實驗的人準備的。