<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:AdvOT10f0fcb1;
        panose-1:0 0 0 0 0 0 0 0 0 0;}
@font-face
        {font-family:"AdvOT92137dae\.I";
        panose-1:0 0 0 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";
        font-weight:bold;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-weight:bold;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-weight:normal">Dear gmx-user,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-weight:normal"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-weight:normal">I would like to simulate an isolated drop of water (without PBC). To prevent water molecules form evaporation I would like to create a spherical constraint around the system. In recent article in PNAS
(Caleman, Hub, van Maaren, van der Spoel, v. 108, 6838 (2011)) where the authors considered similar system it is mentioned that:<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"AdvOT10f0fcb1","sans-serif";font-weight:normal">“To avoid such evaporation, we applied a spherical flat-bottom quadratic potential acting on the water. That potential
was implemented as an additional force </span><span style="font-size:10.0pt;font-family:"AdvOT92137dae.I","sans-serif";font-weight:normal">F
</span><span style="font-size:10.0pt;font-family:"AdvOT10f0fcb1","sans-serif";font-weight:normal">pointing toward the COM of the droplet…”<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-weight:normal"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-weight:normal">Could anybody give a hint how this approach can be realized in GROMACS? I looked through the manual, but did not find clear answer how to implement such a force.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-weight:normal"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-weight:normal">Thank you very much in advance.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-weight:normal">Mikhail<o:p></o:p></span></p>
<p class="MsoNormal"><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:"Arial","sans-serif";font-weight:normal">=============================================<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:"Arial","sans-serif";font-weight:normal">Dr Mikhail Stukan<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Arial","sans-serif";font-weight:normal">Schlumberger Dhahran Carbonate Research Center,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Arial","sans-serif";font-weight:normal">Dhahran Techno Valley - KFUPM,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Arial","sans-serif";font-weight:normal">P.O. Box 39011, Dammam / Doha Camp 31942,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Arial","sans-serif";font-weight:normal">Kingdom of Saudi Arabia</span><span style="font-family:"Arial","sans-serif";font-weight:normal"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:"Arial","sans-serif";font-weight:normal">Tel: +966 3 331 6182<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-family:"Arial","sans-serif";font-weight:normal">Fax:+966 3 330 0845<br>
<a href="mailto:mstukan@slb.com"><span style="color:blue">mstukan@slb.com</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>