<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Hi,<br><br>I think there is still a general issue with optimizing the PME settings<br>(or settings for any electrostatics method) for MD simulations.<br>The question is what exactly one should optimize.<br>If you do a single point calculation, you might only want to minimize<br>the force error. But for MD you might also want good energy conservation.<br>Optimizing the RMS force is certainly not the best choice,<br>since that would lead to a large ewald_rtol, which gives small<br>errors before the cut-off and on most of the long-range part,<br>but large errors just beyond the cut-off. This will lead to bad<br>energy conservation. Even an absolute force error will probably<br>lead to a too large ewald_rtol.<br><br>Berk<br><br><br /><hr />See all the ways you can stay connected <a href='http://www.microsoft.com/windows/windowslive/default.aspx' target='_new'>to friends and family</a></body>
</html>