[gmx-users] yet another question about PBC...
Justin A. Lemkul
jalemkul at vt.edu
Mon Jun 30 17:26:52 CEST 2008
Ramon Crehuet wrote:
> Please don't flame me for asking yet another question about the "out of
> the box" and PBC issue. I've tried to RTFM, RTFML, but could not find
> the question I'm asking.
> I (think) I understand what PBC implies and the fact that "out of the
> box" makes no sense because the system is periodic and gromacs treats
> the periodicity correctly. However, when calculating the rmsd of a
> trajectory with g_rms with respect to the protein I get a warning:
> "cannot make broken molecules whole without a run input file, don't
> worry mdrun doesn't write broken molecules".
> and the rmsd has unrealistic jumps when the chains of the protein (a
> tetramer) leave the box (but if I calculate distances with g_bond, they
> are always ok)
> If I correct the trajectory with trjconf -pbc nojump, the new trajectory
> seems to behave as expected.
> So my questions are, i)should I worry about the warning? ii) why does
> g_rms not treat PBC as expected? Should I use some options?
As far as the warning, it's due to the fact that you specified something
like a .pdb or .gro file for the -s flag, and not the .tpr file. Not so
important in this case, but it would affect, for example, making
molecules whole with trjconv.
As for the second part, perhaps someone else has better insight, but
usually the quick answer is to compensate for PBC using trjconv, as you
have. I know there are a few messages about this in the list archives,
but I think they usually end with "use trjconv."
> Thanks in advance,
> gmx-users mailing list gmx-users at gromacs.org
> 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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
jalemkul[at]vt.edu | (540) 231-9080
More information about the gromacs.org_gmx-users