[gmx-users] problems when using command "trjconv"
Justin A. Lemkul
jalemkul at vt.edu
Thu Jun 10 16:56:54 CEST 2010
mircial at sjtu.edu.cn wrote:
> Dear All:
>
> I encounter some problems when I am using command "trjconv" to remove
> the periodic boundary condition.
>
> I am doing molecular dynamics simulations of a pentamer, and some of the
> subunits "jump" from one box to its adjacent one. Thus, I try to remove
> the periodic boundary conditions by trjconv command.The commands I used
> are as follows:
>
> trjconv -s trpfile -f xtcfile -o new_xtcfile -n index.ndx -pbc nojump
>
> In the results I found the five subunits are put together to form the
> initial pentamer which is what I expected. However, I also found some of
> the subunits are cut off a part (i.e., the subunit is broken into one
> large piece and several small pieces, the large piece forms pentamer
> with other subunits, and the small pieces distribute in the boundary of
> the box).
>
> Then, I try to solve this problem by using the follow command:
>
> trjconv -s trpfile -f new_xtcfile_lastCommand -o another_new_xtcfile -n
> index.ndx -pbc whole
>
> this time, the subunits are not broken, but the pentamer is not
> formed(i.e., some subunits in one box and some in the adjancent one).
>
> Thus, I am writing to find a way in order to remove the periodic
> boundary condition and put the subunits together as a pentamer with out
> break a seperate subunit into several parts.Could anyone be so kind to
> give me some suggestions? Thank you very much for your insight.
>
Another option is to choose a residue near the central interface of the pentamer
and assign it to an index group, then use trjconv -center (choosing this residue
as the group to center).
-Justin
> Best Wishes
>
> R-X Gu
>
>
>
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
More information about the gromacs.org_gmx-users
mailing list