计算化学公社

标题: 关于晶体不同层数对几何优化收敛的影响 [打印本页]

作者
Author: THY 时间: 2025-6-19 22:28
标题: 关于晶体不同层数对几何优化收敛的影响
各位老师，本人使用CP2K对AlN晶体进行几何优化，发现随着AlN层数的增加，几何优化变得十分困难。

本人一共研究了4x4超胞下，单层，两层，三层，四层，四种结构的AlN，发现当层数小于等于两层，十分容易收敛，但是大于两层以后难以收敛，每轮SCF迭代次数基本都超过128了，

本人考虑过k点，diag，OT，smearing，固定层数，使用GC优化，等方法，并按照sob老师博文CP2K中遇到SCF难收敛时的解决方法 - 第一性原理 (First Principle) - 计算化学公社里面的说法都尝试进行解决

但是最终都无法解决，并且发现结构收敛与否好像只与层数有，但是本人想模拟晶体表面结构，所以层数需要在3~4层左右，请问该如何解决。

附上四种层数的cif文件，以及单层和三层的inp文件和out文件，单层的很容易收敛，但是三层的基本每轮SCF迭代次数基本都超过128，但是最终还是收敛了，不知道这种收敛靠不靠谱？

(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 3)

作者
Author: get-it 时间: 2025-6-20 21:06
是否做了结构弛豫？

作者
Author: THY 时间: 2025-6-21 15:54

get-it 发表于 2025-6-20 21:06
是否做了结构弛豫？

老师，除了单层没有做，2-3层结构的都是固定底层，开放表面一层

作者
Author: sobereva 时间: 2025-6-22 07:05
当前结果没有问题
几何优化初期就是容易SCF难收敛，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）这页幻灯片也强调了：
(, 下载次数 Times of downloads: 18)

作者
Author: THY 时间: 2025-6-22 14:49

sobereva 发表于 2025-6-22 07:05
当前结果没有问题
几何优化初期就是容易SCF难收敛，北京科音CP2K第一性原理计算培训班（http://www.keinsc ...

sob老师，您好，像这种每轮SCF迭代次数都超过128次，导致最后out文件有很多个不收敛的Warning，但最后四个收敛标准都达到YES了，这种结果能接受吗？（例如1楼中的"3layer-AlN.out"）还是说应该在inp文件里面上调整迭代次数？

作者
Author: get-it 时间: 2025-6-22 23:36
本帖最后由 get-it 于 2025-6-23 00:03 编辑

THY 发表于 2025-6-21 15:54
老师，除了单层没有做，2-3层结构的都是固定底层，开放表面一层

有文章计算过，纤锌矿AlN的层数比较少时，类似于石墨的纯平结构会更稳定。所以你固定底层很可能就没法让结构弛豫收敛。建议放开所有原子的限制，也放开晶格限制，作cell optimization.

作者
Author: sobereva 时间: 2025-6-23 04:59

THY 发表于 2025-6-22 14:49
sob老师，您好，像这种每轮SCF迭代次数都超过128次，导致最后out文件有很多个不收敛的Warning，但最后四 ...

优化到后期，优化的每一步的SCF最终都收敛了，没什么问题

作者
Author: THY 时间: 2025-6-23 09:27

get-it 发表于 2025-6-22 23:36
有文章计算过，纤锌矿AlN的层数比较少时，类似于石墨的纯平结构会更稳定。所以你固定底层很可能就没法让 ...

谢谢老师，我试一下

作者
Author: THY 时间: 2025-6-23 09:27

sobereva 发表于 2025-6-23 04:59
优化到后期，优化的每一步的SCF最终都收敛了，没什么问题

谢谢老师

作者
Author: Jiangqin 时间: 2025-7-1 11:25

THY 发表于 2025-6-22 14:49
sob老师，您好，像这种每轮SCF迭代次数都超过128次，导致最后out文件有很多个不收敛的Warning，但最后四 ...

你好，请问在.out文件中最后的哪里查看四个收敛标准是YES还是NO哦

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