第 3 章的地圖:一章其實裝了兩門課(書 p.47–48 章首目錄頁)
第 3 章的標題是 Crystal Binding and Elastic Constants。章首目錄頁(書 p.47)照原文印的順序與頁碼是:CRYSTALS OF INERT GASES(p.49,下轄 Van der Waals-London Interaction p.53、Repulsive Interaction p.56、Equilibrium Lattice Constants p.58、Cohesive Energy p.59);IONIC CRYSTALS(p.60,下轄 Electrostatic or Madelung Energy p.60、Evaluation of the Madelung Constant p.64);COVALENT CRYSTALS(p.67);METALS(p.69);HYDROGEN BONDS(p.70);ATOMIC RADII(p.70,下轄 Ionic Crystal Radii p.72);ANALYSIS OF ELASTIC STRAINS(p.73,下轄 Dilation p.75、Stress Components p.75);ELASTIC COMPLIANCE AND STIFFNESS CONSTANTS(p.77,下轄 Elastic Energy Density p.77、Elastic Stiffness Constants of Cubic Crystals p.78、Bulk Modulus and Compressibility p.80);ELASTIC WAVES IN CUBIC CRYSTALS(p.80,下轄 Waves in the [100] Direction p.81、Waves in the [110] Direction p.82);SUMMARY(p.85)。
習題共 14 題,橫跨 p.85–88:1. Quantum solid(85)、2. Cohesive energy of bcc and fcc neon(86)、3. Solid molecular hydrogen(86)、4. Possibility of ionic crystals R⁺R⁻(86)、5. Linear ionic crystal(86)、6. Cubic ZnS structure(86)、7. Divalent ionic crystals(86)、8. Young's modulus and Poisson's ratio(87)、9. Longitudinal wave velocity(87)、10. Transverse wave velocity(87)、11. Effective shear constant(87)、12. Determinantal approach(87)、13. General propagation direction(88)、14. Stability criteria(88)。前七題屬鍵結,後七題屬彈性——習題的切分線和目錄的切分線完全一致。
看這張目錄要抓住一件事:這一章的前半(p.49–72)與後半(p.73–85)在教科書的傳統分類裡其實是兩門課。前半是化學鍵與內聚能,後半是連續介質彈性力學。Kittel 把它們併在一起,理由要到下一個區塊才講得清楚。本段只負責 p.47–52,也就是「章首導論」加上「惰性氣體晶體」剛開頭的那一頁正文(p.49),以及緊接著的三張整頁週期表(p.50 內聚能、p.51 熔點、p.52 體模量與壓縮率)。
為什麼「鍵結」和「彈性常數」會被塞進同一章我的判讀
書上沒有明說,但只要把兩件事擺在一起就看得出來:鍵結決定的是一條 U(R) 曲線——晶體總能量對原子間距 R 的函數;而彈性常數量的是同一條曲線在極小值附近的曲率。
用一維的話說:內聚能是 U 在最低點的深度,平衡晶格常數是最低點的位置,彈性常數(在各向同性的情形就是體模量)是最低點的曲率。三者是同一條曲線的三個問題,所以放同一章。這也解釋了為什麼章首導論會在同一段話裡同時提到 Table 1(內聚能)、Table 2(熔點)與 Table 3(體模量)——它們全是同一條 U(R) 的不同讀數。
這個觀點在本段最後那個「由內聚能估體模量」的補推導裡會被量化:只要知道每個原子佔多大體積、鍵有多深,就能把體模量估到只差 1.4 到 7 倍——而這批材料的體模量本身橫跨了 443 倍,這已經是很強的結論了。
唯一的問題與唯一的答案:全部都是靜電(書 p.49 第一段)
原文用一句話定調整章:“In this chapter we are concerned with the question: What holds a crystal together?” 然後立刻給出答案,而且這個答案的語氣非常強硬:“The attractive electrostatic interaction between the negative charges of the electrons and the positive charges of the nuclei is entirely responsible for the cohesion of solids.”(粗體為我所加)——電子的負電荷與原子核的正電荷之間的吸引性靜電交互作用,完全負責固體的內聚。
接著把其他候選者一次排除:“Magnetic forces have only a weak effect on cohesion, and gravitational forces are negligible.” 磁力對內聚只有微弱影響,重力可以忽略。
最後一句是整段最容易被讀者滑過、但其實最重要的一句:“Specialized terms categorize distinctive situations: exchange energy, van der Waals forces, and covalent bonds.” ——交換能、凡得瓦力、共價鍵這些專有名詞,只是替不同情境貼的分類標籤,不是不同的基本作用力。
這句話的份量在於:初學者讀完固態物理的鍵結章,常常會在腦中建立起一張「五種力」的清單(凡得瓦力、離子鍵、共價鍵、金屬鍵、氫鍵),彷彿宇宙提供了五種膠水。Kittel 在這裡明白否定這件事。所有的膠水都是庫侖力,差別只在於電子在空間中怎麼分布、量子力學(包括包立原理與交換效應)如何限制那個分布。這是接下來整章推導的立論基礎:p.53 起的凡得瓦-倫敦交互作用,會從兩個諧振子的純庫侖耦合項出發;p.60 起的馬德倫能,會從 ±q²/r 的純庫侖求和出發。兩者用的是同一條定律。
把「磁力很弱、重力可忽略」換算成數字(書上只給定性判斷)補推導
原文只寫了 “Magnetic forces have only a weak effect on cohesion, and gravitational forces are negligible.”,沒有給任何數字。可是「弱」和「可忽略」差了多少數量級,其實是理解這一章的第一個標尺——因為它告訴你:如果你算出來的鍵能有幾個 eV,那它一定來自庫侖,不可能來自別的地方。以下用最簡單的兩個估算補上。
算完之後值得注意兩件事。第一,磁偶極—磁偶極作用對應的溫度尺度只有 0.1 K 上下,可是 Fe、Co、Ni 的居里溫度都是幾百到一千 K——差了四個數量級。這正是磁性有序不可能由磁偶極作用造成的證明,真正的來源是交換作用,而交換作用是庫侖作用加上包立原理的產物。第二,重力小到 10⁻³⁶,這是「固態物理完全不需要寫進重力」的量化理由。
- 重力 vs 庫侖(兩個質子):比值 = G m_p² / (k_e e²)。代入 G = 6.674×10⁻¹¹ N·m²/kg²、m_p = 1.673×10⁻²⁷ kg、k_e = 8.988×10⁹ N·m²/C²、e = 1.602×10⁻¹⁹ C。
- 分子 = 6.674×10⁻¹¹ × (1.673×10⁻²⁷)² = 6.674×10⁻¹¹ × 2.799×10⁻⁵⁴ = 1.868×10⁻⁶⁴。分母 = 8.988×10⁹ × (1.602×10⁻¹⁹)² = 8.988×10⁹ × 2.566×10⁻³⁸ = 2.307×10⁻²⁸。
- 比值 = 1.868×10⁻⁶⁴ / 2.307×10⁻²⁸ = 8.1×10⁻³⁷。若改用兩個電子(質量 9.109×10⁻³¹ kg),比值更小:2.4×10⁻⁴³。這就是 negligible 的量化意思。
- 磁偶極 vs 庫侖:兩個波耳磁子 μ_B = 9.274×10⁻²⁴ J/T,相距 r = 2 Å = 2×10⁻¹⁰ m,偶極—偶極能量尺度 U_dd = (μ₀/4π)·μ_B²/r³ = 10⁻⁷ × (9.274×10⁻²⁴)² / (2×10⁻¹⁰)³。
- = 10⁻⁷ × 8.601×10⁻⁴⁷ / 8×10⁻³⁰ = 10⁻⁷ × 1.075×10⁻¹⁷ = 1.075×10⁻²⁴ J = 6.7×10⁻⁶ eV,換成溫度是 1.075×10⁻²⁴ / 1.381×10⁻²³ = 0.078 K。
- 同距離下的庫侖尺度:e²/(4πε₀ r) = 14.4 eV·Å / 2 Å = 7.2 eV。兩者相差 7.2 / 6.7×10⁻⁶ ≈ 1.1×10⁶ 倍。
- 對照 Table 1(書 p.50)的內聚能範圍 0.020 eV(Ne)到 8.90 eV(W):磁偶極作用比最弱的惰性氣體鍵結還要小 3000 倍,比最強的 W 小六個數量級。weak effect 名副其實。
Fig. 1:四型鍵結的差別全出在「最外層電子與離子實的電荷分布」(書 p.48 Fig. 1;引句出自 p.49)
原文把整章的立論壓縮成一句:“The observed differences between the forms of condensed matter are caused in the final analysis by differences in the distribution of the outermost electrons and the ion cores (Fig. 1).” ——凝聚態形式之間觀察到的差異,歸根究柢是最外層電子與離子實的分布不同造成的。
Fig. 1 的圖說逐字如下:“Figure 1 The principal types of crystalline binding. In (a) neutral atoms with closed electron shells are bound together weakly by the van der Waals forces associated with fluctuations in the charge distributions. In (b) electrons are transferred from the alkali atoms to the halogen atoms, and the resulting ions are held together by attractive electrostatic forces between the positive and negative ions. In (c) the valence electrons are taken away from each alkali atom to form a communal electron sea in which the positive ions are dispersed. In (d) the neutral atoms are bound together by the overlapping parts of their electron distributions.”
【Fig. 1 的畫面:四個灰底方框,網頁不放圖,以下為目視描述】
(a) Crystalline argon (van der Waals):中灰底,七個大小完全相同的白色圓,各標 Ar,排成一個中央一顆、外圍六顆的緊密堆積圖樣。圓與圓相切但不重疊。畫面要傳達的是:每顆原子都是一顆閉殼層的硬球,彼此只是靠上去,電子雲沒有被對方搶走,也沒有伸進對方裡面。鍵結來自圖說裡那個關鍵詞——電荷分布的漲落(fluctuations),也就是瞬時偶極,而不是任何永久的電荷不平衡。
(b) Sodium chloride (ionic):中灰底,大圓標 Cl⁻、小圓標 Na⁺,排成棋盤式交錯:上排 Cl⁻、Na⁺、Cl⁻,中排 Na⁺、(大顆)Cl⁻、Na⁺,下排 Cl⁻、Na⁺、Cl⁻。圓仍然相切不重疊。這張圖同時傳達兩件事:電荷被整顆搬走了(所以標 ±),以及離子半徑因此劇烈改變——Na 失去一顆電子後縮成小球,Cl 得到一顆電子後脹成大球。這個大小反轉會在 p.70 的 ATOMIC RADII 一節被量化。
(c) Sodium (metallic):這格的底色明顯比其他三格深,只畫了五個很小的白圓標 Na⁺(四角加中央),彼此離得很開、完全不接觸。深底色本身就是圖:它代表圖說裡的 communal electron sea(共有電子海)。價電子不再屬於任何一顆原子,而是瀰漫整格;留下的 Na⁺ 離子實很小,稀疏地「泡」在電子海裡。這張圖與 (a) 的對比極為刻意——(a) 是大球相切、底色淺;(c) 是小球疏離、底色深。
(d) Diamond (covalent):底色與 (a)、(b) 完全相同(四格之中只有 (c) 明顯加深),五個大白圓標 C(四角加中央),而且中央那顆與四角那顆明顯互相重疊,重疊區域的輪廓線都畫出來了。這就是圖說的 “the overlapping parts of their electron distributions”——電子沒有被搬走(所以不標電荷),也沒有變成公有(所以不畫深底),而是兩顆原子的電子分布在鍵的位置上疊起來、共用。
把四格並排看,Kittel 的論證就成立了:四張圖用的畫筆只有三種——圓的大小(電荷有沒有轉移)、圓有沒有重疊(電子分布有沒有共用)、底色(電子有沒有公有化)。關於底色要特別註明一點:它只有深/淺兩態,不是四級灰階。(a)、(b)、(d) 三格的底色完全相同,只有金屬的 (c) 那一格明顯加深——也就是說這支畫筆是二元的,專門用來把「價電子公有化」的金屬與其餘三型分開,並不區分凡得瓦、離子、共價三者。三個變數的不同組合就生出四型鍵結,而背後的力自始至終只有庫侖。
把 Fig. 1 整理成一張對照表——以及一個必須先打的預防針我的判讀
下表把 Fig. 1 的四格拆成三個可辨識的變數,右欄補上它們各自在本章的哪一節被量化。表中「本章哪裡處理」欄的頁碼取自章首目錄頁(書 p.47),已於上一個區塊逐條核對。
| 類型(Fig. 1 標號) | 最外層電子去了哪裡 | 離子實帶電嗎 | 電子分布有無重疊/公有 | 本章哪裡處理 |
|---|---|---|---|---|
| 凡得瓦(a,固態 Ar) | 留在原子上,只是瞬時漲落 | 中性 | 幾乎不重疊 | CRYSTALS OF INERT GASES,p.49 起 |
| 離子(b,NaCl) | 整顆轉移給鹵素 | 帶 ±q | 不重疊(閉殼層) | IONIC CRYSTALS,p.60 起 |
| 金屬(c,Na) | 被剝離成公有電子海 | 帶 + | 公有化 | METALS,p.69 |
| 共價(d,金剛石) | 留在原子間,兩顆共用 | 中性 | 明顯重疊 | COVALENT CRYSTALS,p.67 |
預防針:這四型是一個連續譜,不是四個互斥的抽屜。Fig. 1 畫得乾淨俐落,很容易讓人以為每個晶體都能被分進其中一格。實際上:GaAs 的鍵同時有共價與離子成分;過渡金屬的鍵同時有金屬性與方向性的 d 電子共價成分;就連離子晶體裡也有凡得瓦貢獻——這一點 Kittel 自己在書 p.61 講得很明白:“The van der Waals part of the attractive interaction in ionic crystals makes a relatively small contribution to the cohesive energy in ionic crystals, of the order of 1 or 2 percent.”
所以正確的讀法是:Fig. 1 給的是四個極端情形,是座標系的四個端點,真實材料落在它們張出的空間裡。這也是為什麼本章要把四型分開處理——只有在極端情形下,模型才簡單到可以解析求解。
定義一:cohesive energy(內聚能)——逐字釘死(書 p.49)
原文的定義句是:“The cohesive energy of a crystal is defined as the energy that must be added to the crystal to separate its components into neutral free atoms at rest, at infinite separation, with the same electronic configuration.”
翻成中文並保留每一個限定詞:晶體的內聚能,定義為「為了把晶體的組成分離成靜止的、無限遠分開的、具有相同電子組態的中性自由原子,所必須加給晶體的能量」。
注意三件事。第一,它是加進去的能量(energy that must be added),所以是正值;很多書寫成「晶體比自由原子低多少能量」,符號相反,物理相同。第二,終態被三個條件完全限定死:at rest(靜止,沒有動能)、at infinite separation(無限遠,沒有殘餘交互作用)、with the same electronic configuration(相同電子組態,沒有電子激發)。第三,也是最關鍵的——終態是 neutral free atoms(中性自由原子)。
這第三點就是下一個定義存在的全部理由。
定義二:lattice energy(晶格能)——只用在離子晶體,終態是自由離子(書 p.49)
原文緊接著寫:“The term lattice energy is used in the discussion of ionic crystals and is defined as the energy that must be added to the crystal to separate its component ions into free ions at rest at infinite separation.”
晶格能這個詞用在離子晶體的討論裡,定義為「為了把晶體的組成離子分離成靜止、無限遠分開的自由離子,所必須加給晶體的能量」。
兩個定義唯一的差別,就是終態選了什麼參考態:
| cohesive energy | lattice energy | |
|---|---|---|
| 終態 | 中性自由原子 | 自由離子(帶電) |
| 適用範圍 | 所有晶體(原文語氣是通用定義) | 原文明說 used in the discussion of ionic crystals |
| 對 NaCl 的實驗值 | 6.4 eV/分子單位(書 p.61 Fig. 8) | 7.9 eV/分子單位(書 p.60 正文與 p.61 Fig. 8) |
這條分界線是本段最該記住的一件事。如果你把 7.9 eV 當成 NaCl 的內聚能寫進報告,或把 6.4 eV 拿去跟馬德倫能比較,後面整段離子晶體的推導都會對不上。下一個區塊把這兩個數字怎麼換算補出來。
定義裡每一個限定詞為什麼都不能少我的判讀
教科書的定義句常常被當成廢話讀過去。這一句不行——它的每個限定詞都對應到一個會算錯的情境。
(1)at rest(靜止)。如果終態的自由原子還帶著動能,你量到的「分離所需能量」就會偏小。更重要的是,這個限定詞在初態那一側留了一個後門:晶體裡的原子在 0 K 時不是靜止的,它們有零點運動。所以「實驗內聚能」與「把原子固定在格點上算出來的理論內聚能」之間必然有一個差額,那個差額就是零點能。這正是 p.49 那條腳註、以及 p.59 那串 28%/10%/6%/4% 的量子修正在處理的事。
(2)at infinite separation(無限遠)。終態不能有殘餘交互作用。對凡得瓦晶體特別要小心:−C/R⁶ 衰減得慢,把「無限遠」近似成「幾個奈米」在數值計算裡會累積可觀誤差。
(3)with the same electronic configuration(相同電子組態)。終態的自由原子必須處在與晶體中相同的電子組態——實務上就是基態。Table 1 的圖說把這件事寫得更死:“Energy required to form separated neutral atoms in their ground electronic state from the solid at 0 K at 1 atm.”(書 p.50;該頁的標題與圖說有文字層,可逐字信任,只有表格數值才是目視判讀)。沒有這個限定,你可以把原子拆成任何激發態,答案就不唯一了。
(4)neutral(中性)。這是與 lattice energy 唯一的分界。它之所以重要,是因為離子晶體的自然計算終點是離子——馬德倫求和 ±q²/r 算出來的東西,天生是對自由離子定義的。要轉成內聚能,必須額外付出(或回收)游離能與電子親和能。下一個式子就是這筆帳。
內聚能與晶格能的換算(以 NaCl 為例)
E_coh = U_latt − I(Na) + A(Cl) = 7.9 − 5.1 + 3.6 = 6.4 eV(每分子單位)
p.61E_coh = 相對於分離的中性原子的內聚能;U_latt = 相對於分離的自由離子的晶格能;I(Na) = 鈉的游離能,書 p.61 Fig. 8 圖上標 5.14 eV,圖說取近似值 5.1;A(Cl) = 氯的電子親和能,Fig. 8 圖上標 3.61 eV,圖說取近似值 3.6。
Table 1 的原文結論:週期表上的四段落差(書 p.49–50)
原文對 Table 1 只寫了四句話,但這四句話是本章的全局座標:
“Values of the cohesive energy of the crystalline elements are given in Table 1. Notice the wide variation in cohesive energy between different columns of the periodic table. The inert gas crystals are weakly bound, with cohesive energies less than a few percent of the cohesive energies of the elements in the C, Si, Ge . . . column. The alkali metal crystals have intermediate values of the cohesive energy. The transition element metals (in the middle columns) are quite strongly bound.”
拆成四條可檢驗的陳述:
(一)不同縱行之間差異很大。
(二)惰性氣體晶體鍵結很弱,內聚能不到 C、Si、Ge…那一行的百分之幾。
(三)鹼金屬晶體居中。
(四)中間各行的過渡金屬鍵結相當強。
用 Table 1 的數字逐條核對(全部為整頁影像目視判讀,書 p.50):
(二)Ar 0.080 eV/C 7.37 eV = 1.1%;Ne 0.020/Si 4.63 = 0.43%;Kr 0.116/Ge 3.85 = 3.0%;Xe 0.16/C 7.37 = 2.2%。原文說 less than a few percent,最大的 Kr/Ge 是 3.0%,勉強貼著上緣,其餘都在 2% 以內——成立。
(三)Na 1.113、K 0.934、Rb 0.852、Cs 0.804、Li 1.63 eV。確實夾在惰性氣體(0.02–0.2)與 C/Si/Ge(3.85–7.37)之間——成立。
(四)中間各行最大值:W 8.90、Os 8.17、Ta 8.10、Re 8.03、Nb 7.57、Ir 6.94 eV。W 的 8.90 eV 是整張表的最大值,比 C 的 7.37 eV 還高——成立,而且比原文的語氣更強。
順帶一提,鹼金屬那一行由上往下是 Li 1.63 → Na 1.113 → K 0.934 → Rb 0.852 → Cs 0.804,單調下降;惰性氣體那一行由上往下是 Ne 0.020 → Ar 0.080 → Kr 0.116 → Xe 0.16 → Rn 0.202,單調上升。兩條相反的趨勢並存在同一張表上,其物理根源完全不同:鹼金屬是價電子離核越來越遠、金屬鍵越來越鬆;惰性氣體是原子極化率隨原子序增大、凡得瓦力越來越強。後者到 p.53 的倫敦色散推導會自然浮現。
Table 1 的判讀說明與一個缺欄(書 p.50,整頁影像目視判讀)我的判讀
本區塊全部數值皆為整頁影像目視判讀。書本 p.50 在這份 PDF 裡沒有文字層,只有標題與圖說抽得出來;表格本體是一整張影像,以下數值是把該頁以 3 倍解析度算繪後放大目視讀出的。
版面:每個元素格子有三列數字,由上而下依序是 kJ/mol、eV/atom、kcal/mol——這三行單位標籤印在表格中央的三條橫向箭頭上,不是印在欄首,很容易漏看。中間那列(eV/atom)有灰色網底,這是版面用來凸顯物理學家最常用單位的手法。
一個必須指出的缺欄:Table 1 沒有 H 與 He。左上角第一格直接從 Li 開始,右上角第一格直接從 B 開始,整個第一週期(H、He)不存在。這不是印刷缺陷,而是有物理理由的:氦在零壓下根本不固化(見 p.49 腳註),沒有「固體的內聚能」可言;氫的固體是分子晶體 H₂ 而非原子晶體,用「每原子內聚能」描述會產生歧義。對照之下,Table 3(p.52)就有 H 與 He 兩格(H 的體模量 0.002、He 的 0.00),因為那裡量的是加壓後的凝聚相。這種「同一章的三張表收錄範圍不同」的細節,在引用時一定要注意。
本段會反覆用到的幾筆(單位:kJ/mol、eV/atom、kcal/mol):
| 元素 | kJ/mol | eV/atom | kcal/mol |
|---|---|---|---|
| Ne | 1.92 | 0.020 | 0.46 |
| Ar | 7.74 | 0.080 | 1.85 |
| Kr | 11.2 | 0.116 | 2.68 |
| Xe | 15.9 | 0.16 | 3.80 |
| Rn | 19.5 | 0.202 | 4.66 |
| Li | 158. | 1.63 | 37.7 |
| Na | 107. | 1.113 | 25.67 |
| Cs | 77.6 | 0.804 | 18.54 |
| Al | 327. | 3.39 | 78.1 |
| Cu | 336. | 3.49 | 80.4 |
| Si | 446. | 4.63 | 106.7 |
| Ge | 372. | 3.85 | 88.8 |
| Fe | 413. | 4.28 | 98.7 |
| C | 711. | 7.37 | 170. |
| W | 859. | 8.90 | 205.2 |
(表中若干格子原書印成「711.」「170.」這種尾隨句點的形式,我照原樣保留,那個點是有效位數的標記而非小數點。)
還有一筆對不上的地方值得記下來:Table 1 給 Ne 的內聚能是 1.92 kJ/mol、Xe 是 15.9 kJ/mol,但書 p.53 的 Table 4(惰性氣體晶體性質,有文字層可抽)給的是 Ne 1.88 kJ/mol、Xe 16.0 kJ/mol(eV 欄則是 0.02 與 0.17)。Ar 的 7.74/0.080 與 Kr 的 11.2/0.116 兩表完全一致。這是兩張表資料來源不同造成的小差異(Table 1 註明資料由 Prof. Leo Brewer 提供),差距約 2%(Ne 是 0.04/1.88 = 2.1%,Xe 只有 0.6%),不影響任何結論,但引用時要說清楚你用的是哪一張表。
內聚能的三種單位換算
1 eV/atom = 96.485 kJ/mol = 23.061 kcal/mol;1 kcal = 4.184 kJ
p.50換算因子 96.485 kJ/(mol·eV) = N_A·e = 6.02214×10²³ mol⁻¹ × 1.602177×10⁻¹⁹ C;23.061 = 96.485 / 4.184。
Table 2(熔點):原文只說 roughly,我們來量「roughly」有多 rough(書 p.51,整頁影像目視判讀)我的判讀
書本 p.51 同樣是一整頁表格影像,文字層只抽得到標題:“Table 2 Melting points, in K. (After R. H. Lamoreaux)”。以下數值為算繪後放大目視判讀。
原文對這張表只有半句話:“The melting temperatures (Table 2) and bulk modulii (Table 3) vary roughly as the cohesive energies.”(書 p.49;注意 Kittel 在這裡把 bulk moduli 拼成 modulii,這是原書的拼法,引用時照抄。)
把 Table 2 的熔點除以 Table 1 的內聚能(eV/atom),得到一個「每 eV 值多少 K 熔點」的比值:
| 元素 | T_m (K) | E_coh (eV) | T_m/E_coh (K/eV) |
|---|---|---|---|
| Ne | 24.56 | 0.020 | 1228 |
| Ar | 83.81 | 0.080 | 1048 |
| Kr | 115.8 | 0.116 | 998 |
| Xe | 161.4 | 0.16 | 1009 |
| Li | 453.7 | 1.63 | 278 |
| Na | 371.0 | 1.113 | 333 |
| K | 336.3 | 0.934 | 360 |
| Cs | 301.6 | 0.804 | 375 |
| Al | 933.5 | 3.39 | 275 |
| Cu | 1358 | 3.49 | 389 |
| Si | 1687 | 4.63 | 364 |
| Ge | 1211 | 3.85 | 315 |
| Fe | 1811 | 4.28 | 423 |
| W | 3695 | 8.90 | 415 |
結論很清楚:內聚能橫跨 445 倍(0.020 → 8.90 eV)、熔點橫跨 150 倍(24.56 → 3695 K),而兩者的比值只在 275 到 1228 之間變動,只有 4.5 倍。這就是 roughly 的量化意義——它不是敷衍,而是一個橫跨三個數量級仍然成立的粗略正比關係。
更有意思的是這 4.5 倍怎麼分布:惰性氣體全部落在 998–1228,其餘全部落在 275–423,兩群之間幾乎不重疊。也就是說,在同一種鍵結類型內部,比值幾乎是常數;跨鍵結類型才會跳一階。這個現象的物理是:熔化比較的是固相與液相的自由能,而液相保留多少鍵能與鍵的方向性有關;凡得瓦鍵完全沒有方向性,熔化時破壞的比例小,所以相對來說要更高的 T/E_coh 才熔得掉。
幾個判讀細節:(1)Table 2 與 Table 1 一樣沒有 H 與 He 兩格。(2)C 的格子是空的——表上有 C 的方框與元素符號,但沒有印任何數字。(3)P 的格子印了兩個數字,前面各帶一個斜體小字母:w 317 與 r 863,顯然是白磷與紅磷兩種晶形。(4)Rn、At、Fr 三格空白。(5)Lu 下方那格印成 Lw(Lawrencium 的舊符號),現行符號是 Lr。
Table 3(體模量與壓縮率):三張表裡資訊量最大的一張(書 p.52,整頁影像目視判讀)我的判讀
書本 p.52 在這份 PDF 裡連標題都沒有文字層——pypdf 對該頁只抽得到一個「52」。以下全部內容來自算繪後目視判讀。
標題與圖說逐字(目視):“Table 3 Isothermal bulk modulii and compressibilities at room temperature / After K. Gschneidner, Jr., Solid State Physics 16, 275–426 (1964); several data are from F. Birch, in Handbook of physical constants, Geological Society of America Memoir 97, 107–173 (1966). Original references should be consulted when values are needed for research purposes. Values in parentheses are estimates. Letters in parentheses refer to the crystal form. Letters in brackets refer to the temperature: [a] = 77 K; [b] = 273 K; [c] = 1 K; [d] = 4 K; [e] = 81 K.”
單位(印在表格中央兩條橫線上):Bulk modulus in units 10¹² dyn/cm² or 10¹¹ N/m²;Compressibility in units 10⁻¹² cm²/dyn or 10⁻¹¹ m²/N。每格上列是體模量 B,下列(灰底)是壓縮率 κ。
三個讀表時必須注意的陷阱:
(1)標題寫 at room temperature,但惰性氣體那幾格根本不可能在室溫。He、Ne、C 標 [d] = 4 K;Ar、Kr 標 [a] = 77 K;N 標 [e] = 81 K;Hg 標 [c] = 1 K;Ga 標 [b] = 273 K。方括號是溫度、圓括號是晶形,這兩種括號長得很像但意思完全不同,而原書偏偏把同一個字母 b 的兩種用法排在相鄰的兩列:P 印的是 (b)(圓括號,依圖說的規則讀作晶形,即黑磷),Ga 印的是 [b](方括號,273 K)。把 Ga 誤讀成晶形註記,正是這張表最典型的陷阱。其餘圓括號的晶形註記還有 S 的 (r)、Sn 的 (g)、Ce 的 (γ)。
(2)圓括號裡的數值是估計值,例如 Os 的 (4.18)/(0.24)、Fr 的 (0.020)/(50.)、Pm 的 (0.35)/(2.85)。引用估計值要註明。
(3)這張表有 H 與 He,而 Table 1、Table 2 沒有。He [d] 的體模量印成 0.00、壓縮率 1168;H [d] 是 0.002 / 500。
本段會用到的幾筆(B 單位 10¹¹ N/m²,κ 單位 10⁻¹¹ m²/N):
| 元素 | B | κ | 附註 |
|---|---|---|---|
| He | 0.00 | 1168 | [d] = 4 K |
| H | 0.002 | 500 | [d] = 4 K |
| Ne | 0.010 | 100 | [d] = 4 K |
| Ar | 0.013 | 79 | [a] = 77 K |
| Kr | 0.018 | 56 | [a] = 77 K |
| Xe | —(空白) | — | |
| Cs | 0.020 | 50. | |
| K | 0.032 | 31. | |
| Na | 0.068 | 14.7 | |
| Li | 0.116 | 8.62 | |
| Ge | 0.772 | 1.29 | |
| Al | 0.722 | 1.385 | |
| Si | 0.988 | 1.012 | |
| Cu | 1.37 | 0.73 | |
| Fe | 1.683 | 0.594 | |
| W | 3.232 | 0.309 | |
| C | 4.43 | 0.226 | [d] = 4 K |
He 那格的 0.00 是這一整段最有故事的一個數字。它不是「氦沒有體模量」,而是四捨五入的結果——下一個區塊會從 κ = 1168 把它還原出來。而它之所以小到印不出來,原因就寫在 p.49 的那條腳註裡。
補算 Table 3 的內部一致性:κ = 1/B,順便把 He 的「0.00」還原成真正的數字補推導
書上沒有提醒讀者 Table 3 的兩列數字並非獨立——壓縮率就是體模量的倒數,κ ≡ −(1/V)(∂V/∂p) = 1/B。這條關係給了我們一個免費的判讀檢查:如果我對影像的目視判讀有錯,兩列就對不起來。
單位要先對齊:B 的單位是 10¹¹ N/m²,κ 的單位是 10⁻¹¹ m²/N,兩者恰好互為倒數的單位,所以只要檢查「上列數字 × 下列數字 = 1」就好。這個單位搭配顯然是原書刻意設計的。
算完之後,最有價值的副產品是 He 那一格:它證明了「0.00」是四捨五入而非真零,而且還原出來的 8.6×10⁷ N/m²(約 856 bar)本身就是一個有物理意義的數字——它比金剛石小了約三個半數量級(5175 倍),這是全表最大的落差。
- Ne:B × κ = 0.010 × 100 = 1.00 ✔
- Ar:0.013 × 79 = 1.027(判讀值只有兩位有效數字,1/79 = 0.01266 四捨五入即 0.013)✔
- Kr:0.018 × 56 = 1.008 ✔
- Li:0.116 × 8.62 = 1.000 ✔ Na:0.068 × 14.7 = 1.000 ✔ Cs:0.020 × 50. = 1.00 ✔
- Si:0.988 × 1.012 = 1.000 ✔ Cu:1.37 × 0.73 = 1.000 ✔ W:3.232 × 0.309 = 0.999 ✔ C:4.43 × 0.226 = 1.001 ✔
- 十組全部通過,誤差都在四捨五入範圍內。這同時驗證了我對 p.52 整頁影像的判讀是正確的。
- 還原 He:由 κ = 1168 × 10⁻¹¹ m²/N,得 B = 1/(1.168×10⁻⁸ m²/N) = 8.56×10⁷ N/m² = 0.000856 × 10¹¹,四捨五入到兩位小數正是表上印的 0.00。
- 還原 H:κ = 500 × 10⁻¹¹ → B = 1/(5.00×10⁻⁹) = 2.00×10⁸ = 0.002 × 10¹¹,與表上印的 0.002 完全一致 ✔
- 落差有多大:C 的 4.43 對 He 的 0.000856,比值 5175 倍;壓縮率則是 1168/0.226 = 5168 倍(兩者應相等,差異來自四捨五入)。整張週期表的剛性跨越三個半數量級,而最軟的那一格正是氦。
補出「熔點與體模量大致隨內聚能變化」背後的量綱論證:B ≈ E_coh / V_atom補推導
原文那句 “The melting temperatures (Table 2) and bulk modulii (Table 3) vary roughly as the cohesive energies.” 只是一個觀察陳述,書上沒有給任何理由。可是理由其實只有一行量綱分析,而且做出來的預測準到令人意外,值得補上。
設定:體模量的定義是 B = V(∂²U/∂V²),量綱是「能量/體積」。晶體裡唯一的能量尺度是每原子的內聚能 E_coh,唯一的長度尺度是原子佔的體積 V_atom。所以量綱上必然有
B = c · E_coh / V_atom,其中 c 是一個無量綱數。
量綱分析能做到的就到這裡——它不能告訴你 c 是多少,只能說「如果內聚能與原子體積是唯一的兩個尺度,那 c 應該是 O(1)」。下面用實際資料把 c 算出來,看看它是不是真的 O(1)、以及它變動多大。
資料來源:E_coh 取自書 p.50 Table 1 的 eV/atom 欄(整頁影像目視判讀);B 取自書 p.52 Table 3(整頁影像目視判讀);晶格常數取自書 p.20 Table 3「Crystal structures of the elements」(同樣是整頁影像,已算繪目視判讀)。單位換算用 1 eV/ų = 1.602×10¹¹ N/m²,剛好就是 Table 3 的單位,很方便。
結果(見下方步驟):c 落在 1.4 到 6.9 之間。在 E_coh/V_atom 橫跨 1440 倍、B 橫跨 443 倍的資料上,c 只變動 5 倍。這就是 “vary roughly as” 的定量內容。
c 為什麼不是 1、而且為什麼惰性氣體偏大?因為 c 其實藏著位能曲線在極小值附近的形狀資訊。對 Lennard-Jones 這種「內側極陡、外側極平」的位能(書 p.58 Fig. 6 的圖說就是這麼描述的:“Notice how steep the curve is inside the minimum, and how flat it is outside the minimum.”),曲率相對於井深偏大,所以 c 偏大;金屬的位能曲線相對平緩,c 偏小。另外惰性氣體那三筆的 B 是在 4 K 或 77 K 量的,而 Ar 的 77 K 已經是它熔點 83.81 K 的 92%、Kr 的 77 K 則是熔點 115.8 K 的 66%,兩者都有程度不等的熱軟化,而熱軟化會壓低 B,所以真實的 c 應該比算出來的更大——這反而讓「惰性氣體 c 偏大」的結論更穩固。
- 換算因子:1 eV/ų = 1.602177×10⁻¹⁹ J / 10⁻³⁰ m³ = 1.602×10¹¹ N/m²,即 1.602(以 Table 3 的 10¹¹ N/m² 為單位)。所以 B_pred = 1.602 × E_coh[eV] / V_atom[ų]。
- 每原子體積:fcc 與 diamond 的立方晶胞邊長為 a 時,fcc 每胞 4 個原子 → V_atom = a³/4;diamond 每胞 8 個 → a³/8;bcc 每胞 2 個 → a³/2。
- Ne(fcc, a = 4.46 Å @4K):V = 4.46³/4 = 88.72/4 = 22.18 ų;B_pred = 1.602×0.020/22.18 = 0.00145。實測 0.010 → c = 6.9
- Ar(fcc, a = 5.31 Å @4K):V = 149.72/4 = 37.43 ų;B_pred = 1.602×0.080/37.43 = 0.00342。實測 0.013 → c = 3.8。(附帶驗證:a/√2 = 5.31/1.4142 = 3.755 Å,與書 p.53 Table 4 印的最近鄰距離 3.76 Å 吻合,說明我對 p.20 的判讀正確。)
- Kr(fcc, a = 5.64 Å @4K):V = 179.41/4 = 44.85 ų;B_pred = 1.602×0.116/44.85 = 0.00414。實測 0.018 → c = 4.3
- Cs(bcc, a = 6.045 Å @5K):V = 220.90/2 = 110.45 ų;B_pred = 1.602×0.804/110.45 = 0.01166。實測 0.020 → c = 1.7
- Na(bcc, a = 4.225 Å @5K):V = 75.41/2 = 37.71 ų;B_pred = 1.602×1.113/37.71 = 0.04728。實測 0.068 → c = 1.4
- Ge(diamond, a = 5.658 Å):V = 181.14/8 = 22.64 ų;B_pred = 1.602×3.85/22.64 = 0.2724。實測 0.772 → c = 2.8
- Si(diamond, a = 5.430 Å):V = 160.10/8 = 20.01 ų;B_pred = 1.602×4.63/20.01 = 0.3706。實測 0.988 → c = 2.7
- Cu(fcc, a = 3.61 Å):V = 47.05/4 = 11.76 ų;B_pred = 1.602×3.49/11.76 = 0.4754。實測 1.37 → c = 2.9
- W(bcc, a = 3.16 Å):V = 31.55/2 = 15.78 ų;B_pred = 1.602×8.90/15.78 = 0.9037。實測 3.232 → c = 3.6
- C(diamond, a = 3.567 Å):V = 45.39/8 = 5.673 ų;B_pred = 1.602×7.37/5.673 = 2.081。實測 4.43 → c = 2.1
- 總結:E_coh/V_atom 從 Ne 的 0.0009 eV/ų 到 C 的 1.30 eV/ų(1440 倍),B 從 0.010 到 4.43(443 倍),而 c 只在 1.4 – 6.9(5 倍)之間。量綱分析抓對了主導行為。
CRYSTALS OF INERT GASES:七句話(含頁尾被切斷的半句)立下的全部前提(書 p.49)
這一節的開頭極短,但每一句都是後面整整十頁推導所依賴的前提。逐字抄錄並逐條拆解:
(一)“The inert gases form the simplest crystals.” 惰性氣體形成最簡單的晶體。——這是選它當第一個例子的理由:沒有電荷轉移、沒有方向性鍵、沒有公有電子,只剩下最純粹的一種效應。
(二)“The electron distribution is very close to that of the free atoms.” 電子分布非常接近自由原子的分布。——這一句是整個 Lennard-Jones 處理的合法性來源。正因為晶體裡的電子分布幾乎沒被扭曲,才可以用氣相測到的 ε 與 σ 去算固體,而書 p.58 明白宣告了這個好處:“Values of ε and σ given in Table 4 can be obtained from gas-phase data, so that calculations on properties of the solid do not involve disposable parameters.”(無可調參數!這在固態物理裡是罕見的奢侈。)
(三)“Their properties at absolute zero are summarized in Table 4.” 絕對零度下的性質整理在 Table 4——那張表在書 p.53,屬於下一段的範圍。
(四)“The crystals are transparent insulators, weakly bound, with low melting temperatures.” 這些晶體是透明絕緣體、鍵結很弱、熔點很低。——「透明絕緣體」這四個字是能隙的直接後果:外殼層填滿代表價帶全滿,能隙又寬到可見光(1.6–3.1 eV)打不動,所以既不導電也不吸光。這是第 7、8 章能帶論的一個現成註腳。
(五)“The atoms have very high ionization energies (see Table 5).” 原子有非常高的游離能。——Table 5 也在後面。這一句與(二)互為因果:游離能高 → 電子抓得緊 → 電子分布不易被鄰居扭曲。
(六)“The outermost electron shells of the atoms are completely filled, and the distribution of electron charge in the free atom is spherically symmetric.” 最外層電子殼層完全填滿,自由原子的電荷分布球對稱。——「球對稱」這四個字是 p.53 起那段推導的起手式:正因為球對稱,永久偶極矩為零,所以鍵結不可能來自靜態偶極,只能來自漲落。p.53 原文把這個邏輯寫得很清楚:“If the charge distributions on the atoms were rigid, the interaction between atoms would be zero, because the electrostatic potential of a spherical distribution of electronic charge is canceled outside a neutral atom by the electrostatic potential of the charge on the nucleus.”
(七)最後半句:“In the crystal the inert gas atoms pack together as closely as possible” ——在晶體中,惰性氣體原子盡可能緊密地堆積。這句話在 p.49 頁尾被切斷(後面接一個腳註記號 1),下半句在 p.53 的開頭。
跨頁的那半句話:fcc、以及「內聚能只有游離能的 1% 或更少」(接續至書 p.53)
書本 p.50、51、52 三頁全部是整頁表格,正文從 p.49 頁尾直接跳到 p.53。那個被切斷的句子完整讀起來是:
“In the crystal the inert gas atoms pack together as closely as possible: the crystal structures (Fig. 2) are all cubic close-packed (fcc), except He³ and He⁴.”(前半在 p.49,後半在 p.53)
——除了 He³ 與 He⁴ 以外,全部是立方最密堆積(fcc)。這個結構結論非常乾淨,而它的例外正好又是氦,指向同一個原因。(附帶對照:書 p.20 的 Table 3「Crystal structures of the elements」也印著 Ne fcc 4.46、Ar fcc 5.31、Kr fcc 5.64、Xe fcc 6.13,而 He⁴ 那格印的是 hcp 3.57/5.83 @2K,與這裡的「除了 He 以外都是 fcc」一致。以上為 p.20 整頁影像目視判讀。)
接著 p.53 問了整節的核心問題並給出量化前提:
“What holds an inert gas crystal together? The electron distribution in the crystal is not significantly distorted from the electron distribution around the free atoms because not much energy is available to distort the free atom charge distributions. The cohesive energy of an atom in the crystal is only 1 percent or less of the ionization energy of an atomic electron. Part of this distortion gives the van der Waals interaction.”
這段話的邏輯鏈值得慢慢讀:可用的能量太少(內聚能只有游離能的 1% 或更少)→ 所以電子分布扭曲不了多少 → 所以可以把原子當成幾乎剛性的球 → 而那一點點扭曲,就是凡得瓦交互作用的全部來源。注意最後一句 “Part of this distortion gives the van der Waals interaction.”——凡得瓦力不是「扭曲之外」的東西,它就是那一點點扭曲本身。這是 p.53 起用兩個耦合諧振子建模的物理依據。
把「1 percent or less」逐一驗算(書上沒給,而且 Xe 其實略微超標)補推導
原文說內聚能只有原子電子游離能的 1 percent or less,但沒有列出任何一個數字。這條陳述是整個「電子分布幾乎不被扭曲」論證的量化根基,值得親自查一次——尤其因為驗算之後會發現它不是每個元素都嚴格成立。
資料:內聚能取自書 p.50 Table 1 的 eV/atom 欄(整頁影像目視判讀);游離能取自書 p.53 Table 4 的 “Ionization potential of free atom, eV” 欄(該頁有文字層,可直接抽取)。
結論先講:Ne、Ar、Kr 都遠低於 1%,Xe 是 1.3%,略微超過。原文的 1 percent or less 應該理解為一個數量級陳述而非嚴格上界。這不影響論證——1.3% 與「幾乎不扭曲」的圖像完全相容——但如果你要把這個數字寫進報告,就該照實寫成「約 1% 或更低,最重的 Xe 約 1.3%」。
- Ne:E_coh = 0.020 eV,游離能 21.56 eV → 0.020/21.56 = 0.093%
- Ar:E_coh = 0.080 eV,游離能 15.76 eV → 0.080/15.76 = 0.51%
- Kr:E_coh = 0.116 eV,游離能 14.00 eV → 0.116/14.00 = 0.83%
- Xe:E_coh = 0.16 eV(Table 1)或 0.17 eV(Table 4),游離能 12.13 eV → 1.32% 或 1.40%,兩者都略高於 1%
- Rn(Table 1 有內聚能 0.202 eV,但 Table 4 沒收 Rn,無法計算)
- 趨勢的物理:由 Ne 到 Xe,比值單調上升約 14 倍。分子(內聚能)上升 8 倍是因為極化率隨原子變大而增加;分母(游離能)下降 1.8 倍是因為最外層電子離核越來越遠。兩個效應同向疊加,所以越重的惰性氣體越「不惰性」,電子分布被扭曲的相對程度越大。這解釋了為什麼 Xe 會形成化合物(XeF₄ 等)而 Ne 不會。
- 對照 He:Table 4 給 He 的游離能 24.58 eV,是全表最高,但 He 在 Table 1 裡沒有內聚能可填(它不固化)。所以這個比值對 He 根本定義不出來——這本身就是下一個區塊要處理的事。
p.49 的那條腳註:整段最重要的六行字(逐字全文)
書本 p.49 頁尾有一條標號 1 的腳註,掛在 “the inert gas atoms pack together as closely as possible” 那句話上。它的字體比正文小,很容易被跳過,但它其實是整段最有資訊密度的一塊,而且直接決定了 p.59、p.60 的量子修正該怎麼理解。逐字全文:
“Zero-point motion of the atoms (kinetic energy at absolute zero) is a quantum effect that plays a dominant role in He³ and He⁴. They do not solidify at zero pressure even at absolute zero temperature. The average fluctuation at 0 K of a He atom from its equilibrium position is of the order of 30 to 40 percent of the nearest-neighbor distance. The heavier the atom, the less important the zero-point effects. If we omit zero-point motion, we calculate a molar volume of 9 cm³ mol⁻¹ for solid helium, as compared with the observed values of 27.5 and 36.8 cm³ mol⁻¹ for liquid He⁴ and liquid He³, respectively.”
(校註:文字層把 temperature 一詞的換行連字號留在中間,抽出來是 “temp- erature”;已還原。上標的 3、4 在文字層是普通數字,已還原為 He³、He⁴。)
把它拆成六條獨立的陳述:
(1)零點運動 = 絕對零度下的動能,是量子效應。這一句給了定義。古典力學裡 0 K 就是所有原子靜止不動;量子力學不允許——把原子局限在一個有限的空間裡,測不準原理就強迫它有動量的不確定性,因而有動能。
(2)它在 He³ 與 He⁴ 裡扮演主導角色(dominant role)。注意用詞是 dominant,不是 important——它不是修正項,它比鍵能還大。
(3)它們在零壓下即使到絕對零度也不會固化。這是(2)的直接後果,也是全書唯一一個「基態是液體」的例子。
(4)He 原子在 0 K 偏離平衡位置的平均漲落,是最近鄰距離的 30 到 40 個百分點。這是一個令人吃驚的數字:原子晃動的幅度是鍵長的三分之一以上。作為對照,一般固體在熔點附近的均方根位移也只有鍵長的 10% 左右(這是 Lindemann 判據的內容,Kittel 在這裡沒提)。He 在絕對零度就已經「晃得比別的固體在熔點還厲害」。
(5)原子越重,零點效應越不重要。這句話給了整個惰性氣體序列的組織原則,也預告了 p.59 那串按質量遞減的量子修正 28%、10%、6%、4%。
(6)如果略去零點運動,算出固態氦的莫耳體積是 9 cm³/mol,而實測液態 He⁴、He³ 分別是 27.5 與 36.8 cm³/mol。這是腳註裡唯一的三個數字,也是唯一可以拿來驗算的地方。下面兩個區塊就做這件事。
補出腳註沒寫的算術之一:9 cm³/mol 是怎麼算出來的補推導
腳註只說 “If we omit zero-point motion, we calculate a molar volume of 9 cm³ mol⁻¹ for solid helium”,完全沒說怎麼算。可是這個「算」其實用的就是本章後面幾頁會建立的整套機器,在這裡先把它跑一遍,有兩個好處:一是驗證我們確實看懂了 Kittel 的模型,二是把 p.53、p.58、p.59 的三個結果串成一條可操作的流程。
需要的三個輸入:
(一)He 的 Lennard-Jones 參數 σ = 2.56 Å(書 p.53 Table 4,該頁有文字層)。注意 Table 4 給 He 的那一列只有三個數字——游離能 24.58 eV、ε = 14×10⁻¹⁶ erg、σ = 2.56 Å——最近鄰距離、內聚能、熔點三欄全部空白,表下還特別註明 “(liquid at zero pressure)”。Kittel 在表裡就已經誠實承認 He 沒有固體資料了。
(二)fcc 的平衡條件 R₀/σ = 1.09(書 p.59 式 (14))。原文說這個值 “the same for all elements with an fcc structure”,而且是純粹由 Lennard-Jones 位能對 fcc 晶格求和後極小化得到的,不含任何元素資訊。
(三)fcc 的幾何:晶胞邊長 a 與最近鄰距離 R₀ 的關係 a = √2 R₀,每個晶胞含 4 個原子。
為什麼算固態氦可以用 fcc?前一個區塊才剛引 p.53 說「除了 He³、He⁴ 以外全是 fcc」,這裡卻用 fcc 幾何算氦,看起來像自打嘴巴,所以要先交代:這個「假想的古典固態氦」本來就不存在,Kittel 算的是「若氦乖乖待在格點上會怎樣」,堆積方式只能照惰性氣體的通例取最密堆積。而且取 fcc 或 hcp 幾乎不影響答案——理想 hcp 與 fcc 都是最密堆積,在相同最近鄰距離下每原子體積完全相同;Lennard-Jones 的晶格求和常數兩者也幾乎一樣(書 p.59:fcc 為 12.13188/14.45392,hcp 為 12.13229/14.45489),所以 R₀/σ = 1.09 這個平衡條件照用無妨。
結果:9.25 cm³/mol,與 Kittel 的 9 吻合到 3%。這同時證明了一件事:腳註裡那個 9,就是「假裝氦是一顆乖乖待在格點上的古典原子」算出來的 Lennard-Jones 預測值。它與現實差了三倍,差距全部由零點運動買單。
一個容易混淆的地方要先澄清:注意不要用 Fig. 6 圖說裡的 R/σ = 2^(1/6) ≈ 1.12。那個 1.12 是兩顆孤立原子的 Lennard-Jones 極小值位置;晶體裡每顆原子有 12 個最近鄰再加上更遠的殼層,吸引項的求和係數(14.45)比排斥項(12.13)大,把平衡位置往內拉到 1.09。用 1.12 會算出 10.1 cm³/mol,與 9 差 12%,明顯偏離。
- 平衡最近鄰距離:R₀ = 1.09 × σ = 1.09 × 2.56 Å = 2.790 Å
- fcc 晶胞邊長:a = √2 × R₀ = 1.41421 × 2.790 = 3.946 Å
- 晶胞體積:a³ = 3.946³ = 61.45 ų
- 每原子體積:V_atom = a³/4 = 61.45/4 = 15.36 ų
- 莫耳體積:V_m = V_atom × N_A = 15.36×10⁻²⁴ cm³ × 6.02214×10²³ mol⁻¹ = 9.25 cm³/mol ✔(Kittel: 9 cm³ mol⁻¹)
- 對照檢查(用錯誤的 1.12 會怎樣):R = 1.12 × 2.56 = 2.867 Å → a = 4.055 → a³/4 = 16.67 ų → 10.04 cm³/mol。與 9 差 12%,可見 Kittel 用的確實是 1.09 而不是 2^(1/6)。
- 反向檢查(這個模型對別的惰性氣體準不準):同樣流程套 Ar,σ = 3.40 Å → R₀ = 3.706 Å,而 Table 4 印的實測最近鄰距離是 3.76 Å,只差 1.4%——書 p.59 稱之為 “The agreement with (14) is remarkable.”。對 Ne,σ = 2.74 → R₀ = 2.987 Å,實測 3.13 Å,差 4.6%,而且是實測偏大。這個「越輕偏差越大、而且一律往外偏」的模式,就是 p.59 說的 “The slight departure of R₀/σ for the lighter atoms from the universal value 1.09 predicted for inert gases can be explained by zero-point quantum effects.” He 只是把這個趨勢推到極端而已。
補出腳註沒寫的算術之二:27.5 與 36.8 cm³/mol 反推出來的膨脹幅度,正好對上「30 到 40 percent」補推導
腳註把兩件事並排寫出來卻沒有把它們連起來:一邊說「漲落是最近鄰距離的 30 到 40 percent」,一邊說「略去零點運動算出 9 cm³/mol,實測是 27.5 與 36.8」。這兩句話其實是同一件事的兩種說法,而書上沒有點破。連起來的方式只是開一次三次方根。
物理圖像:零點運動讓原子在格點附近大幅晃動。晃動的原子等效上佔據了更大的體積——如果硬把它們擠回古典的平衡距離,動能(∝ 1/R²,見下一個公式區塊)會急遽升高。系統的因應方式是把晶格撐開,用位能的一點損失換動能的大幅下降。所以「漲落幅度佔鍵長的比例」與「實際鍵長相對於古典預測的膨脹比例」應該是同一個量級。
結論:算出來是 +44%(He⁴)與 +59%(He³),腳註說漲落是 30–40%。同一個量級,而且順序正確(He³ 比 He⁴ 更誇張)。對一個純粹靠量綱串起來的論證,這已經是很強的一致性。差距(44% vs 40%)也合理:膨脹量還包含了位能與動能競爭後的最佳化結果,不會恰好等於漲落幅度。
順帶解決一個常被問到的問題:為什麼 He³ 的莫耳體積比 He⁴ 大?兩者的電子結構完全相同(同樣是兩顆 1s 電子),所以 ε 與 σ 一模一樣,位能曲線一模一樣。唯一的差別是原子核的質量:3.016 u 對 4.003 u。輕的那個零點動能大,所以撐得更開。這是一個把「同位素效應」孤立到只剩質量一個變數的乾淨案例,也是書 p.60 那個 Ne²⁰ / Ne²² 晶格常數差異(4.4644 Å 對 4.4559 Å)的極端放大版。
- 古典基準:V_m(古典) = 9.25 cm³/mol(上一個區塊算出;Kittel 寫 9)。對應 R₀ = 2.790 Å。
- He⁴ 的體積膨脹比:27.5 / 9.25 = 2.97 倍
- 換成線膨脹比:2.97^(1/3) = 1.437,即最近鄰距離從 2.790 Å 撐到 4.01 Å,膨脹 +44%
- He³ 的體積膨脹比:36.8 / 9.25 = 3.98 倍 → 線膨脹 3.98^(1/3) = 1.585 → R = 4.42 Å,膨脹 +59%
- 與腳註的 30–40% 對照:漲落 30–40%、膨脹 44%(He⁴),同一量級 ✔。若改用 Kittel 原文的 9(而非 9.25)當基準,He⁴ 的線膨脹為 (27.5/9)^(1/3) = 1.450,即 +45%,結論不變。
- 質量比檢查:He³/He⁴ 的線膨脹比 = 1.585/1.437 = 1.103。若零點能 ∝ 1/M 而膨脹量 ∝ (零點能)^(某冪次),質量比 (4.003/3.016) = 1.327 開三次方是 1.099——與 1.103 只差 0.4%。這暗示膨脹量大致隨 M^(−1/3) 變化,與「零點能 ∝ 1/M」的設定自洽。(此為量綱層級的巧合式檢驗,不是嚴格推導。)
- 注意一個誠實的但書:27.5 與 36.8 是液態的莫耳體積(腳註原文寫 liquid He⁴ 與 liquid He³),而 9 是假想固態的值。拿液體跟固體比嚴格說不對等,但零壓下的氦根本沒有固態可比——這正是腳註要說的重點,所以這個比較是唯一可能的比較。
零點能的最簡估計式(Kittel 在 p.60 給的模型,此處補成公式)
E_zp = p²/2M = (h/λ)²/2M ∝ 1/M (λ 由局限原子的邊界決定)
p.60E_zp = 零點動能;M = 原子質量;λ = 量子波長,由局限空間的邊界決定;p = h/λ 是德布羅意關係。三維盒子邊長 L 的基態則為 E_zp = 3h²/(8ML²),對應 λ = 2L。
把腳註的伏筆一次收回:用零點能的標度律解釋 p.59 的 28%/10%/6%/4%,並外推到氦補推導
腳註說「原子越重,零點效應越不重要」,但沒有給任何量化。書 p.59 則在講完 Lennard-Jones 內聚能之後丟出一串數字:“Quantum-mechanical corrections act to reduce the binding by 28, 10, 6, and 4 percent of Eq. (16) for Ne, Ar, Kr, and Xe, respectively.”——同樣沒有解釋這串數字從何而來。
把兩處連起來,就能做一件書上沒做的事:用標度律預測這串數字,再把同一條標度律外推到氦,看看會發生什麼。
兩個候選標度律:
(甲)Kittel p.60 明講的質量標度:p.60 對硬牆盒子模型只下了一句結論——“the quantum zero-point correction to the energy is inversely proportional to the mass”,也就是修正 ∝ 1/M。(完整寫出來其實是 E_zp/(井深) ∝ 1/(M·L²·ε),取 L = R₀ = 1.09σ 就是 1/(Mσ²ε);但下面刻意只保留 Kittel 明講的 1/M,把井深與尺寸的差異略去。這是有意的簡化,因為這裡要裁決的問題只有一個:質量依賴究竟是 1/M 還是 1/√M。附帶一提,若真的照完整的 1/(Mσ²ε) 歸一化到 Ne = 28%,Ar 會掉到 2.75%、Kr 0.84%、Xe 0.32%,連遞減趨勢的量級都保不住,反而測不出質量指數。)
(乙)諧振子(愛因斯坦)模型:零點能 = (3/2)ħω,而 ω = √(k/M),彈性常數 k 的量綱是 ε/σ²,所以 ħω ∝ (ħ/σ)√(ε/M)。修正 = ħω/(井深 ∝ ε) ∝ ħ/(σ√(Mε))。這個無量綱組合在低溫物理裡有專門的名字,但 Kittel 本章沒有引入它,所以這裡只當作一個量綱組合使用。
兩者對質量的依賴不同(1/M 對 1/√M),正好可以用 p.59 的四個數字來裁決。下面把兩者都算出來,歸一化到 Ne = 28%。
裁決結果:(乙)勝出。Ar/Kr/Xe 的預測值 8.8/4.9/3.0 對上實際的 10/6/4,全部落在 12–25% 的誤差內,而且趨勢完全正確;(甲)的 14.1/6.7/4.3 則在 Ar 處高出 41%。這說明零點能確實該用晶格振動的標度(∝ 1/√M)而不是 Kittel 那句話字面上的 1/M 來估。
然後是收網的一步:把(乙)外推到 He。得到的修正超過 100%——零點能大於鍵能。這在物理上只有一個意思:Lennard-Jones 井根本關不住氦原子,零壓下不存在固態氦。這正是腳註第二句 “They do not solidify at zero pressure even at absolute zero temperature.” 的定量版本。而且 He³(146%)比 He⁴(127%)更誇張,與它更大的液態莫耳體積(36.8 對 27.5 cm³/mol)方向一致。
Lennard-Jones 參數全部取自書 p.53 Table 4(有文字層):ε(10⁻¹⁶ erg)He 14、Ne 50、Ar 167、Kr 225、Xe 320;σ(Å)He 2.56、Ne 2.74、Ar 3.40、Kr 3.65、Xe 3.98。原子量取自標準值。
- 計算標度(乙)的分母 σ√(Mε),M 用原子量(u)、ε 用 10⁻¹⁶ erg、σ 用 Å(比值只看相對大小,單位可以任意固定):
- Ne:2.74 × √(20.18 × 50) = 2.74 × √1009 = 2.74 × 31.76 = 87.0
- Ar:3.40 × √(39.95 × 167) = 3.40 × √6672 = 3.40 × 81.68 = 277.7
- Kr:3.65 × √(83.80 × 225) = 3.65 × √18855 = 3.65 × 137.3 = 501.2
- Xe:3.98 × √(131.3 × 320) = 3.98 × √42016 = 3.98 × 204.9 = 815.6
- 歸一化到 Ne = 28%(修正 ∝ 1/分母):Ar = 28 × 87.0/277.7 = 8.8%(實際 10);Kr = 28 × 87.0/501.2 = 4.9%(實際 6);Xe = 28 × 87.0/815.6 = 3.0%(實際 4)
- 對照標度(甲)——只取 Kittel 明講的 1/M,歸一化到 Ne = 28%:Ar = 28 × 20.18/39.95 = 14.1%(實際 10,高估 41%);Kr = 28 × 20.18/83.80 = 6.7%(實際 6);Xe = 28 × 20.18/131.3 = 4.3%(實際 4)。Kr 與 Xe 兩者都行,Ar 這一點把 1/M 淘汰掉。
- 外推到 He⁴:分母 = 2.56 × √(4.003 × 14) = 2.56 × √56.04 = 2.56 × 7.486 = 19.16 → 修正 = 28 × 87.0/19.16 = 127%
- 外推到 He³:分母 = 2.56 × √(3.016 × 14) = 2.56 × √42.22 = 2.56 × 6.498 = 16.63 → 修正 = 28 × 87.0/16.63 = 146%
- (甲)的外推也給同樣的定性結論:He⁴ = 28 × 20.18/4.003 = 141%,He³ = 28 × 20.18/3.016 = 187%。兩個模型都跨過 100%,結論不依賴模型的選擇——這是好的量級論證該有的樣子。
- 物理判讀:修正超過 100% 代表把零點能加進去以後,「固體」的總能量高於分離的自由原子,也就是固態不再是基態。氦因此在零壓下一路到 0 K 都保持液態;要讓它固化必須外加約 25 大氣壓以上的壓力,用 PV 功把原子壓進井裡(此壓力值 Kittel 本章未給,屬課外資訊)。
四個常見誤解(每一個都會讓後面的推導整段翻車)我的判讀
誤解一:把 7.9 eV 當成 NaCl 的內聚能。這是本段最該防的錯。7.9 eV 是晶格能(相對於分離的自由 Na⁺ 與 Cl⁻),6.4 eV 才是內聚能(相對於分離的中性 Na 與 Cl)。兩者差 1.5 eV,來源是 I(Na) − A(Cl) = 5.14 − 3.61。判斷方法很簡單:問自己終態是離子還是原子。馬德倫求和算出來的天然是對離子定義的量,所以離子晶體的計算一路做下來得到的是晶格能;要跟 Table 1 比較,必須再走一次游離能/電子親和能的帳。
誤解二:以為凡得瓦力來自永久偶極。Fig. 1 的圖說寫得很清楚是 “fluctuations in the charge distributions”,而 p.49 又特別強調自由原子的電荷分布是球對稱的——球對稱代表永久偶極矩恰好為零。惰性氣體原子之間如果只有靜態的電荷分布,交互作用會是零(p.53 原文明說 the interaction between atoms would be zero)。鍵結完全來自瞬時漲落誘發的偶極—誘發偶極關聯,這是一個純量子的效應,古典靜電學給不出來。
誤解三:把 Table 1 的數字當成「每莫耳分子」。Table 1 的圖說寫的是 “Energy required to form separated neutral atoms…”,單位列印的是 eV/atom,所以 kJ/mol 那一列是每莫耳原子。對單原子固體沒差,但如果你拿它跟離子晶體「每分子單位 6.4 eV」這種數字並排比較,就會差一個因子。NaCl 的 6.4 eV 是每個 NaCl 分子單位(即每兩個原子),折成每原子是 3.2 eV。
誤解四:以為 Table 3 的數字全部量在室溫。標題確實寫 at room temperature,但 He、Ne、C 標 [d] = 4 K,Ar、Kr 標 [a] = 77 K,Hg 標 [c] = 1 K。惰性氣體在室溫是氣體,根本不可能有固體的體模量。方括號的溫度註記與圓括號的晶形註記長得太像,是這張表最容易出錯的地方;引用 Ar 的 0.013 時如果不註明 77 K,讀者會以為那是室溫值,而 77 K 已經是 Ar 熔點 83.81 K 的 92%,熱軟化非常可觀。
這一段接到哪裡去:後續小節、以及與熱容/Cth 研究的接點我的判讀
本區塊全部標 mine:以下是把 p.47–52 的內容接到本書後續章節的說明,頁碼依據為章首目錄頁(p.47)與我實際讀過的那幾頁正文。
(一)直接續接:p.53 起的凡得瓦-倫敦交互作用。本段建立的三個前提會被逐一兌現——「電荷分布球對稱、永久偶極為零」使得靜態交互作用為零,所以必須做漲落;「電子分布幾乎不被扭曲」使得可以用氣相參數;「內聚能只有游離能的 1%」保證微擾展開收斂。模型是兩個相距 R 的一維諧振子,用庫侖項耦合,最後得到 −C/R⁶。
(二)p.56–59:排斥項、平衡晶格常數、內聚能。Lennard-Jones 位能 U(R) = 4ε[(σ/R)¹² − (σ/R)⁶](式 10)、fcc 晶格求和的兩個常數 12.13188 與 14.45392(式 12)、平衡條件 R₀/σ = 1.09(式 14)、內聚能 U_tot(R₀) = −2.15 × 4Nε(式 16)。本段補算 9 cm³/mol 時已經先用了式 (14)。
(三)p.59–60:量子修正。28%/10%/6%/4%,以及 Ne²⁰ 與 Ne²² 的晶格常數差異(4.4644 對 4.4559 Å)。本段的零點運動腳註就是為這裡鋪路的。
(四)p.60 起:離子晶體。本段釘死的 cohesive energy/lattice energy 分界線,在這裡開始真正發揮作用。
(五)p.73–85:彈性常數。Table 3 的體模量在這裡會被寫成 B = ⅓(C₁₁ + 2C₁₂)(目錄上的 Bulk Modulus and Compressibility,p.80)。本段做的量綱估計 B ≈ c·E_coh/V_atom 是它的「零階版本」。
(六)與熱容/熱電容 Cth 研究的接點。這一條值得多說幾句,因為 Table 3 是整本書裡第一次出現「可以直接換算成聲速、進而換算成 Debye 溫度」的資料。
鏈條是這樣的:體模量 B → 聲速 v ≈ √(B/ρ) → Debye 溫度 Θ_D ∝ (ħv/k_B)(6π²n)^(1/3) → 低溫熱容 C ∝ (T/Θ_D)³。也就是說,Table 3 的一個數字,原則上就決定了該材料低溫熱容的整條曲線。用 Si 試一次:B = 0.988×10¹¹ N/m²,ρ = 2.329 g/cm³ = 2329 kg/m³ → v ≈ √(0.988×10¹¹/2329) = √(4.24×10⁷) = 6.5×10³ m/s。這個數字與 Si 的實際體聲速(縱波約 8.4 km/s、橫波約 5.8 km/s)同量級,說明這條粗鏈條是可用的。
對 FinFET 熱容(Cth)萃取而言,這裡有兩個可以立刻拿去用的觀念。第一,熱容的體積歸一化必須落到「每原子體積」——也就是慣用晶胞體積除以該晶胞內的原子數,而不是拿慣用晶胞體積直接除。本段算 B_pred 時反覆用到的 a³/4(fcc,每胞 4 個原子)、a³/8(diamond,每胞 8 個)、a³/2(bcc,每胞 2 個),算的就是這個量。Si 的 a = 5.430 Å 給出每原子 20.01 ų,即原子密度 5.00×10²² cm⁻³;Dulong–Petit 上限 3Nk_B 直接由它決定,是 Cth 的高溫飽和值。第二,內聚能與體模量同向變化這件事,等價於「鍵越強 → 聲速越高 → Debye 溫度越高 → 室溫下熱容越偏離 Dulong–Petit 上限」。金剛石的 B 是全表最大(4.43),對應的 Θ_D 也是全元素最高,所以它在室溫的熱容只有 Dulong–Petit 值的一小部分——而 Si 的 Θ_D 約 645 K,室溫(300 K)落在 T/Θ_D ≈ 0.47,正好在 T³ 律與 Dulong–Petit 平台之間最難近似的那一段。這是 FinFET 熱模型裡 Cth 對溫度敏感的根本原因,也是為什麼不能直接套常數熱容。(Θ_D 的具體數值與 T³ 律屬第 5 章 PHONONS II 的內容,此處只是指出接點。)
校註:本段引用的每一個數字的來源與可信度我的判讀
這一段有一個結構上的特殊性必須交代清楚:書本 p.47–52 這六頁裡,真正有文字層的正文只有 p.49 一頁。p.47 是章首目錄、p.48 是習題清單加 Fig. 1、p.50/51/52 是三張整頁的表格影像。
可以逐字信任的部分(PDF 文字層直接抽取,已經過自訂字符解碼器):p.47 的目錄條目與頁碼、p.48 的習題清單與 Fig. 1 圖說、p.49 的全部正文與腳註、p.50 Table 1 的標題與圖說、p.51 Table 2 的標題(這兩頁只有標題/圖說有文字層,表格本體沒有)、p.53 的 Table 4(該頁表格有文字層)、p.58–62 的正文與 Table 6。
唯一一處「英文原文也只能目視判讀」的例外:書 p.52 Table 3 的標題與整段圖說。該頁連標題都沒有文字層,本段引用的那段英文是算繪後放大逐字抄錄的,可信度低於上面那一類,引用時請自行複核。
屬於目視判讀的部分(PDF 該頁無文字層,以 PyMuPDF 算繪成 3.2–6 倍解析度影像後放大目視):
- 書 p.50 Table 1(內聚能)——全部數值
- 書 p.51 Table 2(熔點)——全部數值
- 書 p.52 Table 3(體模量與壓縮率)——全部數值,連標題與圖說都只能目視
- 書 p.20 Table 3(元素晶體結構與晶格常數)——本段用來算每原子體積的 a 值
目視判讀做了哪些交叉檢驗(都在前面的推導區塊裡):(1)Table 1 的三列單位互相換算,C、W、Ne 三筆完全吻合;(2)Table 3 的 B 與 κ 互為倒數,十組全部通過;(3)由 p.20 判讀的 Ar 晶格常數 5.31 Å 反算最近鄰距離 3.755 Å,與 p.53 Table 4(有文字層)印的 3.76 Å 吻合。三組獨立的檢驗都通過,判讀可信度高;但只要是上面列出的四張表,本段一律標註「整頁影像目視判讀」,引用時請自行複核。
兩處文字層瑕疵已在對應位置標註校註:p.49 腳註的 “temp- erature” 換行連字號、p.61 Fig. 8 圖說中 (7.9 − 5.1 + 3.6) = 6.4 這一段浮出到段末。