计算化学公社

标题: cp2k中用QMMM+CDFT的结构优化参数设置求助 [打印本页]

作者
Author: 高章磊 时间: 2025-11-13 22:11
标题: cp2k中用QMMM+CDFT的结构优化参数设置求助
我最近在CP2K里面进行使用QMMM模型的CDFT结构优化，主要目的是为了算出电荷转移前后的结构然后利用四点法计算重组能λ和自由能变δG，最终计算出kET。可是目前在进行结构优化时发现我的体系算了接近200步，花了超过一天的时间还没有完成。我怀疑是我输入文件里的Maximum gradient 和RMS gradient 设置太过严格导致的，但是我又不太敢确定调松能调多少，怕算出来结果没意义。
我的体系是MD跑出来的大的盒子，盒子本身有两万多原子，但中心分子如图二只有173个原子不是很多，然后就是以中心为QM其余为MM做的这个QMMM+CDFT的结构优化。
求助：
1. 虽然看图三确实各个项都在逐渐收敛但是截止目前还是没能收敛到目标值，我想问一下有经验的人这里的参数能放宽多少呢？（用于最终算kET）
2. 因为我要算这个盒子里的很多个分子的kET，然后还得算至少两帧的盒子担心会花过多不必要的时间。或者我这个设置的已经没问题了，就是要很大的计算量呢？

作者
Author: sobereva 时间: 2025-11-15 06:15
受力有一定收敛趋势可以再等等
CP2K基于笛卡尔坐标优化分子体系收敛远比Gaussian等量子化学程序慢得多，更何况这还是柔性大体系。一天没跑完很正常
实在搞不定的话也可以尝试考虑改用NWChem的CDFT

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