关于使用双杂化泛函计算C6Li，FeN4在第一步PBE.inp阶段不收敛的问题

牛顿树上摘苹果 · 发表于 Post on 2025-5-10 21:42:07

我在使用双杂化泛函计算LiC6222超胞的能量是，按照cp2k培训班计算尿素的案例进行计算，但是在第一步PBE.inp阶段就遇到了不收敛的问题，计算FeN4体系也遇到了同样的问题，求助各位老师们。

sobereva · 发表于 Post on 2025-5-10 22:44:49

LiC6222、FeN4的gap非常小，本来就不适合双杂化。而且PBE计算也应该开smearing，而且当前晶胞也不够大，不考虑k点说不过去。
双杂化不是什么情况都适合的

牛顿树上摘苹果 · 发表于 Post on 2025-5-11 21:47:48

本帖最后由牛顿树上摘苹果于 2025-5-11 21:49 编辑

感谢sob老师的指导
1.按照您给的建议，我加了smear和d3修正，给LiC6体系扩了444的胞（由于我想看一下双杂化泛函的结果，双杂化泛函又不支持k点，所以没有考虑k点），使用CC-TZ基组仍然遇到了不收敛的问题，之后我换成DZVP-MOLOPT-SR-GTH成功收敛了。之后我使用这步得到的PBE-RESTART.wfn文件，结合CC-TZVP按照DHDF.inp进行双杂化泛函计算是否合理？
2.但是FeN4的体系设置了上述之后仍然无法收敛，我考虑是否是我磁性设置的不对的问题。.inp文件已上传于附件。
3.另外我想请教一下sob老师，在双杂化泛函计算时，第一步PBE.inp的计算是否没有必要设置辅助基组。然后对于像双杂化泛函这样无法考虑k点的计算， -7 Set direction(s) of applying periodic boundary condition, current: XYZ是否应该设置为NONE

北大-陶豫 · 发表于 Post on 2025-5-11 22:18:00

gap非常小，本来就不适合双杂化

想看一下双杂化泛函的结果

双杂化泛函是掺了MP2关联能。你去看一下MP2关联能的公式，分母上是空轨道能级减去占据轨道能级。如果gap非常小，显然MP2关联能就爆了，这种情况不用“想看一下双杂化泛函的结果”，看一下公式就行了。

牛顿树上摘苹果 · 发表于 Post on 2025-5-12 08:34:43

北大-陶豫发表于 2025-5-11 22:18
双杂化泛函是掺了MP2关联能。你去看一下MP2关联能的公式，分母上是空轨道能级减去占据轨道能级。如果ga ...

谢谢老师指导，我后面再找一些适合做双杂化泛函的体系，我想问一下老师对于我计算FeN4，PBE不收敛老师有没有解决的建议？

北大-陶豫 · 发表于 Post on 2025-5-12 10:00:08

牛顿树上摘苹果发表于 2025-5-12 08:34
谢谢老师指导，我后面再找一些适合做双杂化泛函的体系，我想问一下老师对于我计算FeN4，PBE不收敛老师有 ...

CP2K中遇到SCF难收敛时的解决方法
http://bbs.keinsci.com/forum.php ... 6&fromuid=19420
(出处: 计算化学公社)
gap很小的体系难收敛，可以理解为SCF过程中电子空穴会到处乱窜。这时候我会首先选择加大smearing。

牛顿树上摘苹果 · 发表于 Post on 2025-5-12 21:07:28

老师，如果我使用111k点可以认为FeN4是有带隙的吗？

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

[CP2K] 关于使用双杂化泛函计算C6Li，FeN4在第一步PBE.inp阶段不收敛的问题