1. 项目概述：多目标贝叶斯优化在神经形态电子制造中的创新应用
在柔性电子器件制造领域，光固化工艺参数的优化一直是个复杂难题。传统试错法需要耗费大量时间和材料，而单目标优化又难以兼顾多个性能指标的平衡。我们团队开发的这套人机协同多目标贝叶斯优化系统，成功将实验效率提升了3-5倍。这个系统最突出的特点是它能同时处理两类关键信息：定量测量的电学性能数据和定性的人类专家观察评分。
关键创新点：首次将人类专家对薄膜转化状态的直观判断（如图1所示的显微镜观察结果）量化为约束条件，与电学性能指标共同构建优化空间。
2. 核心方法解析：从理论到实现
2.1 高斯过程回归建模原理
高斯过程(GP)作为贝叶斯优化的核心组件，本质上是在函数空间定义的概率分布。对于光固化工艺而言，我们为每个性能指标建立独立的GP模型：
k(x,x') = σ² exp(-||x-x'||²/2l²) + σₙ²δ(x,x')其中x代表工艺参数组合(辐射能量、脉冲数等)，超参数{l,σ,σₙ}通过最大边际似然估计确定。在我们的实验中，使用Matérn 5/2核函数表现出比常规RBF核更好的外推能力。
2.2 多目标优化框架设计
我们采用基于超体积改进的qEHVI（并行期望超体积改进）算法，其采集函数可分解为：
α(x) = ∫[f(x)-τ]·I(f(x)>τ)p(f|D)df
其中τ为当前Pareto前沿，I(·)为指示函数。为平衡探索与开发，我们设计了动态权重机制：前3轮侧重探索(ω=0.7)，后续轮次逐步转向开发(ω=0.3)。
2.3 人机协同机制实现
专家评分系统经过特殊设计：
- -1：明显未转化(图1a)
- 0：部分转化
- 1：完全转化(图1b)
- 特殊标记：烧焦样品(图1c)
这些评分通过sigmoid函数转化为概率约束： Pconstraint = 1/(1+exp(-(f(x)-0.5)/0.2))
3. 实验设计与关键参数
3.1 拉丁超立方采样(LHS)初始设计
初始30组实验参数范围设计基于前期探索实验：
| 参数 | 范围 | 步长 | 物理意义 |
|---|---|---|---|
| 辐射能量(J/cm²) | 1.0-7.0 | 0.2 | 能量输入强度 |
| 脉冲数 | 1-20 | 1 | 处理次数 |
| 脉冲长度(ms) | 1-20 | 1 | 单次作用时间 |
| 微脉冲数 | 1-30 | 1 | 能量调制精细度 |
| 占空比(%) | 20-70 | 5 | 能量间歇比例 |
3.2 性能评价指标
两个关键性能指标的定义和测试方法：
- 介电常数频率色散比(C100Hz/C1MHz)：
- 使用阻抗分析仪(Agilent 4294A)
- 反映介电材料的频率稳定性
- 漏电流对数绝对值(|log(Ileakage)|)：
- 半导体参数分析仪(Keysight B1500A)
- 表征绝缘性能
4. 结果分析与工艺优化
4.1 Pareto前沿演化
经过7轮优化(共65组实验)，获得的Pareto前沿显示：
- 最佳C100Hz/C1MHz：1.10±0.04 (条件#40)
- 最佳|log(Ileakage)|：6.05±0.24 (条件#40)
- 平衡点条件：C100Hz/C1MHz=2.0时，|log(Ileakage)|可达4.5
4.2 SHAP值参数重要性分析
通过SHAP分析揭示各参数的影响规律：
| 参数 | C100Hz/C1MHz影响 | 漏电流影响 |
|---|---|---|
| 辐射能量 | 负相关(-0.43) | 正相关(0.38) |
| 脉冲数 | 弱正相关(0.12) | 强正相关(0.61) |
| 占空比 | 非线性关系 | 阈值效应 |
5. 材料表征与机理研究
5.1 XPS分析揭示转化机制
对Pareto前沿两端样品(#40和#66)的XPS分析显示：
- 金属氧化物峰(530.4eV)占比：
- #40：19%
- #66：6.7%
- 硝酸根残留(407eV)强度：
- #40：0.12cps
- #66：0.35cps
5.2 温度场模拟
使用SimPulse®软件模拟显示：
- #40条件峰值温度：159°C
- #66条件峰值温度：127°C
- 基底温度均<55°C（保障柔性PET不变形）
6. 实操经验与故障排除
6.1 常见问题解决方案
| 问题现象 | 可能原因 | 解决方案 |
|---|---|---|
| 薄膜局部烧焦 | 能量密度过高 | 降低辐射能量10-15% |
| 转化不均匀 | 微脉冲数不足 | 增加微脉冲数至20以上 |
| 介电常数波动大 | 脉冲长度不合适 | 调整至8-12ms范围 |
6.2 参数调整黄金法则
我们发现最优参数往往出现在以下区间：
- 辐射能量：3.5-4.5 J/cm²
- 脉冲数：12-16次
- 占空比：60-70%
- 微脉冲数：18-24个
7. 技术拓展与应用展望
这套方法已成功应用于：
- 柔性忆阻器阵列制造
- 可穿戴生物传感器
- 神经形态计算器件
在典型的氧化物薄膜晶体管制备中，优化周期从传统的3-4周缩短至5-7天，且器件性能标准差降低42%。一个特别有意思的发现是：人类专家对"部分转化"状态的识别准确率(82%)甚至高于初期GP模型的预测准确率(76%)，这凸显了人机协同的价值。
