Q1：朱老师，不同金属的化学势要怎么求？
A：核心是化学势等于原子库的能量，要先考虑清楚原子从哪来/到哪去；
金属/非气体元素：稳定本征晶相的每原子 DFT 能量，0K近似;
气体元素：取单体分子的能量除以原子个数；
总结就是：取单晶/气体中一个原子的能量
Q2：朱老师，同一个体系，一定是电子数越多能量越低吗？
A：不一定
带电超胞由 VASP 用“均匀背景电荷”中和。这样得到的总能包含人造自相互作用/偶极耦合，随超胞尺寸、真空厚度而变，能量可上可下，并不具有可比性；
加电子可能填到反键（能量↑），也可能填到成键/空导带（能量↓）。有没有降低，取决于你填的是哪类态以及结构是否随之重构；
常用 GGA 存在分数电子的凸性误差：电子会被平均化，看起来加一点电子就降低能量，那是数值偏差，不是物理趋势。
Q3：朱老师，计算DOS出现报错？
A：继续用四面体法，但把 k 网格加密到 ≥4 个不可约 k 点；
保持 Γ-only，不要用四面体法，改用小展宽的高斯/费米展宽；
注意金属体系做 DOS 不建议 Γ-only。
Q4：朱老师，能带的布里渊区两边是不是必须要对称才行，两边都是G点的话，对应的能带位置是不是也要一样？
A：不必须对称；但相同的 k 点上，能量必须一样；
画的是一条路径 Γ−𝑋−𝑆−𝑌−Γ。左右两端都是 Γ（𝑘=0）。在同一次非自洽 band 计算里，Γ点的本征能 完全相同；
整张图不需要以中线镜像对称，因为左端邻近段是G-X，右端是Y-G，方向不同，斜率/曲线形状不必相同。
Q5：朱老师，对于计算ORR吉布斯自由能有点疑问，想请教下。下图计算HOOH最后一步应属于脱附计算，但多数文献的HOOH能量台阶最终归向一个值，想问问这里能量是指？
A：图末端的统一值代表产物 H₂O₂（溶液）+ 干净表面相对初态的自由能；
之所以相同，是因为该终态与位点无关；
能量指生成H2O2的标准自由能变，ΔG=-1.40eV.
为确保大家的学习成果得以沉淀和复用，我们特将朱老师本次及往期的VASP理论计算答疑精华，系统汇编为【理论计算答疑精选汇编】。它将成为你科研路上随取随用的“速查手册”，助你高效攻克难题。
【做计算 找华算】
🏅 华算科技提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。 🎯500+博士团队护航，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏
