[gmx-users] Re: Non-equilibraium MD: simulating a moving
Vitaly V. Chaban
vvchaban at gmail.com
Mon Mar 23 13:10:15 CET 2009
Hi and thank you! This cleared the situation. I didn't care about comm-mode
but only set comm-grps...
And let me one more question. Is there any possibility to define a constant
velocity of my rigid wall, not an acceleration?
>Hi, have you set comm-mode = None? If no, you can set it and try again.
>I want to simulate a rigid plate (wall) moving along one direction and
>reflecting all the particles it meets.
>I have made a potential only for repulsion which acts at the distances
>less than 1 nm (and equals zero at bigger ones) and bound it to a
>particle. If the plate is immovable everything goes OK. To make it go
>acc-grps = WAL
>accelerate = 0 0 1
>where the acceleration is in [nm/ps2] (from the gmx manual 4.0).
>To make a plate motion insensitive to the molecules it reflects, its
>mass was set to a very big value. But during the following run the
>plane moved very slowly while all other particles in the system have
>obtained an acceleration at the same direction. Even those ones which
>are to far to interact with a plate. Why? Then if the mass of the
>plate is reduced the above effect is not so pronounced but then the
>plate becomes very sensible to the interactions, so its motion is not
>as we want...
>If the acceleration is given in [nm/ps2], how the mass can influence
>the motion? I thought this value is just added to the value obtained
>from the interations, but it's evidently not true here.
>Please let me know if it's possible in gromacs to make such moving
>plate I described and if so, what is a good practice to simulate this
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-users