<div dir="ltr"><div>Hi,</div><div><br></div><div>Dear Matt,</div><div><br></div><div>Your video and also the videos inside the playlist were very useful and I learned so much and thank you so much for that. Although I modified my input but I have to recompile my cp2k version with libint and elpa as well.</div><div><br></div><div>What I noticed in the videos was that you used SCF_GUESS RESTART and WFN_RESTART_FILE_NAME ${project}-RESTART.wfn. Now I have some other question in my mind:</div><div><br></div><div>** Does using a RESTART file make the calculations and convergence faster and better? How can I obtain such a file before starting AIMD because I found no such file but I have, for example, *.restart.bak.1 file (after running calculations) but not *.wfn file? and I don't know if I can use it in the EX_RESTART section and obtain good convergence.<br></div><div><br></div><div>** As I saw the "RESTART", I remember that in some papers they have equilibriated the system for some pico-seconds with time step of 10fs and then continued the calculations with smaller time steps like 2fs. To this end I guess that I have to make two inputs: one with 10fs time step and the other with a EXT_RESTART file from the previous and 2fs time step. But I don't know what to put for SCF_GUESS for the second input? ATOMIC?</div><div><br></div><div><br></div><div><br></div><br>On Friday, March 29, 2019 at 4:17:26 PM UTC+4:30, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fexercises%3A2016_summer_school%3Ahfx\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHk7IF1iZmbRac3G0ikELVzzWNnZQ';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fexercises%3A2016_summer_school%3Ahfx\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHk7IF1iZmbRac3G0ikELVzzWNnZQ';return true;">execises</a> and some <a href="https://mattatlincoln.github.io/talks/GhentWorkshop/#/" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmattatlincoln.github.io%2Ftalks%2FGhentWorkshop%2F%23%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHhluJ3pPzhY3GkOYMptEx9gQig3g';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmattatlincoln.github.io%2Ftalks%2FGhentWorkshop%2F%23%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHhluJ3pPzhY3GkOYMptEx9gQig3g';return true;">slides</a> maybe this <a href="https://www.youtube.com/watch?v=snG4fbpI0_g&list=PLrmNhuZo9sgYJqeTWUhdhtYu3Ol988prg&index=6&t=0s" target="_blank" rel="nofollow" onmousedown="this.href='https://www.youtube.com/watch?v\x3dsnG4fbpI0_g\x26list\x3dPLrmNhuZo9sgYJqeTWUhdhtYu3Ol988prg\x26index\x3d6\x26t\x3d0s';return true;" onclick="this.href='https://www.youtube.com/watch?v\x3dsnG4fbpI0_g\x26list\x3dPLrmNhuZo9sgYJqeTWUhdhtYu3Ol988prg\x26index\x3d6\x26t\x3d0s';return true;">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></blockquote></div>