CP2K几何优化Fe2O3（001）表面不收敛

Stardust0831 · 发表于 Post on 2025-10-16 01:21:31

喝酸奶的巧乐兹发表于 2025-10-15 18:23
谢谢老师指导！
您所说的“对Fe2O3的(001)表面来说还不太够只用gamma点就能算准的程度吧”，
请问表面 ...

其实我觉得15A一般就够用了，绝大部分体系在25A内能收敛。
想知道具体多少合适的话，可以做收敛性测试。

Stardust0831 · 发表于 Post on 2025-10-16 01:24:49

Uus/pMeC6H4-/キ发表于 2025-10-15 17:09
哦，你是说“小分子和基底统一用gamma点”，那和“小分子和基底统一用更大k点”一样，也符合我说的“含k ...

其实这个问题换一个可以同时考虑k点并支持DFT+U的程序，就不需要纠结这么多了。不知道QE和ABACUS这些免费程序里有没有满足当前需求的。

现在cp2k的困境就是不开DFT+U会scf难收敛，开DFT+U的时候，听楼主的意思似乎没法考虑k点

喝酸奶的巧乐兹 · 发表于 Post on 2025-10-16 10:07:45

Stardust0831 发表于 2025-10-16 01:24
其实这个问题换一个可以同时考虑k点并支持DFT+U的程序，就不需要纠结这么多了。不知道QE和ABACUS这些免费 ...

谢谢您的解答！QE这些平面波软件计算比cp2k还是慢多了，但现在看来确实换软件比较方便了

Stardust0831 · 发表于 Post on 2025-10-16 10:24:33

喝酸奶的巧乐兹发表于 2025-10-16 10:07
谢谢您的解答！QE这些平面波软件计算比cp2k还是慢多了，但现在看来确实换软件比较方便了

QE会慢一些，但是磁性体系很稳健。
ABACUS可不慢，以Si体系为例：

ABACUS通过使用和cp2k相同的赝势，以及由此赝势生成的轨道文件进行EOS计算。在其他参数尽可能对齐的情况下，ABACUS的平面波方法（pw）计算精度最高，但耗时最长；ABACUS的轨道方法（lcao）计算精度略低于平面波（pw）方法但优于CP2K，同时耗时约为cp2k的60%，在精度和效率上均优于cp2k。

详见：ABACUS对比CP2K精度和效率测试 | Si的状态方程（EOS）

喝酸奶的巧乐兹 · 发表于 Post on 2025-10-31 23:52:35

后续：换用DFT+U和OT完美收敛，吸附多种小分子后也很容易收敛
如果不需要计算轨道信息，使用OT确实是个好方法

Assiduitas · 发表于 Post on 2025-11-1 09:39:36

喝酸奶的巧乐兹发表于 2025-10-31 23:52
后续：换用DFT+U和OT完美收敛，吸附多种小分子后也很容易收敛
如果不需要计算轨道信息，使用OT确实是个好 ...

请问你使用DFT+U，这个U的参数值是怎么选择的？具体的数值来源是哪？

喝酸奶的巧乐兹 · 发表于 Post on 2025-11-5 09:54:49

Assiduitas 发表于 2025-11-1 09:39
请问你使用DFT+U，这个U的参数值是怎么选择的？具体的数值来源是哪？

参考了一点文献的u值，3左右，然后用这个u值进行了几何结构优化，得到的晶胞参数与实验值非常吻合，就这么确定的
不过我也不知道这是不是最合理的做法

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

[CP2K] CP2K几何优化Fe2O3（001）表面不收敛

浏览过的版块