计算化学公社

标题: 请教cp2k结构优化任务的续算方法 [打印本页]

作者
Author: 白菜 时间: 2025-8-23 18:57
标题: 请教cp2k结构优化任务的续算方法
老师们好~
因为结构优化跑了四十轮，SCF都没有收敛的趋势，我手动kill结构优化任务，当前路径有：xxx-1.restart、xxx.inp、xxx.out、xxx-pos-1.pdb、xxx-BFGS.Hessian、xxx-RESTART.kp
在检索论坛和网络后发现提及了3种续算方式：
A法. 用xxx-1.restart文件直接做为输入文件，用cp2k.popt -i xxx-1.restart续算
B法. 将xxx-1.restart文件重命名为xxx.inp，用cp2k.popt -i xxx.inp续算
C法. 从xxx-1.restart读取信息，取消注释Multiwfn创建的xxx.inp输入文件中的以下内容，然后用cp2k.popt -i xxx.inp续算：
   &SCF/WFN_RESTART_FILE_NAME xxx-RESTART.kp
   &DFT/SCF_GUESS RESTART
   &MOTION/&CELL_OPT&BFGS/RESTART_HESSIAN T
   &MOTION/&CELL_OPT&BFGS/RESTART_FILE_NAME xxx-BFGS.Hessian
1. 老师们，请问这3种方法续算方法完全相同吗？
2. 如果我想按照CP2K中遇到SCF难收敛时的解决方法 - 思想家公社的门口：量子化学·分子模拟·二次元中各种方法尝试让SCF收敛，可以直接在续算的输入文件里改关键词吗？会不会有什么问题使得计算续不上呢？
谢谢~

作者
Author: ultramanm87 时间: 2025-8-24 01:49
1. 都行
2. 先尽量收敛，有时候一开始的结构尤其是slab模型不好收敛不少见，但是你现在四十轮优化没有收敛应该先解决SCF收敛的事情

作者
Author: sobereva 时间: 2025-8-24 03:05
1 A和B完全一样，建议用B，免得文件名混淆。C原理上更好一些，由于读取了之前的Hessian和波函数，能使得续算时候更好地承接上一步的状态

2 可以。能否兼容取决于具体做法，比如改成OT时不兼容k点因此没法用

作者
Author: 白菜 时间: 2025-8-24 06:24

ultramanm87 发表于 2025-8-24 01:49
1. 都行
2. 先尽量收敛，有时候一开始的结构尤其是slab模型不好收敛不少见，但是你现在四十轮优化没有收敛 ...

嗯嗯好哒，谢谢ultramanm87老师~╰(*°▽°*)╯

作者
Author: 白菜 时间: 2025-8-24 06:25
嗯嗯好滴，谢谢sob老师~

作者
Author: zhubaonian 时间: 2025-11-26 17:02
你好，请问续算后怎么让输出文件.out也接上中断前的数据，因为需要完整的.out文件进行后续的deepmd数据集构建。

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