Kittel·固態物理

全書逐章 / 第 6 章 自由電子費米氣 / 小節 2

02三維自由電子氣

FREE ELECTRON GAS IN THREE DIMENSIONS

書本 p.137-141 · 小節 2 / 6 · 概念 20 推導 9 公式 19

把一維的「依 Pauli 原理逐層填能階」搬到三維:先對照盒中駐波與週期性邊界條件下的行進平面波,確立 k 空間每 (2π/L)³ 恰有一個允許的 k,再把基態看成填滿半徑 k_F 的 Fermi 球,一路推出只依賴電子濃度的 k_F、ε_F、v_F、T_F 與態密度 D(ε) = 3N/2ε,並以 Table 1 列出的電子濃度驗算這些量——ε_F 是幾個 eV、T_F 是幾萬 K,這個懸殊正是下一段電子熱容之所以只有古典值 1% 的根源。

這一段在解決什麼問題:把一維的「填能階」搬到三維,並換成會走動的波我的判讀

p.134–135 的一維問題已經把整套邏輯的骨架立好:一個電子被無限高的位壘關在長度 L 的線段上,波函數是整數個半波長的駐波,能量 ε_n = (ħ²/2m)(nπ/L)²;接著依 Pauli 不相容原理,每個 n 放自旋向上、向下各一個電子,由 2n_F = N 決定最高填滿的能階,得到一維的 Fermi 能 ε_F = (ħ²/2m)(Nπ/2L)²(式 (4))。p.136–137 再用 Fermi-Dirac 分布 (5) 說明溫度升高時,佔據機率只在 ε_F 附近被抹開。

一維模型有兩個不足,逼得我們必須走到三維。第一,態在能量上的分布方式完全不同:後面會算出一維的態密度隨能量遞減(∝ ε^{−1/2}),三維則隨 √ε 遞增;這個斜率的正負,直接決定化學勢隨溫度往哪個方向偏(p.144 Fig. 8:一維往上、三維往下)。第二,一維用的是固定端的駐波,駐波不帶淨動量,無法描述電流;而本章後半(p.147 起)要談導電,需要的是能攜帶動量的行進波。

本段(p.137–141)的推導鏈可以壓成一行:

薛丁格方程 (6) → 邊界條件(駐波 (7) 或週期性 (8)–(11))→ 允許的 k 在 k 空間排成間距 2π/L 的立方格子 → 基態填滿一顆 Fermi 球 (14)–(15) → k_F (16) → ε_F (17) → v_F (18) → 能量 ≤ ε 的軌域數 (19) → 態密度 D(ε) (20) → D = 3N/2ε (21)。

整條鏈只用到兩個輸入:自由粒子的色散 ε = ħ²k²/2m,以及「k 空間每 (2π/L)³ 的體積恰有一個允許點」這個純計數事實。推到最後,唯一留下來的材料參數是電子濃度 n = N/V。這就是為什麼 p.139 的 Table 1 只要給一欄 n,就能一口氣算出 k_F、v_F、ε_F、T_F 四欄——本講解會實際代數字驗算給你看。

讀這一段時,請把它和第 5 章 p.111 的聲子三維態密度並排看:同樣的週期性邊界條件、同樣的 (2π/L)³、同樣的「球體積除以單格體積」。差別只有三處:(i) 色散從聲子的 ω = vK 換成電子的 ε = ħ²k²/2m;(ii) 每個 K 有 3 個聲子偏振,每個 k 有 2 個電子自旋;(iii) 聲子是玻色子、數目不守恆,Debye 模型要用截止頻率去「湊」出正確的模態總數;電子是費米子、數目固定為 N,Fermi 球的大小由 N 自然決定,不需要任何人為截止。

最後,p.140 的 Fig. 5 把 D(ε) 與有限溫度下的佔據 f·D 畫在一起,並在圖說寫明:升溫時電子被熱激發「由區域 1 到區域 2」。這張圖是下一段(p.141 起)電子熱容論證的起點;本段算出的 D(ε_F) 與 T_F 會在那裡被直接代用。

三維自由粒子薛丁格方程

−(ħ²/2m)(∂²/∂x² + ∂²/∂y² + ∂²/∂z²) ψ_k(r) = ε_k ψ_k(r)

(6)  p.137

ħ:約化 Planck 常數;m:電子質量(本段一律用自由電子質量,有效質量要到 p.145 才以「熱有效質量」的形式出現);ψ_k(r):以波向量 k 標記的單電子波函數,Kittel 稱之為軌域(orbital,p.134 的定義:只含一個電子的系統之波動方程的解);ε_k:該軌域的能量;括號內就是 Laplacian ∇²。

出處推導
p.134 一維式 (1) 的三維推廣:Hamiltonian 只保留動能 p²/2m,再以動量算符 p → −iħ∇ 代入。背後是自由電子模型的兩個假設:(a) 不計傳導電子與離子實的交互作用——p.133 說這個交互作用留到下一章處理;(b) 不計電子彼此的交互作用——p.134 指出,把 N 個電子分派到 N 個不同軌域的「軌域模型」,只有在電子之間沒有交互作用時才是精確的。
極限行為
(i) 若金屬內部有一個處處相同的常數位能 V₀(例如把位能井底當能量原點),它只會讓所有 ε_k 同時平移 V₀,波函數與能量差都不變;所以「忽略位能」真正的意思是「忽略位能在空間中的起伏」。(ii) 一旦加入具晶格週期的位能,平面波不再是本徵函數,本徵函數變成「平面波 × 晶格週期函數」的 Bloch 形式(第 7 章 p.167 式 (7)),並在布里淵區邊界打開能隙;第 7 章 p.164 的近自由電子模型就是從本段的 (12)、(10)、(9) 出發。(iii) 自由電子近似對鹼金屬特別好:p.133 指出鈉的離子實只佔晶體體積約 15%;p.155 也指出 Na、K 的精確 Hall 係數與「每原子一個傳導電子」的自由電子計算值吻合得極好。

立方盒中的駐波解

ψ_n(r) = A sin(πn_x x/L) sin(πn_y y/L) sin(πn_z z/L)

(7)  p.137

L:立方盒的邊長,原點放在立方體的一個角上(書中明寫);n_x, n_y, n_z:正整數(書中明寫);A:常數。書中沒有給 A 的值;若要求 ∫|ψ|²d³r = 1,由 ∫₀ᴸ sin²(πnx/L)dx = L/2 可得 A = (2/L)^{3/2}(此歸一化為本講解補充)。

出處推導
電子被關在邊長 L 的立方體中、壁外位能無限大,波函數必須在六個面上為零——這是 p.134 一維邊界條件 ψ(0) = ψ(L) = 0 的三維版本。以分離變數 ψ = X(x)Y(y)Z(z) 代入 (6),每個方向都是一個一維無限位井,解就是 p.134 式 (2) 的正弦函數;三個方向相乘即得 (7)。代回 (6) 得 ε = (ħ²/2m)(π/L)²(n_x² + n_y² + n_z²)(書中此處沒有寫出,它就是 p.134 式 (3) 在三個方向的加總)。
極限行為
(i) 對應的波數是 k_i = πn_i/L,只落在 k 空間的第一卦限、間距 π/L——和下面週期性條件的 2π/L 不同,但數出來的態數相同(見後面的推導區塊)。(ii) 奈米尺度:當 L 只有幾個奈米,最低幾組 (n_x, n_y, n_z) 的能量間距可以大過 k_BT,這時駐波才是正確的物理描述,週期性邊界條件不再適用;第 18 章 p.528 談矩形奈米線時,明白說截面方向的能量與本徵態就是第 6 章討論過的盒中粒子(particle-in-a-box)解。(iii) 巨觀樣品:L = 1 cm 時最低能階只有 ~10⁻¹⁴ eV(見數值欄),能階間距遠小於任何熱能或 ε_F,能譜實際上是連續的——所以後面可以放心把「數點」換成「算體積」。
典型數值
L = 1 cm:ħ²/2m = 3.81 eV·Å² = 3.81×10⁻¹⁶ eV·cm²,(π/L)² = 9.87 cm⁻²,基態 (1,1,1) 的能量 = 3 × 3.81×10⁻¹⁶ × 9.87 ≈ 1.1×10⁻¹⁴ eV;對照 Cu 的 ε_F = 7.00 eV(Table 1),小了約 15 個數量級。

駐波的三個性質:只取正整數、三維才有的簡併、以及它不帶動量我的判讀

(7) 雖然只有一行,卻藏了三件後面都會用到的事。

第一,為什麼 n_x、n_y、n_z 只能是正整數?逐一檢查就知道:n = 0 時 sin(0) ≡ 0,整個波函數消失,那不是一個態;n 取負數時 sin(−πnx/L) = −sin(πnx/L),只是原波函數乘上 −1,而整體相位不可觀測,所以是同一個態。因此每個方向的獨立量子數只有 1, 2, 3, …,對應的波數 k_i = πn_i/L 只落在 k 空間的第一卦限,相鄰間距 π/L。

第二,三維才真正有簡併。把 (7) 代回 (6),每個方向各貢獻 (ħ²/2m)(πn_i/L)²,總能量 ε = (ħ²/2m)(π/L)²(n_x² + n_y² + n_z²)。一維時一個 n 只對應一個空間態(再乘上自旋 2 重);三維時 (1,1,2)、(1,2,1)、(2,1,1) 三組量子數給出同一個能量,是三重的空間簡併,再乘上自旋就是 6 個軌域。p.135 把「同一能量的軌域數」定義為簡併度——這個定義到了三維才開始有實質內容,而「每單位能量有多少軌域」正是本段最後要算的態密度。

第三,駐波不是動量本徵態。一維的 sin(kx) = [e^{ikx} − e^{−ikx}]/2i,是動量 +ħk 與 −ħk 兩個行進波等量的疊加;對它量動量,會以各一半的機率得到 ±ħk,期望值為零,也不帶淨機率流。對「數有幾個態」這個問題,駐波完全夠用;但只要問題牽涉電流——例如 p.147 起讓電場推著整顆 Fermi 球走——就需要本身帶有確定動量的行進波。這是接下來改用週期性邊界條件的物理動機。(p.137 原文只說這樣做比較方便,並指出第 5 章處理聲子時也這麼做;「為了描述電流需要行進波」是本講解補上的理由。)

改用週期性邊界條件(p.137–138 的論述)

給出駐波 (7) 之後,Kittel 立刻換一條路:他說引入滿足週期性邊界條件的波函數比較方便,並點明這和第 5 章處理聲子時的做法相同(p.137)。具體要求是:波函數在 x、y、z 三個方向都以 L 為週期,即式 (8);y、z 兩個方向同理。

在這個條件下,同時滿足自由粒子薛丁格方程 (6) 與週期性的解,是行進平面波 (9),前提是波向量 k 的分量取 (10) 所列的值(k_y、k_z 同理)。原文接著說明:凡是形如 2nπ/L 的分量,都滿足長度 L 上的週期性,n 為正或負整數;並指出 k 的三個分量就是這個問題的量子數,再加上描述自旋方向的量子數 m_s。最後以 (11) 直接代入,驗算 k_x = 2nπ/L 確實滿足 (8)。

p.138 接著把 (9) 代回 (6) 得到軌域能量 (12),並指出波向量大小與波長的關係 k = 2π/λ;再用動量算符 p = −iħ∇ 說明平面波是線動量的本徵函數、本徵值為 ħk(式 (13)),軌域 k 中粒子的速度為 v = ħk/m。

週期性邊界條件

ψ(x + L, y, z) = ψ(x, y, z)

(8)  p.137

L:週期長度,取為樣品邊長。原文註明 y、z 座標同理,所以實際上是三條獨立的條件。

出處推導
不是來自樣品真實的表面,而是一個計數上的約定:把樣品想成三個方向都首尾相接(一維是環,三維是三維環面),因此沒有牆壁、沒有反射。p.137 說明這與第 5 章處理聲子的做法相同;第 5 章 p.110 對一維原子鏈用過完全相同的手法(Fig. 4 把鏈的首尾接成環),並說明對大系統而言,這個方法不會在任何本質上改變問題的物理。
極限行為
(i) 巨觀極限:體性質(每單位體積、每單位能量的軌域數)與邊界條件的選擇無關,週期性條件與盒中駐波的差別只是表面項,相對大小約 1/(k_F L)。(ii) 表面物理與奈米結構:真實的邊界不能再忽略,(8) 失效——第 18 章處理奈米結構時回到盒中駐波與離散能階(p.520、p.528)。(iii) 第 7 章 p.164 的近自由電子模型開頭,就直接沿用「立方體邊長 L 上的週期性邊界條件」與 (10) 的允許 k 值(該處式 (2))。

