计算化学公社

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

[CP2K] cp2k FPMD SCF不收敛求助

[复制链接 Copy URL]

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的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: 0)

最后一轮SCF.png

polyamide_hydration_FPMD_xtb_3000K.inp

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

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

890

帖子

3

威望

1684

eV
积分
2634

Level 5 (御坂)

傻傻的木瓜

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

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的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, 2025-8-15 09:35 , Processed in 0.183403 second(s), 23 queries , Gzip On.

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