<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div><br><div><div>On 19/10/2011, at 12:08 PM, Sapna Sarupria wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><font class="Apple-style-span" size="3">Dear All,</font>
<div><font class="Apple-style-span" size="3"><br></font></div><div><font class="Apple-style-span" 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 class="Apple-style-span" size="3"><br></font></div><div><font class="Apple-style-span" size="3">MDP File </font></div><div><font class="Apple-style-span" size="3"><br></font></div><div><font class="Apple-style-span" size="3"><div>
title = 100% occupancy CO2 hydrate ; a string</div><div>cpp = /lib/cpp ; c-preprocessor</div><div>dt = 0.002 ; time step</div><div>nsteps = 15000000 ; number of steps</div>
<div>nstcomm = 10 ; reset c.o.m. motion</div><div>nstxout = 0000 ; write coords</div><div>nstvout = 0000 ; write velocities</div><div>
nstlog = 25000 ; print to logfile</div>
<div>nstenergy = 500 ; print energies</div><div>xtc_grps = OW_HW1_HW2_CO2 </div><div>nstxtcout = 1000</div><div>nstlist = 10 ; update pairlist</div>
<div>ns_type = grid ; pairlist method</div><div>coulombtype = PME</div><div>rvdw = 1.2 ; cut-off for vdw</div><div>rcoulomb = 1.2 ; cut-off for coulomb</div>
<div>rlist = 1.2 ; cut-off for coulomb</div><div>DispCorr = EnerPres</div><div>Tcoupl = Nose-Hoover</div><div>ref_t = 270 </div><div>tc-grps = System</div>
<div>tau_t = 0.5 </div><div>Pcoupl = Parrinello-Rahman</div><div>Pcoupltype = semiisotropic ; pressure geometry</div><div>tau_p = 1.0 1.0 ; p-coupoling time</div>
<div>compressibility = 4.5e-5 4.5e-5 ; compressibility</div><div>ref_p = 30.5 30.5 ; ref pressure</div><div>gen_vel = yes ; generate initial vel</div>
<div>gen_temp = 260 ; initial temperature</div><div>gen_seed = 372340 ; random seed</div><div>constraint_algorithm = shake</div><div>constraints = all-bonds </div>
<div><br></div></font></div><div><font class="Apple-style-span" size="3">-------</font></div><div><font class="Apple-style-span" size="3"><br></font></div><div><font class="Apple-style-span" size="3"><br></font></div><div>
<font class="Apple-style-span" size="3">Thank you,</font></div><div><font class="Apple-style-span" size="3"><br></font></div><div><font class="Apple-style-span" size="3">Sincerely,</font></div><div><font class="Apple-style-span" size="3">Sapna</font></div>
-- <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/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>
-----<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: <a href="mailto:Itamar.Kass@monash.edu">Itamar.Kass@monash.edu</a><br>============================================<br>
</div>
<br></div></div></body></html>