从透射电子显微镜(TEM)双晶带轴向量测定晶体取向关系(OR)及重合位置点阵(CSL)的算法——基于三极向量(tripolar vectors)与旋转矩阵的新方法

《Micron》：Revisiting Crystalline Orientation Relationship Measurements via Electron Diffraction - Enhancing Efficiency and Accuracy through Post?Processing

　　

摘要：研究人员提出了一种从透射电子显微镜(transmission electron microscope, TEM)实验测得的双晶粒三极向量(tripolar vectors)中精确提取晶体取向关系(orientation relationship, OR)及相应重合位置点阵(coincidence site lattice, CSL) Σ值的系统化算法。该方法利用TEM样品台沿X轴和Y轴的倾转角度将各晶粒的三组低指数带轴(zone axes)重构为晶体坐标系下的单位向量矩阵M与M′，通过变换矩阵T = M′M?1描述两晶粒间仿射变换，并分别采用三种算法获得正交旋转矩阵(rotation matrix)R：(I)直接将T近似为R；(II)基于最小二乘(least-squares, LS)解R = T(TTT)?1/2；(III)基于奇异值分解(singular value decomposition, SVD)令T = UΣVT后取R = UVT。从R中按标准公式提取旋转轴(axis-angle representation)与旋转角θ，并通过CSL理论关联旋转参数与Σ值。研究人员以铝合金中五组相邻晶粒的TEM数据验证了三种算法的可行性，比较了由仿射矩阵差异引起的偏差角，结果表明SVD与最小二乘算法较直接归一化法更贴近理想CSL OR，且SVD避免了对称矩阵开方的计算开销。该方法可从常规TEM倾转数据中无拟合地获取高精度OR与CSL参数，适用于立方晶系低Σ GB的分析。

研究背景与意义

晶体取向关系(orientation relationship, OR)与重合位置点阵(coincidence site lattice, CSL)是表征多晶材料（特别是铝合金等立方晶系）晶界(grain boundary, GB)几何结构与界面性质的核心参数。传统通过TEM电子衍射图谱定性判断OR（如Nishiyama–Wassermann或Kurdjumov–Sachs关系）受限于测量误差及倾转角度容差（可达±5°），难以定量给出旋转轴和旋转角或对应Σ值。已有基于矩阵的方法常忽略仿射变换矩阵的非正交性，或直接依赖商业软件而缺乏透明算法流程。为克服上述不足，研究人员开展了基于TEM倾转角重建三极向量矩阵、通过三种正交化算法获取严格旋转矩阵并从中解析OR与CSL参数的系统研究，成果发表于《Micron》。

主要关键技术方法

研究人员选取铝合金中五个区域五个样品组的相邻晶粒为研究对象，记录每晶粒三组低指数带轴(zone axis) V_ij及对应双轴倾转角度{α_ij, β_ij}，利用倾转几何将底片坐标系向量还原至晶体坐标系并归一化为向量矩阵M（Grain I）与M′（Grain II经同一旋转后），计算仿射变换矩阵 T = M′M^?1，继而分别采用：(I)直接归一化—列向量分别缩放去非正交性得R = TS^?1；(II)最小二乘法—构造对称矩阵 Q = T^TT，R = TQ^?1/2；(III)奇异值分解(singular value decomposition, SVD)—T = UΣV^T，R = UV^T。从所得正交旋转矩阵 R 中提取迹(trace)得 cosθ = (r₁₁+r₂₂+r₃₃?1)/2，反对称部分得旋转轴分量 p_k= (r_ji?r_ij)/(2sinθ)（i≠j），进而对照CSL公式 tan(θ/2) = (n/m)√d 与 Σ = (m2+nd)/α 判定最近似CSL Σ值，并以旋转矩阵差 D = R₂R₁^?1的偏差角评估算法一致性。

研究结果

3.1–3.3 箔朝向(film orientation)与向量重构

研究人员依据TEM双轴倾转几何关系，将实测带轴在荧光屏投影向量结合倾转角 α（绕X轴）、β（绕Y轴）逆推回晶体坐标系方向矢量，经归一化构建3×3向量矩阵M与M′，并验证三向量间夹角与晶体学理论相符。

3.4.1 Algorithm I: No-fitting solution（直接取法）

研究人员将仿射矩阵T的列向量分别除以其模长分离伸缩矩阵S，取R = TS^?1近似为旋转矩阵。由于T偏离正交性（det≠1、列非正交），此法引入最大偏差，但仍可快速给出OR初值。

3.4.2 Algorithm II: least-squares solution（最小二乘解）

研究人员构造 Q = T^TT，通过 Q 的特征值分解求 Q^?1/2，得 R = TQ^?1/2，该R为T的最邻近正交矩阵(proper rotation matrix, det=+1)，消除非正交与非单位模影响，结果较Algorithm I更接近理想CSL OR。

3.4.3 Algorithm III: SVD solution（奇异值分解解）

研究人员对 T 做奇异值分解 T = UΣV^T，取 R = UV^T（必要时调整 det(R)=+1）。此法避免矩阵平方根运算，数值稳定且与Algorithm II所得旋转矩阵及提取的轴角参数在数值误差范围内一致。

3.5 Extracting rotation axis u and angle θ from the rotation matrix

研究人员从正交旋转矩阵 R 按 cosθ = (Tr(R)?1)/2 求旋转角，当 θ≠0° 或 180° 时由反对称部分求归一化旋转轴 u = [p₁p₂p₃]^T，对于Algorithm I 需先分离伸缩矩阵再代入修正公式。

3.6 Find the CLS OR from the experimental data

研究人员将实验所得 (u, θ) 与经典CSL转动公式对照，通过极图叠合判定最匹配 Σ 值，或由轴角参数按 Σ = (m2+nd)/α（m,n互质，α=1,2,4依晶系而定）计算 Σ。

3.7 Comparison between two rotation matrices

研究人员定义差异矩阵 D = R₂R₁^?1，从其迹按相同公式得两OR间的偏差角(deviation angle)，用以量化Algorithm I与II/III结果差异。

4 Results（结果与算例）

以Sample 2为例，Algorithm I得 θ = 141.33°，u ≈ [0.7634 0.5894 0.2641]；Algorithm II与III均得 θ ≈ 139.85°，u ≈ [0.7762 0.5769 0.2543]，二者偏差角约1.5°，表明直接法引入可测偏差。五组样品均显示SVD与LS结果高度吻合并可匹配至特定低 Σ-CSL OR（如 Σ3, Σ9, Σ11等附近），验证了方法的可靠性。

讨论与结论

研究人员指出，传统直接归一化(Algorithm I)虽简便但因忽略仿射矩阵非正交性会产生可测偏差；最小二乘(Algorithm II)与SVD(Algorithm III)均可获得满足正交归一条件的旋转矩阵，其中SVD计算效率更高且数值稳定性好。所提基于TEM三极向量与旋转矩阵的完整流程——倾转重构、仿射变换、正交化(三种可选)、轴角提取及CSL判定——使研究人员能从常规TEM倾转实验数据直接、无拟合地获取定量OR参数与近似CSL Σ值，为晶界特征分布(grain boundary character distribution, GBCD)分析及界面性能关联研究提供可靠的基础数据支撑。该方法适用于立方晶系金属中低 Σ-CSL 晶界的OR测定与验证。