<div dir="ltr">Hi,<div><br></div><div>hybrid functional calculations using MOLOPT basis sets are extremely expensive. You probably need to use the Auxiliary Density Matrix Method approximation. See these <a href="https://www.cp2k.org/exercises:2016_summer_school:hfx">execises</a> and some <a href="https://mattatlincoln.github.io/talks/GhentWorkshop/#/">slides</a> maybe this <a href="https://www.youtube.com/watch?v=snG4fbpI0_g&list=PLrmNhuZo9sgYJqeTWUhdhtYu3Ol988prg&index=6&t=0s">talk</a>.</div><div><br></div><div>Matt</div><div><br></div><div><br><br>On Thursday, March 28, 2019 at 8:15:26 PM UTC, Mohammad Shakiba wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Dear CP2K users,</div><div><br></div><div>I have installed CP2K on our university cluster with intel 2019 and it works properly for tests. However, when I use it for ab initio molecular dynamics with my structure, the first 20 steps converge but the convergence in the next step takes too long and it does not converge. Actually, it becomes close to the convergence limit of 1.0E-6 (like 2.04E-6 which was the closest) but does not converge even in outer loops and oscillates between 2.0E-6 and 1.0E-4.</div><div><br></div><div> Can you help me find out why this happens and how can I fix it? is it because of my basis sets and potentials? or my structure? Is it practical to increase the MAX_SCF to a high number like 1000? What are the normal ranges for this limit?<br></div><div><br></div><div>By the way I tried both CG and DIIS but the problem still exists.<br></div><div><br></div><div>Thanks in advance.<br></div><div><br></div><div>My structure contains Pb ans S atoms -> (PbS)16. Here are some main parts of my attached inputs:</div><div><br></div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT   ! I used MOLOPT bases file<br>    POTENTIAL_FILE_NAME POTENTIAL          ! with GTH pseudopotentials and PBE functionals<br></div><div>    &MGRID<br>       CUTOFF [Ry] 400<br>    &END<br>    &QS<br>       METHOD GPW <br>       EPS_DEFAULT 1.0E-10<br>       EXTRAPOLATION ASPC <br>    &END<br>    &POISSON      <br>      PERIODIC NONE<br>      PSOLVER  MT<br>    &END POISSON<br>    &SCF                          <wbr>    <br>      SCF_GUESS ATOMIC<br>      MAX_SCF 400<br>      EPS_SCF 1.0E-6<br>      &OT<br>        PRECONDITIONER FULL_SINGLE_INVERSE<br>        MINIMIZER DIIS                          <wbr>                  ! I also tried CG<br>        LINESEARCH 3PNT<br>!        STEPSIZE 0.08<br>!        PRECOND_SOLVER INVERSE_CHOLESKY<br>      &END OT<br>      &OUTER_SCF                    <wbr>                             ! Outer SCF loop for 10 times<br>        MAX_SCF 10<br>        EPS_SCF 1.0E-6<br>      &END</div><div>    &XC<br>      &XC_GRID<br>        XC_DERIV SPLINE2<br>        XC_SMOOTH_RHO NN50<br>      &END XC_GRID<br>      &XC_FUNCTIONAL <br>         &PBE<br>         &END<br>      &END XC_FUNCTIONAL<br>    &END XC</div><div><br></div></div></blockquote></div></div>