[gmx-users] pull_geometry=direction

[Lalita Shaki]

Dear gromacs users,

I have been working on a pulling simulation of the unfolding of a protein
and I have found out that went using pull_geometry=direction the code
checks at each time step if the groups distance is  higher than 0.49 each
of the dimensions of the box. If the distance is higher that 0.49 of at
least one of the dimensions the simulations ends in an error. So pretty
much when you want to unfold a protein you would have to use a really large
cubic box.

My question is if someone knows why gromacs check the magnitude of the
distance again the dimension of the box and does not check each of the
components of the distance vector to it correspondent dimension of the box.
I think this will also avoid the problem of miss calculating the
 pull-reference distance and at the same time will permit to run
calculations with a rectangular box.


Thank you

lalita

-- 
Lalita S. Uribe.
European Master in Theoretical Chemistry and Computational Modeling.
PhD student. Johannes Gutenberg-Universität Mainz