行進平面波

ψ_k(r) = exp(ik·r)

(9)  p.137

k = (k_x, k_y, k_z):波向量;r:位置向量。書中以方框標示此式,且沒有附歸一化常數。

出處推導
把 exp(ik·r) = exp(ik_x x)·exp(ik_y y)·exp(ik_z z) 代入 (6):每個方向的二階導數拉下 −k_i²,得 −(ħ²/2m)(−k²)ψ = (ħ²k²/2m)ψ,所以它對任何 k 都是 (6) 的解;再要求它滿足 (8),k 就被限制在 (10) 的離散值上。
極限行為
(i) 它是動量本徵態((13)),也是行進波:等相位面 k·r = 常數 沿 k 的方向推進。(ii) |ψ|² = 1 處處相同——電子均勻分布在整個樣品中,與「電子氣」的均勻電荷圖像一致;若要歸一化到體積 V,取 V^{−1/2}exp(ik·r)(本講解補充)。(iii) 加入晶格週期位能後,第 7 章 p.167 的 Bloch 定理說本徵函數是「平面波 exp(ik·r) × 具晶格週期的函數 u_k(r)」;p.173 並指出,晶格位能消失時 u_k(r) 變成常數,Bloch 函數就回到自由電子的 exp(ik·r)——(9) 正是這個極限。

週期性條件下允許的波向量分量

k_x = 0 ; ±2π/L ; ±4π/L ; …

(10)  p.137

k_x:波向量的 x 分量;原文註明 k_y、k_z 同理。原文並寫成一般形式:任何形如 2nπ/L 的分量都滿足週期性,n 為正或負整數((10) 本身明列 0)。

出處推導
(9) 代入 (8) 得 exp(ik_x L) = 1,因此 k_x L = 2πn(逐步推導見 (11) 之後的推導區塊)。
極限行為
(i) 相鄰允許值的間距恆為 2π/L,與 k 的大小無關:允許的 k 在 k 空間排成晶格常數 2π/L 的簡單立方格子,每點佔 (2π/L)³ = 8π³/V 的體積——這是 (15) 的全部依據。(ii) 自由電子的 k 沒有上限:這裡沒有晶格,k 可以任意大;對照第 5 章 p.111 式 (17),聲子的 K 被截在 Nπ/L,因為原子是離散的。電子要到第 7 章加入晶格後,才會把 k 收進布里淵區。(iii) L → ∞ 時點距 → 0,k 變成連續變數,對 k 的求和可以換成 (V/8π³)∫d³k。
典型數值
L = 1 cm:2π/L = 6.28 cm⁻¹;Cu 的 k_F = 1.36×10⁸ cm⁻¹(Table 1),比點距大約 2×10⁷ 倍——沿著 Fermi 球的一條直徑就排了約 4×10⁷ 個允許點。

允許值滿足週期性的驗算

exp[ik_x(x + L)] = exp[i2nπ(x + L)/L] = exp(i2nπx/L)·exp(i2nπ) = exp(i2nπx/L) = exp(ik_x x)

(11)  p.137

k_x = 2nπ/L,n 為整數;x:位置座標。

出處推導
把 (10) 的 k_x = 2nπ/L 代回 (9) 的 x 因子,檢驗 (8) 是否成立;關鍵的一步是 exp(i2nπ) = cos 2nπ + i sin 2nπ = 1。
極限行為
(i) 若 k_x 不是 2π/L 的整數倍,平移 L 會多出一個相因子 exp(ik_x L) ≠ 1,週期性被破壞——這就是允許的 k 必須離散的原因。(ii) n = 0 時兩邊恆等於 1,對應均勻的 k_x = 0 分量。(iii) 同樣的論證在第 5 章 p.110 用於聲子(K = 0, ±2π/Na, …),只是那裡還多了原子離散性造成的上界。

由 (8) 解出 (10)、看懂 (11) 的驗算,並確認 (9) 真的是 (6) 的解補推導

書上 (10) 與 (11) 之間只有三句話,這裡把「由條件解出允許值」與「把允許值代回去驗算」兩個方向都走一次。

第 3 步提醒的「三條獨立條件」很重要:如果只要求 exp[i(k_x + k_y + k_z)L] = 1 這一條,允許的 k 會是一疊平面而不是一個立方格子,後面「每 (2π/L)³ 一個點」的計數就不成立。

第 6 步回答了一個常見的困惑:週期性邊界條件並沒有「禁止」駐波,它只是讓我們可以選行進波當基底。選哪一組基底不影響態的數目(因而不影響任何熱力學量),但決定了哪些態是動量本徵態——這在討論電流時才有差別。

  1. 第 1 步(代入):把 (9) 寫成 exp(ik_x x)·exp(ik_y y)·exp(ik_z z),代入 (8) 的左邊:exp[ik_x(x + L)]·exp[i(k_y y + k_z z)] = exp(ik_x x)·exp(ik_x L)·exp[i(k_y y + k_z z)]。
  2. 第 2 步(約去):(8) 要求它等於 exp(ik_x x)·exp[i(k_y y + k_z z)]。兩邊共同的因子處處不為零,可以直接約去,剩下 exp(ik_x L) = 1。
  3. 第 3 步(解):複數指數等於 1 的充要條件是幅角為 2π 的整數倍:k_x L = 2πn,n = 0, ±1, ±2, …。所以 k_x = 0, ±2π/L, ±4π/L, …,這就是 (10)。y、z 方向各自獨立地再做一次——三個方向是三條獨立的週期性要求,不是一條加總的條件。
  4. 第 4 步(讀懂 (11)):(11) 是反方向的驗算。令 k_x = 2nπ/L,則 exp[ik_x(x + L)] = exp[i2nπ(x + L)/L] = exp(i2nπx/L)·exp(i2nπ)。關鍵在於 exp(i2nπ) = cos 2nπ + i sin 2nπ = 1,所以結果回到 exp(i2nπx/L) = exp(ik_x x),週期性成立。
  5. 第 5 步((9) 確實解 (6)):∂²/∂x² exp(ik·r) = (ik_x)² exp(ik·r) = −k_x² exp(ik·r),y、z 同理,所以 ∇²ψ_k = −(k_x² + k_y² + k_z²)ψ_k = −k²ψ_k。代入 (6) 左邊:−(ħ²/2m)(−k²)ψ_k = (ħ²k²/2m)ψ_k。這對任何 k 都成立,本徵值 ħ²k²/2m 就是 (12)。
  6. 第 6 步(為什麼不用 sin、cos):在週期性條件下,cos(k_x x) 與 sin(k_x x)(k_x = 2πn/L,n > 0)也是合法的解——它們是環上的駐波。但它們只是 exp(±ik_x x) 的線性組合:每個 n > 0 要嘛數成 {cos, sin} 兩個,要嘛數成 {exp(+ik_x x), exp(−ik_x x)} 兩個,態數完全相同。第 5 章 p.110 對聲子說過同一件事:改用複數形式 exp(isKa) 時,週期性條件給出的模態數不變。選指數形式的理由是它同時是動量本徵態((13))。

週期性邊界條件到底假設了什麼、什麼時候會失效我的判讀

初學者看到 (8) 常會問:真實的金屬塊明明有表面,憑什麼假設 ψ(x + L) = ψ(x)?

它不是在說樣品真的會自我重複,而是一個計數上的約定。一維時可以想成把線段首尾接成一個環(第 5 章 p.110 的 Fig. 4 就是把原子鏈接成環);三維則是把立方體的三組對面兩兩黏合,成為一個沒有邊界的「三維環面」。沒有牆壁就沒有反射,行進波可以一直走下去而不被逼成駐波,+k 與 −k 兩個方向因此都保留下來。

為什麼可以這麼放心?我們關心的體性質——每單位體積、每單位能量有幾個軌域——只取決於 k 空間中允許點的密度,而點的密度只取決於 V。邊界條件的差別只影響緊貼表面的那一層態,其數目正比於表面積 ∝ L²,相對於體積項 ∝ L³ 是 1/(k_F L) 量級的修正。下一個推導區塊會實際示範:駐波與週期性條件數出完全相同的主項。

兩個容易讀過去的細節。其一,p.137 行文說 n 是「正或負整數」,但 (10) 明確列出 k_x = 0——n = 0 當然也允許;三個分量都為零就是均勻的 k = 0 軌域(第 5 章 p.110 的 Fig. 5 也把 K = 0 標為均勻模態)。其二,自旋是作為第四個獨立量子數 m_s 加進來的:(12) 的能量與自旋無關,所以每個允許的 k 都對應兩個能量相同的軌域,這就是 (15) 左邊因子 2 的來源。

什麼時候會失效?當樣品小到 k_F L 只有幾十,表面項就不再可忽略,真實的邊界必須認真處理。第 18 章 p.520 式 (3) 指出有限尺寸系統的態密度是一串 δ 函數;p.528 談矩形奈米線時,截面方向的能量與本徵態直接用第 6 章的盒中粒子解——也就是本段的 (7),而不是 (9)。換句話說,(7) 與 (9) 的分工是:巨觀樣品用 (9) 比較方便;奈米結構在受限的方向上必須回到 (7)。

駐波與行進波數出同樣多的態:一維對回式 (4),三維對回 1/8 卦限補推導

這是本段最值得親手做一次的計數。書中只說「方便起見」改用週期性條件,並沒有示範兩種邊界條件給出同樣的結果;第 5 章 p.110 則在一維聲子鏈上,用 8 個原子的例子示範過同一件事(固定端與週期性條件都給出每個粒子一個模態)。這裡把它推廣到三維電子。

結論有兩層意思。第一,態數(因而所有熱力學量)與邊界條件的選擇無關,至少在巨觀樣品中如此,所以可以挑最方便的那一種。第二,兩者的差別不在計數,而在態的性質:駐波是 ±k 的疊加、不帶動量;行進波各自帶 ±ħk。計數時隨便挑,討論輸運時必須用行進波。

  1. 第 1 步(三維駐波的 k 格子):由 (7),k_i = πn_i/L,n_i = 1, 2, 3, …。允許點構成晶格常數 π/L 的立方格子,每點佔 (π/L)³ = π³/V 的體積,但只在第一卦限(k_x, k_y, k_z > 0)。
  2. 第 2 步(駐波數到半徑 k):半徑 k 的球只有 1/8 落在第一卦限,所以點數 = (1/8)·(4πk³/3)/(π³/V) = (4πk³/3)·V/(8π³) = Vk³/6π²。
  3. 第 3 步(週期性條件的 k 格子):由 (10),k_i = 2πn_i/L,n_i 為任意整數。晶格常數 2π/L,每點佔 (2π/L)³ = 8π³/V,但整顆球都能用。
  4. 第 4 步(週期性條件數到半徑 k):點數 = (4πk³/3)/(8π³/V) = Vk³/6π²——與第 2 步完全相同。
  5. 第 5 步(為什麼剛好相同):駐波的格子每個方向密一倍(π/L 對 2π/L),體積密度是 2³ = 8 倍;但它只用 1/8 的球。物理上,sin(kx) ∝ e^{ikx} − e^{−ikx} 把 +k 與 −k 合併成一個態,所以駐波把 ±k 只數一次,而格點加密剛好補回來。
  6. 第 6 步(一維驗算,對回 p.135 式 (4)):駐波:k = nπ/L,n = 1…n_F,且 2n_F = N,故 k_F = n_Fπ/L = Nπ/2L。週期性:k = 2πn/L,|k| ≤ k_F 的點數是 2k_F/(2π/L) = k_F L/π,乘上自旋 2 得 2k_F L/π = N,故 k_F = Nπ/2L。兩者都給出 ε_F = (ħ²/2m)(Nπ/2L)²,正是式 (4)。
  7. 第 7 步(差在哪裡):駐波排除了 n_i = 0 的三個座標面,週期性條件則把整顆球都算進去;兩者只差在靠近座標面、數目正比於 (kL)² 的一層點,相對於主項 (kL)³ 是 1/(kL) 量級(係數為個位數)。取 Cu 的 k_F = 1.36×10⁸ cm⁻¹、L = 1 cm,1/(k_F L) ≈ 7×10⁻⁹,完全可忽略;若 L = 5 nm,k_F L ≈ 68,這一層的相對份量就到了一成左右的等級,不能再忽略。

