Dear Alexandra<div><br></div><div>If you properly selected the noisy gamma and the equipartition theorem is satisfied, the generated  Langevin dynamics should sample the canonical ensemble.</div><div><br></div><div>Regards</div><div>Marcella<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, June 21, 2022 at 12:44:20 PM UTC+2 a.celi...@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;"><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></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/e7b3665e-5f19-4fa0-9043-11fda28f124bn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/e7b3665e-5f19-4fa0-9043-11fda28f124bn%40googlegroups.com</a>.<br />