贝叶斯优化在低能电子衍射表面结构分析中的应用

1. 低能电子衍射与表面结构分析的挑战

低能电子衍射（LEED）技术自上世纪60年代发展至今，已成为表面科学领域不可或缺的表征手段。作为一名长期从事表面物理研究的实验人员，我深刻体会到LEED在解析表面原子排列方面的独特价值。当一束低能电子（通常20-500eV）入射到晶体表面时，会与表层原子发生弹性散射，形成特定的衍射图案。这些图案的强度随电子能量变化的曲线——即LEED-I(V)曲线，蕴含着丰富的表面结构信息。

1.1 LEED技术的核心优势

与传统X射线衍射相比，LEED具有几个不可替代的特点：

• 表面敏感性：低能电子的平均自由程仅约5-10Å，这意味着信号主要来自最外几层原子，非常适合研究表面重构现象。在我的实验中，即使吸附单层气体分子，也能观察到明显的衍射图案变化。
• 定量解析能力：通过测量多束衍射斑的I(V)曲线并进行动力学计算，可以精确确定表面原子的三维坐标，精度可达0.01Å。这种精度在研究表面弛豫、重构层形成等细微结构变化时至关重要。
• 实验室级设备需求：与同步辐射光源等大型设施相比，LEED设备可以安装在常规实验室，使得表面结构研究更加便捷。我们实验室的LEED系统就与分子束外延（MBE）设备联用，实现了生长过程的原位监测。

1.2 传统分析方法的局限性

尽管LEED在原理上非常强大，但其定量分析过程却异常复杂。主要困难来自：

多重散射问题：与X射线主要经历单次散射不同，低能电子会在表面原子间经历多次散射。这种非线性过程使得I(V)曲线与原子位置的关系变得高度复杂。我曾计算过一个简单的Cu(111)表面，就需考虑数百条散射路径的干涉效应。

逆问题求解的挑战：从实验I(V)曲线反推表面结构属于典型的逆问题求解。传统方法采用"试错法"：

1. 假设一个结构模型
2. 计算理论I(V)曲线
3. 与实验对比（通过R因子评估）
4. 手动调整参数重复过程

这种方法严重依赖操作者的经验。记得在分析一个简单的Ag(100)表面时，即使使用ViPErLEED.calc软件，也需要反复调整20多个参数，整个过程耗时近两周。

参数耦合问题：表面原子位置、热振动振幅（VIBROCC）、入射角等参数相互影响。在Fe2O3(1102)表面的案例中，53个自由参数形成的高维空间存在大量局部极小值，常规优化算法极易陷入其中无法自拔。

2. 贝叶斯优化框架的设计与实现

针对上述挑战，我们开发了基于物理约束的贝叶斯优化（BO）框架。这个方案的核心思想是将LEED逆问题转化为概率推理过程，通过智能采样策略高效探索参数空间。

2.1 整体架构设计

我们的系统架构包含三个关键模块：

物理计算层：

• 保留完整的ViPErLEED和TensErLEED计算引擎，确保散射物理的严格描述
• 每次迭代都进行全动力学计算，避免简化带来的误差
• 支持并行计算，单次I(V)曲线计算时间控制在1-5分钟（取决于表面复杂度）

概率建模层：

class GaussianProcess: def __init__(self, kernel='Matern2.5'): self.kernel = kernel # 使用Matern 2.5核函数 self.noise_level = 1e-4 # 考虑计算噪声 def fit(self, X, y): # 实现高斯过程回归 self.model = SingleTaskGP(X, y, covar_module=self.kernel) def predict(self, X_new): return self.model.predict(X_new)

优化控制层：

• 采用自适应信任域（TR）策略，动态调整搜索范围
• 集成qUCB（并行上置信界）采集函数
• 实现热重启机制，避免局部最优

2.2 关键技术实现细节

