[gmx-users] direction_periodic
Berk Hess
gmx3 at hotmail.com
Thu Oct 1 15:19:03 CEST 2009
Hi,
Just to be 100% I also ran a simulation with pull rate 0.01 nm/ps.
In this case the pull velocity is smaller than the noise in the COM velocity.
But it is still there, the average velocity is + and - 0.005 nm/ps,
as one would expect. And integrating the velocity gives the poper displacement.
Berk
From: gmx3 at hotmail.com
To: gmx-users at gromacs.org
Subject: RE: [gmx-users] direction_periodic
Date: Thu, 1 Oct 2009 14:30:57 +0200
Hi,
I ran the initial part of the 0.1 nm/ps simulation again (since I deleted the 1 GB trr file).
I get fluctuating velocities for both groups around + and - 0.05 nm/ps with an
(initial) amplitude of the oscillation of 0.02 nm/ps.
This is exactly what you would expect.
Berk
> Date: Thu, 1 Oct 2009 14:08:28 +0200
> From: alexander.herz at mytum.de
> To: gmx-users at gromacs.org
> Subject: Re: [gmx-users] direction_periodic
>
> Hi,
>
> did you look at the actual trajectory (trr/xtc) or the velocity of one
> of the diamond slabs (using g_traj) ?
> I do get a sawtooth style COM distance (pullx.xvg) but both the
> veloc.xvg and the trajectory show that the slabs are not actually moving
> much (the speed of the slab fluctuates wildly around 5e-3nm/ps in a
> 2500ps run when I'm pulling with 0.01nm/ps).
>
> Alex
>
> Berk Hess schrieb:
> > Hi,
> >
> > I just ran a simulation with your input files.
> > I only changed the geometry to direction_periodic and increased the
> > rate to 0.1.
> > The distance decreases linearly with a slight waviness due to the
> > harmonic potential.
> > The location of the reference group shows a sawtooth profile,
> > as one would expect in a periodic system.
> >
> > So it seems to work fine for me.
> >
> > Berk
> >
> > > Date: Thu, 1 Oct 2009 10:24:54 +0200
> > > From: alexander.herz at mytum.de
> > > To: gmx-users at gromacs.org
> > > Subject: Re: [gmx-users] direction_periodic
> > >
> > > Hey,
> > >
> > > did you find the time to look at it?
> > >
> > > Thx,
> > > Alex
> > >
> > > Berk Hess schrieb:
> > > > I am not sure I continued the simulation long enough.
> > > > I'll run it again tomorrow.
> > > >
> > > > Berk
> > > >
> > > > > Date: Mon, 28 Sep 2009 16:57:19 +0200
> > > > > From: alexander.herz at mytum.de
> > > > > To: gmx-users at gromacs.org
> > > > > Subject: Re: [gmx-users] direction_periodic
> > > > >
> > > > > Well..they do move forward for a few ps and then they move back
> > (they
> > > > > oscillate).
> > > > > Did you check for a longer period?
> > > > >
> > > > > Alex
> > > > >
> > > > > Berk Hess schrieb:
> > > > > > Strange...
> > > > > > I tested the code on the files you sent me
> > > > > > and the slabs were moving smoothly.
> > > > > >
> > > > > > Berk
> > > > > >
> > > > > > > Date: Mon, 28 Sep 2009 15:40:16 +0200
> > > > > > > From: alexander.herz at mytum.de
> > > > > > > To: gmx-users at gromacs.org
> > > > > > > Subject: Re: [gmx-users] direction_periodic
> > > > > > >
> > > > > > > Hey,
> > > > > > >
> > > > > > > so using the new direction_periodic with some kernels that
> > > > actually work
> > > > > > > and the sources from the git (unmodified)
> > > > > > > gives a nice smooth pull force but no actual mean displacement
> > > > of the
> > > > > > > groups to be pulled (on average they just stay where they
> > are,see
> > > > > > > attached graph). I'm still trying to pull them into opposite
> > > > directions
> > > > > > > though.
> > > > > > >
> > > > > > > I'm using this setup:
> > > > > > > pull = umbrella
> > > > > > > pull_geometry = direction_periodic
> > > > > > > pull_ngroups = 1
> > > > > > > pull_group0 = DIAM
> > > > > > > pull_group1 = DIAM2
> > > > > > > pulldim = Y N N
> > > > > > > pull_k1 = 1000.0
> > > > > > > pull_rate1 = 0.01
> > > > > > > pull_vec1 = -1.0 0.0 0.0
> > > > > > >
> > > > > > > How can I make the two diamond slabs actually move in opposite
> > > > > > > directions rather than just oscillating around their
> > > > > > > initial positions?
> > > > > > >
> > > > > > > Thx,
> > > > > > > Alex
> > > > > > >
> > > > > > > Berk Hess schrieb:
> > > > > > > > Hi,
> > > > > > > >
> > > > > > > > Don't use fortran.
> > > > > > > > We will get rid of it before the 4.1 release.
> > > > > > > >
> > > > > > > > Berk
> > > > > > > >
> > > > > > > > > Date: Mon, 21 Sep 2009 16:17:17 +0200
> > > > > > > > > From: alexander.herz at mytum.de
> > > > > > > > > To: gmx-users at gromacs.org
> > > > > > > > > Subject: Re: [gmx-users] direction_periodic
> > > > > > > > >
> > > > > > > > > Hey,
> > > > > > > > >
> > > > > > > > > ok, I checked this. It works ok with the assembly
> > kernels on my
> > > > > > local
> > > > > > > > > machine, but it fails with fortran kernels
> > > > > > > > > on the HPC and on my local machine.
> > > > > > > > >
> > > > > > > > > Alex
> > > > > > > > >
> > > > > > > > > Berk Hess schrieb:
> > > > > > > > > > Yes.
> > > > > > > > > >
> > > > > > > > > > Berk
> > > > > > > > > >
> > > > > > > > > > > Date: Mon, 21 Sep 2009 14:42:43 +0200
> > > > > > > > > > > From: alexander.herz at mytum.de
> > > > > > > > > > > To: gmx-users at gromacs.org
> > > > > > > > > > > Subject: Re: [gmx-users] direction_periodic
> > > > > > > > > > >
> > > > > > > > > > > With master branch you mean the code I get via
> > > > > > > > > > >
> > > > > > > > > > > |git clone git://git.gromacs.org/gromacs.git
> > > > > > > > > > >
> > > > > > > > > > > right?
> > > > > > > > > > >
> > > > > > > > > > > Alex
> > > > > > > > > > > |
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > Berk Hess schrieb:
> > > > > > > > > > > > I tested on your system and got good results.
> > > > > > > > > > > > There is still a tricky issue: the pull COM is still
> > > > > > determined
> > > > > > > > > > > > in the "standard" way by summing distances from the
> > > > pbcatom.
> > > > > > > > > > > > Therefore atoms should not change nearest image
> > from he
> > > > > > pbcatom.
> > > > > > > > > > > > This would result in nasty noise in the pull COM and
> > > > force.
> > > > > > > > > > > >
> > > > > > > > > > > > I would suggest that you try to run the master
> > branch and
> > > > > > check
> > > > > > > > > > > > if that works.
> > > > > > > > > > > >
> > > > > > > > > > > > Berk
> > > > > > > > > > > >
> > > > > > > > > > > > > Date: Mon, 21 Sep 2009 14:31:34 +0200
> > > > > > > > > > > > > From: alexander.herz at mytum.de
> > > > > > > > > > > > > To: gmx-users at gromacs.org
> > > > > > > > > > > > > Subject: [gmx-users] direction_periodic
> > > > > > > > > > > > >
> > > > > > > > > > > > > Hey,
> > > > > > > > > > > > >
> > > > > > > > > > > > > after some pain of merging the dev branch into
> > our 4.0.5
> > > > > > version
> > > > > > > > > > I got
> > > > > > > > > > > > > the new pull mode "direction_periodic"
> > > > > > > > > > > > > running over the weekend. There's some weird
> > > > rotation of the
> > > > > > > > pulled
> > > > > > > > > > > > > objects going on and pbc seem weird as well (there
> > > > are water
> > > > > > > > > > molecules
> > > > > > > > > > > > > in those positions where I'd expect the periodic
> > > > image of my
> > > > > > > > diamond
> > > > > > > > > > > > > slab which leaves the box at one side). I guess you
> > > > > > tested the
> > > > > > > > > > new pull
> > > > > > > > > > > > > mode somehow, so any ideas what's going on here? I'm
> > > > still
> > > > > > > > trying to
> > > > > > > > > > > > > perform the same experiment for which I send you the
> > > > > > input files
> > > > > > > > > > while
> > > > > > > > > > > > > ago and for which you kindly implemented the new
> > > > pull mode.
> > > > > > > > > > > > >
> > > > > > > > > > > > > Thx for your help,
> > > > > > > > > > > > > Alex
> > > > > > > > > > > > >
> > > > > > > > > > > > > Berk Hess schrieb:
> > > > > > > > > > > > > > I have committed a new pull geometry
> > > > direction_periodic to
> > > > > > > > the git
> > > > > > > > > > > > > > master branch. It is not documented yet.
> > > > > > > > > > > > > > It works the same at direction, but allows
> > > > distances to be
> > > > > > > > larger
> > > > > > > > > > > > > > than half the box and does not add the pull force
> > > > to the
> > > > > > > > virial.
> > > > > > > > > > > > > >
> > > > > > > > > > > > > > Berk
> > > > > > > > > > > > > _______________________________________________
> > > > > > > > > > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > > > > > > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > > > > > > > > > Please search the archive at
> > > > http://www.gromacs.org/search
> > > > > > > > before
> > > > > > > > > > > > posting!
> > > > > > > > > > > > > Please don't post (un)subscribe requests to the
> > list.
> > > > > > Use the
> > > > > > > > > > > > > www interface or send it to
> > > > gmx-users-request at gromacs.org.
> > > > > > > > > > > > > Can't post? Read
> > > > > > http://www.gromacs.org/mailing_lists/users.php
> > > > > > > > > > > >
> > > > > > > > > > > >
> > > > > > > > > >
> > > > > > > >
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > > > > > > > > Express yourself instantly with MSN Messenger! MSN
> > > > Messenger
> > > > > > > > > > > >
> > <http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
> > > > > > > > > > > >
> > > > > > > > > >
> > > > > > > >
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > > > > > > > >
> > > > > > > > > > > > _______________________________________________
> > > > > > > > > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > > > > > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > > > > > > > > Please search the archive at
> > http://www.gromacs.org/search
> > > > > > before
> > > > > > > > > > posting!
> > > > > > > > > > > > Please don't post (un)subscribe requests to the list.
> > > > Use the
> > > > > > > > > > > > www interface or send it to
> > gmx-users-request at gromacs.org.
> > > > > > > > > > > > Can't post? Read
> > > > > > http://www.gromacs.org/mailing_lists/users.php
> > > > > > > > > > >
> > > > > > > > > > > _______________________________________________
> > > > > > > > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > > > > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > > > > > > > Please search the archive at
> > http://www.gromacs.org/search
> > > > > > before
> > > > > > > > > > posting!
> > > > > > > > > > > Please don't post (un)subscribe requests to the list.
> > > > Use the
> > > > > > > > > > > www interface or send it to
> > gmx-users-request at gromacs.org.
> > > > > > > > > > > Can't post? Read
> > > > http://www.gromacs.org/mailing_lists/users.php
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > >
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > > > > > > Express yourself instantly with MSN Messenger! MSN
> > Messenger
> > > > > > > > > > <http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
> > > > > > > > > >
> > > > > > > >
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > > > > > >
> > > > > > > > > > _______________________________________________
> > > > > > > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > > > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > > > > > > Please search the archive at http://www.gromacs.org/search
> > > > before
> > > > > > > > posting!
> > > > > > > > > > Please don't post (un)subscribe requests to the list.
> > Use the
> > > > > > > > > > www interface or send it to gmx-users-request at gromacs.org.
> > > > > > > > > > Can't post? Read
> > > > http://www.gromacs.org/mailing_lists/users.php
> > > > > > > > >
> > > > > > > > > _______________________________________________
> > > > > > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > > > > > Please search the archive at http://www.gromacs.org/search
> > > > before
> > > > > > > > posting!
> > > > > > > > > Please don't post (un)subscribe requests to the list.
> > Use the
> > > > > > > > > www interface or send it to gmx-users-request at gromacs.org.
> > > > > > > > > Can't post? Read
> > http://www.gromacs.org/mailing_lists/users.php
> > > > > > > >
> > > > > > > >
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > > > > Express yourself instantly with MSN Messenger! MSN Messenger
> > > > > > > > <http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
> > > > > > > >
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > > > >
> > > > > > > > _______________________________________________
> > > > > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > > > > Please search the archive at http://www.gromacs.org/search
> > before
> > > > > > posting!
> > > > > > > > Please don't post (un)subscribe requests to the list. Use the
> > > > > > > > www interface or send it to gmx-users-request at gromacs.org.
> > > > > > > > Can't post? Read
> > http://www.gromacs.org/mailing_lists/users.php
> > > > > > >
> > > > > >
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > > Express yourself instantly with MSN Messenger! MSN Messenger
> > > > > > <http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
> > > > > >
> > > >
> > ------------------------------------------------------------------------
> > > > > >
> > > > > > _______________________________________________
> > > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > > Please search the archive at http://www.gromacs.org/search before
> > > > posting!
> > > > > > Please don't post (un)subscribe requests to the list. Use the
> > > > > > www interface or send it to gmx-users-request at gromacs.org.
> > > > > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> > > > >
> > > > > _______________________________________________
> > > > > gmx-users mailing list gmx-users at gromacs.org
> > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > > Please search the archive at http://www.gromacs.org/search before
> > > > posting!
> > > > > Please don't post (un)subscribe requests to the list. Use the
> > > > > www interface or send it to gmx-users-request at gromacs.org.
> > > > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> > > >
> > > >
> > ------------------------------------------------------------------------
> > > > Express yourself instantly with MSN Messenger! MSN Messenger
> > > > <http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
> > > >
> > ------------------------------------------------------------------------
> > > >
> > > > _______________________________________________
> > > > gmx-users mailing list gmx-users at gromacs.org
> > > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > > Please search the archive at http://www.gromacs.org/search before
> > posting!
> > > > Please don't post (un)subscribe requests to the list. Use the
> > > > www interface or send it to gmx-users-request at gromacs.org.
> > > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> > >
> > > _______________________________________________
> > > gmx-users mailing list gmx-users at gromacs.org
> > > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > > Please search the archive at http://www.gromacs.org/search before
> > posting!
> > > Please don't post (un)subscribe requests to the list. Use the
> > > www interface or send it to gmx-users-request at gromacs.org.
> > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> >
> > ------------------------------------------------------------------------
> > See all the ways you can stay connected to friends and family
> > <http://www.microsoft.com/windows/windowslive/default.aspx>
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > gmx-users mailing list gmx-users at gromacs.org
> > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > Please search the archive at http://www.gromacs.org/search before posting!
> > Please don't post (un)subscribe requests to the list. Use the
> > www interface or send it to gmx-users-request at gromacs.org.
> > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
> _______________________________________________
> gmx-users mailing list gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Express yourself instantly with MSN Messenger! MSN Messenger
_________________________________________________________________
What can you do with the new Windows Live? Find out
http://www.microsoft.com/windows/windowslive/default.aspx
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20091001/f1235e67/attachment.html>
More information about the gromacs.org_gmx-users
mailing list