<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta content="MSHTML 6.00.6000.17023" name="GENERATOR">
<style></style><style id="owaParaStyle">P {
        MARGIN-TOP: 0px; MARGIN-BOTTOM: 0px
}
</style>
</head>
<body ocsi="0">
<div style="FONT-SIZE: 13px; COLOR: #000000; DIRECTION: ltr; FONT-FAMILY: Tahoma">
<div>Hi again,</div>
<div>I thought your suggestions would work for my membrane, but it seems like the removal of COM translation of the bilayer and water separately does not stop the system from translating in the box. My new simulations are now soon crashing again since the lipids&nbsp;are
 crossing the pbc. I was using two comm-grps, one containing DPPC_CHOL_MOL and one with water. How is it possible that the system still is able to translate? Is there any other way to do this? Otherwise I have to manually translate the system back to its original
 position in the box after the simulation has crashed and then continue the simulation, but this is not very good since I cannot use the checkpoint file then and the continuation is not exact. Any suggestion what to do? Thanks in advance.
</div>
<div>&nbsp;</div>
<div>Emma</div>
<div><br>
<div>&nbsp;</div>
</div>
<div style="FONT-SIZE: 16px; COLOR: #000000; FONT-FAMILY: Times New Roman">
<div id="divRpF268400" style="DIRECTION: ltr">
<hr tabindex="-1">
<font face="Tahoma" color="#000000" size="2"><b>Från:</b> gmx-users-bounces@gromacs.org [gmx-users-bounces@gromacs.org] för XAvier Periole [x.periole@rug.nl]<br>
<b>Skickat:</b> den 27 april 2010 16:32<br>
<b>Till:</b> Discussion list for GROMACS users<br>
<b>Ämne:</b> Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box<br>
</font><br>
</div>
<div></div>
<div>
<div><br>
</div>
you should remove the water and lipid bilayer COM separately.&nbsp;
<div>I am not sure what you should do with your small molecule though.&nbsp;</div>
<div>Probably best to add it to the bilayer as you constrain their relative&nbsp;</div>
<div>position!</div>
<div><br>
</div>
<div>
<div>
<div>On Apr 27, 2010, at 4:32 PM, ERIKSSON, EMMA wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite"><span class="Apple-style-span">
<div>
<div style="FONT-SIZE: 13px; COLOR: rgb(0,0,0); DIRECTION: ltr; FONT-FAMILY: Tahoma">
<div>Hi again,</div>
<div>Thanks Xavier for your reply. I didn't know that this mdp option existed. However, I read the manual and also checked the mdout.mdp files for my previous simulations, and I understood it as if those are the default settings even if you don't specify any
 of them in the md.mdp file. The default comm_groups is the whole system so I guess if I'm not writing anything there it will take the whole system. In that case I think that in my previous simulations the translational motion should have been removed for the
 whole system, but since it's obviously not remove something is wrong. Or did I misunderstand everything? My system consists of DPPC lipids, cholesterol, water and one small molecule. Should I specify comm_groups as only the lipids? In that case I get a warning
 from grompp. You wrote comm-grps = membrane solvent &#43; ions. Only the water then?</div>
<div>Sorry that I didn't understand you explanation.</div>
<div>&nbsp;</div>
<div>Emma</div>
<div><br>
<div>&nbsp;</div>
</div>
<div style="FONT-SIZE: 16px; COLOR: rgb(0,0,0); FONT-FAMILY: 'Times New Roman'">
<div id="divRpF75352" style="DIRECTION: ltr">
<hr tabindex="-1">
<font face="Tahoma" color="#000000" size="2"><b>Från:</b><span class="Apple-converted-space">&nbsp;</span><a href="mailto:gmx-users-bounces@gromacs.org" target="_blank">gmx-users-bounces@gromacs.org</a><span class="Apple-converted-space">&nbsp;</span>[<a href="mailto:gmx-users-bounces@gromacs.org" target="_blank">gmx-users-bounces@gromacs.org</a>]
 för XAvier Periole [<a href="mailto:x.periole@rug.nl" target="_blank">x.periole@rug.nl</a>]<br>
