JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

徑向分布函數

對玻璃做繞射,你得到的是寬闊的光暈,而非銳利的斑點——這裡沒有晶格能讓布拉格定律咬住。本篇交給你一件工具,讓「結構」在失去晶格之後依然存活:g(r),那條坐在一顆原子上、數它鄰居的曲線,即使面對一個毫無長程有序的固體,也能給你一個真實的鍵長與一個真實的配位數。

從銳利斑點到寬闊光暈

本階前兩篇已把無序態最核心的事實定了案:一個玻璃或一個非晶質固體沒有長程有序——沒有晶格,因此也沒有倒晶格可供拍照——卻仍保有真實的短程有序,在確定的距離上有確定的最近鄰鍵,甚至在那之外還有中程有序。這帶出一個尖銳的實務問題。在繞射那一階,布拉格定律與倒晶格給了我們一套俐落的機械來讀晶體:銳利的斑點、指標化、由位置定出單位晶胞、由強度定出基元。把同一台繞射儀對準二氧化矽玻璃,斑點卻不見了。取而代之的,是寥寥幾個寬闊、模糊的隆起——散布在一大片角度上的漫射光暈。

為什麼是光暈而不是斑點?銳利的布拉格峰是倒晶格獨有的印記,而只有週期性的晶體才擁有倒晶格;玻璃根本沒有倒晶格,斑點無處可落。但是——這正是關鍵——這片漫射圖案並非毫無特徵。一個真正毫無結構的材料,一團彼此完全無關聯的原子理想氣體,只會散射出平坦、平滑、毫無隆起的背景。光暈之所以存在,恰恰是因為「確實有」短程有序:最近鄰坐在一個定義明確的間距上,而來自這些相關聯配對的波干涉,會把散射強度在某些角度堆高、在另一些角度削薄。光暈,就是短程有序在說話。我們只需要一套統計語言來聽懂它說了什麼——而那套語言,就是徑向分布函數。

坐在一顆原子上:對分布函數 g(r)

想法在此,而且美得非常具體。在玻璃裡隨便挑一顆原子,就坐在它身上,朝外看。只問一個問題:在離我 r 的距離處,其他原子有多擁擠,和材料的平均密度相比如何?這個比值,就是對分布函數 g(r)——也叫成對相關函數。把它讀成一支「擁擠計」。g(r) = 1 表示「這裡就是平均密度,沒什麼特別」;g(r) 大於 1 表示「原子喜歡坐在這個距離,比平均擁擠」;g(r) 小於 1 表示「原子避開這個距離,比平均空曠」。由於每一顆原子都同樣適合當作坐鎮之處,g(r) 其實是一次對所有原子取的平均——它是個統計量,而不是某個凍結排列的快照。

現在往外走,看著 g(r) 一路展開。在極小的 r 處它平貼於零:原子是硬的,無法重疊或互相穿透,所以沒有任何鄰居會靠得比一個鍵長還近。接著 g(r) 猛地躍起,形成一個高而相當銳利的第一峰,峰心落在 r1——也就是最近鄰鍵長——你的鍵結鄰居全擠在那裡。過了它,g(r) 陷入一道谷,再升成一個較寬、較矮的第二峰(第二鄰殼層),又鼓成一個更平緩的第三峰,而這些起伏持續衰減。在幾個鍵長之內,g(r) 就已壓平到 1 並停在那裡:距離一旦夠遠,玻璃就忘了你在哪,每個距離都只顯示平均密度。這些振盪褪為 1,正是「沒有長程有序」的數學面容——相較之下,在晶體裡,峰會保持銳利、永不消亡,一個接一個精確地標出殼層,直到無窮遠。

再做一步修飾,g(r) 就變成這個領域真正拿來畫圖的量。介於 r 與 r + dr 之間那層薄球殼裡的原子數,是 4 乘 pi 乘 r^2 乘 rho_0 乘 g(r) 乘 dr,其中 rho_0 是平均數密度(單位體積內的原子數)。函數 RDF(r) = 4 pi r^2 rho_0 g(r) 才是嚴格意義上的徑向分布函數——字面上就是在每個距離找到的鄰居數目。注意那個 4 pi r^2 的球殼因子:即使在 g(r) 已安定為 1 的遠處,RDF 仍會隨 r^2 一路攀升,單純因為越大的球殼裝得下越多原子。所以 g(r) 是擁擠計(相對平均的密度),而 RDF 是把 g(r) 重新加權成「每單位距離的老實人頭數」。兩者攜帶同樣的資訊;你要畫哪一個,只是方便與否的問題。