高斯过程建模： 我们选择Matérn 2.5核函数构建协方差矩阵，其数学形式为： k(x,x') = (1 + √5d + 5d²/3)exp(-√5d)，其中d为归一化距离

这种核函数对参数空间的非平滑性具有较好的适应性。在实际编码中，我们使用BoTorch库的实现，并对超参数施加适度的先验约束，防止过拟合。

信任域管理： 信任域的动态调整是算法成功的关键。我们设计了一套启发式规则：

• 当连续5次迭代R因子改善>5%，扩大TR半径20%
• 当连续10次迭代无显著改善，收缩TR半径30%
• 最小TR半径设为参数范围的1%

并行采样策略： 对于高维问题（如Fe2O3案例），我们采用Thompson采样与多TR结合：

1. 维护4个独立TR区域
2. 每个区域使用不同的随机种子初始化
3. 通过评估点的质量动态分配计算资源

3. 应用案例与性能分析

3.1 Ag(100)-(1×1)表面重构

作为概念验证，我们首先测试了简单的Ag(100)表面。这个体系虽然只有13个优化参数，但足以展示算法的核心优势。

优化过程解析： 图2展示了典型的收敛轨迹：

• 初期（0-100步）：算法在大TR内广泛探索，R因子从0.27快速降至0.15
• 中期（100-280步）：TR逐步收缩，聚焦有希望的区域
• 后期（>280步）：精细调节VIBROCC参数，R因子收敛至0.077

特别值得注意的是，原子在xy平面的位移主导了早期的R因子下降（图2c），这与DFT计算的能量降低趋势一致（图3a）。这种一致性验证了算法的物理合理性。

关键发现：

• 表面热振动参数不可忽略：固定VIBROCC为体材料值会使R因子恶化5倍（图3c Case C）
• 存在近简并解：多个结构差异微小（<0.01eV）的模型都能很好拟合实验数据
• 高阶衍射束更敏感：(41)束的R因子降幅最大，达50%（图2d）

3.2 Fe2O3(1102)-(1×1)表面重构

这个复杂案例涉及53个自由参数，包括：

• 27个原子的位移（z>0.25）
• 2种原子的VIBROCC
• 电子束入射角θ

算法应对策略：

1. 采用4个并行TR区域
2. 定义有效采样半径Reff监控探索程度（公式8）
3. 在约100步时检测到局部最优，通过扩大Reff成功逃脱（图4b）

性能指标：

• 最终R因子：0.1977（优于手动优化的0.201）
• 总计算次数：约500次（相当于4天计算时间）
• 入射角收敛值：0.874°（与实验估计一致）

4. 技术优势与扩展应用

4.1 与传统方法的对比

我们总结了BO框架的几大优势：

效率提升：

可靠性增强：

• 通过TR机制避免物理不合理区域
• 热重启策略降低陷入局部最优风险
• 不确定性量化指导采样重点

4.2 潜在应用扩展

这套框架可推广到其他表征技术的逆问题求解：

X射线衍射：

• 适用于薄膜/界面结构分析
• 需要调整散射因子计算模块
• 可处理异常散射等特殊情形

中子衍射：

• 用于磁性材料结构确定
• 需考虑核散射与磁散射的干涉
• 对氢原子定位特别有效

电子显微镜：

• 解析原子分辨率STEM图像的原子位置
• 需开发新的正向模型替代LEED计算
• 可结合深度学习进行特征提取

5. 实操建议与经验分享

基于数十次实际测试，总结出以下关键经验：

5.1 参数设置指南

GP超参数选择：

• 长度尺度先验：设为参数范围的20-50%
• 噪声水平：初始设为1e-4，根据收敛情况调整
• 采集函数：低维用qUCB（β=0.3），高维用Thompson采样

计算资源分配：

• 每个TR区域分配至少4个CPU核心
• 设置检查点每50步保存进度
• 内存需求约2GB/参数（53参数案例需100GB）

5.2 常见问题排查

收敛缓慢：

• 检查TR半径是否过小
• 确认GP超参数不会导致过度平滑
• 考虑增加并行TR数量

物理不合理结果：

• 检查参数边界约束是否恰当
• 验证正向计算（如ViPErLEED）的输入文件
• 查看是否有异常噪声干扰实验数据

数值不稳定：

• 标准化所有输入参数到[0,1]范围
• 对R因子加小偏移（如1e-6）避免除零
• 使用双精度浮点运算

这套系统已在GitHub开源（项目地址需审核），包含完整的案例教程和数据处理脚本。对于希望自主开展研究的同行，建议从Ag(100)案例入手，逐步过渡到更复杂体系。在实际操作中保持耐心，因为即使采用BO方法，复杂表面的完全收敛也可能需要数天时间。