[gmx-users] PBC correction to visualize a protein-membrane structure
Juan Munoz-Garcia
juan.munoz-garcia at bioch.ox.ac.uk
Fri May 16 18:58:18 CEST 2014
Dear Justin,
thank you. I’ve tried the following but neither of them worked, I get the same result.
trjconv -f input.gro -o output.gro -s .tpr -trans 0 0 z_box/2 -pbc mol -ur compact
trjconv -f input.gro -o output.gro -s tpr -trans x_box/2 y_box/2 z_box/2 -pbc mol -ur compact
This is the last line of my input.gro file
19.51671 19.51671 13.80040 0.00000 0.00000 0.00000 0.00000 9.75836 9.75836
Juan C.
Your system was probably centered at z=0 instead of z = z_box/2 (commands used
for building the system would also help - I guess I should have been more
clear), hence it's getting split across PBC. trjconv -trans in conjunction with
-pbc mol should be what you need.
-Justin
--
Begin forwarded message:
From: Juan Carlos Munoz Garcia <juan.munoz-garcia at bioch.ox.ac.uk<mailto:juan.munoz-garcia at bioch.ox.ac.uk>>
Subject: Re: PBC correction to visualize a protein-membrane structure
Date: 16 May 2014 14:08:50 BST
To: <gromacs.org_gmx-users at maillist.sys.kth.se<mailto:gromacs.org_gmx-users at maillist.sys.kth.se>>
Thank you Justin,
please find a dropbox link to the image below.
I’ve used
trjconv_mpi -f NPT.trr -o NPT_2.trr -s NPT.tpr -pbc whole
trjconv_mpi -f NPT_2.trr -o NPT_PBC.trr -s NPT.tpr -pbc mol -ur compact -center
and different combinations of those
https://www.dropbox.com/s/az1gvb299y6wh8h/Screen%20Shot%202014-05-06%20at%2021.15.07.png
Thank you.
Regards.
Juan C.
More information about the gromacs.org_gmx-users
mailing list