問題的提出:古典理論預測的 ³⁄₂Nk_B 到哪裡去了
本段一開頭(p.141),Kittel 就把它定位為早期金屬電子論裡最棘手的一個問題:傳導電子的熱容。
古典統計力學的結論很直接:一個自由粒子的熱容是 ³⁄₂k_B。如果金屬裡 N 個原子各把一個價電子交給電子氣,而且這些電子可以自由移動,那麼電子部分的熱容就應該是 ³⁄₂Nk_B,和單原子氣體的原子一模一樣。(本站補充:³⁄₂k_B 來自三個平移自由度,依能量均分定理平均能量是 ³⁄₂k_BT;而「電子可以自由移動」正是用電子氣解釋導電所需要的前提。)
實驗卻不支持。Kittel 指出,室溫下觀測到的電子熱容通常不到這個值的 0.01。這個落差困擾了 Lorentz 等早期研究者:電子在導電過程中表現得像是可以自由移動,為什麼在吸收熱量時卻幾乎不參與?同一群電子,在電導實驗裡是自由的,在量熱實驗裡卻好像被凍住了。
Kittel 說這個問題要等到包立不相容原理與費米分布函數被發現之後才得到解答,並引述 Fermi 本人的結論:比熱在絕對零度時消失(原文 “the specific heat vanishes at absolute zero”),而在低溫下正比於絕對溫度。這兩個特徵——T → 0 時趨於零、低溫時線性於 T——正是本段接下來要從定性(式 (22)–(23))與定量(式 (24)–(36))兩個層次推出來的結果,也是下一小節「金屬熱容的實驗」要檢驗的對象。
為什麼這個矛盾這麼要命:和第 5 章的晶格熱容擺在一起看我的判讀
把這個矛盾放回第 5 章的脈絡,更能看出它的份量。第 5 章 p.117 把晶格熱容的高溫極限 3Nk_B 稱為 Dulong–Petit 值,每莫耳約 3R ≈ 24.9 J/mol-K。如果每個一價金屬原子還額外帶著一個古典自由電子,高溫熱容就應該是 3R + ³⁄₂R = 4.5R ≈ 37 J/mol-K,比 Dulong–Petit 值多出一半。這不是細微修正,而是量熱實驗一眼就能看到的差異——然而實際上看不到。
換句話說,古典電子氣模型在同一個框架內自相矛盾:它用「電子自由移動」來解釋歐姆定律(本章下一段,p.147 起),卻因此預測了一個不存在的熱容。解法不是修補古典模型,而是承認電子是費米子、必須服從包立原理。這是本書中「費米子的量子統計在巨觀可量測量上造成數量級效應」的第一個例子。量子統計本身並不是第一次登場:第 5 章已用 Planck 分布(p.107)處理聲子熱容,Einstein 模型也解釋了鑽石熱容在低溫被壓低(p.117 Fig. 11);這裡輪到的是費米–狄拉克統計與包立原理。
本段的論證方式之後會原封不動地再出現。第 11 章 p.316 的包立自旋磁化率用的是同一套邏輯:古典理論預測的順磁磁化率被壓低,因為只有費米面附近 k_BT 範圍內的電子能翻轉自旋,比例只有 T/T_F。把本段讀懂,那一段幾乎不用再學。
定性解答:只有費米面附近約 k_BT 寬的電子能被加熱(Fig. 5)
Kittel 的定性解答只用了一句話的物理:把樣品從絕對零度加熱時,並不是每個電子都像古典預期那樣多得到 ~k_BT 的能量,只有落在費米能階附近 k_BT 範圍內的軌域上的電子,才會被熱激發。
本站補充:為什麼只有這些電子?關鍵在包立不相容原理。T = 0 時電子從最低能量開始,每個軌域一個電子,填到 ε_F 為止。現在給系統 ~k_BT 的熱能。一個深埋在 ε_F 以下、距離費米面好幾個 k_BT 的電子,想吸收 ~k_BT 往上跳,但它上方 k_BT 範圍內的軌域全都被佔滿了——它無處可去。只有離費米面在 ~k_BT 以內的電子,往上跳一點點就能碰到空軌域,才有辦法參與熱激發。對一般金屬,室溫 k_BT ≈ 0.026 eV,而 ε_F 是幾個 eV(p.139 Table 1:K 為 2.12 eV、Cu 為 7.00 eV),能動的只是費米海最表層薄薄的一層。
Fig. 5(p.140)在說什麼。這張圖畫的是三維自由電子的態密度 D(ε):一條從原點出發、像 ε^{1/2} 那樣上升並逐漸變緩的實線。T = 0 時被佔據的軌域是實線下方從 0 到 ε_F 的灰色區域,右緣在 ε_F 處垂直切斷。虛線是有限溫度(但 k_BT ≪ ε_F)下的佔據密度 f(ε,T)D(ε):遠低於 ε_F 時與實線重合,接近 ε_F 時從實線下方偏離,在 ε_F 附近平滑滑落,並在 ε_F 以上拖出一條很快趨近於零的尾巴。圖上把 ε_F 以下、實線與虛線之間被掏空的一小塊標為「1」,把 ε_F 以上、虛線下方新長出的一小塊標為「2」,並在費米能附近用箭頭標出寬度 k_BT。圖說的意思是:升溫時電子由區域 1 被熱激發到區域 2,平均能量因此增加;而灰色區域的主體完全不動。
有了這個圖像,估計就很自然。若總電子數為 N,能被熱激發的只佔比例約 T/T_F,因為只有它們落在分布頂端 ~k_BT 的範圍內(k_BT/ε_F = T/T_F,T_F ≡ ε_F/k_B 是 p.140 定義的費米溫度)。這 NT/T_F 個電子每個得到 ~k_BT 的熱能,於是得到式 (22) 與 (23)。
電子氣熱能的數量級估計
U_el ≈ (NT/T_F) k_BT
(22) p.141U_el:電子氣因升溫而得到的熱動能(Kittel 稱之為總電子熱動能;嚴格說是相對於 T = 0 的增量,後文記為 ΔU)。N:傳導電子總數。T:溫度。T_F ≡ ε_F/k_B:費米溫度(p.140 定義)。k_B:波茲曼常數。NT/T_F 是「能被熱激發的電子數」,k_BT 是「每個被激發電子得到的能量」。
電子熱容的數量級估計
C_el = ∂U/∂T ≈ Nk_B (T/T_F)
(23) p.141C_el:傳導電子對熱容的貢獻;其餘同式 (22)。
「約為古典值的 0.01 倍」的數值核對,以及數量級論證漏掉的 π²/2補推導
Kittel 在 p.141–142 給的「0.01 或更小」,是由定性式 (23) 得到的數量級陳述。這裡把它和本段最後的精確結果、以及實測值並排,看數量級論證到底差多少。
結論:「小約百分之一」在數量級上是對的——電子熱容確實只有古典值的百分之一到百分之五——但若把 (23) 當成定量公式,會系統性低估約 5 倍。物理直覺抓的是「哪個無因次比值控制現象」(這裡是 T/T_F),精確計算補上的是數值係數(這裡是 π²/2,源自積分 (33))。兩者都重要,但回答的是不同問題。第 5 步還顯示,若直接拿低溫實測 γ 外推,結果比自由電子的精確值再大 25%(K)到 38%(Cu)——這個差距正是後面熱有效質量 m_th/m 要處理的東西;不過這等於把低溫的 m_th/m 原封不動搬到 300 K,其中電子–聲子的部分在 T ≳ θ 時會消退,所以第 5 步的 0.050 與 0.017 應視為上限。
- 數量級版本:由 (23),C_el/(³⁄₂Nk_B) ≈ (Nk_BT/T_F)/(³⁄₂Nk_B) = ⅔·(T/T_F)。T_F = 5×10⁴ K、T = 300 K 時為 ⅔ × 0.006 = 0.004。
- 精確版本:用本段最後推出的 (36) C_el = ½π²Nk_BT/T_F,比值變成 C_el/(³⁄₂Nk_B) = (π²/3)(T/T_F)。同樣條件下 = 3.29 × 0.006 ≈ 0.020。
- 兩者相差 (36)/(23) = ½π² ≈ 4.93 倍(若用 (22) 直接微分得到的 2Nk_BT/T_F 來比,則差 π²/4 ≈ 2.47 倍)。這就是數量級論證的代價:T/T_F 這個標度抓對了,卻漏掉一個約 5 的數值係數。
- 代入具體金屬(T_F 取 p.139 Table 1):K 的 T_F = 2.46×10⁴ K → (π²/3)(300/24600) ≈ 0.040;Cu 的 T_F = 8.12×10⁴ K → ≈ 0.012。
- 改用實測 γ(p.146 Table 2,目視判讀:K 為 2.08、Cu 為 0.695 mJ/mol-K²):300 K 時 γT 分別為 0.624 與 0.209 J/mol-K,除以 ³⁄₂R = 12.47 J/mol-K,得 0.050 與 0.017。(300 K 仍遠低於 T_F,所以自由電子部分的線性律照樣成立;但實測 γ 還包含電子–聲子作用造成的質量增強(p.147 列出的第二個效應)。書外補充:這個增強只在 T ≪ θ 時完整存在,T ≳ θ 時會消退;K 的 θ 只有 91 K(p.116),300 K 遠高於它,Cu 的 θ = 343 K,300 K 已與之相當。因此這兩個數字是把低溫 γ 外推到 300 K 的略高估上限,「線性律照樣成立」只對去掉電子–聲子增強後的 γ 才較恰當。)
定量計算的起點:ΔU 的定義與三個要素
定性論證給了 C_el ∝ T 以及 T/T_F 的標度,但數值係數要靠真正的計算。Kittel 從 p.142 開始推導一個在低溫 k_BT ≪ ε_F 下成立的定量表達式。
出發點是總能量的增量 ΔU ≡ U(T) − U(0):N 個電子的系統從 0 加熱到 T 時,總能量增加了多少。U(T) 是把每個軌域的能量 ε 乘上它在溫度 T 被佔據的機率 f(ε)、再乘上單位能量的軌域數 D(ε),對所有能量積分;U(0) 是基態能量,此時 f 是一個在 ε_F 處從 1 掉到 0 的階梯,所以積分只到 ε_F,被積函數也不需要 f。兩者相減就是式 (24)。
三個要素先講清楚:
f(ε):費米–狄拉克分布(p.136 式 (5)),給出能量為 ε 的軌域在熱平衡下被佔據的機率。它含有化學勢 μ,而 μ 是溫度的函數,必須選得讓總粒子數恰好是 N(p.136)。T = 0 時 μ = ε_F。
D(ε):單位能量範圍內的軌域數(態密度)。三維自由電子由 p.140 式 (20) 給出,∝ ε^{1/2}。Kittel 在該頁註腳說明,嚴格來說它是單粒子態密度,或稱軌域密度。
ε_F:T = 0 時最高被佔據軌域的能量,是一個與溫度無關的固定參考能量。這一點後面非常關鍵:ε_F 不隨 T 變,μ 會。
升溫造成的總能量增量
ΔU = ∫₀^∞ dε ε D(ε) f(ε) − ∫₀^{ε_F} dε ε D(ε)
(24) p.142ΔU ≡ U(T) − U(0);ε:軌域能量;D(ε):態密度(含自旋,單位能量的軌域數);f(ε):溫度 T 時的佔據機率;ε_F:費米能。第一項是 U(T),第二項是 U(0)。
費米–狄拉克分布(式 (24a),附校註)
f(ε, T, μ) = 1 / {exp[(ε − μ)/k_BT] + 1}
(24a) p.142μ:化學勢(溫度的函數,T = 0 時等於 ε_F);k_BT:熱能;ε:軌域能量。f 是軌域被佔據的機率,介於 0 與 1 之間。
粒子數守恆恆等式
N = ∫₀^∞ dε D(ε) f(ε) = ∫₀^{ε_F} dε D(ε)
(25) p.142N:電子總數(固定)。中間的積分是溫度 T 時的總電子數,右邊是 T = 0 時的總電子數。
乘上 ε_F 後的恆等式
(∫₀^{ε_F} + ∫_{ε_F}^∞) dε ε_F f(ε) D(ε) = ∫₀^{ε_F} dε ε_F D(ε)
(26) p.142符號同 (25);ε_F 在此是乘上去的常數。
(24) 減去 (26) 得到 (27):書上省略的重組步驟補推導
(24) 的兩個積分都很大而且幾乎相消。Kittel 的策略是「加一個零」:用粒子數守恆 (25) 造出恆等於零的組合 (26),從 ΔU 中減掉。書上只用一句話帶過(用 (26) 把 (24) 改寫成 (27)),中間的重組完全省略。
減去 ε_F × N 的物理意義是把能量零點移到費米能。把一個電子從 ε_F 以下的 ε′ 搬到 ε_F 以上的 ε,系統能量增加 ε − ε′,可以拆成 (ε − ε_F) + (ε_F − ε′):前者是「把電子放到 ε_F 以上」的代價,後者是「在 ε_F 以下挖一個洞」的代價。(27) 正是這個拆法的積分版本。這個重組之所以合法,是因為被搬上去的電子數與挖出的洞數必須相等(粒子數守恆),所以 ε_F ×(上去的數目 − 洞的數目)= 0,改用 ε_F 當參考點不會改變總能量。若粒子數不守恆,這個「零」就不是零了。
- 由 (26),下面這個量恆等於零:Z ≡ ∫₀^{ε_F} dε ε_F f D + ∫_{ε_F}^∞ dε ε_F f D − ∫₀^{ε_F} dε ε_F D = 0。所以 ΔU = ΔU − Z,數值不變。
- 把 (24) 的第一個積分也拆成兩段:∫₀^∞ dε ε D f = ∫₀^{ε_F} dε ε D f + ∫_{ε_F}^∞ dε ε D f。
- 寫出 ΔU − Z,依積分範圍分組。ε_F 以上只有兩項:∫_{ε_F}^∞ dε ε D f − ∫_{ε_F}^∞ dε ε_F f D = ∫_{ε_F}^∞ dε (ε − ε_F) f(ε) D(ε)。這就是 (27) 的第一項。
- ε_F 以下有四項:∫₀^{ε_F} dε [ε f − ε − ε_F f + ε_F] D。
- 方括號因式分解:ε f − ε − ε_F f + ε_F = −ε(1 − f) + ε_F(1 − f) = (ε_F − ε)[1 − f(ε)]。
- 所以 ε_F 以下的部分 = ∫₀^{ε_F} dε (ε_F − ε)[1 − f(ε)] D(ε),即 (27) 的第二項。兩部分相加得 (27)。
- 檢查正負:第一項積分範圍 ε > ε_F,所以 ε − ε_F > 0,而 f ≥ 0、D ≥ 0;第二項範圍 ε < ε_F,所以 ε_F − ε > 0,而 1 − f ≥ 0。兩項都非負——這就是 Kittel 說「兩項對能量的貢獻都為正」的依據。
ΔU 改寫為激發電子與空缺兩項
ΔU = ∫_{ε_F}^∞ dε (ε − ε_F) f(ε) D(ε) + ∫₀^{ε_F} dε (ε_F − ε) [1 − f(ε)] D(ε)
(27) p.142第一項:f(ε)D(ε)dε 是能量 ε 附近 dε 範圍內「被提升上去的電子數」,(ε − ε_F) 是把一個電子從 ε_F 搬到 ε 所需的能量。第二項:1 − f(ε) 是軌域 ε 上的電子「已被移走」的機率,(ε_F − ε) 是把一個電子從 ε 帶到 ε_F 所需的能量。這兩個解讀都是 p.142 書上的說法。
Fig. 6:ΔU 隨溫度的完整曲線
Kittel 在 p.142 說 ΔU 畫在 Fig. 6(p.143)。這張圖是三維無交互作用費米子氣體的能量隨溫度的變化:縱軸是無因次的 ΔU/Nε_F,從 0 到 1.0;橫軸是 k_BT/ε_F,也從 0 到 1.0,格線每 0.2 一格。
曲線從原點出發,起初非常平坦、向上彎,像一段拋物線——這是低溫時 ΔU ∝ T² 的行為。大約過了 k_BT/ε_F ≈ 0.2(此處目視 ΔU/Nε_F 約 0.09)之後逐漸變直,成為一條斜率近乎固定的直線,經過約 (0.4, 0.3)、(0.6, 0.55)、(0.8, 0.82),並在 k_BT/ε_F ≈ 0.94 左右碰到圖框上緣 1.0。以上數值都是本站算繪 PDF 後的目視判讀。
Fig. 6 的定量檢查:低溫端 (π²/4)t² 與高溫端 ³⁄₂t − ⅗補推導
Fig. 6 是一條沒有附公式的曲線。這裡用本段自己的結果檢查它的兩端,順便看出低溫近似能用到哪裡。第 3、5 步的「精確數值」是本站以完整費米–狄拉克積分計算的,不是書上的數字;它與目視判讀在 0.2–0.8 之間都吻合到 0.01 以內,表示 Fig. 6 畫的是精確曲線而非示意圖。
要注意 Fig. 6 的橫軸範圍極大:k_BT/ε_F = 1 相當於把鉀加熱到約 2.5 萬 K。真實金屬在室溫下 k_BT/ε_F 約 0.01(Fig. 8 圖說),即使一路加熱到熔點也只到約 0.02 的量級(本站以 Table 1 的 T_F 與熔點估算:Cu 約 0.017、Au 約 0.021),也就是圖最左邊、幾乎貼著橫軸的那一小段;這張圖讓你看見整個「由簡併到古典」的過渡,金屬實際所處的只是曲線剛起步、呈二次方的那一角。
解讀:Fig. 6 同時畫了兩種物理——左端的 T² 是簡併費米氣體,右端的線性是古典氣體——T_F 就是兩者的交界尺度。第 3 步說明低溫展開可以放心用到 k_BT/ε_F ~ 0.1(誤差約 2%),對真實金屬(即使到熔點也只約 0.02)更是綽綽有餘。
- 低溫端:把本段最後的結果 (34) C_el = ⅓π²D(ε_F)k_B²T 對溫度積分(T = 0 時 ΔU = 0):ΔU = ∫₀^T C_el dT′ = ⅙π²D(ε_F)(k_BT)²。
- 代入 (35) D(ε_F) = 3N/2ε_F:ΔU = (π²/4)N(k_BT)²/ε_F,即 ΔU/Nε_F = (π²/4)t²,其中 t ≡ k_BT/ε_F,π²/4 ≈ 2.47。
- t = 0.2 時 (π²/4)(0.04) ≈ 0.099;Fig. 6 目視約 0.09;本站以完整費米積分(同時解 μ(T))數值計算得 0.092。t = 0.1 時低溫式 0.0247 對精確值 0.0243,誤差約 2%;t = 0.4 時低溫式 0.39 對精確值 0.30,已明顯高估——曲線在 t ≳ 0.2 以後「變直」就在這裡。
- 高溫端:古典極限 U → ³⁄₂Nk_BT,而 U(0) = ⅗Nε_F(第 6 章習題 1,p.157 式 (60))。所以 ΔU/Nε_F → ³⁄₂t − ⅗:斜率 ³⁄₂、截距 −0.6 的直線。
- t = 1 時古典漸近式給 0.9,精確數值約 1.10。差額來自 k_BT ~ ε_F 時仍然明顯的簡併效應:包立原理迫使粒子佔據較高的能量,使能量高於古典值。
- 金屬的實際工作點:K 在 300 K 時 t = 0.0122,ΔU/Nε_F = 2.47 × 1.49×10⁻⁴ ≈ 3.7×10⁻⁴,即 ΔU/N ≈ 7.8×10⁻⁴ eV——在 Fig. 6 上根本看不出離開了原點。
(27) → (28):微分後兩個積分為何能合成一個補推導
(27) 的兩項各以 ε_F 為積分邊界,看起來被切成兩段,但微分之後可以漂亮地合回一個積分。書上只說「整理各項」就得到 (28)。
(28) 和直接微分 (24) 的差別很值得體會。直接微分 (24) 會得到 ∫ ε(df/dT)D dε;其中與 ε_F 同量級的大部分,要靠 ∫(df/dT)D dε = dN/dT = 0 才消掉。(28) 則把這個消去預先做好了:兩者數學上相等(相差 ε_F·dN/dT = 0),但 (28) 的被積函數在遠離 ε_F 處直接為零,物理上透明,數值上也穩定得多。這是 (27)「以 ε_F 為零點」那一步留下的紅利。
- 對 (27) 求 d/dT。ε_F 由基態定義,與溫度無關;D(ε) 是自由電子的態密度,只與 ε 及體積有關(忽略熱膨脹,也就是算等容熱容)。所以唯一依賴溫度的是 f(ε)——這是 Kittel 在 (28) 之前的陳述。
- 第一項:d/dT ∫_{ε_F}^∞ dε (ε − ε_F) f D = ∫_{ε_F}^∞ dε (ε − ε_F)(df/dT) D。積分限與 T 無關,微分可以移進積分號。
- 第二項:d/dT ∫₀^{ε_F} dε (ε_F − ε)[1 − f] D = ∫₀^{ε_F} dε (ε_F − ε)(−df/dT) D = ∫₀^{ε_F} dε (ε − ε_F)(df/dT) D。負號與 (ε_F − ε) 合併成 (ε − ε_F)。
- 兩項的被積函數現在完全相同,都是 (ε − ε_F)(df/dT)D(ε),積分範圍分別是 [ε_F, ∞) 與 [0, ε_F],接起來正好是 [0, ∞)。得 C_el = dU/dT = ∫₀^∞ dε (ε − ε_F)(df/dT) D(ε),即式 (28)。
- 注意:這裡的 df/dT 是全導數。f 透過兩條路徑依賴 T:顯式的 k_BT,以及化學勢 μ(T)。(28) 在這一步仍是嚴格的;μ(T) 的處理留到 (30) 之前。
電子熱容的積分表示
C_el = dU/dT = ∫₀^∞ dε (ε − ε_F) (df/dT) D(ε)
(28) p.142df/dT:佔據機率對溫度的全導數;其餘同 (27)。
為什麼能把 D(ε) 提出來:Fig. 3 與「費米面附近的峰」
Kittel 在 p.142 說,在金屬感興趣的溫度下 τ/ε_F < 0.01(τ ≡ k_BT 已在第 5 章 p.107 式 (3) 定義,本章在 p.144 式 (30) 前重申;見後面的符號警示),並請讀者看 Fig. 3:(ε − ε_F)df/dT 在 ε_F 附近有很大的正峰。既然被積函數只在費米面附近一小段能量內不為零,而 D(ε) 在這一小段內幾乎不變,就可以用它在 ε_F 的值代替並提出積分號,得到式 (29)。
Fig. 3(p.136)在說什麼。這張圖畫的是費米–狄拉克分布 f(ε) 在多個溫度下的形狀,設定 T_F = 50,000 K、三維氣體,總粒子數固定。橫軸是 ε/k_B,以 10⁴ K 為單位,從 0 到 9;縱軸是 f,從 0 到 1.2。標註的溫度有 500 K、5,000 K、1×10⁴ K、2.5×10⁴ K、5×10⁴ K 與 10×10⁴ K。500 K 的曲線幾乎是一個在 ε/k_B = 5×10⁴ K 處垂直落下的階梯;5,000 K 在費米能附近變成一段平滑的下坡;1×10⁴ K 與 2.5×10⁴ K 更緩;500 K、5,000 K 與 1×10⁴ K 三條曲線大致交會在 ε/k_B ≈ 5、f 略小於 0.5 處。到了 5×10⁴ K(即 T = T_F),曲線在 ε = 0 處只從約 0.5 起跳;10×10⁴ K 時更只剩約 0.2——連最低的軌域都不再確定被佔據,階梯的痕跡完全消失。(以上讀值為本站目視判讀。)圖說提醒:各溫度的 μ 可以從 f = 0.5 的位置讀出。
Fig. 3 畫的是 f 而不是 df/dT。可以這樣理解:比較不同溫度的 f 可以看出,在金屬實際的溫度下,只有費米面附近的佔據機率會隨溫度改變,所以 df/dT 只在那裡顯著,而 (ε − ε_F)df/dT 在 ε_F 兩側都是正的。這個「峰」到底長什麼樣、有多寬,書上沒算,下一個推導區塊補上。
(ε − ε_F)df/dT 的精確形狀:兩個峰,位在 ε_F ± 2.4k_BT補推導
Kittel 說 (ε − ε_F)df/dT 在 ε_F 附近有 large positive peaks(複數),但沒有給出形狀與位置;以下算出兩峰位在 ±2.40k_BT。這關係到「D 能不能提出積分號」這個近似有多好,值得算清楚。
有兩個值得注意的細節。第一,這兩個峰對稱地位在費米面上下,而 ε_F 本身的被積函數是零:恰在費米面上的電子被激發時幾乎不需要能量,對熱容沒有貢獻。第二,窗口比定性論證說的「k_BT 範圍內」寬得多——要到 ±5k_BT 才涵蓋 85%。定性論證只在數量級上成立,這也是它漏掉 π²/2 係數的原因之一。
第 6 步說的一次項相消:若 D(ε) ≈ D(ε_F) + D′(ε_F)(ε − ε_F),一次項帶來 ∫ x³ eˣ/(eˣ + 1)² dx,被積函數是奇函數,積分為零。所以 D 的修正從二次項才開始,相對大小 ~ (k_BT/ε_F)²。
- 暫時取 μ = ε_F(理由見後面 Figs. 7–8 的討論),令 x = (ε − ε_F)/k_BT,則 f = 1/(eˣ + 1)。利用後面 (30) 的結果與 df/dT = k_B·df/dτ,可得 (ε − ε_F)(df/dT) = k_B · x² eˣ/(eˣ + 1)²。
- 函數 g(x) ≡ x² eˣ/(eˣ + 1)² 可改寫為 x²/[4cosh²(x/2)],是 x 的偶函數,在 x = 0 為零。
- 求極大:d ln g/dx = 2/x − tanh(x/2) = 0,即 x·tanh(x/2) = 2,數值解 x ≈ ±2.40,峰值 g ≈ 0.44。
- 半高:g = 0.22 的兩點在 |x| ≈ 1.08 與 |x| ≈ 4.50,所以每個峰從約 1k_BT 延伸到約 4.5k_BT。
- 權重分布(佔全積分 π²/3 的比例,本站數值積分):|x| < 1 只佔 4%;|x| < 3 約 50%;|x| < 5 約 85%;|x| < 8 約 98%。
- 結論:對熱容有貢獻的能量窗口大約是 ε_F ± 8k_BT。300 K 時為 ±0.2 eV;K 的 D ∝ ε^{1/2} 在這個窗口內的相對變化約 ½ × 0.2/2.12 ≈ 5%。而且 D 的一次項因對稱性相消(見下文),把 D(ε_F) 提出來的實際誤差遠小於 5%。
把態密度提出積分號
C_el ≅ D(ε_F) ∫₀^∞ dε (ε − ε_F) (df/dT)
(29) p.143D(ε_F):費米能處的態密度(單位能量的軌域數,含自旋);其餘同 (28)。
把 μ 換成 ε_F:Figs. 7 與 8 在說什麼
(29) 的 df/dT 還藏著一個麻煩:f 依賴化學勢 μ,而 μ 隨溫度改變。Kittel 在 p.143–144 的處理是:檢視 Figs. 7 與 8 中 μ 隨 T 的變化可知,當 k_BT ≪ ε_F 時可以忽略 μ 在分布函數裡的溫度依賴,直接用常數 ε_F 取代 μ。這一步書上是「看圖決定」的,所以兩張圖要讀仔細。
Fig. 7(p.143):三維無交互作用費米子氣體的化學勢 μ 對溫度(以 k_BT 表示)作圖。圖說特別註明,為了作圖方便,μ 與 k_BT 的單位都是 0.763ε_F(書上沒有說明這個數字的來歷)。橫軸 k_BT 從 0 到 3,格線每 0.5 一格;縱軸 μ 從 −6 到 +2。曲線在 k_BT = 0 處從 μ 約 +1.3 出發——這正是 ε_F 換算到此單位的值,ε_F/(0.763ε_F) ≈ 1.31。從 0 到 0.5 的區間塗成灰色,箭頭標為「簡併量子氣體區」,曲線在這段幾乎水平。之後曲線開始下彎,目視約在 k_BT ≈ 1.3–1.4(此單位下)穿過 μ = 0,然後近乎直線地一路下降,到 k_BT = 3 時約為 −4.5。右上方以箭頭標出「古典氣體區」。以上數值都是本站算繪後的目視判讀。
Fig. 8(p.144):放大看低溫區。縱軸是 μ/ε_F,只畫 0.95 到 1.05;橫軸是 τ/ε_F(τ ≡ k_BT),從 0 畫到 0.2 多一點。兩條曲線都從 (0, 1.00) 出發,起初都貼著 1.00 走得非常平,接著分岔:標示「One dimension」的往上彎,標示「Three dimensions」的往下彎。本站以像素量測判讀:τ/ε_F = 0.1 時一維約 1.008、三維約 0.991;τ/ε_F = 0.2 時一維約 1.036(約 1.04)、三維約 0.964(約 0.965)。圖說交代兩件事:一般金屬室溫下 τ/ε_F ≃ 0.01,所以 μ 非常接近 ε_F;這些曲線是由系統粒子數積分的級數展開算出來的。
Kittel 的論據就在這裡:在金屬實際所處的 τ/ε_F ≈ 0.01,兩條曲線都還緊貼著 1.00,以 Fig. 8 的解析度根本看不出偏離,所以把 μ 換成 ε_F。
為什麼一維往上、三維往下:D(ε) 的斜率決定 μ 的走向我的判讀
先補一句 Fig. 7 的解讀。μ 變成很大的負值代表什麼?回到 (24a):若 μ 遠低於 0,則對所有 ε ≥ 0 都有 ε − μ ≫ k_BT,f ≈ exp[(μ − ε)/k_BT],整個分布都落在 p.137 所說的波茲曼(馬克士威)尾巴裡——那就是古典氣體。所以 Fig. 7 同時是一張「簡併—古典」相圖。
Fig. 8 最耐人尋味的是兩條曲線方向相反,書上沒有解釋原因。這裡給一個只用 f 的對稱性與粒子數守恆的論證。
先假設升溫時 μ 停在 ε_F 不動。由 f(ε_F + δ) = 1 − f(ε_F − δ),費米面上方距離 δ 處一個軌域被佔據的機率,恰等於下方距離 δ 處一個軌域被清空的機率。所以:上方增加的電子數 ∝ D(ε_F + δ)·f,下方減少的電子數 ∝ D(ε_F − δ)·f。兩者是否相等,完全取決於 D 在費米面兩側是否一樣大。
三維自由電子 D ∝ ε^{1/2} 是遞增函數,D(ε_F + δ) > D(ε_F − δ):上方能容納的比下方釋出的多,若 μ 不動,總電子數會憑空增加。這違反粒子數守恆 (25),所以 μ 必須往下調,把分布整體往低能量推一點。這就是三維曲線往下彎。
一維自由電子的 D ∝ ε^{−1/2}(由 p.134 式 (3) ε_n = (ℏ²/2m)(nπ/L)²,每個 n 容納兩個自旋,能量 ≤ ε 的軌域數 2n ∝ ε^{1/2},再對 ε 微分),是遞減函數,情形剛好相反,μ 往上調。二維的 D 是常數(第 6 章習題 3,p.157),上下完全平衡,μ 在低溫下幾乎不動。習題 3 給的精確結果 μ = k_BT ln[exp(πnℏ²/mk_BT) − 1],用二維的 ε_F = πnℏ²/m(由 n = Dε_F)可改寫為 μ = ε_F + k_BT ln[1 − exp(−ε_F/k_BT)],偏離量是 exp(−ε_F/k_BT) 級,比任何冪次都小。
所以 μ(T) 的低溫走向是在「量」費米能附近態密度的斜率。同一個邏輯在第 8 章會以另一種面貌出現:p.207 式 (47) 給出本質半導體的 μ = ½E_g + ¾k_BT ln(m_h/m_e)——化學勢會往有效質量較小、也就是態密度較小的那一個帶偏移,理由與這裡完全相同。
化學勢的低溫展開(書外補充)
μ(T) ≈ ε_F − (π²/6)(k_BT)² · D′(ε_F)/D(ε_F)
p.144D′(ε_F) ≡ dD/dε 在 ε_F 的值;D′/D 是態密度在費米能處的對數斜率。
忽略 dμ/dT 為什麼安全:修正只有 (T/T_F)² 級補推導
Kittel 以「看 Figs. 7 和 8」作為把 μ 換成 ε_F 的理由。但圖只告訴你 μ 本身幾乎不變;熱容是導數,還需要確認 dμ/dT 的貢獻也可以忽略。這裡把論證補上。
真正讓這個近似成立的是對稱性:μ 的改變只是把整個 f 的「台階」平移,而平移對 (ε − ε_F) 加權的積分,只有在台階中心已經偏離 ε_F 時才有貢獻。因為偏離本身是 T² 級,它對熱容的影響就降到 T³ 級。這種「一次項因對稱相消」的結構在本段出現了三次——(27) 兩項相等、D′ 項消失、dμ/dT 項消失——都源自 f(ε_F + δ) = 1 − f(ε_F − δ)。
- f 對溫度的全導數分成兩部分:df/dT = (∂f/∂T)_μ + (∂f/∂μ)_T·(dμ/dT)。把 μ 換成常數 ε_F,等於丟掉第二部分。
- 第二部分代入 (29):D(ε_F)(dμ/dT)∫dε (ε − ε_F)(∂f/∂μ)。其中 ∂f/∂μ = −∂f/∂ε,是以 μ 為中心、面積為 1 的鐘形偶函數。
- 把 ε − ε_F 寫成 (ε − μ) + (μ − ε_F):前者是奇函數,積分為零;後者是常數,積分就等於 μ − ε_F。所以這一部分 = D(ε_F)(μ − ε_F)(dμ/dT)。
- 由上一個區塊,μ − ε_F ∝ T²,所以 dμ/dT ∝ T,這一部分 ∝ T³。與主項 ⅓π²D(ε_F)k_B²T 相比,相對大小 ~ (k_BT/ε_F)² = (T/T_F)²。
- 同理,第一部分裡若保留 μ 與 ε_F 的差或 D 的斜率,也只帶來 (T/T_F)² 級的相對修正。所以 C_el = γT[1 + O((T/T_F)²)];對室溫金屬,這個修正約 10⁻⁴。
符號警示:本段的 τ 不是碰撞時間;T_F 不是電子溫度;γ 不只一個我的判讀
τ 的兩種身分。Kittel 在 p.144 式 (30) 之前重申 τ ≡ k_BT(首次定義在第 5 章 p.107 式 (3))——本段的 τ 是以能量為單位的溫度。這是 Kittel 沿用 Kittel & Kroemer《Thermal Physics》的記號:附錄 D(p.652)稱 τ = k_BT 為基本溫度;第 5 章 p.114 的 Einstein 模型能量式也寫成 e^{ℏω/τ}。然而下一段(電導率與歐姆定律,p.147)在寫完 (39)–(41) 之後,以「If the collision time is τ」把 τ 用作碰撞時間,p.156 熱導率式裡的 τ 也是碰撞時間。相隔三頁(p.144 → p.147),τ 的量綱就從能量變成時間。讀 (30)–(33) 時,請一律把 τ 讀作 k_BT。
順帶一提:τ 在 p.142(τ/ε_F < 0.01 那句)就已經出現,比 p.144 的重申早兩頁,但這並不是先用後定義:p.142 的 τ/ε_F 沿用的是第 5 章 p.107 的同一個定義,即 k_BT/ε_F。
T_F 不是電子的溫度。p.145 開頭 Kittel 再次提醒:T_F 雖然叫費米溫度,卻不是電子的溫度,只是一個方便的參考記號(p.140 也講過一次)。電子氣的溫度就是晶格的溫度 T;T_F = ε_F/k_B 只是把費米能換成溫度單位,用來組成 T/T_F 這個比值——它衡量電子氣有多簡併。金屬的 T_F 是 10⁴–10⁵ K(Table 1),遠高於熔點,所以固態金屬中的電子氣在任何溫度下都是高度簡併的。
γ 至少有三個意思。本段的 γ 是 Sommerfeld 參數(電子熱容係數,p.145)。第 9 章緊束縛法(p.235)的 γ 則是重疊積分(一個能量),p.243 軌道量子化條件裡的 γ 又是相位修正。偏偏 Kittel 在重費米子處正要你去看第 9 章的緊束縛,兩個 γ 會在同一個論證中相遇,而且方向相反:重疊積分 γ 小 → 有效質量大 → Sommerfeld γ 大。
分布函數對 τ 的導數
df/dτ = [(ε − ε_F)/τ²] · exp[(ε − ε_F)/τ] / {exp[(ε − ε_F)/τ] + 1}²
(30) p.144τ ≡ k_BT:以能量為單位的溫度(p.144 在此式前重申,首次定義見第 5 章 p.107 式 (3);不是碰撞時間);f = 1/{exp[(ε − ε_F)/τ] + 1},已把 μ 換成 ε_F。
(30) 的微分,以及 df/dT = k_B·df/dτ補推導
書上直接寫出 (30),而且從 (29) 的 df/dT 換成 (30) 的 df/dτ 時沒有交代兩者差一個 k_B。這一步沒注意到的話,到 (32) 就會對不上 k_B 的冪次。
- 令 u ≡ (ε − ε_F)/τ,則 f = (eᵘ + 1)⁻¹。
- 鏈鎖律:df/dτ = −(eᵘ + 1)⁻² · eᵘ · (du/dτ)。
- du/dτ = −(ε − ε_F)/τ²(ε 與 ε_F 都與 τ 無關)。
- 代入:df/dτ = [(ε − ε_F)/τ²] · eᵘ/(eᵘ + 1)²,兩個負號相乘為正,即式 (30)。
- (29) 要的是 df/dT。由 τ = k_BT,dτ/dT = k_B,所以 df/dT = k_B·df/dτ。這個多出來的 k_B,就是 (32) 中 k_B² 的其中一個 k_B。
無因次能量變數
x ≡ (ε − ε_F)/τ
(31) p.144x:以熱能 τ = k_BT 為單位、從費米能起算的能量。
(29)+(30)+(31) → (32):k_B²T 從哪裡來,以及下限換成 −∞ 的誤差補推導
這是把 (29) 化為一個純數積分的關鍵一步。書上只說由 (29) 與 (30) 得到 (32),中間的換變數與 k_B 的來源都沒寫。
第 3 步是整個計算最漂亮的地方:所有溫度依賴都從被積函數中消失,全部集中到前面的 k_B²T。這保證 C_el ∝ T,而且比例常數 = D(ε_F) × k_B² × 一個普適純數。溫度的一次方可以拆成兩個來源:df/dT 帶來的 k_B(每單位溫度的能量尺度),以及 dε = τdx 帶來的 τ(有貢獻的能量窗口寬度 ∝ T)——後者正是定性論證中「被激發的電子數 ∝ T」的精確版本。
- 從 (29) C_el = D(ε_F)∫₀^∞ dε (ε − ε_F)(df/dT),代入 df/dT = k_B·df/dτ:C_el = k_B D(ε_F)∫₀^∞ dε (ε − ε_F)(df/dτ)。
- 換變數 (31):ε = ε_F + τx,dε = τ dx;ε = 0 ↔ x = −ε_F/τ,ε → ∞ ↔ x → ∞。
- 被積函數:(ε − ε_F)(df/dτ) = (τx) · (τx/τ²) · eˣ/(eˣ + 1)² = x² eˣ/(eˣ + 1)²——τ 完全消掉。
- 再乘上 dε = τ dx:C_el = k_B D(ε_F) · τ ∫_{−ε_F/τ}^∞ dx x² eˣ/(eˣ + 1)²。
- τ = k_BT,所以前面的係數 k_B·τ = k_B²T:C_el = k_B²T D(ε_F)∫_{−ε_F/τ}^∞ dx x² eˣ/(eˣ + 1)²,即式 (32)。
- 下限換成 −∞ 的誤差:x → −∞ 時被積函數 ≈ x²eˣ。被丟掉的部分 ∫_{−∞}^{−a} x² eˣ dx = (a² + 2a + 2)e^{−a},其中 a = ε_F/τ。a = 100 時約 1.0×10⁴ × 3.7×10⁻⁴⁴ ≈ 4×10⁻⁴⁰,相對於積分值 π²/3 ≈ 3.29 完全可忽略。這是 Kittel 說「eˣ 因子在 x = −ε_F/τ 已經可以忽略」的定量版本。
換變數後的電子熱容
C_el = k_B² T D(ε_F) ∫_{−ε_F/τ}^∞ dx x² eˣ/(eˣ + 1)²
(32) p.144x:由 (31) 定義;τ = k_BT;D(ε_F):費米能處的態密度。
∫ x² eˣ/(eˣ + 1)² dx = π²/3 的完整計算補推導
Kittel 只給結果 (33),沒有推導。這個積分會在固態物理裡反覆出現(熱電、熱導、Wiedemann–Franz 常數),值得完整走一遍。唯一用到的外部結果是 Σ_{n≥1} 1/n² = π²/6。
這和第 5 章 Debye T³ 律的積分(p.114:∫₀^∞ x³/(eˣ − 1) dx = π⁴/15)是同一套技巧——幾何級數展開後逐項積分。差別在於:玻色子的 1/(eˣ − 1) 展開後各項同號,得到 ζ 函數本身;費米子的 1/(eˣ + 1) 展開後各項交錯,得到交錯版本 (1 − 2^{1−s})ζ(s),所以多出第 6 步的那個 ½。這個 π²/3 會一路傳到 Wiedemann–Franz 定律的 Lorenz 數 L = (π²/3)(k_B/e)²(p.157 式 (59))——那個 π²/3 就是這裡的 (33)。
- 偶函數:eˣ/(eˣ + 1)² = 1/(e^{x/2} + e^{−x/2})² = 1/[4cosh²(x/2)],在 x → −x 下不變,x² 也是偶函數。所以 I ≡ ∫_{−∞}^∞ x² eˣ/(eˣ + 1)² dx = 2∫₀^∞ x² eˣ/(eˣ + 1)² dx。
- 分部積分:注意 eˣ/(eˣ + 1)² = −d/dx[1/(eˣ + 1)]。於是 ∫₀^∞ x² · {−d/dx[1/(eˣ + 1)]} dx = [−x²/(eˣ + 1)]₀^∞ + ∫₀^∞ 2x/(eˣ + 1) dx。
- 邊界項:x = 0 時 x² = 0;x → ∞ 時 x²/(eˣ + 1) → 0。邊界項為零,所以 ∫₀^∞ x² eˣ/(eˣ + 1)² dx = 2∫₀^∞ x/(eˣ + 1) dx。
- 幾何級數:x > 0 時 e⁻ˣ < 1,1/(eˣ + 1) = e⁻ˣ/(1 + e⁻ˣ) = Σ_{n=1}^∞ (−1)^{n+1} e^{−nx}。
- 逐項積分:∫₀^∞ x e^{−nx} dx = 1/n²。所以 ∫₀^∞ x/(eˣ + 1) dx = Σ_{n=1}^∞ (−1)^{n+1}/n² = 1 − 1/4 + 1/9 − 1/16 + …
- 交錯級數與 Σ1/n² = π²/6 的關係:Σ(−1)^{n+1}/n² = Σ1/n² − 2Σ_{n 偶}1/n² = Σ1/n² − 2·(1/4)Σ1/m² = ½Σ1/n² = ½ · π²/6 = π²/12。
- 回代:∫₀^∞ x² eˣ/(eˣ + 1)² dx = 2 · π²/12 = π²/6,I = 2 · π²/6 = π²/3。即式 (33)。
費米積分
∫_{−∞}^{∞} dx x² eˣ/(eˣ + 1)² = π²/3
(33) p.144x:無因次能量 (ε − ε_F)/k_BT。
電子熱容的一般式
C_el = ⅓π² D(ε_F) k_B² T
(34) p.144D(ε_F):費米能處的態密度(含自旋;若 D 以每莫耳計,C_el 即為莫耳熱容)。
自由電子在費米能處的態密度
D(ε_F) = 3N/2ε_F = 3N/2k_BT_F
(35) p.144N:傳導電子總數;k_BT_F ≡ ε_F。
自由電子氣的熱容
C_el = ½π² Nk_B T/T_F
(36) p.144N:傳導電子總數;T_F ≡ ε_F/k_B:費米溫度。
γ 真正量的是 D(ε_F):本段結果在全書的回響我的判讀
把 (34) 寫成 γ = ⅓π²k_B²D(ε_F),就會看出本段最重要的觀念轉換:電子熱容量的不是「有多少電子」,而是「費米能處每單位能量有多少軌域」。N 個電子裡絕大多數被包立原理凍結在費米海深處,對熱容毫無貢獻;決定反應的只有費米面上的可用態數 D(ε_F)。這個觀點會在後面反覆出現:
熱導率與 Wiedemann–Franz 定律(本章 p.156–157):Kittel 把 (36) 代入 K = ⅓Cvℓ 得到電子熱導率 K_el = π²nk_B²Tτ/3m(這裡的 τ 是碰撞時間),再除以電導率得到 Lorenz 數 L = (π²/3)(k_B/e)²。那個 π²/3 就是本段的積分 (33)。
包立自旋順磁性(第 11 章 p.316–317):同樣的論證——只有費米面附近 k_BT 範圍內的電子能翻轉自旋,比例 T/T_F——並直接引用本章的 D(ε_F) = 3N/2ε_F = 3N/2k_BT_F。Kittel 在 p.317 還指出,過渡金屬的磁化率偏高,暗示其態密度異常地高,這與電子熱容的量測一致。
電子–電子碰撞(第 14 章 p.417–419):碰撞率被包立原理壓低 (k_BT/ε_F)² 倍,理由也是只有費米面附近的電子有空的末態可去。
超導(第 10 章 p.266–267):鎵的正常態熱容寫成 C/T = 0.596 + 0.0568T²,其中 0.596 正是 Table 2 裡鎵的實測 γ;進入超導態後,電子熱容不再線性於 T,而是隨能隙呈指數變化。線性律的消失本身就是費米面上打開能隙的證據。
合金與過渡金屬(第 22 章 p.635):Kittel 以鎳為例,指出 3d 帶的平均態密度比 4s 帶高一個數量級,並說這大致反映了非鐵磁性過渡金屬的電子熱容相對於簡單一價金屬的增強。
實驗:C = γT + AT³ 與 C/T 對 T² 作圖(Fig. 9)
理論預測有了,接著是實驗(p.145「Experimental Heat Capacity of Metals」)。在同時遠低於 Debye 溫度 θ 與費米溫度 T_F 的溫度下,金屬的熱容可以寫成電子與聲子兩項之和:C = γT + AT³,γ 與 A 是材料常數。電子項線性於 T,在溫度夠低時主宰。
Kittel 建議把數據畫成 C/T 對 T²:兩邊除以 T 得 C/T = γ + AT²(式 (37)),若以 T² 為橫軸,數據點應落在一條直線上,斜率是 A,截距是 γ。
Fig. 9(p.145)在說什麼。這是鉀的實驗數據(引自 W. H. Lien 與 N. E. Phillips)。縱軸 C/T,單位 mJ/mol-K²,範圍 2.0 到 3.0;橫軸 T²,單位 K²,範圍 0 到 0.3。圖上有二十多個黑點,分布在 T² 約 0.06 到 0.3 K² 之間(也就是 T 約 0.25 K 到 0.55 K,本站目視判讀),全部整齊地落在一條直線上;直線向左延伸到 T² = 0 處。圖內標注擬合式 C/T = 2.08 + 2.57T²。換言之,鉀的 γ = 2.08 mJ/mol-K²,A = 2.57 mJ/mol-K⁴。
Kittel 把觀測到的 γ 稱為 Sommerfeld 參數,各金屬的值列在 Table 2。
低溫金屬熱容的線性化形式
C/T = γ + AT²
(37) p.145C:金屬的總熱容(每莫耳);γ:電子熱容係數(Sommerfeld 參數),理論上 γ = ⅓π²D(ε_F)k_B²;A:聲子 T³ 項的係數,Debye 模型下每莫耳原子 A = (12π⁴/5)R/θ³。
從 Fig. 9 反推:鉀的 Debye 溫度與電子/聲子交叉溫度補推導
C/T 對 T² 作圖是一個標準的線性化技巧:γ 不必從一條彎曲的 C(T) 曲線去猜,而是由直線外插到 T² = 0 直接讀出;數據是否共線,本身也在檢驗 C = γT + AT³ 在該溫區是否成立。Fig. 9 的直線同時給出截距與斜率。截距 γ 是本段主角;斜率 A 則連回第 5 章的 Debye 模型。這裡用斜率做一個跨章一致性檢查,並算出電子項何時勝出——這也解釋了為什麼 Fig. 9 的量測溫度只有半個 K 左右。
鉀的電子與聲子熱容(mJ/mol-K),由 Fig. 9 的擬合式計算:
| T (K) | γT | AT³ | 電子項占比 |
|---|---|---|---|
| 0.3 | 0.62 | 0.07 | 90% |
| 0.5 | 1.04 | 0.32 | 76% |
| 1 | 2.08 | 2.57 | 45% |
| 2 | 4.16 | 20.6 | 17% |
| 4 | 8.32 | 164 | 5% |
第 3 步的吻合說明 Fig. 9 的直線確實能拆成物理上獨立的兩部分:電子的 γT 與聲子的 AT³。第 4、5 步則點出實驗上的現實:鉀的聲子很「軟」(θ 只有 91 K),聲子熱容大,要降到 1 K 以下電子項才冒出頭;θ 較高的銅,交叉溫度提高到約 4 K。這也說明為什麼 γ 一定要在液氦溫區甚至更低溫量。
- Debye T³ 律(第 5 章 p.114 式 (32)),每莫耳原子:C_ph = (12π⁴/5)R(T/θ)³,所以 A = (12π⁴/5)R/θ³ ≈ 1944 J/mol-K ÷ θ³。
- 解 θ:θ = [(12π⁴/5)R/A]^{1/3} = [1943.8 / (2.57×10⁻³)]^{1/3} K = (7.56×10⁵)^{1/3} K ≈ 91 K。
- 對照第 5 章 Table 1(p.116;該表在 PDF 中是無文字層的影像,本站算繪後目視判讀):鉀的 θ 為 91 K。從熱容直線的斜率反推出的 θ 與表列值一致。
- 交叉溫度:γT* = AT*³ ⇒ T* = (γ/A)^{1/2} = (2.08/2.57)^{1/2} K ≈ 0.90 K。低於約 0.9 K,電子熱容大於聲子熱容。
- 作為對照,銅的 θ = 343 K(同表目視判讀),A = 1943.8/343³ ≈ 4.8×10⁻⁵ J/mol-K⁴ = 0.048 mJ/mol-K⁴;配合 Table 2 的 γ = 0.695 mJ/mol-K²,T* ≈ (0.695/0.048)^{1/2} ≈ 3.8 K。
重現 Table 2 的「自由電子 γ」:每莫耳原子、T_F 取自 Table 1補推導
Kittel 說自由電子 γ 是用 (17) 與 (34) 算出的,但沒說 Table 2 用的是哪一組 ε_F,也沒說「mol」指什麼。這裡用 p.139 Table 1 的 T_F 直接重算。
完全吻合告訴我們兩件事:(1) Table 2 的「mol」是每莫耳原子,所以多價金屬的自由電子 γ 要乘上價數 z;(2) Table 2 與 Table 1 用的是同一組自由電子參數,兩張表可以對照著讀。也就是說,Table 2 的第二列不是獨立資料,而是把 Table 1 的 T_F 代進 (36) 的結果;真正的實驗資訊只在第一列。
- 由 (36),每莫耳原子(含 N = zN_A 個傳導電子,z 為價數):γ(free) = π²zN_Ak_B/(2T_F) = π²zR/(2T_F),R = 8.314 J/mol-K。
- K(z = 1,T_F = 2.46×10⁴ K):π² × 8.314/(2 × 24600) = 1.668 mJ/mol-K²。Table 2 目視:1.668。
- Cu(z = 1,8.12×10⁴ K):0.505,表列 0.505。Li(5.48×10⁴ K):0.749,表列 0.749。Na(3.75×10⁴ K):1.094,表列 1.094。
- 多價金屬要乘 z:Be(z = 2,16.41×10⁴ K):0.500,表列 0.500;Al(z = 3,13.49×10⁴ K):0.912,表列 0.912;Pb(z = 4,10.87×10⁴ K):1.510,表列 1.509。
- 第三列 m_th/m 也逐一驗算為「第一列 ÷ 第二列」:K 為 2.08/1.668 = 1.247 → 1.25;Cu 為 0.695/0.505 = 1.376 → 1.38;Li 為 1.63/0.749 = 2.18;Be 為 0.17/0.500 = 0.34。本站對 Table 2 中三列都有數值的 22 個金屬全部驗算過,捨入後一致——這同時是對本站目視判讀的交叉檢查。
熱有效質量的定義
m_th/m ≡ γ(observed)/γ(free)
(38) p.145m_th:熱有效質量(thermal effective mass);m:自由電子質量;γ(observed):實測 Sommerfeld 參數;γ(free):以自由電子質量 m 與同一電子濃度,由 (17) 與 (34) 算出的值。
為什麼 γ ∝ m:兩種看法補推導
Kittel 用一句話說明 (38) 的形式很自然:ε_F 與電子質量成反比,所以 γ ∝ m。這裡把兩種等價的理由寫開。
看法二比較有物理:固定電子濃度就固定了費米球的大小 k_F;質量愈大,ε = ℏ²k²/2m 愈平,同一個 k 空間殼層對應的能量間隔愈小,單位能量裡塞進的軌域就愈多。「重」的電子 = 平坦的能帶 = 高的 D(ε_F) = 大的 γ。這條等價鏈是理解重費米子的關鍵:γ 大兩三個數量級,意思就是費米能處的能帶極其平坦。
- 看法一(Kittel p.145 的說法):由 (17) ε_F = (ℏ²/2m)(3π²N/V)^{2/3},在電子濃度 N/V 固定時 ε_F ∝ 1/m,因此 T_F ∝ 1/m。
- 代入 (36):γ = ½π²Nk_B/T_F ∝ 1/T_F ∝ m。
- 看法二(更根本):把 ε_F = ℏ²k_F²/2m 代入 (20),D(ε_F) = (V/2π²)(2m/ℏ²)^{3/2}(ℏk_F/√(2m)) = Vmk_F/(π²ℏ²)。k_F = (3π²N/V)^{1/3} 只由電子濃度決定((16)),與質量無關,所以 D(ε_F) ∝ m。
- 代入 (34):γ = ⅓π²k_B²D(ε_F) ∝ m。
- 所以 m_th/m = γ(obs)/γ(free) = D_obs(ε_F)/D_free(ε_F):熱有效質量其實是費米能處態密度的比值,只是用「質量」的語言表達。
Table 2 的內容(p.146,影像表,目視判讀)
Table 2 的標題是金屬電子熱容常數 γ 的實驗值與自由電子值;註明數據來自 N. Phillips 與 N. Pearlman 提供的彙編,熱有效質量依式 (38) 定義。表格依週期表排列,每個元素格內有三列:第一列是觀測 γ(mJ mol⁻¹ K⁻²),第二列(灰底)是自由電子計算 γ(同單位),第三列是 m_th/m = 觀測 γ / 自由電子 γ。
關於資料來源的說明:p.146 在這份 PDF 中是一張旋轉 90° 的影像,沒有文字層(文字抽取只得到頁碼)。以下數值是本站以 300 dpi 算繪、轉正後目視判讀的結果;有三列數值的金屬都已用「第三列 = 第一列 ÷ 第二列」驗算自洽。
以下選錄三列齊全的簡單金屬(價數欄取自 p.139 Table 1,不在 Table 2 內):
| 金屬 | 價數 | 觀測 γ | 自由電子 γ | m_th/m |
|---|---|---|---|---|
| Li | 1 | 1.63 | 0.749 | 2.18 |
| Na | 1 | 1.38 | 1.094 | 1.26 |
| K | 1 | 2.08 | 1.668 | 1.25 |
| Rb | 1 | 2.41 | 1.911 | 1.26 |
| Cs | 1 | 3.20 | 2.238 | 1.43 |
| Cu | 1 | 0.695 | 0.505 | 1.38 |
| Ag | 1 | 0.646 | 0.645 | 1.00 |
| Au | 1 | 0.729 | 0.642 | 1.14 |
| Be | 2 | 0.17 | 0.500 | 0.34 |
| Zn | 2 | 0.64 | 0.753 | 0.85 |
| Al | 3 | 1.35 | 0.912 | 1.48 |
| Ga | 3 | 0.596 | 1.025 | 0.58 |
| Pb | 4 | 2.98 | 1.509 | 1.97 |
過渡金屬在表中只有觀測值,第二、三列留空。它們的 γ 多半明顯偏大,例如 Sc 10.7、Y 10.2、V 9.26、Pd 9.42、Ni 7.02、Pt 6.8;但也有不大的,如 Cr 1.40、Mo 2.0、W 1.3。半金屬則極小:As 0.19、Sb 0.11、Bi 0.008。B、C、N、Si、P、Ge 的格子是空的。
偏離 1 的三個來源
Kittel 說 m_th/m 偏離 1 牽涉三個不同的效應(p.145–147):
(1) 能帶效應:傳導電子與剛性晶格的週期位勢交互作用。電子在這個位勢中的有效質量稱為能帶有效質量(band effective mass)。
(2) 電子–聲子交互作用:電子會把鄰近的晶格極化或扭曲,於是移動中的電子試圖拖著附近的離子一起走,提高了電子的有效質量。
(3) 電子–電子交互作用:移動中的電子在周圍電子氣中引起慣性反作用,同樣提高有效質量。
用這三個效應讀 Table 2我的判讀
這三項正好是自由電子氣刻意忽略的三件事:晶格的靜態週期位勢、晶格的動態響應(聲子)、電子彼此之間的庫侖作用。本段只列出名目。
跨章對照。三個效應在書中後面都有正式處理。能帶有效質量:第 8 章 p.199 指出能量隨 k 變化很小時有效質量就很大,並說第 9 章的緊束縛近似能快速說明窄能帶如何形成;第 9 章 p.235 給出簡單立方的緊束縛結果——能量被限制在寬 12γ 的能帶內(此處 γ 是重疊積分),帶底有效質量 m* ≅ ℏ²/2γa²,重疊愈弱、能帶愈窄、有效質量愈大。電子–聲子:第 14 章 p.420 的極化子(polaron)正是「電子拖著離子走而變重」的圖像,並給出極化子質量與耦合常數 α 的關係。電子–電子:第 14 章 p.417 以 Landau 費米液體理論處理,Kittel 在那裡用了和 p.147 幾乎相同的描述,並說在鹼金屬中這項質量增加大約是 25% 的量級。
以下是本站對 Table 2 的解讀,屬於教學性的推論,不是書上的陳述。
鹼金屬(Na、K、Rb 約 1.25):第 9 章 p.230 說鹼金屬的費米面只佔布里淵區體積一半、離區邊界很遠,Na 的費米面非常接近球形,Cs 則偏離球形約 10%。能帶效應既小,m_th/m ≈ 1.25 就與 p.417 的「電子–電子效應約 25%」量級相符;Cs 的 1.43 較大,也與它較明顯的非球形一致。但要小心:m_th 是三個效應的總和,電子–聲子同樣有貢獻,光憑熱容無法歸因。Li 的 2.18 明顯異常,說明即使同是鹼金屬,也不能假設能帶效應一定小。
貴金屬(Cu 1.38、Ag 1.00、Au 1.14):Ag 與 Au 的電子濃度幾乎相同(Table 1:5.85 vs 5.90×10²² cm⁻³),m_th/m 卻是 1.00 vs 1.14,提示費米面附近的能帶細節對 D(ε_F) 有實質影響。
小於 1 的例子(Be 0.34、Ga 0.58、Zn 0.85):D(ε_F) 比自由電子還小。Kittel 在 p.230 說 Be、Mg 的晶格作用弱、費米面近乎球形,但因為各有兩個價電子,費米球的體積恰等於一個區,必然跨出第一布里淵區。區邊界處打開的能隙會改變費米能附近的態密度;Be 的 0.34 顯示這類效應在熱容上一點也不小。至於各金屬實際由哪個機制主導,本段的資訊不足以判斷。
過渡金屬(Sc、Y、V、Pd 約 9–11):比一價的貴金屬大一個數量級以上。這與書中兩處說法一致:第 11 章 p.317 說過渡金屬的高磁化率暗示態密度異常地高,且與電子熱容量測相符;第 22 章 p.635 以鎳說明 3d 帶的平均態密度比 4s 帶高一個數量級。窄的 d 帶 = 平坦 = 高 D(ε_F) = 大 γ。
半金屬(Bi 0.008、Sb 0.11、As 0.19):γ 極小,因為傳導載子極少。第 8 章 p.186 的 Fig. 1 把 Bi、Sb、As 放在金屬與半導體之間的載子濃度區,p.215 說明半金屬的導帶底只比價帶頂略低,兩帶只有很小的重疊。載子少、D(ε_F) 小,γ 就小。
重費米子
p.147 最後一小段談一類極端情形。有幾種金屬化合物被發現具有巨大的電子熱容常數 γ——比一般值高兩到三個數量級。這些重費米子化合物包括 UBe₁₃、CeAl₃ 與 CeCu₂Si₂。Kittel 用保留的語氣寫道:有人提出,這些化合物中的 f 電子可能具有高達 1000m 的慣性質量,原因是相鄰離子上 f 電子波函數的重疊很弱;並請讀者參見第 9 章的緊束縛。
重費米子的數量級:m* = 1000m 意味著什麼我的判讀
p.147 那段話雖短,卻把本段的線索全部串起來:γ ∝ D(ε_F) ∝ m_th;而重疊弱 → 能帶窄 → 有效質量大(第 9 章緊束縛)→ D(ε_F) 大 → γ 大。它可以看成「三個效應」中能帶效應被推到極端的例子。
用本段的公式做幾個數量級估計,可以感受「1000m」的份量。以下是本站的示意計算,不是這些化合物的真實模型。
γ 的量級。一般金屬的 γ ~ 1 mJ/mol-K²(Table 2 的簡單金屬多在 0.6–3 之間)。大兩三個數量級,就是 ~10²–10³ mJ/mol-K²,也就是每莫耳每 K² 可接近 1 J。
費米溫度的量級。電子濃度固定時 T_F ∝ 1/m(由 (17))。拿鉀的 T_F = 2.46×10⁴ K 當基準,質量變成 1000 倍時 T_F* ~ 25 K。後果很重要:(36) 的線性律要求 T ≪ T_F,對重費米子而言意味著要到幾 K 甚至更低才看得到 C = γT;在幾十 K,這些電子已經不能當作簡併費米氣體處理。
能帶寬度的量級。用第 9 章 p.235 簡單立方緊束縛的 m* ≅ ℏ²/2γa²(γ 為重疊積分):取 a = 5 Å,要 m* = 1000m,需要 γ = (ℏ²/2m)/(1000a²) = 3.81 eV·Å² / (1000 × 25 Ų) ≈ 1.5×10⁻⁴ eV,能帶寬 12γ ≈ 1.8 meV,相當於約 21 K。與上面的 T_F* 同量級——兩種估計殊途同歸:質量大 1000 倍,就是整條能帶窄到只有幾十 K。
對 (29) 的意涵。把 D(ε) 提出積分號的前提,是 D 在 ~k_BT 範圍內不變。當整條能帶只有幾十 K 寬,這個前提只在極低溫成立;這也是重費米子的 γ 必須在很低溫下量的原因。
順帶釐清:「重費米子」的「重」不是電子本身變重,而是費米能處的能譜極其平坦,用 (38) 的語言表達成很大的熱有效質量。
常見誤解我的判讀
| 常見誤解 | 正確理解 |
|---|---|
| 費米面以下的電子沒有動能,所以不貢獻熱容 | 它們的動能很大(平均 ⅗ε_F,數 eV),只是周圍軌域全滿、無法改變狀態,所以不能吸熱 |
| T_F 是電子的溫度 | p.140 與 p.145 兩度強調不是;電子溫度就是晶格溫度 T,T_F 只是 ε_F 的溫度單位寫法 |
| (30)–(33) 裡的 τ 是碰撞時間 | 本段 τ ≡ k_BT(p.144);碰撞時間 τ 從下一段 p.147 才出現 |
| 只有恰好位於 ε_F 的電子貢獻熱容 | 貢獻集中在 ε_F ± 2–5k_BT,兩個峰在 ±2.4k_BT;ε_F 本身的貢獻為零 |
| 電子熱容正比於電子總數 N | 正比於 D(ε_F);過渡金屬的 γ 可比貴金屬大一個數量級以上 |
| m_th 就是能帶有效質量 | m_th 由 (38) 定義,同時包含能帶、電子–聲子、電子–電子三種效應 |
| 在室溫量熱容就能得到 γ | 室溫聲子熱容約 3R;電子部分在簡單金屬約 1–3%,過渡金屬可達約 10%,無論哪一類,都不能從室溫熱容直接分離出 γ;γ 要在遠低於 θ 的溫度用 C/T 對 T² 外插 |
| 升溫時 μ 一定下降 | 取決於 D′(ε_F) 的正負:三維下降、一維上升、二維幾乎不變(Fig. 8、習題 3) |
| (23) 與 (36) 的差別只是嚴謹度,數值差不多 | 兩者差 ½π² ≈ 4.9 倍;(23) 只能用來估數量級 |
Cth 研究接點:室溫下電子熱容可以忽略到什麼程度延伸知識
對 FinFET 自熱熱容(Cth)的研究,本段給出一個可以直接引用的數量級結論:室溫下,Cu、W、Al 這類金屬中電子對熱容的貢獻只有約 1–2%;過渡金屬(Ti、Ru、Co、Ta、Pt、Ni、V、Pd 等)以低溫 γ 估算則可達約 4–11%。元件熱容仍以晶格(聲子)為主。
以 Table 2 的實測 γ(目視判讀)估算 300 K 的電子熱容,並與 Dulong–Petit 值 3R ≈ 24.9 J/mol-K 比較(直接沿用低溫 γ,並忽略 C_p 與 C_v 的差):Cu 的 γT = 0.695×10⁻³ × 300 ≈ 0.21 J/mol-K,約 0.8%;W 為 1.3×10⁻³ × 300 ≈ 0.39 J/mol-K,約 1.6%;Al 為 1.35×10⁻³ × 300 ≈ 0.41 J/mol-K,約 1.6%。但 FinFET 接觸與阻障層常用的過渡金屬不能套用同一個數字:以同樣算法,Table 2 中這些純金屬的 γ 給出 Ti(3.35)與 Ru(3.3)約 4.0%、Co(4.73)約 5.7%、Ta(5.9)約 7.1%、Pt(6.8)約 8.2%、Ni(7.02)約 8.4%,V(9.26)與 Pd(9.42)更達約 11%(括號內為 γ,單位 mJ mol⁻¹ K⁻²)。這些是以低溫 γ 外推的估計:低溫 γ 含電子–聲子增強,而這個增強在 300 K 會消退,所以實際值會再低一些,但一般仍高於 Cu、W、Al 的 1–2%。至於矽,載子濃度比金屬低好幾個數量級(即使重摻雜的源/汲極或強反轉層中的電子可能是簡併的),電子熱容仍然微不足道。
實務意涵:在 Cth 萃取中,把熱容建模為「晶格比熱 × 體積」大體是有根據的;但「電子熱容的修正小於一般量測不確定度」只對 Cu、W、Al 這類金屬成立,過渡金屬層的電子項不宜一概忽略。只有在極短時間尺度、電子與晶格尚未達成熱平衡(電子溫度 ≠ 晶格溫度)時,才需要把電子熱容單獨拿出來考慮——那已超出 Kittel 本段的範圍,屬於書外內容。