Kittel·固態物理

全書逐章 / 第 5 章

05

聲子 II:熱性質

Phonons II. Thermal Properties

書本 p.105–1296 個小節62 條公式53 段推導約 100,962 字熱容關聯 直接

①主要內容與研究問題

第 5 章要解決的問題,Kittel 在 p.107 章首用一句話講完:We discuss the heat capacity of a phonon gas and then the effects of anharmonic lattice interactions on the phonons and on the crystal. 這一句同時也是全章的目錄——前半講熱容,後半講非諧性帶來的後果(熱膨脹與熱導率)。

前半的核心問題只有一個:把第 4 章交出的「模式清單」換算成量熱計讀得到的數字。 第 4 章給的是 ω(K),一條純古典的色散曲線,跟溫度毫無關係。要變成熱容,還缺兩塊:每個模式裡裝多少能量(Planck 分布,式 (2)),以及總共有多少模式、頻率怎麼分布(簡正模計數與態密度 D(ω))。把兩塊乘起來積分得到 U_lat,對 T 微分就是 C_lat。p.107 到 p.118 這十二頁,從頭到尾就是在鋪這一條流水線,Debye 模型與 Einstein 模型只是同一條流水線上兩種不同的 D(ω) 輸入。認清這個結構,你會發現這一章的公式比看起來少得多。

它為什麼被放在這個位置? 因為它是第 4 章唯一的出口。第 4 章幾乎通篇是古典力學(只有 p.99 半頁的 QUANTIZATION OF ELASTIC WAVES 是量子的),它交出頻率清單就下班了;量子統計要到第 5 章才接上去。同時它必須排在第 6 章之前——第 6 章的電子熱容整節都是拿來跟聲子熱容對照的,p.145 那條 C = γT + AT³ 把兩章的結果寫在同一條式子裡,而 C/T 對 T² 作圖分離兩項,是低溫量熱的標準手法。

但這一章裡有一道很多人沒注意到的接縫。 p.107 到 p.118 全部站在諧近似上:簡正模互不作用、頻率與溫度無關、熱容可以只靠「數模式」算出來。翻到 p.119,Kittel 當場把這個假設拆掉,並列出它的五個後果,其中兩條直接判了諧近似死刑——晶體不會熱膨脹,而且兩個晶格波互不作用,單一波永不衰減(於是平均自由程無限大、熱導率無限大)。所以本章後半不是前半的延伸,而是前半那套方法失效之後的補救:前半是計數問題,後半是碰撞問題。把兩半當成同一件事讀,就會在「熱容大所以導熱好」這類地方栽跟頭——Table 1 裡 Cu 的 θ = 343 K、κ = 4.01 W cm⁻¹K⁻¹,Si 的 θ = 645 K、κ = 1.48,排序完全相反(Table 1 在 p.116,該頁 PDF 無文字層,數值為算繪後目視判讀)。

還有一個容易被跳過的前提:p.107 第一句就聲明本章算的是 C_V,而實驗量的是 C_p(原文逐字:which is what the experiments determine),腳註給了 C_p − C_V = 9α²BVT。固體裡這個差通常很小,但它提醒你:教科書的熱容與你量到的熱容,嚴格說不是同一個量。

這一章如果只能帶走一條公式,帶走 p.118 的式 (37):D(ω) = [V/(2π)³]∮dS_ω/v_g。 Debye 模型只是它在「ω = vK、等頻面是球」這個特例下的結果;T³ 律也只是低頻 D(ω) ∝ ω² 的必然推論,而低頻的線性色散是真實晶體本來就有的性質,不是 Debye 硬塞的——所以 T³ 律在 T → 0 是漸近精確的,Debye 近似的粗糙不影響它。式 (37) 的推導從頭到尾只用了「K 空間態密度均勻」與「梯度垂直於等值面」兩件事,沒有任何一步用到「這是晶格振動」;Kittel 自己在 p.119 只用一句話帶過這個事實:We can use this result also in electron band theory. 分母的群速為零處態密度出現奇異(Van Hove 奇異點,同頁),這套機制在聲子譜上與在電子能帶上是同一件事。本章教的其實是一種「把 K 空間的幾何翻譯成能譜上的權重」的通用技術,熱容只是它的第一個應用。

