求助：Al2O3计算激发能时候不收敛

周志博

在对Al2O3进行电子激发的计算时候，按照靛青(indigo)的激发态计算例子进行设置，计算的时候发现在第一步自洽场计算的时候就无法收敛，converge达到几百。其中理论方法设为PBE0，基组设为6-31G*，计算的态数设为20，周期性设置为NONE。
我尝试了使用multiwfn中得其他全电子基组，皆无法收敛，converge甚至有达到10的7次方，理论方法有换过B3LYP，也同样不收敛，还尝试过cutoff提升到800，打开色散矫正，开启smearing，皆无济于事。想请问为什么会如此？
注：Al2O3.out文件是6-31G*基组的计算结果，Al2O3(1).out是converge达到上万时候的结果。

sobereva

你算的是团簇还是晶体？算小团簇根本没必要用CP2K

并且把下文仔细认真看了，免得做出用SCREEN_ON_INITIAL_P T却不读像样波函数当初猜的典型低级错误
CP2K做杂化泛函计算的关键要点和简单例子
http://sobereva.com/690（http://bbs.keinsci.com/thread-41087-1-1.html）

周志博

sobereva 发表于 2024-5-15 22:40
你算的是团簇还是晶体？算小团簇根本没必要用CP2K

并且把接下来仔细认真研究，免得做出用SCREEN_ON_INITIA ...

老师，我计算的任务是晶体，只有十个原子。
我读了您的这篇文章，CP2K做杂化泛函计算的关键要点和简单例子http://sobereva.com/690（http://sobereva.com/690）/bbs.keinsci.com/thread-41087-1-1.html）
觉得我的问题可能是初猜结构的不合理，然后尝试先使用泛函PBE进行初步计算，以得到比较合理的初猜结构（我的基组选择是6-31G*），然后再进行杂化泛函PBE0的计算，但是使用PBE计算时仍然无法收敛。我继续尝试几个能使scf收敛的方法，仍然无果。最好的仅收敛到0.3左右，然后便开始在此范围内震荡。

关于这个报错
*** WARNING in pw/ps_wavelet_methods.F:238 :: Density non-zero on the ***
*** edges of the unit cell: wrong results in WAVELET solver           ***
这个报错是由于初步猜测结构不合理导致的报错么？我在网上搜索过，有人回答是在inp文件加入字节&TOPOLOGY&CENTER_COORDINATES，不过似乎并无作用，multiwfn产生的inp文件中已经有这字节

sobereva

周志博 发表于 2024-5-17 17:47
老师，我计算的任务是晶体，只有十个原子。
我读了您的这篇文章，CP2K做杂化泛函计算的关键要点和简单例 ...

晶胞不合理，当前根本就不是周期性体系
而且都知道是晶体，周期性哪能用NONE，此时Multiwfn给你默认用的PSOLVER完全不适合当前情况

没周期性体系计算常识的话最好系统性学学再做，要不然胡算瞎算完全是白耽误工夫，就算计算没报错结果也可能根本没意义。北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）

周志博

sobereva 发表于 2024-5-18 03:21
晶胞不合理，当前根本就不是周期性体系
而且都知道是晶体，周期性哪能用NONE，此时Multiwfn给你默认用的 ...

明白了，谢谢老师

周志博

sobereva 发表于 2024-5-18 03:21
晶胞不合理，当前根本就不是周期性体系
而且都知道是晶体，周期性哪能用NONE，此时Multiwfn给你默认用的 ...

老师，我使用xyz周期性之后，可以正常收敛。然后我在计算激发态的时候，先用PBE计算了一次，正常收敛，然后用优化后的restart结构做PBE0的计算，SCF却不收敛。
emm，两者同为PBE泛函，做PBE0泛函的时候不收敛，属于正常现象么？这种情况还需要用杂化泛函计算其激发态么？

sobereva

周志博 发表于 2024-5-23 22:07
老师，我使用xyz周期性之后，可以正常收敛。然后我在计算激发态的时候，先用PBE计算了一次，正常收敛，然 ...

把此文仔细看了，免得犯低级错误
CP2K做杂化泛函计算的关键要点和简单例子
http://sobereva.com/690（http://bbs.keinsci.com/thread-41087-1-1.html）