<div class="gmail_quote"><div>Hello,<br><br>&nbsp;Thanks for your response. I have used the same tpr files for both the simulations. I compiled the tpr file in gromacs 3.2.1 and then used that tpr file to run simulations using the two versions of gromacs (
i.e. only the mdrun was from different versions). I am not using any dispersion correction. There is also no switch or shift being used. I have pasted the run.mdp file below. The energies are different in the two cases. The LJ(SR) and the Coulomb are both different, resulting in different potential energies (PE) for the two. The PE of gromacs 
v3.3.1 is higher than that obtained by gromacs v3.2.1 (i.e v3.2.1 is more negative). In the pressure virials, I find that the diagonal terms are the same for the two versions (Vir-XX,YY and ZZ) however the off-diagonal are significantly different (magnitude of 
v3.2.1 &gt; v3.3.1),given below. Does this information give any indication as to what may be wrong? If not I will submit the tpr file to bugzilla.<br><br>Thank you<br><br>Sapna<br><br>quantity&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; gromacs v3.2.1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; gromacs 
v3.3.1<br>Pressure (bar)&nbsp; &nbsp; 0.596018&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.838094&nbsp;&nbsp;&nbsp; <br>Vir-XX&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp; 26640.5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 26661.2&nbsp;&nbsp;&nbsp; <br>Vir-XY&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6865.12&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.36518&nbsp;&nbsp;&nbsp; <br>Vir-XZ&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6865.21&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; 0.0121227
&nbsp;&nbsp;&nbsp; <br>Vir-YX&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; 6865.8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 1.48257&nbsp;&nbsp;&nbsp; <br>Vir-YY&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; 26652.3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 26642.1&nbsp;&nbsp;&nbsp; <br>Vir-YZ&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; 6886.45&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 3.93327&nbsp;&nbsp;&nbsp; <br>Vir-ZX&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6861.23&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp; -0.0529923
&nbsp;&nbsp;&nbsp; <br>Vir-ZY&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; 6882.24&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 3.87739&nbsp;&nbsp;&nbsp; <br>Vir-ZZ&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 26665&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; 26647.2&nbsp;&nbsp;&nbsp; <br><br>Run.mdp file:<br>title&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; Water simulation&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; a string<br>cpp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; /lib/cpp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; c-preprocessor
<br>dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.002&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; time step (ps)<br>nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; number of steps<br>nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; reset c.o.m. motion<br>nstxout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 500&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; write coords
<br>nstvout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 500&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; write velocities<br>nstlog&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 100&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; print to logfile<br>nstenergy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 500&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; print energies<br>nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; =&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; update pairlist
<br>ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; =&nbsp; grid&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; pairlist method<br>coulombtype&nbsp;&nbsp;&nbsp; =&nbsp; PME<br>rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; cut-off for vdw<br>rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; cut-off for coulomb
<br>rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; cut-off for coulomb<br>Tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; Berendsen<br>tc-grps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; System<br>tau_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; =&nbsp; 0.5<br>ref_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; =&nbsp; 300<br>Pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; Berendsen
<br>Pcoupltype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; isotropic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; pressure geometry<br>tau_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; p-coupoling time<br>compressibility&nbsp; =&nbsp; 4.5e-5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; compressibility<br>ref_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 
1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; reference pressure<br>gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; yes&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; generate initial velocities<br>gen_temp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 300&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; initial temperature<br>gen_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 2432384243&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; random seed
<br><br><br><br><br><br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Date: Tue, 18 Dec 2007 08:28:52 -0500<br>From: &quot;sapna sarupria&quot; &lt;
<a href="mailto:sapna.sarupria@gmail.com">sapna.sarupria@gmail.com</a>&gt;<br>Subject: [gmx-users] Different system volumes from gromacs versions<br> &nbsp; &nbsp; &nbsp; &nbsp;3.2.1 &nbsp; and 3.3.1<br>To: <a href="mailto:gmx-users@gromacs.org">
gmx-users@gromacs.org</a><br>Message-ID:<br> &nbsp; &nbsp; &nbsp; &nbsp;&lt;<a href="mailto:3e794d50712180528p1500aefdpb080cbcacf1e537d@mail.gmail.com">3e794d50712180528p1500aefdpb080cbcacf1e537d@mail.gmail.com</a>&gt;<br>Content-Type: text/plain; charset=&quot;iso-8859-1&quot;
<br><br>Dear users,<br><br> &nbsp; I had posted this question earlier too but did not get any response. I am<br>re-posting it hoping someone can clarify this for me.<br><br> &nbsp; I performed NPT simulations of few different systems using both gromacs
<br>v3.2.1 and v3.3.1. I use exactly the same parameters, starting files,<br>topology files etc for both the versions. The systems are pure water system<br>(spc/e 10684 water molecules) and a protein-water system. However, the
<br>average volume of the system that I get is from the simulations is different<br>for the two versions. In case of pure water, gromacs v3.2.1 gives me a lower<br>volume (~311.497 nm^3) than v3.3.1 (~324.604 nm^3) and the difference is
<br>significantly higher than the fluctuations in the volume. The pressure is<br>maintained at 1 bar and the temperature is maintained at 300 K using<br>Berendsen barostat and thermostat respectively. I used SETTLE to maintain
<br>water geometry. I was wondering if someone could tell me why this happens?<br>Has there been a modification in the pressure virial calculation or<br>something related to the same that leads to these differences?<br><br>
Thank you<br><br><br>Regards<br>Sapna<br><br>--<br>Sapna Sarupria<br>Ph.D. Student - Chemical Engineering<br>Rensselaer Polytechnic Institute<br>Troy, New York 12180<br>U.S.A.<br>Ph#: (518)276-3031<br>Life isn&#39;t about finding yourself. Life is about creating yourself.
<br>George Bernard Shaw.<br>Dare to Dream<br><br><br><br>--<br>Sapna Sarupria<br>Ph.D. Student - Chemical Engineering<br>Rensselaer Polytechnic Institute<br>Troy, New York 12180<br>U.S.A.<br>Ph#: (518)276-3031<br>Life isn&#39;t about finding yourself. Life is about creating yourself.
<br>George Bernard Shaw.<br>Dare to Dream<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://www.gromacs.org/pipermail/gmx-users/attachments/20071218/f2801895/attachment-0001.html" target="_blank">
http://www.gromacs.org/pipermail/gmx-users/attachments/20071218/f2801895/attachment-0001.html</a><br><br>------------------------------<br><br>Message: 7<br>Date: Tue, 18 Dec 2007 15:13:16 +0100<br>From: David van der Spoel &lt;
<a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</a>&gt;<br>Subject: Re: [gmx-users] Different system volumes from gromacs<br> &nbsp; &nbsp; &nbsp; &nbsp;versions 3.2.1 &nbsp;and 3.3.1<br>To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">
gmx-users@gromacs.org</a>&gt;<br>Message-ID: &lt;<a href="mailto:4767D57C.30005@xray.bmc.uu.se">4767D57C.30005@xray.bmc.uu.se</a>&gt;<br>Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br><br>sapna sarupria wrote:
<br>&gt;<br>&gt; Dear users,<br>&gt;<br>&gt; &nbsp; &nbsp;I had posted this question earlier too but did not get any response.<br>&gt; I am re-posting it hoping someone can clarify this for me.<br>&gt;<br>&gt; &nbsp; &nbsp;I performed NPT simulations of few different systems using both
<br>&gt; gromacs v3.2.1 and v3.3.1. I use exactly the same parameters, starting<br>&gt; files, topology files etc for both the versions. The systems are pure<br>&gt; water system (spc/e 10684 water molecules) and a protein-water system.
<br>&gt; However, the average volume of the system that I get is from the<br>&gt; simulations is different for the two versions. In case of pure water,<br>&gt; gromacs v3.2.1 gives me a lower volume (~311.497 nm^3) than v3.3.1
<br>&gt; (~324.604 nm^3) and the difference is significantly higher than the<br>&gt; fluctuations in the volume. The pressure is maintained at 1 bar and the<br>&gt; temperature is maintained at 300 K using Berendsen barostat and
<br>&gt; thermostat respectively. I used SETTLE to maintain water geometry. I was<br>&gt; wondering if someone could tell me why this happens? Has there been a<br>&gt; modification in the pressure virial calculation or something related to
<br>&gt; the same that leads to these differences?<br>&gt;<br>Have you checked the energy differences?<br><br>Have you compared tpr files?<br><br>Is the dispersion correction on and are you using some kind of<br>shift/switch? (This has changed between 
3.2 and 3.3).<br><br>If you can not find the reason for this and it is reproducible then<br>please submit a bugzilla with a 3.2.1 tpr file that reproduces the problem.<br>&gt; Thank you<br>&gt;<br>&gt;<br>&gt; Regards<br>
&gt; Sapna<br>&gt;<br>&gt; --<br>&gt; Sapna Sarupria<br>&gt; Ph.D. Student - Chemical Engineering<br>&gt; Rensselaer Polytechnic Institute<br>&gt; Troy, New York 12180<br>&gt; U.S.A.<br>&gt; Ph#: (518)276-3031<br>&gt; Life isn&#39;t about finding yourself. Life is about creating yourself.
<br>&gt; George Bernard Shaw.<br>&gt; Dare to Dream<br>&gt;<br>&gt;<br>&gt;<br>&gt; --<br>&gt; Sapna Sarupria<br>&gt; Ph.D. Student - Chemical Engineering<br>&gt; Rensselaer Polytechnic Institute<br>&gt; Troy, New York 12180
<br>&gt; U.S.A.<br>&gt; Ph#: (518)276-3031<br>&gt; Life isn&#39;t about finding yourself. Life is about creating yourself.<br>&gt; George Bernard Shaw.<br>&gt; Dare to Dream<br>&gt;<br>&gt;<br>&gt; ------------------------------------------------------------------------
<br>&gt;<br>&gt; _______________________________________________<br>&gt; gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>&gt; <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">
http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>&gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>&gt; Please don&#39;t post (un)subscribe requests to the list. Use the
<br>&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">
http://www.gromacs.org/mailing_lists/users.php</a><br><br><br>--<br>David van der Spoel, Ph.D.<br>Molec. Biophys. group, Dept. of Cell &amp; Molec. Biol., Uppsala University.<br>Box 596, 75124 Uppsala, Sweden. Phone: &nbsp;+46184714205. Fax: +4618511755.
<br><a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</a> &nbsp; &nbsp;<a href="mailto:spoel@gromacs.org">spoel@gromacs.org</a> &nbsp; <a href="http://folding.bmc.uu.se" target="_blank">http://folding.bmc.uu.se</a><br><br><br>------------------------------
<br><br></blockquote></div><br>Sapna Sarupria<br>Ph.D. Student - Chemical Engineering<br>Rensselaer Polytechnic Institute<br>Troy, New York 12180<br>U.S.A.<br>Ph#: (518)276-3031<br>Life isn&#39;t about finding yourself. Life is about creating yourself. 
<br>George Bernard Shaw. <br>Dare to Dream