自由電子軌域的能量(色散關係)

ε_k = (ħ²/2m)k² = (ħ²/2m)(k_x² + k_y² + k_z²)

(12)  p.138

k = |k|:波向量大小,與波長的關係為 k = 2π/λ(原文同頁);ε_k:軌域 k 的能量。

出處推導
把 (9) 代入 (6)(見上方推導第 5 步)。
極限行為
(i) 能量只依賴 |k|、不依賴方向:等能面是以原點為中心的球面——這就是基態被佔據區域為一顆球(Fig. 4)的原因。(ii) 拋物線色散是自由電子的特徵;在半導體帶底附近仍可寫成拋物線,只是換成有效質量並把能量原點移到帶邊:第 8 章 p.205 式 (36) ε_k = E_c + ħ²k²/2m_e。因此本段 (20) 的 √ε 形式會在半導體中保留下來。(iii) 相對論極限:當 ħk 接近 mc,(12) 失效;p.158 習題 4(b) 改用 ε = ħkc。
典型數值
ħ²/2m = 3.81 eV·Å²;k = 1 Å⁻¹ 時 ε = 3.81 eV。代入 Cu 的 k_F = 1.36 Å⁻¹(Table 1 的三位有效數字)得 7.05 eV;用較精確的 1.3576 Å⁻¹(k_F² = 1.843 Å⁻²)得 7.02 eV,表列 7.00 eV。

平面波是線動量的本徵函數

pψ_k(r) = −iħ∇ψ_k(r) = ħkψ_k(r)

(13)  p.138

p = −iħ∇:線動量算符;ħk:本徵值。原文隨後說平面波 ψ_k 是線動量的本徵函數、本徵值 ħk,並在行文中寫出軌域 k 的粒子速度 v = ħk/m(該式無編號,見下一個區塊)。

出處推導
對 (9) 作用梯度:∇exp(ik·r) = ik·exp(ik·r),再乘上 −iħ:(−iħ)(ik) = ħk。
極限行為
(i) 駐波 (7) 不是本徵態:它是 ±ħk 的疊加,動量量測會以各一半的機率得到 ±ħk,期望值為零。(ii) 基態 Fermi 球對原點對稱,每個 k 都有 −k 配對,總動量 Σħk = 0,沒有電流;p.147 起讓電場把整顆 Fermi 球平移,才產生淨動量與電流。(iii) 在晶體中,第 7 章 p.173 把 ħk 稱為電子的晶體動量(crystal momentum):k 出現在碰撞過程的守恆律中(書中說這些守恆律其實是躍遷的選擇定則),例如電子吸收波向量 q 的聲子時 k + q = k′ + G,可以多出一個倒晶格向量 G——這是它與自由電子的真實動量 ħk 不同之處。
典型數值
Cu:ħk_F = 1.0546×10⁻²⁷ erg·s × 1.358×10⁸ cm⁻¹ = 1.43×10⁻¹⁹ g·cm/s。

軌域 k 中的粒子速度

v = ħk/m

p.138

v:處於軌域 k 之電子的速度(向量)。此式在原文中寫於行文內,沒有式號。

出處推導
由 (13) 的動量本徵值 p = ħk 與非相對論關係 p = mv。等價地,也可由群速度得到:第 8 章 p.191 式 (1) 寫 v = ħ⁻¹∇_k ε,對 (12) 取 k 梯度得 ħ⁻¹·ħ²k/m = ħk/m,兩者一致。
極限行為
(i) 平面波的相速度 (ε/ħ)/k = ħk/2m 只有粒子速度的一半——波峰移動的速度不是電子的速度;電子速度是群速度。這對自由電子已經成立,對能帶電子則是唯一正確的寫法。(ii) 在 Fermi 面上 |v| = v_F,即 (18)。(iii) 能帶電子的速度由 ħ⁻¹∇_k ε(k) 給出,一般不再平行於 k(第 8 章)。
典型數值
Cu:v_F = 1.57×10⁸ cm/s(Table 1),約為光速的 0.5%,非相對論處理是安全的。

平面波帶著多少電流:歸一化、機率流與 Fermi 球的總動量補推導

(13) 與其後的 v = ħk/m 在書上只佔兩行;這裡補上把「動量本徵態」翻譯成「電流」的幾步。這幾步在本章後半的電導率(p.147–148)與 Hall 效應都會默默用到。

第 5 步還說明了一件後面很重要的事:電場雖然讓球內每個電子都平移同樣的 δk,但佔據分布的淨改變只出現在 Fermi 面附近——球的前緣多了一層被佔據的態、後緣少了一層,球內深處的電子,其 ±k 對手都在,彼此抵消。

  1. 第 1 步(歸一化,書中未寫):(9) 沒有歸一化常數。在週期性條件下 |exp(ik·r)|² = 1,∫_V |ψ_k|² d³r = V,所以歸一化的軌域是 ψ_k = V^{−1/2}exp(ik·r):一個電子均勻攤在整個樣品上。
  2. 第 2 步(機率流):量子力學的機率流密度 j = (ħ/2mi)(ψ*∇ψ − ψ∇ψ*)。代入 ψ_k = V^{−1/2}e^{ik·r}:ψ*∇ψ = ik/V,ψ∇ψ* = −ik/V,相減得 2ik/V,所以 j = (ħ/2mi)(2ik/V) = ħk/(mV)。
  3. 第 3 步(解讀):j = (1/V)·(ħk/m)——密度為 1/V 的「一個電子」以速度 ħk/m 流動。這就是 p.138「軌域 k 中的粒子速度 v = ħk/m」的機率流版本。
  4. 第 4 步(駐波的機率流為零):對實函數 ψ = A sin(πn_x x/L)…,ψ* = ψ,ψ*∇ψ − ψ∇ψ* ≡ 0,所以 j = 0。駐波不帶電流,與前面的論證一致。
  5. 第 5 步(Fermi 球的總動量):基態填滿以原點為中心的球;對每個被佔據的 k,−k 也在球內且同樣被佔據,所以 Σħk = 0、Σj = 0:平衡態的金屬沒有淨電流。要有電流,佔據分布就必須不再對原點對稱——p.147 起的做法正是讓電場把整顆 Fermi 球平移 δk,平移後總動量為 Nħδk。
  6. 第 6 步(相速度 ≠ 粒子速度):平面波的時間因子是 exp(−iεt/ħ),相速度 = (ε/ħ)/k = ħk/2m,只有 ħk/m 的一半;群速度 dω/dk = (1/ħ)dε/dk = ħk/m 才等於粒子速度——這正是第 8 章 p.191 式 (1) 以 v = ħ⁻¹∇_k ε 定義電子速度的理由。

基態是一顆 Fermi 球(p.138 與 Fig. 4)

p.138 把 N 個自由電子的基態畫成 k 空間中的一個幾何物體:被佔據的軌域可以表示為 k 空間中一個球內部的點;球面上的能量就是 Fermi 能,位於 Fermi 面上的波向量大小為 k_F,滿足 (14)。

接著是計數:由 (10) 可見,k 空間中每 (2π/L)³ 的體積元恰有一個允許的波向量——也就是一組不同的量子數 (k_x, k_y, k_z)。於是在體積 4πk_F³/3 的球內,軌域總數由 (15) 給出,左邊的因子 2 來自每個允許的 k 都有兩個自旋量子數值。由 (15) 解出 (16),Kittel 特別指出 k_F 只依賴粒子濃度。

Fig. 4(p.138)在說什麼(網頁無圖,以文字還原):三條互相垂直的座標軸 k_x、k_y、k_z 交於一點,以此點為中心畫了一顆灰色的實心球(表面帶有淡淡的斜向虛線陰影)。一支箭頭從球心指向球面,標為 k_F;右側一條引線指向球面,標注「Fermi 面,位於能量 ε_F」。圖說的意思是:在 N 個自由電子系統的基態,被佔據的軌域填滿半徑 k_F 的球,而 ε_F = ħ²k_F²/2m 是波向量為 k_F 之電子的能量。

為什麼是球?為什麼整個故事的主角是 Fermi 面?我的判讀

為什麼是球。基態就是總能量最低的狀態。Pauli 原理規定每個軌域(一個 k 加一個自旋值)至多一個電子,所以要放進 N 個電子、又要總能量最小,唯一的辦法是從能量最低的軌域開始往上填,直到 N 個填完。由 (12),自由電子的能量只依賴 |k|、不依賴方向,所以「能量低於某值的所有 k」恰好構成一顆以原點為中心的實心球。在 k 空間的解析度(2π/L)下,球的邊界其實是鋸齒狀的,但對 N ~ 10²³,這種鋸齒完全看不見。

如果能量依賴方向呢?在晶體中 ε(k) 通常與方向有關,等能面不再是球,被佔據區域的邊界也跟著變形——但「佔據區與空區之間有一個清楚的邊界曲面」這個概念保留下來,仍叫 Fermi 面。第 9 章 p.249 以銅為例:自由電子 Fermi 球的直徑是 9.80/a,而布里淵區兩個六角形面之間的距離是 10.88/a,球碰不到區界;但真實銅的 Fermi 面長出頸部去接觸六角形面。

為什麼主角是 Fermi 面,不是整顆球。在絕對零度,球內每個軌域都滿、球外都空。任何微小的擾動——一點熱能、一點電場——想改變某個電子的狀態,都必須把它送進一個空的軌域。球深處的電子周圍全是滿的,沒有地方可去;只有 Fermi 面附近的電子才找得到鄰近的空位。這個觀察貫穿整章:

  • 熱容:只有 Fermi 面附近約 k_BT 範圍內的電子能被熱激發(p.141,對應 Fig. 5 的區域 1 → 2);
  • 導電:電場讓整顆球平移,但佔據分布的淨改變只發生在 Fermi 面上(p.147–148);
  • 散射:p.134 說傳導電子很少被其他傳導電子散射,是 Pauli 原理的後果——同樣是因為能讓電子「落腳」的空位只在 Fermi 面附近。

所以 k_F、v_F、D(ε_F) 這三個「Fermi 面上的量」,會比 ε 的整個分布更常出現在後面的公式裡。

Fermi 能與 Fermi 波向量

ε_F = (ħ²/2m)k_F²

(14)  p.138

ε_F:Fermi 能——基態下最高填滿軌域的能量(p.135 的定義;p.137 並指出絕對零度時化學勢 μ 等於 ε_F);k_F:Fermi 波向量的大小,即 Fermi 球的半徑。

出處推導
把 (12) 用在 Fermi 面上(|k| = k_F)。它本身只是色散關係的一個特例;真正的內容在下一步——k_F 由電子數決定。
極限行為
(i) ε_F 按定義是 T = 0 的量,不隨溫度改變;T > 0 時變動的是化學勢 μ(p.144 Fig. 8:三維略降),圖說指出一般金屬在室溫 τ/ε_F 約 0.01,所以 μ 與 ε_F 非常接近。(ii) 能帶電子的 ε(k) 不再是 ħ²k²/2m:Fermi 能仍定義為最高填滿能量,但 (14) 的形式失效;帶底附近可用 ħ²k²/2m* 延續。(iii) 與一維式 (4) 對照:形式完全相同,差別只在 k_F 與 N 的關係(一維 k_F = Nπ/2L,三維見 (16))。
典型數值
Cu:由 n 算得 k_F = 1.3576 Å⁻¹,ε_F = 3.810 eV·Å² × 1.843 Å⁻² = 7.02 eV;Table 1 列 7.00 eV。

(15) 的每一個因子:從「一點一格」數到 N = Vk_F³/3π²,再解出 (16)補推導

(15) 在書上是一行等式,三個等號各藏了一個動作:最左邊「數點」、中間「化簡」、最右邊「令它等於電子數」。

