计算化学公社

标题: 求助vasp优化金属氧化物吸附模型不能收敛 [打印本页]

作者
Author: 卡卡西 时间: 2023-3-18 21:34
标题: 求助vasp优化金属氧化物吸附模型不能收敛
本帖最后由卡卡西于 2023-4-15 16:54 编辑

各位老师，我对CoAl2O4金属氧化物模型切面优化，之后计算优化模型对糠醛分子的吸附，计算了两周一直不收敛

，尝试过各种吸附构型，离子步都难以收敛，现在不确定是模型问题还是INCAR参数设置问题，求指教。（因为体系是三斜晶系，在优化切面时，按照VASPKIT给出的K-POINTS 1*2*1进行计算，系统一直报错，扩大到2*4*2依旧不行，后发现调整为1*1*1不再报错，也完成了优化，所以对于吸附体系K-POINTS也取1*1*1，而且晶胞比较大，有187个原子，也无法再增大K-POINTS，不知道和这个有没有关系？）INCAR和吸附模型如下图所示

作者
Author: 卡卡西 时间: 2023-3-18 21:36
对了，由于计算耗时过长，我曾尝试将INCAR中的IBRION改为1，但是并没有改善

作者
Author: 含光君 时间: 2023-3-19 00:18
1. 分清SCF收敛和结构优化收敛。69次收敛的是SCF，这69步加起来只是结构优化中的一步。原子受力离收敛还早。
2. 对于目前这个体系的原子数与k点，使用的核数太少，算得慢是自然的事。可加上NCORE=20关键词。或者降低k点精度/使用更多核心数。改IBRION和电子步耗时没有关系。
3. 合理。

作者
Author: 卡卡西 时间: 2023-3-19 13:44

含光君发表于 2023-3-19 00:18
1. 分清SCF收敛和结构优化收敛。69次收敛的是SCF，这69步加起来只是结构优化中的一步。原子受力离收敛还早 ...

嗯嗯，谢谢老师

作者
Author: 卡卡西 时间: 2023-4-15 13:42
刚更新了问题，目前吸附模型一直不收敛，也没找到原因，求各位老师指教

作者
Author: Deepast 时间: 2023-4-15 16:12
没输出文件看不懂
ediff至少1e-5，不然受力算不准很难收敛
吸附计算最好加d3
开始优化的时候可以先用ibrion=1加速收敛，跑一段时间后换成2
这种体系最好做好计算速度和收敛性测试，特别是scf和并行参数，用vasp_gam跑。正常的话prec=normal就够了。

作者
Author: 卡卡西 时间: 2023-4-15 16:48
本帖最后由卡卡西于 2023-4-15 16:55 编辑

Deepast 发表于 2023-4-15 16:12
没输出文件看不懂
ediff至少1e-5，不然受力算不准很难收敛
吸附计算最好加d3

嗯嗯，谢谢老师，您的话我再研究一下，有些东西不是很懂。ediff和ediffg是想先分别用1e-4和-0.05，等收敛后再调整为1e-5和-0.02，这样可以提高计算速度，这是从网上学到的一个方法，另外我又上传了没有收敛的outcar文件，麻烦您看看再指导一下，感激不尽

作者
Author: Deepast 时间: 2023-4-15 17:11

卡卡西发表于 2023-4-15 16:48
嗯嗯，谢谢老师，您的话我再研究一下，有些东西不是很懂。ediff和ediffg是想先分别用1e-4和-0.05，等收敛 ...

才40多步离子步，不收敛很正常。你这个体系构型不好跑个两三百步，甚至500步以上也不奇怪

作者
Author: 卡卡西 时间: 2023-4-15 17:18

Deepast 发表于 2023-4-15 17:11
才40多步离子步，不收敛很正常。你这个体系构型不好跑个两三百步，甚至500步以上也不奇怪

哦哦，不过这只是其中的一个模型，我看有的模型算了100步力越来越来越不收敛了（如下面的数据所示），以为是我的设置有问题。
1 4.61392
2 9.27844
3 1.49262
4 1.32597
5 0.858911
6 1.03015
7 0.739568
8 1.20341
9 0.690736
10 1.09181
11 0.743849
12 1.0785
13 0.589998
14 0.924558
15 0.606664
16 1.43953
17 1.10657
18 0.75559
19 0.462345
20 1.42316
21 1.15758
22 0.733403
23 1.72404
24 0.836611
25 0.734007
26 1.22099
27 1.0866
28 1.05094
29 1.07985
30 0.699213
31 1.43029
32 1.08542
33 1.25734
34 2.65432
35 1.18112
36 1.07343
37 1.07084
38 1.30015
39 1.22891
40 1.29947
41 1.89894
42 1.38701
43 1.42366
44 1.49066
45 3.85597
46 1.59936
47 1.37259
48 1.2188
49 1.12481
50 0.720393
51 0.972002
52 0.865274
53 1.54973
54 1.24909
55 0.80393
56 0.833437
57 1.49504
58 1.19848
59 1.03934
60 1.43378
61 1.31055
62 1.2431
63 2.0176
64 1.37032
65 0.994056
66 1.08977
67 0.792219
68 1.17692
69 0.437346
70 0.498278
71 0.30776
72 0.58411
73 0.40984
74 1.11952
75 0.870625
76 0.644503
77 1.31966
78 0.777327
79 0.3832
80 0.71844
81 0.517036
82 0.697519
83 0.527567
84 0.882983
85 0.496279
86 0.685212
87 2.07539
88 1.15243
89 0.773913
90 0.773377
91 1.39333
92 0.884774
93 0.912574
94 0.99271
95 0.953947
96 0.85841
97 0.746021
98 1.1194
99 0.779905

作者
Author: Deepast 时间: 2023-4-15 17:31
不太好说，这要看结构有没有发生较大的改变，比如说糠醛分子优化过程中是不是转了个角度什么的
在我看来这都还好，主要是初猜问题
另外这种体系一定要有耐心，先跑个几百步再说

作者
Author: Deepast 时间: 2023-4-15 17:33
还有12核实在太小了，不建议算这么大的体系

作者
Author: 卡卡西 时间: 2023-4-15 18:57

Deepast 发表于 2023-4-15 17:33
还有12核实在太小了，不建议算这么大的体系

好的好的，谢谢老师

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