<br><br>
<div>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Thanks Pedro,<br><br>if I did not misunderstand you I should go for something like this:<br><br>pcoupl = berendsen
<br>pcoupltype = surface-tension<br>ref_p = 600 1.025<br>tau_p = 5.0 5.0<br>compressibility = 9.8e-5 4.53e-5<br><br>Is it correct?</blockquote>
<div>&nbsp;</div>
<div>Hmmm... &quot;Surface tension&quot; is a tricky subject&nbsp;at the microscopic level. </div>
<div>Besides, I&#39;m no expert, so I won&#39;t give you a &quot;correct&quot; answer, but only suggestions.&nbsp;</div>
<div>In the first place, I think your 600 bar.nm value is high for a bilayer. With a little calculation and approximations, it would be equivalent to -20 atm of &quot;lateral pressure&quot; (or, more precisely, lateral components of the pressure tensor -&gt; Pxx and Pyy).
</div>
<div>Also, I would use the 9.8e-5 for the zz compressibility. In fact, this is the most important value, as you can see in pages 26-27 of the manual.</div>
<div>&nbsp;</div>
<div>Regards.</div>
<div>&nbsp;</div>
<div>Pedro.</div><br>&nbsp;</div><br>