CP2K用QM/MM方法结合TDDFT跑MD模拟报错

yjb

各位老师好，最近在用CP2K模拟荧光蛋白酶激发态模拟，想用QM/MM结合TDDFT的方法模拟酶的构象变化。先用QM/MM MD的方法跑了一段平衡（用的方法是B3LYP-D3(BJ)），平衡之后想取最后一帧结构进行QM/MM + TDDFT MD模拟，来查看酶的构象变化。遇到如下报错：


*******************************************************************************
*   ___                                                                       *
*  /   \                                                                      *
* [ABORT]                                                                     *
*  \___/                             CPASSERT failed                          *
*    |                                                                        *
*  O/|                                                                        *
* /| |                                                                        *
* / \                                               xc/xc_rho_set_types.F:309 *
*******************************************************************************


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

           11 xc_calc_2nd_deriv_analytical
           10 tddfpt_apply_xc_analytic
            9 fhxc_kernel
            8 tddfpt_compute_Aop_evects
            7 tddfpt_davidson_solver
            6 tddfpt
            5 qs_energies_properties
            4 qs_energies
            3 qs_forces
            2 qs_mol_dyn_low
            1 CP2K



想请教各位老师，这样做是不是合理？

xinlanyh

请问下你是赵老师组的学生不？这体系看着好眼熟啊