RADIAL DISTRIBUTION   g(r) = local density / average density

CRYSTAL (long-range order): sharp spikes at exact shell radii, forever
  g(r)|   |       |         |            |             |
      |   |       |         |            |             |
    1 |---|-------|---------|------------|-------------|---
      +---+-------+---------+------------+-------------+--> r

GLASS / LIQUID (short-range order only): peaks broaden and damp to 1
  g(r)|    __
    2 |   /  \      __
      |  /    \   _/  \__          _____________________
    1 |-/------\-/------\___.---''                        g -> 1
      |/        V           (no correlation: average density)
    0 |__
      +----+---------+----------+-----------+-----------+--> r
         r1        ~2 r1
         ^
  first peak:  centre = nearest-neighbour bond length
               area   = coordination number
同一條曲線一眼就分辨得出晶體與玻璃。晶體的 g(r) 是一排在精確殼層半徑上、永不消亡的銳利尖峰(長程有序)。玻璃的 g(r) 有一個銳利的第一峰(短程有序完好無缺),但外側諸峰在幾個鍵長之內就變寬、衰減到 1——週期性沒了。第一峰的峰心就是鍵長;它的面積就是配位數。

讀第一個峰:鍵長與配位數

第一峰是座小金礦,因為它有三個特徵,各自回報一件不同的物理事實。它的「峰心」是最近鄰鍵長——你的鍵結鄰居所坐落的平衡間距。它的「寬度」量的是這些鍵有多分散:一部分來自靜態無序(玻璃是帶著一整段略有差異的距離凍結下來的,而非某個精確值),一部分來自熱振動,把每顆原子在其平均位置附近抹糊。而它的「面積」——RDF 第一峰底下所收攏的原子數——就是配位數,最近鄰的平均數目。這正是你在密堆積那一階為晶體所認識的同一個配位數(面心立方給 12,鑽石立方矽給 4),只不過現在是為一個「毫無晶格」的固體量出來的。

用非晶質矽(a-Si)給它填上真實數字。它量到的 g(r) 把第一峰放在 r1 = 2.35 埃——在峰寬的誤差內,與晶態鑽石立方矽中的 Si–Si 鍵一模一樣。把 RDF 的第一峰積分到它的第一個極小處,面積約為 4。所以玻璃裡的每一顆矽原子,仍有 4 個最近鄰排成一個四面體,和完美晶體裡完全相同:短程有序基本上原封不動。失去的東西住在更外側。第二鄰峰明顯變寬,第三鄰峰被抹向 g -> 1 的背景,再往外曲線就平了。這一條曲線,就是整個非晶質態的一幅全景——短程有序完好、長程有序不存——而它,正是任何 a-Si 的連續無規網絡模型必須重現的標靶,也就是下一篇的主題。

結構因子 S(Q):偵測器真正看見的

一個尷尬的事實在此:沒有任何儀器直接量到 g(r)。偵測器記錄下來的,是散射強度隨角度的分布圖——那些寬闊的光暈對角度的變化。經過打包與歸一化,這張圖就是靜態結構因子 S(Q),寫成散射向量長度的函數,Q = 4 pi sin theta / lambda。(這個 Q 不過是你在繞射那一階認識的散射向量長度 s = 2 sin theta / lambda 的 2 pi 倍——同一支箭頭,物理學家的單位。)對玻璃而言,S(Q) 是一條平滑曲線,帶著幾個寬闊的隆起——又是那些光暈——它們振盪,然後在高 Q 處安定到 S(Q) -> 1。這個高 Q 極限的 1,是實空間裡 g -> 1 極限在倒空間的孿生子:兩者都只是材料在說「這麼遠已經沒有任何關聯留下」。

現在來看美妙的部分,而它是一位穿了新衣的老朋友。g(r) 與 S(Q) 是一對傅立葉變換對——正是那把曾將晶體變成其倒晶格的「實空間對倒空間」二重性,如今為一個沒有晶格的材料效力。一次正弦傅立葉變換,就把量到的 S(Q) 轉成實空間的徑向分布函數:示意地說,G(r) 正比於「Q 乘 [S(Q) 減 1] 乘 sin(Q r)」對 Q 的積分。於是實驗住在 Q 空間,你在那裡只看得到寬闊模糊的光暈,而一次變換就把 r 空間的圖景交到你手上——銳利的鍵長、可數的配位殼層——那是你真能拿來推理的東西。這整套流程「量總散射、再傅立葉變換成 g(r)」,正是對分布函數分析,玻璃、液體與奈米結構物質的日常工具。

  1. 收集總散射。用 X 光或中子照射玻璃或液體,記錄強度對角度——一條由寥寥幾個寬闊光暈構成的平滑曲線,任何地方都沒有銳利峰。
  2. 把角度換算成 Q = 4 pi sin theta / lambda,再做校正與歸一化(扣掉背景、除去原子散射),得到結構因子 S(Q)。
  3. 把 S(Q) 減 1 做正弦傅立葉變換到實空間,得到 g(r),再乘上 4 pi r^2 rho_0 得到 RDF,即每個距離上的鄰居數。
  4. 讀第一峰:它的峰心是最近鄰鍵長,它的寬度是鍵距的分散程度(靜態無序加上熱振動)。
  5. 把 RDF 的第一峰積分到它的第一個極小——那面積就是配位數,即最近鄰的平均數目。