第 2 步的近似值得停下來想一次:它要求球內點數夠多。這在巨觀金屬中毫無問題(1 cm³ 的 Cu 有 8.45×10²² 個傳導電子),但對只有幾十、幾百個電子的金屬團簇就不成立——第 18 章 p.547 描述的「magic numbers」團簇殼層結構,正是球體積近似失效、離散能階的細節浮上檯面的例子。

  1. 第 1 步(格子):由 (10),允許的 k = (2π/L)(n_x, n_y, n_z),n_i 為任意整數。這些點在 k 空間排成晶格常數 2π/L 的簡單立方格子;把空間切成以每個點為中心的小立方格,每格恰含一點,體積 (2π/L)³ = 8π³/L³ = 8π³/V。
  2. 第 2 步(把數點換成算體積):當 k_F 遠大於格距 2π/L(巨觀樣品必然如此),球內點數 ≈ 球體積 ÷ 單格體積;誤差只來自被球面切開的那一層格子,相對大小 ~ 2π/(k_F L)。
  3. 第 3 步(k 點數):(4πk_F³/3) ÷ (8π³/V) = (4πk_F³/3)·V/(8π³) = Vk_F³/6π²。其中 4/(3·8) = 1/6,π/π³ = 1/π²。
  4. 第 4 步(乘上自旋):每個 k 有 m_s = ±½ 兩個軌域,每個軌域至多一個電子,基態全部填滿:N = 2 × Vk_F³/6π² = Vk_F³/3π²。這就是 (15):2·(4πk_F³/3)/(2π/L)³ = (V/3π²)k_F³ = N。
  5. 第 5 步(解出 k_F):k_F³ = 3π²N/V,所以 k_F = (3π²N/V)^{1/3},即 (16)。
  6. 第 6 步(反過來用):n = N/V = k_F³/3π²;Fermi 球的赤道截面積為 πk_F²。第 9 章 p.249 就是這樣用的:由 Table 1 取金的 k_F = 1.2×10⁸ cm⁻¹,算出極值軌道面積約 4.5×10¹⁶ cm⁻²,與磁矩振盪週期的實驗值大致相符(π × (1.2×10⁸)² = 4.5×10¹⁶ ✓)。
  7. 第 7 步(量綱檢查):N/V 的單位是 cm⁻³,開立方得 cm⁻¹,正是波向量的單位;3π² 是純數。

Fermi 球內的軌域總數

2 · (4πk_F³/3)/(2π/L)³ = (V/3π²)k_F³ = N

(15)  p.138

4πk_F³/3:Fermi 球的體積;(2π/L)³:每個允許 k 所佔的 k 空間體積;2:每個允許 k 的兩個自旋量子數值(原文明寫);V = L³:樣品體積;N:電子總數。

出處推導
(10) 給出的「每 (2π/L)³ 一個允許 k」+ Pauli 原理(每個軌域至多一個電子)+ 基態從最低能量往上填滿(逐步推導見上一個區塊)。
極限行為
(i) 「球體積 ÷ 單格體積」只在球內點數很多時成立;對 N ~ 10²³ 的巨觀樣品相對誤差 ~N^{−1/3},完全可忽略;對只含數百個電子的團簇,會被殼層結構取代(第 18 章 p.547)。(ii) 若自旋被完全極化(每個 k 只放一個電子),因子 2 變為 1,同樣的 N 需要大 2^{1/3} ≈ 1.26 倍的 k_F。(iii) 長方體樣品 L_x×L_y×L_z 的單格體積是 8π³/(L_xL_yL_z) = 8π³/V,結果不變——(15) 只依賴 V,與樣品形狀無關。
典型數值
Cu,V = 1 cm³:k_F = 1.358×10⁸ cm⁻¹,k_F³ = 2.502×10²⁴ cm⁻³,Vk_F³/3π² = 2.502×10²⁴/29.61 = 8.45×10²²,回到表列的電子數(自洽)。

Fermi 波向量

k_F = (3π²N/V)^{1/3}

(16)  p.138

N/V = n:傳導電子濃度。原文在此式後指出 k_F 只依賴粒子濃度。

出處推導
解 (15):k_F³ = 3π²N/V。
極限行為
(i) 式中沒有 m 也沒有 ħ——k_F 是純粹的計數(幾何)結果,質量只在把 k_F 換成 ε_F、v_F 時才出現。(ii) k_F r_0 = (9π/4)^{1/3} ≈ 1.92(r_0 為每個電子分到的球半徑,推導見 r_s 區塊):Fermi 波長 λ_F = 2π/k_F 與電子間距同一尺度。(iii) 其他維度:一維 k_F = Nπ/2L(p.135 式 (4) 中的組合);二維同法可得 k_F = (2πn)^{1/2}(n 為面密度,本講解補充)。(iv) 濃度差很多時 k_F 只差立方根:濃度 → 0 時 k_F → 0,但真正拿來和 k_BT 比較的是 ε_F(見 (17))。
典型數值
Table 1:Cs 0.64×10⁸、Cu 1.36×10⁸、Be 1.93×10⁸ cm⁻¹。Be 與 Cs 的濃度差 24.2/0.91 = 26.6 倍,k_F 只差 26.6^{1/3} ≈ 2.98 倍(表列 1.93/0.64 = 3.0)。

k_F 只看濃度:L 為什麼消失,以及 k_F 為什麼和晶格同一個尺度我的判讀

p.138 在 (16) 之後只補了一句:k_F 只依賴粒子濃度。這句話值得拆成三層來讀。

第一層:L 為什麼會消失。球內的點數正比於 V(因為單格體積 8π³/V 反比於 V),電子數 N 也正比於 V;兩者相等時 V 被約掉,只剩比值 N/V。k_F、ε_F、v_F、T_F 因此全是強度量,與樣品大小、形狀無關;而 N 與後面的 D(ε) 是外延量,正比於 V。這個分類在後面檢查公式時很有用:凡是把強度量與外延量等號相連的式子,一定漏了 V 或 N。

第二層:k_F 不含質量。(16) 裡沒有 m,也沒有 ħ。質量只在把 k_F 翻譯成能量與速度時才出現:ε_F ∝ 1/m、v_F ∝ 1/m。這個分工在後面很有用:p.145 用「熱有效質量」m_th 描述實測電子熱容係數 γ 與自由電子值的比,理由正是 ε_F 與電子質量成反比、因而 γ ∝ m;被質量修正改寫的是 ε_F 與 D(ε_F),而(對球形 Fermi 面)k_F 仍由濃度決定。

第三層:k_F 與電子間距是同一個尺度。若把每個電子分到的體積寫成半徑 r_0 的球,(4π/3)r_0³ = 1/n,代入 (16) 得 k_F r_0 = (9π/4)^{1/3} ≈ 1.92,也就是 λ_F = 2π/k_F ≈ 3.27 r_0。以 Cu 為例,λ_F = 2π/(1.358×10⁸ cm⁻¹) ≈ 4.6 Å;若依 p.249 對單價 fcc 金屬的寫法 n = 4/a³,由 Table 1 的 n 反推 a = (4/n)^{1/3} ≈ 3.62 Å——Fermi 波長與晶格常數同數量級。這不是巧合:每個原子貢獻約一個電子,電子間距自然就是原子間距。其後果是 Fermi 球的大小必然與布里淵區相當(p.249:Cu 的 k_F = 4.90/a,最近的六角形區界面則在 10.88/(2a) = 5.44/a 處),所以晶格位能對 Fermi 面附近的電子不可能完全沒有影響——這是第 7–9 章能帶論的伏筆。

Fermi 能與電子濃度

ε_F = (ħ²/2m)(3π²N/V)^{2/3}

(17)  p.140

符號同 (14)、(16)。書中以方框標示此式,並評論它把 Fermi 能與電子濃度 N/V 聯繫起來。

出處推導
(14) 與 (16) 相結合(原文一句帶過;代入步驟見下方推導區塊)。
極限行為
(i) ε_F ∝ n^{2/3}:Table 1 中 Be 與 Cs 濃度差 26.6 倍,ε_F 應差 26.6^{2/3} ≈ 8.9 倍,表列 14.14/1.58 = 8.95。(ii) ε_F ∝ 1/m:p.145 以此說明為什麼電子熱容係數 γ ∝ m,進而以 γ(實測)/γ(自由) 定義熱有效質量 m_th/m。(iii) 稀薄極限:n = 10¹⁶ cm⁻³(輕摻雜半導體的量級)時 ε_F ≈ 1.7×10⁻⁴ eV(自由電子質量估計),遠小於室溫 k_BT ≈ 0.026 eV,電子氣變成古典的 Boltzmann 氣體——就是 p.137 所說的高能尾極限。(iv) 極高密度:ε_F 大到 v_F 接近光速時 (12) 失效,p.158 習題 4 改用 ε = ħkc,Fermi 能約為 ħc(N/V)^{1/3}。(v) 二維:由 p.157 習題 3 附註的面態密度 m/πħ²,得 ε_F = πħ²n/m,正比於 n 的一次方;一維則 ∝ (N/L)²(式 (4))。
典型數值
實用式 ε_F = 1.693 eV × (n/10²² cm⁻³)^{2/3}(推導見下)。Cu:1.693 × 8.45^{2/3} = 1.693 × 4.149 = 7.02 eV;表列 7.00 eV。

Fermi 速度

v_F = (ħk_F/m) = (ħ/m)(3π²N/V)^{1/3}

(18)  p.140

v_F:位於 Fermi 面上之電子的速度大小。

出處推導
軌域 k 的粒子速度 v = ħk/m(p.138)取 |k| = k_F,再代入 (16)。
極限行為
(i) 與溫度無關,0 K 時依然存在:它是 Pauli 原理強加的「零點運動」,不是熱運動。Cu 的 v_F 比 300 K 的古典熱速度 (3k_BT/m)^{1/2} ≈ 1.2×10⁷ cm/s 大約 13 倍。(ii) v_F ∝ n^{1/3}/m。(iii) v_F/c:Cu 約 0.5%,表中最大的 Be(2.23×10⁸ cm/s)約 0.7%,金屬中都可用非相對論處理;p.158 習題 4 的白矮星與脈衝星則是 v_F 逼近 c 的反例。(iv) 凡是「電子多快跑到下一次碰撞」的問題,速度都是 v_F:p.156 以 ℓ = v_Fτ 寫平均自由徑,p.159 習題 10 以 τ ≈ d/v_F 估計薄膜的表面散射。
典型數值
Cu:1.57×10⁸ cm/s = 1.57×10⁶ m/s。實用式 v_F = 7.72×10⁷ cm/s × (n/10²² cm⁻³)^{1/3}。

由 (14)(16) 到 (17)(18),以及三條可以心算的實用式補推導

這一區塊的目的,是讓你以後看到任何金屬的電子濃度,都能在腦中直接估出 ε_F、k_F、v_F。第 5 步的三條式子只用到 (16)–(18) 與物理常數,沒有任何新物理;它們在自由電子質量、兩個自旋、單一球形 Fermi 面的前提下成立。

  1. 第 1 步((17)):把 (16) 平方:k_F² = (3π²N/V)^{2/3}。代入 (14) ε_F = (ħ²/2m)k_F²,即得 ε_F = (ħ²/2m)(3π²N/V)^{2/3}。書中只以一句「由 (14) 與 (16)」帶過,實際就是這一步代入。
  2. 第 2 步((18)):把 |k| = k_F 代入 v = ħk/m,再用 (16):v_F = ħk_F/m = (ħ/m)(3π²N/V)^{1/3}。
  3. 第 3 步(兩者的關係):由 (14) 與 (18),ε_F = ħ²k_F²/2m = ½m(ħk_F/m)² = ½mv_F²。p.156 推導金屬的電子熱導率時正是用這個關係,把熱容與 ℓ = v_Fτ 組合起來。
  4. 第 4 步(數值常數):ħ²/2m = 3.810 eV·Å²,(3π²)^{1/3} = 3.094,(3π²)^{2/3} = 9.571。以 10²² cm⁻³ = 10⁻² Å⁻³ 為 n 的單位,(10⁻²)^{1/3} = 0.2154 Å⁻¹,(10⁻²)^{2/3} = 0.04642 Å⁻²。
  5. 第 5 步(三條實用式):k_F = 3.094 × 0.2154 × (n/10²²)^{1/3} Å⁻¹ = 0.666 Å⁻¹ × (n/10²² cm⁻³)^{1/3};ε_F = 3.810 × 9.571 × 0.04642 × (n/10²²)^{2/3} eV = 1.693 eV × (n/10²² cm⁻³)^{2/3};v_F = (ħ/m)k_F,ħ/m = 1.158 cm²/s,得 v_F = 7.72×10⁷ cm/s × (n/10²² cm⁻³)^{1/3}。
  6. 第 6 步(驗一下):Cu,n = 8.45:8.45^{1/3} = 2.037,8.45^{2/3} = 4.149。k_F = 1.357 Å⁻¹,ε_F = 7.02 eV,v_F = 1.572×10⁸ cm/s,與 Table 1 的 1.36×10⁸ cm⁻¹、7.00 eV、1.57×10⁸ cm/s 一致。
  7. 第 7 步(標度律):k_F ∝ n^{1/3};ε_F ∝ n^{2/3}/m;v_F ∝ n^{1/3}/m。濃度加倍,ε_F 只增加 2^{2/3} ≈ 1.59 倍。

