<pre>Dear Gromacs users,</pre>I have do simulation of surfactant at oil/water interface for a long time. But I find it too difficult to get a reasonable interfacial tension. The key question is the pcoupltype.<br><pre>My system: </pre>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Box size:  3.0 * 3.0 * 13 <br>Molecule Number: decane 230, water 2041<br>
pcoupltype = semiisotropic </blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Gromacs version: 4.5.5</blockquote>
<blockquote style="margin:0 0 0 40px;border:none;padding:0px"></blockquote> <div>1. When I set compressibility =  4.5e-5 4.5e-5.  Surfactant number is 32 . The box will be shrink at x and y. At last simualion will end for a fatal error:</div>
<div><br></div><div><div>One of the box vectors has become shorter than twice the cut-off length or box_yy-|box_zy| or box_zz has become smaller than the cut-off.</div><div>For more information and tips for troubleshooting, please check the GROMACS.</div>
<div><br></div><div>2. When I add the number of surfactants to 40. compressibility = 4.5e-5 4.5e-5. The interfacial tension is -20bar/nm for two interfaces.  It is not reasonable because it should be positive. </div><div>
<br></div><div>3. When I compressibility = 0 4.5e-5.  Surfactant number is 32. The interfacial tension is 794bar/nm for two interfaces.</div><div><br></div><div><br></div><div>My question is:</div><div><br></div><div>What is the right way to simulate for this system?  My work confused me too much.  Very appreciate for you help. </div>
<div><br></div><div>Hugh</div><div><br></div><div><pre><br></pre><pre>  </pre></div></div>