它撐起後面哪幾章? 第 6 章拿同一條 K = ⅓Cvℓ 去做電子熱導與 Wiedemann–Franz 定律(p.156);第 18 章把本章的連續體圖像縮到奈米線,在 p.561–562 給出彈道熱導 G_th,並指出完美穿透時它只由基本常數與絕對溫度決定——那是本章 Fourier 定律的反面極限;第 19 章更是直接沿用,p.580 逐字寫著 As in Chapter 5, the expression for the thermal conductivity K has the form K = ⅓cvℓ,並用熔融石英的低溫熱容當作 Debye T³ 律的反例。這是全書處理晶格熱性質的唯一一章,二十五頁,後面三章談熱的地方全部回指這裡。

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

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

01聲子熱容:普朗克分布、簡正模計數、一維態密度PHONON HEAT CAPACITY — Planck distribution / Normal mode enumeration / Density of states in one dimension書本 p.105-110 · 概念 19 推導 9 公式 17把第 4 章算出的色散關係 ω(K) 換算成可以跟量熱計比對的晶格熱容 C_lat:需要「每個模式裝幾個聲子」(Planck 分布)與「總共有幾個模式」(簡正模計數與態密度)兩塊拼圖,這一段把兩塊都備齊。02聲子熱容:三維態密度與 Debye 模型PHONON HEAT CAPACITY — Density of states in three dimensions / Debye model for density of states書本 p.111-113 · 概念 17 推導 12 公式 15前一節已經把「每個聲子模態該貢獻多少能量」用 Planck 分布解決了,剩下的是一個純計數問題:三維晶體到底有多少個振動模態、頻率怎麼分布?這一段先用週期性邊界條件把 K 空間的格點數算乾淨(得到與任何色散關係無關的恆等式 D(ω)=(VK²/2π²)(dK/dω)),再用 Debye 的 ω=vK 把它具體化成 Vω²/2π²v³,並以「N 個原胞恰好用光 N 個聲學模態」定出截止頻率 ω_D 與 Debye 溫度 θ,最後積出一條只依賴 T/θ 的普適熱容曲線——高溫回到古典的 3Nk_B,低溫則預告下一節的 T³ 定律。03聲子熱容:Debye T³ 定律、Einstein 模型、D(ω) 的一般結果PHONON HEAT CAPACITY — Debye T3 law / Einstein model of the density of states / General result for D(ω)書本 p.114-119 · 概念 15 推導 10 公式 12Debye 積分在 T ≪ θ 時可以真的積出來——本段把 ∫₀^∞ x³/(eˣ−1)dx = π⁴/15 算到底,得到 C_V ≅ 234 Nk_B(T/θ)³;再用 Einstein 的單頻模型示範「如果晶體沒有低頻聲學模會發生什麼事」;最後把態密度從 Debye 的球面特例推廣成 D(ω) = (V/8π³)∮dS_ω/v_g,一條連電子能帶論都通用的公式。04非諧晶體交互作用與熱膨脹ANHARMONIC CRYSTAL INTERACTIONS / Thermal expansion書本 p.119-121 · 概念 22 推導 4 公式 6把前面整章默默使用的「諧近似」攤開來檢查:它到底許諾了哪五件事,這五件事在真實晶體怎麼一條條破功,以及用最小的一個模型(U = cx² − gx³ − fx⁴ 加上 Boltzmann 平均)算出熱膨脹係數從哪裡來。05熱導率與聲子氣的熱阻THERMAL CONDUCTIVITY / Thermal resistivity of phonon gas書本 p.121-124 · 概念 16 推導 4 公式 7熱容問「存得住多少熱」,這一段問「熱怎麼流過去」——先說明 Fourier 定律 j_U = −K dT/dx 這個形式本身就已經宣稱了熱傳是擴散而非彈道,再用氣體動力論推出 K = ⅓Cvℓ 並用它回推真實材料的聲子自由程,最後論證:光是把自由程限制住還不夠,必須有能同時重分配能量又鬆弛晶體動量的碰撞,否則熱阻是零。06Umklapp 過程、缺陷散射與習題Umklapp processes / Imperfections / PROBLEMS書本 p.125-129 · 概念 19 推導 14 公式 5K = ⅓Cvℓ 裡的 ℓ 到底由什麼決定——這一段先證明只有需要倒晶格向量 G 的 umklapp 碰撞才產生熱阻、再說明當 U 過程在低溫被 exp(−θ/2T) 凍結之後,接手的是試片尺寸與同位素亂序,最後用五道習題把零點能、維度與非諧性收束成 Grüneisen 常數。

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

  • 「Debye 溫度 θ 是材料常數。做熱分析要用的時候,從 Table 1 抄一個下來就對了——矽是 645 K,這是查得到的數字。」
    θ 不是一個數,是一族數。書本 p.116 Table 1 的表頭逐字是 Low temperature limit of θ, in Kelvin——它明白告訴你這是低溫極限值,是配 (32) 的 T³ 律用的那一個。真實晶體的 D(ω) 不是拋物線,不同溫區的熱容對聲子譜的不同部位加權,所以從高溫熱容擬合、從彈性常數計算、從熔點經驗式估計得到的 θ 彼此都不相等(矽:熱容低溫極限 645 K,由彈性常數算出的約 696 K——後者用的是手冊彈性常數,非書中數值)。更隱蔽的是約定問題:θ 的定義經過式 (23) 綁在「N 個原胞恰好用光 N 個聲學模態」上,但 Table 1 依慣例是用 N = 原子數(每原子 3 個模態)定的。矽是金剛石結構、一個原胞兩個原子(書 p.16–17),若改用原胞數去算,同一組彈性常數給出的 θ 會從約 696 K 掉到約 553 K(差一個 2^⅓ 因子)。兩套約定的 θ 差 26%,而 Debye 熱容對 θ 敏感——300 K 的矽 dlnC_V/dlnθ = −0.42(我算的),換算下來 C_V(300 K) 會差 8.7%(我算的)。抄一個數字進去,你其實抄進了一個沒寫出來的約定。
    引用 θ 之前一定要問三件事:(1) 它是低溫極限、高溫擬合,還是由彈性常數算的?(2) 分母的 N 是原子數還是原胞數?(3) 它要跟哪個公式配對使用?可操作的檢核:拿你手上的 θ 去算高溫極限,看它有沒有回到正確的 Dulong–Petit 值。以矽為例,用 θ = 645 K、N = 原子數、n = 8/a³ = 4.997×10²² cm⁻³(a = 5.430 Å,書 p.17)算出的古典上限是 3nk_B = 2.070 J cm⁻³K⁻¹(我算的),這個數必須與你用原子數密度直接算的一致;若你的 θ 配上的是原胞數,上限會只剩一半,馬上現形。另一個檢核:算 T/θ。同一個 300 K,對鑽石(θ = 2230 K)是 T/θ = 0.135、C 只有古典值的 17%;對鉛(θ = 105 K)是 T/θ = 2.86、已經 99% 飽和。同一個溫度、同一個元件裡的不同材料,可能一個早已飽和、另一個還深在量子區。(θ 與 κ 的數值取自 Table 1,該頁在 PDF 中無文字層,為算繪後目視判讀。)
  • 「熱容和熱導是同一套聲子物理的兩面。K = ⅓Cvℓ 就把它們綁在一起了,所以量到其中一個、再估個聲速,另一個就出得來。熱容大的材料導熱也好。」
    K = ⅓Cvℓ 確實把兩者寫在同一條式子裡,但那條式子裡有三個因子,而 C 與 ℓ 背後是完全不相干的兩套物理。C 只由 D(ω) 與 Planck 分布決定,是一個純計數問題:有多少模式、每個裝多少能量,不涉及任何碰撞。ℓ 則完全由碰撞決定——非諧耦合、Umklapp 過程、邊界、同位素亂序、雜質,這些東西對模式計數毫無影響。書裡就有現成的反例。Table 1(p.116)把 θ 與 300 K 熱導率印在同一格,但它們的排序完全相反:Cu 的 θ = 343 K、κ = 4.01 W cm⁻¹K⁻¹;Si 的 θ = 645 K、κ = 1.48。銅的晶格熱容並不比矽大,它熱導高是因為熱由電子搬運(第 6 章 p.156)。更強的反例在低溫:Table 2(p.122)顯示石英從 0 °C 降到 −190 °C 時,C 從 2.00 掉到 0.55 J cm⁻³K⁻¹(少了 3.6 倍),K 卻從 0.13 漲到 0.50 W cm⁻¹K⁻¹(多了 3.8 倍)。C 與 K 反向而行,差額全部由 ℓ 吸收(40 → 540 Å)。若 C 與 K 真是同一件事的兩面,這張表不可能長成這樣。
    把 C 與 K 當成兩個獨立的實驗問題,用兩套獨立的方法量、兩套獨立的模型解釋。可操作的判準是「對尺寸的敏感度」:C 是三維模態計數的結果,只要結構還是三維的、D(ω) 就不變,所以 C 對尺寸幾乎免疫;ℓ 卻會被幾何截斷(p.126:the thermal conductivity becomes a function of the dimensions of the specimen),所以 K 對尺寸極度敏感。這個不對稱直接決定了萃取策略:在奈米元件裡,C_th 是三個熱參數中最良好定義的一個,而 R_th 是最不良好定義的;τ = R_th·C_th 隨結構縮小而拉長,那個拉長幾乎全部來自 R_th。實務上:看到有人拿塊材 κ 除以塊材 C 得到 D 再套進奈米結構,要立刻警覺——分子被尺寸效應壓低、分母沒有,D 會被系統性高估。
  • 「非諧性是聲子散射的來源,所以非諧性就是熱阻的來源。三聲子碰撞越頻繁,熱導率越低。」
    這句話聽起來完全順理成章,而且 p.123 前半也確實是這樣鋪陳的(With anharmonic lattice interactions, there is a coupling between different phonons which limits the value of the mean free path)。但 p.123–124 整整兩頁的目的,就是把這個直覺打掉。Kittel 的原話很直接:It is not sufficient to have only a way of limiting the mean free path, but there must also be a way of establishing a local thermal equilibrium distribution of phonons. 然後他給出那個「subtle reason」:滿足 K₁ + K₂ = K₃ 的三聲子碰撞(N 過程)守恆聲子氣的總晶體動量 J(式 46),一團以漂移速度移動的平衡聲子分布,撞再多次也不會被減速——熱阻是零。Fig. 16a 的無摩擦管就是這件事的圖像:管裡的分子撞來撞去,熱阻率照樣是零。只有需要倒晶格向量的 umklapp 過程 K₁ + K₂ = K₃ + G(式 47,p.125)才會翻轉聲子流。p.126 最後一句把結論寫死:the phonon mean free path which enters (42) is the mean free path for umklapp collisions between phonons and not for all collisions between phonons. 而且這件事反過來也成立:靜態缺陷與晶界散射會改變動量,卻不改變聲子能量(p.123:the frequency of the scattered phonon is equal to the frequency of the incident phonon),所以它們能限制 ℓ、卻無法自己建立局域熱平衡。
    判準不是「有沒有碰撞」,是「總晶體動量 J 有沒有被改變」,外加「能量有沒有被重新分配」。 拿這兩把尺去量任何一個散射機制,就能立刻分類:N 過程——重分配能量、不改變 J → 不產生熱阻;U 過程——兩者都做 → 產生內稟體熱阻;邊界/缺陷/同位素——改變 J、不重分配能量 → 產生熱阻但不熱化,而且給出的是尺寸相依的 Casimir 型熱阻,不是材料常數。數值上的後果:U 過程需要波向量約 ½G 的高能聲子,其數目按 exp(−θ/2T) 凍結(p.126 逐字,注意分母那個 2 不能掉)。對矽(θ = 645 K)算一下:300 K 時 exp(θ/2T) = 2.93,77 K 時暴增到 65.9(我算的)——低溫下 U 過程幾乎關掉,此時你量到的熱阻幾乎全是界面與邊界的貢獻,改材料不會有用,改界面才會。
  • 「室溫遠高於大多數材料的 Debye 溫度,所以 300 K 的矽早就到 Dulong–Petit 了。做自熱分析時把 C_th 當成常數,誤差可以忽略。」

    先算一個數再說。Table 1 給矽 θ = 645 K,所以 300 K 對應 T/θ = 0.465——這連 θ 都沒到,怎麼會「遠高於」。對式 (30) 的 Debye 積分做數值積分(我算的),得到 300 K 時 C_V/3Nk_B = 0.802:室溫矽的晶格熱容只有古典飽和值的八成,而且還在往上爬。換成工程單位(用 a = 5.430 Å,書 p.17,得 n = 4.997×10²² cm⁻³、V_m = 12.052 cm³/mol):

    T (K)T/θC_V/3Nk_BC_vol (J cm⁻³K⁻¹)
    770.1190.1230.254
    3000.4650.8021.660
    3500.5430.8491.757
    4000.6200.8811.824
    古典上限 3nk_B—1.0002.070

    (θ 與 n 的來源如上;C 欄全部是我對書中式 (30) 做數值積分算出來的。1.660 與工程手冊常用的 1.63–1.66 一致,可交叉驗證。)300 → 400 K 之間 C_vol 上升 9.8%,而自熱造成的溫升正好就落在這個區間。這不是小數點後的問題。

    先算 T/θ,再決定能不能當常數。 可操作的門檻:要讓 C 走到古典值的 99% 以上,大約需要 T/θ ≳ 3;對矽那是 1900 K 以上,元件永遠到不了。所以矽的 C_th 在任何工作溫度都不是常數。更要命的是這個誤差的形狀。 它不是隨機誤差,而是與溫升相關的系統性誤差:功率越大、溫升越高、C_th 被低估得越多。在 Z_th(f) 或脈寬掃描的資料上,這會偽裝成「非線性熱阻」或「功率相依的時間常數」,讓你去追一個不存在的物理。做法:把 C_th(T) 保留成溫度的函數(用 θ = 645 K 的 Debye 積分即可,誤差在百分之幾以內),或至少在你的溫升區間取中點值而不是起點值。低溫量測更要小心——77 K 時 C_vol 只有室溫的 15%(0.254 vs 1.660),若用室溫值去反推熱阻,誤差是數倍而不是數個百分點。
  • 「Fourier 定律 j_U = −K dT/dx 只是熱導率的定義,定義當然永遠成立。所以不管結構多小,量到熱流和溫差就能定義一個有效熱導率,那就是材料性質。」
    p.121 整段就是在反駁「這只是個定義」。Kittel 寫完 (41) 之後沒有往下算,而是先審問這條式子:This form implies that the process of thermal energy transfer is a random process. 接著給反證法:If the energy were propagated directly through the specimen without deflection, then the expression for the thermal flux would not depend on the temperature gradient, but only on the difference in temperature ΔT between the ends of the specimen, regardless of the length of the specimen. 關鍵在於長度 L 出不出現:dT/dx ≈ ΔT/L,寫成梯度就等於宣稱「棒子加長一倍、通量減半」。這是一個可以實驗否證的物理宣稱,不是記號約定。當 ℓ ≳ L 時它會被否證,此時硬算出來的「有效 K」會隨試片尺寸漂移——那不是材料性質,是幾何與界面的性質。這件事書裡其實接了回來:第 18 章 p.561–562 給出完美穿透(純彈道)時的熱導 G_th,並指出它只由基本常數與絕對溫度決定,根本沒有材料參數的位置。
    判準是熱傳的 Knudsen 數 Kn ≡ ℓ/L(ℓ 是能量載子平均自由程,L 是熱流方向的特徵尺寸)。Kn ≪ 1 擴散、(41) 成立、R = L/(KA) 正比於 L;Kn ≫ 1 彈道、(41) 失效、熱阻趨近與 L 無關的界面熱阻;Kn ~ 1 兩者都不對,要解 Boltzmann 輸運方程。矽的 ℓ 是多少?書上沒有替矽算過,但書上示範過算法——p.580 就是拿熔融石英的 c、K、⟨v⟩ 反解出 ℓ ≈ 6×10⁻⁸ cm。照樣做:ℓ = 3κ/(C_vol·v),取 κ = 1.48 W cm⁻¹K⁻¹(Table 1)、C_vol = 1.660 J cm⁻³K⁻¹(我算的)、v = 5×10⁵ cm/s(Table 2 表題指定的代表聲速),得 ℓ_Si(300 K) ≈ 53 nm(我算的)。對照 FinFET 的鰭寬數奈米、閘長數十奈米——Kn ≳ 1,室溫就成立。p.126 逐字:This effect was discovered by de Haas and Biermasz.——他們當年要在低溫才看得到的尺寸效應(Fig. 18 的橫軸只到 100 K),在今天的元件裡是室溫常態。最好的診斷是尺寸序列:做一組只改長度(或膜厚)的結構,畫 R 對 L。過原點的直線 → 擴散;直線但有截距 → 擴散加上固定界面熱阻;在小 L 端彎曲趨於平台 → 你已經在彈道/過渡區,報出來的 K 不是材料常數。這一步該排在萃取流程的最前面,而不是當成補充討論。

