计算化学公社

标题: 求助CP2K计算两层Au（111）收敛问题 [打印本页]

作者
Author: zhouchen 时间: 2025-4-1 14:01
标题: 求助CP2K计算两层Au（111）收敛问题
大家好，我计算了两层Au（111）,在XY方向具有周期性。利用PBE0计算，加入DFT-D3，使用diag对角化，开启smearing。但是收敛趋势显示，一开始收敛正常，但是后面收敛值就开始增大。请问这是为什么？是我的两层Au（111）建模有问题还是计算设置有问题？

作者
Author: Stardust0831 时间: 2025-4-1 15:27
本帖最后由 Stardust0831 于 2025-4-1 15:28 编辑

是否读取了PBE下收敛的波函数当初猜？
参考：《CP2K做杂化泛函计算的关键要点和简单例子》（http://sobereva.com/690）

当前金属的厚度也有些薄了。

作者
Author: zhouchen 时间: 2025-4-1 15:35

Stardust0831 发表于 2025-4-1 15:27
是否读取了PBE下收敛的波函数当初猜？
参考：《CP2K做杂化泛函计算的关键要点和简单例子》（http://sobere ...

您好，我已经使用PBE的收敛波函数作为初猜了。
请问金属的厚度也会影响体系的收敛吗？我当时想的是本来就不容易收敛，所以尽可能的减少金原子数

作者
Author: sobereva 时间: 2025-4-1 22:12
两层根本什么都不是，就算能算完，审稿人一看这结构也会质疑。一般至少用4层原子
而且当前用PBE0的意义不明。

作者
Author: zhouchen 时间: 2025-4-2 10:42

sobereva 发表于 2025-4-1 22:12
两层根本什么都不是，就算能算完，审稿人一看这结构也会质疑。一般至少用4层原子
而且当前用PBE0的意义不 ...

老师，相同的设置情况下，我用了2，3层Au(111)大概200原子左右，都不收敛（收敛振荡），用了4层之后，原子数更多，288个左右，这次反而能收敛，而且收敛趋势非常明显，请问这是因为什么？

作者
Author: sobereva 时间: 2025-4-2 14:40

zhouchen 发表于 2025-4-2 10:42
老师，相同的设置情况下，我用了2，3层Au(111)大概200原子左右，都不收敛（收敛振荡），用了4层之后，原 ...

层数太少的情况边界效应、有限尺寸效应太明显，电子结构非常诡异、偏离真实的界面的电子结构状态都很远

作者
Author: zhouchen 时间: 2025-4-7 21:49

sobereva 发表于 2025-4-2 14:40
层数太少的情况边界效应、有限尺寸效应太明显，电子结构非常诡异、偏离真实的界面的电子结构状态都很远

老师，同样的四层金原子，我在上面吸附了一个很小的分子，利用cp2k的PBE进行结构优化之后，再用PBE0去算，还是收敛振荡。

作者
Author: sobereva 时间: 2025-4-7 23:50

zhouchen 发表于 2025-4-7 21:49
老师，同样的四层金原子，我在上面吸附了一个很小的分子，利用cp2k的PBE进行结构优化之后，再用PBE0去算 ...

CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

用PBE0算这个没有额外优势

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