[gmx-users] Dimer g_rms

[XAvier Periole]

On Jun 30, 2010, at 10:34 AM, Carla Jamous wrote:

> Dear all,
>
> I'm running my first simulation of a dimer.
> When I run g_rms on Calpha (lsq fit and RMSD calculation on Calpha)  
> to get the RMSD of the whole dimer, the graph is ascending, till  
> reaching a value of 8 A.
> in this case, I take a .tpr file, a .xtc file and a .ndx file of the  
> whole dimer.
>
> However, when I run g_rmsd on Calpha, but this time, to get the RMSD  
> of one chain (one monomer), the graph is constant with value around  
> 1 A. and this for both chains, when taken alone.
> In this case, I take a .tpr file, a .xtc file and a .ndx file of the  
> monomer.
>
> So is this an error of using g_rms?
Well this is not an error. Evidently your monomers are stable  
(rmsd~0.1nm)
and you dimer is not.
This might be result of either an actual deformation of the dimer or
one of your monomer is crossing the periodic boundaries and thereby
the rmsd is not representative.
try to trjconv your trajectory using the -pbc nojump and run g_rms  
again.
>
> Thank you.
>
> Carla
> -- 
> gmx-users mailing list    gmx-users at 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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php