计算化学公社

标题: CP2K计算MnFe2O4的(111)表面不收敛 [打印本页]

作者
Author: zhaoquanyou 时间: 2025-2-10 19:20
标题: CP2K计算MnFe2O4的(111)表面不收敛
我想用cp2k对MnFe2O4的(111)表面做几何优化和计算能量，后边算对H2O2的吸附能和差分等数据，结果从开头就卡住了，我使用卢老师提供的方法比如说设置Mn和Fe为5的初猜磁矩和设置BS，对MnFe2O4晶体能做计算能收敛，几何优化和能量计算都没有问题，切面之后用这两种方法都试了试，结果都是不收敛的，我想的原因是表面Mn和Fe的化学环境发生了变化，导致不是+2和+3价了，可是我又不知道表面的他们是什么化学价态。我还有一个问题，最近学了一些cp2k和VASP的计算，包括之前做的CASTEP的计算发现，每种算出来的能量都是不一样的，cp2k算出来的特别特别大，这个是正常的吗？

作者
Author: sobereva 时间: 2025-2-10 20:31
对表面原子尝试不同的自旋磁矩的设置进行计算。如果发现QE等程序更好收敛，也可以拿QE算出来的表面原子自旋磁矩信息当CP2K计算的初猜。或者搜索计算过此体系的相关文献参考其结果进行初猜的设置

不同程序算的绝对能量没有可比性

别乱找地方发帖，CP2K的问题必须发到第一性原理版块的CP2K分类里。这次给你移动了，下次删帖扣分处理

不要24小时内在思想家公社QQ群和计算化学公社论坛里问同一问题，注意看群规（http://sobereva.com/QQrule.html），以及论坛置顶的论坛新社员必读贴。下次发现删帖扣分处理

作者
Author: zhaoquanyou 时间: 2025-2-10 20:46

sobereva 发表于 2025-2-10 20:31
对表面原子尝试不同的自旋磁矩的设置进行计算。如果发现QE等程序更好收敛，也可以拿QE算出来的表面原子自旋 ...

好的，谢谢老师

作者
Author: zhaoquanyou 时间: 2025-2-10 22:29

sobereva 发表于 2025-2-10 20:31
对表面原子尝试不同的自旋磁矩的设置进行计算。如果发现QE等程序更好收敛，也可以拿QE算出来的表面原子自旋 ...

老师，我还有一个问题，我之前拿VASP和CASTEP计算的计算磁矩可以作为参考吗，还是说最好用QE啊

作者
Author: sobereva 时间: 2025-2-10 23:13

zhaoquanyou 发表于 2025-2-10 22:29
老师，我还有一个问题，我之前拿VASP和CASTEP计算的计算磁矩可以作为参考吗，还是说最好用QE啊

可以

作者
Author: zhaoquanyou 时间: 2025-3-4 10:58
老师，我按照CASTEP算出来的这个表面磁矩，Mn设置为4，Fe根据情况设置为-2和3，作为CP2K的初猜，发现还是不收敛，是我的设置有问题吗？

作者
Author: zhaoquanyou 时间: 2025-3-10 08:44

sobereva 发表于 2025-2-10 23:13
可以

老师，我用CASTEP算出来的磁矩是这样，Mn设置为4，Fe根据情况设置为-2和3，作为CP2K的初猜，发现还是不收敛，是我的设置有问题吗？

作者
Author: sobereva 时间: 2025-3-10 22:56

zhaoquanyou 发表于 2025-3-10 08:44
老师，我用CASTEP算出来的磁矩是这样，Mn设置为4，Fe根据情况设置为-2和3，作为CP2K的初猜，发现还是不 ...

贴图方式不对，仔细看置顶的社员必读贴了解怎么贴图

SCF不收敛问题下文该说的都说了
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

作者
Author: zhaoquanyou 时间: 2025-3-11 09:14

sobereva 发表于 2025-3-10 22:56
贴图方式不对，仔细看置顶的社员必读贴了解怎么贴图

SCF不收敛问题下文该说的都说了

好的好的，谢谢老师

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