<div dir="ltr">Dear All,<br><br>I am trying to run the single point energy calculation using wB97X(-rV) functional. However I got an error and I don't know how to deal with it:<br> *******************************************************************************<br> *   ___                                                                       *<br> *  /   \                                                                      *<br> * [ABORT]                                                                     *<br> *  \___/                                                                      *<br> *    |                                                                        *<br> *  O/|                                                                        *<br> * /| |                                                                        *<br> * / \                                                    hfx_admm_utils.F:695 *<br> *******************************************************************************<br><br><br> ===== Routine Calling Stack ===== <br><br>            8 hfx_admm_init<br>            7 qs_ks_build_kohn_sham_matrix<br>            6 rebuild_ks_matrix<br>            5 qs_ks_update_qs_env<br>            4 scf_env_do_scf_inner_loop<br>            3 scf_env_do_scf<br>            2 qs_energies<br>            1 CP2K<br><br>Probably the error is for the initialization with Auxiliary Density Matrix Methods. I tried to change different ADMM_PURIFICATION_METHOD but it doesn't help.<br>Could you please give me any suggestion on what I should change in the input file?<br><br>Thanks,<br>Huy<br><br>Please see the input below (only for wB97X)<br><br>&FORCE_EVAL<br>  METHOD Quickstep<br>  &DFT<br>    BASIS_SET_FILE_NAME BASIS_MOLOPT <br>    BASIS_SET_FILE_NAME BASIS_ADMM <br>    BASIS_SET_FILE_NAME BASIS_ADMM_MOLOPT<br>    POTENTIAL_FILE_NAME GTH_POTENTIALS <br><br>    &MGRID<br>      CUTOFF 600<br>    &END MGRID<br><br>    &QS<br>      EPS_DEFAULT 1.0E-13<br>      EXTRAPOLATION ASPC <br>      EXTRAPOLATION_ORDER 2<br>    &END QS<br>   <br>    &AUXILIARY_DENSITY_MATRIX_METHOD<br>      ADMM_PURIFICATION_METHOD CAUCHY<br>    &END<br><br>    &SCF <br>       SCF_GUESS  RESTART <br>       EPS_SCF 3.0E-7 <br>       MAX_SCF 500 <br>       ADDED_MOS 500 <br>       CHOLESKY INVERSE <br>       &SMEAR  ON <br>         METHOD FERMI_DIRAC <br>         ELECTRONIC_TEMPERATURE [K] 300 <br>       &END SMEAR <br>       &DIAGONALIZATION <br>          ALGORITHM STANDARD <br>       &END DIAGONALIZATION<br>       &MIXING <br>           METHOD  DIRECT_P_MIXING<br>           ALPHA   0.3 <br>           BETA    1.5 <br>           NBROYDEN  16 <br>       &END <br><br>     &END SCF <br><br>   &XC   <br>     &XC_FUNCTIONAL<br>        &BECKE97<br>            PARAMETRIZATION  wB97X-V<br>            SCALE_X 1.0<br>            SCALE_C 1.0<br>        &END <br>     &END XC_FUNCTIONAL<br>      &HF<br>        &SCREENING<br>          EPS_SCHWARZ 1.0E-7<br>          SCREEN_ON_INITIAL_P TRUE<br>        &END<br>        &MEMORY<br>          MAX_MEMORY 10000<br>          EPS_STORAGE_SCALING 0.1<br>        &END<br>        &INTERACTION_POTENTIAL<br>          OMEGA 0.3<br>          SCALE_LONGRANGE 0.833<br>          SCALE_COULOMB 0.167<br>          POTENTIAL_TYPE MIX_CL<br>          CUTOFF_RADIUS 6.0<br>          T_C_G_DATA t_c_g.dat<br>        &END<br>        FRACTION 0.167<br>      &END<br><br>    &END XC<br>  &END DFT<br><br>  &SUBSYS<br>    &CELL<br>      ABC 12.3250085893  12.3250085893  25.0<br>      ALPHA_BETA_GAMMA 90.0 90.0 120.0<br>    &END CELL<br><br>    &TOPOLOGY<br>      COORD_FILE_NAME file.xyz<br>      COORD_FILE_FORMAT XYZ<br>    &END<br><br>    &KIND O<br>      BASIS_SET DZVP-MOLOPT-GTH<br>      AUX_FIT_BASIS_SET  FIT3<br>      POTENTIAL GTH-PBE-q6 <br>    &END KIND<br><br>    &KIND H<br>      BASIS_SET DZVP-MOLOPT-GTH<br>      AUX_FIT_BASIS_SET  FIT3<br>      POTENTIAL GTH-PBE-q1<br>    &END KIND<br>    &KIND C<br>      BASIS_SET DZVP-MOLOPT-GTH<br>      AUX_FIT_BASIS_SET  FIT3<br>      POTENTIAL GTH-PBE-q4<br>    &END KIND<br><br>  &END SUBSYS<br>&END FORCE_EVAL<br><br>&GLOBAL<br>  PROJECT WG<br>  RUN_TYPE ENERGY <br>  PRINT_LEVEL medium <br>&END GLOBAL<br><br><br></div>