计算化学公社

标题: 求助VASP计算O2在N掺杂碳/纯碳材料上吸附的模型结果与文献差异较大，应该如何优化？ [打印本页]

作者
Author: zrf12138 时间: 2023-11-21 19:42
标题: 求助VASP计算O2在N掺杂碳/纯碳材料上吸附的模型结果与文献差异较大，应该如何优化？
大家好，小白在计算中有一些问题求教：
做了一个O2在N掺杂碳/纯碳材料上吸附的模型，使用VASP计算，参考文献大多都是Gaussian 09做的非周期性计算
参考文献中结合能大多为负，我使用周期性计算得到的结果均＞0，非周期性计算结果不收敛
这种情况应该怎么办呢，是计算有错误，还是软件不同导致的差异，很头疼。。

作者
Author: Stars 时间: 2023-11-21 23:48
结构不合理。将你的周期性做平移，与文献结构并不相同（比如氢在周期性体系中有的不应存在））

作者
Author: zrf12138 时间: 2023-11-22 09:28

Stars 发表于 2023-11-21 23:48
结构不合理。将你的周期性做平移，与文献结构并不相同（比如氢在周期性体系中有的不应存在））

嗯，我明白这个问题。但是现在用周期性模型去模拟文献结构对我来说有些困难，如果直接计算文献中的非周期结构又并不收敛，应该怎么处理呢，能给一些建议吗？谢谢了

作者
Author: Stars 时间: 2023-11-23 18:18

zrf12138 发表于 2023-11-22 09:28
嗯，我明白这个问题。但是现在用周期性模型去模拟文献结构对我来说有些困难，如果直接计算文献中的非周期 ...

看看你非周期性计算的模型和输入。

作者
Author: zrf12138 时间: 2023-11-24 00:01

Stars 发表于 2023-11-23 18:18
看看你非周期性计算的模型和输入。

这是POSCAR和INCAR

作者
Author: 白皛白 时间: 2023-11-24 15:27
主要的问题是INCAR中IVDW=0这一项，吸附模型是需要使用DFT-D方法的，建议使用DFT-D3。另外使用ISPIN=2的话，就需要把磁矩的初猜加上，否则默认所有原子的磁矩都是1，这显然是不合理的，容易导致SCF步骤不收敛的问题。

作者
Author: zrf12138 时间: 2023-11-24 16:41

白皛白发表于 2023-11-24 15:27
主要的问题是INCAR中IVDW=0这一项，吸附模型是需要使用DFT-D方法的，建议使用DFT-D3。另外使用ISPIN=2的话 ...

感谢解答，我再重新试试

作者
Author: yxzinc 时间: 2023-11-24 16:42
感觉你两个模型差别有点大，一个是有限的分子一个是无限的一维材料，电子结构大不相同。vasp也可以做近似做非周期计算，只要把三个方向的真空区域都设置的大一些，但是高斯和VASP所用的基组和泛函差别太大了，得到的结果可能完全无法比较

欢迎光临计算化学公社 (http://bbs.keinsci.com/)