<div dir="ltr"><b>Hello,<br><br>I have a question regarding dispersion corrections for the SCC DFTB method implemented in CP2K.<br><br>I tell the program to use Dispersion correction, but the line in the output shows only zeroes like:</b><br><br><span style="font-family: courier\ new, monospace;">  Core Hamiltonian energy:                                     -4.17194433927250<br>  Repulsive potential energy:                                   0.07921577780740<br>  Electronic energy:                                            0.02107155371499<br>  Dispersion energy:                                            0.00000000000000<br><br>  Total energy:                                                -4.07165700775011<br><br> ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):               -4.071657007750285<br><br><font face="arial,sans-serif"><b>My input for a simple water molecule looks like:</b><br><br><span style="font-family: courier\ new, monospace;">&GLOBAL<br>  PROJECT H2O<br>  RUN_TYPE GEO_OPT<br>  PRINT_LEVEL LOW<br>&END GLOBAL<br><br>&FORCE_EVAL<br>  METHOD QS<br>  &DFT<br>    &MGRID<br>      CUTOFF 500<br>    &END MGRID<br>    &QS<br>      METHOD DFTB<br>      &DFTB<br>        SELF_CONSISTENT    T<br>        DO_EWALD           T<br>        DISPERSION         <br>        &PARAMETER<br>          PARAM_FILE_PATH  /home/Luyj/test/cp2k510/data/DFTB/scc/<br>          PARAM_FILE_NAME  scc_parameter<br>          DISPERSION_PARAMETER_FILE /home/Luyj/test/cp2k510/data/dftd3.dat<br>          DISPERSION_TYPE D3<br>          #UFF_FORCE_FIELD  uff_table<br>        &END PARAMETER<br>      &END DFTB<br>    &END QS<br>    &SCF<br>      MAX_SCF 30<br>      EPS_SCF 1.0E-6<br>      SCF_GUESS ATOMIC<br>      &OT<br>        MINIMIZER DIIS<br>        PRECONDITIONER FULL_SINGLE_INVERSE<br>      &END<br>    &END SCF<br>  &END DFT<br>  &SUBSYS<br>    &CELL<br>      ABC 12.4138 12.4138 12.4138<br>    &END CELL<br>    &COORD<br>      O      12.235322       1.376642      10.869880<br>      H      12.415139       2.233125      11.257611<br>      H      11.922476       1.573799       9.986994<br>    &END COORD<br>  &END SUBSYS<br>&END FORCE_EVAL<br><br>&MOTION<br>  &GEO_OPT<br>    TYPE MINIMIZATION<br>    MAX_DR    1.0E-03<br>    MAX_FORCE 1.0E-03<br>    RMS_DR    1.0E-03    <br>    RMS_FORCE 1.0E-03<br>    MAX_ITER 200<br>    OPTIMIZER CG<br>    &CG<br>      MAX_STEEP_STEPS  0<br>      RESTART_LIMIT 9.0E-01<br>    &END CG<br>  &END GEO_OPT<br>  &CONSTRAINT<br>    &FIXED_ATOMS<br>      COMPONENTS_TO_FIX XYZ<br>      LIST 1<br>    &END FIXED_ATOMS<br>  &END CONSTRAINT<br>&END MOTION </span><br></font></span><br><br><b>Is something wrong with my input or is parameter set not supported?</b><br><br></div>