<div>Hi Victor,</div><div>from your setup the main things are:</div><div><ul><li>the cutoff is too small - should be set like a standard DFT calculation - something around 400 Ry and the REL_CUTOFF to 60 Ry or similar.</li><li>the timestep is very aggressive. If you plot your energy as function of time I think you will see the conserved quantity drift badly at some point in your simulation and then it explodes. Try something like <br /><ul><li><span style="font-family: Courier New;">stepsize [au_t] 0.1</span></li></ul></li></ul></div><div>Matt<br /></div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, 19 May 2023 at 07:48:32 UTC+1 Victor Volkov 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></div><div>Good day.<br></div><div>I am exploring molecular dynamics under field.</div><div></div><div></div><div>I set glycine at the center of a box,</div><div>and suggested Ehrenfest dynamics with a time step of about 0.04 fs,</div><div>while the custom field file sets the external field to ZERO. <br></div><div></div><div></div><div>The simulation proceeds alright, till the time-frame 500, after which</div><div>there is a peculiar structural "expansion".</div><div></div><div><b>Q1: Would you comment this?</b></div><div></div><div>Also, I attach a fragment of the momentum file:</div><div>CP2K prints reference point position to be "chaotic",</div><div>while I select 
REFERENCE COAC: since the atomic charges of the glycine, <br></div><div>the reference point should be next to the C-alpha atom. <br></div><div></div><div>
<b>Q2:  would you comment on the usage of the setting I selected?  </b><b><br></b></div><div></div><div>Thank you.</div><div>Victor<br> </div><div><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/c389ccfe-3201-4f0a-83e7-e931000e2366n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/c389ccfe-3201-4f0a-83e7-e931000e2366n%40googlegroups.com</a>.<br />