<div dir="ltr">Hi,<div><br></div><div>I did a simulation with NVE ensemble. The energy plot is shown in figure (Energy.jpg).</div><div>The plot should be flat one in NVE simulation, but my plot has a fluctuation.<br>Is there any problem in my input file (attached file)?   </div><div>Is the energy converged? Is the fluctuation reasonable?</div><div><br></div><div>Could anyone give me some suggestions?</div><div>Thank you very much.</div><div><br></div><div>Best wishes,</div><div>Zuoheng</div><div><br>在 2018年9月13日星期四 UTC+2下午4:23:29，zuoh...@gmail.com写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><code>Hello everyone,<br><br>Recently, I want to confirm the suitable time step, when I use core shell model for OH- ions and SPC/Fw for water molecules. <br>I have tested if I use core shell model for water, 1 fs of time step is too large. The energy of NVE simulation is not converged significantly.<br>However, after I transfered the SPC/Fw to replace the core-shel model for water, the energy of NVE simulation seemed to be converged (shown in figure, blue line for the energy plot and orange line for the averaged energy).<br>There are 1 Fe3+, 3 OH-(core shell model) and 28 water molecules(SPC/Fw) in the system. <br>I am not sure if the energy is converged? or the fluctuation is too strong? <br>If it is converged, can I use 1 fs of time step to do next simulation?<br>Could anyone give me some suggestions?<br><br>Thank you very much!!<br><br>Best,<br>Zuoheng  <br><br></code><p style="text-align:center;clear:both"><a style="clear:left;margin-bottom:1em;float:left;margin-right:1em"><img src="https://groups.google.com/group/cp2k/attach/867783c4b3b00/Fe3OH3.png?part=0.1&view=1&authuser=0" alt="Fe3OH3.png" width="320" height="211" border="0"></a></p></div></blockquote></div></div>