Kittel·固態物理

全書逐章 / 第 5 章 聲子 II:熱性質 / 小節 5

05熱導率與聲子氣的熱阻

THERMAL CONDUCTIVITY / Thermal resistivity of phonon gas

書本 p.121-124 · 小節 5 / 6 · 概念 16 推導 4 公式 7

熱容問「存得住多少熱」,這一段問「熱怎麼流過去」——先說明 Fourier 定律 j_U = −K dT/dx 這個形式本身就已經宣稱了熱傳是擴散而非彈道,再用氣體動力論推出 K = ⅓Cvℓ 並用它回推真實材料的聲子自由程,最後論證:光是把自由程限制住還不夠,必須有能同時重分配能量又鬆弛晶體動量的碰撞,否則熱阻是零。

這一段在解決什麼問題,以及本文的引用體例我的判讀

第 5 章的前三分之二在算熱容:Planck 分布 (2)、模態計數、Debye 模型、Einstein 模型、D(ω)。熱容回答的是「這塊晶體存得住多少熱」。接著 p.119–121 的 ANHARMONIC CRYSTAL INTERACTIONS 講非諧性,回答的是「原子間的位能不是純拋物線會怎樣」,結果是熱膨脹。

到了 p.121 的 THERMAL CONDUCTIVITY,問題換成第三個、也是工程上最現實的一個:熱怎麼「流」過去。這一段要交出三樣東西:

第一,熱導率 K 的定義本身要成立的前提是什麼。Kittel 在寫下 j_U = −K dT/dx 之後,立刻花了一整段解釋「這個式子的形式本身就已經假設了某件事」——熱傳是隨機(擴散)過程,不是彈道直穿。這一段是整本書裡 ballistic 與 diffusive 熱傳對比的原始出處,對做 Cth(熱容/熱阻)萃取的人是最關鍵的接點,本文會逐句拆。

第二,K 的微觀表達式 K = ⅓Cvℓ 從哪來。書上用半頁的氣體動力論推完,但中間跳了好幾步:½ 從哪來、⅓ 從哪來、為什麼兩者最後不會互相干擾。本文把每一步補齊。

第三,什麼樣的散射機制才「算數」。這是 Thermal resistivity of phonon gas 這個子節的靈魂:光是把平均自由程 ℓ 限制住還不夠,還必須能建立局域熱平衡。Kittel 用一個非常漂亮的反例說明——三聲子 N 過程 K₁ + K₂ = K₃ 明明是不折不扣的碰撞,卻完全不產生熱阻。Fig. 16a–c 三張圖就是為了把這件事講死。

引用體例。本文中標「逐字」的英文句子,都取自書本該頁,並已用 PyMuPDF 把該頁算繪成影像目視核對過(這份 PDF 的文字層會把行間公式浮到段末、把下標順序打亂,不核對會出錯)。還原只做三件不改變任何文字的事:去掉跨行連字號(pro- ceed → proceed)、把浮出的行間公式與上下標放回紙面位置、把被拆成「1 3」的分數還原成 ⅓。除此之外一字不改。凡標「校註」的地方,都會說明我的判讀依據。

熱導率的定義:四個被寫進定義裡的限定詞

p.121 的開場只有一句話,但每個限定詞都有份量:

「The thermal conductivity coefficient K of a solid is defined with respect to the steady-state flow of heat down a long rod with a temperature gradient dT/dx:」

拆開來看:

「steady-state」(穩態)。K 不是一個瞬時量。定義它的實驗是:把一根棒子的兩端固定在兩個熱庫上,等到棒子裡每一點的溫度都不再隨時間變化,才去量熱流。這排除了所有暫態效應——也就是說,熱容 C 在這個定義裡被刻意排除掉了。∂T/∂t = 0,所以熱擴散方程 ∂T/∂t = (K/C)∂²T/∂x² 退化成 ∂²T/∂x² = 0,溫度沿棒子線性分布。這一點在實驗上很重要:量 K 的實驗與量 C 的實驗是兩種不同的實驗,而量熱擴散率 D = K/C 的實驗(如 3ω 法、時域熱反射 TDTR)又是第三種。

「a long rod」(長棒)。長,是為了讓一維近似成立(側向散熱可忽略),也是為了讓下一段要講的「擴散」假設有意義——棒子要比平均自由程長很多,才會有夠多次碰撞。

「with a temperature gradient dT/dx」(有溫度梯度)。注意 Kittel 說的是梯度(局域的、逐點定義的量),不是兩端的溫差。這個選擇本身就是一個物理宣稱,下一段會攤開來講。

「where j_U is the flux of thermal energy, or the energy transmitted across unit area per unit time.」j_U 的單位是 W cm⁻²(能量/面積/時間),所以 K 的單位是 W cm⁻¹ K⁻¹。下標 U 代表 energy(內能)。第 6 章 p.156 會再用 K = ⅓Cvℓ 這條式子處理金屬的電子熱導與 Wiedemann–Franz 定律,但該處只寫 K、K_el、σ,並沒有再用 j_U 這個記號;j_U 下一次出現要到第 8 章 p.215 定義 Peltier 係數 j_U = Π j_q 的地方。

(本區塊的英文引句逐字取自 p.121,並已用 PyMuPDF 算繪該頁目視核對;「四個限定詞」的拆解方式、以及穩態量測與熱擴散率量測(3ω 法、TDTR)的對照,是本講解為釐清定義前提所補的說明。)

熱導率的定義(Fourier 定律的一維形式)

j_U = −K dT/dx

(41)  p.121

j_U:熱能通量,即單位時間內通過單位面積傳輸的能量,單位 W cm⁻²(原文:「the flux of thermal energy, or the energy transmitted across unit area per unit time」);K:熱導率係數,單位 W cm⁻¹ K⁻¹;dT/dx:溫度梯度,單位 K cm⁻¹。負號表示熱從高溫流向低溫——梯度為正(右邊熱)時通量向左,j_U 為負。

出處推導
這不是從任何更基本的式子推來的,它是 K 的定義式,前提由 p.121 明白列出:(i) 穩態(steady-state)、(ii) 長棒(long rod)、(iii) 存在良好定義的局域溫度梯度 dT/dx。第三個前提最強:它預設晶體內部每一點都能指定一個局域溫度,也就是預設了局域熱平衡存在。p.123 整段就是在追問「什麼機制保證這件事」。
極限行為
ℓ ≫ L(彈道極限)時這個形式失效。原文對此講得很直接:如果能量不經偏折直穿,「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」——通量只看兩端溫差 ΔT,與試片長度無關。此時寫成 j_U ∝ ΔT 才對,硬套 (41) 會得到一個隨長度增大的「表觀 K」,那不是材料性質。ℓ ≪ L(擴散極限)時 (41) 成立,K 才是材料的本徵性質。三維推廣是 j_U = −K·∇T,各向異性晶體的 K 是二階張量——Table 2 的石英特別註明「Parallel to optic axis」正是因為這個緣故。
典型數值
0 °C(273 K,接近室溫)的石英(平行光軸)K = 0.13 W cm⁻¹ K⁻¹(Table 2, p.122)。若在 1 cm 長的棒上加 100 K 溫差,dT/dx = 100 K cm⁻¹,則 j_U = 13 W cm⁻²。

「這個形式本身就蘊含隨機過程」:原文的完整論證

這是 p.121 最該讀慢的一段。Kittel 寫完 (41) 之後沒有直接跳去算 K,而是先回頭審問這條式子:

「This form implies that the process of thermal energy transfer is a random process.」

(這個形式蘊含了熱能傳遞是一個隨機過程。)

然後他解釋什麼叫「不是隨機過程」:

「The energy does not simply enter one end of the specimen and proceed directly (ballistically) in a straight path to the other end, but diffuses through the specimen, suffering frequent collisions.」

(能量並不是單純地從試片一端進入、直線彈道式地直達另一端,而是在試片中擴散,經歷頻繁的碰撞。)

接著是反證法的核心:

「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.」

(如果能量毫無偏折地直接穿過試片,那麼熱通量的表達式將不依賴於溫度梯度,而只依賴於試片兩端的溫差 ΔT,與試片長度無關。)

最後收尾:

「The random nature of the conductivity process brings the temperature gradient and, as we shall see, a mean free path into the expression for the thermal flux.」

(傳導過程的隨機性,把溫度梯度、以及我們稍後會看到的平均自由程,帶進了熱通量的表達式。)

這段論證的邏輯結構是這樣的。梯度 dT/dx 與溫差 ΔT 的差別,全部在於長度 L 出不出現:dT/dx ≈ ΔT/L。如果通量正比於梯度,那麼同樣的兩端溫差下,棒子加長一倍、通量減半;如果通量只正比於 ΔT,那麼加長棒子完全不影響通量。這是一個可以直接做的實驗:拿同一材料、不同長度的樣品,固定兩端溫差量熱流。

擴散輸運下,熱流會隨長度下降(因為中間的碰撞次數 ∝ L,每次碰撞都在「阻擋」);彈道輸運下,熱流與長度無關(沒有碰撞,中間那段長度什麼事也沒做)。所以「j_U 依賴梯度」這件事不是數學上的方便,而是一個關於微觀機制的實質宣稱:能量載子在棒內走的是隨機漫步,不是直線。

換句話說:Fourier 定律的合法性,等價於系統處於擴散區。這句話是整段的精髓,也是後面 Table 2 為什麼要去回推 ℓ 的理由——ℓ 與 L 誰大誰小,決定了你寫下的 K 到底是不是材料性質。

(四段英文引句全部逐字取自 p.121 並經算繪目視核對;「梯度與溫差的差別全在長度出不出現」這個形式化說法,以及據此設計長度序列實驗的建議,是本講解補的解讀。)

Ballistic 對 Diffusive:把 p.121 的論證變成可操作的判準我的判讀

上一段是 Kittel 的原話。這一段把它翻成做熱量測時真正用得上的形式。這是本節對 Cth/熱阻萃取研究最重要的接點,值得寫長一點。

判準只有一個無因次數。令 ℓ 為能量載子(在介電固體裡就是聲子)的平均自由程,L 為熱流方向上的特徵尺寸(樣品長度、薄膜厚度、FinFET 的鰭高、接觸孔的孔徑),定義

