hklstrain：一种基于劳埃晶体学方法计算取向依赖微应变（micro-strain）的算法

《Crystals》：hklstrain: An Algorithm for Orientation-Dependent Micro-Strain Calculation Using Laue Crystallography Methodology

　　

利用劳埃晶体学方法精确表征结晶材料内部微应变分布具有重要意义。基于Green–Lagrange变形矩阵的应变张量在表示hkl晶格应变张量时往往精度不足且评价标准存在缺陷，尤其对于非立方晶系。为此，研究人员设计了一种新型算法hklstrain。该算法依据理论标定的斑点位置获取无应变晶面间距d0，并根据劳埃X射线反射原理与取向信息通过全局拟合方法获取有应变晶面间距dmes，进而将所有关注晶面的微应变映射到样品坐标系。hklstrain应变分辨率可达2×10?4，可定量确定取向依赖的微应变场，为结晶材料表征提供直观认识，已成功应用于上海同步辐射装置（SSRF）BL03HB光束线进行微应变分析。与Green–Lagrange应变张量法相比，hklstrain能给出可靠的取向依赖晶面间距变化微应变图，在劳埃晶体学结晶材料应变张量分析中具有显著优势。

研究背景方面，结晶材料作为航空航天、能源与高性能电子系统的关键构件，其内部应变状态直接影响力学性能、服役寿命与失效行为。传统应变测量手段如电阻应变片仅获宏观平均信号，电子背散射衍射（EBSD）只能辅助分析，扫描电镜（SEM）关联取向与应变但需样品预处理，透射电镜（TEM）虽达原子尺度却制备繁琐且限于微小体积，数字图像相关（DIC）缺取向信息且仅表面测量，拉曼光谱受限于强拉曼信号材料。劳埃晶体学方法具备无损、原位、多环境适用优势，但现有LaueTools、XMAS、XtalCAMP、PYXIS等软件均基于Green–Lagrange变形矩阵，将该张量映射至hkl晶格应变对非立方晶系困难且不便。因此研究人员在SSRF BL03HB光束线开发开源算法hklstrain，以劳埃原理与数值优化克服上述缺陷。

关键技术方法上，研究人员以正交晶系LaAlO₃(100)单晶、镍基单晶高温合金叶片及MgO单晶为样本，在BL03HB使用Pilatus 2M探测器与五维样品台进行固定几何角二维面扫描（步长50μm×50μm，曝光0.2s），上游Laueprocess 3.9输出理论/实测斑点位与初始取向欧拉角；hklstrain分两支处理：一支以globalfitting在单胞体积不变约束下拟合实测衍射方向(θ,φ)_mes得d_mes,i，另一支由hkldcal算d_hkl,i、hklecal与hklmicroecal由ε_hkl,i=(d_mes,i?d_hkl,i)/d_hkl,i得单斑应变，再将晶面法向分解至样品系用最小二乘求ε_xx/ε_yy/ε_zz（要求斑点i>5），取向由欧拉角经旋转矩阵R得ORI_[100]/ORI_[010]/ORI_[001]并以K-means聚类可视化；波长标定另用四晶单色器能量扫描2500帧交叉验证。

研究结果部分，Miller Indices Assignment and Validation小节中，LaAlO₃(100)劳埃图全部斑点经Laueprocess标定hkl，hklstrain理论波长与能量扫描中心波长一致，证实标定可靠。Orientation Classification小节中，由方程(7)计算的ORI_[100]呈中下部聚类与局部梯度，ORI_[010]左下象限有明显方向梯度，ORI_[001]近乎均匀，说明空间异质主要来自[010]轴变化且[001]向应变极低。Micro-Strain Analysis小节显示，总应变ε图全局近零；ε_xx在中心区有压缩带（负应变），ε_yy左下象限有更高幅拉伸带，ε_zz全图均匀蓝（近零），证实即便整体无变形，[100]/[010]向仍存在加工残余应力导致的局域微应变分量。Compare with Different Methods小节指出，LaueTools/XMAS等因商业闭源或标定失败未生成Green–Lagrange矩阵；hklstrain对MgO单晶达?2.21×10^?4应变、RMS 1.9×10^?5 rad（分辨率2×10^?4），对带涂层镍基叶片得应变1.395×10^?3、RMS 8.95×10^?4 rad，验证原位适用性。Further Development小节阐明局限：仅适结晶材料、表亚表层数毫米、多相重叠斑需人工干预、严重宽化斑待后续版本处理。

讨论与结论翻译：研究人员开发并应用于SSRF BL03HB的开源hklstrain算法具以下特征——(1)以劳埃X射线反射与数值优化测d_mes；(2)易适配其他光束线或软件；(3)单张劳埃衍射图即得全微应变张量；(4)利用所有匹配劳埃斑评估微应变；(5)以残差均方根（RMS）作数值优化判据；(6)应变分辨率达2×10^?4；(7)滤除伪斑，适用于带涂层单晶；(8)检出整体Green–Lagrange变形矩阵会遗漏的hkl晶面间距变化；(9)将应变数据关联取向信息，助研究人员直观识别取向差热点并理解应变产生物理机制。该文发表于《Crystals》。