计算化学公社

 找回密码 Forget password
 注册 Register
Views: 4141|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] vasp,结构优化和自洽计算结果不一致问题

[复制链接 Copy URL]

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位好,有时候做完结构优化后,用弛豫完的结构做自洽计算,得到的结果和结构优化最后的离子步结果不一致,差距在零点几~几eV,而有时又十分吻合,请问,这是什么原因造成的?无论是结构优化还是自洽计算,电子步都是收敛的。

689

帖子

2

威望

4101

eV
积分
4830

Level 6 (一方通行)

2#
发表于 Post on 2022-5-15 12:50:45 | 只看该作者 Only view this author
优化晶胞参数了吗? 如果晶胞体积有变化的话 两次计算的平面波数量会不完全一样 导致误差;当然这个问题不大

269

帖子

0

威望

2058

eV
积分
2327

Level 5 (御坂)

3#
发表于 Post on 2022-5-15 14:29:22 | 只看该作者 Only view this author
优化还有自洽的时候电子和离子步的收敛线是一致的么
努力做计算的小菜鸡

3

帖子

0

威望

166

eV
积分
169

Level 3 能力者

4#
发表于 Post on 2022-5-15 14:34:24 | 只看该作者 Only view this author
你所指的结果不一致是能量不一致吗?如果是的话,结构优化是离子步中嵌套电子步,当最后一个离子步收敛时,计算就停止了,认为达到收敛,但此时电子步并未收敛,所以需要静态自洽(单点能计算),单点计算时结构不变(NSW=0/1或 IBRION=-1),只是电子步优化,收敛后即最准确的能量。两者一般来说相差很小。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
乐平 + 5 我很赞同

查看全部评分 View all ratings

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-5-15 15:31:15 | 只看该作者 Only view this author
jiangning198511 发表于 2022-5-15 12:50
优化晶胞参数了吗? 如果晶胞体积有变化的话 两次计算的平面波数量会不完全一样 导致误差;当然这个问题不 ...

优化了,明白您的意思了。这样的话,在结构弛豫完后,最后再做一次自洽计算比较好吧

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2022-5-15 15:31:52 | 只看该作者 Only view this author
吴双 发表于 2022-5-15 14:29
优化还有自洽的时候电子和离子步的收敛线是一致的么

基本一致

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-5-15 15:34:23 | 只看该作者 Only view this author
Miles 发表于 2022-5-15 14:34
你所指的结果不一致是能量不一致吗?如果是的话,结构优化是离子步中嵌套电子步,当最后一个离子步收敛时, ...

明白您的意思了,我隐约有这种感觉,您这一解释我就清晰了。看来,每一次结构优化完毕后,再做一次自洽计算比较好

3

帖子

0

威望

166

eV
积分
169

Level 3 能力者

8#
发表于 Post on 2022-5-15 16:41:39 | 只看该作者 Only view this author
Shana 发表于 2022-5-15 15:34
明白您的意思了,我隐约有这种感觉,您这一解释我就清晰了。看来,每一次结构优化完毕后,再做一次自洽计 ...

本着严谨的态度,计算形成能、吸附能、分离能等能量时,都应该做单点计算得到电子步收敛时的能量。当然有的偷懒只用结构优化的能量来表示,差异也不明显

评分 Rate

参与人数
Participants 2
eV +8 收起 理由
Reason
晨晨晨嘛 + 3
乐平 + 5 我很赞同

查看全部评分 View all ratings

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2022-5-15 17:27:41 | 只看该作者 Only view this author
Miles 发表于 2022-5-15 16:41
本着严谨的态度,计算形成能、吸附能、分离能等能量时,都应该做单点计算得到电子步收敛时的能量。当然有 ...

好的,明白了

2

帖子

0

威望

37

eV
积分
39

Level 2 能力者

10#
发表于 Post on 4 day ago | 只看该作者 Only view this author
转个计算老司机关于该问题的帖子:https://www.bilibili.com/opus/497026122446620883

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 01:00 , Processed in 0.202400 second(s), 22 queries , Gzip On.

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