“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 664|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] cp2k FPMD SCF不收敛求助

[复制链接 Copy URL]

261

帖子

0

威望

1467

eV
积分
1728

Level 5 (御坂)

实验组内的DFT计算、第一性原理、MD模拟爱好者

我cp2k中使用了gfn1-xtb搭配d3校正做涉及断键的FPMD模拟过程,使用OT技术,步长0.2fs,在3000K\1bar下跑NPT20ps,体系只有C/H/O/N四种元素,约580个原子。过程中出现scf不收敛,提示“SCF run NOT converged. To continue the calculation regardless, please set the keyword IGNORE_CONVERGENCE_FAILURE”。我看了最后一轮SCF收敛情况,发现震荡。我是直接启用IGNORE_CONVERGENCE_FAILURE继续运行,还是修改OT算法和OT极小化方法呢

最后一轮SCF.png (63.83 KB, 下载次数 Times of downloads: 20)

最后一轮SCF.png

polyamide_hydration_FPMD_xtb_3000K.inp

34.98 KB, 下载次数 Times of downloads: 4

目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

1295

帖子

6

威望

2756

eV
积分
4171

Level 6 (一方通行)

傻傻的木瓜

2#
发表于 Post on 2025-3-15 20:19:39 | 只看该作者 Only view this author
CP2K里GFN1-xTB做涉及变胞的任务可能有问题,之前http://bbs.keinsci.com/thread-51341-1-1.html有讨论,注意观察轨迹是否有类似那帖说的收缩问题。可以考虑改用DFT跑AIMD,比如常规的PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH组合;如果计算资源有限算不动如此大体系,或许可以缩小点。
√546=23.36664289109

261

帖子

0

威望

1467

eV
积分
1728

Level 5 (御坂)

实验组内的DFT计算、第一性原理、MD模拟爱好者

3#
 楼主 Author| 发表于 Post on 2025-3-16 15:02:33 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-3-15 20:19
CP2K里GFN1-xTB做涉及变胞的任务可能有问题,之前http://bbs.keinsci.com/thread-51341-1-1.html有讨论,注 ...

好的,感谢。我也正有此意
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-4-13 10:05 , Processed in 1.119827 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list