Hi,<div><br></div><div>I assume you'll add a KPOINTS section in the future, the most obvious improvement is to use a higher electronic temperature. Probably the most commonly used figure in literature is 0.2 eV which is ~2300 K, you may of course be able to use a lower one. You might also see some improvement from a smaller alpha (0.01 - 0.1), these settings can require a bit of experimentation to optimize. I prefer Broyden mixing over Pulay - but I sometimes get jobs that hang when using Broyden for cell optimizations and switch to Pulay or Kerker instead. My 'preference' is based solely on Broyden tending to require fewer SCF cycles and being more consistent, it's not a large improvement but adds up over 1000's of MD steps.</div><div><br></div><div>Marcella recommended this recently for a calculation on bulk Au (inside the &SCF section),</div><div><br></div><div><div>      ADDED_MOS 200</div><div>      &DIAGONALIZATION T</div><div>        ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &MIXING T</div><div>        METHOD BROYDEN_MIXING</div><div>        ALPHA 0.01</div><div>        BETA 0.5</div><div>        NBUFFER 8</div><div>      &END MIXING</div><div>      &SMEAR</div><div>        METHOD FERMI_DIRAC</div><div>        ELECTRONIC_TEMPERATURE 500</div><div>      &END SMEAR</div></div><div><br></div>Another recommendation from Matthias for bulk Cu, includes an example KPOINT section which you can adapt to your calculation,<br><br>&DFT<br> BASIS_SET_FILE_NAME BASIS_MOLOPT <br> POTENTIAL_FILE_NAME GTH_POTENTIALS<br> &KPOINTS<br>  SCHEME MONKHORST-PACK 2 2 2<br>  FULL_GRID yes<br>  SYMMETRY yes<br>  VERBOSE yes<br>  PARALLEL_GROUP_SIZE -1<br> &END KPOINTS<br> &MGRID<br>  NGRIDS 5<br>  CUTOFF 400.0<br>  REL_CUTOFF 60.0<br> &END MGRID<br> &QS<br>  EPS_DEFAULT 1.0E-12<br>  EXTRAPOLATION use_prev_p<br> &END QS<br> &SCF<br>  ADDED_MOS 60<br>  EPS_SCF 1.0E-8<br>  MAX_SCF 300<br>  SCF_GUESS restart<br> &DIAGONALIZATION yes<br>  ALGORITHM STANDARD<br> &END DIAGONALIZATION<br> &MIXING yes<br>  ALPHA 0.4<br>  BETA 1.0<br>  METHOD broyden_mixing<br>  NBROYDEN 8<br> &END MIXING<br> &SMEAR on<br>  METHOD FERMI_DIRAC<br>  ELECTRONIC_TEMPERATURE [K] 2000.0<br> &END SMEAR<div><br></div><div><div>-T<br><br></div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, January 27, 2021 at 1:21:41 PM UTC-5 GENG YUAN 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;">Hello everyone and hope everything goes well!<div><br></div><div>I started learning CP2K recently and trying to run some DFT simulations. I am currently working on performing a cell optimization of bulk Ni, but I'm not getting an inner loop SCF convergence. I saw some similar problems in this discussion forum and modified my input file based on some suggestions from other discussions, such as adding smearing, changing the SCF convergence criterion, etc. Nothing has worked so far. Moreover, I thought my simulations won't take quite long as each of them only has 32 Ni atoms, but none of them finished within 30 minutes (I set up the total wall time for 30 minutes, so all simulations will be terminated if they didn't finish in 30 mins).</div><div><br></div><div>I attached both my input file and xyz file here for your reference, my bulk Ni model contains 2X2X2 unit cell (the lattice constant I use for Ni in my work is 3.52 Angstrom, the structure of my Ni is FCC), therefore I made #ABC 7.04 7.04 7.04 in my input file under &cell subsection. Any suggestions on how to get convergence are welcome!</div><div><br></div><div>Many thanks in advance,</div><div>Sincerely,</div><div>Geng Yuan</div><div><br></div><div><br></div></blockquote></div>