计算化学公社

标题: cp2k计算结构优化时总是卡死，卡在 Extrapolation method: initial_guess，但无报错 [打印本页]

作者
Author: 汪杰 时间: 2025-12-23 19:18
标题: cp2k计算结构优化时总是卡死，卡在 Extrapolation method: initial_guess，但无报错
如题，在计算一个结构优化的任务时，莫名的卡在 Extrapolation method: initial_guess
inp文件与计算日志文件log在附件中 (, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 3)

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-12-23 19:58
这个FeB和C3N4的异质结倾斜而细长，不好观察结构是否合理，和我在http://bbs.keinsci.com/thread-54317-1-1.html五楼后两个链接里的求助帖是类似问题。可以上传从数据库下载的两种晶体的原始结构cif，并确认一下是否用FeB的这个晶面形成异质结。

作者
Author: aoj195 时间: 2025-12-23 20:05
这个问题感觉像软件bug，我之前在做变胞优化的时候也遇到过这个问题，运行固定步数后就会卡在initial_guess，只能不停的运行restart文件直到任务结束

作者
Author: 汪杰 时间: 2025-12-23 20:07

Uus/pMeC6H4-/キ发表于 2025-12-23 19:58
这个FeB和C3N4的异质结倾斜而细长，不好观察结构是否合理，和我在http://bbs.keinsci.com/thread-54317-1-1 ...

感谢回答，该异质结使用了vaspkit造出来的，失配度在1%以内。
这个异质结也是在尝试是否可行。从目前结果来看的话，可能时结构不行导致优化难以为继？

作者
Author: 汪杰 时间: 2025-12-23 20:14

aoj195 发表于 2025-12-23 20:05
这个问题感觉像软件bug，我之前在做变胞优化的时候也遇到过这个问题，运行固定步数后就会卡在initial_guess ...

cp2k很考验技术，这个体系具有磁性也是一个扰动点。
CG算法就会出现卡住的问题，BFGS倒是不会出现卡住，但是很难收敛，十几步一直在震荡。
要想收敛还得做许多尝试

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-12-23 20:34

汪杰发表于 2025-12-23 20:07
感谢回答，该异质结使用了vaspkit造出来的，失配度在1%以内。
这个异质结也是在尝试是否可行。从目前结 ...

我链接里面的相关内容你有看到吗？异质结建模只看失配度还不够，对一个晶格矢量夹角偏小的周期性盒子计算既不方便前后处理观察，也可能从算法层面上给积分格点的划分与运算带来未知的数值问题（这点不针对任何特定第一性原理程序），所以我让上传原始结构就是想看能不能改用更好的异质结结构（既包括调整夹角，也可能可以减小原子数量，但是只有尝试了才知道）。这里g-C3N4的结构还好说，主要是不知道这FeB是哪个数据库的哪种晶型，以及用了什么晶面。

作者
Author: 汪杰 时间: 2025-12-23 21:22

Uus/pMeC6H4-/キ发表于 2025-12-23 20:34
我链接里面的相关内容你有看到吗？异质结建模只看失配度还不够，对一个晶格矢量夹角偏小的周期性盒子计算 ...

明白了，FeB那个我也是受人请求，他给的结构。不知道出处。
此外，vaspkit构建异质结时，已经遍历了所有可能的结构。我从中挑选了一个原子数量可以接受的
夹角也稍微能接受的一个构型。其他夹角更大的，原子数要到大几百去了。

作者
Author: sobereva 时间: 2025-12-26 04:39
变胞优化卡住这事本身是非常老生常谈的问题，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里专门讲了：
(, 下载次数 Times of downloads: 0)

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