VASP中1×1×1优化结构做4×4×4 DOS时残余力很大怎么办？

RYL

各位老师好，想请教一个 VASP 结构优化与 DOS 后处理中的 k 点设置问题。
我有一个体相体系，结构优化时采用：
PBEGamma-centered 1×1×1 KPOINTSEDIFFG = -0.03 eV/ÅISIF = 3IBRION = 2结构优化正常结束，OUTCAR 中显示 reached required accuracy，最终最大力已经小于 0.03 eV/Å。因此，在 1×1×1 这个结构优化设置下，结构满足了设定的力收敛标准。
之后我用优化得到的 CONTCAR 作为固定结构，进行 DOS 计算。DOS 计算中没有再更新离子位置，而是将 k 点加密为：
Gamma-centered 4×4×4IBRION = -1NSW = 0ICHARG = 11DOS 计算电子步可以正常收敛，但 OUTCAR 中打印的 Hellmann-Feynman 力明显变大，最大力约为 2 eV/Å。
我想请教的是：

• 在这种情况下，4×4×4 非自洽 DOS 结果是否仍可以作为“固定 1×1×1 优化结构上的电子态密度”来分析？
• DOS 静态计算中较大的残余力，是否意味着必须先用 4×4×4 重新优化结构，再计算 DOS？
• 更准确的表述是否应为：该结构在 1×1×1 优化设置下是收敛的，但尚未证明其在 4×4×4 设置下也处于力平衡状态？
  

我的理解是，DOS 非自洽计算本身不进行离子优化，主要判断电子步和态密度采样是否合理；但 OUTCAR 中打印的力也反映了该结构在更密 k 点设置下的受力情况。因此不确定这种 DOS 是否适合用于正式讨论，还是必须重新在更密 k 点下优化结构。
请各位老师指教，谢谢。

RYL

现在发现用111做结构优化，设定力的收敛标准为-0.03ev，结构优化收敛以后 用CONTCAR和111k点密度做非自洽DOS计算 第一步跑完以后输出的最大受力还是0.4几 这种情况到底是怎么回事啊

阮文升-河内国大

你算 DOS 的时候用的是 ICHARG = 11？这节约计算量但是会导致精度下降。
（1）根据我个人经验，我觉得你遇到的问题可以不用管，当然严谨起见最好做 k 点收敛测试。但建议你给出更多关于体系性质的信息，帮助大家判断。
（2）如果算得动的话，建议用 ICHARG=1 或者读波函数代替。
（3）如果计算资源有限或者需要和之前的计算保持统一，用 ICHARG = 11 也可以，但是这时候波函数不是自洽的（电子密度不是来自当前波函数，而是固定为和一开始读的电子密度保持一致），此时得到的受力可能缺乏意义。