<br><br>----- Original Message -----<br>From: Ramachandran G &lt;gtrama@gmail.com&gt;<br>Date: Thursday, July 1, 2010 11:42<br>Subject: [gmx-users] amber/charmm force field and HB lifetime at low        temperature<br>To: gmx-users@gromacs.org<br><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Hi gromacs users:<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;&nbsp;&nbsp;&nbsp; Using amber/charmm force field we simulated a solvated protein <br><br>I've never heard of such a forcefield. What do you mean?<br><br>&gt; system at different temperature(300,250 and 200K).<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>I used TIP5P water model and applied NVT ensemble . We analyzed the hydrogen bond life time correlation function(HBCF) to study the protein water interactions. At 300K, the result looks good and comparable with the literature data. But at 250 and 200K, during the initial 0-2ps, HBCF decay is very fast before it start to relax.<br><br>Forcefields are typically parameterized on data around laboratory conditions, so trying to make observations on systems whose conditions are&nbsp; significantly far from those is doubtful. I don't have any real idea what "decay is very fast before it start to relax" means, either.<br><br>&gt; This problem does not arise when i tried with other force field which came along with gromacs package(like GROMOS,OPLS..) I don't know where i am making mistakes; patching the amber and charmm force fields to gromacs package lead to this problem? <br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;"></font><br>Mark