计算化学公社

标题: Cp2k XTB CPASSERT failed,SCF 无法收敛问题求助，谢谢 [打印本页]

作者
Author: yinxing 时间: 2022-3-30 20:46
标题: Cp2k XTB CPASSERT failed,SCF 无法收敛问题求助，谢谢
本帖最后由 yinxing 于 2022-3-31 12:30 编辑

参考https://www.cp2k.org/howto:biochem_qmmm 的例子，进行蛋白与镧离子复合物QMMM模拟，MM与QM区域耦合采用COULOMB方法，QS 采用XTB+D3色散校正，SCF_GUESS 设置为 MOPAC， MAX_SCF 与OUT_SCF分别设置200与100，OUT_ NVT系综下采用NOSE TIMECON [fs] 100控温方法进行模拟。输出文件出现如下错误：
Total energy:                                           -389.29857543363477

  outer SCF iter =  201 RMS gradient = 0.18E+01 energy =    -389.2985754336
  outer SCF loop FAILED to converge after  201 iterations or **** steps

*** WARNING in qs_scf.F:576 :: SCF run NOT converged ***

Calculate PDOS at iteration step                               10

            ---- PDOS: start iteration on the KS states ---

ENERGY| Total FORCE_EVAL ( QMMM ) energy [a.u.]:          7971.008685129578225

*******************************************************************************
* ___                                                                      *
*  / \                                                                   *
* [ABORT]                                                                   *
*  \___/                            CPASSERT failed                         *
* |                                                                      *
*  O/|                                                                      *
* /| |                                                                      *
* / \                            motion/thermostat/thermostat_utils.F:1846 *
*******************************************************************************

===== Routine Calling Stack =====

         2 qs_mol_dyn_low
         1 CP2K
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

作者
Author: Eudaimonia 时间: 2022-3-31 09:00
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚

作者
Author: yinxing 时间: 2022-3-31 12:30

Eudaimonia 发表于 2022-3-31 09:00
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚

明白，谢谢

作者
Author: www1 时间: 2022-7-7 09:53
我也出现了一样的报错，请问楼主解决了吗~

作者
Author: yinxing 时间: 2022-9-12 14:41

www1 发表于 2022-7-7 09:53
我也出现了一样的报错，请问楼主解决了吗~

应该是QM区域的原子与Link原子出现错误，需要仔细核对。

作者
Author: tbw@zyc 时间: 2022-12-7 09:19
本帖最后由 tbw@zyc 于 2022-12-7 09:38 编辑

SCF steps 200, 提示201都没converge 可以直接double 看看，也有可能是你的QM体系原子数太多体系的配置参数不合适

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