计算化学公社

标题: CP2K计算晶体表面吸附孤立原子不收敛问题 [打印本页]

作者
Author: THY 时间: 2025-5-19 22:19
标题: CP2K计算晶体表面吸附孤立原子不收敛问题
各位老师，本人使用CP2K想要计算在ZnO晶体表面吸附了一个孤立的O原子和一个孤立的N原子，根据电子轨道分布，给O原子设置磁矩2μB，给N原子设置磁矩3μB，但是体系一直振荡，不收敛，请问该如何解决？

(, 下载次数 Times of downloads: 7) 附件为inp文件

作者
Author: sobereva 时间: 2025-5-20 06:22
该说的都说了
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）
当前晶胞尺寸不大，不考虑k点不合理

另外，孤立原子带的单电子状态和化学吸附后的状态并没有可比性。

作者
Author: THY 时间: 2025-5-20 08:25

sobereva 发表于 2025-5-20 06:22
该说的都说了
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/threa ...

老师，您的意思是, 不应该设置磁矩吗？

作者
Author: sobereva 时间: 2025-5-21 04:03

THY 发表于 2025-5-20 08:25
老师，您的意思是, 不应该设置磁矩吗？

没说不该设置，只是说盲目套用原子孤立状态的磁矩进行设置没意义。

作者
Author: THY 时间: 2025-5-21 10:36

sobereva 发表于 2025-5-21 04:03
没说不该设置，只是说盲目套用原子孤立状态的磁矩进行设置没意义。

谢谢老师，假如孤立状态的O原子磁矩为2μB，孤立状态N原子磁矩为3μB，那吸附态的O和N的磁矩应该是小于2μB和3μB吧？这种情况下，怎么在inp里面设置初始磁矩呢，是需要一点一点的试吗？

作者
Author: sobereva 时间: 2025-5-22 03:59

THY 发表于 2025-5-21 10:36
谢谢老师，假如孤立状态的O原子磁矩为2μB，孤立状态N原子磁矩为3μB，那吸附态的O和N的磁矩应该是小于2 ...

是

不需要特别讲究，可以试几个（设得不理想可能SCF难收敛，此时再试其它的），反正SCF过程自旋磁矩也会自发往实际合理方向调整

作者
Author: THY 时间: 2025-5-22 09:45

sobereva 发表于 2025-5-22 03:59
是

不需要特别讲究，可以试几个（设得不理想可能SCF难收敛，此时再试其它的），反正SCF过程自旋磁矩也 ...

谢谢老师

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