<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
On 17/10/2011 11:41 PM, Mikhail Stukan wrote:
<blockquote
cite="mid:757CC336D2D08E4BB8D36B3AD5960B874722DADF@NL0230MBX06N1.DIR.slb.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<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]-->
<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><br>
</p>
</div>
</blockquote>
<br>
I think the only way to do something like this in GROMACS is with a
COM virtual site and such a distance restraint from each water
molecule to that virtual site.<br>
<br>
Mark<br>
</body>
</html>