Kn ≡ ℓ / L (熱傳的 Knudsen 數)

  • Kn ≪ 1:擴散區。載子在跨越試片的過程中碰撞很多次,局域熱平衡成立,(41) 成立,K 是材料常數,熱阻 R = L/(KA) 與長度成正比。
  • Kn ≫ 1:彈道區。載子從一端直穿到另一端都不碰撞,(41) 失效。熱阻不再與 L 成正比,而是趨近一個與 L 無關的界面/接觸熱阻。此時若還硬要報一個「有效熱導率」K_eff,它會隨 L 變小而變小(因為 R 不縮小、L 縮小,K_eff = L/(RA) 跟著縮小),量出來的數字不是材料性質,是幾何與界面的性質。
  • Kn ~ 1:過渡區。兩者都不對,要解 Boltzmann 輸運方程(BTE)或用 Casimir/Fuchs–Sondheimer 這類修正。

為什麼這件事在奈米元件上會炸開。Table 2 的數字馬上就告訴你尺度:0 °C(273 K,接近室溫)的 NaCl ℓ 約 23 Å,石英約 40 Å——這種尺度下,任何微米級的樣品都穩穩在擴散區,Kittel 寫 (41) 完全安全。但

  • 矽在室溫的聲子平均自由程分布極寬,大部分熱是由 ℓ 遠大於 100 nm 的長波聲子攜帶的;
  • FinFET 的鰭寬只有數奈米、閘極長度數十奈米;
  • 而 Table 2 也已經示範了:降溫會讓 ℓ 暴增(石英從 40 Å 變成 540 Å,約 13.5 倍;NaCl 從 23 Å 變成書上印的 100 Å,約 4.3 倍——不過這一筆與 (44) 對不回來,用書上自己的 C、K 重算是 162 Å、約 7 倍,詳見後文驗算)。低溫量測會把你從擴散區推進過渡甚至彈道區。

所以在奈米元件的熱萃取裡,「先驗證自己在哪一區」是第一步,不是補充說明。如果你的 Cth/Rth 萃取流程預設了 Fourier 定律(絕大多數緊緻模型、RC 熱網路、有限元熱模擬都預設了),那麼 Kn ≳ 1 時萃取出來的「熱導率」會隨量測結構的尺寸漂移,而這個漂移不是量測誤差,是模型用錯了。

尺寸相依性是最好的診斷。這正是 p.121 那段反證法的直接應用:做一組只改長度(或膜厚)的樣品,畫 R 對 L 的圖。

  • 通過原點的直線 → 擴散,斜率給 1/(KA)。
  • 直線但截距不為零 → 擴散加上一個固定的界面熱阻。
  • 曲線、或 R 在小 L 時明顯偏離線性而趨於平台 → 你進入彈道/過渡區了。

Kittel 自己在第 18 章接了回來。p.562(第 18 章 Nanostructures)討論奈米線的熱輸運時寫道:「Note that in the case of perfectly transmitted phonons through the channel, the thermal conductance (78) is determined only by fundamental constants and the absolute temperature.」——完美穿透(即純彈道)時,熱導只由基本常數與絕對溫度決定,也就是熱導量子。那是 Kn → ∞ 的極端;p.121 這一段是 Kn → 0 的極端。兩端之間的整片灰色地帶,就是今天奈米尺度熱輸運研究的主戰場。

(本段的 Knudsen 數判準、尺寸相依性診斷與元件尺度的討論是本講解補充;Kittel 在 p.121 只給了「彈道時通量與長度無關」這個原則性的對比。)

氣體動力論的熱導率公式

K = ⅓ C v ℓ

(42)  p.122

C:單位體積的熱容(原文:「the heat capacity per unit volume」),單位 J cm⁻³ K⁻¹;v:粒子平均速度(「the average particle velocity」);ℓ:粒子兩次碰撞之間的平均自由程(「the mean free path of a particle between collisions」)。注意 C 是體積熱容不是莫耳熱容,也不是每顆粒子的熱容——每顆粒子的熱容在推導裡寫作小寫 c,兩者的關係是 C ≡ nc。

出處推導
由 (43)→(44) 的初等氣體動力論推導得到,再與定義式 (41) 比對係數。原文明說這是「From the kinetic theory of gases」,而且是先給結果再補推導(「We give the elementary kinetic theory which leads to (42).」)。推導中用到三個假設:(i) 粒子以平均自由程 ℓ = vτ 為尺度隨機漫步;(ii) 速度分布各向同性,因而 ⟨v_x²⟩ = ⅓⟨v²⟩;(iii) 速率為定值(「If, as for phonons, v is constant」),這一條對聲子在 Debye 近似下成立,對真實氣體只是近似。
極限行為
這是一條「單一弛豫時間、灰體」的公式,它把所有模態塞進同一個 v 與同一個 ℓ。真實情況要寫成對模態的積分 K = ⅓∫C(ω)v(ω)ℓ(ω)dω,不同頻率的聲子有完全不同的自由程;把它縮成單一 ℓ,得到的是一個熱流加權平均,這就是 Table 2 那些 ℓ 的真正身分。高溫極限(T ≫ θ_D):C 趨於 Dulong–Petit 的定值 3nk_B,而 ℓ ∝ 1/T(p.123),所以 K ∝ 1/T。低溫極限(T ≪ θ_D):C ∝ T³(Debye),而 ℓ 因 umklapp 過程凍結而暴增,最後被晶界尺寸截斷(ℓ → 樣品尺寸,即 Casimir 極限),於是 K ∝ T³ 隨溫度下降而下降。兩個趨勢相反,中間必然有一個熱導率極大值——這正是純淨介電晶體在數十 K 附近會出現 K 峰的原因。K → ∞ 的情形:ℓ → ∞(純諧振晶體、或只有 N 過程)時 K 發散,熱阻為零,這是 p.123–124 整段的主題。
典型數值
用 Table 2 的石英 0 °C 資料反過來檢查:C = 2.00 J cm⁻³K⁻¹、v = 5×10⁵ cm/s、ℓ = 40 Å = 4×10⁻⁷ cm,得 K = ⅓ × 2.00 × 5×10⁵ × 4×10⁻⁷ = 0.133 W cm⁻¹K⁻¹,與表上的 0.13 相符。

補完 (43):從 ½n⟨|v_x|⟩ 到 j_U = −n⟨v_x²⟩cτ dT/dx補推導

推導前要先擺正三件事。第一,這是氣體動力論,Kittel 是先在分子氣體的語言裡把式子推出來,再整批搬到聲子氣。第二,全程用的是單一弛豫時間近似:每顆粒子走一段 ℓ = vτ 就徹底忘掉來歷,跟當地重新熱平衡——這一步偷偷用掉了「局域熱平衡存在」這個假設,而 p.123 整節就是在檢討這個假設什麼時候真的成立。第三,全程是線性響應:泰勒展開只留一階,所以只對「ΔT 遠小於 T 本身」的小梯度成立。

推導後最值得記住的是係數的來歷。很多人記得 K = ⅓Cvℓ 卻說不清 ⅓ 哪來,因為推導中確實同時出現了 ½(單向粒子流)和 2(自由程兩端相距 2ℓ_x),兩者相乘為 1,於是唯一倖存的數值因子是來自 ⟨v_x²⟩ = ⅓⟨v²⟩ 的 ⅓。這個 ⅓ 是維度數的倒數:一維鏈裡不會有這個 ⅓,二維膜裡會是 ½。當你在文獻上看到 K = ½Cvℓ 卻沒有錯,通常就是因為那是二維系統(或者作者用了不同的自由程定義)。

另一條等價的記帳法。有些教科書不走「兩端相距 2ℓ_x」的路,而是說:向右的粒子帶來 x − ℓ_x 的能量、向左的帶走 x + ℓ_x 的能量,兩股各 ½n⟨|v_x|⟩,一步到位。也有人乾脆寫「單向通量 n⟨|v_x|⟩ 乘上一個自由程的溫降」,跳過 ½ 和 2。三種寫法算出來一樣,但只有把 ½ 和 2 都寫明,才看得出 ⅓ 是幾何來的、不是碰巧湊出來的。

  1. 設定(原文逐字):「The flux of particles in the x direction is ½n⟨|v_x|⟩, where n is the concentration of molecules; in equilibrium there is a flux of equal magnitude in the opposite direction. The ⟨···⟩ denote average value.」也就是:n 是粒子數密度,向 +x 走的粒子通量是 ½n⟨|v_x|⟩,向 −x 走的也是 ½n⟨|v_x|⟩。
  2. 為什麼是 ½。平衡時速度分布對 v_x → −v_x 對稱,所以恰好一半的粒子 v_x > 0、一半 v_x < 0。單位時間內穿過單位面積、由左往右的粒子數,就是「v_x > 0 的那一半密度」乘上「這一半的平均 |v_x|」= (n/2)·⟨|v_x|⟩。淨通量為零(因為兩個方向抵消),但能量淨通量不為零,因為兩個方向來的粒子攜帶的能量不同。
  3. 每顆粒子交出多少能量(原文逐字):「If c is the heat capacity of a particle, then in moving from a region at local temperature T + ΔT to a region at local temperature T a particle will give up energy c ΔT.」即單顆粒子從 T + ΔT 處移動到 T 處,交出 cΔT 的能量。這裡的 c 是單顆粒子的熱容(小寫),不是 (42) 裡的體積熱容 C。
  4. ΔT 是誰的溫差(原文式子,p.122 有排版但未編號):ΔT = (dT/dx)·ℓ_x = (dT/dx)·v_x τ,其中「τ is the average time between collisions」。關鍵在於這個 ΔT 是一條自由程兩端的溫差(「ΔT between the ends of a free path of the particle」),而自由程在 x 方向的投影是 ℓ_x = v_x τ。
  5. 把兩個方向合起來(這是書上只用一句「from both senses of the particle flux」帶過、最容易卡住的一步)。向 +x 走的粒子,是從 x − ℓ_x 那個較「上游」的位置出發、帶著那裡的溫度資訊過來的;向 −x 走的粒子,帶的是 x + ℓ_x 的資訊。兩股粒子流的大小相同(都是 ½n⟨|v_x|⟩),但攜帶的能量差了 c·[T(x−ℓ_x) − T(x+ℓ_x)]。
  6. 把溫度做泰勒展開:T(x−ℓ_x) − T(x+ℓ_x) = −2ℓ_x (dT/dx) + O(ℓ_x³)。注意這裡出現一個因子 2——它來自「兩端相距 2ℓ_x」,也就是原文所說的 both senses。
  7. 組裝:j_U = ½n⟨|v_x|⟩ · c · [T(x−ℓ_x) − T(x+ℓ_x)] = ½n⟨|v_x|⟩ · c · (−2 ℓ_x dT/dx) = −n⟨|v_x| ℓ_x⟩ c dT/dx。½ 與 2 在這裡相乘抵消,只剩下 1。這就是為什麼最後的答案裡看不到 ½。
  8. 把 ℓ_x 換成 v_x τ:|v_x| ℓ_x = |v_x| · |v_x| τ = v_x² τ(取絕對值是因為兩個方向都算過了,符號已由上一步的溫差處理掉)。於是 j_U = −n⟨v_x²⟩ c τ dT/dx,正是 (43) 的第一個等號。
  9. 各向同性把 x 換成全速率:v² = v_x² + v_y² + v_z²,各向同性下 ⟨v_x²⟩ = ⟨v_y²⟩ = ⟨v_z²⟩,故 3⟨v_x²⟩ = ⟨v²⟩,即 ⟨v_x²⟩ = ⅓⟨v²⟩。代回得 j_U = −⅓ n⟨v²⟩ c τ dT/dx,這是 (43) 的第二個等號。⅓ 就是從這裡來的,而且只從這裡來。