④本章最重要的結論

聲子熱容的計算歸結為求態密度 D(ω):Debye 模型(聲速常數、截斷於 ω_D)給出高溫的古典值 3Nk_B 與低溫的 T³ 定律 C_V ≅ 234 Nk_B(T/θ)³,與實驗(如固態氬)吻合;Einstein 模型雖然低溫以指數衰減而不符實驗,但常用來近似光學聲子部分。真實晶體的熱膨脹、聲子間散射與有限熱傳導都無法由諧近似解釋,必須引入非諧項。熱傳導以 K = ⅓Cvℓ 描述,而能造成熱阻的是會改變聲子總動量的 umklapp 過程 K₁ + K₂ = K₃ + G,不是 N 過程;低溫下 umklapp 失效,平均自由程改由試樣尺寸與同位素、雜質等缺陷散射決定,此時 K 隨 C 而呈 T³ 變化。

書中逐字

"This function is called the density of modes or, more often, density of states." (p.108)

"These processes, discovered by Peierls, are called umklapp processes." (p.125)

"There is no thermal expansion." (p.119)

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

承 Ch4 的色散關係,接上 Planck 分佈得到聲子熱容;本章後半的熱導率與 Ch6 的電子熱導形成對照(並在 Ch6 p.156 由 Wiedemann–Franz 定律接起來)。Ch19 討論非晶低溫熱容時直接回指本章的 Debye T³ 律並指出偏離;Ch18 把本章的連續體圖像縮到奈米尺度。

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

