[gmx-users] Martini simulation problem in recentering trajectory so that the bilayer is at the center
tsjerkw at gmail.com
Fri Nov 6 13:25:32 CET 2009
Hmm, illustrating again that it pays off sometimes to spend a few more
moments thinking about things before writing something. You're right
again. Removel of jumps is performed before centering, but, for the
next frame, the centered system will be used to remove the jumps. So
that's bound to give garbled results for anything but the first frame.
Solution: store the shift that is applied and use it when determining
whether atoms jumped...
On Fri, Nov 6, 2009 at 1:14 PM, XAvier Periole <x.periole at rug.nl> wrote:
> Best Tsjerk,
>>> Well you can see the bilayer as protein complex where each monomer is
>>> a leaflet. Applying -pbc nojump on the trajectory would simply keep the
>>> leaflets together. Then the system can be centered using the bilayer as
>>> reference: this would put the bilayer at the center of the z axis. The
>>> "out of the box" can then be corrected using -pbc mol.
>>> It is trivial is the reference structures are correct.
>> Ah, my imagination got stuck a bit... Yes, that would work. And
>> removal of jumps is applied before centering, so these can be combined
>> in one command. That would make:
>> trjconv -center -pbc nojump
>> trjconv -pbc mol
> Actually I believe the mixing nojump and centering never worked for me!
> I generally do a trial first :))
>> That may be the quickest work-around and should indeed solve the COM
>> trend in the z-direction. But it may give COM jitter or COM trend in
>> the (xy) plane.
>> Tsjerk A. Wassenaar, Ph.D.
>> Junior UD (post-doc)
>> Biomolecular NMR, Bijvoet Center
>> Utrecht University
>> Padualaan 8
>> 3584 CH Utrecht
>> The Netherlands
>> P: +31-30-2539931
>> F: +31-30-2537623
>> 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
> 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 thewww interface
> or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Tsjerk A. Wassenaar, Ph.D.
Junior UD (post-doc)
Biomolecular NMR, Bijvoet Center
3584 CH Utrecht
More information about the gromacs.org_gmx-users