Bridging Atomistic and Continuum Descriptions of Nanoscale Dislocation Loops in Tungsten

核融合反應爐的防護裝甲預計需在極端輻射下運行 2.5 年，表面積高達 850 平方公尺。然而，科學家目前用來研究材料輻射損傷的原子級模擬，極限僅能處理數千個原子與極短的時間尺度。為了解決這個巨大的尺度落差，英國華威大學與原子能管理局的研究團隊提出了一種跨尺度模型，成功將奈米級的原子位移數據與連續體力學結合，證實連續體模型在距離缺陷核心兩倍半徑處，即可精準預測鎢金屬的應力衰減率。

核融合裝甲面臨的跨尺度運算挑戰

鎢金屬（Tungsten）因具有極高的熔點、極低的物理濺射率與良好的導熱性，被公認為核融合反應爐中面對電漿設備防護裝甲的首選材料。它在遭受中子活化後的半衰期相對較短，約在 100 年左右，提高了未來回收利用的可行性。然而，長時間暴露在高能輻射下，鎢金屬的物理與機械性質無可避免地會隨著時間顯著劣化。

在微觀尺度上，輻射損傷主要表現為 Frenkel pair（法蘭克爾缺陷：晶體中原子離開原位形成空位與間隙原子的配對）的累積。這些高能碰撞級聯產生的點缺陷簇，隨後會鬆弛重組成能量更穩定的位錯環與空洞。當位錯數量增加，它們開始相互作用並形成複雜的網路結構。過去科學家常依賴 DFT（密度泛函理論：以量子力學計算多電子系統電子結構的方法）與分子動力學來研究這些現象。

真實的電漿防護設備運作時間以年為單位，表面積廣闊，這遠遠超出了分子動力學所能負載的運算極限。為了預測長時間受照材料的演變，科學家必須依賴能夠處理宏觀現象的連續體模型。然而，連續體模型必須先透過微觀觀察來進行校準，以確保所預測的巨觀物理特性符合現實。

連續體線性彈塑性的力偶極子模型

古典位錯理論從 CLE（連續體線性彈塑性：將材料視為連續介質來分析受力變形的理論）出發，假設總應變可以分解為彈性與塑性兩個部分。研究團隊利用 Volterra 積分方程式，推導出材料內部任何一點受位錯影響而產生的位移場。

為了解決龐大的運算問題，研究團隊針對距離位錯環極遠的「遠場」（far-field）狀態進行漸進式近似展開。當觀測距離遠大於位錯環半徑時，模型將每一個位錯環視為介質中的點狀力偶極子。計算結果顯示，位移場與位錯環所圍成的面積成正比，這個面積本質上反映了產生位錯環所需的額外體積。

將位錯環轉換為力偶極子，大幅簡化了傳統離散位錯模擬中需要沿著整條位錯線進行積分的繁重計算負擔。在這個遠場近似中，總應變等同於彈性應變，且與彈性格林函數的梯度直接相關。由於彈性格林函數會隨著距離的倒數衰減，模型預測遠場位移會依循平方衰減，而應力場則會隨距離的三次方快速衰減。

0.1 dpa 輻射損傷下的遠場預測閾值

模型所推導出的遠場預測，前提是各個位錯環之間的距離必須夠遠。在輻射損傷程度達到 0.1 dpa（每原子位移數）的高強度受損環境中，鎢金屬內部的位錯環密度大約為 $2.8 \times 10^{22} \, \text{m}^{-3}$。在此極端密度下，將位錯環視為系統中相互作用的粒子，可以推算出它們之間的平均間距約為 330 埃。

觀測數據指出，這些位錯環的半徑大多落在 25 埃 到 50 埃 之間。為避免受到缺陷核心微觀物理效應的干擾，模型與原子模擬數據的交互檢驗，必須在距離核心大約 250 埃 處進行。

這個距離確立了連續體模型與原子模擬之間的「分離長度尺度」。只有在這個距離之外，連續體理論才能有效且安全地接手，用以預測大規模缺陷群集對宏觀材料性質的影響，這正是跨尺度銜接成功的關鍵點。

LAMMPS 原子模擬與多重邊界條件測試

為驗證連續體力學預測的準確度，團隊使用原子模擬環境（ASE）庫建構了 BCC（體心立方：一種常見的金屬晶體結構）鎢晶格。模擬過程全程在 LAMMPS 軟體中進行，並採用了多種精準的嵌入原子模型（EAM）勢函數。針對最大型的系統結構，若常規的共軛梯度法無法達到小於 $10^{-6} \, \text{eV}/\text{Å}$ 的力學容差，則會額外啟用 FIRE 演算法強制降級放鬆結構。

模擬設定為一個球體，半徑在 50 埃 到 500 埃 之間變化。為了解決有限尺寸帶來的表面效應干擾，團隊測試了三種外殼邊界條件：固定在原始晶格位置、強制位移至連續體預測位置，以及完全無牽引的自由放鬆狀態。

透過不同尺寸與邊界條件的交叉比對，研究團隊找出了介於位錯核心異常區與外殼表面干擾區之間的「信心區域」。只有在這個純淨的線性彈性區域內，原子尺度的實際位移才會被採納，並與連續體模型產生的應變場數據進行嚴格回歸比對。

證實遠場應力衰減率與近核非物理極端值

結果驗證了連續體模型的預測：流體靜壓與 J2 剪切不變量兩項關鍵應力指標，確實在遠場呈現精準的距離三次方衰減規律。即使鎢金屬具有異向性彈性張量，但因其齊納異向性比率接近 1.0，整體應力場仍圍繞著穿過位錯環中心的軸線呈現優美的圓柱對稱性。

這項對比同時揭示了單一連續體模型的物理侷限。在距離核心極近的區域（半徑小於 12 埃），連續體模型推算出的應力會達到極度誇張的 1000 GPa。然而，實驗測得的鎢金屬體積模量約為 160 GPa，剪切模量約為 310 GPa。

近核區域出現的非物理數值，明確標示了材料連續性假設破裂的界線。跨尺度模型的價值在於，它清楚界定了原子級數據應該在何處終止，而連續體力學應該在何處接手，從而建立起一套計算資源可負擔、同時又精準的物理特性預測系統。

跨尺度模型成功界定了微觀原子與宏觀連續體的交界，這是未來預測大面積核融合設備壽命不可或缺的關鍵運算橋樑。