直接

全書的熱容主場。Debye 模型、T³ 律、Dulong–Petit、Einstein 模型、非諧性與熱膨脹、熱導率 K = ⅓Cvℓ 全在這一章。必讀 p.107–118 與 p.121–128。

關聯程度:直接——這是全書的熱容主場,而且是唯一的一場

先把搜尋證據攤開。對全書 700 頁做逐頁文字檢索,下列詞彙出現次數皆為 0:self-heating、thermal time constant、thermal impedance、Kapitza、boundary resistance、thermal diffusivity、transient、Foster、Cauer。thermal resistance 全書只出現 1 次(p.127 Fig. 19 圖說,講 Casimir 的邊界散射)。所以 Z_th(f)、脈寬掃描、RC 熱網路、Foster/Cauer 拆解這整套語言,在 Kittel 的框架裡完全沒有對應物——以下凡屬元件層級的橋接,一律是【延伸知識】。

但它們需要的材料層輸入,本章全部給齊了,而且不只給了數字,還給了「什麼時候這些數字不能用」的判準。下面分四步走。


第一步:把書上的數字換成你要的 C_th(T)

書中直接可用的量只有兩個:θ = 645 K 與 κ(300 K) = 1.48 W cm⁻¹K⁻¹(Table 1,書 p.116;該頁在 PDF 中沒有文字層,數值為算繪放大後目視判讀)。原子數密度不必查表——書 p.17 逐字給了矽的晶格常數 a = 5.430 Å,金剛石結構每個慣用立方胞 8 個原子(書 p.16–17),所以 n = 8/a³ = 4.997×10²² cm⁻³(我算的),V_m = N_A/n = 12.052 cm³/mol(我算的)。

