<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>Gromacs has the half as a pre-factor in the potential.<br><br>Note that in 4.0.5 there is a g_wham program which can automatically<br>do the WHAM analysis on the mdrun output.<br><br>Berk<br><br>&gt; Date: Tue, 1 Sep 2009 13:19:32 +0800<br>&gt; From: mircial@sjtu.edu.cn<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: [gmx-users] Force constant for umbrella sampling potential<br>&gt; <br>&gt; Hello everyone,<br>&gt; <br>&gt; I would like to know does the harmonic biasing potential for running<br>&gt; umbrella sampling in gromacs have a 0.5 in front or is the 0.5 term<br>&gt; incorporated into the force constant value that we specify?<br>&gt; On other words is the potential V=0.5*k(x-xo)^2 or V=k(x-xo)^2.<br>&gt; <br>&gt; I am using Grossfieldã] WHAM<br>&gt; algorithm and it assumes that the biasing potential used has a 0.5 in front,<br>&gt; but I know for example CHARMM does not include the 0.5 and I just need to<br>&gt; know if this is also the case with gromacs
 <br>&gt; <br>&gt;   Thank you in advance.<br>&gt; <br>&gt; R-X  Gu<br>&gt; <br>&gt; <br>&gt; _______________________________________________<br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br><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>