<div>Dear cp2k-users,</div><div>thanks to Marcell's suggestions and to Onofrio:</div><div><br></div><div><div>-Always start sgcpmd from a well converged wfn obtained along a standard BOMD run with diagonalisation and smearing</div><div>-  Let the first steps run on many scf iterations if needed, i.e., set MAX_SCF large enough</div><div>-
  But use a relative large eps_scf. The wfn is not going to be optimised
 as in the BOMD, energy convergence  around 10^-5 is desirable</div><div>- Use FULL_ALL preconditioner + ASPC</div><div>- Use a small time step</div><div>- Tune properly GAMMA and NOISY_GAMMA</div></div><div><br></div><div>Now my calculations are converging! Thanks a lot! Per MD step roughly 5 iterations.<br></div><div>My system is water/Pt(111),  therefore I used the GAMMA and NOISY_GAMMA from Lan, Jinggang, Jürg Hutter, and Marcella Iannuzzi. "First-principles simulations of an aqueous CO/Pt (111) interface." <i>The Journal of Physical Chemistry C</i> 122.42 (2018): 24068-24076.</div><div><br></div><div>My question is regarding the cons. Q. In ENSEMBLE NVT the cons. Q is equal to the total energy of my system plus energy of thermostat. The energy of the thermostat can be also printed (for only NVT ENSEMBLE) and I was able to verify it. In ENSEMBLE LANGEVIN, which is also at constant temperature, the cons. Q is equal to the total energy of my system. Is that correct? <br></div><div>Maybe I am missing something, because I would interpret it as my system is NVE instead of NVT. Can anybody enlighten me, please?</div><div><br></div><div>I attached the plots of the energies (Hartree units), E_T = cons. Q - Ek - Ep in case of Langevin and NVT ensemble. I know that is a very short simulations but it is just for getting some understanding.</div><div><br></div><div>Many thanks in advance,</div><div>Alexandra<br><em></em></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/6842685a-07a4-4b1f-856f-da07e61f1230n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6842685a-07a4-4b1f-856f-da07e61f1230n%40googlegroups.com</a>.<br />