用 θ = 645 K 對書中式 (30) 的 Debye 積分做數值積分(以下 C 欄全部是我算的):

T (K)T/θC_V/3Nk_BC_mol (J mol⁻¹K⁻¹)C_vol (J cm⁻³K⁻¹)相對 300 K
770.1190.1233.060.25415%
3000.4650.80220.011.660100%
3500.5430.84921.171.757106%
4000.6200.88121.981.824110%
∞(Dulong–Petit)—1.00024.942.070125%

交叉驗證:1.660 J cm⁻³K⁻¹ 與工程手冊常用的 1.63–1.66 相符(手冊值為外部來源,非書中)。

直接後果:C_th 不是常數,300 → 400 K 上升 9.8%,而這正是自熱溫升的典型範圍。把 C_th 當常數擬合,會產生一個與功率相關的系統性偏差,在資料上偽裝成「非線性熱阻」。


第二步(最有價值的一步):本章的數字直接告訴你「ns 脈寬到底在量誰的熱容」

把書中的 κ 與上表的 C_vol 相除,得到熱擴散率(我算的):

D = κ / C_vol = 1.48 / 1.660 = 0.891 cm² s⁻¹(矽,300 K;手冊值約 0.88,相符)

時域的擴散長度 √(Dt) 與頻域的週期加熱穿透深度 √(D/πf)(兩欄都是我算的):

