Experimentally-validated multi-slice simulation of electron diffraction patterns
研究團隊首次將多層法模擬導入第 5 階泰勒展開，並與多晶鋁鎂合金實驗數據對齊，打破傳統完美晶格的技術限制。
- 採用多層法 (MS) 的第 5 階泰勒展開，成功在運算資源與高精度繞射圖樣產出間取得平衡。
- 搭配專屬的等向性畸變校正模型，消除硬體幾何變形，使 MS5 的精準度可媲美主流 BW 演算法。
- 首次透過實驗數據驗證多層法模擬，為未來次微米尺度下包含各種缺陷的晶體表徵提供新工具。
高解析度電子背散射繞射 (HR-EBSD) 技術在次微米尺度的彈性應變測量中扮演著關鍵角色。近期一篇發布於 arXiv 的論文首度將多層法 (Multi-slice, MS) 的第 5 階泰勒展開應用於實際多晶鋁鎂合金的實驗數據比對，成功打破過去該演算法僅限於純理論探討的瓶頸，讓非完美晶體的精確繞射圖樣模擬成為現實。
高解析度 EBSD 應變測量面臨的完美晶格限制
在現代材料科學中，高解析度電子背散射繞射 (HR-EBSD) 技術的快速進步，大幅提升了研究人員在次微米空間解析度下測量彈性應變與評估差排密度 (dislocation density) 的能力。為了從電子顯微鏡擷取到的實驗影像中提取極微小的應變資訊，分析系統需要仰賴高質量的動態模擬圖樣，並透過高精度的交叉相關 (cross-correlation，一種計算兩組數據相似度的數學方法) 演算法來進行圖樣匹配。目前學界與業界最廣泛使用的圖樣模擬技術是 Bloch Wave (BW) 方法。
BW 方法基於量子力學中的布洛赫定理，將電子波視為平面波的疊加，透過求解特徵值方程式來決定波函數的振幅與波向量。對於擁有數百億個原子的完美晶體，這種方法能夠極其準確地預測出 Kikuchi 極 (Kikuchi poles，繞射帶的交會點) 與帶狀圖形的幾何位置與相對亮度，且運算效率極高。然而，BW 方法的數學核心建立在理想狀態的假設上，這意味著它在本質上被嚴格侷限於完美週期性的晶體結構。當材料內部存在微觀缺陷、晶界或局部原子位移時，特徵值方程式的邊界條件將被破壞，BW 算法便無法真實反映這些非週期性結構對電子散射產生的擾動，成為高階材料表徵的一大技術限制。
追蹤電子波逐層傳播的多層法演算法特性
為了解決完美晶格假設帶來的限制，另一種稱為多層法 (Multi-slice, MS) 的模擬方案逐漸受到重視。不同於 BW 方法直接求解整個晶體域的特徵值，MS 方法採用了一種更貼近物理傳遞過程的演算法邏輯。它將待測樣品沿著電子束的入射方向，虛擬切割成無數個厚度僅達原子尺度的極薄切片。當入射電子波穿透材料時，系統會計算波函數在穿越原子層時與靜電位勢發生的相位交互作用，隨後再利用自由空間傳播函數模擬波函數在相鄰薄片間的真空傳遞。
這種逐層迭代的計算方式，優勢在於它能夠動態跟隨電子波在材料內部的演化軌跡。因為每一層薄片的位勢可以被獨立定義，這種機制可以輕鬆引入空位 (vacancies)、間隙原子或複雜晶界等局部不均勻結構，並能呈現出比傳統方法更豐富的繞射細節。不過，由於需要對每一層進行大量的二維傅立葉轉換與卷積運算，MS 方法過去一直伴隨著極高的計算成本。這使得它長期以來被視為偏向理論開發的工具，鮮少有研究將其生成的模擬圖樣直接與真實世界的 EBSD 實驗數據進行系統性的對齊與檢驗。
捨棄高能假設並導入第 5 階泰勒展開式優化
為了讓這項具備潛力的演算法能夠實際應用於實驗數據的解析，Xinke Xiao 等人組成的研究團隊針對傳統的 MS 架構進行了核心數學理論的改良。過去在穿透式電子顯微鏡 (TEM) 的模擬中，MS 通常會採用高能假設 (high-energy hypothesis) 來簡化計算；但 EBSD 涉及的是背散射電子的大角度繞射，低能量電子的波長較長，這種簡化會導致波函數的相位計算出現嚴重失真。
研究團隊決定放棄這項妥協，轉而採用更高階的泰勒展開式 (higher-order Taylor expansions) 來精準逼近單向薛丁格方程式 (forward-only Schrodinger equation)。單向方程式忽略了背向散射波的干涉以降低計算複雜度，而透過高階展開則能彌補截斷誤差。經過反覆的模擬測試，研究人員發現採用第 5 階展開的演算法，即本研究定義的 MS5 模型，可以更精確地捕捉波函數在強烈位勢變化下的非線性相位擾動，避免了低階展開常遇到的數值發散問題。MS5 成功在龐大的計算資源消耗與高精度的繞射圖樣產出之間，找到了最具實用價值的平衡點。
多晶鋁鎂合金實驗驗證與等向性畸變校正模型
確立了高效能的 MS5 模型後，研究團隊隨即面臨如何證明數值模擬能對齊物理觀測的挑戰。團隊選用了在工業界具備代表性的多晶鋁鎂合金 (polycrystal Al-Mg alloys) 作為實驗載體，擷取其高質量的 EBSD 實驗圖樣來作為標定基準。然而，真實電子背散射的探測器是一個平面螢光幕，當球面的繞射訊號投影到平面時，必然會產生徑向的空間變形，直接比對會產生極大的殘差。
為了解決這個硬體層面的光學限制，團隊特別搭配了一套量身打造的等向性畸變校正模型 (isotropic distortion correction model)。這套模型透過多項式擬合，將實驗圖樣中的每一像素重新映射到理想的球面座標系上，精確消除了探測器邊緣的非線性拉伸。隨後再結合標準立體三角形重建 (standard stereographic triangle reconstruction) 演算法，將多晶體的取向分佈進行標準化處理。經過這系列嚴密的空間校正後，MS5 生成的動態模擬圖樣不僅與實驗影像高度吻合，其在帶狀對比度與極點定位上的精度，已經完全具備與傳統 BW 方法一較高下的實力。
解鎖次微米尺度缺陷晶體表徵應用的新可能性
從材料檢測與晶體學分析的宏觀角度來看，這項研究達成了該領域內一項重要的基礎里程碑：首次為 MS EBSD 模擬方法提供了堅實的實驗數據背書。長期以來，從事先進金屬材料或半導體磊晶層研究的科學家，往往受限於傳統模擬工具無法處理局部微觀畸變，導致在分析高應力或高差排密度區域時出現不可忽略的誤差。
MS5 演算法的成功驗證與落地，意味著未來的分析軟體可以不再妥協於完美晶格的預設立場。研究人員將有能力針對特定的奈米級缺陷結構建立模型，並直接生成對應的預測繞射圖樣，進而透過反向運算找出實際材料中的真實缺陷型態。這不僅填補了主流演算法在非完美晶體上的技術空白，更為下一代以缺陷動態表徵為核心的材料分析技術，鋪平了底層的運算道路。
透過捨棄高能假設並採用第 5 階泰勒展開，研究團隊成功將多層法模擬與多晶合金實驗數據精準對齊，為未來處理複雜缺陷晶體的 EBSD 分析奠定了紮實的演算法基礎。