电化学材料dft计算

结合案例图分析，第一张案例图是三种催化剂的HER极化曲线，横轴为相对可逆氢电极（RHE）的电势，纵轴为电流密度。灰色undoped曲线在−0.3 V左右才出现明显电流，且达到−40 mA cm⁻²需要更负的电势；蓝色N-doped曲线整体左移，说明N掺杂已经激活了部分活性位；红色defect-rich曲线左移最多，在−0.19 V附近就达到了−10 mA cm⁻²，对应的过电位η₁₀=190 mV。曲线在高电流区都呈现类似的Tafel斜率趋势，说明速率决定步骤没有改变，改变的是氢吸附自由能ΔG_H*的强度。

第二张案例图是d轨道与s轨道的投影态密度（PDOS）。红色实线代表d states，蓝色虚线代表s states。图中标注d带中心ε_d=−1.5 eV，相对于费米能级位置适中。按照Sabatier原理，ε_d过于接近费米能级会过强地吸附H，过于远离则吸附太弱；−1.5 eV这个位置通常对应中等强度的H*吸附，有利于HER。图内阴影区表示费米能级附近的态密度贡献，d states在−2 eV附近有主峰，说明活性位点的d电子主导了与H 1s轨道的杂化。

从这两组数据可以清晰看出，这是一个电化学材料的DFT计算项目——客户给了未掺杂、N掺杂和富缺陷三种二维材料样品，想从理论端解释为什么富缺陷样品在10 mA cm⁻²电流密度下的过电位最低。极化曲线是实验端的宏观表现，PDOS是计算端的微观解释，两者对照才能讲清楚活性差异的电子结构来源。

项目计算流程上，我用VASP做了几何优化、静态自洽和态密度计算。表面模型用2×2超胞，真空层设为18 Å，K点采6×6×1，交换关联泛函用PBE。HER自由能用ΔG_H* = ΔE_ZPE − TΔS近似，ZPE从频率计算得到。缺陷模型我没有直接去掉一个原子就结束，而是先做了位形熵估算，再比较了空位、棱位和边缘位的形成能，选形成能最低且配位数变化最大的构型作为富缺陷代表。

踩坑经验有两条。第一个是溶剂化效应：早期我用真空层模型算出的η₁₀只有130 mV，加了隐式溶剂模型VASPsol后，由于界面电场重排，过电位被修正到190 mV，与实验趋势更接近。第二个是赝势选择：该体系含过渡金属，我一开始用标准PAW赝势没加DFT+U，d带位置偏浅约0.3 eV，导致ΔG_H*预测过强；加了U=3.5 eV后ε_d才落到−1.5 eV，与文献对照合理。

综合来看，这张图的核心价值在于把宏观极化曲线与微观PDOS对应起来，定量给出了缺陷工程对HER活性的提升机制。如果要做进一步诊断，我建议把极化曲线拆成Tafel斜率和交换电流密度两个指标，同时结合微观动力学模拟给出turnover frequency。关于电化学材料的DFT计算，我在电化学材料dft计算栏目里整理了更多关于HER、OER和CO₂RR的理论分析经验。