脈寬 t√(Dt)頻率 f√(D/πf)
100 ps94 nm1 GHz168 nm
1 ns299 nm100 MHz533 nm
10 ns944 nm10 MHz1.68 µm
100 ns2.99 µm1 MHz5.33 µm

現在把這個對上元件尺度。取一個示意的鰭:10 nm(寬)× 30 nm(高)× 100 nm(長)= 3.0×10⁻¹⁷ cm³。

結論很硬:在 ns 量級,你量到的 C_th 幾乎不含通道矽的資訊。 鰭在 10 ps 內就已經內部等溫了,ns 級的 Z_th 曲線描述的是 STI、埋氧層、基板與內連線。要「看見」鰭本身,脈寬得壓到 10 ps 以下。

反過來說,這也給了 Z_th(f) 一個乾淨的物理讀法:每一個頻率對應一個體積,穿透深度 √(D/πf) 就是那個體積的尺度,而該層的熱容就是「體積 × C_vol」。這是 Foster/Cauer 階梯網路能成立的物理根據,也是為什麼擬合出來的第一個 RC 級距不該被解讀成「通道熱容」——除非你的量測頻寬真的到了 10 GHz 以上。


第三步:本章的三個警告,每一個都會咬到萃取流程

(1)Fourier 定律可能在你的結構裡已經破了。 書上沒有替矽算過平均自由程,但書上示範過算法——p.580 就是拿熔融石英的 c、K、⟨v⟩ 反解 ℓ ≈ 6×10⁻⁸ cm。照做:ℓ = 3κ/(C_vol·v),取 κ = 1.48(Table 1)、C_vol = 1.660(我算的)、v = 5×10⁵ cm/s(Table 2 表題指定的代表聲速),得

