<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Hi,<br><br>There are no settings that matter in particular for pulling compared to a normal simulation.<br>All the settings that matter normally matter here as well.<br>NVT pulling should give nearly identical results no NPT,<br>as long as you set the volume correctly.<br><br>Berk<br><br><hr id="stopSpelling">To: gmx-users@gromacs.org<br>Subject: Re : [gmx-users] PULL CODE AND NpT ensemble<br>Date: Thu, 27 Aug 2009 08:03:56 -0400<br>From: gigitang5@aol.fr<br><br>
<div> <font face="Arial, Helvetica, sans-serif">Hi Berk,<br>
Thanks for your answer. I just want to be sure that the procedure is correct. No I do not look at the pressure/volume fluctutation. But I do want to make a clever choice of the ensemble so that no significant extra artifact are generated, since the number of contraint to be calculated is not negligible. I previously used no P coupling and I was wondering is that feature could explain that the average force plot is not smooth function. Other criteria such as Cut-off value for LJ and coulomb may also be checked out?<br>
Geraldine <br>
</font></div>
<div> <br>
</div>
<div> <br>
</div>
-----E-mail d'origine-----<br>
De : Berk Hess <gmx3@hotmail.com><br>
A : Discussion list for GROMACS users <gmx-users@gromacs.org><br>
Envoyé le : Jeudi, 27 Août 2009 13:45<br>
Sujet : RE: [gmx-users] PULL CODE AND NpT ensemble<br>
<br>
<div id="EC_AOLMsgPart_3_25200f8e-ab4d-4c8b-a490-2ab0f7c5763d">
<style>
.ExternalClass #EC_AOLMsgPart_3_25200f8e-ab4d-4c8b-a490-2ab0f7c5763d .EC_hmmessage P
{padding:0px;}
.ExternalClass #EC_AOLMsgPart_3_25200f8e-ab4d-4c8b-a490-2ab0f7c5763d body.EC_hmmessage
{font-size:10pt;font-family:Verdana;}
</style>
Hi,<br>
<br>
I would use v-rescale and Berendsen pressure coupling.<br>
Although Berendsen pressure coupling, strictly speaking, does not give the correct ensemble,<br>
in practice it produces negligible artifacts, unless you are looking at pressure/volume fluctuations.<br>
<br>
Berk<br>
<br>
<hr id="EC_stopSpelling">To: <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
Date: Thu, 27 Aug 2009 07:20:00 -040
0<br>
From: <a href="mailto:gigitang5@aol.fr">gigitang5@aol.fr</a><br>
Subject: [gmx-users] PULL CODE AND NpT ensemble<br>
<br>
<font face="Arial, Helvetica, sans-serif">Hi all,<br>
I would like to use pull code to evaluate the free energy of pulling system1 from system2. Free energy calculation should be done for NpT ensemble. All tutorials found for calculating free energy used <br>
1) berensen for T and P coupling or <br>
2) berensen for T coupling AND Parrinelo for P coupling or<br>
3)Vrescale for T and no P coupling<br>
<br>
I am a bit lost in which one should lead to real ensemble.<br>
Is any one could help?<br>
Regards,<br>
geraldine, phD <br>
</font><br>
<hr>Express yourself instantly with MSN Messenger! <a href="http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/">MSN Messenger</a>=
</div>
<div id="EC_AOLMsgPart_4_25200f8e-ab4d-4c8b-a490-2ab0f7c5763d" style="font-family: Tahoma,Verdana,Arial,Sans-Serif; font-size: 12px; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);">
<pre style="font-size: 9pt;"><tt>_______________________________________________<br><br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br><br>Please search the archive at <a href="http://www.gromacs.org/search">http://www.gromacs.org/search</a> before posting!<br><br>Please don't post (un)subscribe requests to the list. Use the <br><br>www in<br>terface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br><br>Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a><br><br></tt></pre>
</div><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>