<div>Hi Jonathan,</div><div><br></div><div>Did you try to use smearing? It's very useful for metallic systems. See the final of the tutorial of Energy and Force here https://www.cp2k.org/howto:static_calculation. There's shows how to use it.<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Em segunda-feira, 29 de março de 2021 às 07:02:56 UTC-3, Matt W escreveu:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Couple of comments. This is actually a rather difficult calculation to get correct - the full solid state calcs will probably me much easier to converge.<div>(i) your grid cuttoff is probably far to low for a q30 pseudo. You could also try uing GAPW method (METHOD GAPW in the QS section)</div><div>(ii) Your cell is probably too small for the MT poisson solver try increasing the size of the cell. WAVELET solver may be more efficient for atomic calcs.</div><div>(iii) ALPHA might be too big - try reducing it. Maybe then restart and increase a little, or try OT method.</div><div>Matt<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, March 22, 2021 at 2:08:15 PM UTC <a href data-email-masked rel="nofollow">jona...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br><div>Dear cp2k developers,</div><div> I am trying to test the calculation using a single atom of cerium. The coordinates should be ok, I am using a single atom at the origin of a box of 5 A of side length. Also both basis set and pseudopotential should be ok, I have found them in the data directory. I do not understand why the SCF does not converge. Surely, there is an error of mine but I cannot see it. I have put the muliplicity to 3, because in its electronic configuration Ce has 2 unpaired. electrons. This is my input.</div><div>&GLOBAL<br>  PROJECT Ce<br>  RUN_TYPE ENERGY<br>  PRINT_LEVEL LOW<br>&END GLOBAL<br>&FORCE_EVAL<br>  METHOD QS<br>  &DFT<br>    BASIS_SET_FILE_NAME BASIS_MOLOPT_UCL<br>    POTENTIAL_FILE_NAME GTH_POTENTIALS<br>    CHARGE 0<br>    MULTIPLICITY 3<br>    &MGRID<br>      CUTOFF 550<br>      REL_CUTOFF 80<br>    &END MGRID<br>    &QS<br>      EPS_DEFAULT 1.0E-12<br>      EXTRAPOLATION PS<br>      EXTRAPOLATION_ORDER 3    <br>    &END QS<br>   &SCF<br>      SCF_GUESS RESTART<br>      EPS_SCF 1.0E-7<br>      MAX_SCF 500<br>      &MIXING T<br>        METHOD BROYDEN_MIXING<br>        ALPHA 0.4<br>        NBROYDEN 8<br>      &END MIXING<br>      <br>      &DIAGONALIZATION<br>        ALGORITHM STANDARD<br>      &END DIAGONALIZATION<br>    &END SCF<br>   &XC<br>       DENSITY_CUTOFF     1.0000000000000000E-10<br>       GRADIENT_CUTOFF     1.0000000000000000E-10<br>       TAU_CUTOFF     1.0000000000000000E-10<br>       &XC_FUNCTIONAL  NO_SHORTCUT<br>        &PBE<br>            PARAMETRIZATION ORIG<br>        &END PBE<br>       &END XC_FUNCTIONAL<br>     &END XC  <br>    &POISSON<br>      PERIODIC NONE<br>      POISSON_SOLVER MT<br>    &END POISSON<br>  &END DFT<br>  &SUBSYS<br>    &CELL<br>      ABC 5. 5. 5.<br>      PERIODIC NONE<br>    &END CELL<br>    &COORD<br>    Ce 0.00 0.00 0.00<br>    &END COORD<br>&KIND Ce<br>      ELEMENT   Ce<br>      BASIS_SET DZVP-MOLOPT-SR-GTH-q30<br>      POTENTIAL GTH-PBE-q30<br>    &END KIND<br>&END SUBSYS<br>&END FORCE_EVAL<br></div><div>Have a nice day,</div><div>Jonathan<br></div></blockquote></div></blockquote></div>