[gmx-users] yet another question about PBC...

[Justin A. Lemkul]

Ramon Crehuet wrote:
> Hi,
> 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."

-Justin


> Thanks in advance,
> Ramon
>
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.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-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
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================