关于用OT结构优化之后，算单点能用对角化读取OT收敛波函数Convergence 越来越大的问题

ddddnight

老师，由于需要结构的态密度，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

sobereva

对角化时开了smearing，和OT时没有可比性，不属于相同计算条件

ddddnight

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

老师，不用smearing还是一样

ddddnight

之后我增大了cutoff还是不行，这种还有其他办法吗，老师

ball2006

OT下增大cutoff重新计算一下，然后再以此得到wfn文件对角化进行计算试试

ddddnight

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

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

ball2006

算单点能的时候增大cutoff，至于增加到多少得看具体研究的体系。我看你模拟的体系含有Na，cutoff可以取到1000试试

ddddnight

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

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

ball2006

是的

ddddnight

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

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

ball2006

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

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

ddddnight

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

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

ball2006

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

建议重新优化一下

ddddnight

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

好的，我试试，感谢