<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'><div style="text-align: left;">Hi,<br><br>This issue is caused by the use of domain decomposition instead<br></div>particle decomposition.<br><br>However, this does not have much to do with the "nojump" behavior<br>you are talking about.<br>The pbcatom is only used for determining the COM of one group.<br>The distance between COM's was always done with full pbc,<br>which means that the vector between the two closest images is used.<br><br>Note that the choice of pbcatom only has an effect for groups with<br>and x, y or z size that is more than half the box length.<br>(and it does not at all affect how the distance between groups<br>is determined)<br><br>Berk<br><br><br><hr id="stopSpelling">&gt; Date: Thu, 30 Oct 2008 13:16:56 -0400<br>&gt; From: chris.neale@utoronto.ca<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: [gmx-users] pull_pbcatom1 and pull code in gromacs 4<br>&gt; <br>&gt; Hello,<br>&gt; <br>&gt; The pull-code .mdp option pull_pbcatom1 is defined in the online manual as:<br>&gt; <br>&gt; "The reference atom for the treatment of periodic boundary conditions <br>&gt; For determining the COM, all atoms are put at their periodic image which <br>&gt; is closest to pull_pbcatom1. A value of 0 means that the middle atom <br>&gt; (number wise) is used."<br>&gt; <br>&gt; I suppose that this is necessary since mdrun now writes broken <br>&gt; molecules. However, this is suboptimal for replica exchange umbrella <br>&gt; sampling (pull = umbrella, pull_geometry = distance) in which one does <br>&gt; not necessarily want to utilize different periodic images of pull_group0 <br>&gt; and pull_group1.<br>&gt; <br>&gt; a) Perhaps I don't understand. Is it possible that although a particular <br>&gt; image of each pulled group is chosen, the pulling is actually still done <br>&gt; along the vector between the closest two images?<br>&gt; <br>&gt; b) If not, is there a simple code modification to utilize a pbc=nojump <br>&gt; type definition here during dynamics?<br>&gt; <br>&gt; Thanks,<br>&gt; Chris.<br>&gt; _______________________________________________<br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://www.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>