關於 S(Q) 給了什麼、又藏了什麼,用兩則誠實的收尾。最靠低 Q 的第一個隆起——常被稱作「第一銳利繞射峰」——是中程有序的指紋:那是伸展到最近鄰之外的關聯,例如二氧化矽中共角四面體那鬆散的網絡重複,大約出現在 Q 約 1.5 每埃處。那是第一殼層之外的真實結構,而 S(Q) 看得見它。但現在來看那個深刻的極限:由於玻璃是各向同性的,你量到的永遠只是一條經球面平均的、一維的曲線。變換到 g(r) 會丟掉每一絲方向資訊——鍵「角」、環的統計、真正的三維堆積。真正不同的三維排列,可以共用完全相同的 g(r)。所以 g(r) 與 S(Q) 緊緊「約束」一個結構,卻永遠無法唯一地「釘死」它——這是你在晶體那裡遇過的相位問題,在無序世界的表親。(還有個較小的實務皺褶:你只能把 Q 量到某個有限的上限,在那裡截斷變換,會把一些小小的假漣漪撒進 g(r) 裡,稱為截斷漣漪。)

每個玻璃模型的裁判

正因為 g(r) 是我們真能量到的東西,它就成了裁判:任何無序固體的結構模型,唯有重現觀測到的 g(r) 與 S(Q),才算可信。接下來兩篇會建起兩大模型家族,而兩者都在這條標準線前受審。對二氧化矽這類共價玻璃,查哈里亞森的連續無規網絡取用剛性的 SiO4 四面體——矽扮演網絡形成者——並透過共用的橋氧,角對角地把它們連成一張非週期的三維網。每個四面體內部的 O–Si–O 角維持在約 109.5 度的銳利值,但連接四面體之間的 Si–O–Si 角則廣泛分散。正是這個組合,讓 g(r) 呈現一個極銳利的第一峰、卻是衰減的外側諸峰。摻進一個像 Na2O 這樣的網絡修飾者,它會扯斷其中一些橋,製造出非橋氧,使網絡鬆動——這就是為什麼鈉鈣玻璃遠比純二氧化矽更易熔化、更好加工背後的結構化學。

至於金屬玻璃,勝出的圖景全然不同:硬球的無規密堆積,原子像倒進袋子、搖實後的滾珠般彼此卡緊,配位數接近 12 到 13,帶著濃厚的二十面體局部有序氣味——但任何地方都沒有晶面。它的 g(r) 帶著一個洩底的印記:一個特徵性地「分裂」的第二峰,那種「一肩一鼓」的形狀,是密集無規堆積會顯示、而簡單晶體不會有的。當你在量到的曲線裡看見那個分裂的第二峰,你正直視著一個密集無規堆積的金屬玻璃,而非細晶粒的晶體——這是純粹從 g(r) 的形狀下的診斷,也就是第五篇的全部故事。

RDF 鋪好、卻無法獨力回答的最後一問:這一切究竟為何會凍結下來?一個被冷卻得夠快、閃過結晶的液體,會成為過冷液體;隨著它變冷,黏度飆升,而在玻璃轉變處,原子已無法在實驗的時間內重新排列——整個結構在動力學上鎖死。這是一種動力學的凍結,而非熱力學的相變:玻璃轉變溫度會隨冷卻速率漂移,而自由體積(原子要相互擠身而過所需的那一點點餘裕)在它以下被凍結。金屬結晶得太積極,以致最早的金屬玻璃需要每秒約 10^6 克耳文的快速淬火。把勝算撥向玻璃的,是阻挫:當一個液體最偏愛的局部堆積——二十面體,連同它被禁止的五重對稱——根本無法週期性地鋪滿空間時,結晶就在幾何上受挫,熔體於是滑進玻璃,而非找到一個晶格。