這一節在解決什麼問題我的判讀
第 3 章開頭把問題定得很乾脆：是什麼把晶體黏在一起？書上 p.49 的答案是「電子的負電荷與原子核的正電荷之間的吸引性靜電交互作用完全負責固體的內聚」，磁力只有微弱影響，重力可以忽略。既然如此，剩下的工作就是分門別類：不同的凝態形式之所以差別這麼大，照書上的說法，最終是因為最外層電子與離子實（ion cores）的分布不同。
前一節處理惰性氣體晶體：那裡的原子是中性的、閉殼層的，吸引力只剩下微弱的感應偶極—偶極（van der Waals）作用，隨距離以 1/R⁶ 衰減，內聚能只有每原子零點幾個 eV。這一節把電荷打開：讓原子先交換一個電子變成 ±q 的離子，吸引力立刻升級成 1/r 的庫侖作用。1/r 比 1/R⁶ 長程得多，也強得多，結果就是鹼金屬鹵化物的晶格能比惰性氣體晶體的內聚能大上數十到數百倍（NaCl 每分子單元 7.9 eV；Table 4 的惰性氣體是每原子 Ne 0.02 eV、Xe 0.17 eV，所以對 Ne 約四百倍、對 Xe 只有約 46 倍——注意兩者的計量單位一個是每分子單元、一個是每原子）。
但「1/r 很強」這句話單獨拿出來沒有用，因為晶體裡不只有一對離子：參考離子周圍有異號的最近鄰（吸引），也有同號的次近鄰（排斥），第三層又是異號……整個晶體的靜電能是一個正負交錯的無窮級數。這個級數的和能不能收斂、收斂到多少，就是 Madelung 常數 α 要回答的問題。這一節的全部技術內容可以壓成一句話：
把 1/r 的兩體庫侖作用在整個晶格上求和，得到 −Nαq²/R；再加上一個射程極短的排斥項把晶體撐開，由平衡條件定出 R₀，最後得到晶格能 −(Nαq²/R₀)(1 − ρ/R₀)。
後面的每一個式子都只是在把這句話說清楚。讀這一節時建議一直把兩個問題放在手邊：（a）α 為什麼會收斂、收斂到 1.75 這種「略大於 1」的數？（b）排斥項為什麼只削掉 10% 而不是 50%？這兩個問題的答案分別對應第 (21)(25) 式與第 (24) 式。
起點：離子都是閉殼層，電荷分布近似球對稱
書上開宗明義：「Ionic crystals are made up of positive and negative ions. The ionic bond results from the electrostatic interaction of oppositely charged ions.」兩種最常見的離子晶體結構——氯化鈉結構與氯化銫結構——在第 1 章已經出現過。
接著是這一整節得以成立的關鍵物理前提，書上 p.60 逐字寫的是：「The electronic configurations of all ions of a simple ionic crystal correspond to closed electronic shells, as in the inert gas atoms.」書上拿氟化鋰當例子，按照書前扉頁的週期表：
| 物種 | 電子組態 | 等電子於 |
|---|---|---|
| Li（中性） | 1s² 2s | — |
| F（中性） | 1s² 2s² 2p⁵ | — |
| Li⁺ | 1s² | He |
| F⁻ | 1s² 2s² 2p⁶ | Ne |
這件事為什麼重要？因為惰性氣體原子有閉殼層，而閉殼層的電荷分布是球對稱的。書上由此推論：「We expect that the charge distributions on each ion in an ionic crystal will have approximately spherical symmetry, with some distortion near the region of contact with neighboring atoms.」
球對稱這個性質是後面所有計算的隱形支柱（以下的 Gauss 定律論證為本站補充）。原因是：對於球對稱的電荷分布，在分布之外的任何一點，它產生的電位與「把全部電荷集中在球心的點電荷」完全一樣（這是 Gauss 定律的標準結論）。所以把離子當成點電荷 ±q 放在格點上，不是粗糙的近似，而是在殼層不重疊的前提下嚴格成立的。這也解釋了為什麼排斥項只需要考慮最近鄰：只有最近鄰才真正發生殼層重疊，一旦重疊，球對稱的點電荷等效性就失效，需要另外一個經驗項來描述。
注意書上用了「approximately」與「with some distortion」。真實晶體裡離子並非完美球形，接觸區有變形；這個變形正是後面第 16 章談介電常數與極化率、以及殘餘的 van der Waals 貢獻（1–2%）的來源。（原書 p.56 談電子極化率時把讀者指向「Chapter 15」，但第 8 版的 Dielectric Constant and Polarizability 一節實際在第 16 章。）
Fig. 7：X 光電子密度圖是「球對稱」的實驗證據
書上說這幅球對稱的圖像「is confirmed by x-ray studies of electron distributions (Fig. 7)」。網頁不放圖，以下用文字把 Fig. 7 描述清楚（等高線的數值來自文字層，圖形拓樸為算繪後目視判讀）。
圖說逐字是：「Electron density distribution in the base plane of NaCl, after x-ray studies by G. Schoknecht. The numbers on the contours give the relative electron concentration.」
畫面內容：這是 NaCl 底面（base plane）上的一張等高線圖，也就是切過一層離子的正方形截面。四個角落與中心交錯排列著 Na⁺ 與 Cl⁻——左上 Na⁺、右上 Cl⁻、左下 Cl⁻、右下 Na⁺，形成棋盤式的交替。每個離子周圍是一圈一圈幾乎同心的圓形等高線，由內而外標著 100、50、10、2、1、0.5、0.1 等相對電子濃度。圖上還有兩條水平與兩條垂直的細線，分別通過上下兩對、左右兩對離子的中心，構成一個網格，標示晶軸方向。
三個要從圖上讀出來的訊息：
第一，等高線幾乎是圓的。這就是「球對稱」的直接證據。如果離子之間有明顯的共價鍵（電子雲沿鍵方向拉長、在鍵中央堆積），等高線會在連線方向上明顯變形成啞鈴狀——圖上沒有。
第二，離子之間有一片很低的谷。兩個離子連線的中點附近，等高線降到 0.1 以下；四個離子圍成的正方形中央有一條封閉等高線標著 0.05（標籤印在這條等高線的右下側，不在正方形的幾何中心），它圍出的中心區域是全圖密度最低處。也就是說，離子之間的電子密度掉了三個數量級。這正是「離子晶體＝一堆彼此幾乎不共享電子的帶電球」這個圖像的實驗依據，也是把交互作用寫成兩體庫侖項 ±q²/r 的正當性來源。
第三，Cl⁻ 的等高線比 Na⁺ 伸得遠。在圖上 Cl⁻ 的 0.1 等高線半徑明顯大於 Na⁺ 的。這與離子半徑的常識一致：陰離子得到一個電子、核電荷不足以束縛，半徑膨脹；陽離子失去一個電子、核電荷相對過剩，半徑收縮。後面 Table 7 的最近鄰距離 R₀（例如 LiF 2.014 Å vs. RbI 3.671 Å）就是這兩個半徑相加的結果。
一個容易誤讀的地方：等高線上的數字是相對電子濃度，不是絕對的 electrons/ų，所以不要拿它去積分算電子數。它只能用來比較「哪裡多、哪裡少」。
數量級論證：先確認方向沒走錯
在動手做任何級數求和之前，書上先做了一個 30 秒就能算完的估計，目的是建立信心。書上 p.60 逐字寫的是：「A quick estimate suggests that we are not misguided in looking to electrostatic interactions for a large part of the binding energy of an ionic crystal.」
估計的內容是：NaCl 晶體中，一個正離子到最近的負離子的距離是 2.81 × 10⁻⁸ cm，「the attractive coulomb part of the potential energy of the two ions by themselves is 5.1 eV」。注意「by themselves」這三個字——這是把其他所有離子都拿掉、只剩孤零零一對離子時的庫侖吸引能。
這個 5.1 eV 要拿去跟什麼比？書上接著說：「This value may be compared (Fig. 8) with the experimental value of 7.9 eV per molecular unit for the lattice energy of crystalline NaCl with respect to separated Na⁺ and Cl⁻ ions.」
這個對照是這一節最漂亮的一步棋（以下兩點的解讀為本站補充）。一對離子就已經給了 5.1 eV，整個晶體給的是 7.9 eV——比值 7.9/5.1 ≈ 1.55。這說明兩件事：
- 靜電作用確實撐起了絕大部分的結合能。如果庫侖項只佔一小部分，5.1 和 7.9 之間會差一個數量級，而不是差 55%。既然量級對得上，「往靜電交互作用去找答案」這個方向就是對的——這正是書上「we are not misguided」的意思。
- 但整個晶體的結果不是「一對」的簡單放大。7.9 不等於 5.1 × 6（六個最近鄰），也不等於 5.1。這中間的差距，恰恰就是後面 Madelung 常數 α 與排斥修正因子 (1 − ρ/R₀) 要負責解釋的。先記住這個 1.55，等到 (24) 式算完，你會看到它被拆成 α × (1 − ρ/R₀) = 1.7476 × 0.886 ≈ 1.548。
書上這句話之後只寫了一句：「We now calculate the energy more closely.」然後就進入 Electrostatic or Madelung Energy 小節。
兩個離子之間的長程靜電交互作用
U_coul(r) = ± q² / r （CGS） U_coul(r) = ± q² / (4πε₀ r) （SI）
p.60q：離子電荷的大小（對鹼金屬鹵化物是一個基本電荷 e）；r：兩離子中心的距離；上方的 + 號用於同號離子（排斥），− 號用於異號離子（吸引）；ε₀：真空電容率。
cohesive energy 與 lattice energy：第 3 章一開頭埋的分界線
這裡要把第 3 章第一段埋下的定義串回來。書上 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.」
「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.」
差別只有一個詞：neutral free atoms 還是 free ions。但這個差別在離子晶體裡會造成好幾個 eV 的落差，因為從「中性原子」走到「離子」要付游離能、又會拿回電子親和力。（以下的對照表與延伸討論為本站補充。）
| 名稱 | 終態 | NaCl 的實驗值 |
|---|---|---|
| lattice energy（晶格能） | 分離的 Na⁺ 與 Cl⁻ 自由離子 | 7.9 eV / 分子單元 |
| cohesive energy（內聚能） | 分離的中性 Na 與 Cl 自由原子 | 6.4 eV / 分子單元 |
注意在惰性氣體那一節，這兩個詞不需要區分（原子本來就是中性的），所以第一次讀的時候很容易把它當成同義詞滑過去。到了離子晶體，這條線才真正咬人。書上在 p.63 又重申了一次晶格能的定義：「The total lattice energy is defined as the energy required to separate the crystal into individual ions at an infinite distance apart.」——請注意這裡說的是 individual ions。
換句話說：這一整節從 (17) 到 (24) 算出來的 U_tot，是晶格能，不是內聚能。(24) 式的 −Nαq²/R₀(1 − ρ/R₀) 的參考態是分離的自由離子。如果你要的是相對於中性原子的內聚能，還得再做一次 Fig. 8 的帳。這是考試與論文裡最常見的數量級錯誤之一。
Fig. 8 的能量收支：6.4 eV 是怎麼來的補推導
書上只把答案 (7.9 − 5.1 + 3.6) = 6.4 印在圖說裡，中間每一項的正負號為什麼那樣配，得自己想清楚。以上是補出來的逐步帳。
這筆帳的結構其實就是Born–Haber 循環的骨架（書上這一節沒有用這個名字），它在實驗上非常重要：晶格能無法直接量測（你沒辦法真的把晶體一顆一顆拆成氣相離子），但游離能、電子親和力、內聚能（可由昇華熱等熱化學數據得到）都量得到。於是實務上是反過來用這條式子：
晶格能 = 內聚能 + 游離能 − 電子親和力 = 6.4 + 5.14 − 3.61 = 7.9 eV
書上 Table 7 的 Experimental 欄位（NaCl 182.6 kcal/mol）就是這樣得到的——所謂「實驗晶格能」其實是由一串熱化學量拼出來的，不是一次量出來的。這也是為什麼書上特別加一句「All values on the figure are experimental」：Fig. 8 三個數字都是量得到的，只有它們拼起來的邏輯是我們加的。
最後注意一個術語上的小陷阱：Fig. 8 第三列的 7.9 eV 在圖上被標成「Cohesive energy」，但圖說又說「The lattice energy with respect to separated ions is 7.9 eV」。以 p.49 的嚴格定義來看，7.9 eV 是晶格能；圖上那個標籤是寬鬆用法。真正的內聚能是圖說算式給的 6.4 eV。
步驟 0：先把 Fig. 8 用文字畫出來。這張圖由三列「化學反應式」組成，每一列左邊是初態、右邊是終態，箭頭上下標著能量的名稱。內容如下（算繪後目視確認）：
第一列：
Na（Gas） + 5.14 eV → Na⁺（Gas） + e，5.14 eV 下方標 Ionization energy。第二列：
e + Cl（Gas） → Cl⁻（Gas） + 3.61 eV，3.61 eV 下方標 Electron affinity。第三列：
Na⁺（Gas） + Cl⁻（Gas） → Na⁺Cl⁻（Crystal） + 7.9 eV，7.9 eV 下方標 Cohesive energy。圖說逐字是：「The energy per molecule unit of a crystal of sodium chloride is (7.9 − 5.1 + 3.6) = 6.4 eV lower than the energy of separated neutral atoms. The lattice energy with respect to separated ions is 7.9 eV per molecule unit. All values on the figure are experimental.」並註明游離能見 Table 5、電子親和力見 Table 6。
步驟 1：規定符號。把「系統釋出能量」記為負的能量變化 ΔE < 0（終態比初態低），「系統吸收能量」記為 ΔE > 0。我們要算的是
ΔE_total = E(晶體) − E(分離的中性 Na 原子 + 分離的中性 Cl 原子)
一路走三步：中性原子 → 離子 → 晶體。
步驟 2：游離 Na。這一步要花錢。書上圖中寫的是「Na + 5.14 eV → Na⁺ + e」，也就是必須注入 5.14 eV 才能把電子拔掉。
ΔE₁ = +5.14 eV
步驟 3：Cl 抓電子。這一步賺回來。「e + Cl → Cl⁻ + 3.61 eV」，系統放出 3.61 eV。Table 6 的說明逐字是「The electron affinity is positive for a stable negative ion.」——電子親和力為正就代表 Cl⁻ 是穩定的。
ΔE₂ = −3.61 eV
步驟 4：離子凝聚成晶體。這一步賺最多。「Na⁺ + Cl⁻ → 晶體 + 7.9 eV」。
ΔE₃ = −7.9 eV
這個 7.9 eV 正是晶格能——它的初態剛好就是「分離的自由離子」，完全符合 p.49 對 lattice energy 的定義。
步驟 5：加總。
ΔE_total = ΔE₁ + ΔE₂ + ΔE₃ = (+5.14) + (−3.61) + (−7.9) = −6.37 eV
書上圖說寫成 (7.9 − 5.1 + 3.6) = 6.4 eV，是把「系統降低了多少」寫成正數，並且把 5.14 捨成 5.1、3.61 捨成 3.6。兩者是同一筆帳，只差一個整體符號與四捨五入。
ΔE_total = −6.4 eV，所以 NaCl 晶體的內聚能是 6.4 eV / 分子單元。
警告：兩個「5.1 eV」不是同一件事我的判讀
這是本節最容易踩的地雷，而且書上完全沒有提醒，因為它是一個純粹的數字巧合。
在 p.60，書上說「一對孤立 Na⁺Cl⁻ 離子的庫侖吸引位能是 5.1 eV」。
在 p.61 的 Fig. 8 圖說，算式 (7.9 − 5.1 + 3.6) = 6.4 裡面也有一個 5.1。
這兩個 5.1 eV 毫無關係。
- p.60 的 5.1 eV = q²/R = 14.4 / 2.81，是靜電位能，量綱來自庫侖定律。
- 圖說裡的 5.1 eV = Na 的第一游離能 5.14 eV 四捨五入，是原子物理的量，跟晶體結構、跟 R 一點關係都沒有。
它們碰巧都接近 5.1，純屬偶然。Na 的游離能 5.14 eV 是 Fig. 8 給的（原始資料在 Table 5，該表在本站讀取的頁面範圍之外，故不引用其他元素的數值）；而 e²/R₀ = 14.3996/2.81 = 5.12 eV 是由 Table 7 的幾何算出來的。兩個數字用的輸入完全沒有交集，卻在小數第一位撞在一起。
怎麼防呆？看到 5.1 eV 出現時，問一句「它的初態和終態是什麼」：
- 庫侖那個 5.1：初態是兩個相距無窮遠的離子，終態是相距 2.81 Å 的一對離子。
- 游離能那個 5.1：初態是中性 Na 原子，終態是 Na⁺ 加一個無窮遠的電子。
完全不同的兩個過程。寫報告或答題時把這兩個數字混用，是會被直接扣分的那種錯。
參考離子 i 的總交互作用能
U_i = Σ′_j U_ij
(17) p.62U_i：把離子 i 與晶體中所有其他離子的交互作用能全部加起來的結果；U_ij：離子 i 與離子 j 之間的交互作用能；求和號上的一撇 ′ 表示排除 j = i 這一項（書上逐字：「where the summation includes all ions except j = i」），也就是不算自己跟自己的交互作用（那會發散）。
兩體交互作用能：Born–Mayer 形式
U_ij = λ exp(−r_ij / ρ) ± q² / r_ij （CGS） U_ij = λ exp(−r_ij / ρ) ± q² / (4πε₀ r_ij) （SI）
(18) p.62λ：排斥作用的強度參數，量綱是能量；ρ：排斥作用的射程參數，量綱是長度；r_ij：離子 i 與 j 的距離。書上逐字：「where the + sign is taken for the like charges and the – sign for unlike charges.」——庫侖項的符號：同號取 +（排斥），異號取 −（吸引）。注意 λ exp(−r/ρ) 這一項永遠是正的（排斥），不帶 ±。
為什麼離子晶體改用 λ exp(−r/ρ)，而不沿用惰性氣體的 R⁻¹²
惰性氣體那一節用的是 Lennard-Jones 位能，排斥項是 R⁻¹²。到了離子晶體，書上換成指數形式，並且給了兩個理由。原文（p.62）逐字如下：
「We treat the strength λ and range ρ as constants to be determined from observed values of the lattice constant and compressibility; we have used the exponential form of the empirical repulsive potential rather than the R⁻¹² form used for the inert gases. The change is made because it may give a better representation of the repulsive interaction. For the ions, we do not have gas-phase data available to permit the independent determination of λ and ρ.」
拆開來看是三件事：
（一）換形式的理由：「it may give a better representation of the repulsive interaction」。書上的用詞是 may，相當保守——這是一個經驗上的選擇，不是定理。注意書上並沒有在這裡給出微觀的理由（那屬於本站補充，見下一段）。
（二）參數怎麼定：由晶格常數與壓縮率反推。λ 和 ρ 是兩個未知數，所以需要兩個實驗輸入。書上選的是觀測到的晶格常數（給出 R₀）與壓縮率（給出曲率，即體模量 B）。Table 7 就是這樣建起來的：它同時列出 R₀、B、zλ、ρ 四欄，前兩欄是輸入，後兩欄是輸出。
（三）為什麼不能像惰性氣體那樣獨立定參數：「we do not have gas-phase data available」。這一句是本段最有洞察力的地方。惰性氣體的 ε 與 σ 可以從氣相的第二維里係數、黏滯係數等數據獨立定出來，然後拿去預測晶體性質——那是一個真正的預測。但離子沒有這條路：你在氣相裡沒有一罐自由的 Na⁺ 和 Cl⁻ 讓你量它們的兩體位能曲線。所以 λ、ρ 只能從晶體本身的數據回推。
這一點決定了 Table 7 的「Calculated」欄該怎麼解讀。因為 R₀ 已經被用來定參數了，模型當然會重現 R₀——那不是預測。真正有預測力的是晶格能：R₀ 與 B 是輸入，晶格能是輸出。Table 7 最後一欄的 Calculated 就是某個以 R₀、B 為輸入的模型輸出（書上沒有交代 Tosi 用的是哪一個模型，本站不臆測），而書上 p.65 說計算值與觀測值「in exceedingly good agreement」。這才是模型的成績單。
（順帶一提排斥項的射程假設：書上在 (19) 式只在最近鄰之間放排斥項，其他所有距離只留庫侖項。這是另一個獨立的簡化，理由是 ρ ≈ 0.1R₀ 使得次近鄰的排斥能小到可以丟掉——見後面的討論。）
ρ 的物理意義，以及指數形式背後的微觀理由我的判讀
書上對 ρ 只給了一句精確的定義性描述（p.62）：「We note that ρ is a measure of the range of the repulsive interaction; when r = ρ, the repulsive interaction is reduced to e⁻¹ of the value at r = 0.」
把這句話寫成式子就很清楚：排斥項在 r = 0 的值是 λ e⁰ = λ；在 r = ρ 的值是 λ e⁻¹ = λ/2.718 = 0.368λ。所以 ρ 就是排斥位能衰減 1/e 所需的距離——它是一個「衰減長度」，和 RC 電路的時間常數、或指數衰減的半衰期是同一類的量。從 Table 7 看，ρ 落在 0.29–0.37 Å 之間，也就是不到半個埃。
以下是書上沒有給、但值得知道的微觀理由（本站補充）：
為什麼是指數，不是冪次？排斥的來源是兩個閉殼層的電子雲重疊（Pauli 原理不許同自旋電子佔同一軌域，被迫升到更高能階；書上 p.85 的章末摘要對排斥作用的一般性說法是「the electrostatic repulsion of overlapping charge distributions and the Pauli principle, which compels overlapping electrons of parallel spin to enter orbitals of higher energy」）。而束縛態的電子波函數在原子核外以指數衰減：ψ ~ exp(−κr)，其中 ℏ²κ²/2m 對應束縛能。兩個原子的波函數重疊積分因此 ~ exp(−κd)，而排斥能大致正比於重疊的某個冪次，仍然是指數。所以 exp(−r/ρ) 在函數形式上有物理根據，而 R⁻¹² 純粹是為了計算方便挑的（12 = 2 × 6，平方項好算）。ρ 的倒數 1/ρ 大致就是價電子波函數的衰減常數 2κ 的量級：ρ ≈ 0.3 Å 對應 1/ρ ≈ 3.3 Å⁻¹，與閉殼層外緣的波函數衰減率同量級，數量級上說得通。
這對後面章節有什麼影響？先破除一個常見的誤解：指數項每微分一次帶出一個 1/ρ ≈ 3.1 Å⁻¹（NaCl），而 R⁻¹² 每微分一次帶出 12/R₀ ≈ 4.3 Å⁻¹——兩者其實是同一個量級，不能說指數排斥的導數「比冪次項大得多」。真正的差別在於位能井的形狀：排斥項的陡峭程度決定了位能井底附近的二次項係數（→ 體模量、聲速）與三次項係數（→ 非諧性 → 熱膨脹與熱阻）。第 5 章談熱膨脹與聲子—聲子散射時，用的就是井底的三次項。
為什麼是 U_tot = N U_i 而不是 2N U_i補推導
這一步的算術很簡單，但概念上是初學者最常出錯的地方之一，所以值得把三個數字的角色講清楚。先補一個對照：同一件事在惰性氣體那一節也出現過，只是長相不同。書上 p.58 的 (11) 式寫的是 U_tot = ½N(4ε)[…]，並在 p.59 解釋「The factor ½ occurs with the N to compensate for counting twice each pair of atoms.」那裡 N 是原子數，所以要明寫 ½；這裡 N 是分子單元數（離子數的一半），½ 已經被吸收進 N 裡，所以看不到 ½。兩處是同一個防重複計數的機制。
三個數字的角色：
- 2N 是離子的個數。
- N 是分子單元（formula unit）的個數，也是「鍵」或「離子對」的有效個數。
- U_i 是單一個離子與其他全部離子的交互作用能。
錯誤寫成 U_tot = 2N U_i 的人，其實是把每一對交互作用算了兩次，會得到剛好兩倍大的晶格能。檢查方法很簡單：量綱檢查不會抓到這個錯，只有數值對照能抓到。NaCl 的實驗晶格能是 7.9 eV / 分子單元；如果你算出 15.8 eV，回頭找這個因子 2。
另一個相關的記帳陷阱：「每分子單元」與「每離子」差兩倍。Table 7 的晶格能欄位單位是 kcal/mol，這個 mol 指的是分子單元的莫耳數（NaCl 一莫耳 = 6.022 × 10²³ 個 Na⁺Cl⁻ 對），不是離子的莫耳數。把 182.6 kcal/mol 換成 eV 時，用 1 eV = 23.061 kcal/mol，得 7.92 eV——是每分子單元，與 Fig. 8 的 7.9 eV 相符。
最後注意書上在 (17) 的 U_i 裡用的是 Σ′（排除 j = i），而 U_tot = N U_i 這一步用的是「U_i 對所有離子相同」——這兩件事都依賴忽略表面效應。對一個 1 cm³ 的晶體，表面離子佔比約 10⁻⁷，完全可以忽略；但對奈米顆粒（比如 5 nm 的晶粒），表面離子可能佔 20% 以上，這時 U_tot = N U_i 就不再準確。這是奈米材料的表面能為什麼重要的起點。
- 設定。一個含 N 個「分子單元」的離子晶體共有 2N 個離子（N 個陽離子 + N 個陰離子）。書上 p.63 逐字寫的是「the total lattice energy U_tot of a crystal composed of N molecules or 2N ions as U_tot = N U_i. Here N, rather than 2N, occurs because we must count each pair of interactions only once or each bond only once.」書上只給了這一句話，以下把它拆成算式。
步驟 1：把總能量寫成所有「無序配對」之和。系統的真正總能量是對所有相異的離子對 {i, j} 各算一次：
U_tot = Σ_{pairs {i,j}} U_ij
這裡 {i, j} 與 {j, i} 是同一對，只出現一次。
步驟 2：改寫成「對每個離子求和」。如果改成先固定 i、對所有 j ≠ i 求和，再對所有 i 求和：
Σ_i Σ′_j U_ij = Σ_i U_i
這個雙重求和把每一對 {i, j} 數了兩次——一次是 i 當參考離子、j 當鄰居，另一次是 j 當參考離子、i 當鄰居。
步驟 3：除以 2。
U_tot = ½ Σ_i U_i
步驟 4：用「U_i 與參考離子種類無關」這個性質。前面 (17) 的討論已經說過，在 NaCl 結構裡，無論參考離子是正離子還是負離子，U_i 的值都一樣（忽略表面效應）。於是 Σ_i U_i = (離子總數) × U_i = 2N U_i。
U_tot = ½ × 2N × U_i = N U_i
把交互作用按「是不是最近鄰」分兩種情形
U_ij = { λ exp(−R/ρ) − q²/R （最近鄰） { ±(1/p_ij)(q²/R) （其他）
(19) p.63R：晶體中的最近鄰距離（nearest-neighbor separation）；p_ij：無量綱的距離比，定義為 r_ij ≡ p_ij R，所以最近鄰的 p_ij = 1；λ、ρ 同 (18)。第二行的 ± 依照同號／異號取號。
從 (19) 走到 (20)：把求和拆成兩堆補推導
這一段是全節最容易在符號上翻船的地方，所以把它拆得比書上細。書上從 (19) 直接跳到 (20)，中間「合併兩堆」與「符號約定」兩步完全沒寫。
要記住的結構：(20) 只有兩項，而且兩項的角色完全不同——
- zλ e^(−R/ρ)：正的、短程的、只跟最近鄰有關。裡面的 z 是局域的配位數。
- −αq²/R：負的、長程的、跟整個晶格有關。裡面的 α 是把全部無窮多個離子加起來的結果。
換句話說，排斥項認得的是「鄰居」，庫侖項認得的是「整個晶體」。這個不對稱是離子鍵與共價鍵最大的差別，也是為什麼離子晶體的性質（晶格能、體模量）對結構型態（NaCl vs. CsCl vs. ZnS）敏感，而共價晶體幾乎只看配位數。
步驟 1：把 (17) 的求和依照 (19) 的兩種情形拆開。
U_i = Σ′_j U_ij = Σ_{最近鄰 j} [λ exp(−R/ρ) − q²/R] + Σ_{其餘 j} [±(1/p_ij)(q²/R)]
步驟 2：處理第一堆（最近鄰）。設 z 為任一離子的最近鄰個數（書上逐字：「where z is the number of nearest neighbors of any ion」）。這 z 項每一項都完全相同（都在距離 R 處、都是異號），所以
Σ_{最近鄰} [λ exp(−R/ρ) − q²/R] = z λ exp(−R/ρ) − z (q²/R)
步驟 3：處理第二堆（其餘）。把共同因子 q²/R 提出來：
Σ_{其餘} ±(1/p_ij)(q²/R) = (q²/R) Σ_{其餘} (±)/p_ij
步驟 4：合併兩堆的庫侖部分。注意步驟 2 留下的 −z(q²/R) 可以寫成 (q²/R) × (−z)，而最近鄰的 p_ij = 1、符號為負（異號離子），所以 −z 正好就是「最近鄰對 Σ′(±)/p_ij 的貢獻」。於是兩堆的庫侖項可以合成一個完整的求和：
U_i = z λ exp(−R/ρ) + (q²/R) Σ′_j (±)/p_ij
書上對這一點有一句話提醒（p.63）：「The sum should include the nearest-neighbor contribution, which is just z.」——意思就是這裡的 Σ′ 必須把最近鄰那 z 項算進去，別因為前面已經把它們挑出來處理排斥項就忘了它們的庫侖部分。
步驟 5：把 ± 的約定一次講死，再定義 Madelung 常數。這是整段唯一容易翻船的地方，所以先釐清符號。(18) 裡的 ±q²/r 是「同號取 +、異號取 −」（那是位能本身的物理符號）；但 (21) 裡的 (±) 是另一套約定——書上 p.64 明講：「If we take the reference ion as a negative charge, the plus sign will apply to positive ions and the minus sign to negative ions.」也就是異號取 +、同號取 −，剛好相反。
為什麼要翻過來？因為書上要讓 α 是正的：「For (20) to give a stable crystal it is necessary that α be positive.」把負號提到 α 前面統一寫成 −αq²/R，α 本身就變成正數，讀起來直觀。於是定義
α ≡ Σ′_j (±)/p_ij （異號取 +、同號取 −）
這就是 (21) 式。最近鄰全是異號、p_ij = 1，貢獻 +z（NaCl 是 +6）；次近鄰同號、p = √2，貢獻 −12/√2 = −8.49；第三層異號、p = √3，貢獻 +8/√3 = +4.62……正負交錯抵消，最後收斂到 +1.7476。
步驟 5′：把步驟 3、4 的求和用 α 改寫。步驟 4 得到的庫侖總和是「以 (18) 的物理符號」寫的，也就是異號帶負、同號帶正，正好是 α 的相反數：
Σ′_j (±)_{(18)}/p_ij = −α
所以
U_i = z λ exp(−R/ρ) + (q²/R)(−α) = z λ exp(−R/ρ) − α q²/R
負號代表淨吸引，符合物理。
步驟 6：得到 (20)。
U_i = z λ exp(−R/ρ) − α q²/R
再用前面證過的 U_tot = N U_i：
U_tot = N ( zλ e^(−R/ρ) − αq²/R )
總晶格能（平衡前，R 仍為變數）
U_tot = N U_i = N ( z λ e^(−R/ρ) − α q² / R ) （CGS）
(20) p.63N：分子單元（離子對）個數，晶體共有 2N 個離子；z：任一離子的最近鄰個數（NaCl 結構 z = 6，CsCl 結構 z = 8，立方 ZnS 結構 z = 4）；λ、ρ：Born–Mayer 排斥參數；α：Madelung 常數；R：最近鄰距離，此處仍是自由變數。SI 要把 q² 換成 q²/4πε₀。
Madelung 常數的定義
α ≡ Σ′_j (±) / p_ij ≡ Madelung constant
(21) p.63p_ij：由 r_ij ≡ p_ij R 定義的無量綱距離比（最近鄰 p_ij = 1）；(±)：以參考離子為準，異號離子取 +、同號離子取 −；Σ′：對所有 j ≠ i 求和。α 是純數，沒有量綱。
平衡條件
N dU_i/dR = − (N z λ / ρ) exp(−R/ρ) + N α q² / R² = 0 （CGS）
(22) p.63各符號同 (20)。這是 dU_tot/dR = 0 展開後的結果。
決定平衡距離 R₀ 的超越方程
R₀² exp(−R₀/ρ) = ρ α q² / (z λ) （CGS）
(23) p.63R₀：平衡最近鄰距離；其餘同前。SI 要把 q² 換成 q²/4πε₀。
平衡時的總晶格能
U_tot = − (N α q² / R₀) ( 1 − ρ / R₀ ) （CGS）
(24) p.63U_tot：2N 個離子在平衡距離 R₀ 時的總晶格能（負值代表比分離的自由離子低）；−Nαq²/R₀ 這一項書上明確命名為 Madelung 能量（「The term −Nαq²/R₀ is the Madelung energy.」）；(1 − ρ/R₀)：排斥項造成的修正因子。
(20) → (22) → (23) → (24) 的完整代數補推導
書上把 (20)(22)(23)(24) 四條式子排在一起，只用「may be written, using (20) and (23)」一句話帶過中間的代數。真正卡住讀者的是步驟 5 的技巧：看到 (23) 是超越方程，直覺會想「那要怎麼解 R₀？」——但其實根本不用解。把平衡條件當成一條代數恆等式用，直接消去指數項，答案就出來了。
這個技巧在固態物理裡反覆出現，值得抽象出來記住：
若 u(R) = A f(R) + B g(R)，平衡條件 A f′ + B g′ = 0 給出 A f = −B (g′/f′) f。把它代回 u，就能把 u(R₀) 只用一種能量表示。
在 Lennard-Jones 的情形（惰性氣體那一節），同樣的操作給出 U(R₀) = −(排斥冪次與吸引冪次的比值決定的因子) × 吸引能；在這裡指數排斥給出的因子是 (1 − ρ/R₀)。章末習題 5（Linear ionic crystal）要你對 A/Rⁿ 型排斥重做一遍，答案是 U(R₀) = −(2Nq² ln 2/R₀)(1 − 1/n)——注意那裡的因子是 (1 − 1/n)，對 n = 9 是 0.889，與指數模型的 0.886 幾乎一樣，這解釋了為什麼兩種排斥形式給出的晶格能差別很小。
還有一個檢查值得做：(24) 對 R₀ 的依賴幾乎是純 1/R₀（因為 ρ 的變化比 R₀ 小）。所以同一個結構型態的一系列晶體，晶格能應該大致隨 1/R₀ 變化。用 Table 7 檢查：LiF（R₀ = 2.014 Å，242.3 kcal/mol）與 RbI（R₀ = 3.671 Å，144.9 kcal/mol）。1/R₀ 的比是 3.671/2.014 = 1.823，晶格能的比是 242.3/144.9 = 1.672。兩者相差 9%，來源正是修正因子不同：LiF 的 ρ/R₀ = 0.291/2.014 = 0.144，RbI 的 ρ/R₀ = 0.348/3.671 = 0.0948，(1−0.144)/(1−0.0948) = 0.945，把 1.823 × 0.945 = 1.723，與 1.672 差 3%。剩下的差距來自 van der Waals、零點能與實驗誤差。
步驟 1：寫下要微分的對象。每分子單元的能量
u(R) ≡ U_tot/N = z λ e^(−R/ρ) − α q² / R
步驟 2：微分。
du/dR = zλ · d/dR[e^(−R/ρ)] − αq² · d/dR[R⁻¹]
第一項：d/dR e^(−R/ρ) = e^(−R/ρ) · (−1/ρ) = −(1/ρ)e^(−R/ρ)。
第二項：d/dR R⁻¹ = −R⁻², 所以 −αq² × (−R⁻²) = +αq²/R²。
du/dR = −(zλ/ρ) e^(−R/ρ) + αq²/R²
乘上 N 就是書上的 (22)。
步驟 3：令 R = R₀ 使其為零，得平衡條件。
(zλ/ρ) e^(−R₀/ρ) = α q² / R₀² ……（★）
步驟 4：整理成 (23)。（★）兩邊同乘 ρ R₀² / (zλ)：
R₀² e^(−R₀/ρ) = ρ α q² / (z λ) ✔ 這就是 (23)。
步驟 5：這一步是關鍵——不要去解 (23)，而是用（★）把指數項換掉。由（★）直接得
z λ e^(−R₀/ρ) = ρ α q² / R₀²
（把（★）兩邊同乘 ρ 即可。）注意左邊正是 (20) 裡的排斥項在 R = R₀ 的值。這是全段最漂亮的一招：我們不需要知道 R₀ 到底是多少，只要知道它滿足平衡條件，就能把排斥項用庫侖項表示出來。
步驟 6：代回 (20)。
u(R₀) = zλ e^(−R₀/ρ) − αq²/R₀
= ρ α q² / R₀² − α q² / R₀
提出共同因子 −αq²/R₀：
u(R₀) = − (α q² / R₀) [ 1 − ρ/R₀ ]
步驟 7：乘上 N。
U_tot = − (N α q² / R₀)( 1 − ρ / R₀ ) ✔ 這就是 (24)。
(1 − ρ/R₀) 到底削掉多少：用 Table 7 逐一檢查補推導
書上對這件事只寫了一句「ρ is of the order of 0.1R₀」。上面把它用 Table 7 的數據完整驗證，並且抽出趨勢。
這個結論的重要性在於它讓「純靜電圖像」變得可用。如果修正因子是 0.5，你就不能說「離子晶體的結合基本上是靜電的」；但既然它是 0.89，你可以說。這正是書上 p.61 那句斷言的量化基礎：「The main contribution to the binding energy of ionic crystals is electrostatic and is called the Madelung energy.」
順便把三個小數字放在一起對照，它們常被搞混：
| 效應 | 對內聚能的相對大小 | 出處 |
|---|---|---|
| Madelung（靜電）能量 | 100%（基準） | (24) 的 −Nαq²/R₀ |
| 排斥項的修正 | −9.5% 到 −14% | (24) 的 (1 − ρ/R₀)，本站由 Table 7 算出 |
| van der Waals | +1% 到 +2% | 書上 p.61 明講 |
也就是說：排斥比 van der Waals 重要一個數量級，這是為什麼 (18) 一定要有排斥項、卻可以完全不寫 van der Waals 項。書上把後者明講成「a relatively small contribution to the cohesive energy in ionic crystals, of the order of 1 or 2 percent」。
- 步驟 1：書上的宣稱。p.63 逐字：「We shall find that ρ is of the order of 0.1R₀, so that the repulsive interaction has a very short range.」這句話說的是一個數量級的宣稱，可以直接用 Table 7 查證。
步驟 2：從 Table 7 取 R₀ 與 ρ，算出比值 ρ/R₀ 與修正因子 (1 − ρ/R₀)。
晶體 R₀ (Å) ρ (Å) ρ/R₀ 1 − ρ/R₀ LiF 2.014 0.291 0.144 0.856 LiCl 2.570 0.330 0.128 0.872 LiI 3.000 0.366 0.122 0.878 NaF 2.317 0.290 0.125 0.875 NaCl 2.820 0.321 0.114 0.886 NaI 3.237 0.345 0.107 0.893 KCl 3.147 0.326 0.104 0.896 KI 3.533 0.348 0.0985 0.902 RbCl 3.291 0.323 0.0981 0.902 RbI 3.671 0.348 0.0948 0.905 - 步驟 3：讀出結論。16 個鹼金屬鹵化物的 ρ/R₀ 全部落在 0.0948（RbI）至 0.144（LiF） 之間，中位數約 0.11。書上說「of the order of 0.1R₀」完全成立。修正因子 (1 − ρ/R₀) 落在 0.856 至 0.905，也就是排斥項削掉靜電吸引的 9.5% 到 14%。
- 步驟 4：注意一個系統性趨勢。ρ 本身變化不大（0.29–0.37 Å，只變了 25%），但 R₀ 從 2.01 變到 3.67 Å（變了 82%）。所以 ρ/R₀ 主要由 R₀ 決定：晶體越大，排斥的相對重要性越低。最極端的是 LiF（最小的晶體，0.144）與 RbI（最大的，0.0948）。
- 步驟 5：為什麼 ρ 幾乎是常數？ρ 是閉殼層電子雲的指數衰減長度。所有這些離子都是 He、Ne、Ar、Kr、Xe 型的閉殼層，最外層電子的束縛能都在 20–40 eV 等級（不是價電子的幾個 eV），對應的波函數衰減長度自然都在 0.3 Å 附近。ρ 是離子的性質，R₀ 是晶體的性質——這是為什麼兩者變化幅度差這麼多。（此步為本站補充的物理詮釋。）
Madelung 常數：歷史、定義與正負號約定
Evaluation of the Madelung Constant 這個小節一開頭是三句歷史與方法的交代（p.64，逐字）：
「The first calculation of the coulomb energy constant α was made by Madelung. A powerful general method for lattice sum calculations was developed by Ewald and is developed in Appendix B. Computers are now used for the calculations.」
三件事：（1）α 由 Madelung 首先算出，常數因此得名。（2）Ewald 發展了一套通用的晶格求和方法，書上放在附錄 B。（3）現在都用電腦算。
接著書上重述 (21) 的定義，然後交代了關鍵的符號約定：
「For (20) to give a stable crystal it is necessary that α be positive. If we take the reference ion as a negative charge, the plus sign will apply to positive ions and the minus sign to negative ions.」
這一句要讀慢一點。它其實是兩個獨立的陳述：
第一，α 必須為正。這不是自然而然的，而是一個穩定性要求。回看 (20)，庫侖項寫成 −αq²/R；只有 α > 0 才代表吸引、才可能有束縛態。如果某個假想結構算出 α < 0，那個結構在靜電上就是不穩定的，根本不會形成。
第二，符號怎麼配。以負離子為參考時，正離子（異號、吸引）取 +，負離子（同號、排斥）取 −。反過來以正離子為參考時，結論一樣（負離子取 +），因為 (21) 描述的是「異號取 +」。
接著書上給了一個等價但更好用的定義，就是 (25) 式。
Madelung 常數的等價定義（用真實距離而非距離比）
α / R = Σ′_j (±) / r_j
(25) p.64r_j：第 j 個離子到參考離子的真實距離（有長度量綱）；R：最近鄰距離；(±)：異號取 +、同號取 −。左邊 α/R 有 1/長度 的量綱，兩邊量綱一致。
一維交錯離子鏈：α = 2 ln 2 的每一步補推導
書上這個一維例子寫得很完整（步驟 1 到 5 都在書上），本站補的是步驟 6 與表格化的整理。但即使書上寫得完整，有兩個地方值得停下來想：
（一）為什麼 α = 1.386 而不是 2？如果只算兩個最近鄰，α 會是 2。但次近鄰是同號的、要減掉 2/2 = 1；第三層加回 2/3；第四層減 2/4……一路抵消到 1.386。這個「比配位數小很多」的現象在三維同樣發生：NaCl 的 z = 6，α 卻只有 1.7476。Madelung 常數不是配位數，它是配位數經過無窮多層正負抵消之後的殘值。
（二）為什麼一維算得出來，三維算不出來？因為一維的幾何是平凡的：第 n 層永遠恰好 2 個離子、距離恰好 nR、符號恰好交替。三維沒有這種規律——距離 √2R 有 12 個、√3R 有 8 個、2R 有 6 個……而且同一層裡正負離子可能混合（在某些結構裡），符號不再是簡單交替。書上下一段開頭就是這個意思：「In three dimensions the series presents greater difficulty. It is not possible to write down the successive terms by a casual inspection.」
一維的結果 2 ln 2 在章末習題 5 會再用到：那題要你對一條含 2N 個交錯離子、排斥項為 A/Rⁿ 的鏈，證明 U(R₀) = −(2Nq² ln 2/R₀)(1 − 1/n)。注意那裡的 2 ln 2 正是這裡算出的 α，而 (1 − 1/n) 是冪次排斥版的修正因子（對應本節指數排斥的 (1 − ρ/R₀)）。
設定（Fig. 9）。書上 Fig. 9 畫的是一條無窮長的直線，線上交錯排列著正負離子，相鄰兩離子的距離都是 R，圖上標出一個「Reference ion」。書上逐字：「As an example, we compute the Madelung constant for the infinite line of ions of alternating sign in Fig. 9. Pick a negative ion as reference ion, and let R denote the distance between adjacent ions.」
把參考離子（負）放在原點，往右數：+ 在 R、− 在 2R、+ 在 3R、− 在 4R……往左完全對稱。
步驟 1：列出每一項的距離與符號。用 (25) 的形式 α/R = Σ′(±)/r_j。
第幾層 距離 r_j 該層離子種類 相對參考離子（負） 符號 每層個數 1 R 正 異號（吸引） + 2 2 2R 負 同號（排斥） − 2 3 3R 正 異號 + 2 4 4R 負 同號 − 2 n nR 交替 交替 (−1)^(n+1) 2 「每層 2 個」是這題唯一的幾何輸入，書上逐字解釋：「the factor 2 occurs because there are two ions, one to the right and one to the left, at equal distances r_j.」
步驟 2：把求和寫出來。
α/R = 2[ 1/R − 1/(2R) + 1/(3R) − 1/(4R) + … ]
這就是書上印出來的那一行。
步驟 3：把 1/R 提出來，得到純數。
α = 2[ 1 − 1/2 + 1/3 − 1/4 + … ]
這也是書上印出來的。右邊那個級數叫交錯調和級數（alternating harmonic series）。
步驟 4：用書上給的展開式求和。書上逐字：「We sum this series by the expansion
ln(1 + x) = x − x²/2 + x³/3 − x⁴/4 + …」
對照可知：只要令 x = 1，右邊就變成 1 − 1/2 + 1/3 − 1/4 + …，正是我們要的級數。於是
1 − 1/2 + 1/3 − 1/4 + … = ln(1 + 1) = ln 2
步驟 5：得到答案。
α = 2 ln 2 = 1.386294…
書上結語：「Thus the Madelung constant for the one-dimensional chain is α = 2 ln 2.」
步驟 6（書上沒寫，但數學上必要的一步）：x = 1 恰好在收斂半徑的邊界上。ln(1+x) 的 Taylor 級數收斂半徑是 |x| < 1，x = 1 是端點。在端點上級數仍然收斂（交錯級數判別法：項的絕對值單調遞減趨零），而且收斂到 ln 2——這需要 Abel 定理（若冪級數在端點收斂，則其和等於函數在該點的極限）才嚴格成立。書上直接代 x = 1 沒有出錯，但這一步不是「顯然」的。
這件事和三維的困難是同一個根源：交錯調和級數是條件收斂（把絕對值加起來就是發散的調和級數），條件收斂級數的和依賴求和順序（Riemann 重排定理保證可以重排成任何值）。書上在下一頁對三維提出的警語，講的正是這件事。
三維的警語：為什麼不能「隨手逐項寫下來」
這是整節最該原樣傳達、也最常被跳過的一段。書上 p.65 開頭逐字：
「In three dimensions the series presents greater difficulty. It is not possible to write down the successive terms by a casual inspection. More important, the series will not converge unless the successive terms in the series are arranged so that the contributions from the positive and negative terms nearly cancel.」
兩個層次的困難，書上用「More important」明確分了輕重：
困難一（較輕）：寫不出通項。「It is not possible to write down the successive terms by a casual inspection.」一維可以一眼看出第 n 層是 2 個離子、距離 nR、符號 (−1)^(n+1)。三維的 NaCl 結構，你得先枚舉所有整數座標 (l, m, n)，算出距離 √(l²+m²+n²) × (a/2)，再判斷 l+m+n 的奇偶決定正負——這已經不是「casual inspection」能做的。
困難二（更重要）：級數根本不會收斂，除非你把正負項排成幾乎互相抵消。這句話的措辭要小心讀：書上說的不是「收斂得很慢」，而是 will not converge。
這裡的物理圖像是（以下為本站補充的詮釋）：如果你笨拙地一層一層往外加（例如先把所有正離子加完再加負離子，或按某種讓正負不平衡的順序），每加一層就相當於在晶體外面裹上一層淨電荷，而一個帶淨電荷的球殼在中心產生的電位不隨半徑衰減——於是級數不收斂。只有當你的求和方式讓每一個「批次」都是電中性的（或幾乎電中性的），外層的貢獻才會快速衰減。
這正是 Evjen 法（把每個離子按它被幾個立方體共用而拆成分數電荷，使每個立方殼層恰好電中性；Evjen 法本身並未出現在原書，此處為本站補充）與 Ewald 法（把求和拆成實空間的快收斂部分與倒空間的快收斂部分）背後的想法。書上把 Ewald 法放在附錄 B，並說「Computers are now used for the calculations.」
實務上的意義（本站補充）：如果你自己寫程式算 Madelung 常數，直接用 for l,m,n in range(-L,L) 加一遍 (-1)^(l+m+n)/sqrt(l²+m²+n²)，你會得到一個隨 L 震盪、不收斂的數列——這不是程式有 bug，而是書上警告的這件事。要嘛用 Ewald 法，要嘛至少用球對稱的中性化（例如對截斷邊界的離子加權）。
為什麼條件收斂會讓答案依賴求和順序我的判讀
書上把「the series will not converge unless…」講清楚了，但沒有說為什麼。以下是本站補充的解釋，分數學與物理兩層。
數學層：條件收斂與 Riemann 重排定理。
一個級數 Σ aₙ 若 Σ|aₙ| 發散但 Σ aₙ 收斂，稱為條件收斂。Riemann 重排定理說：條件收斂級數可以藉由重排項的順序，收斂到任何你指定的實數，甚至發散。Madelung 級數正是這種級數——把所有項取絕對值就是 Σ 1/r_j，而三維空間中距離 r 到 r+dr 的殼層裡有 ∝ r² dr 個離子，所以 Σ|1/r_j| ~ ∫ r²(1/r) dr = ∫ r dr 發散得很厲害。所以「Madelung 級數的和」這個說法嚴格來講是不完整的，必須同時指定求和方式。
物理層：為什麼不同求和方式對應不同的物理情境。
更深一層的理解是：不同的求和順序對應不同的表面電荷分布。想像你把晶體截斷成一個有限大的塊：
- 如果截斷面是電中性的（陰陽離子一樣多），外部沒有淨電荷，中心的電位收斂到「正確」的值。
- 如果截斷面帶了淨電荷（例如你把一個 NaCl 立方體切在全是 Na⁺ 的面上），表面就有一層電荷，它在中心產生一個與晶體大小無關的額外電位。
所以「求和順序」在物理上就是「你怎麼選晶體的表面」。這不是數學家的吹毛求疵——真實的極性表面（polar surface，例如 NaCl 的 (111) 面）確實有這個問題，稱為極性災難（polar catastrophe），真實晶體會透過表面重構或吸附來中和它。這是表面科學與氧化物介面（例如著名的 LaAlO₃/SrTiO₃ 二維電子氣）的核心議題之一。
實作上的建議（本站補充）：
- 最簡單可靠的方法是 Ewald 求和（書上附錄 B）：把 1/r 拆成 erfc(ηr)/r（實空間快收斂）與 erf(ηr)/r（其傅立葉變換在倒空間快收斂），兩邊各自絕對收斂，再減掉自作用項。
- 快速驗證程式對不對，就拿 NaCl 結構跑一遍，看能不能得到 1.747565（書上給到小數第六位，正好可以當作六位精度的測試基準）。
- 一維的 2 ln 2 = 1.386294 也是很好的單元測試。
三個典型 Madelung 常數，與跨書比對時最大的踩雷點
書上 p.65 給了一張小表，前面逐字寫的是：「Typical values of the Madelung constant are listed below, based on unit charges and referred to the nearest-neighbor distance:」
| Structure | α |
|---|---|
| Sodium chloride, NaCl | 1.747565 |
| Cesium chloride, CsCl | 1.762675 |
| Zinc blende, cubic ZnS | 1.6381 |
兩個前提條件寫在表前那句話裡，一個都不能漏：
（一）based on unit charges。這些 α 是對 ±1 價離子算的。對 ±2 價的晶體（例如 MgO、BaO，也是 NaCl 結構），α 的數值不變（因為 α 只是幾何量），但 (20)(24) 裡的 q² 要換成 (2e)² = 4e²，晶格能一口氣變成四倍。章末習題 7（Divalent ionic crystals）就是在練這件事：比較假想的 Ba⁺O⁻ 與 Ba⁺⁺O⁻⁻。
（二）referred to the nearest-neighbor distance。這就是書上 p.64 明講的那個陷阱：
「The value given for α will depend on whether it is defined in terms of the nearest-neighbor distance R or in terms of the lattice parameter a or in terms of some other relevant length.」
這是跨書比對 Madelung 常數時最常見的錯誤來源。有些教科書（尤其是化學系的結晶學教材）習慣用晶格參數 a 當基準，那麼同一個 NaCl 結構的 α 會印成 3.49513 而不是 1.747565。兩個數字都對，只是定義不同。
怎麼判斷你看到的 α 是哪一種？看它是否落在 1.6–1.8 附近。以最近鄰距離為基準的 α，對常見的三維結構都在這個窄帶裡（下一段解釋為什麼）；如果你看到 2.5、3.5 這類數字，幾乎可以確定是以 a 為基準的。
為什麼三個結構的 α 差這麼少？這是一個值得想的問題（以下為本站補充）。NaCl 的 z = 6、CsCl 的 z = 8、ZnS 的 z = 4，配位數差了兩倍，α 卻只在 1.64 到 1.76 之間（差 7%）。原因是：配位數越大的結構，最近鄰越多（有利），但次近鄰的同號離子也越近（不利），兩者幾乎抵消。CsCl 結構裡，次近鄰（6 個同號 Cs）的距離只有 2R/√3 = 1.155R，非常近；而 NaCl 結構的次近鄰在 √2 R = 1.414R。所以 CsCl 雖然有 8 個最近鄰，淨得益只比 NaCl 多 0.9%。
這個「只差 0.9%」有很重要的後果：既然靜電上 NaCl 與 CsCl 結構幾乎一樣好，決定實際採用哪個結構的就不是 Madelung 能量，而是排斥項——也就是離子半徑比。這呼應了書上 p.61 那句話：「The ions arrange themselves in whatever crystal structure gives the strongest attractive interaction compatible with the repulsive interaction at short distances between ion cores.」注意「compatible with」三個字——結構是吸引與排斥妥協的結果，不是吸引單獨決定的。
換基準長度時 α 怎麼變：NaCl 與 CsCl 的實算補推導
書上只寫了一句警語（「The value given for α will depend on whether it is defined in terms of the nearest-neighbor distance R or in terms of the lattice parameter a or in terms of some other relevant length.」），沒有給換算的例子。上面把三個結構都算了一遍。
實務上怎麼防呆？三個檢查，做完再引用任何一本書的 α：
- 看數值範圍。1.6–1.8 → 以 R 為基準；2–4 → 大概是以 a 為基準。
- 看有沒有「per ion pair」還是「per ion」。有些文獻把 α 定義成包含因子 2（對應 U_tot = −(αq²/2R) × 2N ions 的寫法），數值會差一倍。
- 自己算一次晶格能對回實驗值。這是最終極的檢查：用你手上的 α 與該書的公式算 NaCl，看能不能得到 ≈ 7.9 eV / 分子單元（= 182 kcal/mol = 762 kJ/mol）。算出 15.8 或 3.95，就是差了因子 2；算出別的，就是基準長度弄錯。
一個延伸提醒（本站補充）：這個「無量綱常數依賴你選什麼尺度」的問題不是 Madelung 常數獨有的。第 6 章的 Fermi 能量、第 5 章的 Debye 溫度、第 2 章的結構因子歸一化，都有類似的慣例分歧。養成習慣：引用任何無量綱常數之前，先確認它的分母是什麼。
步驟 1：釐清什麼是不變的。由 (25)，真正的物理量是
S ≡ Σ′_j (±)/r_j
這是「參考離子位置的靜電位（除以 q）」，只依賴晶體結構與尺寸，不依賴你選什麼長度當基準。Madelung 常數是把它無量綱化的結果：選一個長度 L，定義 α_L ≡ L × S。
步驟 2：換基準的規則。
α_L / α_{L′} = L / L′
也就是說，基準長度越大，α 越大，而且成正比。這條規則只需要 (25) 就能得到，一行代數。
步驟 3：NaCl 結構。在 NaCl 結構裡，立方晶格參數 a 與最近鄰距離 R 的關係是 R = a/2（Na⁺ 與 Cl⁻ 沿立方邊交替排列，相鄰兩離子相隔半個邊長）。於是
α_a = α_R × (a/R) = 1.747565 × 2 = 3.49513
這就是某些書上印的 NaCl Madelung 常數。
步驟 4：CsCl 結構。CsCl 結構是簡單立方的 Cl⁻ 網格，Cs⁺ 在體心。最近鄰距離是體對角線的一半：R = (√3/2) a，即 a/R = 2/√3 = 1.154701。於是
α_a = 1.762675 × 1.154701 = 2.03528
步驟 5：立方 ZnS（閃鋅礦）結構。Zn 與 S 各自構成 fcc，相對位移 (¼,¼,¼)a，所以最近鄰距離是體對角線的四分之一：R = (√3/4) a，a/R = 4/√3 = 2.309401。
α_a = 1.6381 × 2.309401 = 3.78303
步驟 6：檢查「α 落在 1.6–1.8」這個判別法。
結構 α（基準 R） R/a α（基準 a） NaCl 1.747565 1/2 3.49513 CsCl 1.762675 √3/2 = 0.8660 2.03528 立方 ZnS 1.6381 √3/4 = 0.4330 3.78303 以 R 為基準時三個值擠在 1.64–1.76 的窄帶；以 a 為基準時散落在 2.04–3.78。所以「看到 1.6–1.8 就是以最近鄰距離為基準」這個判別法是可靠的。
Fig. 10：KCl 的能量分解圖，逐條讀
書上圖說逐字：「Energy per molecule of KCl crystal, showing Madelung (coulomb) and repulsive contributions.」內文引用處寫的是：「The Madelung and repulsive contributions to the binding of a KCl crystal are shown in Fig. 10.」
網頁不放圖，以下把 Fig. 10 完整描述（曲線形狀為算繪後目視判讀，數值標註來自文字層與目視確認）。
座標軸。橫軸標「R, in 10⁻⁸ cm」（也就是 Å），刻度 1 到 6。縱軸標「Energy, in eV」，從 −14 到 +14，每 2 eV 一格。橫軸畫在 y = 0 的位置，所以正負能量分居上下。
曲線一：Repulsive energy，標註 (2.4 × 10⁴) exp(−R/0.30) eV。位於橫軸上方（正能量）。它從左上方陡降：在 R ≈ 2.2 Å 已高達 14 eV 以上（超出圖框），到 R = 3 附近降到約 1 eV，R = 4 之後幾乎貼著橫軸。這條曲線的形狀就是「射程極短」四個字的圖示——它在 R 從 2.2 增加到 4.0（不到兩倍）的過程中掉了三個數量級以上。
曲線二：Coulomb energy，標註 (25.2/R) eV。位於橫軸下方（負能量）。它從左下的約 −12.9 eV（曲線左端點約在 R ≈ 1.95）單調上升，經過 R = 3 附近的約 −8.4 eV，到圖右端 R ≈ 6 時約 −4 eV。校註：圖上的標籤寫「(25.2/R) eV」沒有帶負號，但畫出來的曲線明確是負的；正確的表達應是 −25.2/R eV。這是原書標註的疏漏，不影響圖形本身。
曲線三：Total energy（虛線）。只畫在能看出結構的區間（約 R = 2.4 到 3.6）。它從左邊的高處（約 −3.7 eV，因為排斥項把總能推高）下降，在 R ≈ 3.1 Å 附近觸底，最低點約 −7.3 eV，然後往右上與庫侖曲線合流（因為排斥項已經消失）。圖上有一個標註「Equilibrium position」，用一條引線指向這個最低點。
這裡有一個非常常見的誤讀必須點破：平衡位置是總能量曲線的極小值，不是兩條曲線的交點。事實上這兩條曲線根本不相交——排斥曲線永遠在橫軸上方，庫侖曲線永遠在下方。圖上那個「Equilibrium position」的引線明確指向虛線的谷底。這一點與 (22) 式 dU_tot/dR = 0 完全一致：平衡由導數為零定義，與兩個分項的大小關係無關。
從圖上該讀出的三件事：
- 總能量的極小值很淺、很寬。井底附近曲率不大，這對應 KCl 的體模量只有 1.74 × 10¹¹ dyn/cm²（Table 7）＝ 17.4 GPa，比金剛石之類的共價晶體（約 443 GPa）小一個多數量級（約 25 倍）。
- 左邊的牆很陡、右邊的尾很長。這個不對稱就是非諧性：壓縮比拉伸難得多。第 5 章談熱膨脹時，熱膨脹係數為正正是因為這個不對稱——溫度升高、振幅變大，平均位置往「軟」的那一側（大 R）偏移。
- 在平衡點，排斥能只有庫侖能的約 10%。從圖上讀：R₀ ≈ 3.1 處，庫侖約 −8.1 eV，總能約 −7.3 eV，差額 0.8 eV 就是排斥能，佔 10%。這與 (24) 的 ρ/R₀ = 0.30/3.1 ≈ 0.097 完全吻合。
驗算 Fig. 10：25.2 這個係數從哪來、平衡點在哪補推導
書上只把 Fig. 10 的兩條曲線公式印在圖上，沒有說明它們和 Table 7 的關係。上面把這條線索補完，順便驗證了三件事：Fig. 10 的 25.2 確實等於 αe²、Fig. 10 的平衡點確實是總能極小值、Table 7 的參數確實能重現實驗晶格能。
兩個要留意的地方：
（一）Fig. 10 與 Table 7 的參數不一致。這在教科書裡很常見（圖和表可能來自不同年代的文獻）。做作業或引用時，用 Table 7——它註明了來源是 M. P. Tosi, Solid State Physics 16, 1 (1964)，而 Fig. 10 沒有註明來源。
（二）我用 (24) 算出的 165.3 接近 Table 7 的「Experimental」165.8，而不是「Calculated」161.6。這值得誠實說明：Table 7 的 Calculated 欄與 (24) 式並不是同一個計算——書上沒有交代 Tosi 的 Calculated 欄具體用了什麼模型（可能包含 van der Waals 項、零點能修正等）。本站不臆測其細節。可以確定的是：用 (24) 加上 Table 7 的 R₀ 與 ρ，能在 0.5% 以內重現實驗晶格能，這已經足以支持書上 p.65 的結論「The calculated values of the lattice energy are in exceedingly good agreement with the observed values.」
步驟 1：驗證庫侖項的係數 25.2。KCl 是 NaCl 結構，α = 1.747565。庫侖能（每分子單元）是 −αq²/R；取 q = e，CGS 下 e² = 14.3996 eV·Å。
α e² = 1.747565 × 14.3996 = 25.164 eV·Å
所以庫侖能 = −25.16/R eV（R 以 Å 為單位）。書上圖中標的 25.2/R 就是這個數四捨五入。✔
步驟 2：驗證排斥項的量級。圖上標 2.4 × 10⁴ eV 與 ρ = 0.30 Å。對照 Table 7 的 KCl：zλ = 2.05 × 10⁻⁸ erg，ρ = 0.326 Å。把 zλ 換成 eV：
2.05 × 10⁻⁸ erg ÷ (1.602 × 10⁻¹² erg/eV) = 1.28 × 10⁴ eV
Fig. 10 用的是 (2.4 × 10⁴, 0.30)，Table 7 用的是 (1.28 × 10⁴, 0.326)——兩組參數不同。這不是矛盾：λ 與 ρ 是強相關的擬合參數（ρ 稍微變小，λ 就要變大才能在 R₀ 附近給出同樣的排斥能），Fig. 10 顯然用的是一組較早、較粗的擬合值。
步驟 3：用 Fig. 10 的參數求平衡點。解 du/dR = 0：
(zλ/ρ) e^(−R/ρ) = α e²/R²
(2.4 × 10⁴ / 0.30) e^(−R/0.30) = 25.2/R²
8.0 × 10⁴ e^(−R/0.30) = 25.2/R²
試 R = 3.10：左 = 8.0 × 10⁴ × e^(−10.333) = 8.0 × 10⁴ × 3.25 × 10⁻⁵ = 2.60；右 = 25.2/9.61 = 2.62。幾乎相等。
試 R = 3.09：左 = 2.69；右 = 2.64。已經反號。
R₀ ≈ 3.10 Å。與圖上虛線谷底的位置一致。✔（Table 7 的實驗值是 3.147 Å，差 1.5%。）
步驟 4：算出谷底深度。代 R₀ = 3.10：
庫侖項 = −25.2/3.10 = −8.13 eV
排斥項 = 2.4 × 10⁴ × e^(−10.333) = 2.4 × 10⁴ × 3.25 × 10⁻⁵ = +0.78 eV
總能 = −7.35 eV
與圖上虛線最低點約 −7.3 eV 相符。✔
步驟 5：用 Table 7 的參數重算，對照。R₀ = 3.147 Å，ρ = 0.326 Å，用 (24)：
庫侖項 = −25.164/3.147 = −7.996 eV
修正因子 = 1 − 0.326/3.147 = 0.8964
U/N = −7.996 × 0.8964 = −7.17 eV
換算 kcal/mol：7.17 × 23.061 = 165.3 kcal/mol。Table 7 的 KCl 實驗值是 165.8（0 K 值 169.5），差 0.3%。✔
Table 7：16 個鹼金屬鹵化物，怎麼讀這張表
書上 p.66 的 Table 7 標題是「Properties of alkali halide crystals with the NaCl structure」，說明文字逐字：「All values (except those in square brackets) at room temperature and atmospheric pressure, with no correction for changes in R₀ and U from absolute zero. Values in square brackets at absolute zero temperature and zero pressure, from private communication by L. Brewer.」資料來源註明「Data from various tables by M. P. Tosi, Solid State Physics 16, 1 (1964).」
| 晶體 | R₀ (Å) | B (10¹¹ dyn/cm²) | zλ (10⁻⁸ erg) | ρ (Å) | 晶格能 實驗 | 晶格能 計算 |
|---|---|---|---|---|---|---|
| LiF | 2.014 | 6.71 | 0.296 | 0.291 | 242.3 [246.8] | 242.2 |
| LiCl | 2.570 | 2.98 | 0.490 | 0.330 | 198.9 [201.8] | 192.9 |
| LiBr | 2.751 | 2.38 | 0.591 | 0.340 | 189.8 | 181.0 |
| LiI | 3.000 | (1.71) | 0.599 | 0.366 | 177.7 | 166.1 |
| NaF | 2.317 | 4.65 | 0.641 | 0.290 | 214.4 [217.9] | 215.2 |
| NaCl | 2.820 | 2.40 | 1.05 | 0.321 | 182.6 [185.3] | 178.6 |
| NaBr | 2.989 | 1.99 | 1.33 | 0.328 | 173.6 [174.3] | 169.2 |
| NaI | 3.237 | 1.51 | 1.58 | 0.345 | 163.2 [162.3] | 156.6 |
| KF | 2.674 | 3.05 | 1.31 | 0.298 | 189.8 [194.5] | 189.1 |
| KCl | 3.147 | 1.74 | 2.05 | 0.326 | 165.8 [169.5] | 161.6 |
| KBr | 3.298 | 1.48 | 2.30 | 0.336 | 158.5 [159.3] | 154.5 |
| KI | 3.533 | 1.17 | 2.85 | 0.348 | 149.9 [151.1] | 144.5 |
| RbF | 2.815 | 2.62 | 1.78 | 0.301 | 181.4 | 180.4 |
| RbCl | 3.291 | 1.56 | 3.19 | 0.323 | 159.3 | 155.4 |
| RbBr | 3.445 | 1.30 | 3.03 | 0.338 | 152.6 | 148.3 |
| RbI | 3.671 | 1.06 | 3.99 | 0.348 | 144.9 | 139.6 |
晶格能單位是 kcal/mol（欄位標題逐字：「Lattice energy compared to free ions, in kcal/mol」）；體模量欄位標題是「Bulk modulus B, in 10¹¹ dyn/cm² or 10¹⁰ N/m²」；LiI 的 B 印成 (1.71)，但原書對這組圓括號未置一詞（Table 7 的說明只解釋了方括號＝絕對零度、零壓；圓括號或許表示該值較不確定，這是本站的推測）。
書上對這張表的評語只有一句，但份量很重（p.65）：「The calculated values of the lattice energy are in exceedingly good agreement with the observed values.」
讀表的四條線索（以下四條為本站補充的讀表心得；原書對 Table 7 只留了「in exceedingly good agreement」一句評語）：
（一）欄位的邏輯是「兩進兩出」。R₀ 與 B 是輸入（實驗量），zλ 與 ρ 是由它們擬合出來的參數，最後一欄的 Calculated 晶格能則是某個以 R₀、B 為輸入的模型輸出（書上未交代 Tosi 用的是哪一個模型，本站不臆測——後面用 (24) 重算 NaCl 會落在 Experimental 欄上，見該節說明）。所以前兩欄重現得完美不算本事，最後一欄才是模型的成績。
（二）逐行比較 Experimental 與 Calculated。最好的是 LiF（242.3 vs 242.2，差 0.04%）、NaF（214.4 vs 215.2，差 0.4%）、KF（189.8 vs 189.1，差 0.4%）、RbF（181.4 vs 180.4，差 0.6%）——全是氟化物。最差的是 LiI（177.7 vs 166.1，差 6.5%）與 LiBr（189.8 vs 181.0，差 4.6%）。規律很清楚：陰離子越大、陽離子越小，模型越差。這是可以理解的：大的、易極化的陰離子（I⁻、Br⁻）在小陽離子的強電場下會明顯變形，違背了「球對稱剛性離子」的前提；同時 van der Waals 貢獻也隨極化率上升。
（三）R₀ 與 B 成反比關係。LiF 的 R₀ 最小（2.014）而 B 最大（6.71）；RbI 的 R₀ 最大（3.671）而 B 最小（1.06）。這符合直覺：晶體越緊實越難壓縮。粗略地，B ∝ R₀⁻ⁿ，n 約 3–4。
（四）ρ 幾乎是常數。16 個值全部落在 0.290–0.366 Å，只變動 26%，而 R₀ 變動 82%。前面已經討論過為什麼。
用 Table 7 的參數重算 NaCl 晶格能：兩條路殊途同歸補推導
書上沒有做這個算例，但它是把整節串起來的最好方式：一路從 (20) 算到實驗值，中間沒有任何自由參數可以調。
三個要帶走的換算常數（研究所階段值得背）：
| 常數 | 值 | 用途 |
|---|---|---|
| e²（CGS） | 14.3996 eV·Å | 距離用 Å，庫侖能直接是 eV |
| 1 eV/分子 | 23.061 kcal/mol = 96.485 kJ/mol | 對回化學文獻 |
| 1 erg | 6.242 × 10¹¹ eV | Table 7 的 zλ 欄位 |
一個提醒關於「Calculated」欄。上面算出 182.4–182.6，落在 Table 7 的 Experimental（182.6）上，而不是 Calculated（178.6）。前面已經說過，書上沒有交代 Tosi 的 Calculated 欄用的是什麼模型，本站不臆測。如果你的作業要求「用書上的公式算 NaCl 晶格能並與表比較」，正確的做法是：算出約 182 kcal/mol，並註明它接近 Experimental 欄而非 Calculated 欄，且說明兩者可能來自不同模型——誠實交代差異比硬湊到 178.6 好。
- 目標。驗證 (20) 與 (24) 在 NaCl 上給出一致且正確的答案。輸入全部取自 Table 7 與正文：α = 1.747565，R₀ = 2.820 Å，ρ = 0.321 Å，zλ = 1.05 × 10⁻⁸ erg，z = 6，e² = 14.3996 eV·Å。
步驟 1：單位換算。zλ = 1.05 × 10⁻⁸ erg。1 eV = 1.602 × 10⁻¹² erg，所以
zλ = 1.05 × 10⁻⁸ / 1.602 × 10⁻¹² = 6554 eV
（順帶：單一鍵的 λ = 6554/6 = 1092 eV。）
路線 A：直接用 (20) 在 R = R₀ 計算。
排斥項：zλ e^(−R₀/ρ) = 6554 × exp(−2.820/0.321) = 6554 × exp(−8.785)
exp(−8.785) = 1.530 × 10⁻⁴
排斥項 = 6554 × 1.530 × 10⁻⁴ = +1.003 eV
庫侖項：−α e²/R₀ = −1.747565 × 14.3996 / 2.820 = −25.164/2.820 = −8.923 eV
U/N = 1.003 − 8.923 = −7.920 eV
路線 B：用 (24) 的封閉式。
U/N = −(α e²/R₀)(1 − ρ/R₀) = −8.923 × (1 − 0.321/2.820) = −8.923 × (1 − 0.11383) = −8.923 × 0.88617 = −7.907 eV
步驟 2：比較兩條路線。−7.920（路線 A）vs −7.907（路線 B），相差 0.16%。
為什麼不完全相同？因為 (24) 的推導用了平衡條件 (23)，也就是假設 Table 7 的 (R₀, ρ, zλ) 三者精確滿足 (23)。實際上它們是分別擬合、各自四捨五入的，所以只近似滿足。檢查一下：平衡條件要求排斥項 = (ρ/R₀) × (αe²/R₀) = 0.11383 × 8.923 = 1.0158 eV，而路線 A 直接算出的是 1.003 eV——差 1.3%，正是這個 1.3% 造成上面 0.16% 的差異。這是一個有用的自我檢查：兩條路線算出來差太多，就表示你抄錯了某個參數。
步驟 3：換成 kcal/mol 對回 Table 7。1 eV/分子單元 = 23.061 kcal/mol。
路線 A：7.920 × 23.061 = 182.6 kcal/mol
路線 B：7.907 × 23.061 = 182.4 kcal/mol
Table 7 的 NaCl 實驗值是 182.6 kcal/mol（0 K 值 185.3）。
- 步驟 4：對回 Fig. 8。182.6 kcal/mol ÷ 23.061 = 7.92 eV，與 Fig. 8 圖說的「The lattice energy with respect to separated ions is 7.9 eV per molecule unit」一致。✔ 整條鏈（(17) → (24) → Table 7 → Fig. 8）自洽。
這一節刻意沒算的東西：van der Waals、電子親和力表、與零點能
一個好的模型要知道自己丟掉了什麼。書上在這一節明確點名了三樣被忽略或被外部輸入的東西。
（一）van der Waals：1–2%，被丟掉了。書上 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.」
注意這個 1–2% 的大小剛好與前面算出的 Experimental vs Calculated 差距（NaCl 是 2.2%、LiF 是 0.04%）同量級。所以「模型算出來偏低幾個百分比」有一部分就是缺了這一項。還要注意 van der Waals 是額外的吸引，所以加上它會讓晶格能變大（更負），方向是對的。
（二）電子親和力：Table 6，由實驗輸入。書上 p.62 的 Table 6 標題是「Electron affinities of negative ions」，說明逐字：「The electron affinity is positive for a stable negative ion.」資料來源是 H. Hotop and W. C. Lineberger, J. Phys. Chem. Ref. Data 4, 539 (1975)。完整數值如下：
| 原子 | 電子親和力 (eV) | 原子 | 電子親和力 (eV) |
|---|---|---|---|
| H | 0.7542 | Si | 1.39 |
| Li | 0.62 | P | 0.74 |
| C | 1.27 | S | 2.08 |
| O | 1.46 | Cl | 3.61 |
| F | 3.40 | Br | 3.36 |
| Na | 0.55 | I | 3.06 |
| Al | 0.46 | K | 0.50 |
值得注意的規律：鹵素（Cl 3.61、F 3.40、Br 3.36、I 3.06）遠高於其他元素，這正是它們特別容易形成穩定陰離子、因而特別容易形成離子晶體的原因。另外 Cl 的 3.61 比 F 的 3.40 大——雖然 F 電負度更高。原因是 F⁻ 的 2p 殼層太小，額外電子擠進去要付很大的電子—電子排斥代價。
還有一個容易被忽略的事實：表上全都是「加第一個電子」的親和力，而且都是正的（放熱）。加第二個電子通常是吸熱的。章末習題 7 給了 O 的數據：「the electron affinities of the first and second electrons added to the neutral oxygen atom are 1.5 and −9.0 eV」——第二個電子親和力是 −9.0 eV，也就是 O²⁻ 在自由空間中根本不穩定。它之所以能存在於 MgO、BaO 之類的晶體裡，完全是靠周圍陽離子提供的 Madelung 電位把它穩住。這是 Madelung 能量最深刻的一個後果。
（三）零點能：這一節完全沒提。前一節（惰性氣體）花了不少篇幅談零點動能的量子修正，甚至提到 Ne²⁰ 與 Ne²² 的晶格常數差異。到了離子晶體，書上在這一節沒有討論零點修正。從 Table 7 的方括號可以間接看到它的量級：NaCl 室溫 182.6、0 K 185.3，差 1.5%；這個 2.7 kcal/mol 的差額混合了熱膨脹與零點能兩種效應。離子晶體的離子比惰性氣體原子重、束縛又強得多（書上前一節說「The heavier the atom, the smaller the quantum correction」），所以量子修正在這裡不是主要問題。
這一節最容易出錯的七件事我的判讀
整理成清單，考前掃一遍。全部都是前面已經討論過的內容，這裡集中列出。
1. 把 U_tot 寫成 2N U_i。正確是 N U_i，因為每一對交互作用只能算一次。錯了會得到剛好兩倍的晶格能。→ 檢查：NaCl 應該是 7.9 eV/分子單元，不是 15.8。
2. 把 (24) 算出的東西叫「內聚能」。它是晶格能（參考態是分離的自由離子）。NaCl 的內聚能是 6.4 eV，晶格能是 7.9 eV，差在游離能與電子親和力。
3. 把 p.60 的庫侖 5.1 eV 與 Fig. 8 的游離能 5.1 eV 混為一談。純粹的數字巧合，兩者物理內容完全不同。
4. 引用 α 時不確認基準長度。NaCl 可以是 1.747565（基準 R）或 3.49513（基準 a）。看到 1.6–1.8 才是以最近鄰距離為基準。
5. 以為 Fig. 10 的平衡位置是兩條曲線的交點。平衡位置是總能量曲線的極小值；那兩條曲線一正一負，根本不相交。平衡由 (22) 的 dU/dR = 0 定義。
6. 用 (19) 時忘了 Σ′ 要包含最近鄰的庫侖項。書上特別提醒：「The sum should include the nearest-neighbor contribution, which is just z.」排斥項只算最近鄰，但庫侖項所有離子都要算——包括最近鄰。
7. 寫程式直接暴力加 Madelung 級數。會得到不收斂的震盪結果。這不是 bug，是書上明講的「the series will not converge unless the successive terms in the series are arranged so that the contributions from the positive and negative terms nearly cancel」。要用 Ewald 法（附錄 B）或中性化的求和方式。
額外一條關於單位制：這一節書上全程用 CGS，所以庫侖項是 q²/r 沒有 4πε₀。書上兩次提醒要換 SI（p.62 的「In SI units the coulomb interaction is ±q²/4πε₀r」與 p.63 的「For SI, replace q² by q²/4πε₀」）。寫程式時建議一律用 e² = 14.3996 eV·Å 這個組合常數，它已經把單位制吃進去了，距離給 Å、答案就是 eV，不會出錯。
往後怎麼用：章內習題與後續章節的接點我的判讀
章末摘要怎麼收這一節。書上 p.85 的 SUMMARY 用兩行收束整節：「Ionic crystals are bound by the electrostatic attraction of charged ions of opposite sign. The electrostatic energy of a structure of 2N ions of charge ±q is (CGS) U = −Nαq²/R = −N Σ(±)q²/r_ij, where α is the Madelung constant and R is the distance between nearest neighbors.」注意摘要裡的式子是純 Madelung 能量（沒有 (1 − ρ/R₀) 修正），也就是硬球極限。
第 3 章內部的習題接點：
| 習題 | 練的是什麼 | 與本節的關係 |
|---|---|---|
| 4. Possibility of ionic crystals R⁺R⁻ | 用 Table 5、6 評估假想的 Na⁺Na⁻ 晶體相對於金屬鈉的穩定性（給定 Na 的電子親和力 0.78 eV） | Fig. 8 的能量收支反過來用 |
| 5. Linear ionic crystal | 2N 個交錯離子、排斥 A/Rⁿ，證明 U(R₀) = −(2Nq² ln 2/R₀)(1 − 1/n)，再算壓縮功 | 一維 α = 2 ln 2 加上冪次排斥版的 (24) |
| 6. Cubic ZnS structure | 用 Table 7 的 λ、ρ 與正文的 α，算 KCl 在立方 ZnS 結構下的內聚能，與 NaCl 結構比較 | α 與 z 同時改變的效果 |
| 7. Divalent ionic crystals | BaO（NaCl 結構，R₀ = 2.76 Å），估計 Ba⁺O⁻ 與 Ba⁺⁺O⁻⁻ 相對於中性原子的內聚能 | q → 2e 使庫侖能變四倍；O 的第二電子親和力是 −9.0 eV |
習題 5 與 7 特別值得做完：5 讓你把 (24) 的推導技巧完整複製一次；7 讓你親眼看到「二價離子晶體之所以存在，完全是 Madelung 能量硬把不穩定的 O²⁻ 撐住」。
往後章節的接點：
- 第 3 章後半（Analysis of Elastic Strains 起）：(24) 的位能曲線在 R₀ 附近展開到二次項，就是彈性常數。Table 7 的體模量 B 正是把 λ、ρ 定出來的第二個輸入。
- 第 4 章（聲子）：把 (18) 的兩體位能在 R₀ 附近展開，二次項係數就是力常數 C，決定色散關係 ω(K)。離子晶體特有的縱光學（LO）與橫光學（TO）聲子分裂，根源正是長程庫侖力——這是 Madelung 求和困難在動力學上的對應版本。
- 第 5 章（聲子熱性質）：位能曲線的三次項（Fig. 10 那條總能曲線的左右不對稱）給出非諧性，決定熱膨脹與聲子—聲子散射，進而決定熱導率。
- 第 16 章（介電體與鐵電體）：這一節假設離子是剛性的球；一旦允許離子極化，就進入介電常數與極化率（Clausius–Mossotti 關係）的領域。Fig. 7 提到的「some distortion near the region of contact」就是極化率的胚胎。（LST 關係不在這一章，它在第 14 章 Plasmons, Polaritons, and Polarons。）
- 第 17 章（表面與介面物理）：Madelung 級數的條件收斂直接關聯到極性表面的穩定性問題。
與熱容／熱傳研究的接點（本站補充）。如果你關心的是熱容或熱導（例如元件層級的熱阻抽取），這一節提供的是最底層的輸入：位能曲線的二階導數給聲速與 Debye 溫度（→ 熱容的溫度標度），三階導數給 Grüneisen 參數（→ 熱膨脹與 Umklapp 散射率 → 熱導）。具體地說，(24) 的 (1 − ρ/R₀) 因子之所以只有 0.89，是因為 ρ/R₀ ≈ 0.1；而 ρ/R₀ 小正意味著排斥項極陡、位能井左壁非常硬，這直接造就了離子晶體相對於前一節的惰性氣體晶體明顯較硬的位能井壁，因而 Debye 溫度較高、非諧性較弱。基準要說清楚：換成金剛石、Si 這類共價晶體當對照，鹼金屬鹵化物反而是 Debye 溫度低得多、非諧性強得多的那一方。想從第一原理估計一個介電材料的熱容曲線時，這條鏈（位能參數 → 力常數 → 聲子譜 → C_v(T)）就是起點。