Table 1(p.139):自由電子 Fermi 面參數表

p.140 說,Table 1 列出部分金屬的 k_F、v_F、ε_F 計算值,另外也列出定義為 ε_F/k_B 的 T_F。表本身印在 p.139(橫排整頁,有文字層,數值可直接擷取並與算繪影像核對)。

標題與條件。表題是金屬在室溫下的自由電子 Fermi 面參數「計算值」;括號註明例外:Na、K、Rb、Cs 取 5 K、Li 取 78 K 的數據。

欄位(由左至右):價數;金屬;電子濃度(cm⁻³);半徑參數 r_n(上標 a 指向表底註腳);Fermi 波向量(cm⁻¹);Fermi 速度(cm s⁻¹);Fermi 能(eV);Fermi 溫度 T_F ≡ ε_F/k_B(K)。全表用 CGS 單位。

讀表的小陷阱:10 的冪次只印在第一列(Li)——4.70×10²²、1.11×10⁸、1.29×10⁸、5.48×10⁴——以下各列只印係數,冪次沿用同欄第一列。所以 Cu 列的「8.45、1.36、1.57、8.12」要讀成 8.45×10²² cm⁻³、1.36×10⁸ cm⁻¹、1.57×10⁸ cm/s、8.12×10⁴ K;ε_F 欄本來就以 eV 為單位,r_n 欄無因次,都沒有冪次。

分組:價數 1:Li、Na、K、Rb、Cs、Cu、Ag、Au;價數 2:Be、Mg、Ca、Sr、Ba、Zn、Cd;價數 3:Al、Ga、In;價數 4:Pb、Sn(w)。

數值範圍(全表的兩端都是 Cs 與 Be):

量Cs(最稀)Be(最密)
電子濃度 n(10²² cm⁻³)0.9124.2
r_n5.631.88
k_F(10⁸ cm⁻¹)0.641.93
v_F(10⁸ cm/s)0.752.23
ε_F(eV)1.5814.14
T_F(10⁴ K)1.8316.41

註腳 a:無因次半徑參數定義為 r_n = r_0/a_H,其中 a_H 是第一 Bohr 半徑,r_0 是含有一個電子的球的半徑。

Cu 列逐項驗算:把 n = 8.45×10²² cm⁻³ 代進 (16)–(18)補推導

Table 1 同一列的五個計算量(k_F、v_F、ε_F、T_F、r_n)全部只由一個輸入 n 決定。逐項代入後,k_F、v_F、r_n 三位有效數字全對;ε_F 與 T_F 比重算值低約 0.3%。

這個小差異在全表幾乎每一列都以同一方向出現(下一區塊),最可能的原因是編表時採用的物理常數或中間值捨入不同,而不是誤植——因為表中 T_F/ε_F 的比值始終在 11,600 K/eV 左右,表的內部是自洽的。

  1. 常數(CGS):ħ = 1.0546×10⁻²⁷ erg·s,m = 9.109×10⁻²⁸ g,k_B = 1.3807×10⁻¹⁶ erg/K,1 eV = 1.6022×10⁻¹² erg,a_H = 0.5292×10⁻⁸ cm。
  2. k_F(式 (16)):3π² = 29.609;3π²n = 29.609 × 8.45×10²² = 2.502×10²⁴ cm⁻³;開立方:(2.502)^{1/3} = 1.3576,(10²⁴)^{1/3} = 10⁸,所以 k_F = 1.358×10⁸ cm⁻¹。表列 1.36×10⁸ ✓。
  3. v_F(式 (18)):ħ/m = 1.0546×10⁻²⁷ / 9.109×10⁻²⁸ = 1.1577 cm²/s;v_F = 1.1577 × 1.3576×10⁸ = 1.572×10⁸ cm/s。表列 1.57×10⁸ ✓。
  4. ε_F(式 (17),用 ε_F = ½mv_F² 最省事):v_F² = 2.470×10¹⁶ cm²/s²;½mv_F² = 0.5 × 9.109×10⁻²⁸ × 2.470×10¹⁶ = 1.125×10⁻¹¹ erg;除以 1.6022×10⁻¹² erg/eV 得 7.02 eV。表列 7.00 eV(低 0.3%)。
  5. T_F:ε_F/k_B = 1.125×10⁻¹¹ / 1.3807×10⁻¹⁶ = 8.15×10⁴ K。表列 8.12×10⁴ K(低 0.3%,與 ε_F 的差同源)。
  6. r_n:(4π/3)r_0³ = 1/n,所以 r_0 = [3/(4πn)]^{1/3} = [3/(1.0619×10²⁴)]^{1/3} = (2.825×10⁻²⁴)^{1/3} = 1.414×10⁻⁸ cm;r_n = 1.414/0.5292 = 2.67。表列 2.67 ✓。
  7. 交叉檢查:T_F/ε_F 應等於 1/k_B = 11,604 K/eV;表中 8.12×10⁴/7.00 = 11,600 K/eV ✓。又 v_F/c = 1.57×10⁸/3.00×10¹⁰ = 0.52%,非相對論近似成立。

全表 20 列驗算:誤差都在 1% 以內補推導

用與上一區塊相同的常數,對 Table 1 全部 20 列由 n 重算 k_F、v_F、ε_F、T_F、r_n。下表列出 10 個代表性金屬(每格為「表列/重算」):

金屬n(10²² cm⁻³)k_F(10⁸ cm⁻¹)v_F(10⁸ cm/s)ε_F(eV)T_F(10⁴ K)r_n
Li4.701.11/1.1161.29/1.2924.72/4.755.48/5.513.25/3.25
Na2.650.92/0.9221.07/1.0683.23/3.243.75/3.763.93/3.93
K1.400.75/0.7460.86/0.8632.12/2.122.46/2.464.86/4.86
Cs0.910.64/0.6460.75/0.7481.58/1.591.83/1.845.63/5.61
Cu8.451.36/1.3581.57/1.5727.00/7.028.12/8.152.67/2.67
Au5.901.20/1.2041.39/1.3945.51/5.536.39/6.413.01/3.01
Be24.21.93/1.9282.23/2.23214.14/14.1616.41/16.431.88/1.88
Mg8.601.37/1.3661.58/1.5817.13/7.108.27/8.242.65/2.66
Al18.061.75/1.7492.02/2.02411.63/11.6513.49/13.522.07/2.07
Pb13.201.57/1.5751.82/1.8249.37/9.4510.87/10.972.30/2.30

結論:

  • k_F 與 v_F:20 列全部在 ±1% 以內(最大為 Cs 的 k_F,−0.9%),絕大多數只是最後一位的捨入差。
  • ε_F 與 T_F:只有 Mg、K、Sr 三列的表列值略高(Mg 約 +0.3%,另兩列約 +0.1%);其餘 17 列都比重算值略低,方向一致,多在 0.03–0.6% 之間;偏差最大的是 Pb(約 −0.9%),其次 Ga 與 Cs(T_F 約 −0.8%)、Ba(ε_F −0.7%)。
  • r_n:除 Cs 差 0.015 外,其餘 19 列與重算值的差都不超過 0.01。
  • 內部自洽:各列 T_F/ε_F 都落在 11,580–11,625 K/eV,與 1 eV/k_B = 11,604.5 K 一致。

這些差異都遠小於自由電子模型本身的不確定性——p.145–147 列出三種會讓實測量偏離自由電子值的效應(週期晶格位能造成的能帶有效質量、電子與聲子的交互作用、電子彼此的交互作用)。沒有任何一列像是誤植,可以把 Table 1 放心當成「由 n 經 (16)–(18) 算出」的表來用。

讀 Table 1 的四個注意事項:價數是假設、量測溫度、Sn(w)、以及書中他處的引用我的判讀

1. 「價數」欄是自由電子模型的假設,不是量測結果。電子濃度 = 價數 × 原子濃度。例如 Al 的 18.06×10²² cm⁻³ = 3 × 6.02×10²²,意思是假設每個 Al 原子的三個價電子全部成為自由電子。對單價的鹼金屬,這個假設很好(p.155:Na、K 的 Hall 係數與每原子一個電子的計算值吻合極佳);對多價金屬則不然——p.156 指出三價的 Al、In,其 Hall 係數對應的是每原子一個正電荷載子,符號與大小都和「三個負電荷載子」不符,Be 也出現正號;這要到第 7–9 章的能帶論才能解釋。所以多價金屬那幾列,應讀成「若價電子全部自由,Fermi 參數會是多少」的名目值。

2. 為什麼要註明量測溫度。n 取決於原子間距,而原子間距隨熱膨脹改變,所以同一金屬在不同溫度下的 n、以及由它算出的 k_F、ε_F 會略有差異。表題註明 Na、K、Rb、Cs 用 5 K、Li 用 78 K 的數據,其餘為室溫;書中沒有說明為何這幾種金屬特別採用低溫數據。

3. Sn(w)。表中錫記作 Sn(w),書中沒有解釋 (w),全書文字層也搜不到 white tin 字樣。一般文獻以此表示白錫(金屬性的 β 相)——這是書外知識,僅供參考。

4. 書中他處引用本表時的一個小出入。第 9 章 p.239 估計鈉的內聚能時,說 Fermi 能是 3.1 eV、取自「Table 6.1」,但本表 Na 列是 3.23 eV,相差約 4%。書中沒有說明差異來源(本講解推測可能引用了不同溫度或不同來源的濃度),使用時以本表為準即可。同章 p.249 引用金的 k_F = 1.2×10⁸ cm⁻¹,則與本表一致。

無因次半徑參數 r_n(習題 8 記作 r_s)

r_n = r_0/a_H

p.139

r_0:含有一個電子的球之半徑;a_H:第一 Bohr 半徑(0.529 Å;p.159 習題 8 寫作 a_H = ħ²/e²m,CGS)。表中記作 r_n;p.159 習題 8 以完全相同的定義記作 r_s。

出處推導
Table 1 註腳 a 的定義。「(4π/3)r_0³ = V/N = 1/n」是把「含一個電子的球」翻成算式(每個電子分到的體積 1/n 等於該球的體積),註腳本身沒有寫出這一步。
極限行為
(i) r_n 小 = 電子密(Be 1.88);r_n 大 = 電子稀(Cs 5.63)。(ii) 以 r_s 表示時,每電子平均動能 ∝ 1/r_s²、庫侖能 ∝ 1/r_s(p.159 習題 8:2.21/r_s² 與 −1.80/r_s,單位 Ry),兩者之比 ∝ r_s:r_s 越大,交互作用相對越重要,「自由」電子的圖像越不自然。(iii) 習題 8 用這兩項求總能極小,得平衡值 r_s = 2.45,落在 Table 1 金屬的範圍(1.88–5.63)之內。
典型數值
Cu:r_0 = [3/(4π × 8.45×10²²)]^{1/3} = 1.414 Å,r_n = 1.414/0.529 = 2.67(表列 2.67)。

用 r_s 改寫 k_F、ε_F、v_F,並接上習題 8 的 2.21/r_s²補推導

