<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.6944.0">
<TITLE>RE: [gmx-users] Water density</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->

<P><FONT SIZE=2>Hi Viswanadham Sridhara and Michael Patra,<BR>
<BR>
Thank you very much for your reply. I have solved the problem. It was just spc water model which dose not give good estimate of density at room temperature. I am now using spce which is giving correct density even in vacuum.<BR>
<BR>
Thanking you all,<BR>
<BR>
<BR>
Abu Naser<BR>
<BR>
School Of Life Sciences<BR>
Heriot-Watt University<BR>
Edinburgh EH14 4AS<BR>
Email: mn2@hw.ac.uk<BR>
Phone: +44(0)1314518265<BR>
Fax : +44(0) 131 451 3009<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
<BR>
-----Original Message-----<BR>
From: gmx-users-bounces@gromacs.org on behalf of Viswanadham Sridhara<BR>
Sent: Mon 20/02/2006 5:11 PM<BR>
To: Discussion list for GROMACS users<BR>
Subject: Re: [gmx-users] Water density<BR>
<BR>
This paper gives lot of information dealing with SPC and other water models.<BR>
Hope it helps you.<BR>
David van der Spoel, Paul J. van Maaren and Herman J.C. Berendsen<BR>
A systematic study of water models for molecular simulation: Derivation of<BR>
water models optimized for use with a reaction field.<BR>
J Chem Phys 108 pp. 10220-10230 (1998)<BR>
<BR>
<BR>
<BR>
On 2/20/06, David van der Spoel &lt;spoel@xray.bmc.uu.se&gt; wrote:<BR>
&gt;<BR>
&gt; Naser, Md Abu wrote:<BR>
&gt; &gt; Thanks david for reply. I have been trying to fix density of a water<BR>
&gt; &gt; box at 300k temp. At first I scaled&nbsp; the density of spc216.gro using<BR>
&gt; &gt; editconf and then I have created spc water box using genbox and scaled<BR>
&gt; &gt; the water box using editconf -density 1000.&nbsp; Finally I allowed<BR>
&gt; &gt; some vacuum in z direction using editconf again and run for a short<BR>
&gt; &gt; time( 2ns) under no pressure and checked using g_density which giving<BR>
&gt; &gt; me less density than what I start off. Is it because of vaccum?<BR>
&gt; Of course. Same number of waters, larger volume. Run with NPT to test<BR>
&gt; the density. Compare to literature.<BR>
&gt; &gt;<BR>
&gt; &gt; regards,<BR>
&gt; &gt; Abu Naser<BR>
&gt; &gt; School Of Life Sciences<BR>
&gt; &gt; Heriot-Watt University<BR>
&gt; &gt; Edinburgh EH14 4AS<BR>
&gt; &gt; Email: mn2@hw.ac.uk &lt;<A HREF="mailto:mn2@hw.ac.uk">mailto:mn2@hw.ac.uk</A>&gt;<BR>
&gt; &gt; Phone: +44(0)1314518265<BR>
&gt; &gt; Fax : +44(0) 131 451 3009<BR>
&gt; &gt;<BR>
&gt; &gt;<BR>
&gt; &gt; ------------------------------------------------------------------------<BR>
&gt; &gt; *From:* gmx-users-bounces@gromacs.org on behalf of David van der Spoel<BR>
&gt; &gt; *Sent:* Mon 20/02/2006 2:02 PM<BR>
&gt; &gt; *To:* Discussion list for GROMACS users<BR>
&gt; &gt; *Subject:* Re: [gmx-users] Water density<BR>
&gt; &gt;<BR>
&gt; &gt; Naser, Md Abu wrote:<BR>
&gt; &gt; &gt; Hi All,<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt; When I run a water box after fixing its density(1kg/L) using editconf,<BR>
&gt; &gt; &gt; the density of the box always is going down to 955gm/L. I have been<BR>
&gt; &gt; &gt; wondering why this is happening. Can anyone give some clue?<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; details please.<BR>
&gt; &gt;<BR>
&gt; &gt; &gt; Abu Naser<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt; School Of Life Sciences<BR>
&gt; &gt; &gt; Heriot-Watt University<BR>
&gt; &gt; &gt; Edinburgh EH14 4AS<BR>
&gt; &gt; &gt; Email: mn2@hw.ac.uk<BR>
&gt; &gt; &gt; Phone: +44(0)1314518265<BR>
&gt; &gt; &gt; Fax : +44(0) 131 451 3009<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt; Abu Naser<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt; School Of Life Sciences<BR>
&gt; &gt; &gt; Heriot-Watt University<BR>
&gt; &gt; &gt; Edinburgh EH14 4AS<BR>
&gt; &gt; &gt; Email: mn2@hw.ac.uk<BR>
&gt; &gt; &gt; Phone: +44(0)1314518265<BR>
&gt; &gt; &gt; Fax : +44(0) 131 451 3009<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt;<BR>
&gt; ------------------------------------------------------------------------<BR>
&gt; &gt; &gt;<BR>
&gt; &gt; &gt; _______________________________________________<BR>
&gt; &gt; &gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<BR>
&gt; &gt; &gt; <A HREF="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</A><BR>
&gt; &gt; &gt; Please don't post (un)subscribe requests to the list. Use the<BR>
&gt; &gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt; &gt; &gt; Can't post? Read <A HREF="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</A><BR>
&gt; &gt;<BR>
&gt; &gt;<BR>
&gt; &gt; --<BR>
&gt; &gt; David.<BR>
&gt; &gt; ________________________________________________________________________<BR>
&gt; &gt; David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,<BR>
&gt; &gt; Dept. of Cell and Molecular Biology, Uppsala University.<BR>
&gt; &gt; Husargatan 3, Box 596,&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 75124 Uppsala, Sweden<BR>
&gt; &gt; phone:&nbsp; 46 18 471 4205&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fax: 46 18 511 755<BR>
&gt; &gt; spoel@xray.bmc.uu.se&nbsp;&nbsp;&nbsp; spoel@gromacs.org&nbsp;&nbsp; <A HREF="http://folding.bmc.uu.se">http://folding.bmc.uu.se</A><BR>
&gt; &gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<BR>
&gt; &gt; _______________________________________________<BR>
&gt; &gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<BR>
&gt; &gt; <A HREF="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</A><BR>
&gt; &gt; Please don't post (un)subscribe requests to the list. Use the<BR>
&gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt; &gt; Can't post? Read <A HREF="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</A><BR>
&gt; &gt;<BR>
&gt; &gt; ------------------------------------------------------------------------<BR>
&gt; &gt;<BR>
&gt; &gt; _______________________________________________<BR>
&gt; &gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<BR>
&gt; &gt; <A HREF="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</A><BR>
&gt; &gt; Please don't post (un)subscribe requests to the list. Use the<BR>
&gt; &gt; www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt; &gt; Can't post? Read <A HREF="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</A><BR>
&gt;<BR>
&gt;<BR>
&gt; --<BR>
&gt; David.<BR>
&gt; ________________________________________________________________________<BR>
&gt; David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,<BR>
&gt; Dept. of Cell and Molecular Biology, Uppsala University.<BR>
&gt; Husargatan 3, Box 596,&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 75124 Uppsala, Sweden<BR>
&gt; phone:&nbsp; 46 18 471 4205&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fax: 46 18 511 755<BR>
&gt; spoel@xray.bmc.uu.se&nbsp;&nbsp;&nbsp; spoel@gromacs.org&nbsp;&nbsp; <A HREF="http://folding.bmc.uu.se">http://folding.bmc.uu.se</A><BR>
&gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<BR>
&gt;<BR>
&gt;<BR>
&gt; _______________________________________________<BR>
&gt; gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<BR>
&gt; <A HREF="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</A><BR>
&gt; Please don't post (un)subscribe requests to the list. Use the<BR>
&gt; www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt; Can't post? Read <A HREF="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</A><BR>
&gt;<BR>
<BR>
<BR>
<BR>
--<BR>
Viswanadham Sridhara,<BR>
Graduate Research Assistant,<BR>
Old Dominion University, &quot;VIRGINIA&quot;.<BR>
<BR>
<BR>
</FONT>
</P>

</BODY>
</HTML>