ℓ_Si(300 K) ≈ 53 nm(我算的;巧的是這與 Table 2 裡石英 −190 °C 的 540 Å 幾乎同一個數量級)

對上鰭寬數奈米、閘長數十奈米,Kn = ℓ/L ≳ 1 在室溫就成立。p.121 的反證法在此完全適用:此時報出來的「有效熱導率」會隨結構尺寸漂移,那是幾何與界面的性質,不是材料常數。診斷方法是尺寸序列(R 對 L 作圖,見誤解五)。

(2)好消息:C_th 對尺寸效應的免疫力遠高於 R_th。 C 是模態計數的結果,只要結構還是三維的、D(ω) 就不變;ℓ 卻會被幾何截斷(p.126)。所以在奈米元件裡R_th 被強烈放大、C_th 幾乎不變,τ = R_th·C_th 的拉長幾乎全部記在 R_th 頭上。這對你是策略性的好消息:C_th 是三個熱參數裡最良好定義的一個,用它去反推「有效受熱體積」比用 R_th 去反推「有效熱導率」可靠得多。

(3)但 C_th 的免疫力有下限,而且下限在低溫。 式 (17) 的 K 空間格點間距是 2π/L。鰭寬 10 nm → 2π/L = 6.28×10⁸ m⁻¹,只比 K_D = (6π²n)^⅓ = 1.436×10¹⁰ m⁻¹ 小 22.8 倍(我算的),受限方向只剩約 18 個允許 K 點,D(ω) 的低頻端開始出現子帶結構。室溫下這個修正還小,因為 Debye 譜是「頭重腳輕」的:ω_D/2 以下的模態只佔 12.5%,0.8ω_D 以上佔 48.8%(我算的),室溫熱容由高頻模主導,而高頻模對邊界不敏感。但到 77 K 以下,熱容的貢獻整個移到低頻端,侷限效應就不能忽略了。所以室溫做 C_th 萃取比低溫安全。 第 18 章 p.561 正是處理這個尺度,p.561–562 給出一維聲子子帶的熱容與彈道熱導 G_th。

