<p>Dear CP2K community,</p>
<p>I am currently using CP2K to compute the interaction energy of a chloroform dimer with BSSE correction, following the counterpoise method. However, I am encountering serious SCF convergence issues when using the SCAN meta-GGA functional.</p><p>Has anyone experienced similar issues with SCAN functional in dimer calculations? Is there any recommended strategy for improving SCF convergence in these cases? Any advice on practical workarounds would be greatly appreciated.</p><p>Thank you in advance for your help.</p><p>Here’s my input file:</p>&GLOBAL<br />  PROJECT D_dimer_DISTANCE<br />  RUN_TYPE BSSE<br />&END GLOBAL<br /><br />&FORCE_EVAL<br />  METHOD Quickstep<br />  &DFT<br />    BASIS_SET_FILE_NAME  BASIS_MOLOPT<br />    POTENTIAL_FILE_NAME  POTENTIAL<br />    CHARGE 0<br />    MULTIPLICITY 1<br /><br />    &MGRID<br />      CUTOFF 700<br />    &END MGRID<br /><br />    &QS<br />      METHOD GPW<br />      EPS_DEFAULT  1.0E-8<br />      EPS_PGF_ORB  1.0E-8<br />    &END QS<br /><br />    &SCF<br />      SCF_GUESS ATOMIC<br />      EPS_SCF 1.0E-7<br />      MAX_SCF 3000<br />      &OT<br />        MINIMIZER       DIIS<br />        SAFE_DIIS       T <br />        PRECONDITIONER  FULL_SINGLE_INVERSE<br />        ENERGY_GAP      0.002<br />        STEPSIZE        0.05<br />      &END<br />      &OUTER_SCF<br />        EPS_SCF 1.0E-7<br />        MAX_SCF 500<br />      &END<br />    &END SCF<br />    <br />    &POISSON            <br />      PERIODIC NONE<br />      PSOLVER  IMPLICIT   <br />    &END POISSON<br /><br />    &XC<br />      &XC_FUNCTIONAL <br />        &MGGA_X_SCAN<br />        &END MGGA_X_SCAN<br />        &MGGA_C_SCAN<br />        &END MGGA_C_SCAN<br />      &END XC_FUNCTIONAL<br />    &END XC<br />  &END DFT<br /><br />  &SUBSYS<br />    &CELL<br />      ABC 16.0 10.0 10.0<br />      PERIODIC NONE           ! Non periodic calculation.<br />    &END CELL<br /><br />    &COORD<br />C    0.0000000    0.0000000    0.0000000<br />H    1.0817310    0.0000000    0.0000000<br />Cl  -0.5404300    1.6724230    0.0000000<br />Cl  -0.5404300   -0.8362120   -1.4483610<br />Cl  -0.5404300   -0.8362120    1.4483610<br />C    3.0000000    0.0000000    0.0000000<br />H    4.0818090    0.0000000    0.0000000<br />Cl   2.4467380    1.6699320    0.0000000<br />Cl   2.4467380   -0.8349660   -1.4462040<br />Cl   2.4467380   -0.8349660    1.4462040<br />    &END COORD<br /><br />    <br />    &KIND C_ghost<br />      ELEMENT C<br />      GHOST T<br />      BASIS_SET TZVP-MOLOPT-GTH-q4<br />      POTENTIAL GTH-SCAN-q4<br />    &END KIND<br /><br />    &KIND H_ghost<br />      ELEMENT H<br />      GHOST T<br />      BASIS_SET TZVP-MOLOPT-GTH-q1<br />      POTENTIAL GTH-SCAN-q1    <br />    &END KIND<br /><br />    &KIND Cl_ghost<br />      ELEMENT Cl<br />      GHOST T<br />      BASIS_SET TZVP-MOLOPT-GTH-q7<br />      POTENTIAL GTH-SCAN-q7    <br />    &END KIND<br /><br />    &KIND C<br />      ELEMENT C<br />      BASIS_SET TZVP-MOLOPT-GTH-q4<br />      POTENTIAL GTH-SCAN-q4<br />    &END KIND<br /><br />    &KIND H<br />      ELEMENT H<br />      BASIS_SET TZVP-MOLOPT-GTH-q1<br />      POTENTIAL GTH-SCAN-q1    <br />    &END KIND<br /><br />    &KIND Cl<br />      ELEMENT Cl<br />      BASIS_SET TZVP-MOLOPT-GTH-q7<br />      POTENTIAL GTH-SCAN-q7    <br />    &END KIND<br />  &END SUBSYS<br /><br />  &BSSE<br />    &FRAGMENT<br />      LIST 1..5  !  monomer A<br />    &END FRAGMENT<br /><br />    &FRAGMENT<br />      LIST 6..10  !  monomer B<br />    &END FRAGMENT<br /><br />    &CONFIGURATION<br />      GLB_CONF 1 1    ! First & Second FRAGMENT are being used.<br />      SUB_CONF 1 1    ! E_AB, dimer energy.<br />      CHARGE 0<br />      MULTIPLICITY 1<br />    &END CONFIGURATION<br /><br />    &CONFIGURATION<br />      GLB_CONF 1 0    ! Only First FRAGMENT is being used.<br />      SUB_CONF 1 0    ! E_A, B as ghost, monomer A energy.<br />      CHARGE 0<br />      MULTIPLICITY 1<br />    &END CONFIGURATION<br /><br />    &CONFIGURATION<br />      GLB_CONF 0 1    ! Only Second FRAGMENT is being used.<br />      SUB_CONF 0 1    ! E_B, A as ghost, monomer B energy.<br />      CHARGE 0<br />      MULTIPLICITY 1<br />    &END CONFIGURATION<br /><br />    &CONFIGURATION<br />      GLB_CONF 1 1    ! First & Second FRAGMENT are being used.<br />      SUB_CONF 1 0    ! E_AB(A), B as ghost, monomer A + ghost B<br />      CHARGE 0<br />      MULTIPLICITY 1<br />    &END CONFIGURATION<br /><br />    &CONFIGURATION<br />      GLB_CONF 1 1    ! First & Second FRAGMENT are being used.<br />      SUB_CONF 0 1    ! E_AB(B), A as ghost, monomer B + ghost A<br />      CHARGE 0<br />      MULTIPLICITY 1<br />    &END CONFIGURATION<br />    <br />    &FRAGMENT_ENERGIES <br />    &END FRAGMENT_ENERGIES<br />    <br />    ! E_BSSE = E_AB - (E_A + E_B) + (E_AB(A) + E_AB(B))<br />   <br />    #&PRINT<br />    <br />    #&END PRINT<br />  &END BSSE<br />&END FORCE_EVAL<br /><p><br /></p><p><br /></p><p>

</p><p>Best regards,</p><p>
Sam</p>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/8eac5297-87a8-4833-86fa-8396a6222360n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/8eac5297-87a8-4833-86fa-8396a6222360n%40googlegroups.com</a>.<br />