r_s 形式的好處,是把所有 Fermi 參數化成「一個純數 ÷ r_s 的冪次」,材料之間的差別全部集中在 r_s 這一個數上。第 4 步用到的 a_H 與 Ry 的定義都取自 p.159 習題 8,沒有引入書外的定義。

  1. 第 1 步(r_0 與 n):「含一個電子的球」的體積等於每個電子分到的體積:(4π/3)r_0³ = 1/n,所以 n = 3/(4πr_0³)。
  2. 第 2 步(k_F r_0):代入 (16):k_F³ = 3π²n = 3π² · 3/(4πr_0³) = 9π/(4r_0³),所以 k_F r_0 = (9π/4)^{1/3} = 1.919,是一個與材料無關的純數。
  3. 第 3 步(k_F 以 r_s 表示):r_0 = r_s a_H,a_H = 0.5292 Å,所以 k_F = 1.919/(r_s a_H) = 3.63 Å⁻¹/r_s。Cu(r_s = 2.67):1.36 Å⁻¹ ✓。
  4. 第 4 步(ε_F 以 r_s 表示):ε_F = (ħ²/2m)k_F² = [ħ²/(2m a_H²)] · (9π/4)^{2/3}/r_s²。p.159 習題 8 給 a_H = ħ²/e²m、Ry = me⁴/2ħ²,代入得 ħ²/(2m a_H²) = ħ²·m²e⁴/(2m·ħ⁴) = me⁴/2ħ² = 1 Ry。又 (9π/4)^{2/3} = 3.683,所以 ε_F = 3.683 Ry/r_s² = 50.1 eV/r_s²(1 Ry = 13.61 eV)。
  5. 第 5 步(對回習題 8(a)):p.157 習題 1 要證明基態總動能 U₀ = (3/5)Nε_F(式 (60)),所以每電子平均動能 = (3/5) × 3.683 Ry/r_s² = 2.21 Ry/r_s²——正是 p.159 習題 8(a) 要你證明的結果。兩道習題透過本段的 (17) 串在一起。
  6. 第 6 步(v_F 以 r_s 表示):v_F = ħk_F/m = [ħ/(m a_H)] · 1.919/r_s。ħ/(m a_H) = e²/ħ = 2.19×10⁸ cm/s,所以 v_F = 4.20×10⁸ cm/s ÷ r_s。Cu:4.20/2.67 = 1.57×10⁸ cm/s ✓。
  7. 第 7 步(拿表的兩端驗算):Cs,r_n = 5.63:ε_F = 50.1/5.63² = 50.1/31.7 = 1.58 eV(表列 1.58);Be,r_n = 1.88:ε_F = 50.1/3.53 = 14.2 eV(表列 14.14)。

r_s 在說什麼:稀薄度,以及動能與庫侖能的拉鋸我的判讀

Table 1 多給了一欄 r_n,正文卻沒有說明它的用途;它的意義要從 p.159 的習題 8 讀出來。

它是一把無因次的「密度尺」。r_0 是每個電子分到的球半徑,a_H 是氫原子的特徵長度;兩者相除得到一個不帶單位的數,告訴你電子彼此相隔幾個 Bohr 半徑。Table 1 裡的金屬都落在約 1.9–5.6:Be 最密(1.88),Cs 最稀(5.63);鹼金屬群集在大的一端(3.25–5.63),貴金屬與其餘二至四價金屬多在約 2.1–3.0 之間(鹼土金屬 Ca、Sr、Ba 例外,約 3.3–3.7)。

為什麼用 Bohr 半徑當單位。一旦用 a_H 量長度、用 Rydberg 量能量,電子氣的兩個主要能量就都變成只含 r_s 的簡單式子(前一區塊):平均動能 2.21/r_s² Ry;庫侖能(習題 8 (b)、(c) 兩項相加)−1.80/r_s Ry。動能 ∝ 1/r_s² 只是 (17) 的 n^{2/3} 換個寫法;庫侖能 ∝ 1/r_s 則因為電荷間距 ∝ r_0。

拉鋸的意義。兩者之比 |庫侖能|/動能 = (1.80/r_s)/(2.21/r_s²) ≈ 0.81 r_s。代入 Cu 的 r_s = 2.67,比值約 2.2——庫侖能其實比動能還大。換句話說,把傳導電子當成「自由」的,並不是因為交互作用很小。p.133 對此的說法很節制:自由電子模型最有用的,是那些主要取決於傳導電子運動(動理)性質(kinetic properties)的物理量;p.134 並指出傳導電子很少被其他傳導電子散射,是 Pauli 原理的後果。本段的 Fermi 球圖像正好解釋了後者:只有 Fermi 面附近的電子有空位可以散射進去。

記號提醒。同一個量在 Table 1 記作 r_n、在習題 8 記作 r_s,定義完全相同;本講解在推導中沿用 r_s。

Fermi 溫度

T_F ≡ ε_F/k_B

p.140

k_B:Boltzmann 常數;T_F:以溫度單位表示的 Fermi 能。此式在 p.140 行文中定義,無式號;Table 1 欄標題亦寫 T_F ≡ ε_F/k_B。

出處推導
純定義,沒有推導。p.144 式 (35) 以 k_BT_F ≡ ε_F 的形式再次使用。
極限行為
(i) T ≪ T_F:簡併電子氣(所有金屬在室溫),熱只能激發 ε_F 附近 ~k_BT 內的電子,比例 ~T/T_F。(ii) T ≳ T_F:古典 Boltzmann 氣體(p.143 Fig. 7 標為古典氣體區);金屬在熔化前到不了這一區,稀薄的半導體載子則室溫就在這一區。(iii) 熱容:p.144 式 (36) C_el = (1/2)π²Nk_BT/T_F,T_F 直接成為電子熱容的溫度尺度。
典型數值
Cu:7.00 eV × 11,604 K/eV = 8.12×10⁴ K(表列 8.12×10⁴ K);300 K 時 T/T_F = 3.7×10⁻³,ε_F/k_BT ≈ 270。

書中的警告:T_F 與電子氣的溫度無關(p.140、p.145)

在介紹 Table 1 的同一段裡,Kittel 定義 T_F ≡ ε_F/k_B,並立刻加上一句括號警告:「The quantity T_F has nothing to do with the temperature of the electron gas!」(p.140)。到了 p.145 他又提醒一次:T_F 雖然叫做 Fermi 溫度,卻不是電子的溫度,只是一個方便的參考記號。

T_F 其實在本章更早就出現過:p.136 的 Fig. 3 畫 Fermi-Dirac 分布時,就是取 T_F = 50,000 K,把 500 K、5,000 K、1×10⁴ K、2.5×10⁴ K、5×10⁴ K、10×10⁴ K 各溫度的曲線畫在一起;圖說並說明化學勢可由 f = 0.5 的位置讀出。

到了 p.141–142,T_F 的角色就明朗了:升溫時只有約 T/T_F 比例的電子能被熱激發,所以電子熱容比古典值 (3/2)Nk_B 小了一個約 0.01 或更小的因子(p.142 以 T_F ~ 5×10⁴ K 估計)。

T_F 為什麼高達數萬 K,以及它真正在量什麼我的判讀

高,是因為 Pauli,不是因為熱。要把 N 個電子放進 k 空間,每個軌域只能放一個,就必須一路填到 ε_F;ε_F 的大小完全由濃度決定((17)),與溫度無關。金屬的電子濃度約 10²² cm⁻³,塞進這麼多費米子所需的「最高能量」就是幾個 eV——換成溫度單位便是幾萬 K。Table 1 的 T_F 從 1.83×10⁴ K(Cs)到 16.41×10⁴ K(Be),比任何金屬能以固態存在的溫度都高得多。

T/T_F 是「量子程度」的刻度。T ≪ T_F 的電子氣稱為簡併的(degenerate):分布接近 T = 0 的階梯,只在 ε_F 附近 k_BT 的範圍內被抹開。T ≳ T_F 時才回到古典的 Boltzmann 氣體;p.143 的 Fig. 7 把這兩個區域分別標為簡併量子氣體區與古典氣體區。室溫下金屬的 T/T_F 只有 0.002–0.016(Table 1 的兩端),電子氣永遠處在深度簡併區。反之,輕摻雜半導體的載子濃度低了 6 個數量級以上,T_F 只有幾 K(n = 10¹⁶ cm⁻³ 時約 2 K,自由電子質量估計),室溫下是古典氣體——同一套公式,兩種完全相反的極限。

用 Fig. 3 看 T/T_F 的效果(p.136,目視判讀)。該圖以 T_F = 5×10⁴ K 作圖:500 K(T/T_F = 0.01)的曲線在 ε_F 處幾乎是垂直的階梯;5,000 K(0.1)已看得出明顯的圓角;到 T = T_F = 5×10⁴ K 時,連 ε = 0 處的佔據機率都只剩約 0.5,也就是化學勢已降到能量原點附近——分布完全不像階梯了。大多數金屬在室溫落在比 500 K 那條曲線更陡的位置(Cu 的 T/T_F ≈ 0.004);最稀的 K、Rb、Cs(0.012–0.016)則略緩於它,但仍遠比 5,000 K 那條陡。

與聲子對比。第 5 章的 Debye 溫度多在 10²–10³ K,室溫下聲子已接近古典極限;電子的 T_F 卻是數萬 K,室溫下是深度量子氣體。這個對比決定了 p.145 的低溫熱容形式 C = γT + AT³:在同時遠低於 θ 與 T_F 的溫度下,電子項線性、聲子項立方,溫度夠低時線性的電子項會勝出。

態密度的定義與註腳 1(p.140–141)

在 Table 1 之後,Kittel 轉向本段的最後一個量:每單位能量範圍內的軌域數 D(ε),稱為態密度(density of states)。他在這個詞上加了註腳 1:嚴格說來,D(ε) 是單粒子態的密度,也就是軌域密度(density of orbitals)。Fig. 5 的縱軸也正是這樣標的。

做法是:用 (17) 求出能量 ≤ ε 的軌域總數,得 (19);對 ε 微分,得態密度 (20)(書中以方框標示)。p.141 接著比較 (19) 與 (20),在能量 ε 處得到更簡潔的 (21),並評論:在 Fermi 能處,每單位能量的軌域數在一個量級為 1 的因子之內,就等於傳導電子總數除以 Fermi 能——正如我們所預期的。

「軌域」一詞的意義要回到 p.134–135:軌域是只含一個電子的系統之波動方程的解;Kittel 用它來區分「N 個交互作用電子系統的精確量子態」與「把 N 個電子分派到 N 個不同軌域所構成的近似態」。在一維的例子中,軌域以 n 與 m_s 兩個量子數標記,每個軌域至多容納一個電子。

能量 ≤ ε 的軌域總數

N = (V/3π²)(2mε/ħ²)^{3/2}

(19)  p.140

此處的 N 應讀作 N(ε):能量不超過 ε 的軌域總數(兩個自旋都算、對整個體積 V);ε:任意能量。

出處推導
原文說用 (17) 求能量 ≤ ε 的軌域總數。做法是把 (17) 中的 ε_F 換成任意 ε、再解出 N:若把能量 ≤ ε 的軌域全部填滿,所需的電子數就是讓 ε 恰好成為「Fermi 能」的那個 N(逐步推導見下一區塊)。
極限行為
(i) ε → 0:N ∝ ε^{3/2} → 0。(ii) ε = ε_F:N(ε_F) 等於實際的電子總數。(iii) (19) 是「盒子」的性質,與電子實際填到哪裡無關——它數的是存在的軌域,不是被佔據的軌域;T > 0 時佔據數要再乘上 Fermi-Dirac 分布 f。(iv) N(ε) ∝ V·m^{3/2}·ε^{3/2}。
典型數值
V = 1 cm³、ε = 7.00 eV:2mε/ħ² = 7.00/3.810 = 1.837 Å⁻²,其 3/2 次方 = 2.490 Å⁻³,除以 3π² = 29.61 得 0.0841 Å⁻³ = 8.41×10²² cm⁻³——回到 Cu 的 n(差 0.5% 來自表列 ε_F 的捨入)。

(17) → (19) → (20):兩條路,同一個 √ε補推導

書上從 (17) 到 (19) 到 (20),每一步都只用一句話;這裡把兩個「顯然」的步驟補齊,並給出不依賴 (17) 的第二條路。

