Hi Onofrio<div><br></div><div>Question: have you set this up and run a more standard simulation to work out the convergence first?</div><div><br></div><div>If I guess, it's probably Cu as a transition metal with d-orbitals that is difficult (I am suffering with Au).<br><br></div><div><div>Kind Regards</div><div><br></div><div>Sam </div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, June 15, 2022 at 3:20:20 PM UTC+2 onofr...@gmail.com 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;">Dear CP2K community,<br><br>I'm new to CP2K and I'm trying to set up a second-generation Car-Parrinello (SGCP) MD simulation. I have started with a geometrically optimized system of toluene adsorbed on a 3-layer Cu slab (8x6). I would to know what sections/keywords in the input file must be used to correctly set up a SGCPMD simulation and, possibly, the best settings to improve the speed of such simulations since I have not found any related tutorial (but only few things in the cp2k google group). <br>If I well understood, it has to be turn on the Langevin dynamics section in combination with the OT method. For example:<br><br>&FORCE_EVAL<br>  ...<br>  &DFT<br>    ...<br>    &QS<br>      METHOD GPW<br>      EXTRAPOLATION ASPC<br>      EXTRAPOLATION_ORDER 3<br>      EPS_DEFAULT 1.0E-10<br>    &END QS<br>    ...<br>    &SCF<br>    ...<br>      &OT<br>        PRECONDITIONER FULL_SINGLE_INVERSE   ! or FULL_KINETIC<br>        MINIMIZER DIIS<br>        N_DIIS 7<br>      &END OT<br>      &OUTER_SCF<br>        MAX_SCF 3<br>        EPS_SCF 1.0E-05<br>      &END OUTER_SCF    <br>    ...<br>    &END SCF<br>    ...<br>  &END DFT<br>  ...<br>&END FORCE_EVAL<br><br>&MOTION<br>  &MD<br>    ENSEMBLE LANGEVIN<br>    ...<br>    &LANGEVIN<br>      GAMMA 0.001<br>    &END LANGEVIN<br>    ...<br>  &END MD<br>  ...<br>&END MOTION<br><br>Is it correct? To follow what I just mentioned, I prepared an input file (attached below) based on sections/keywords defined just above and I got from the output file (attached below) only <br>    <b><i>"Leaving inner SCF loop after reaching   200 steps" </i></b><div>(200 steps because I set MAX_SCF 200)<br>resulting in a very bad convergence. <br>As a result, it seems that I have not been running any SGCPMD simulation but a simple bad BOMD. Am I missing something?<br>How can I correctly set up my input file for SGCPMD?<br><br>Any help would be greatly appreciated.</div><div>Thanks.<br><br>Onofrio<br></div></blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/35d7f287-058c-4f36-a5d1-410ed43c7131n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/35d7f287-058c-4f36-a5d1-410ed43c7131n%40googlegroups.com</a>.<br />