<b>Skickat:</b><span class="Apple-converted-space">&nbsp;</span>den 27 april 2010 12:41<br>
<b>Till:</b><span class="Apple-converted-space">&nbsp;</span>Discussion list for GROMACS users<br>
<b>Ämne:</b><span class="Apple-converted-space">&nbsp;</span>Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box<br>
</font><br>
</div>
<div></div>
<div>
<div><br>
</div>
The crash seems to be expected as by crossing the pbc the distance will change significantly
<div>and in way the algorithm can not handle.&nbsp;</div>
<div><br>
</div>
<div>Note that the overall translational motion of your system should always be removed.&nbsp;</div>
<div><br>
</div>
<div>The removal of the COM motion of your bilayer should be sufficient to prevent the overall</div>
<div>motion of the bilayer. have a look at the following option in the mdp:</div>
<div>
<div>; mode for center of mass motion removal</div>
<div>comm-mode &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= Linear</div>
<div>; number of steps for center of mass motion removal</div>
<div>nstcomm &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div>
<div>; group(s) for center of mass motion removal</div>
<div>comm-grps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= membrane solvent&#43;ions</div>
<div><br>
</div>
</div>
<div><br>
<div>
<div>On Apr 27, 2010, at 12:07 PM, ERIKSSON, EMMA wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div>
<div style="FONT-SIZE: 13px; COLOR: rgb(0,0,0); DIRECTION: ltr; FONT-FAMILY: Tahoma">
<div>
<div>Hi all!</div>
<div>&nbsp;</div>
<div>I have a problem regarding lipid bilayer simulations in Gromacs 4. During some of my simulations the whole system is moving in the z direction in the box,&nbsp;meaning that&nbsp;after some time the lipids&nbsp;are going out in the bottom of the box and coming in in the
 top of the box, since I'm using periodic boundary conditions. This&nbsp;doesn't matter&nbsp;(I think) when running non-constrained simulations, however when I'm constraining the distance between the lipids (pull_geometry=cylinder) and a molecule in the system the system
 explodes and the simulations crash when the lipids are starting to cross over to the other side. The fact that the system is moving in the box must be the problem since the system explodes exactly when the first lipid&nbsp;passes over to the other side and nothing
 like this ever happens when the bilayer is not moving in the box. Is there any way to freeze the cylindrical COM of the lipids or something like that so that they stay more or less in the middle of the box all the time? I don't want to use freezegrps and freezedim
 = N N Y as this freezes the lipids completely in the z direction, and that's not what I want, I want them to be free to move as before but I want to stop the whole system from moving too much in the z direction. Anyone that has experienced a similar problem
 or know how to go about with this? I would really appreciate any help I can get. Thanks.</div>
<div>&nbsp;</div>
<div>Emma</div>
</div>
</div>
--<span class="Apple-converted-space">&nbsp;</span><br>
gmx-users mailing list &nbsp;&nbsp;&nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<span class="Apple-converted-space">&nbsp;</span><a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a><span class="Apple-converted-space">&nbsp;</span>before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<span class="Apple-converted-space">&nbsp;</span><br>
www interface or send it to<span class="Apple-converted-space">&nbsp;</span><a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read<span class="Apple-converted-space">&nbsp;</span><a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a></div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
--<span class="Apple-converted-space">&nbsp;</span><br>
gmx-users mailing list &nbsp;&nbsp;&nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<span class="Apple-converted-space">&nbsp;</span><a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a><span class="Apple-converted-space">&nbsp;</span>before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<span class="Apple-converted-space">&nbsp;</span><br>
www interface or send it to<span class="Apple-converted-space">&nbsp;</span><a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read<span class="Apple-converted-space">&nbsp;</span><a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a></div>
</span></blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</body>
</html>