计算化学公社

标题: vasp算自由能修正，在频率计算这一步为什么不收敛？ [打印本页]

作者
Author: 将军金甲夜不脱 时间: 2020-4-11 18:29
标题: vasp算自由能修正，在频率计算这一步为什么不收敛？
本帖最后由将军金甲夜不脱于 2020-4-11 18:39 编辑

vasp算自由能修正，在频率计算这一步，输入grep reach OUTCAR ,不会出现'reached required accuracy - stopping structural energy minimisation'字样，但是在前一步结构优化时都轻松收敛了。
比如，我算了乙醇的频率计算，发现无法收敛，再看离子步，正好算了55个离子步。在查看频率，27个频率也都算出来了。下面是我的结果图和INCAR,KPOINTS.现在就是很纳闷，到底我这个频率算出来了吗
？我又算了水分子，发现也是这种情况，有没有大神指导啊

作者
Author: 卡开发发 时间: 2020-4-11 19:25
声子谱的情况NSW=1（参考手册6.61.2的例子），要不然NSW=500算500次？

作者
Author: 将军金甲夜不脱 时间: 2020-4-11 22:24
本帖最后由将军金甲夜不脱于 2020-4-11 22:25 编辑

卡开发发发表于 2020-4-11 19:25
声子谱的情况NSW=1（参考手册6.61.2的例子），要不然NSW=500算500次？

我没算过声子谱，这个就是单纯结构优化后算个频率。我看网上的教程说NSW>0就行哇大神

作者
Author: 卡开发发 时间: 2020-4-12 10:06

将军金甲夜不脱发表于 2020-4-11 22:24
我没算过声子谱，这个就是单纯结构优化后算个频率。我看网上的教程说NSW>0就行哇大神

1、这两个概念等同，可以把晶格振动假想成一个粒子，不同振动能级就对应着这些粒子的能级，这个粒子就是声子。
2、NSW好像是>0就行，确实记错了。做声子谱（频率）计算的时候，如果是IBRION=5，则是计算结构的单点+6N个有限移动，比如乙醇是9个原子，则总离子步数应该是1+9*6=55，55步完成就自然结束了，得到每一个结构的力在用于差分近似来得到动力学矩阵。计算过程不牵扯到极小化这样的迭代步骤，因此不会提示stopping structural energy minimisation。

作者
Author: 将军金甲夜不脱 时间: 2020-4-12 18:01
本帖最后由将军金甲夜不脱于 2020-4-12 18:57 编辑

卡开发发发表于 2020-4-12 10:06
1、这两个概念等同，可以把晶格振动假想成一个粒子，不同振动能级就对应着这些粒子的能级，这个粒子就是 ...

谢谢大神，你不仅把我问的问题解决了，还给我普及了声子谱的知识，真是听君一席话胜读十年书啊

。
我还有2个问题，1：就是我算水分子的频率时，虚频特别大（见下图），INCAR还是上面乙醇那个。我看有人说虚频不大于100都可以忽略，但是这个也太大了吧。。。
2：算乙醇的频率用不用加vdw=11这一项呢？

作者
Author: 卡开发发 时间: 2020-4-12 19:00

将军金甲夜不脱发表于 2020-4-12 18:01
谢谢大神，你不仅把我问的问题解决了，还给我普及了声子谱的知识，真是听君一席话胜读十年书啊。
...

优化和算声子谱（振动频率）应当采用同一计算精度。另外虚频来源要么是结构优化没做好（包括精度太低引起的情况），要么是边界取得太小。

作者
Author: BADEGG 时间: 2023-6-8 22:57

卡开发发发表于 2020-4-12 10:06
1、这两个概念等同，可以把晶格振动假想成一个粒子，不同振动能级就对应着这些粒子的能级，这个粒子就是 ...

老师您好，请教一个问题。
计算过程中虽然不牵扯极小化过程，但是如果计算时电子步不收敛，最终还是会影响我频率的计算的结果。是吗？
另外，我用vasp在计算过渡态频率时遇到第一个电子步不收敛该如何解决呢？是我跑的过渡态就是有问题的吗？
表面反应，原子数78，频率仅计算反应涉及原子的频率，过渡态计算和频率计算的收敛标准都是：能量1E-6，力0.02。

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