Physics-aware global Rietveld refinement for high-energy X-ray diffraction microscopy with application to reconstructing intragranular orientation and strain fields

为了了解固体材料的内在生命，科学家们通常会观察它们是如何由内而外构建起来的。大多数金属和陶瓷并非平滑、均匀的整体；它们是由无数微小的晶体（称为“晶粒”）紧密堆积而成的马赛克。当这些材料受到挤压、拉伸或加热时，应力并不会均匀分布。相反，应力会集中在特定位置，通常就在一个晶粒与另一个晶粒交界的地方。这些微观层面的相互作用决定了一座桥梁能否屹立不倒，一个涡轮叶片能否幸存，或者一个陶瓷杯是否会破碎。几十年来，研究人员一直利用强大的 X 射线在不破坏材料的情况下窥视其内部，创建晶粒取向的三维图谱。然而，这种传统的视角在细节处理上显得模糊不清。它能告诉科学家晶粒的平均方向，却无法清晰地看到单个晶粒内部发生的细微形状和应力变化。其结果是，得到的图像往往过于粗糙，无法解释最关键的行为，例如裂纹是如何开始形成的。

加州理工学院的一个研究小组开发了一种新的方法来锐化这种图像，将模糊的快照转变为应力与应变的超高清电影。他们的这种方法（在最近的一项研究中有所描述）将整个材料视为一个单一的、相互连接的系统，而不是一系列分离的部分。他们不再试图逐个晶粒地拟合数据，而是使用一种方法在计算机上模拟整个 X 射线实验，然后不断调整材料的内部图谱，直到模拟结果与现实世界的数据完美匹配。至关重要的一点是，他们在计算机模拟中直接植入了物理定律。这确保了重建的应力和应变场不仅是数学上的猜测，而且是遵循固体变形普遍规则的、物理上可能的各种状态。通过这样做，他们现在可以以前所未有的清晰度观察到单个晶粒内部隐藏的应力景观，为观察铝酸氮氧化物等脆性材料的力学行为提供了一个新的窗口。

该团队面临的挑战在于，分析这些 X 射线图像的标准方法往往会丢弃有价值的信息。传统技术通常将探测器上光亮和暗点的复杂模式转换为简单的二元数据——本质上是决定一个点是否存在——然后忽略光的强度。这就像是通过仅仅计算有多少种乐器在演奏，却忽略了每种音符的音量和音调，来试图理解一场交响乐。研究人员意识到，X 射线信号的完整强度（即在探测器上表现出的亮度变化）包含了理解材料内部细微变化的钥匙。为了解锁这一点，他们创建了一种新的算法，将整个模拟的 X 射线模式与整个实验模式进行逐像素对比。他们使用了一种被称为“最优传输”（optimal transport）的复杂数学工具，来衡量两个模式之间的差异。该工具能够找到将一种模式转化为另一种模式的最有效方式，从而让计算机能够精确计算出如何调整材料的内部图谱，使模拟结果与现实相符。

这种方法的独特之处在于，它不仅调整晶粒的取向，还调整它们之间的边界。在以往提高分辨率的尝试中，研究人员通常假设晶粒之间的边界是固定的，这限制了最终图谱的准确性。这种新方法将晶粒的形状视为一个可以变化的变量。它使用一种虚拟力量来推拉晶粒边界，移动它们，直到模拟的 X 射线模式与实验模式完美对齐。这类似于一位雕塑家在精炼一座雕像时，不仅通过平滑表面，还通过改变底层结构来获得正确的比例。通过在演化晶粒边界的同时演化内部应力场，该方法确保了最终的重建既符合 X 射线数据，又符合弹性物理定律。

研究人员通过两条路径测试了这种新方法：首先是使用完全已知的计算机生成数据，然后是使用现实世界的实验数据。在计算机模拟中，他们创建了一个虚拟的铝酸氮氧化物块（一种坚硬且易碎的陶瓷），并对其施加了压缩载荷。由于他们预先知道标准答案，因此可以测量新方法接近真相的程度。结果令人瞩目。新方法恢复内部应力和取向场的精度显著高于现有技术。它正确识别了晶粒边界的位置以及每个晶粒内应力的细微变化，而旧方法在这些领域经常失败或产生误差。研究表明，通过在重建过程中强制执行物理定律，该方法可以找回那些在传统分析的噪声中实际上已经丢失的细节。

当团队将该方法应用于真实的实验数据时，结果同样令人信服。他们使用了从先进光子源（Advanced Photon Source）实验室中采集的压缩铝酸氮氧化物样本的数据。该样本是一个大约 1.4 毫米宽的小块，其中钻有一个孔以模拟应力集中器。研究人员使用聚焦 X 射线束扫描样本，在样本旋转的过程中捕捉了数千张图像。利用他们的新算法，他们重建了材料的三维状态。生成的图谱揭示了晶粒内应力和应变的复杂网络，展示了载荷是如何分布在复杂的微观结构上的。重建过程能够解析出此前不可见的特征，提供了高保真度的视图，展示了材料在压力下是如何发生形变的。

该研究最重要的发现之一是，新方法足够鲁棒，可以纠正初始图谱中的重大错误，例如取向错误的晶粒或形状略有偏差的晶粒。然而，该方法有一个特定的局限性：如果初始图中遗漏了某些晶粒，它无法创造出新的晶粒。算法会演化现有晶粒的边界使其到达正确位置，但如果一个晶粒在初始猜测中完全不存在，则无法在不修改框架的情况下将其恢复。这意味着初始猜测必须足够密集，以包含所有的晶粒，即使它们的形状或取向最初是不完美的。一旦正确的晶粒存在，通过将模拟与真实数据进行比较的迭代过程，模型就可以实现自我修正，移动边界到真实位置，并调整内部应力场以符合物理现实。这种韧性意味着该方法可以应用于广泛的材料和实验条件，而无需为每一个细节都准备一个完美的起点，使其成为研究复杂微观结构的实用工具。

研究还强调了利用完整 X 射线信号强度的重要性。通过利用来自探测器的每一比特信息，而不仅仅是斑点的存在或缺失，研究人员能够提取出以前无法获得的细节水平。这种方法使他们能够高精度地绘制出晶粒内的取向和应力变化。虽然研究中的合成测试假设每个晶粒内的取向是统一的（零莫赛克扩散）以验证核心方法，但能够在较大体积而非仅仅是极小切片内映射这些变化的能力，为理解多晶材料在极端条件下的行为开辟了新的可能性。研究人员证明，他们的方法可以处理包含数百万个数据点的集合，证明了其计算成本对于现实世界的应用是可控的。

这项工作的意义不仅在于更好的图像。通过提供一种准确重建材料内部状态的方法，该方法为设计从航空航天部件到医疗植入物等一切产品的工程师和科学家提供了新工具。准确了解应力如何在材料中集中，有助于预测失效的时间和位置，从而实现更安全、更耐用的设计。能够以高分辨率和物理准确性进行三维观察，正在改变我们对材料科学的理解。它将该领域从对结构的静态观察转变为对材料如何响应现实世界力量的动态理解。

研究人员向公众发布了名为 PARA-X 的软件，允许他人使用并基于其工作进行开发。这种透明度确保了科学界可以针对不同的材料和实验设置测试该方法，进一步验证其效用。研究结论指出，这种具备物理感知能力的方法代表了高能 X 射线衍射显微技术领域的重大进步。它架起了实验室收集的原始数据与材料内部复杂物理现实之间的桥梁，为微观世界提供了一个清晰、详细且符合物理逻辑的视角。随着该技术的完善及其在更复杂场景中的应用，它有望揭示关于固体基本力学的全新见解，帮助解决材料工程中最具挑战性的问题之一。