Cp2k XTB CPASSERT failed,SCF 无法收敛问题求助，谢谢

yinxing

参考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.
--------------------------------------------------------------------------

Eudaimonia

在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚

yinxing

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

明白，谢谢

www1

我也出现了一样的报错，请问楼主解决了吗~

yinxing

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

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

tbw@zyc

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