求助，用CP2K计算YBCO超导体系1100K下的nvt，每一步SCF的电子步迭代要好几十次

bentty_lee

如题，有没有什么方法，可以像别的体系一样，每轮SCF迭代个不超过10次就收敛了。一轮几十次，这计算量也太大了。

我查看了sober老师的CP2K中遇到SCF难收敛时的解决方法 - 第一性原理 (First-principles) - 计算化学公社这个帖子，感觉最主要的可能的改善方式是改磁性设置。

但是YBaCuO7是超导体，YBaCuO6有反铁磁性，我这个体系是用SQS掺杂了O6.0625、6.1875、3125等几个浓度吧，很难确定Cu2O层和CuO层具体的电子组态。我个人觉得由于1100K是很高的温度了，可以不考虑磁性，所以我的INPUT里面也没考虑磁性。

这个材料本身没有带隙，所以也不能用OT。基组我是参考文献里的选的（如图）。

为了计算的快一点，我已经把截断能设置的比较小了，smearing也用上了。不知道还有没有什么其他方法能改善我现在的参数啊，我想大概算一万步离子步的话，时间控制在一周算完。

Assiduitas

OT本身就是非常有利于AIMD计算。材料没有带隙不能使用OT这个句话是否正确，需要斟酌一下。

bentty_lee

Assiduitas 发表于 2025-9-15 16:29
OT本身就是非常有利于AIMD计算。材料没有带隙不能使用OT这个句话是否正确，需要斟酌一下。

直接上OT很难收敛，金属算不了空轨道。我用对角化的波函数当初猜再算OT倒是很快收敛了，就是不知道这样会不会有问题。