计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sunhao
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 用不同的ALGO参数优化出来的结构能量差别很大。

[复制链接 Copy URL]

38

帖子

0

威望

217

eV
积分
255

Level 3 能力者

16#
 楼主 Author| 发表于 Post on 2019-12-4 09:30:42 | 只看该作者 Only view this author
卡开发发 发表于 2019-12-4 08:50
1、后续如果吸附的体系比较大的话可以适当放宽一点。
2、做收敛性测试的目的之一就是保证能量和力计算可 ...

如图,我在大师兄教程上看到的。他说的意思是只要不是变胞优化,就可以统一使用正常的相对较低的ENCUT值。所以我在想是不是可以用大的ENCUT值优化好的晶胞参数固定住,再用小的ENCUT值优化一下,后面的吸附统一用小的ENCUT值去计算,这样能量有可比性。但不知道这种想法对不对。。主要是ENCUT设置成800计算起来还是很吃力的。。。后续吸附计算都用800比较浪费时间。希望老师能够指导指导。

1575422609(1).jpg (55.13 KB, 下载次数 Times of downloads: 65)

1575422609(1).jpg

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

17#
发表于 Post on 2019-12-4 09:54:21 | 只看该作者 Only view this author
sunhao 发表于 2019-12-4 09:30
如图,我在大师兄教程上看到的。他说的意思是只要不是变胞优化,就可以统一使用正常的相对较低的ENCUT值 ...

无论是晶格常数计算还是一般的构型优化,stress和force对动能截断都是有依赖的,测试的目的就是为了获得可靠合理的参数,高或低应当有个定量的标准。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

38

帖子

0

威望

217

eV
积分
255

Level 3 能力者

18#
 楼主 Author| 发表于 Post on 2019-12-4 10:11:30 | 只看该作者 Only view this author
卡开发发 发表于 2019-12-4 09:54
无论是晶格常数计算还是一般的构型优化,stress和force对动能截断都是有依赖的,测试的目的就是为了获得 ...

好的,我明白了。非常感谢老师的指导!

38

帖子

0

威望

217

eV
积分
255

Level 3 能力者

19#
 楼主 Author| 发表于 Post on 2019-12-11 10:44:03 | 只看该作者 Only view this author
卡开发发 发表于 2019-12-4 09:54
无论是晶格常数计算还是一般的构型优化,stress和force对动能截断都是有依赖的,测试的目的就是为了获得 ...

老师您好,我先用ISIF=3优化好了晶格,然后用ISIF=2固定住晶格之后在优化结构。用ISIF=2优化好的CONTCAR复制成POSCAR,想检验一下是不是优化好了。但是发现不是一个离子步结束,而且能量差别依然很大,磁矩也差别很大。请问老师这可能是什么原因呢?我的体系比较复杂,我设置的力的收敛标准是-0.05.

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

20#
发表于 Post on 2019-12-11 11:19:46 | 只看该作者 Only view this author
sunhao 发表于 2019-12-11 10:44
老师您好,我先用ISIF=3优化好了晶格,然后用ISIF=2固定住晶格之后在优化结构。用ISIF=2优化好的CONTCAR ...

1、一直用ISIF=3就行,没必要切换到ISIF=2。
2、看看是否有可能磁矩猜错。你可以LORBIT=11把磁矩print出来,然后按照那个趋势去做初猜,看看能量是否有较大差异。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

38

帖子

0

威望

217

eV
积分
255

Level 3 能力者

21#
 楼主 Author| 发表于 Post on 2019-12-11 11:31:37 | 只看该作者 Only view this author
卡开发发 发表于 2019-12-11 11:19
1、一直用ISIF=3就行,没必要切换到ISIF=2。
2、看看是否有可能磁矩猜错。你可以LORBIT=11把磁矩print出 ...

好的老师,我试试。谢谢您

38

帖子

0

威望

217

eV
积分
255

Level 3 能力者

22#
 楼主 Author| 发表于 Post on 2019-12-11 11:43:10 | 只看该作者 Only view this author
卡开发发 发表于 2019-12-11 11:19
1、一直用ISIF=3就行,没必要切换到ISIF=2。
2、看看是否有可能磁矩猜错。你可以LORBIT=11把磁矩print出 ...

老师您能加我一下好友吗?我想单独请教你一下。我的微信号:wxid_cjt40nkp3h9712

62

帖子

0

威望

744

eV
积分
806

Level 4 (黑子)

23#
发表于 Post on 2022-1-4 17:26:31 | 只看该作者 Only view this author
卡开发发 发表于 2019-11-30 03:25
1、调整ALGO没有用处啊,ALGO是控制SCF中对角化的方法的。很难收敛的情况RMM-DIIS确实可能不稳定,相对而言 ...

老师,结构优化的时候是不是一般选择ALGO=Normal比较好?
Life was like a box of chocolates. You never know what you're going to get.

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

24#
发表于 Post on 2022-1-4 23:36:57 | 只看该作者 Only view this author
Adair 发表于 2022-1-4 17:26
老师,结构优化的时候是不是一般选择ALGO=Normal比较好?

是,一般建议ALGO=N,体系很大的话,如果数值还比较稳定,ALGO=V或F也行。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 13:28 , Processed in 0.173415 second(s), 22 queries , Gzip On.

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