VASP | 理论计算常见问题解答-80

A：这是典型的对角化阶段程序崩溃，ZHEGV 是 LAPACK 里解广义本征值问题 𝐴𝑥=𝜆𝐵𝑥Ax=λBx的子程序，要求重叠矩阵 𝐵B 正定。

①检查POSCAR文件是否有原子过近；

②删除 WAVECAR/CHGCAR；

③提高计算精度；

④不要用旧的 NPAR；改用 KPAR+NCORE，或者都不设置。

Q2：朱老师，如果在钙钛矿MDABCOKI₃将钾离子替换成NH₄⁺，是直接在MS中将K换成N，然后加H，导出POSCAR进行后续计算吗？这样做正确吗？

A：一、不用固定，不改电荷；

二、NH_₄^⁺在中心是可转动的，一个取向可能卡局域极小。建议至少试 3–6 种起始取向；

三、若能找到 [MDABCO]NH_₄I_₃ 的实验 CIF（或文献给的坐标），直接据此建模是最省心的。

Q3：朱老师，我对铈金属单质催化剂优化，得到的结构原子位置没有发生大的变化，但是吸附了我的物质后，催化剂原子位置有较大变化，导致吸附能有8点多，我还计算了铈掺杂别的元素，原子位置也有变化，那现在铈单质催化剂怎么去优化呢？

A：可以少放开一些表面原子，但不能全部固定；

看看ce的potcar没有4f；

得到 8 eV 级别的吸附能其实包含了表面重构，部分氧化的形成能，不是分子在金属上的弱吸附能。。

Q4：朱老师，对于自由离子而言，Fe³⁺的d带中心应该是多少啊？

A：自由离子 Fe^³⁺没有“d 带中心”这个量;

“d-band center/ d 带中心”是针对固体/表面里连续的 d 带（有参考的费米能级 ​）定义的；

自由离子只有离散 3d 轨道，不存在带也没有费米能级​，因此这个描述子在物理上不适用，也没有唯一数值。

Q5：朱老师，下面三幅图能用vasp算吗？

A：可以做，具体是：

用 VASP 得到电荷密度/静电势文件 → 后处理在分子表面上着色（ESP）→ 对 ESP 做球谐多极展开即可。

VASP本身不画图、不做球谐分析，但配合 VASPKIT / Multiwfn / VMD 就能复现图 b、c。

【做计算 找华算】

🏅 华算科技提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。

🎯500+博士团队护航，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏

为确保大家的学习成果得以沉淀和复用，我们特将朱老师本次及往期的VASP理论计算答疑精华，系统汇编为【理论计算答疑精选汇编】。它将成为你科研路上随取随用的“速查手册”，助你高效攻克难题。