Recovering full crystallographic orientation in specimen coordinates from polarised μ-FTIR spectra

技术摘要：从极化 $\mu$-FTIR 光谱中恢复完整的晶体学取向

问题背景
傅里叶变换红外（FTIR）光谱法是定量分析名义无水矿物（NAMs）中痕量羟基（OH）的标准方法，因其易于获取且速度较快。然而，双折射晶体会各向异性地吸收红外光，这意味着测得的吸光度取决于晶体相对于极化矢量的取向。虽然极化 FTIR 理论上可以同时解析 OH 含缺陷的浓度和取向，但这需要预先知道晶体相对于仪器的取向。

目前的实践面临显著局限性：

1. 非极化 FTIR：大多数研究使用随机取向晶体的非极化光谱。这种方法丢弃了取向信息，并会导致系统性地低估 OH 浓度，除非测量数百个颗粒，但这通常并不实际。
2. 预取向：对预取向晶体进行极化光谱测量需要测量前对齐样品，这对于细粒材料而言是一个障碍。
3. EBSD 耦合：将 FTIR 与电子背散射衍射（EBSD）结合使用解决了取向问题，但需要两台仪器，在转移样品时会引入对准误差（通常为 $\pm 8^\circ$），并且是平均化取向数据而非保持其与光谱的内部一致性。
4. 现有的光谱方法：以往尝试直接从光谱中恢复取向的方法（例如 Asimow 等人，2006 年）存在局限性：它们没有将晶体置于试样参考系中（使得与其它技术的比较变得困难），依赖于可能受到检测器非线性影响的固定光谱范围，且缺乏对其他极化技术的通用性。

方法论
作者提出了一种利用极化微米级 FTIR ($\mu$-FTIR) 光谱，在试样坐标系下恢复各向异性晶体完整晶体学取向（以欧拉角表示）的方法。其核心创新是“波长截面法”（wavelength-section method），该方法重建了晶体各向异性响应的二维截面。

• 理论框架：该方法基于透射率方程（Asimow 等人，2006 年），该方程将透射率与晶体主轴、样品厚度和极化矢量的取向联系起来。不同于以往解决极化矢量相对于晶体轴向的方法，本方法求解的是晶体相对于固定外部（显微镜）坐标系的取向。
• 参考系：作者定义了三个耦合的参考系：样品系、显微镜系（其中极化器固定在 x 轴方向）和晶体系（由晶胞定义）。通过寻找能够使测得的透射率数据与源自三个主轴方向的理论“透射率包络线”相匹配的旋转（使用 Bunge $zxz$ 欧拉角）来确定取向。
• 数据采集：该方法需要在单个诊断波长下，收集在不同极化器角度 ($\alpha$) 下的极化光谱。每个 $(T, \alpha)$ 对被转换为透射率包络线上的一个笛卡尔点。
• 优化：算法通过迭代应用试探旋转，使测量点与理论包络线之间的失配最小化。为了处理取向代价曲面的非凸性质（存在多个局部极小值），作者测试并建议使用带有多个随机初始猜测值的梯度优化器或微分进化全局优化器。
• 波长选择：该方法要求选择一个透射率包络线具有足够各向异性（主透射率差异显著）且样品厚度能使透射率保持在检测器线性响应范围内（通常为 10–90% T）的单一波长。作者通过展示橄榄石在二氧化硅倍频区（1500–2100 $\text{cm}^{-1}$）的应用证明了这一点。

主要贡献

1. 完整取向恢复：该方法恢复的是试样坐标系下的完整晶体学取向（欧拉角），而不仅仅是极化矢量相对于晶体轴的取向。这使得直接与 EBSD 和其他结构数据进行比较成为可能。
2. 单波长灵活性：通过在单一波长下运行，该方法避免了对宽光谱范围的需求，允许用户避开检测器非线性或基线问题。
3. 开源实现：作者提供了 FTIRkit，这是一个实现了该方法的开源 Python 工具，依赖于标准的科学库（NumPy, Pandas, SciPy, Matplotlib）。
4. 验证：通过蒙特卡洛模拟和已知取向（通过 EBSD 测量）的橄榄石晶体实验数据对该方法进行了验证。

结果

• 模拟性能：在无噪声模拟中，该方法达到了极高的精度。使用带有 20 个初始猜测值的梯度优化器，在 >99.75% 的情况下，恢复的取向与真实值在 1° 以内。中位角度精度约为 0.004°。
• 噪声敏感性：引入高斯噪声后，精度随测量次数（样本量）的增加而单调提高。该方法对透射率噪声比对角度 ($\alpha$) 噪声更敏感。降低透射率噪声可显著提高取向恢复的精度。
• 实验验证（橄榄石）：该方法在来自 Kilbourne Hole 陨石圈的橄榄石颗粒上进行了测试，涉及两种厚度（447.5 $\mu\text{m}$ 和 244 $\mu\text{m}$）。
  • 对于较薄的样品，70% 的颗粒（16/23）在与 EBSD 相比的 25° 失配阈值内被成功恢复。
  • 对于较厚的样品，50% 的颗粒（15/30）达到了该阈值。
  • 成功与特定波长和样品厚度的组合有关，以确保透射率保持在检测器线性响应范围内且包络线形状清晰。
  • 失败通常归因于无法用理论模型重现其测得光谱的颗粒（表明可能存在参考标准、样品夹杂物或厚度变化的问题），而非算法本身的失效。

意义与声明
作者将这项工作定位为解决极化光谱学中“取向问题”的方案。他们声称，通过利用 FTIR 信号本身的各向异性，该方法消除了对第二台仪器（EBSD）的需求，并保持了光谱与取向之间的内部一致性。

作者谦虚地指出，该方法：

• 适用于任何具有主方向参考光谱且拥有相关光谱-取向模型的双折射晶相。
• 可以扩展到其他极化技术（拉曼、紫外-可见光谱、激光直接红外光谱）。
• 为光谱范围法提供了一个更灵活的替代方案，允许用户避开退化的波数区域。

作者并未声称该方法在所有实验条件下都是完美的；他们明确指出，成功取决于选择合适的波长和厚度以维持检测器线性响应，并且现实中的挑战（如样品移动、环氧树脂变形）可能会影响与 EBSD 的实际对比精度。主要贡献在于证明了仅通过极化光谱即可恢复完整的取向，从而消除了常规极化 $\mu$-FTIR 分析的一个主要障碍。