<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">You got pretty close to what the density of SPC water is at 298K<div><br></div><div>According to "Temperature dependence of TIP3P, SPC, and TIP4P water from NPT Monte Carlo simulations: Seeking temperatures of maximum density", William L. Jorgensen, Corky Jenson, J. Comp. Chem, Volume 19 Issue 10, Pages 1179 - 1186 (1998)<div>it should be 0.985g/ml, and TIP4P will give you 1.001g/ml. </div><div><br></div><div>Remember that you're simulating SPC water, not real water. </div><div><br></div><div>Sander</div><div><br></div><div><br></div><div><br><div><div>On Feb 24, 2010, at 13:23 , Sunil Thapa wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><table cellspacing="0" cellpadding="0" border="0" style="position: static; z-index: auto; "><tbody><tr><td valign="top" style="font: inherit;"><div>Respectable Experts</div>
<div>Grettings</div>
<div>I am studying diffusion of an oxygen molecule in 255 molecule of SPC water. I need to determine diffusion coefficient at temperature 298 K and 1 atm pressure. In order to do that</div>
<div>using editconf I scaled the density of the system including oxygen to 1g/l. I checked the density after energy minimization (steepest descent) which gave the reading of 0.76 g/l. After that I equilibrated the system to pressure of 1 bar and temperature of 298 K using NPT simulation (berendsen). The density reading was 0.98. When I checked the density table at 298 K, it was about 0.9927 g/l. </div>
<div>I am afraid if my density is not consistent with experimental density of the system. Is this some bug of gromacs 3.2.1 or my understanding of thermodynamics is very poor for pressure and liquid density. Please help me on this. </div>
<div>Should I use the same NPT simulation to derive the experimental diffusion coefficient ?</div>
<div> </div>
<div>Looking forward to hearing from you</div>
<div>Neal, Nepal</div><font face="Arial" size="2"></font></td></tr></tbody></table><br>
-- <br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at http://www.gromacs.org/search before posting!<br>Please don't post (un)subscribe requests to the list. Use the <br>www interface or send it to gmx-users-request@gromacs.org.<br>Can't post? Read http://www.gromacs.org/mailing_lists/users.php</blockquote></div><br></div></div></body></html>