计算化学公社

 找回密码 Forget password
 注册 Register

vasp算自由能修正,在频率计算这一步为什么不收敛?

查看数: 13269 | 评论数: 6 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2020-4-11 18:29

正文摘要:

本帖最后由 将军金甲夜不脱 于 2020-4-11 18:39 编辑 vasp算自由能修正,在频率计算这一步,输入grep reach OUTCAR ,不会出现'reached required accuracy - stopping structural energy minimisation'字样,但是 ...

回复 Reply

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

老师您好,请教一个问题。
计算过程中虽然不牵扯极小化过程,但是如果计算时电子步不收敛,最终还是会影响我频率的计算的结果。是吗?
另外,我用vasp在计算过渡态频率时遇到第一个电子步不收敛该如何解决呢?是我跑的过渡态就是有问题的吗?
表面反应,原子数78,频率仅计算反应涉及原子的频率,过渡态计算和频率计算的收敛标准都是:能量1E-6,力0.02。
卡开发发 发表于 Post on 2020-4-12 19:00:18
将军金甲夜不脱 发表于 2020-4-12 18:01
谢谢大神,你不仅把我问的问题解决了,还给我普及了声子谱的知识,真是听君一席话胜读十年书啊。
...

优化和算声子谱(振动频率)应当采用同一计算精度。另外虚频来源要么是结构优化没做好(包括精度太低引起的情况),要么是边界取得太小。
卡开发发 发表于 Post on 2020-4-12 10:06:18
将军金甲夜不脱 发表于 2020-4-11 22:24
我没算过声子谱,这个就是单纯结构优化后算个频率。我看网上的教程说NSW>0就行哇大神

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

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

将军金甲夜不脱 发表于 Post on 2020-4-11 22:24:23
本帖最后由 将军金甲夜不脱 于 2020-4-11 22:25 编辑
卡开发发 发表于 2020-4-11 19:25
声子谱的情况NSW=1(参考手册6.61.2的例子),要不然NSW=500算500次?

我没算过声子谱,这个就是单纯结构优化后算个频率。我看网上的教程说NSW>0就行哇大神
卡开发发 发表于 Post on 2020-4-11 19:25:54
声子谱的情况NSW=1(参考手册6.61.2的例子),要不然NSW=500算500次?

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 14:47 , Processed in 0.222359 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list