動力論的能量淨通量

j_U = −n⟨v_x²⟩ c τ dT/dx = −⅓ n⟨v²⟩ c τ dT/dx

(43)  p.122

n:粒子數密度(cm⁻³);c:單顆粒子的熱容(J K⁻¹);τ:兩次碰撞之間的平均時間(「the average time between collisions」);⟨v_x²⟩:速度 x 分量平方的平均;⟨v²⟩:速率平方的平均。原文引言:「The net flux of energy (from both senses of the particle flux) is therefore」。

出處推導
由 ½n⟨|v_x|⟩(x 方向單向粒子通量)、cΔT(單顆粒子跨越一條自由程所交出的能量)、ΔT = (dT/dx)ℓ_x = (dT/dx)v_xτ(自由程兩端的溫差)三者組合而成。第二個等號用的是速度分布各向同性:⟨v_x²⟩ = ⅓⟨v²⟩。書上沒寫中間的組裝步驟,本文前一個推導區塊補齊。
極限行為
梯度要小:泰勒展開只保留一階,所以要求一條自由程上的溫降遠小於當地溫度,ℓ|dT/dx| ≪ T。這個條件在強脈衝加熱或次微米熱點附近會被破壞。單一 τ:這裡假設所有粒子共用一個弛豫時間;真實聲子氣的 τ 強烈依賴頻率與極化,要做 K = ⅓∫C(ω)v(ω)²τ(ω)dω。注意 ⟨v²⟩ ≠ ⟨v⟩²:(43) 裡站著的是 ⟨v²⟩,只有在下一步假設 v 為定值時才能寫成 v²,這正是原文接下來那句「If, as for phonons, v is constant」的用意。對馬克士威分布的真實氣體,⟨v²⟩/⟨v⟩² = 3π/8 ≈ 1.18,所以 ⅓Cvℓ 對氣體只是量級估計。

聲子氣的熱通量(v 為定值時的 (43))

j_U = −⅓ C v ℓ dT/dx

(44)  p.122

C ≡ nc:單位體積熱容;ℓ ≡ vτ:平均自由程。原文逐字:「with ℓ ≡ vτ and C ≡ nc. Thus K = ⅓Cvℓ.」

出處推導
由 (43) 在「v 為定值」的條件下改寫。原文逐字:「If, as for phonons, v is constant, we may write (43) as」。代數上是:−⅓n⟨v²⟩cτ = −⅓(nc)(v²τ) = −⅓C·v·(vτ) = −⅓Cvℓ。把 v²τ 拆成 v×(vτ) 這一步就是 ℓ ≡ vτ 的定義所在。
極限行為
「v is constant」對聲子為什麼合理:在 Debye 近似下所有聲學模的色散是 ω = vK,群速度等於相速度等於聲速,與波向量無關,因此 ⟨v²⟩ = v² 嚴格成立。但這個近似在布里淵區邊界附近會崩潰——第 4 章的單原子鏈色散告訴我們,K → π/a 時 dω/dK → 0,群速度趨近於零。高溫下被大量激發的正是這些區邊界附近的短波聲子,它們幾乎不輸運能量。光學支更極端:光學支多半色散很弱、群速度遠小於聲學支,所以光學聲子對 K 的直接貢獻通常很小(它們對 K 的主要影響是間接的——當成散射通道,吃掉聲學聲子的自由程)。這是 (44) 這種「單一 v」寫法最大的系統性偏差來源。與 (41) 比對即得 (42) K = ⅓Cvℓ,這就是整個推導的目的。
典型數值
NaCl 0 °C:C = 1.88 J cm⁻³K⁻¹、v = 5×10⁵ cm/s、ℓ ≈ 23 Å,得 j_U/(dT/dx) = −⅓(1.88)(5×10⁵)(2.3×10⁻⁷) ≈ −0.072 W cm⁻¹K⁻¹,與表上的 K = 0.07 一致。

Debye 的搬家:從分子氣體到聲子氣體

(42) 原本是一條氣體的公式,Kittel 在 p.122 交代了它怎麼變成固體的公式:

「This result was applied first by Debye to describe thermal conductivity in dielectric solids, with C as the heat capacity of the phonons, v the phonon velocity, and ℓ the phonon mean free path.」

(這個結果最早由 Debye 應用於描述介電固體中的熱導率,其中 C 取聲子的熱容、v 取聲子速度、ℓ 取聲子平均自由程。)

這一句話很短,但做了三次概念替換,每一次都值得停下來想:

n(分子數密度)→ 聲子數密度。分子數是守恆的,聲子數不守恆——升溫會憑空產生聲子,降溫會湮滅聲子。這個差別在 (42) 裡被藏起來了,因為 (42) 裡根本沒有 n,只有 C = nc。把 n 和 c 綁成 C,正好繞開了「聲子數不守恆」這個麻煩:不管有幾顆聲子、每顆熱容多少,你要的只是總體積熱容,而那是第 5 章前半已經用 Debye 模型算好的東西。這是一個非常聰明的迴避。

v(分子熱運動速率)→ 聲速。分子氣體的 v 是熱速度,正比於 √(T/m),會隨溫度變;聲速幾乎不隨溫度變(只透過非諧性與熱膨脹有微弱的溫度依賴)。這是聲子氣與分子氣最大的不同之一,也是 Table 2 可以對兩個相差 190 K 的溫度用同一個 v = 5×10⁵ cm/sec 的理由。

ℓ(分子間碰撞的自由程)→ 聲子平均自由程。分子的碰撞是實打實的兩體碰撞,機制不需要解釋;聲子在諧振近似下根本不碰撞——簡正模是嚴格的本徵態,各自獨立演化到天荒地老。所以「聲子的平均自由程」這個概念本身就需要非諧性才能存在。這正是下一節 Thermal resistivity of phonon gas 的開場白要處理的事。

適用範圍的提醒:Kittel 說的是「dielectric solids」(介電固體)。金屬裡熱主要由電子攜帶,聲子熱導只是小配角,那要等到第 6 章的 Wiedemann–Franz 定律才處理。

(引句逐字取自 p.122;三次「概念替換」的拆法與其後的討論是本講解補充,Kittel 在此處只寫了那一句話。)

Table 2(p.122)原表重現

表題:Table 2 Phonon mean free paths

