一种用于测试生产级赝势并构建具化学精度赝势的合成赝势基准

《Physical Chemistry Chemical Physics》：A synthetic pseudopotential benchmark for testing production pseudopotentials and constructing pseudopotentials with chemical accuracyOpen Access

　　

范数守恒（norm-conserving）赝势保证全电子与赝波函数产生的静电单极项一致。研究人员引入该概念的一般化形式，同时保证由轨道杂化产生的高阶多极项近似守恒。通过纳入浅核态，研究人员可在密度泛函理论（DFT）下使原子化能误差始终小于所谓化学精度（chemical accuracy）1 kcal mol?1。为展示新赝势精度并评估其他赝势，研究人员生成了高度精确全电子计算得到的合成原子化能基准，该数据库含前六周期65种元素共669个分子。与其他库不同，该基准面向平面波或小波等系统基组，且分子测试集结构多样性高于周期性块体测试集。基于高阶多极守恒赝势，研究人员提出一种将异核分子原子化能误差分解到各元素单个赝势的方案，从而可识别低质量元素赝势。用该方法测试ABINIT、QE等主流代码默认赝势后发现，某些元素误差相当大，且在很多情况下超过PBE交换关联泛函本身带来的误差。

论文解读：《Physical Chemistry Chemical Physics》发表的合成赝势基准研究

研究背景方面，现代计算材料研究高度依赖密度泛函理论（DFT）。Materials Project、AflowLib、OQMD等大型数据库均使用赝势（PSP）替代芯电子以提升效率。然而预测性计算化学要求能量达到化学精度（chemical accuracy，1 kcal mol^?1），现有赝势仍引入超出该限的误差。当前三类活跃基准均只含固体块体体系，以方程-of-状态、晶格常数、体模量及Δ测度比较全电子线性增广平面波加局域轨道（LAPW+LO）参考，难以反映分子多样化学环境中的赝势表现。因此研究人员开展一项兼具化学多样性与全电子精度的合成分子基准，并构建多极守恒赝势（MPC PSP）以逼近化学精度。

主要关键技术方法包括：采用LAPW+LO作全电子参考并在5×5×7 ?至7×7×7 ?小周期盒中放置669个中性分子以兼顾孤立极限与计算成本；构造双空间高斯型赝势并通过拟合原子Kohn–Sham（KS）本征值、单极项、高阶多极项及组态能差实现多极守恒；用原子零阶正则近似（AZORA）与无限阶相对论近似（IORA）处理标量相对论效应；以Γ点采样、密度/势对称化或分数占据避免破缺对称解；用Quantum ESPRESSO（QE）与平面波、BigDFT与小波、MRChem与多小波交叉验证。样本队列为H至Po（省惰性气体）前六周期65元素669分子，其中d元素取假设结构，轻元素几何来自NIST已知实验结构。

研究结果各小节如下。1 多极守恒赝势：普通范数守恒只保单极，研究人员令参与成键的n、?轨道杂化交叉项多极Q_?m在全部电子与赝体系中一致，用拟合流程生成H–Br非相对论与AZORA标量相对论MPC PSP，芯/价划分较生产赝势更深，KS本征值误差≤1 mHa、组态能差≤1 kcal mol^?1。2 偶极矩：对18个小分子，LAPW参考与NIST高斯型轨道（GTO）平均误0.5%、最大3.5%；MPC PSP平均0.4%、最大1.5%，显著优于GBRV、ONCV、HGH、kjpaw、rrkjus、JTH。3 分子集：669分子仅与一、二行元素成键以覆盖宽氧化态与低配位数，小盒周期像间距≥3 ?，Γ点即可获得可与全电子比较结果。4 全电子参考：LAPW+LO取r_MT 0.7–1.2 Bohr、K_max=10 Bohr^?1，与多小波MRChem在周期/自由边界下误差一致，最大差0.02 kcal mol^?1，证明小盒与Γ点充分。5 相对论近似：AZORA与IORA对前三周期分子原子化能差≤0.10 kcal mol^?1，对I、Po最高氧化态亦≤0.1 kcal mol^?1，故标量相对论模型可用于基准与赝势生成。6 原子化能：6.1 MPC PSP在QE中以E_cut^wf=1500 Ry、E_cut^pot=6000 Ry计算，前四周期最大误1.1 kcal mol^?1（CoOF₃中Co⁵⁺），其余低于化学精度；轻元素子集75%误差<0.10 kcal mol^?1。6.2 开源生产赝势：GBRV、ONCV、HGH、kjpaw、rrkjus、JTH、SSSP等超化学精度分子占比20%–59%，且异核误差难归因；研究人员用MPC参考将总误按元素分解，发现某些元素赝势误差大于PBE泛函误差，ABINIT/QE默认集亦有大问题元素。

讨论与结论部分总结：研究人员指出已有固体基准的Δ测度优化不足以保障分子原子化能精度；MPC PSP通过保高阶多极与浅核态把赝化误差压到化学精度内，但代价是平面波截断能远高于常规PSP（如Br需800 Ry对比ONCV 80 Ry）。元素贡献分解方案可诊断生产赝势中特定元素缺陷，显示主流开源族在高位氧化态与重元素上仍不足。结论为：该合成基准填补了分子化学多样性空白，MPC PSP在PBE下对H–Po多数体系达化学精度，且误差分解框架证明某些默认赝势误差已超过交换关联泛函误差，对未来高通量材料筛选与预测计算化学的赝势选型具有重要校正意义。