低溫還有第二個陷阱:77 K 時 C_vol 只剩室溫的 15%,同時 U 過程按 exp(−θ/2T) 凍結(矽:300 K 時 exp(θ/2T) = 2.93,77 K 時 65.9,我算的),ℓ 暴增、Kn 惡化。Table 2 已經示範過:石英的 ℓ 從 0 °C 的 40 Å 漲到 −190 °C 的 540 Å,13.5 倍(書中數字)。


第四步:哪些部分可以略過(誠實版)

書中段落對 C_th 萃取的用處建議
p.107–113 Planck 分布、態密度、Debye 模型直接產出 C_vol(T)必讀,而且要會自己算積分
p.116 Table 1θ 與 κ 的來源必讀(注意是低溫極限的 θ)
p.114–115 T³ 律、Fig. 9/10室溫元件用不到(T³ 只在 T ≲ θ/50 ≈ 13 K 乾淨)略讀,知道它的適用界線即可
p.117 Einstein 模型光學支的近似;矽的光學支室溫仍大半凍結略讀
p.119–121 非諧性與熱膨脹與 C_th 無關,除非你要做熱機械應力或 Grüneisen 分析可以跳過
p.121–122 Fourier 定律的前提、K = ⅓Cvℓ、Table 2給出 Kn 判準與 ℓ 的算法必讀
p.123–128 Umklapp、缺陷、同位素只影響 R_th,不影響 C_th若目標純粹是 C_th,讀結論即可

最後一條誠實的限定

Kittel 這二十五頁給的全部是塊材、單一材料、平衡態聲子的物理。真實的 Z_th(f) 裡還有本章完全沒碰的東西:串聯在散熱路徑上的 STI/埋氧/基板/內連線熱容、界面熱阻(Kapitza,全書 0 次)、以及電子—聲子能量交換不夠快造成的非平衡熱載子效應(全書未處理)。本章是骨架,不是全部。

但骨架裡的三樣東西——C_vol(T) 的 Debye 算法、D = κ/C_vol 給出的「頻率 ↔ 體積」對應、以及 Kn = ℓ/L 的合法性判準——是任何 C_th 萃取流程都繞不開的。先把這三樣釘死,再去談網路拆解。