VASP问题答疑系列

Q1：朱老师，a方向晶格常数不失配，但b方向失配了，还要继续调b吗？ A：①需要，直到都匹配； ②做supercell配对； ③施加应变； ④换取向或换切面方向 Q2：朱老师，我粘…

Q1：朱老师，算了100步没有收敛，这个算是成功了么？ A：电子自洽（SCF）在达到设定的最大电子步数 NELM=100 时仍没满足 EDIFF 收敛条件； 这一离子步的电子结构没…

Q1：朱老师，β型的NiMoO4的cif怎么找呀？在materials project上只有这两个? A：MP给的两个NiMoO4，严格说都不是多数电化学文献里常说的高温β-NiM…

Q1：朱老师，LDH材料OER计算，很多文献只给了DFT+U的U值和最后一步的POSCAR，但是Ni-Fe的mag数值没有明确给出，我想做OER计算是不是INCAR当中我只写ISP…

Q1：朱老师，请问我做AIMD的体系跑MOF材料，在NVT下设置的初始和末态温度为100K，0.5fs的步长，跑200步，OSZICAR中的T就从100开始掉到50左右，请问我应该…

Q1：朱老师，不同金属的化学势要怎么求？ A：核心是化学势等于原子库的能量，要先考虑清楚原子从哪来/到哪去； 金属/非气体元素：稳定本征晶相的每原子 DFT 能量，0K近似; 气体…

Q1：朱老师，有时候算出来的态密度图，-10eV左边，最高的这个峰有用吗？ A：①这是内层的s电子的贡献； ②多数情况下没用，可以直接裁掉； ③什么时候有意义？赝势/价电子选择检查…

Q1：朱老师，想看圈起来那三个原子的差分电荷二维图，应该切哪个面呢？ A：只能先算整体3D差分，再截面获取2D差分图； Edit → Plane → New → By 3 poin…

Q1：朱老师，这个是什么问题，poscar的吗？ A：这是典型的对角化阶段程序崩溃，ZHEGV 是 LAPACK 里解广义本征值问题 ??=???Ax=λBx的子程序，要求重叠矩阵…

Q1：朱老师，怎么算氢吸附自由能？ A：没有问题，这样可以直接拿去优化； 标准流程：基态结构→构建 H*→几何优化→能量处理→自由能修正→考虑电位与 pH; 研究“水解离/Volm…

Q1：朱老师，这样哪里有问题需要改？ A：这是 VASP 在 电子步数(NELM)用完仍未自洽 时的通用告警。几何、力、能量此步都不可靠。 ①从干净的初始起步：删除 WAVECAR…

Q1：朱老师，请问在计算吸附能的过程中是不是IVDW这个参数也一直保留，在仅计算金属催化剂表面的能量时发现加IVDW和不加有显著差异？ A：要么全程都开，要么全程都关；绝不要混用；…

Q1：朱老师，neb计算时，这个nebef.pl之后第二列全是0？ A：一、可能是没有编译好VTST； 二、在计算过渡态之前需要检查原子是否一一对应； 三、第二个这个截图是正常的。…

Q1：朱老师，计算FE110表面H2的吸附时，需要考虑色散校正开启IVDW=12吗？ A：需要开启，H₂在 Fe (110) 表面的吸附通常属于化学吸附，而非纯粹物理吸附； 研究分…

Q1：朱老师，我算出来的声子谱是这样的是什么原因，有非常明显的虚频？ A：若前期结构优化不彻底，原子会存在 自发位移的趋势，声子谱会直接将这种趋势转化为明显虚频，具体表现为： 一、…

Q1：朱老师，这是我做的二维钙钛矿材料，沿（0 0 1）切面后加了个真空层。主要是想后面做缺陷，再吸附分子做缺陷钝化，这个切面步骤是有必要的吧？ A：①切面是聚焦表面 / 界面研究…

Q1：朱老师，想研究一个二维材料和一个衬底slab模型异质结界面，结构优化是用isif=3还是2？ A：优先选ISIF=3，无原子固定用ISIF=3，并且固定真空层方向不做晶胞优化…

Q1：朱老师，这个报错：“WARNING: Sub-Space-Matrix is not hermitian in DAV”怎么解决？ A：金属表面试试ISMEAR=2 SIGM…

Q1：朱老师，请问我这个过渡态计算最有可能是过渡态的是哪个点啊 A：如果结构没问题，能量最高点是过渡态 Q2：朱老师，我是要按照这里显示的去下载对应版本的python吗？但我下载好…

  Q1：朱老师，我同一个模型切同一面后，在算功函数时，如何确定那个暴露面最好呢？是通过什么数值来确定的呢？ A：根据实验结果哪个是暴露面 Q2：朱老师，这里不同缺陷原子对net …