第二條路更值得記:它把態密度拆成「k 空間的幾何(k²)」乘上「色散關係(dk/dε)」,換任何色散關係都能直接套用——例如把 m 換成有效質量、把能量原點移到導帶底,就得到第 8 章 p.205 的式 (37)。√ε 的來源也在第 6 步看得一清二楚:k² 項給出 ε 的一次方,dk/dε 給出 ε^{−1/2},相乘剩下 ε^{1/2}。

  1. 第 1 步(重新詮釋 (17)):(17) 說:N 個電子填滿時,最高能量是 ε_F = (ħ²/2m)(3π²N/V)^{2/3}。反過來讀:對任意能量 ε,能量 ≤ ε 的軌域數 N(ε) 就是滿足 ε = (ħ²/2m)(3π²N(ε)/V)^{2/3} 的那個數。
  2. 第 2 步(解出 N(ε)):兩邊乘 2m/ħ²:2mε/ħ² = (3π²N/V)^{2/3}。兩邊取 3/2 次方:(2mε/ħ²)^{3/2} = 3π²N/V。所以 N(ε) = (V/3π²)(2mε/ħ²)^{3/2},即 (19)。
  3. 第 3 步(微分):把 (19) 寫成 N = (V/3π²)(2m/ħ²)^{3/2}·ε^{3/2},對 ε 微分:dN/dε = (V/3π²)(2m/ħ²)^{3/2}·(3/2)ε^{1/2}。
  4. 第 4 步(化簡係數):(1/3π²)·(3/2) = 1/(2π²)。所以 D(ε) = dN/dε = (V/2π²)(2m/ħ²)^{3/2}ε^{1/2},即 (20)。
  5. 第 5 步(第二條路:直接數 k 空間球殼):半徑 k 到 k + dk 的球殼體積 4πk²dk,含 k 點 4πk²dk/(8π³/V) = Vk²dk/2π² 個,乘上自旋 2 得 Vk²dk/π² 個軌域。
  6. 第 6 步(把 dk 換成 dε):由 ε = ħ²k²/2m,k = (2mε)^{1/2}/ħ,dk/dε = m/(ħ²k) = (1/2)(2m/ħ²)^{1/2}ε^{−1/2}。所以 D(ε) = (V/π²)k²·(dk/dε) = (V/π²)(2mε/ħ²)·(1/2)(2m/ħ²)^{1/2}ε^{−1/2} = (V/2π²)(2m/ħ²)^{3/2}ε^{1/2}——與第 4 步相同。
  7. 第 7 步(對照第 5 章):第 5 章 p.112 式 (20) 的聲子一般式是 D(ω) = (VK²/2π²)(dK/dω)(每個偏振)。把 ω 換成 ε、再乘上電子的 2 個自旋,就是第 6 步的 D(ε) = (Vk²/π²)(dk/dε)。計數公式完全相同,只有色散關係與內部自由度的數目不同。

三維自由電子的態密度

D(ε) ≡ dN/dε = (V/2π²)·(2m/ħ²)^{3/2}·ε^{1/2}

(20)  p.140

D(ε):每單位能量範圍內的軌域數,對整個樣品體積 V、兩個自旋都算在內(註腳 1:嚴格說是單粒子態密度,即軌域密度)。書中以方框標示此式。

出處推導
對 (19) 微分;或直接數 k 空間球殼(兩條路見上一區塊)。
極限行為
(i) ε → 0 以 √ε 趨近零:帶底附近的態很少。(ii) D 正比於 V(外延量);除以 V 得每單位體積的態密度——第 8 章 p.205 式 (37) 正是「由 (6.20)」寫出單位體積版本 (1/2π²)(2m_e/ħ²)^{3/2}(ε − E_c)^{1/2},把 m 換成有效質量、能量原點移到導帶底。(iii) 維度:二維為常數(p.157 習題 3 附註:每單位面積 m/πħ²),一維 ∝ ε^{−1/2}(由式 (4) 可推,見 (21) 前面的推導區塊「(19)(20) → (21)」)。(iv) 有限尺寸:第 18 章 p.520 式 (3) 指出有限系統的態密度是一串 δ 函數,(20) 是其平滑化的巨觀極限。(v) D 的斜率為正,因此三維的化學勢隨溫度略降(p.144 Fig. 8);斜率為負的一維則相反。
典型數值
Cu 在 ε_F:每單位體積 (1/2π²)(2m/ħ²)^{3/2}ε_F^{1/2} ≈ 1.8×10²² eV⁻¹cm⁻³;等價地 3n/2ε_F = 3 × 8.45×10²²/(2 × 7.00) = 1.81×10²² eV⁻¹cm⁻³。換成每個原子為 0.214 eV⁻¹(Cu 單價,每原子一個傳導電子;兩個自旋都算)。

D(ε) 的計數約定:整個 V、兩個自旋——以及全書其他地方怎麼換約定我的判讀

「D(ε) 對整個體積 V、已含自旋因子」這句話看似平凡,卻是全書最常見的因子 2 與因子 V 錯誤的源頭。

本段的約定。D(ε) 數的是軌域,而在 Kittel 的用法中,一個軌域是「一個空間態加一個自旋值」(p.135:每個 n 有自旋向上、向下兩個軌域;p.138:(15) 所說的「軌域總數」已經乘上自旋因子 2)。所以 (19)、(20) 的 N(ε) 與 D(ε) 兩個自旋都算在內,而且對整個樣品體積 V——D(ε) 的單位是「個/能量」,是外延量。註腳特別強調 density of orbitals,用意是提醒:這裡數的是單電子軌域,不是 N 電子系統的多體量子態(p.134 所區分的那兩種「態」)。

全書其他地方換了約定:

出處寫法約定
本段 p.140 式 (20)D(ε) = (V/2π²)(2m/ħ²)^{3/2}ε^{1/2}整個 V、兩個自旋
p.144 式 (34)、(35)D(ε_F) = 3N/2ε_F同上(書中寫明由 (21) 得來)
第 8 章 p.205 式 (37)D_e(ε) = (1/2π²)(2m_e/ħ²)^{3/2}(ε − E_c)^{1/2}每單位體積、兩個自旋;有效質量;能量原點在導帶底(書中寫明由 (6.20) 得來)
第 10 章 p.278BCS 的 T_c 判據明寫用「一個自旋」的 D(ε_F)
第 11 章 p.317 式 (43)½D(ε + μB)明寫這是「一個自旋方向」的軌域密度;D(ε_F) 本身仍引用第 6 章的 3N/2ε_F
第 18 章 p.520 式 (3)一串 δ 函數之和有限系統的離散態密度

實用守則:每次看到 D,先問三件事——(a) 是對整個 V、每單位體積,還是每個原子?(b) 兩個自旋都算,還是只算一個?(c) 是每單位能量,還是每單位頻率(第 5 章聲子的 D(ω))?三者任何一個弄錯,就會差一個因子 2、一個 V 或一個 ħ。

(19)(20) → (21):3/2 從哪裡來,以及任意維度的一般式補推導

(21) 看起來只是 (19)、(20) 的比值,但它把態密度從一個含 V、m、ħ 的公式,變成只需要「電子數」與「能量」的無因次關係——這就是它在下一段熱容計算中特別好用的原因(p.144 直接寫「由 (21)」代入)。

第 6、7 步順便說明了:一維、二維、三維態密度的差別,完全來自 k 空間的維度。斜率的正負(一維遞減、二維平坦、三維遞增)正是 p.144 Fig. 8 中一維與三維化學勢往相反方向走的根源(二維的情形見習題 3)。

  1. 第 1 步:把 (19) 寫成 N(ε) = Cε^{3/2},其中 C = (V/3π²)(2m/ħ²)^{3/2} 與 ε 無關。
  2. 第 2 步:(20) 就是 D(ε) = dN/dε = (3/2)Cε^{1/2}。
  3. 第 3 步:相除:D/N = (3/2)Cε^{1/2}/(Cε^{3/2}) = 3/(2ε),所以 D(ε) = 3N(ε)/2ε,即 (21)。所有材料常數(V、m、ħ)都消掉了。
  4. 第 4 步(在 Fermi 能取值):N(ε_F) 就是電子總數 N,所以 D(ε_F) = 3N/2ε_F;再用 ε_F = k_BT_F 得 D(ε_F) = 3N/2k_BT_F——這就是 p.144 的式 (35)。
  5. 第 5 步(3/2 的來源):取對數微分,d ln N/d ln ε = 3/2。N ∝ k³(三維球體積),k ∝ ε^{1/2}(拋物線色散),所以 N ∝ ε^{3/2}。一般地,d 維中 N ∝ k^d ∝ ε^{d/2},D = (d/2)·N/ε。
  6. 第 6 步(一維驗算):由 p.135 式 (4) 解出 N = (2L/πħ)(2mε)^{1/2},微分得 D₁ = (L/πħ)(2m/ε)^{1/2} = N/2ε,符合 d = 1。一維態密度 ∝ ε^{−1/2},隨能量遞減。
  7. 第 7 步(二維驗算):二維 N = 2·πk²/(2π/L)² = Ak²/2π(A = L²),代入 k² = 2mε/ħ² 得 N = Amε/πħ²,所以 D₂ = Am/πħ² = N/ε,是常數——與 p.157 習題 3 附註「二維自由電子氣的軌域密度與能量無關、每單位面積為 m/πħ²」一致。

態密度與累積軌域數的關係

D(ε) ≡ dN/dε = 3N/2ε

(21)  p.141

此處的 N 是 (19) 的 N(ε),不是固定的電子總數;只有在 ε = ε_F 時才等於電子總數。原文寫的是比較 (19) 與 (20),在能量 ε 處得到此式。

出處推導
(20) 除以 (19):D/N = (3/2)/ε(逐步推導見上一區塊)。
極限行為
(i) ε = ε_F:D(ε_F) = 3N/2ε_F = 3N/2k_BT_F——p.144 式 (35) 寫明由 (21) 得到,並代入熱容式 (34)。(ii) 原文的物理解讀(p.141):在 Fermi 能處,每單位能量的軌域數在一個量級為 1 的因子之內,等於傳導電子總數除以 Fermi 能。那個因子就是 3/2:N/ε_F 是 [0, ε_F] 區間的平均態密度,而 D ∝ √ε 遞增,端點值比平均值大 3/2 倍。(iii) 倒數 1/D(ε_F) = 2ε_F/3N 是 Fermi 能附近相鄰能階的平均間距——第 18 章 p.547 式 (55) 正是用 (6.21) 寫出 Δε = 1/D(ε_F) = 2ε_F/3N。(iv) d 維推廣 D = (d/2)N/ε,3/2 就是 d/2。
典型數值
Cu:D(ε_F)/N = 3/(2 × 7.00 eV) = 0.214 eV⁻¹(每個電子,也就是每個 Cu 原子;兩個自旋都算)。Au 奈米粒子 R = 2 nm:N = 5.90×10²² cm⁻³ × (4π/3)(2×10⁻⁷ cm)³ ≈ 1.98×10³,Δε = 2 × 5.51 eV/(3 × 1.98×10³) ≈ 1.9 meV,與 p.547 所說「約 2 meV」一致。

Fig. 5(p.140)文字重建:D(ε) 與有限溫度下的佔據我的判讀

Fig. 5 是本段唯一的函數圖,也是下一段的論證基礎。網頁不放圖,以下依原書算繪影像逐項還原(形狀細節屬目視判讀)。

座標軸。橫軸是能量 ε(軸下標注 Energy, ε,附向右箭頭),軸上只標了一個刻度 ε_F;縱軸是 D(ε),軸旁文字為「Density of orbitals, relative scale」——縱軸只是相對尺度,沒有數值。

實線:D(ε)。從原點出發,起初幾乎垂直地竄升,然後逐漸變緩、向右上方延伸,越過 ε_F 之後仍繼續上升——這就是 (20) 的 √ε 形狀(原點處切線垂直,因為 dD/dε ∝ ε^{−1/2} 在 ε → 0 時發散)。

灰色陰影:T = 0 的佔據。從 0 到 ε_F、實線以下的區域全部塗灰,右邊界是位於 ε_F 的一條垂直線。圖說寫明陰影面積代表絕對零度時被填滿的軌域。因為 T = 0 時 f 是階梯函數,被佔據的態密度 f·D 就是「D 截到 ε_F 為止」。

虛線:有限溫度的 f(ε,T)D(ε)。圖說指明虛線是有限溫度下被填滿軌域的密度,而且溫度低到 k_BT 遠小於 ε_F。低能處虛線與實線重合;接近 ε_F 時虛線先微微拱起、隨即向下彎離實線,在 ε = ε_F 處大約落在 D(ε_F) 的一半高度(與 f(μ) = ½、μ ≈ ε_F 一致),越過 ε_F 後拖著一條迅速衰減的尾巴,貼近橫軸向右延伸。

區域 1 與區域 2。「1」標在 ε_F 左側、實線與虛線之間的灰色楔形內:這些軌域在 T = 0 時是滿的,升溫後部分被騰空。「2」標在 ε_F 右側、虛線之下的白色區域內:這些軌域在 T = 0 時是空的,升溫後部分被填入。圖說的結論是:溫度由 0 升到 T 時平均能量增加,因為電子被熱激發由區域 1 跑到區域 2。

k_BT 標尺。圖的上方有一支左右雙向箭頭,標為 k_BT,兩端各拉一條垂直細線,大致跨在 ε_F 兩側——表示被「抹開」的能量範圍只有 k_BT 的量級。(這是示意圖:虛線開始偏離實線的位置,畫得比標尺左緣還早一些,不必拿尺去量。)

由區域 1 激發到區域 2:下一段熱容論證的起點我的判讀

Fig. 5 的圖說只有兩句,卻已經包含電子熱容問題的全部答案;p.141 起的下一段就是把它翻成算式。這裡先把邏輯走一遍,讓你讀下一段時知道每一步對應圖上的哪裡。

古典預期錯在哪。p.141 指出,古典統計力學給每個自由粒子 (3/2)k_B 的熱容,N 個電子就是 (3/2)Nk_B;但室溫實測的電子貢獻通常不到這個值的 1%。古典圖像默認每個電子升溫時都能多拿到 ~k_BT 的能量。

Pauli 原理怎麼改寫它。位於 Fermi 海深處(ε_F − ε ≫ k_BT)的電子,吸收 ~k_BT 的能量後想去的軌域都已經滿了——它無處可去。只有離 ε_F 在 k_BT 以內的電子,才找得到 k_BT 範圍內的空軌域。這正是圖中虛線只在 ε_F 附近、寬約 k_BT 的範圍內偏離階梯的原因。

數量級估計。被激發的電子數 ≈ 區域 2 的面積 ≈ D(ε_F)·k_BT。用 (21),D(ε_F) = 3N/2ε_F,所以被激發的電子數 ≈ (3N/2ε_F)k_BT = (3/2)N(T/T_F)——只是總數的 T/T_F 量級。每個被激發的電子多拿到 ~k_BT,總激發能 U ~ N(T/T_F)k_BT,對 T 微分得 C ~ Nk_B(T/T_F)。這就是 p.141 的 (22)、(23);精確係數要到 p.144 才算出:式 (34) C_el = (1/3)π²D(ε_F)k_B²T,代入 (35) D(ε_F) = 3N/2k_BT_F,得 (36) C_el = (1/2)π²Nk_BT/T_F。從頭到尾,唯一需要的電子結構資訊就是 D(ε_F)——本段 (21) 的用處就在這裡。

代入 Cu 的數字。室溫 T = 300 K,T/T_F = 300/81,200 ≈ 3.7×10⁻³。被激發的電子約佔 (3/2)(T/T_F) ≈ 0.55%。C_el 與古典值 (3/2)Nk_B 之比為 (π²/3)(T/T_F) ≈ 0.012——正是 p.141–142 所說「約 0.01 或更小」的量級。

容易被忽略的細節:區域 1 與區域 2 的面積必須相等。f·D 曲線下的總面積就是電子數 N,而 N 不隨溫度改變(p.136 Fig. 3 圖說也強調粒子總數與溫度無關),所以從區域 1 騰出的電子數必須等於進入區域 2 的電子數。但三維的 D(ε) 是遞增的:若 μ 固定在 ε_F,f 對 ε_F 的對稱抹開會讓上方(態較多)得到的電子多於下方失去的。為了守住 N,化學勢必須略往下移——這正是 p.144 Fig. 8 中三維曲線 μ/ε_F 隨 τ/ε_F 下降的原因;一維的態密度遞減,所以一維曲線反而上升。常見金屬的 τ/ε_F 約 0.01,這個偏移非常小,所以 p.144 可以放心把 μ 換成 ε_F。

與 p.142 式 (27) 的對應。p.142 把能量增量 ΔU 拆成兩個積分:第一個是把電子從 ε_F 送到 ε > ε_F 的軌域所需的能量(對應區域 2),第二個是把電子從 ε_F 以下的軌域帶到 ε_F 所需的能量(對應區域 1 被騰空的部分);兩者都為正。Fig. 5 的兩塊區域就是這兩個積分的幾何圖像。

本段最容易犯的六個錯誤我的判讀

1. 「T_F 是電子的溫度,所以金屬裡的電子很熱。」

錯在把能量尺度當成溫度。T_F 只是把 ε_F 以溫度單位表示,p.140 與 p.145 兩度強調它與電子氣的溫度無關。平衡時電子與晶格同溫;電子「能量高」是因為 Pauli 原理迫使它們往上填,不是因為熱。正確理解:T_F 是判斷簡併程度的標尺,T ≪ T_F 時電子氣是高度簡併的量子氣體。

2. 「ε_F 是電子的平均能量。」

錯。ε_F 是最高填滿軌域的能量;基態每電子平均能量是 (3/5)ε_F(p.157 習題 1,式 (60):U₀ = (3/5)Nε_F)。快速驗證:平均能量 = ∫₀^{ε_F} εD dε ÷ ∫₀^{ε_F} D dε,D ∝ ε^{1/2},分子 ∝ (2/5)ε_F^{5/2},分母 ∝ (2/3)ε_F^{3/2},比值 (3/5)ε_F。因為 D 隨 ε 增加、電子偏重高能端,平均值比 ε_F/2 大。

3. 「(21) 說 D = 3N/2ε,所以態密度隨能量遞減。」

錯在把 (21) 中的 N 當成固定的電子總數。(21) 是 (20) 與 (19) 的比,其中 N 是 N(ε)——能量 ≤ ε 的軌域數,本身 ∝ ε^{3/2}。代回去 D ∝ ε^{3/2}/ε = ε^{1/2},是遞增的。只有在 ε = ε_F 時,N(ε_F) 才等於電子總數,得到 p.144 的 D(ε_F) = 3N/2ε_F。

4. 「D(ε) 的因子 2 與體積因子全書一致,公式可以直接搬。」

不一致。本段的 D 是整個 V、兩個自旋;第 8 章 p.205 是每單位體積;第 10、11 章明寫改用一個自旋(見前面的對照表)。直接搬公式,最常錯一個 2 或一個 V。

5. 「駐波與週期性邊界條件給出不同的物理,其中一個是錯的。」

兩者對巨觀樣品給出完全相同的態數(Vk³/6π²),差別只是 1/(k_F L) 量級的表面項。它們真正的差別在態的性質:駐波不帶動量,行進波帶 ±ħk。計數時可任選;討論電流時用行進波;處理奈米尺度的受限方向時,必須回到真實邊界(第 18 章)。

6. 「v_F 是熱運動速度;降到 0 K 電子就停下來。」

錯。v_F 由濃度決定((18)),與溫度無關,0 K 時依然存在——它是 Pauli 原理強加的「零點運動」。Cu 的 v_F = 1.57×10⁸ cm/s,而 300 K 的古典熱速度 (3k_BT/m)^{1/2} 只有約 1.2×10⁷ cm/s,小了一個數量級以上。後面 p.156 電子熱導率所用的平均自由徑 ℓ = v_Fτ,速度也是 v_F 而非熱速度。

接到 Cth:這一段的量在元件熱分析裡出現在哪裡延伸知識

(本區塊為延伸內容:Kittel 書中沒有元件熱阻或 FinFET 的討論。以下只做數量級連結,引用的書內公式都標明頁碼。)

1. 電子對熱容的貢獻可以略去——Cth 裡的「C」幾乎全是晶格。本段給出的 T_F ~ 10⁴–10⁵ K,意味著室溫下金屬電子的熱容只有古典值的 ~1%。以 Cu 為例,用 p.144 式 (36) 估計:C_el = (π²/2)Nk_B(T/T_F) ≈ 4.93 × 3.7×10⁻³ Nk_B ≈ 0.018Nk_B,而室溫下晶格熱容接近古典的 3Nk_B(第 5 章 Debye 模型的高溫極限),兩者之比約 0.6%。所以在萃取元件或內連線的熱容網路時,金屬中的電子熱容可以安全略去,熱容由聲子決定——那是第 5 章的領域。

2. 但金屬裡的熱主要是電子在搬——v_F 直接進入熱導。p.156 式 (56) 給出電子熱導率 K_el = π²nk_B²Tτ/3m,推導中用了 ℓ = v_Fτ 與 ε_F = ½mv_F²;Wiedemann–Franz 定律(p.156–157)再把它與電導率綁在一起。FinFET 的接觸、via 與 BEOL 銅導線,導熱主要靠傳導電子。當線寬縮小到與平均自由徑可比時,表面散射把 τ 截短——p.159 習題 10 正是用 τ ≈ d/v_F 估計薄膜的最大面電阻;同樣的截短也會壓低 K_el。v_F 就是這個尺寸效應的速度尺度。

3. 通道中的載子通常不簡併。把本講解的實用式 ε_F = 1.693 eV × (n/10²² cm⁻³)^{2/3}(自由電子質量、單一能谷)用到半導體:n = 10¹⁸ cm⁻³ 時 ε_F ≈ 3.6 meV、T_F ≈ 40 K ≪ 300 K,電子氣是非簡併的,應使用 p.137 的高能尾 Boltzmann 極限(第 8 章 p.205 式 (35) 也是這樣處理);n ~ 10²⁰ cm⁻³ 的重摻雜源/汲極則 ε_F ≈ 0.08 eV、T_F ≈ 900 K,屬於簡併區。矽真實的有效質量與多能谷會改變這些數字(書外知識),但「通道非簡併、重摻雜區簡併」的分界不會變。無論哪一種,載子熱容都遠小於晶格熱容,Cth 中的熱容仍以矽晶格為主。

4. 奈米尺度的邊界條件。fin 寬度只有數奈米時,受限方向的量子化必須用盒中駐波 (7) 處理(第 18 章 p.528 對奈米線的做法),而不是週期性條件 (9)。這主要影響載子能階與遷移率;對 Cth 的影響是間接的(改變焦耳熱在何處、以何種速率經電子–聲子耦合交給晶格),已超出本段範圍。

誠實的結論:本段對 Cth 的直接貢獻是兩個數量級判斷——電子熱容可忽略(T ≪ T_F)、金屬導熱的速度尺度是 v_F。更定量的部分(熱導率、Wiedemann–Franz、尺寸效應)在本章 p.156–159。

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

往前接:

  • 第 5 章 p.110–112:週期性邊界條件、(2π/L)³ 的計數(p.111 式 (18)(19))、以球殼體積換態密度(p.112 式 (20))——本段把整套工具從聲子搬到電子(p.137 原文就說與第 5 章處理聲子的做法相同)。
  • 本章 p.134–135:一維盒中電子、軌域的定義、Pauli 填法、式 (4)。
  • 本章 p.136–137:Fermi-Dirac 分布 (5) 與 T = 0 時 μ = ε_F。

本章稍後直接使用:

  • p.141–144 電子熱容:Fig. 5 的區域 1 → 2 論證;(34) 需要 D(ε_F);(35) 直接引用 (21);(36) 以 T_F 表示。
  • p.147–148 電導與 Ohm 定律:Fermi 球在電場中整體平移(需要平面波的動量本徵值 ħk)。
  • p.156–157 熱導與 Wiedemann–Franz 定律:用到 ε_F = ½mv_F² 與 ℓ = v_Fτ。
  • 習題(p.157–159):習題 1 的 U₀ = (3/5)Nε_F;習題 2 以 Table 1 估計鉀的電子氣體積模數;習題 3 的二維態密度 m/πħ²;習題 4 的相對論 Fermi 氣體(白矮星、脈衝星);習題 5 液態 He³ 的 ε_F 與 T_F;習題 8 以 r_s 表示平均動能 2.21/r_s² Ry。

往後被引用:

  • 第 7 章 p.164:近自由電子模型的出發點式 (1)–(3),就是本段的 (12)、(10)、(9)。
  • 第 8 章 p.205 式 (37):半導體導帶態密度寫明由 (6.20) 得來,改成每單位體積與有效質量。
  • 第 9 章 p.239:鈉的內聚能估計引用「Table 6.1」的 Fermi 能,並用平均動能 = 0.6 ε_F(該處寫 3.1 eV,與本表 Na 列 3.23 eV 有約 4% 的出入)。
  • 第 9 章 p.249:銅的自由電子 Fermi 球直徑 9.80/a 與布里淵區比較;金的 k_F 引用「Table 6.1」算極值軌道面積。
  • 第 10 章 p.278:BCS 的 T_c 判據使用單一自旋的 D(ε_F)。
  • 第 11 章 p.317:Pauli 順磁性,D(ε_F) = 3N/2ε_F 寫明取自第 6 章,並以 ½D 表示單一自旋。
  • 第 14 章 p.404–405:Thomas–Fermi 近似以 (6.17) 寫局域 Fermi 能,並用自由電子的 D(ε_F);書中給 Cu(n₀ = 8.5×10²² cm⁻³)的屏蔽長度為 0.55 Å。
  • 第 18 章 p.520、p.528、p.547:有限系統的 δ 函數態密度、奈米線截面的盒中粒子解,以及以 (6.21) 估計金屬奈米粒子在 Fermi 能處的平均能階間距(Au、R = 2 nm 時約 2 meV)。