优化NiFe2O4晶胞，cp2k几何优化不收敛，但SCF收敛

yiranfengbai · 发表于 Post on 2024-5-14 00:09:43

老师您好，我用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文件的设置是否有错误呢？该怎么修改呢？

sobereva · 发表于 Post on 2024-5-14 07:00:54

步步为营，不变胞优化成功再说变胞的事、不+U成功再说+U的事。而且文献里的Ueff参数普遍是对平面波程序优化的，根本没法直接给CP2K用，严重偏大

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

ljh123 · 发表于 Post on 2024-11-16 19:02:43

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

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

sobereva · 发表于 Post on 2024-11-17 02:04:58

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 优化NiFe2O4晶胞，cp2k几何优化不收敛，但SCF收敛