一个领域科学家指挥AI智能体写出量子化学程序:分子HF能量与独立程序吻合到纳哈特里量级
核心概要
vibe-qc 是一个面向分子与固体的开源电子结构程序,其代码几乎全部由 AI 智能体在一位领域科学家的指导下编写。它把分子与周期体系放进同一套 Python 接口和 C++17 数值核心,并附带一套可复用的验证集:分子哈特里—福克能量与独立程序吻合到纳哈特里量级,常规 CCSD(T) 也在小分子上做了对照。验证集把独立程序比对、已发表基准和内洽检验分开记录,未达容差的比较同样保留。
深度剖析
一个通用电子结构程序现在可以主要由 AI 智能体编写,而由领域科学家负责选择物理模型、判断结果和批准发布。 电子结构方法的实现与维护长期依赖学生和博士后手工完成,程序往往比写它的人存活更久,因此社区库被用来减少重复实现。 首个哈特里—福克原型在一天后即与独立程序做了对照,较宽的初始实现约用十周完成,此后验证与修订持续进行。
程序把分子与周期计算放在同一套接口下,分子侧覆盖平均场、微扰、耦合簇、活性空间与激发态,周期侧只做到哈特里—福克和 Kohn–Sham DFT。 数值核心用 C++17 承担积分驱动的 Fock 构建、晶格求和、导数核与交换关联网格,Python 负责计算设置、文件读写与流程控制。 积分、泛函、晶体对称性、有效核势积分和傅里叶变换分别调用 libint、libxc、spglib、libecpint 与 FFTW,分子比对以 orca 6.1.1 和 openmolcas 为参照,周期高斯计算以 pyscf 和 crystal 23 为参照。
启示与展望
已演示的范围是分子哈特里—福克、Kohn–Sham DFT、MP2、CCSD(T)、活性空间与线性响应激发态的小体系比对,以及周期哈特里—福克和 DFT 的 GDF 与 BIPOLE 路线;GPW/GAPW、变胞优化、公开驱动下的周期有效核势和相关周期方法不在已验证范围内。
作者明确说明约十周的开发进度只是单个项目的观察,不是对生产率的受控估计;局部相关仍是试点实现,其截断误差单独报告;同机计时中 vibe-qc 在过氧化氢的 DLPNO-CCSD(T) 计算上比 orca 6.1.1 慢,作者称实现仍需优化。
