计算化学公社

标题: 计算含金属大体系GEO_OPT不收敛 [打印本页]

作者
Author: ma455173220 时间: 2024-3-12 07:02
标题: 计算含金属大体系GEO_OPT不收敛

请问老师使用CP2K计算大体系时结构优化GEO_OPT不收敛，反而四个收敛参数以及结构能量越来越大一般有什么解决办法呢？SCF收敛没有问题，用的是OT算法。结构优化使用的BFGS算法，并且即使把TRUST_RADIUS设到0.05也依旧是一样的问题。输入文件请见附件中的cp2k-OT.inp。
体系中含金属，所以理论上来讲应该使用SMEAR，但是使用SCF对角化算法时一直提示非正交，不知会不会造成结果的不准确？输入文件和输出的WARNING请见附件中的cp2k-DIA.inp和cp2k-DIA.png。
当使用OT算法时，LBFGS算法倒是能收敛，但是感觉结构有点儿不靠谱，毕竟四个收敛参数基本都是NO，CG算法也测试过，同样不收敛，明不是很复杂的体系，只是掺杂一个原子而已。

作者
Author: sobereva 时间: 2024-3-12 12:27
这种情况必须用对角化+smearing

我重新创建了输入文件，用CP2K 2024.1没出现那个warning

(, 下载次数 Times of downloads: 50)

(, 下载次数 Times of downloads: 14)

作者
Author: ma455173220 时间: 2024-3-12 13:51
感谢老师的回复，我用超算自带的CP2K发现好像也没有WARNING，感觉有可能是我自己编译的ELPA的问题......哎...不过结构优化不收敛倒是依旧如常...

作者
Author: ma455173220 时间: 2024-3-12 18:01
问题解决了，如果之后有碰到同样问题的人可以一试：解决方法很简单，把TRUST_RADIUS 设成0.02就可以收敛了......我是真没想到0.05居然都不够小......

作者
Author: lywbanner 时间: 2024-10-12 12:58
感谢，探索者！

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