<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Berk Hess wrote:
<blockquote cite="midCOL113-W52E898FAB16DC44AEF36F18E3C0@phx.gbl"
 type="cite">
  <style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style><br>
  <br>
&gt; Date: Mon, 1 Mar 2010 15:44:28 +0100<br>
&gt; From: <a class="moz-txt-link-abbreviated" href="mailto:r.friedman@bioc.uzh.ch">r.friedman@bioc.uzh.ch</a><br>
&gt; To: <a class="moz-txt-link-abbreviated" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt; Subject: Re: [gmx-users] NVE of water<br>
&gt; <br>
&gt; Mark Abraham wrote:<br>
&gt; &gt; On 2/03/2010 12:39 AM, Ran Friedman wrote:<br>
&gt; &gt;&gt; Hi,<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; I would still argue that double precision is important.<br>
&gt; &gt;<br>
&gt; &gt; Oh? The discussion of Table 4 of<br>
&gt; &gt; <a class="moz-txt-link-freetext" href="http://pubs.acs.org/doi/abs/10.1021/ct700301q">http://pubs.acs.org/doi/abs/10.1021/ct700301q</a> (2008 GROMACS 4
JCTC<br>
&gt; &gt; paper) suggested to me that single-precision NVE could be
done well in<br>
&gt; &gt; GROMACS. Am I missing something?<br>
&gt; &gt;<br>
&gt; &gt; Mark<br>
&gt; "The presented benchmarks were performed in the NVT ensemble"
(section IX).<br>
&gt; Or am I missing something?<br>
  <br>
No, but everything that affects energy conservation in NVE also affects
it in NVT,<br>
in addition the thermostat affect the integration accuracy in NVT (and
in NVT<br>
you do not measure energy conservation from the total energy, but from
the<br>
conserved energy quantity).<br>
  <br>
Double precision can be important for energy conservation, but often
other<br>
factors deteriorate the energy conservation orders of magnitude from
what<br>
can be reached in single precision already. Double precision is only
required<br>
for testing or when you really need to generate an NVE ensemble.<br>
  <br>
Berk<br>
</blockquote>
Thanks for clarifying Berk - I was referring to the latter (i.e., when
you really need to generate an NVE ensemble). I guess this is the main
reason to run with NVE. Since most people are not running NVE with
large systems or for production simulations, the added accuracy often
justifies the use of double precision.<br>
<br>
Best,<br>
Ran<br>
<br>
</body>
</html>