Hi all,<div>I&#39;d like to calculate the PMF for a small molecule moving from a membrane-bound protein into the membrane itself.  I tried using the &quot;cylinder&quot; option for pull_geometry, but when I tried to analyze the results with g_wham it said that &quot;cylinder&quot; isn&#39;t a supported option.  Is this correct?  Is there anything I can do to still use this geometry?  </div>

<div><br></div><div>If it won&#39;t work, I thought maybe I could instead use &quot;distance&quot; for pull_geometry, and then restrict pull_dim to the z axis (N N Y, right?).  I&#39;m worried that if I do this and the centers of mass of my molecule and the bilayer don&#39;t align, then the umbrella potential will still be directed along the vector connecting the two COMs and the Z-component of the potential will therefore be something smaller (dependent on the angle between the COM-COM vector and the z axis).  Does anyone know how this is implemented?  </div>

<div><br></div><div>Thanks,</div><div>Gard</div>