表下說明(逐字):「[Calculated from (44), taking v = 5 × 10⁵ cm/sec as a representative sound velocity. The ℓ's obtained in this way refer to umklapp processes.]」

(由 (44) 計算,取 v = 5 × 10⁵ cm/sec 作為代表性的聲速。以此方式得到的 ℓ 指的是 umklapp 過程。)

CrystalT, °CC, in J cm⁻³K⁻¹K, in W cm⁻¹K⁻¹ℓ, in Å
Quartz ᵃ02.000.1340
−1900.550.50540
NaCl01.880.0723
−1901.000.27100

ᵃ Parallel to optic axis.(平行於光軸。)

讀這張表要抓住三件事。

第一,ℓ 不是量出來的,是算出來的。表題說明白了「Calculated from (44)」——實驗量到的是 C 和 K,ℓ 是用 ℓ = 3K/(Cv) 反推的。所以這一欄的可信度完全建立在 (44) 這個模型上;如果模型不對(例如樣品其實在彈道區),這個 ℓ 就沒有物理意義。

第二,所有四筆都用同一個 v = 5×10⁵ cm/sec。這是一個粗糙但合理的代表值(典型固體聲速 3–7 km/s = 3–7 ×10⁵ cm/s)。它同時被用在石英與 NaCl、也被用在 0 °C 與 −190 °C——後者尤其重要,因為聲速幾乎不隨溫度變,所以這個做法是站得住的。但這也意味著表中的 ℓ 帶有一個共同的系統性不確定度:v 差 30%,ℓ 就差 30%。

第三,這些 ℓ 指的是 umklapp 過程。這句註記(「The ℓ's obtained in this way refer to umklapp processes.」)是一個很強的物理宣稱,而且它是前瞻的——p.122 講這句話的時候,umklapp 這個詞在書上還沒有定義;在此之前它只在第 5 章章首目錄(p.105)以「Umklapp processes 125」的形式出現過,正式定義要到 p.125 的 Umklapp Processes 小節。它的意思是:在這兩個溫度、這兩種(相對純淨、相對大塊的)樣品裡,限制熱流的主要機制既不是晶界、也不是雜質、更不是 N 過程,而是 U 過程。這個宣稱的理由,正是 p.123 到 p.125 這三頁要論證的東西:只有 U 過程能讓聲子氣的動量鬆弛,因此只有 U 過程能產生體熱阻。

(校註:兩個數值上的觀察)表中「Quartz」與表下註記「Parallel to optic axis」——石英是三方晶系,熱導率是張量,所以必須註明方向;表中沒有給垂直光軸的數值。另一個觀察與數字有關,放在下一個區塊處理。

(表格、表題說明與註腳逐字重現自 p.122,並已算繪該頁目視核對欄位對應;「讀這張表要抓住三件事」的歸納與其後的解讀是本講解補充。)

親手回推 Table 2:量綱檢查與四筆數字的逐一驗算補推導

為什麼要自己算一遍。因為 Table 2 的 ℓ 欄本來就是算出來的,不是量出來的——表題白紙黑字寫著「Calculated from (44)」。既然如此,讀者就有權(也有義務)把這個計算重跑一次。這件事只需要一台計算機,卻能同時確認三件事:你有沒有搞懂 C 是體積熱容、你的單位換算對不對、以及書上的數字彼此自洽不自洽。

量綱那一步不要跳過。K = ⅓Cvℓ 這條式子最常見的錯誤,是把 C 誤當成莫耳熱容(J mol⁻¹K⁻¹)或比熱(J g⁻¹K⁻¹)。這兩種誤用在量綱檢查時會立刻現形:用莫耳熱容算出來的「長度」會帶著 mol⁻¹,根本不是長度。所以每次代數值之前先把單位串一遍,是最便宜的保險。

校註:NaCl −190 °C 這一筆對不回來。我用 PyMuPDF 把書本 p.122 算繪成影像放大目視確認過,表上印的確實是 1.00 / 0.27 / 100,不是文字層抽取錯位。四筆裡有三筆能在 3% 以內重現,唯獨這一筆差了 62%。就本頁所印的數字與表題給的 v 而言,(44) 給出的是 162 Å 而不是 100 Å。本文不去猜測原因(可能是 K 值取自不同來源、可能是印刷誤植、也可能是原始文獻對這一筆用了不同的 v),只忠實報告:依表題的算法與表內的 C、K,這一筆算出 162 Å。

這個發現本身就是一堂課:教科書上的數字表不一定內部自洽,而唯一能發現這件事的方法,就是自己算一遍。若在研究中要引用 Table 2 的 ℓ,建議用 C 與 K 自行重算,而不要直接抄 ℓ 欄。

順便看一個趨勢。把四筆的 ℓ 排開:石英 40 → 540 Å(13.5 倍),NaCl 23 → 162 Å(依重算值,7 倍;依書上印的 100 Å 則是 4.3 倍)。不論用哪一組數字,結論都一樣:從 0 °C 降到 −190 °C,聲子平均自由程增加了將近一個數量級以上。這個現象要解釋,見下一節。

  1. 反解公式:由 K = ⅓Cvℓ 得 ℓ = 3K / (C v)。
  2. 量綱檢查(一定要先做):[K]/([C][v]) = (W cm⁻¹K⁻¹) ÷ (J cm⁻³K⁻¹ × cm s⁻¹) = (W cm⁻¹K⁻¹) ÷ (J cm⁻²K⁻¹s⁻¹) = (W cm⁻¹K⁻¹) × (cm² K s / J) = W·cm·s / J。因為 W·s = J,所以整串化為 cm。量綱正確,答案是長度。
  3. 單位換算備忘:1 cm = 10⁸ Å。所以算出來的 cm 值乘 10⁸ 就是 Å。
  4. 石英 0 °C:ℓ = 3 × 0.13 / (2.00 × 5×10⁵) = 0.39 / 1.00×10⁶ = 3.9×10⁻⁷ cm = 39 Å。表上印 40 Å。✔ 相符(四捨五入)。
  5. 石英 −190 °C:ℓ = 3 × 0.50 / (0.55 × 5×10⁵) = 1.50 / 2.75×10⁵ = 5.45×10⁻⁶ cm = 545 Å。表上印 540 Å。✔ 相符。
  6. NaCl 0 °C:ℓ = 3 × 0.07 / (1.88 × 5×10⁵) = 0.21 / 9.40×10⁵ = 2.23×10⁻⁷ cm = 22.3 Å。表上印 23 Å。✔ 相符(K 只給到兩位小數,若原始值是 0.072 就正好得 23.0 Å)。
  7. NaCl −190 °C:ℓ = 3 × 0.27 / (1.00 × 5×10⁵) = 0.81 / 5.00×10⁵ = 1.62×10⁻⁶ cm = 162 Å。表上印 100 Å。✘ 對不上。
  8. 對不上的程度:162 vs 100,差 62%,遠超過四捨五入能解釋的範圍(其他三筆的偏差都在 3% 以內)。反推需要什麼條件才能得到 100 Å:若 C = 1.00、v = 5×10⁵ 不動,則需要 K ≈ 0.167 W cm⁻¹K⁻¹;若 K = 0.27、v 不動,則需要 C ≈ 1.62 J cm⁻³K⁻¹;若 C、K 都不動,則需要 v ≈ 8.1×10⁵ cm/s。三者都與表上印的數字或表題說明不符。

為什麼低溫下 ℓ 會暴增:兩個相反的趨勢在打架我的判讀

Table 2 最值得玩味的,是 C、K、ℓ 三欄的走向彼此矛盾:

量0 °C → −190 °C 的變化原因
C(體積熱容)下降(石英 2.00 → 0.55;NaCl 1.88 → 1.00)高頻模態被凍結,往 Debye 的 T³ 走
K(熱導率)上升(石英 0.13 → 0.50;NaCl 0.07 → 0.27)與 C 相反
ℓ(自由程)暴增(一個數量級上下)必須夠大,才能同時解釋上面兩列

邏輯是這樣的。K = ⅓Cvℓ,v 基本不變。C 掉了 3.6 倍(石英),而 K 反而漲了 3.8 倍,所以 ℓ 必須漲 (2.00/0.55) × (0.50/0.13) ≈ 14.0 倍才收得平;書上印的是 540/40 = 13.5,略低的原因是 ℓ 欄四捨五入(由 (44) 重算的 ℓ 比值是 545/39 ≈ 14.0)。換句話說,表裡 ℓ 的暴增,不是獨立的實驗事實,而是「C 下降但 K 上升」這個觀測事實在 (44) 框架下的必然推論。

物理上為什麼 ℓ 會暴增。這一節(p.122)還沒有給答案,答案分兩段出現在後面:

第一段在 p.123:高溫下「the total number of excited phonons is proportional to T」,碰撞頻率正比於可碰撞的聲子數,所以 ℓ ∝ 1/T。單看這一條,0 °C(273 K)降到 −190 °C(83 K)只會讓 ℓ 漲 273/83 ≈ 3.3 倍。不夠。

第二段在 p.125–126 的 Umklapp Processes:真正造成熱阻的不是所有三聲子碰撞,而是需要 K₁ + K₂ = K₃ + G 的 umklapp 過程;而 U 過程要求參與的聲子波向量夠大(約為布里淵區尺度的一半以上),這種高能聲子的數目在低溫下按 exp(−θ_D/2T) 這種指數形式凍結,遠比 1/T 陡峭。指數凍結才是 ℓ 暴增的主因。Table 2 的表題註記「The ℓ's obtained in this way refer to umklapp processes」正是在預告這件事。

這個暴增不會無限持續。當 ℓ 長到與樣品尺寸相當時,聲子開始主要撞晶界而不是撞彼此,ℓ 被幾何截斷在樣品尺寸(Casimir 極限)。此後 ℓ 不再增加,K = ⅓Cvℓ 就完全由 C ∝ T³ 主導,K 隨溫度下降。所以純淨介電晶體的 K(T) 曲線是一個倒 U 形:低溫端 T³ 上升、高溫端 1/T 下降、中間(典型在 θ_D/20 附近)有一個尖銳的極大。Table 2 的 −190 °C 那一欄,大致就落在這個峰的高溫側。

注意這裡有一個尺度上的自我檢查。ℓ = 540 Å = 54 nm。如果樣品只有幾百奈米厚,那麼 −190 °C 時 Kn = ℓ/L 已經接近 0.1–1,前面講的彈道/擴散判準就要拿出來用了。這就是為什麼「低溫」與「小尺寸」這兩件事在熱輸運裡總是一起出問題。

(本區塊中「C 下降 × K 上升 = ℓ 必須暴增」的反推、指數凍結與 Casimir 極限的討論,屬於本講解的補充與跨節整合;Kittel 在 p.122 只給了表格數字,在 p.123 給了 ℓ ∝ 1/T,在 p.125 起才處理 umklapp。)

Thermal Resistivity of Phonon Gas(p.123):ℓ 由誰決定

子節一開場就把答案分成兩類:

「The phonon mean free path ℓ is determined principally by two processes, geometrical scattering and scattering by other phonons.」

(聲子平均自由程 ℓ 主要由兩種過程決定:幾何散射,以及被其他聲子散射。)

第一類:幾何散射。接下來這句同時解釋了幾何散射什麼時候會獨大,也解釋了聲子-聲子散射為什麼需要非諧性:

「If the forces between atoms were purely harmonic, there would be no mechanism for collisions between different phonons, and the mean free path would be limited solely by collisions of a phonon with the crystal boundary, and by lattice imperfections.」

(如果原子間的力是純諧振的,那麼不同聲子之間就沒有碰撞的機制,而平均自由程將只受限於聲子與晶體邊界的碰撞,以及晶格缺陷。)

「There are situations where these effects are dominant.」(有些情況下這些效應是主導的。)

這句話的分量要抓住。純諧振晶體裡,簡正模是哈密頓量的嚴格本徵態,各模之間沒有任何耦合:你今天放進去一顆 K = K₀ 的聲子,一萬年後它還是 K = K₀。沒有碰撞就沒有有限的自由程,ℓ = ∞,K = ⅓Cvℓ = ∞——純諧振晶體的熱導率是無限大的。這不是數學病態,這是一個真實而深刻的陳述:熱阻完全是非諧性的產物。第 5 章前半用諧振近似算熱容算得非常成功(Debye T³ 定律精確得不得了),但同一個近似用來算熱導率會得到無限大——這是諧振近似最戲劇性的失敗,也是為什麼這一章的後半非講非諧性不可。

Kittel 說「There are situations where these effects are dominant」,指的就是上一節提過的低溫、小樣品的情形:聲子-聲子散射凍結掉之後,剩下的只有晶界與缺陷,ℓ 被樣品尺寸釘住。

第二類:聲子-聲子散射,需要非諧耦合。

「With anharmonic lattice interactions, there is a coupling between different phonons which limits the value of the mean free path. The exact states of the anharmonic system are no longer like pure phonons.」

(有了非諧的晶格交互作用,不同聲子之間就有了耦合,限制了平均自由程的值。非諧系統的精確態不再像純粹的聲子。)

最後那句話很有意思,值得多想一秒:「精確態不再像純粹的聲子」。這是在提醒讀者,「聲子互相碰撞」這個說法本身就是一種近似語言。嚴格說,非諧晶體的本徵態根本不是聲子;我們之所以還能說「聲子碰撞」,是因為非諧項小,可以把諧振簡正模當成好的零階基底,再把非諧項當微擾處理——三次非諧項(∝ x³)給出三聲子過程,四次項給出四聲子過程。第 5 章 p.120 的式 (38) U(x) = cx² − gx³ − fx⁴ 就是這個展開,而同一頁明說「Three-phonon processes are caused by third-order terms in the lattice potential energy.」(三聲子過程由晶格位能中的三次項造成。)——p.120 用這個 −gx³ 項算出熱膨脹 (40),p.123 用同一個項造出聲子碰撞。熱膨脹與熱阻是同一個非諧項的兩個後果,這是本章前後兩半的真正連結。

(本區塊英文引句逐字取自 p.123(連結段落的引句取自 p.120);「純諧振晶體的 K 是無限大」這個推論、以及把它認定為諧振近似最戲劇性失敗的評價,是本講解補充。)

高溫下的聲子平均自由程

ℓ ∝ 1/T

p.123

ℓ:聲子平均自由程;T:絕對溫度。此關係在原文中以行內式印出,沒有給式號,所以本條不列 eqNumber。

出處推導
兩個來源。理論:原文逐字「The theory of the effect of anharmonic coupling on thermal resistivity predicts that ℓ is proportional to 1/T at high temperatures, in agreement with many experiments.」(非諧耦合對熱阻效應的理論預測,在高溫下 ℓ 正比於 1/T,這與許多實驗相符。)定性論證(原文逐字):「We can understand this dependence in terms of the number of phonons with which a given phonon can interact: at high temperature the total number of excited phonons is proportional to T. The collision frequency of a given phonon should be proportional to the number of phonons with which it can collide, whence ℓ ∝ 1/T.」論證鏈是三環:高溫下被激發的聲子總數 ∝ T → 給定聲子的碰撞頻率 ∝ 可與之碰撞的聲子數 ∝ T → 自由程 ∝ 1/碰撞頻率 ∝ 1/T。
極限行為
「高溫」的意思是 T ≫ θ_D(Debye 溫度)。此時每個模態的佔據數 ⟨n⟩ 都遠大於 1,且 ⟨n⟩ ≈ k_BT/ℏω 正比於 T。低溫(T ≪ θ_D)時這個關係完全失效:不但係數變了,連函數形式都變了——能參與 umklapp 過程的高能聲子數按 exp(−θ_D/aT)(a 為 2 附近的常數)指數凍結,ℓ 的增長遠快於 1/T。這正是 Table 2 中 ℓ 從 40 Å 跳到 540 Å(13.5 倍,而 273/83 只有 3.3 倍)所要求的。再更低溫,ℓ 撞上樣品尺寸而飽和(Casimir 極限),與 T 無關。推論到 K:高溫下 C 已飽和在 Dulong–Petit 值、v 不變,故 K ∝ 1/T——這是純淨介電晶體高溫熱導率最著名的實驗規律。
典型數值
用 Table 2 檢驗:0 °C = 273 K、−190 °C = 83 K,比值 3.29。若純由 ℓ ∝ 1/T 支配,石英的 ℓ 應從 40 Å 變成約 132 Å;實際是 540 Å。差了 4 倍——這個落差就是 umklapp 凍結(指數項)的貢獻,也直接證明 83 K 對石英而言已經不算「高溫」(石英的 Debye 溫度量級在數百 K)。

補完「高溫下激發聲子總數 ∝ T」這一步補推導

為什麼要補這一步。Kittel 在 p.123 只寫「at high temperature the total number of excited phonons is proportional to T」,把它當成已知。它確實是第 5 章前半的直接推論,但書上沒有把「Planck 分布 → 高溫展開 → 總數 ∝ T」這條線明確接起來,讀者容易停在「好像記得有這回事」的層次。上面七步把它接完整。

要注意的一個陷阱:τ 在這一章有兩個意思。p.107 定義 τ ≡ k_B T(一個能量,Kittel 稱之為 fundamental temperature,全書熱力學部分都這樣用);而 p.122 的 (43) 裡的 τ 是碰撞間的平均時間(「τ is the average time between collisions」)。同一章、相隔十五頁、同一個符號、兩個完全不同的物理量。這在讀 Kittel 時一定要靠上下文分辨,尤其是在做這種跨節推導、把 Planck 分布代進動力論公式的時候。本文在需要區分時一律把前者寫成 k_BT。

這條論證鏈的極限在哪裡。它有一個被藏起來的假設:「碰撞頻率正比於可碰撞的聲子數」。這句話對 N 過程與 U 過程一視同仁——但 p.123–125 的整個重點正是兩者在熱阻上的效果天差地別。所以嚴格說,ℓ ∝ 1/T 描述的是「聲子-聲子散射的總頻率」,而能產生熱阻的只是其中的 U 過程部分。高溫下 U 過程沒有被凍結(能量夠高的聲子到處都是),兩者的溫度依賴恰好一致,所以這個粗糙論證在高溫下剛好答對。低溫下它就錯了,錯的方式正是上一條 limits 說的:U 過程被指數凍結,總散射率與熱阻散射率脫鉤。

  1. 出發點:第 5 章 p.107 的 Planck 分布(式 (2),p.108 重述為 (7))給出單一模態的平衡佔據數 ⟨n⟩ = 1 / (exp(ℏω/τ) − 1),其中 p.107 定義 τ ≡ k_B T。
  2. 高溫展開:T ≫ θ_D 時,對所有被激發的模態都有 ℏω/τ = ℏω/(k_BT) ≪ 1。用 exp(x) ≈ 1 + x + x²/2 展開分母:exp(ℏω/k_BT) − 1 ≈ ℏω/(k_BT) + ½(ℏω/k_BT)² + …
  3. 取領頭項:⟨n⟩ ≈ k_BT / (ℏω)。這就是古典(Rayleigh–Jeans)極限——每個模態的能量 ⟨n⟩ℏω ≈ k_BT,正是能量均分定理。
  4. 對所有模態求和:聲子總數 N_ph = Σ_{K,p} ⟨n_{K,p}⟩ ≈ Σ_{K,p} k_BT/(ℏω_{K,p}) = (k_BT/ℏ) · Σ_{K,p} 1/ω_{K,p}。
  5. 關鍵觀察:Σ_{K,p} 1/ω_{K,p} 只依賴晶體的色散關係與模態總數(3N 個),完全不含溫度。所以 N_ph ∝ T。這就是原文「at high temperature the total number of excited phonons is proportional to T」那句話的完整內容。
  6. 接到碰撞頻率:三聲子過程的散射率在費米黃金律下正比於末態/中間態的佔據因子;高溫下這些佔據數都 ∝ T,所以碰撞頻率 1/τ ∝ N_ph ∝ T。
  7. 接到自由程:ℓ = vτ,v 與 T 無關,故 ℓ ∝ τ ∝ 1/T。∎

關鍵轉折:限制自由程還不夠,必須能建立局域熱平衡

p.123 的後半是整段最精采、也最容易讀過去的部分。Kittel 先拋出一個看似多餘的要求:

「To define a thermal conductivity there must exist mechanisms in the crystal whereby the distribution of phonons may be brought locally into thermal equilibrium. Without such mechanisms we may not speak of the phonons at one end of the crystal as being in thermal equilibrium at a temperature T₂ and those at the other end in equilibrium at T₁.」

(要定義熱導率,晶體中必須存在某種機制,使聲子的分布能夠被帶到局域的熱平衡。沒有這樣的機制,我們就不能說晶體一端的聲子處在溫度 T₂ 的熱平衡、另一端處在 T₁ 的熱平衡。)

這句話直接回頭咬住 (41)。(41) 裡寫著 dT/dx——要寫得出「x 處的溫度 T(x)」,就必須承認 x 處有一個良好定義的溫度;而溫度只對熱平衡態有定義。所以 Fourier 定律不只需要「碰撞夠多」(擴散條件),還需要「碰撞能讓分布鬆弛到平衡形式」。這是兩個獨立的要求。

接著 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.」

(光有限制平均自由程的辦法是不夠的,還必須有建立聲子局域熱平衡分布的辦法。)

然後給出第一個反例:靜態缺陷與晶界。

「Phonon collisions with a static imperfection or a crystal boundary will not by themselves establish thermal equilibrium, because such collisions do not change the energy of individual phonons: the frequency ω₂ of the scattered phonon is equal to the frequency ω₁ of the incident phonon.」

(聲子與靜態缺陷或晶體邊界的碰撞,本身並不能建立熱平衡,因為這類碰撞不改變個別聲子的能量:散射後聲子的頻率 ω₂ 等於入射聲子的頻率 ω₁。)

為什麼「不改變頻率」就不夠?因為熱平衡分布 ⟨n(ω)⟩ = 1/(exp(ℏω/k_BT) − 1) 規定的是能量在不同頻率之間怎麼分配。靜態散射體只改變聲子的方向(K → K',|K| 不變),完全不在頻率之間搬運能量。所以如果你丟進去一個頻率分布很怪的聲子群(例如全部集中在某個 ω₀),純彈性散射把它們的方向打亂之後,它們的頻率分布還是那麼怪——你可以說它們「動量平衡」了,但不能說它們「熱平衡」了,因為根本沒有一個溫度能描述這個分布。

這是一個容易被忽略的區分:動量鬆弛 ≠ 能量再分配 ≠ 熱化。靜態散射做到了第一件(部分),沒做到第二、三件。接下來 Kittel 要給的第二個反例,則恰恰相反——三聲子 N 過程做到了能量再分配,卻沒做到動量鬆弛。

(三段英文引句逐字取自 p.123 並經算繪目視核對——特別是 T₂ 與 T₁ 的先後順序,PDF 文字層會把下標浮出並打亂;「為什麼不改變頻率就不夠」的說明,以及「動量鬆弛 ≠ 能量再分配 ≠ 熱化」這個三分法,是本講解補充。)

三聲子正常過程(N process)

K₁ + K₂ = K₃

(45)  p.123

K₁、K₂:碰撞前兩顆聲子的波向量;K₃:碰撞後聲子的波向量。這是三聲子過程中波向量完全相加的那一類,後文(p.125)稱為 normal process(N 過程)。原文引言:「It is rather remarkable also that a three-phonon collision process (45) will not establish equilibrium, but for a subtle reason: the total momentum of the phonon gas is not changed by such a collision.」

出處推導
非諧晶格交互作用的三次項(U(x) 展開中的 x³ 項)在微擾論下產生的三聲子躍遷,配合晶體中的波向量選擇律。晶體中嚴格成立的選擇律是「總波向量守恆,至多相差一個倒晶格向量 G」(第 4 章 p.100 的 (30)(31)(32) 已用同一原理處理 X 射線光子的彈性散射,以及伴隨聲子產生/吸收的非彈性散射)。(45) 是 G = 0 的特例。同一頁之後 p.125 的 (47) K₁ + K₂ = K₃ + G 才是 G ≠ 0 的 umklapp 情形。
極限行為
能量守恆也必須同時成立:p.125 對 Fig. 17 的說明中明說「For all processes, N or U, energy must be conserved, so that ω₁ + ω₂ = ω₃.」所以 (45) 只是兩條守恆律中的一條。這類碰撞做到什麼:它在不同頻率之間搬運能量(ω₁ + ω₂ → ω₃),因此能把聲子分布熱化成 Planck 形式。它做不到什麼:它不改變聲子氣的總動量 J(見 (46)),因此不能讓一個帶淨漂移的分布停下來。結論就是原文所說的:只有 N 過程時,熱阻為零、熱導率無限大。極端情形:K₂ = 0 沒有物理意義(沒有第二顆聲子);K₁ = K₂ 時是兩顆同樣的聲子併成一顆 K₃ = 2K₁,只要 2K₁ 還落在第一布里淵區內就是 N 過程,一旦超出區邊界就必須折回,那就變成 U 過程了——N 與 U 的分界完全由布里淵區的幾何決定。

聲子氣的總動量(晶體動量)

J = Σ_K n_K ℏK

(46)  p.123

J:聲子氣的總(晶體)動量;n_K:具有波向量 K 的聲子數(原文:「Here n_K is the number of phonons having wavevector K.」);ℏK:單顆聲子攜帶的晶體動量;求和跑遍所有 K。

出處推導
把單顆聲子的晶體動量 ℏK(第 4 章 p.100 的 PHONON MOMENTUM 一節)按佔據數加總。原文的用法是:對 (45) 型碰撞,「the phonon momentum (46) is conserved, because on collision the change in J is K₃ − K₂ − K₁ = 0」——碰撞中 J 的改變量正是 K₃ − K₂ − K₁,由 (45) 知它等於零。
極限行為
J = 0 是真正的熱平衡:晶體色散恆滿足 ω(K) = ω(−K)(時間反演/倒易對稱的結果,各向異性晶體也成立,不必假設 ω 只依賴 |K|),所以平衡時 ⟨n_K⟩ = ⟨n_−K⟩,Σ_K n_K ℏK 中每一對 ±K 互相抵消,求和後 J = 0,沒有淨熱流。J ≠ 0 是「漂移平衡」:分布仍是平衡形式,但整體被推向某個方向,聲子氣像一陣風一樣沿棒子吹。原文:「An equilibrium distribution of phonons at a temperature T can move down the crystal with a drift velocity which is not disturbed by three-phonon collisions of the form (45).」此時有淨能量輸運卻不需要溫度梯度,代入 (41) 相當於 dT/dx = 0 而 j_U ≠ 0,即 K → ∞。J 不是真動量:它是晶體動量,只在模 ℏG 的意義下守恆——這正是為什麼 (47) 的 umklapp 過程能改變 J 而不違反任何守恆律。第 4 章 p.100 已經說過真動量是嚴格守恆的,差額由整塊晶體帶走。U 過程下 J 的改變量是 −ℏG,是一個布里淵區尺度的巨量,一次碰撞就足以把聲子流的方向翻轉(p.125 對 Fig. 17b 的說明:「in the umklapp process the direction of the x-component of the phonon flux has been reversed」)。

為什麼「J 守恆」就等於「熱阻為零」:把這一步的邏輯寫死補推導

推導前先把三個「平衡」分清楚,因為 p.123 的整段論證全靠這個區分:

名稱分布長什麼樣誰能建立它有沒有熱流
完全熱平衡Planck 分布,J = 0需要既改變能量、又改變 J 的過程(U 過程)沒有
漂移平衡Planck 分布被整體推移,J ≠ 0N 過程可以維持它有,而且不需要梯度
非熱分布頻率分布不是 Planck 形式只有彈性散射時會停在這裡溫度根本沒定義

Kittel 用的詞是「complete thermal equilibrium」——那個 complete 就是在區分第一列與第二列。N 過程能把你從第三列帶到第二列(它重新分配能量),但永遠帶不到第一列(它動不了 J)。

補充一個書上沒有寫、但把這件事講到底需要的形式。(45) 型碰撞所守恆的量有兩個:總能量 Σ n_K ℏω_K 與總晶體動量 J = Σ n_K ℏK。統計力學的一般結論是:碰撞算符的零空間由守恆量張成,所以 N 過程碰撞算符的穩態解不是只有一個 Planck 分布,而是一整族帶漂移速度 u 的分布

n_K = 1 / { exp[ ℏ(ω_K − u·K) / k_B T ] − 1 }

這叫做位移的(drifted)Bose–Einstein 分布,u 就是原文所說的 drift velocity。u 是自由參數,代表 J 這個守恆量可以取任意值而碰撞不加干涉。這正是「thermal resistivity is zero」在數學上的準確說法。(此形式為本講解補充;Kittel 在 p.123 用的是文字敘述「can move down the crystal with a drift velocity」,並未寫出這個分布函數。)

這件事在物理上的類比在哪裡:它與超流、與無雜質金屬中的電子氣(電子-電子碰撞守恆總動量,因此純電子-電子散射也不產生電阻)是同一個結構。凡是「碰撞守恆了那個攜帶流的守恆量」,該項流的阻抗就是零。這是輸運理論裡最反覆出現的定理之一,值得記牢。

  1. 前提一(原文 (46) 的後續):對所有形如 (45) 的三聲子碰撞,J 的改變量是 K₃ − K₂ − K₁ = 0,故 J 守恆。
  2. 前提二(原文 p.123):「For a distribution with J ≠ 0, collisions such as (45) are incapable of establishing complete thermal equilibrium because they leave J unchanged.」(對於 J ≠ 0 的分布,(45) 這類碰撞無法建立完全的熱平衡,因為它們讓 J 保持不變。)
  3. 推論(原文 p.123 末至 p.124,此句跨頁):「If we start a distribution of hot phonons down a rod with J ≠ 0, the distribution will propagate down the rod with J unchanged. Therefore there is no thermal resistance.」(若我們讓一個 J ≠ 0 的熱聲子分布沿棒子出發,這個分布會帶著不變的 J 沿棒子傳播下去。因此沒有熱阻。)
  4. 把它翻成穩態輸運的語言:熱阻的微觀定義是「維持一個穩定的熱流需要多大的溫度梯度」。如果一個攜帶熱流的分布(J ≠ 0)可以在沒有任何驅動的情況下自己維持下去,那麼維持這個熱流所需的梯度是零。
  5. 代回 (41):j_U ≠ 0 而 dT/dx = 0 ⟹ K = −j_U / (dT/dx) → ∞。熱導率發散,熱阻率 1/K = 0。
  6. 另一個看法:J 是一個守恆量,而守恆量沒有弛豫時間。要有有限的熱阻,就必須有機制讓「攜帶熱流的那個量」衰減。在聲子氣裡攜帶熱流的正是 J(熱流 ∝ Σ n_K ℏω_K v_K,在 Debye 近似 ω = vK 下與 J 同方向且成正比)。N 過程守恆 J,所以它對熱流的弛豫完全沒有貢獻——它會把分布內部攪均勻,但攪不掉整體的漂移。
  7. 因此真正需要的是「不守恆 J 的碰撞」。靜態缺陷與晶界確實不守恆 J(它們把動量交給了晶體整體),但它們不改變聲子頻率,所以不能熱化;U 過程既改變 J(改變 −ℏG)又重新分配能量,兩件事都做。這就是為什麼 Table 2 的表題要說那些 ℓ「refer to umklapp processes」。

Fig. 16a(p.124):無摩擦管裡的漂移平衡氣體

網頁上沒有圖,所以先把圖本身描述出來。

圖在畫什麼。一根很長的水平管子,上下管壁畫成斜線剖面(表示是實體管壁,但設定為無摩擦)。管子兩端都是開口的:左端標「Hot gas in」(熱氣進入),右端標「Hot gas out」(熱氣流出)。管內畫滿密密麻麻的小箭頭,代表氣體分子的速度向量——它們指向各種方向(隨機熱運動),但整體明顯偏向右。管子上方另有一個大箭頭指向右,標示「Mass flow of gas」(氣體的質量流)。

圖說逐字。「Flow of gas molecules in a state of drifting equilibrium down a long open tube with frictionless walls. Elastic collision processes among the gas molecules do not change the momentum or energy flux of the gas because in each collision the velocity of the center of mass of the colliding particles and their energy remain unchanged. Thus energy is transported from left to right without being driven by a temperature gradient. Therefore the thermal resistivity is zero and the thermal conductivity is infinite.」

(氣體分子在無摩擦管壁的長開口管中處於漂移平衡狀態下的流動。氣體分子之間的彈性碰撞過程不改變氣體的動量或能量通量,因為在每次碰撞中,碰撞粒子的質心速度與它們的能量都保持不變。因此能量從左傳到右,而不需要由溫度梯度驅動。所以熱阻率為零,熱導率為無限大。)

這張圖是整段論證的「基準案例」。要抓住的是三個設計選擇,每一個都在對應一個物理要件:

管壁無摩擦 ⟺ 沒有外部動量匯。真實管子的壁會把分子的縱向動量吸走(這就是黏滯阻力),氣流會慢下來。設成無摩擦,就是明確地把「動量可以流失到系統外」這條路封死。這對應到晶體裡的「沒有晶界散射、沒有缺陷散射」。

兩端開口、有質量流 ⟺ J ≠ 0 的初始條件。左端持續灌入熱氣,右端讓它流走,所以管內維持一個穩定的淨漂移。這對應到 Fig. 16c 裡的「聲子源與聲子槽」。

只有彈性的兩體碰撞 ⟺ 只有守恆總動量的內部碰撞。圖說裡「the velocity of the center of mass of the colliding particles and their energy remain unchanged」把這件事講得非常精確:兩體彈性碰撞守恆質心速度(等價於守恆總動量)與總能量。這正是 (45) 型 N 過程在分子氣體中的對應物。

結論的推理路徑。內部碰撞守恆總動量 → 漂移速度不會衰減 → 熱氣可以一路輸運能量到右端 → 產生能量通量卻不需要任何溫度梯度 → 由 (41),K = −j_U/(dT/dx) 在 dT/dx → 0 而 j_U 有限時發散 → 熱阻率為零。

這張圖的功能,是讓讀者在一個完全熟悉的古典場景裡先接受「碰撞很多但熱阻為零」這件反直覺的事,然後才有辦法接受聲子氣裡同樣的事。

(圖形描述依 p.124 原圖算繪後目視判讀,圖說逐字取自 p.124;「三個設計選擇」的對應關係與結論的推理路徑是本講解補充。)

Fig. 16b(p.124):封閉管——熱導率的「正規」定義場景

圖在畫什麼。同樣一根水平管子,但這次兩端都被封死(左右兩端畫上與上下管壁相同的斜線剖面,形成一個密閉盒子)。左端外標「Hot end」(熱端),右端外標「Cold end」(冷端)。管內同樣畫滿小箭頭,但這次箭頭的方向是真正隨機的,沒有整體偏向——左右、上下都有,看不出淨流。與 Fig. 16a 對比,上方那個「Mass flow of gas」的大箭頭不見了。

圖說逐字。「The usual definition of thermal conductivity in a gas refers to a situation where no mass flow is permitted. Here the tube is closed at both ends, preventing the escape or entrance of molecules. With a temperature gradient the colliding pairs with above-average center of mass velocities will tend to be directed to the right, those with below-average velocities will tend to be directed to the left. A slight concentration gradient, high on the right, will be set up to enable the net mass transport to be zero while allowing a net energy transport from the hot to the cold end.」

(氣體中熱導率的通常定義,指的是不允許有質量流的情況。這裡管子兩端封閉,防止分子逸出或進入。在溫度梯度下,質心速度高於平均的碰撞對傾向於指向右方,低於平均的傾向於指向左方。管中會建立起一個微小的濃度梯度、右側較高,使得淨質量輸運為零,同時仍允許從熱端到冷端有淨能量輸運。)

這張圖在做的事,是把 (42) 的適用場景界定清楚。Fig. 16a 說明「有質量流時 K 發散」,那麼我們平常說的「氣體的熱導率」到底是什麼?答案:是禁止質量流這個額外約束下的量。這不是技術細節,這是定義的一部分——「The usual definition ... refers to a situation where no mass flow is permitted」。

最精微的一句話在中間。「the colliding pairs with above-average center of mass velocities will tend to be directed to the right, those with below-average velocities will tend to be directed to the left」——注意他講的是碰撞對的質心速度,不是個別分子。為什麼?因為熱端的分子平均動能大、速度大,它們組成的碰撞對質心速度也大;這些「快的對」在碰撞後仍保有原來的質心速度(彈性碰撞守恆質心速度),於是繼續往右跑。反過來,從冷端來的「慢的對」往左跑。結果是:快的往右、慢的往左,淨能量向右(熱端 → 冷端),但這樣一來質量會淨往右搬。

所以需要那個「微小的濃度梯度」。為了讓淨質量輸運歸零,氣體會自己在右側(冷端)堆高一點密度,用「多一點慢分子往左」去抵消「少一點快分子往右」。這是一個自洽的穩態:密度梯度不是外加的,是系統為了滿足「無質量流」這個邊界條件自己長出來的。這在熱輸運裡有個名字叫熱擴散/Soret 效應的靜態對應。

對應到晶體。這裡有一個重要的不對稱:氣體分子守恆(所以需要濃度梯度去平衡),但聲子數不守恆(聲子可以在晶體內任意產生與湮滅)。所以聲子氣沒有「不准質量流」這個約束,也不需要長出濃度梯度。聲子氣裡對應的約束不是「質量流為零」,而是「J 必須能鬆弛掉」——而那需要 U 過程。Kittel 把 Fig. 16b 放在這裡,正是為了讓讀者看清楚:氣體用「封住兩端」解決的問題,晶體必須用另一種辦法解決。

(圖形描述依 p.124 原圖算繪後目視判讀,圖說逐字取自 p.124;「為什麼需要那個濃度梯度」的展開、以及聲子數不守恆的對照,是本講解補充。)

Fig. 16c(p.124):聲子源與聲子槽——把類比搬進晶體

圖在畫什麼。一個水平長方框代表晶體。左側標「Phonon source」(聲子源),右側標「Phonon sink」(聲子槽)。框內畫滿小箭頭代表聲子,方向有些散亂但整體明顯偏右;框的中央標著「N processes」(N 過程)。框的上方有一個大箭頭指向右,標「Net flow of phonons」(聲子的淨流動)。

整張圖的視覺結構與 Fig. 16a 刻意做成一樣——大箭頭、偏向右的小箭頭、兩端的進出口。這是圖形上的明示:這兩張圖說的是同一件事。

圖說逐字。「In a crystal we may arrange to create phonons chiefly at one end, as by illuminating the left end with a lamp. From that end there will be a net flux of phonons toward the right end of the crystal. If only N processes (K₁ + K₂ = K₃) occur, the phonon flux is unchanged in momentum on collision and some phonon flux will persist down the length of the crystal. On arrival of phonons at the right end we can arrange in principle to convert most of their energy to radiation, thereby creating a sink for the phonons. Just as in (a) the thermal resistivity is zero.」

(在晶體中我們可以安排主要在一端產生聲子,例如用一盞燈照射左端。從該端會有一道朝向晶體右端的淨聲子通量。如果只發生 N 過程 (K₁ + K₂ = K₃),聲子通量在碰撞中動量保持不變,因而會有一部分聲子通量一路延續到晶體的另一端。當聲子到達右端時,我們原則上可以安排把它們的大部分能量轉換成輻射,從而為聲子創造一個匯。正如 (a) 中一樣,熱阻率為零。)

這張圖把 Fig. 16a 的每一個元件都換成了晶體裡的對應物:

Fig. 16a(氣體)Fig. 16c(晶體)物理角色
左端「Hot gas in」用燈照射左端產生聲子源:提供 J ≠ 0 的初始條件
分子間彈性碰撞N 過程 K₁ + K₂ = K₃守恆總動量的內部碰撞
無摩擦管壁(隱含)沒有 U 過程、沒有晶界與缺陷散射沒有外部動量匯
右端「Hot gas out」把聲子能量轉成輻射匯:讓穩態能成立
熱阻率為零熱阻率為零同一結論

「用燈照射左端」這個設定要特別留意。它不是隨口一說。它的作用是提供一個不靠溫度梯度的聲子源——燈光在左端被吸收、產生聲子,這個過程與棒子內部的溫度分布無關。有了它,我們就能製造出一個「有淨聲子流但沒有溫度梯度」的狀態,從而讓 K → ∞ 這個結論不只是數學推論,而是一個原則上可以設計出來的實驗。

「原則上」這三個字也要留意。「we can arrange in principle to convert most of their energy to radiation」——Kittel 很誠實:這是一個思想實驗。真實晶體裡 U 過程永遠存在(只是在低溫下被凍結),所以純 N 過程的情形不可能完全實現。但這個思想實驗的價值不在於能不能做,而在於它把「什麼機制才產生熱阻」這個問題的答案隔離出來了。

(橋接)Fig. 16 還有第 (d) 張,但它印在下一頁 p.125,屬於 Umklapp Processes 小節的開頭,本段範圍只到 p.124。它的圖說是:「In U processes there is a large net change in phonon momentum in each collision event. An initial net phonon flux will rapidly decay as we move to the right. The ends may act as sources and sinks. Net energy transport under a temperature gradient occurs as in (b).」——U 過程中每一次碰撞事件都有很大的淨聲子動量改變,初始的淨聲子通量會隨著往右移動而迅速衰減。這就是 (a)(b)(c) 三張圖鋪陳了整整一頁之後要交出的答案。

(圖形描述依 p.124 原圖算繪後目視判讀;Fig. 16c 圖說逐字取自 p.124,Fig. 16d 圖說逐字取自 p.125。對應表與其後的討論是本講解補充。)

把 Fig. 16a–c 統一起來看:熱阻的存在條件清單我的判讀

三張圖看完,如果只記住一句話,應該是這一句:熱阻不是「有碰撞」造成的,是「碰撞破壞了某個守恆量」造成的。

把 p.123–124 的所有討論整理成一張檢查表,一個熱傳機制要能產生有限的、可定義的熱導率,必須同時滿足:

要件為什麼需要靜態缺陷/晶界N 過程U 過程
限制平均自由程 ℓ否則 K = ⅓Cvℓ 發散;也否則系統是彈道的,(41) 根本不適用✔✔✔
在頻率之間重新分配能量(建立 Planck 分布,讓「溫度」有定義)(41) 裡的 T(x) 需要局域熱平衡✘(ω₂ = ω₁)✔✔
鬆弛總晶體動量 J否則漂移分布不衰減,熱流不需要梯度就能維持✔(動量給了晶體整體)✘(J 守恆)✔(J 改變 −ℏG)

三欄看下來,只有 U 過程三項全過。這就是為什麼 Table 2 的表題會寫「The ℓ's obtained in this way refer to umklapp processes」——那句在 p.122 看起來突兀的註記,到這裡才完全說得通。

兩個反例各自缺了不同的那一塊,這個對稱性很漂亮,值得記住:

  • 靜態缺陷/晶界:能鬆弛 J,不能熱化(不改變 ω)。
  • N 過程:能熱化(重分配 ω),不能鬆弛 J。
  • 兩者單獨都不足以定義熱導率。

一個容易混淆的追問:那把兩者加在一起夠不夠?在有限樣品裡,晶界散射確實提供了 J 的鬆弛通道,N 過程提供了熱化通道,兩者並存時確實會有有限的 K。但這時的 K 依賴樣品尺寸(因為 J 的鬆弛尺度是樣品尺寸),所以它不是材料的本徵熱導率——這正是低溫 Casimir 區的情形,也再度回到 p.121 那段「彈道時通量只看 ΔT、與長度無關」的警告。要得到一個真正只屬於材料的 K,必須讓 U 過程主導,也就是要求樣品尺寸 ≫ U 過程的自由程。

這張檢查表在研究上的用法。當你面對一個「熱導率量出來怪怪的」的系統(合金、超晶格、奈米線、多晶薄膜、非晶質),問這三個問題通常就能定位問題:ℓ 被誰限制?能量能不能在模態之間重分配?動量的鬆弛尺度是材料內稟的還是幾何給的?第三個問題的答案若是「幾何」,那麼你手上那個數字就不是材料常數。

(本區塊是對 p.121–124 內容的整合與表格化,屬本講解補充;三項要件本身全部來自原文的敘述。)

五個常見誤解我的判讀

誤解一:「K = ⅓Cvℓ 裡的 C 是莫耳熱容。」

不是。原文寫得很清楚:「C is the heat capacity per unit volume」,單位是 J cm⁻³ K⁻¹。Table 2 的欄位標頭也印著 J cm⁻³K⁻¹。用莫耳熱容代進去,量綱馬上就錯(會多出一個 mol⁻¹)。每次代數字之前先串一遍單位,這個錯誤就不會發生。

誤解二:「⅓ 是從三個維度來的,½ 是從雙向來的,所以應該是 ⅙。」

不是。推導裡確實同時出現 ½(單向粒子通量 ½n⟨|v_x|⟩)和一個因子 2(自由程兩端相距 2ℓ_x),這兩個相乘等於 1,互相抵消。唯一存活下來的數值因子是 ⟨v_x²⟩ = ⅓⟨v²⟩ 的那個 ⅓。這一點書上沒寫,本文的推導區塊逐步補出來了。

誤解三:「τ 是 k_BT。」

在這一節不是。第 5 章 p.107 定義 τ ≡ k_BT(一個能量),但 p.122 的 (43) 裡 τ 是「the average time between collisions」(一個時間)。同一章、同一個符號、兩個意思。做跨節推導(例如把 Planck 分布代進 (43))時,這個衝突會直接害你算錯量綱。本文一律把溫度寫成 k_BT 以避開它。

誤解四:「碰撞越多,熱阻越大。」

這是 p.123–124 整段要打掉的直覺。N 過程是不折不扣的碰撞,碰撞率再高也不產生任何熱阻,因為它守恆 J。Fig. 16a 的無摩擦管更是把這件事畫成了圖:管裡的分子撞來撞去,熱阻率照樣是零。產生熱阻的是「破壞 J」的碰撞,不是「碰撞」本身。

誤解五:「Table 2 的 ℓ 是實驗量到的聲子平均自由程。」

不是。表題明說「Calculated from (44)」——ℓ 是拿實驗量到的 C 與 K、配上一個假設的 v = 5×10⁵ cm/sec 反推出來的。它是一個模型相依的推導量,而且是把所有頻率的聲子強行平均成單一數值的結果。真實的聲子自由程譜 ℓ(ω) 跨越好幾個數量級,用單一數字代表它會嚴重低估長波聲子的貢獻。本文的驗算還發現,四筆裡有一筆(NaCl −190 °C)連用書上自己的公式與數字都對不回來。

與 Cth/熱萃取研究的接點:這一段到底給了什麼我的判讀

如果你的研究主題是元件熱容/熱阻的萃取,p.121–124 這四頁給了三樣直接可用的東西,值得逐一點名。

第一,一個「你的模型合不合法」的判準。整套 Cth/Rth 緊緻模型(RC 熱網路、有限元熱模擬、3ω 法、TDTR 的擬合)全部建立在 Fourier 定律 (41) 之上。p.121 明確告訴你,(41) 的合法性等價於「熱是擴散的、不是彈道的」,而檢驗方法是看通量依賴梯度還是依賴兩端溫差。轉成實驗設計就是:做尺寸序列。熱阻對長度(或膜厚)的圖若不過原點、若在小尺寸端彎曲,你的 Fourier 假設已經在漏了。這是第一個該做、卻最常被跳過的檢查。

第二,一組數量級的直覺。Table 2 給了「0 °C(273 K,接近室溫)、一般介電晶體,ℓ 是幾十埃」這個錨點。這個數字讓你知道 Kittel 那個年代的塊材量測為什麼從來不用擔心彈道效應(樣品是毫米級,Kn ~ 10⁻⁵)。它同時也讓你知道為什麼今天要擔心:現代元件的特徵尺寸已經逼近甚至低於聲子自由程,而且矽的聲子自由程分布遠比 Table 2 的「幾十埃」寬——這正是 Table 2 那種「單一 ℓ」描述的失效之處。更重要的是 Table 2 示範了降溫會把 ℓ 推高一個數量級,所以低溫量測特別容易踩進過渡區。

第三,也是最有價值的一樣:熱阻的來源清單。p.123 的三段論證告訴你,熱阻需要「限制自由程 + 建立局域熱平衡 + 鬆弛晶體動量」三件事同時成立。這張清單直接可以拿來拆解一個奈米結構的熱阻是從哪來的:

  • 界面/晶界貢獻的是動量鬆弛但不熱化的部分——它給你尺寸相依的 Casimir 型熱阻,不是材料常數。
  • 體內 U 過程貢獻的是真正的材料內稟熱阻。
  • 而如果你的結構小到 U 過程來不及發生(Kn ≫ 1),那麼你量到的「熱阻」幾乎全是界面的,改材料不會有用,改界面才會。

這個拆解對元件設計的意義很直接:知道熱阻卡在哪一項,才知道該改什麼。

還要補一句 Kittel 這幾頁沒講、但做元件必須知道的限制。這裡談的全是介電固體中的聲子熱傳(p.122 明說「dielectric solids」)。真實元件裡還有電子熱傳(金屬互連、重摻雜區,第 6 章 p.156 的 Wiedemann–Franz 定律)、界面熱阻(Kapitza 阻抗)、以及電子-聲子之間的能量交換不夠快所造成的非平衡(熱載子效應)。p.121–124 是整個熱輸運圖像的骨架,不是全部。但骨架抓對了,後面所有修正才有地方掛。

(本區塊是跨領域連結與研究方法上的補充,不是原文內容;所引用的 Table 2 數字與 p.121、p.123 的論證則來自原文。)

本段回顧:四頁書的邏輯骨架我的判讀

把 p.121–124 壓縮成一條線:

1. 定義(p.121)。j_U = −K dT/dx,(41)。四個限定詞:穩態、長棒、有梯度、通量是每單位面積每單位時間的能量。

2. 定義蘊含的物理(p.121)。寫成「正比於梯度」而非「正比於兩端溫差」,就已經宣稱了熱傳是隨機(擴散)過程。彈道傳輸下通量與長度無關,(41) 失效。這是全書 ballistic vs diffusive 的原始出處。

3. 微觀公式(p.122)。氣體動力論給 K = ⅓Cvℓ,(42);推導路徑 (43) → (44)。½(單向流)與 2(自由程兩端)相乘抵消,⅓ 唯一地來自 ⟨v_x²⟩ = ⅓⟨v²⟩。Debye 首先把它搬到介電固體,C 取聲子熱容、v 取聲速、ℓ 取聲子平均自由程。

4. 數字(Table 2, p.122)。石英與 NaCl 在 0 °C 與 −190 °C 的 C、K,以及由 (44) 反推的 ℓ。0 °C 時幾十埃,−190 °C 時增加數倍到十餘倍。表題註明這些 ℓ 指的是 umklapp 過程。(自行驗算發現 NaCl −190 °C 那一筆用書上的公式與數字算出 162 Å 而非印出的 100 Å。)

5. ℓ 由誰決定(p.123)。兩類:幾何散射(晶界、缺陷)與聲子-聲子散射。純諧振晶體裡聲子不互撞,ℓ 只被幾何限制——熱阻完全是非諧性的產物。非諧耦合的理論預測高溫下 ℓ ∝ 1/T,論證鏈是:激發聲子總數 ∝ T → 碰撞頻率 ∝ 可碰撞聲子數 → ℓ ∝ 1/T。

6. 但限制 ℓ 還不夠(p.123)。要定義熱導率,必須有機制建立局域熱平衡,否則「這一端溫度 T₂、那一端 T₁」這句話沒有意義。靜態缺陷與晶界不改變聲子頻率(ω₂ = ω₁),不能熱化。

7. 更微妙的反例(p.123–124)。三聲子 N 過程 K₁ + K₂ = K₃ (45) 也不行,因為聲子動量 J = Σ n_K ℏK (46) 在這類碰撞中的改變量是 K₃ − K₂ − K₁ = 0。J ≠ 0 的漂移平衡分布可以不受阻礙地一路傳下去,熱阻為零、熱導率無限大。

8. 三個類比(Fig. 16a–c, p.124)。無摩擦開口管中的漂移氣體(熱阻為零)/兩端封閉的管(熱導率的正規定義場景,用微小濃度梯度維持零質量流)/晶體中以燈為聲子源、以輻射為聲子槽(只有 N 過程時熱阻同樣為零)。

9. 懸而未決的問題,答案在下一頁:什麼機制既改變聲子頻率、又改變 J?Umklapp 過程 K₁ + K₂ = K₃ + G(p.125 的 (47))。Fig. 16d 與 Fig. 17 就是為它畫的。