计算化学公社

标题: 优化NiFe2O4晶胞，cp2k几何优化不收敛，但SCF收敛 [打印本页]

作者
Author: yiranfengbai 时间: 2024-5-14 00:09
标题: 优化NiFe2O4晶胞，cp2k几何优化不收敛，但SCF收敛
老师您好，我用CP2K优化NiFe2O4晶胞，泛函PBE，使用了2X2X2超胞，磁矩设置为：Fe(tet): -4，Fe(oct): 4，Ni(oct): 2，MULTIPLICITY设置为：1+2*8*8=129 {参考论文中e(tet): -3.98，Fe(oct): 4.14，Ni(oct): 1.68 }。U值中Fe: 3.3（相似结构CoFe2O4文献中Fe: 3.3），Ni设为2（参考文献中为6，但测试发现U值超过4会出错，所以准备测试2~4）。之前计算设置SCF的收敛限1E-6时，运行了24h发现，从第30步开始能量上升，SCF不收敛。把SCF的收敛限设为1E-5时，发现每一步的SCF收敛，但几何优化步数达到400步时却不收敛。
请问老师，inp文件的设置是否有错误呢？该怎么修改呢？

作者
Author: sobereva 时间: 2024-5-14 07:00
步步为营，不变胞优化成功再说变胞的事、不+U成功再说+U的事。而且文献里的Ueff参数普遍是对平面波程序优化的，根本没法直接给CP2K用，严重偏大

几何优化不收敛常规解决方式在北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里我都讲过

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

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

作者
Author: ljh123 时间: 7 day ago

sobereva 发表于 2024-5-14 07:00
步步为营，不变胞优化成功再说变胞的事、不+U成功再说+U的事。而且文献里的Ueff参数普遍是对平面波程序优化 ...

sob老师，我是打算对一个金属络合物和水的溶液跑AIMD。跑AIMD需要对该提体系作预结构优化。请问这种溶液体系，需要考虑k点和cutoff吗？

作者
Author: sobereva 时间: 6 day ago

ljh123 发表于 2024-11-16 19:02
sob老师，我是打算对一个金属络合物和水的溶液跑AIMD。跑AIMD需要对该提体系作预结构优化。请问这种溶液 ...

算什么都得恰当设定CUTOFF
是否需要考虑k点完全看当前盒子尺寸，必须把关键信息交代出来

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