<div dir="ltr"><div>Dear Satyanarayana B,<br></div><div><br></div><div>  Indeed the energy does not "drift", it explodes. What did happen in an NVT simulation with otherwise the same input, did the simulation work fine?</div>
<div><br></div><div>  To me the cut-off energy of 350 Ry looks low, even with the PBE functional (smoother than BLYP); did you test that it is sufficient?</div><div><br></div><div>    Greetings from Zurich,</div><div><br>
</div><div>       apsi</div><div><br></div><div>PS Does some one know, are the methods in &XC_GRID ... &END XC_GRID consistent with an NPT simulation?</div><div class="gmail_extra"><br><br><div class="gmail_quote">
2013/9/19 Matt W <span dir="ltr"><<a href="mailto:MattWa...@gmail.com" target="_blank">MattWa...@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>Hi,<br></div><div><br></div><div>well your system is clearly exploding. Your cell_ref parameters look inconsistent with your actual cell - maybe the stress tensor is very bad!</div><div class="im"><div>
<br></div><div>  &SUBSYS</div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>    &CELL</div><div>      ABC 21.8620 16.2099 43.6020</div>
<div>      ALPHA_BETA_GAMMA 90.0 92.7540 90.0</div><div>     &CELL_REF</div><div>      ABC 22.8620 17.2099 44.6020</div><div>      ALPHA_BETA_GAMMA 90.0 126.0 90.0</div><div>     &END CELL_REF</div><div>    </div>
<div>    &END CELL</div><div><br></div></div></div></blockquote></div></div><div class="HOEnZb"><div class="h5">

<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%2Bun...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-<br>  Ari P Seitsonen / <a href="mailto:Ari.P.S...@iki.fi">Ari.P.S...@iki.fi</a> / <a href="http://www.iki.fi/~apsi/">http://www.iki.fi/~apsi/</a><br>
  Physikalisch-Chemisches Institut der Universität Zürich<br>  Tel: +41 44 63 55 44 97  /  Mobile: +41 79 71 90 935
</div></div>