<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Correct.<br>Make sure than your section A spreads over 90% (or 95%) of the x-range.<br>Then you are sure that the distance in x between the pbc-atom and every<br>other atom in A is not more than 0.45*boxlength plus a bit of fluctation.<br><br>Berk<br><br>> Date: Fri, 2 Oct 2009 11:49:37 +0200<br>> From: alexander.herz@mytum.de<br>> To: gmx-users@gromacs.org<br>> Subject: Re: [gmx-users] direction_periodic<br>> <br>> Hi,<br>> <br>> thx for suggesting the workaround, just to make sure I understand the<br>> workaround correctly:<br>> <br>> I devide each slab in 2 sections, one section contains 90% of the atoms<br>> (section A), the other the rest (section B).<br>> Now I set the pbc_atom to the center atom of section A and pull it using<br>> constraints?<br>> <br>> Alex<br>> <br>> Berk Hess schrieb:<br>> > Hi,<br>> ><br>> > This is due to the periodicity issue of the COM calculation that I<br>> > have mentioned before.<br>> > The problem is that some atoms are at close to half a box length from<br>> > the pbc_atom<br>> > and those can jump back and forth between two periodic images in the<br>> > distance calculations.<br>> > This "noise" will make the constraint code crash, since it tries to<br>> > correct for it immediately.<br>> > The umbrella potential smooths this effect out somewhat and can<br>> > therefore run stably,<br>> > but you would still see the noise in the results.<br>> ><br>> > Fixing this is difficult, since it would require storing the periodic<br>> > image of each atom<br>> > in the pull group at the start of the simulation and determining the<br>> > image for the COM<br>> > calculation using the initial image. Also the checkpoint file should<br>> > store the images.<br>> > I would not like to put such code in the standard version.<br>> ><br>> > A workaround could be to pull on 90% of the slabs, leaving a part of<br>> > 10% along x out.<br>> > Then define the pbc_atom in the middle of the 90%. The 10% will just<br>> > move along.<br>> > This should solve your problem and should not produce artifacts.<br>> ><br>> > Berk<br>> ><br>> > > Date: Fri, 2 Oct 2009 10:25:21 +0200<br>> > > From: alexander.herz@mytum.de<br>> > > To: gmx-users@gromacs.org<br>> > > Subject: Re: [gmx-users] direction_periodic<br>> > ><br>> > > You did suggest that, but it crashes at step 0, as said :)<br>> > ><br>> > > Alex<br>> > ><br>> > > Berk Hess schrieb:<br>> > > > Hi,<br>> > > ><br>> > > > Then you should use constraint pulling instead of an umbrella<br>> > potential.<br>> > > > I have been wondering all the time why you are using a potential<br>> > > > and not a constraint. I think I had suggested using a constraint some<br>> > > > time ago, but I am not sure that I did.<br>> > > ><br>> > > > Berk<br>> > > ><br>> > > > > Date: Fri, 2 Oct 2009 09:45:42 +0200<br>> > > > > From: alexander.herz@mytum.de<br>> > > > > To: gmx-users@gromacs.org<br>> > > > > Subject: Re: [gmx-users] direction_periodic<br>> > > > ><br>> > > > > Hey,<br>> > > > ><br>> > > > > I want to pull the diamonds with more or less exactly the speed<br>> > I give<br>> > > > > for the pulling without the massive noise<br>> > > > > (in gromacs 3.3 this worked using afm type pulling, after some<br>> > > > > equilibration the slabs settled to an almost constant<br>> > > > > pull speed without soo much fluctuation using the same<br>> > parameters I use<br>> > > > > noew for umbrella).<br>> > > > > So I tried using a constraint rather than umbrella, which alwas<br>> > crashes<br>> > > > > at the first step (even if I did some equilibration<br>> > > > > using umbrella pulling beforehands).<br>> > > > ><br>> > > > > So how can I achieve a pull speed which is about constant 0.01<br>> > nm/ps and<br>> > > > > with fluctuations less than 10% of this speed?<br>> > > > ><br>> > > > > Thx,<br>> > > > > Alex<br>> > > > > _______________________________________________<br>> > > > > gmx-users mailing list gmx-users@gromacs.org<br>> > > > > http://lists.gromacs.org/mailman/listinfo/gmx-users<br>> > > > > Please search the archive at http://www.gromacs.org/search before<br>> > > > posting!<br>> > > > > Please don't post (un)subscribe requests to the list. Use the<br>> > > > > www interface or send it to gmx-users-request@gromacs.org.<br>> > > > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>> > > ><br>> > > ><br>> > ------------------------------------------------------------------------<br>> > > > See all the ways you can stay connected to friends and family<br>> > > > <http://www.microsoft.com/windows/windowslive/default.aspx><br>> > > ><br>> > ------------------------------------------------------------------------<br>> > > ><br>> > > > _______________________________________________<br>> > > > gmx-users mailing list gmx-users@gromacs.org<br>> > > > http://lists.gromacs.org/mailman/listinfo/gmx-users<br>> > > > Please search the archive at http://www.gromacs.org/search before<br>> > posting!<br>> > > > Please don't post (un)subscribe requests to the list. Use the<br>> > > > www interface or send it to gmx-users-request@gromacs.org.<br>> > > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>> > ><br>> > > _______________________________________________<br>> > > gmx-users mailing list gmx-users@gromacs.org<br>> > > http://lists.gromacs.org/mailman/listinfo/gmx-users<br>> > > Please search the archive at http://www.gromacs.org/search before<br>> > posting!<br>> > > Please don't post (un)subscribe requests to the list. Use the<br>> > > www interface or send it to gmx-users-request@gromacs.org.<br>> > > Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>> ><br>> > ------------------------------------------------------------------------<br>> > What can you do with the new Windows Live? Find out<br>> > <http://www.microsoft.com/windows/windowslive/default.aspx><br>> > ------------------------------------------------------------------------<br>> ><br>> > _______________________________________________<br>> > gmx-users mailing list gmx-users@gromacs.org<br>> > http://lists.gromacs.org/mailman/listinfo/gmx-users<br>> > Please search the archive at http://www.gromacs.org/search before posting!<br>> > Please don't post (un)subscribe requests to the list. Use the <br>> > www interface or send it to gmx-users-request@gromacs.org.<br>> > Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>> <br>> _______________________________________________<br>> gmx-users mailing list gmx-users@gromacs.org<br>> http://lists.gromacs.org/mailman/listinfo/gmx-users<br>> Please search the archive at http://www.gromacs.org/search before posting!<br>> Please don't post (un)subscribe requests to the list. Use the <br>> www interface or send it to gmx-users-request@gromacs.org.<br>> Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>                                            <br /><hr />See all the ways you can stay connected <a href='http://www.microsoft.com/windows/windowslive/default.aspx' target='_new'>to friends and family</a></body>
</html>