<div dir="ltr">Hello,<br><br>I have done umbrella sampling with &quot;pull=umbrella&quot; and I found that the pulling group has high fluctuations and sometimes moving out of the periodic box. I think that the harmonic potential is not properly applied and thus the pulling group is not retained with the specified COM distance between the reference and pulling group. Can we use pull=constraint option to retain the pulling group within the COM distance? my pulling code is as bellow with restrain at 0.5nm distance from ref group. I just want to get some idea about this pull code modification.  <br>

<br>pull            = constraint<br>pull_geometry    = distance<br>pull_dim        = N N Y<br>pull_start      = no <br>pull_ngroups    = 1<br>pull_group0     = Chain_A <br>pull_group1     = Chain_B <br>pull_init1      = 0.50        <br>

pull_rate1      = 0.0<br>pull_k1         = 1000      <br>pull_nstxout    = 500       <br>pull_nstfout    = 500       <br><br><br>Regards,<br>Vijay.<br>
</div>