<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Hi Sapna,</div><div><br></div><div>It is a good question, but such of personal preferences. &nbsp;I found out that my simulations tend to crash more often when I am using 4.5.5. It might be the General-reaction -field I use or something else, I am not sure. Also, it seems I will not gain much from upgrading, so I stick to 4.0.7.</div><div><br></div><div>All this is not to say that 4.5.5 is not a good simulation package. It is a good and useful MD package, and I highly appreciate all the time and effort put into it.</div><div><br></div><div>Best,</div><div>Itamat</div><div><br></div><div>PS, there are still people around using 3.3.3 or even 3.2.1, so I am not that behind.&nbsp;</div><br><div><div>On 19/10/2011, at 1:14 PM, Sapna Sarupria wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Dear Itamar,<div><br></div><div>Thank you for your response. Why did you prefer to stick to 4.0.7? Is there something particularly difficult to transfer from 4.0.7 version to 4.5 version?</div><div><br></div><div>Sapna<br>

<br><div class="gmail_quote">On Tue, Oct 18, 2011 at 9:36 PM, Itamar Kass <span dir="ltr">&lt;<a href="mailto:itamar.kass@monash.edu">itamar.kass@monash.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<div style="word-wrap:break-word">Hi&nbsp;Sapna,<div><br></div><div>I had the same problem few weeks ago.</div><div><br></div><div>Apparently, some parameters are different (eg.&nbsp;nstpcouple).</div><div><br></div><div>Cheers,</div>

<div>Itamar<font color="#888888"><br></font><div><br></div><div>PS&nbsp;I found out that I prefer to stick to 4.0.7.</div><div><br></div><div><br></div><div></div></div></div><br><div style="word-wrap:break-word"><div><div><br>

<div><div>On 19/10/2011, at 12:08 PM, Sapna Sarupria wrote:</div><br><blockquote type="cite"><font size="3">Dear All,</font>
<div><font size="3"><br></font></div><div><font size="3">I have a simulation set up for a mixture of carbon dioxide and water system which runs perfectly on Gromacs 4.0.5 version. I have run these simulations at different temperatures and pressures, from the same starting configuration and the simulation proceeds smoothly in the 4.0.5 version. However, when I use exactly the same files to run the simulation in 4.5.4 version, the system crashes and writes out several pdb files. I was wondering if there is anything specific that has been changed between the two versions. I have pasted the mdp file of the simulations below. Any input will be much appreciated.</font></div>



<div><font size="3"><br></font></div><div><font size="3">MDP File&nbsp;</font></div><div><font size="3"><br></font></div><div><font size="3"><div>

title &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;100% occupancy CO2 hydrate ; a string</div><div>cpp &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;/lib/cpp &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; c-preprocessor</div><div>dt &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;0.002 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; time step</div><div>nsteps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;15000000 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; number of steps</div>



<div>nstcomm &nbsp; &nbsp; &nbsp;= &nbsp;10 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; reset c.o.m. motion</div><div>nstxout &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;0000 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; write coords</div><div>nstvout &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;0000 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; write velocities</div><div>


nstlog &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;25000 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; print to logfile</div>
<div>nstenergy &nbsp; &nbsp; &nbsp; = &nbsp;500 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; print energies</div><div>xtc_grps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;OW_HW1_HW2_CO2&nbsp;</div><div>nstxtcout &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;1000</div><div>nstlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;10 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; update pairlist</div>



<div>ns_type &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;grid &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; pairlist method</div><div>coulombtype &nbsp; &nbsp;= &nbsp;PME</div><div>rvdw &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;1.2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; cut-off for vdw</div><div>rcoulomb &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;1.2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; cut-off for coulomb</div>



<div>rlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;1.2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; cut-off for coulomb</div><div>DispCorr &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;EnerPres</div><div>Tcoupl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;Nose-Hoover</div><div>ref_t &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 270&nbsp;</div><div>tc-grps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;System</div>



<div>tau_t &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;0.5&nbsp;</div><div>Pcoupl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;Parrinello-Rahman</div><div>Pcoupltype &nbsp; &nbsp; &nbsp;= &nbsp;semiisotropic &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; pressure geometry</div><div>tau_p &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= &nbsp;1.0 &nbsp; 1.0 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; p-coupoling time</div>



<div>compressibility &nbsp; = &nbsp;4.5e-5 &nbsp;4.5e-5 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; compressibility</div><div>ref_p &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;30.5 &nbsp;30.5 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; ref pressure</div><div>gen_vel &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;yes &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; generate initial vel</div>



<div>gen_temp &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;260 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; initial temperature</div><div>gen_seed &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;372340 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; random seed</div><div>constraint_algorithm = shake</div><div>constraints &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp;all-bonds&nbsp;</div>



<div><br></div></font></div><div><font size="3">-------</font></div><div><font size="3"><br></font></div><div><font size="3"><br></font></div><div>

<font size="3">Thank you,</font></div><div><font size="3"><br></font></div><div><font size="3">Sincerely,</font></div><div><font size="3">Sapna</font></div>

-- <br>gmx-users mailing list &nbsp;&nbsp;&nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>

Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>Please don't post (un)subscribe requests to the list. Use the <br>

www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a></blockquote>

</div><br><div>
-----<br>"In theory, there is no difference between theory and practice. But, in practice, there is." - Jan L.A. van de Snepscheut<br><br>===========================================<br>| Itamar Kass, Ph.D.<br>| Postdoctoral Research Fellow<br>

|<br>| Department of Biochemistry and Molecular Biology<br>| Building 77 Clayton Campus<br>| Wellington Road<br>| Monash University,<br>| Victoria 3800<br>| Australia<br>|<br>| Tel: <a href="tel:%2B61%203%209902%209376" value="+61399029376" target="_blank">+61 3 9902 9376</a><br>

| Fax: <a href="tel:%2B61%203%209902%209500" value="+61399029500" target="_blank">+61 3 9902 9500</a><br>| E-mail:&nbsp;<a href="mailto:Itamar.Kass@monash.edu" target="_blank">Itamar.Kass@monash.edu</a><br>============================================<br>


</div>
<br></div></div></div><br>--<br>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></blockquote></div><br><br clear="all"><div><br></div>
</div>
-- <br>gmx-users mailing list &nbsp;&nbsp;&nbsp;<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/Support/Mailing_Lists/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/Support/Mailing_Lists</blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; ">-----<br>"In theory, there is no difference between theory and practice. But, in practice, there is." - Jan L.A. van de Snepscheut<br><br>===========================================<br>| Itamar Kass, Ph.D.<br>| Postdoctoral Research Fellow<br>|<br>| Department of Biochemistry and Molecular Biology<br>| Building 77 Clayton Campus<br>| Wellington Road<br>| Monash University,<br>| Victoria 3800<br>| Australia<br>|<br>| Tel: +61 3 9902 9376<br>| Fax: +61 3 9902 9500<br>| E-mail:&nbsp;<a href="mailto:Itamar.Kass@monash.edu">Itamar.Kass@monash.edu</a><br>============================================<br></span>
</div>
<br></body></html>