zhenlong li wrote:
Dear all,

In my simulation, when the simulation box shrinks, the center of box always shift towards the origin even if the center motion is set to removed every timestep. I just wonder if there is an anticipated center position for a given simulation box in gromacs. The box type is dodecahedron in my simulation.


Yes, this is expected. By default, a unit cell built by Gromacs will be centered at (x/2, y/2, z/2).

My problem is how to make sure a group of atoms, such as proteins, sit at the center of the simulation box after the simulation box become stable. Current, as the center of the box floats, it is very difficult to maintain it by position restraint of the center group.


If you are applying periodic boundary conditions, there is no such thing as a center (since the system is infinite). While it may be a useful convention for visualization purposes, there is no effect on the calculation itself.

-Justin

--
========================================

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

========================================
--
gmx-users mailing list    [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to [email protected].
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to