Hi, I am trying to look at the free energy differences as a function of end-to-end distance of peptides in water. I am running a set of simulations with the two atoms of interest joined by a type 2 constraint, with the length varying from 0.5 to 3 nm. However, these simulations all explode, even though I am starting with conformations from an unconstrained simulation that are taken through another round of energy minimization and relaxtion before running the constrained simulations.
There is an mdrun option, -rcon, that I have not been setting. I read this as the distance over which atoms interacting through the constraint will be communicated between processors for the P-LINC algorithm. since some of the end-to-end distances can be relatively large, I was wondering if, in principle, I should be setting -rcon to the constarint distance +1 nm to ensure that the constraint is not randomly disappearing from consideration when the two ends of the peptide are in different domains. Thanks for any insights. Warren Gallin-- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php