<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Hi,<br><br>I hope you are using 4.0.5, I fixed several bug in the pull code for older 4.0 versions.<br><br>The problems you are seing could be due to&nbsp; pbc.<br>Do you have pbc in Z, and what is the height of your box?<br><br>A safer setup is:<br>pull_geometry = direction<br>pull_vec1 = 0 0 1<br>But they should give the same answers if you do not have pbc issues.<br><br>Berk<br><br>&gt; Date: Fri, 31 Jul 2009 12:13:07 +0200<br>&gt; From: alexander.herz@mytum.de<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: [gmx-users] pulling<br>&gt; <br>&gt; Hey,<br>&gt; <br>&gt; I appear to have serious trouble understanding how to set up the pulling<br>&gt; properly.<br>&gt; <br>&gt; I have many configurations of a protein partially adsorbed to a froozen<br>&gt; surface (the configs differ<br>&gt; in the amount of the protein that has been desorbed).<br>&gt; Now I want the pulling to keep the distance of the desorbed end of the<br>&gt; protein to the surface using the harmonic pot.<br>&gt; Now the documentation is not very clear how this all works so I ran<br>&gt; several experiments to figure it out but I failed.<br>&gt; I use the following options:<br>&gt; <br>&gt; ;PULLING<br>&gt; pull                     = umbrella<br>&gt; pull_geometry            = distance<br>&gt; pull_dim                 = N N Y<br>&gt; pull_nstxout             = 1000<br>&gt; pull_nstfout             = 1000<br>&gt; pull_ngroups             = 1<br>&gt; pull_group0              = GLD<br>&gt; pull_group1              = ASN<br>&gt; pull_vec1                = 0.0 0.0 0.0<br>&gt; pull_init1               =   5.27778<br>&gt; pull_rate1               = 0.0<br>&gt; pull_k1                  = 100<br>&gt; <br>&gt; <br>&gt; where gld is the surface and asn is the end residue of the protein and<br>&gt; pull_init1 is set to the desired COM distance of the two<br>&gt; groups (gld is froozen). I use the same settings for all runs, only<br>&gt; changing pull_init1 to get the desired distance.<br>&gt; Now for some reason using this setup either pulls the ASN end of the<br>&gt; protein completely onto the surface or very far away from it depending<br>&gt; on the value I use for pull_init1.<br>&gt; So the distance between what and what shall I put for pull_init1? What<br>&gt; else is wrong?<br>&gt; <br>&gt; Thx,<br>&gt; Alex<br>&gt; _______________________________________________<br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>