VASP计算CO₂在Cu(111)面吸附能——一个催化剂筛选项目的实战复盘

这个项目接到的任务是筛选一系列Cu基催化剂对CO₂还原的活性，第一步就要算CO₂在各晶面上的吸附能。听起来是个常规活——几何优化、算能量、做差——但实际跑起来，吸附位点的选择、范德华修正的开与关、以及溶剂化效应要不要加，三个变量叠在一起，结果能差出去0.3 eV以上。

0.3 eV对吸附能意味着什么？足以把”强吸附”翻成”弱吸附”，整套筛选逻辑直接崩掉。

模型搭建阶段：吸附位点的选择比想象中更关键

Cu(111)面是fcc金属最密排面，常规思路是遍历四个高对称位点：top、bridge、fcc hollow、hcp hollow。项目初期就是这么做的——3×3超胞、4层slab、底部两层固定、15 Å真空层——在四个位点上各放一个CO₂分子，跑结构优化。

前两个位点的收敛很顺利。top位吸附能-0.12 eV，bridge位-0.18 eV，都在弱物理吸附区间。但fcc hollow跑出了-0.55 eV——这个数值如果成立，CO₂在Cu(111)上的活化就比文献普遍报道的值高出一大截。

问题出在哪里？把优化后的结构拉出来看才发现：fcc hollow的CO₂在弛豫过程中偏离了初始位置，C原子滑到了bridge和hollow之间的过渡区，O-C-O键角从180°弯到了162°。这不是fcc hollow吸附——这是一个非对称的bridge-like构型。

关键教训：高对称点位初始放置不能盲目信。CO₂是线形分子，在金属表面的吸附构型远比原子吸附复杂，初始取向（平躺vs竖直）和C端朝下vs O端朝下都要分别测试。这个项目后来补了8种初始构型，最终确认最稳定的吸附构型是C端朝下的桥位倾斜构型，吸附能-0.21 eV。

能量计算：范德华修正加不加？

VASP的DFT-D3修正对吸附能的影响在Cu(111)这种金属表面上一直是个灰色地带。PBE泛函本身对范德华相互作用的描述很差，但DFT-D3的阻尼函数在金属表面是否过度修正，不同文献的结论并不一致。

这个项目做了两组对照：

差距最大的一组（PBE vs PBE-D3 zero）差了0.15 eV。对CO₂这种弱吸附体系，这个差异是实质性的——无修正的结果几乎可以认为不吸附，加D3后变成明确吸附。

项目最终的选择是PBE-D3(BJ)，依据来自Grimme原始文献中BJ阻尼在金属表面的验证数据集——Cu(111)上苯的吸附能，BJ阻尼与实验值的偏差在0.05 eV以内。optB86b-vdW虽然物理图像更漂亮，但VASP中的实现版本对slab模型的偶极修正兼容性有问题，在这个项目里放弃了。

溶剂化效应：要不要用VASPsol？

CO₂电还原在实际反应中是在水溶液里进行的。隐式溶剂模型VASPsol可以在DFT层面引入介电环境，但它的参数选择——尤其是表面张力参数和溶质腔半径——对吸附能的影响甚至比范德华修正还大。

这个项目对比了三组：真空、VASPsol（ε=78.4，水）、VASPsol（ε=78.4 + 0.1 M K⁺）。结果出乎意料——加溶剂化后吸附能从-0.21 eV变成了-0.08 eV，几乎又推回了弱吸附区间。

这个现象在文献里是有先例的：J. Phys. Chem. C上有一篇2021年的工作专门讨论过，VASPsol对带电表面（如施加外电场的electrochemical界面）的描述尚可，但对中性slab上的中性分子吸附，隐式溶剂模型的极化响应偏强，低估了吸附能。在这个项目里，VASPsol被限定在”定性趋势判断”层面使用——比较不同掺杂体系在溶剂中的吸附能相对高低——绝对值的锚定还是回到真空DFT-D3的结果。

回过头看

CO₂在Cu(111)上的吸附能，这个项目最后报出的值是-0.21 ± 0.05 eV。±0.05 eV这个区间不是为了好看——它是吸附位点（±0.03 eV）、范德华修正方法（±0.04 eV）和slab厚度测试（±0.02 eV）三组收敛性测试的累积结果。

真正有价值的不是那个-0.21，而是这套测试流程：8种初始构型遍历→4组vdW方法对比→3种溶剂环境测试。催化剂筛选不是算出吸附能就结束了——你得知道这个数字是在什么条件下成立的，以及换了条件它能偏多远。