Dear Mark,<br><br>I have checked this possibility and it turned out that the positional restraints and constraints did not lead to MD instabilities when run in serial mode. Unfortunately, the system is pretty large (it consists of 27255 atoms) so it runs quite slowly in serial mode. That is why, I do want to run mdrun in parallel mode and make an advantage of our cluster. Thus, I would like to ask you for some advice how to get around the problem.<br>
<br>Philip<br><span class="gmail_quote"></span>University of Sofia<br>Faculty of Chemistry<br>Department of Physical Chemistry<br>Quantum and Computational Chemistry Lab<br><br><br><div><span class="gmail_quote">2008/1/28, Mark Abraham <<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Philip Shushkov wrote:<br>> Dear GMX-USERS,<br>><br>> I would like to ask you for some advice because I have serious problem<br>> in setting up a PR MD in GROMACS of DPPC monolayer. Prior to actual MD<br>> simulation I have done exhaustive energy minimization up to 10 kJ/mol*nm<br>
> convergence criterion. Firstly, I minimized the system with restraints<br>> applied to the heavy atoms of the lipid molecules and then the whole<br>> system consisting of 81 DPPC molecules and 5575 water molecules.<br>
> Afterward, in order to avoid some unexpected behavior I minimized the<br>> system with less stringent restraints on the heavy lipid atom which I<br>> would keep in the course of the equilibration, a final step before the<br>
> equilibrating MD run. However, though the three stage minimization<br>> procedure which surely relaxed all bad contacts, the first step of the<br>> MD integration gave high positive energy due to the position restraints<br>
> contribution. Some trial and error steps allowed me to find that the<br>> problem was in the grompp program and may be the parallelization of<br>> GROMACS. The position restraints were OK with *.tpr file constructed<br>
> without the -shuffle option of grompp. Unfortunately, now the LINCS<br>> algorithm crashes and obviously the systems tends to explode. But if I<br>> again turn off the restraints and prepare the *.tpr with -shuffle using<br>
> the same starting configuration as before, everything looks very well.<br>> In conclusion, I would like to humbly ask you where the problem could be<br>> and are there any problems reported concerning PR constraint molecular<br>
> dynamics in parallel mode in GROMACS.<br><br>Position restraints and analytical constraints on the same atoms are<br>likely to lead to numerical problems, particularly if your initial<br>structure (to which you are presumably restraining) doesn't satisfy the<br>
constraints (which get re-inforced every integration step).<br><br>You can eliminate the parallel code as a source of problems by doing a<br>test run on one processor.<br><br>Mark<br>_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at <a href="http://www.gromacs.org/search">http://www.gromacs.org/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/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a><br>
</blockquote></div><br>