<div class="gmail_quote">Hello!</div>
<div class="gmail_quote"> </div>
<div class="gmail_quote">On Mon, Mar 10, 2008 at 8:53 PM, Diederica <<a href="mailto:diederic...@gmail.com" target="_blank">diederic...@gmail.com</a>> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><br><br><br>---------- Forwarded message ----------<br>From: Teodoro Laino <<a href="mailto:teodor...@gmail.com" target="_blank">teodor...@gmail.com</a>><br>
Date: 10 mrt, 19:03<br>Subject: nose thermostat<br>To: cp2k<br><br><br>Dear Diederica,<br><br>Let's go step by step:<br><br>(1) First and most important.. there should not be differences  <br>between old and new versions of cp2k. If this is the case it may  <br>
easily be<br>that we introduced a bug during the changes. I would highly  <br>appreciate if you could send us the whole input file including<br>geometry.<br>In this way we can identify the problem and fix it.<br></blockquote>

<div> </div>
<div>These are added, as an attachment this time</div>
<div> </div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><span></span><br>(2) The drift that you observe depends very probably on the semi-<br>empirical (Ben can comment on that.. we spent a couple of<br>
days trying to identify something similar but for NVE runs). Are you  <br>doing a condensed phase? or an isolated molecule?<br>How large is your molecule compared to the box?<br>Most important, if you switch to NVE do you also observe the drift?<br>
</blockquote>
<div> </div>
<div>I'm modelling an isolated molecule, it should fit in the box. Using NVE without temperature rescaling, with an analogue structure, an energy drift was observed, that got worse after about 20 ps, though no non-convergence was observed. The input file and an energy-plot are added.</div>

<div> </div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><span></span><br>(3) Nose has been tested extensively and it conserves the energy. So  <br>must necessarily be the quality of your forces.</blockquote>

<div> </div>
<div>The potential and kinetic energy are conserved indeed, but the conserved quantity is not</div>
<div> </div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><span></span><br><br>For the next time, can you please put tables and input files as  <br>attachment (maybe .tgz) so that the formatting of mail readers are  <br>
not changed and the input files not digested ? Thanks..<br><br>Ciao,<br>teo</blockquote>
<div> </div>
<div>Thanks for helping!</div>
<div> </div>
<div>Greetz,</div>
<div>Diederica</div></div>
<p><br> </p>