[gmx-users] trjconv -pbc

David spoel at xray.bmc.uu.se
Mon Mar 10 20:15:17 CET 2003


On Mon, 2003-03-10 at 17:19, Andre Farias de Moura wrote:
> 
> hello david,
> 
> the genconf program worked just fine and I could view my
> micelles whole again, thank you.
> 
> nonetheless, I've got another related question: the analyses
> programs (g_gyrate, for instance) consider those molecules
> in the central simulation box or do they apply periodic
> boundary conditions and the minimum image convention?
> 
> I'm not pretty sure, but it seems to me that it is not
> enough to apply PBC and minimum image convention to get
> the micelle whole again (this procedure would locate the
> images nearer the simulation box origin), and it's not
> correct to use those molecules in the central box either.
> in both cases, the micellar center of mass and radius of
> gyration would not be correct, they would change every
> time a lipid molecule jumps across the boundary. is that
> right?
> 
> I think these programs should find out which molecules are
> actually in the aggregate before any structural calculation
> is performed... I just don't know how to do it yet.

I have implemented such an option in trjconv:
trjconv -pbc cluster

it is in CVS, but you want to make a separate tree for this, since mdrun
can not be trusted, use it only for trjconv and maybe some analysis
tools (g_cluster, g_clustsize)
> 
> thanks in advance,
> 
> andre
> 
> 
> On 9 Mar 2003, David wrote:
> 
> > On Sun, 2003-03-09 at 09:58, Andre Farias de Moura wrote:
> > >
> > > hello rainer,
> > >
> > > my tpr-file is ok, lipid molecules are not broken after
> > > "-pbc whole" option is applied. any other clue?
> >
> > Since you have a pbc system, lipids will have to be either on one side
> > of the box, or on the other. If they move they may be transported to the
> > other side. What people usually do for animations of bilayers is to
> > extend the pbc (using e.g. genconf)
> >
> >
> > --
> > Groeten, David.
> > ________________________________________________________________________
> > Dr. David van der Spoel, 	Dept. of Cell and Molecular Biology
> > Husargatan 3, Box 596,  	75124 Uppsala, Sweden
> > phone:	46 18 471 4205		fax: 46 18 511 755
> > spoel at xray.bmc.uu.se	spoel at gromacs.org   http://xray.bmc.uu.se/~spoel
> > ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
> > _______________________________________________
> > gmx-users mailing list
> > gmx-users at gromacs.org
> > http://www.gromacs.org/mailman/listinfo/gmx-users
> > Please don't post (un)subscribe requests to the list. Use the
> > www interface or send it to gmx-users-request at gromacs.org.
> >
> 
> _______________________________________________
> gmx-users mailing list
> gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-request at gromacs.org.
-- 
Groeten, David.
________________________________________________________________________
Dr. David van der Spoel, 	Dept. of Cell and Molecular Biology
Husargatan 3, Box 596,  	75124 Uppsala, Sweden
phone:	46 18 471 4205		fax: 46 18 511 755
spoel at xray.bmc.uu.se	spoel at gromacs.org   http://xray.bmc.uu.se/~spoel
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++



More information about the gromacs.org_gmx-users mailing list