计算化学公社

标题: 关于用OT结构优化之后，算单点能用对角化读取OT收敛波函数Convergence 越来越大的问题 [打印本页]

作者
Author: ddddnight 时间: 2024-3-11 23:06
标题: 关于用OT结构优化之后，算单点能用对角化读取OT收敛波函数Convergence 越来越大的问题
老师，由于需要结构的态密度，cp2k计算晶体表面单点能，OT收敛之后，对角化读取OT收敛的.wfn文件，为什么这个convergence这么大啊，还越来越大，我把NBROYDEN增加到20，ALPHA减小到0.2还是这样，是不是结构有问题。相关输入输出文件已经上传。
OT收敛的.wfn文件：链接：https://pan.baidu.com/s/1rVZdP3qcM__6PeuIToXM3A?pwd=jd16 提取码：jd16
--来自百度网盘超级会员V5的分享
SCF WAVEFUNCTION OPTIMIZATION

  Step    Update method    Time Convergence       Total energy Change
  ------------------------------------------------------------------------------
   1 NoMix/Diag. 0.20E+00 14.5    0.08900804 -32985.3128026740 -3.30E+04
   2 Broy./Diag. 0.20E+00 16.8    0.26219850 -32982.8964797727  2.42E+00
   3 Broy./Diag. 0.20E+00 16.7    0.47372682 -33005.7831703441 -2.29E+01
   4 Broy./Diag. 0.20E+00 15.7    0.93622278 -32962.3350946594  4.34E+01
   5 Broy./Diag. 0.20E+00 16.9    0.97599664 -33060.9102296500 -9.86E+01
   6 Broy./Diag. 0.20E+00 16.7    1.04081235 -32949.1200677070  1.12E+02
   7 Broy./Diag. 0.20E+00 16.6    1.06575081 -33081.6233511831 -1.33E+02
   8 Broy./Diag. 0.20E+00 16.9    1.09276795 -32945.3188796000  1.36E+02
   9 Broy./Diag. 0.20E+00 17.3    1.45037019 -33089.4454640624 -1.44E+02
10 Broy./Diag. 0.20E+00 16.9    2.10990966 -32977.1729669421  1.12E+02
11 Broy./Diag. 0.20E+00 17.1    4.35576353 -33086.9334353722 -1.10E+02
12 Broy./Diag. 0.20E+00 16.9 15.40540776 -32560.8669316402  5.26E+02
13 Broy./Diag. 0.20E+00 17.2 372.25764144 -32220.8023891155  3.40E+02
14 Broy./Diag. 0.20E+00 17.6 375.05245040 -30275.8035381751  1.94E+03
15 Broy./Diag. 0.20E+00 17.0 981.48718729 -30402.5034684129 -1.27E+02
16 Broy./Diag. 0.20E+00 17.6 981.48704745 -30002.8649452983  4.00E+02
17 Broy./Diag. 0.20E+00 17.5 939.63033302 -30571.1371907507 -5.68E+02
18 Broy./Diag. 0.20E+00 17.3 939.63068430 -29961.1095845925  6.10E+02
19 Broy./Diag. 0.20E+00 17.8 802.54788221 -30434.3854612933 -4.73E+02
20 Broy./Diag. 0.20E+00 17.6 802.54751832 -29998.2290568626  4.36E+02
21 Broy./Diag. 0.20E+00 18.1 262.69819343 -30262.9669030778 -2.65E+02
22 Broy./Diag. 0.20E+00 18.0 262.69991927 -30868.1922164356 -6.05E+02
23 Broy./Diag. 0.20E+00 17.7 712.53990112 -29972.9445784307  8.95E+02
24 Broy./Diag. 0.20E+00 18.0 712.53699649 -29968.7708891473  4.17E+00
25 Broy./Diag. 0.20E+00 18.0 249.69885143 -30435.8471216474 -4.67E+02
26 Broy./Diag. 0.20E+00 17.6 249.69945758 -29991.7404388678  4.44E+02
27 Broy./Diag. 0.20E+00 18.4 275.33962100 -30509.9334460865 -5.18E+02
28 Broy./Diag. 0.20E+00 17.4 275.33888122 -29719.4765644519  7.90E+02
29 Broy./Diag. 0.20E+00 17.9 981.59055941 -30860.6673017609 -1.14E+03
30 Broy./Diag. 0.20E+00 18.0 981.59305124 -29761.4534647838  1.10E+03

作者
Author: sobereva 时间: 2024-3-12 10:43
对角化时开了smearing，和OT时没有可比性，不属于相同计算条件

作者
Author: ddddnight 时间: 2024-3-12 13:12

sobereva 发表于 2024-3-12 10:43
对角化时开了smearing，和OT时没有可比性，不属于相同计算条件

老师，不用smearing还是一样

作者
Author: ddddnight 时间: 2024-3-15 22:49
之后我增大了cutoff还是不行，这种还有其他办法吗，老师

作者
Author: ball2006 时间: 2024-3-16 20:11
OT下增大cutoff重新计算一下，然后再以此得到wfn文件对角化进行计算试试

作者
Author: ddddnight 时间: 2024-3-17 20:28

ball2006 发表于 2024-3-16 20:11
OT下增大cutoff重新计算一下，然后再以此得到wfn文件对角化进行计算试试

感谢您的解答，请问OT下增大cutoff是结构优化还是说算一下能量，cutoff大概增加到多少呀，OT重新算的时候需要读取之前结构优化的wfn文件吗？

作者
Author: ball2006 时间: 2024-3-19 19:29
算单点能的时候增大cutoff，至于增加到多少得看具体研究的体系。我看你模拟的体系含有Na，cutoff可以取到1000试试

作者
Author: ddddnight 时间: 2024-3-19 20:02

ball2006 发表于 2024-3-19 19:29
算单点能的时候增大cutoff，至于增加到多少得看具体研究的体系。我看你模拟的体系含有Na，cutoff可以取到10 ...

我是这样的，计算过程中Na用全电子基组，您的意思是OT算单点增加cutoff，然后diag算读取OT的.wfn是吧

作者
Author: ball2006 时间: 2024-3-20 22:58
是的

作者
Author: ddddnight 时间: 2024-3-22 17:25

ball2006 发表于 2024-3-20 22:58
是的

很遗憾啊，还是不行，同样一个体系，La收敛了，Ce就不行

作者
Author: ball2006 时间: 2024-3-31 21:45

ddddnight 发表于 2024-3-22 17:25
很遗憾啊，还是不行，同样一个体系，La收敛了，Ce就不行

Ce的4f有单电子，要不要考虑一下自旋

作者
Author: ddddnight 时间: 2024-4-1 22:47

ball2006 发表于 2024-3-31 21:45
Ce的4f有单电子，要不要考虑一下自旋

考虑自旋的话是不是要重新优化

作者
Author: ball2006 时间: 2024-4-2 19:19

ddddnight 发表于 2024-4-1 22:47
考虑自旋的话是不是要重新优化

建议重新优化一下

作者
Author: ddddnight 时间: 2024-4-2 19:32

ball2006 发表于 2024-4-2 19:19
建议重新优化一下

好的，我试试，感谢

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