引言从炒菜试错到先算后做传统材料研发的流程有点像炒菜凭经验和直觉配出配方烧一炉测性能不行再来。一个新型电池正极材料从实验室到量产十年是常态。问题在于候选空间太大——仅无机晶体的可能组合就超过 $10^{10}$ 种靠人力一炉一炉烧永远烧不完。机器学习切入的位置很直接用模型预测成分 结构 → 性能把 99% 的候选在虚拟筛选阶段就淘汰掉实验资源集中在前 1%。Google DeepMind 的 GNoME 在 2023 年一次性预测出 220 万种稳定晶体结构相当于把人类已知的稳定无机材料库扩大了一个数量级。这篇文章聊聊这条 pipeline 的技术细节模型怎么选、数据从哪来、以及实战里哪些坑最常见。三条技术路线筛选、生成、自动实验高通量虚拟筛选是最成熟的路线。流程是从结构数据库里拉出候选 → 用 ML 模型预测目标性能形成能、带隙、离子电导率等→ 按预测值排序 → 只对头部候选跑 DFT 精算或实验验证。核心收益是把 DFT 单次数小时的计算压到毫秒级推理。生成式设计更进一步不从已有库筛选而是让模型直接画出新结构。扩散模型是当下主流——MatterGen、CDVAE 这类工作把晶体表示为原子坐标 晶格参数的连续体从噪声出发逐步去噪生成满足目标性能的结构。难点不在生成而在生成的东西能不能合成所以评估时必须叠加合成可行性打分。自动实验闭环是终极目标模型推荐候选 → 机器人实验室自动合成与表征 → 结果回流更新模型。伯克利的 A-Lab 已经在氧化物体系上跑通了这套闭环41 天内自主合成了 41 种新化合物。这类系统的瓶颈往往不在算法而在自动化表征设备的可靠性。模型怎么选图神经网络是主力晶体的自然表示是图原子是节点近邻原子之间的键是边。GNN 因此成为材料 ML 的标配。不同架构的关键差异在如何编码几何信息| 模型 | 几何信息利用 | 特点 | 适用场景 | | --- | --- | --- | --- | | CGCNN | 仅距离 | 开山之作简单稳定 | 数据量小的基线 | | SchNet | 连续径向基 | 平滑处理距离 | 分子/能量面 | | MEGNet | 加入全局状态 | 可预测温度相关性质 | 多性质联合预测 | | M3GNet | 三体相互作用 通用势 | 可做结构弛豫 | 高通量筛选 优化 | | CHGNet / MACE | 等变消息传递 | 精度接近 DFT | 分子动力学、弛豫 |近几年的一个明显趋势是通用原子间势universal interatomic potentialM3GNet、CHGNet、MACE 这类模型在覆盖元素周期表大部分元素的大规模数据上预训练下游可以零样本做结构弛豫、相稳定性判断把原本需要 DFT 的环节整个替掉。代码示例用 PyG 训练形成能预测模型下面用 PyTorch Geometric 搭一个最小可用的晶体图网络预测材料的形成能eV/atomimport torch from torch_geometric.nn import CGConv, global_mean_pool from torch_geometric.data import Data class CrystalGNN(torch.nn.Module): def __init__(self, node_dim92, edge_dim40, hidden128): super().__init__() self.embed torch.nn.Embedding(100, node_dim) # 原子序数嵌入 self.conv1 CGConv(node_dim, edge_dim) self.conv2 CGConv(node_dim, edge_dim) self.head torch.nn.Sequential( torch.nn.Linear(node_dim, hidden), torch.nn.Softplus(), torch.nn.Linear(hidden, 1), ) def forward(self, data: Data): x self.embed(data.z) # z: 原子序数 x self.conv1(x, data.edge_index, data.edge_attr).softplus() x self.conv2(x, data.edge_index, data.edge_attr).softplus() x global_mean_pool(x, data.batch) # 图级读出 return self.head(x).squeeze(-1) # 构图要点: 以 8Å 为截断半径连边, 边特征用高斯基展开原子间距离 # 训练: loss L1(pred, formation_energy), 数据集用 MatBench 基准这套结构在 Matbench 形成能任务上能拿到约 0.03 eV/atom 的 MAE已经接近 DFT 本身的误差水平。要进一步提升可以换 M3GNet 预训练权重做微调。数据从哪来主流材料数据库对比| 数据库 | 规模 | 数据类型 | 访问方式 | | --- | --- | --- | --- | | Materials Project | 15 万 晶体 | DFT 计算性质 | 免费 APIpymatgen | | OQMD | 100 万 条目 | DFT 热力学稳定性 | 免费 API | | AFLOW | 350 万 条目 | 高通量 DFT | REST API | | NOMAD | 1 亿 条目 | 计算 实验混合 | 开放存档 | | MatBench | 基准集合 | 13 个标准任务 | 直接下载 |给新手的建议先用 Materials Project 的 API 拉数据练手pymatgen 库已经把结构