[AMBER] EPTOT shows big difference between the different force fields

From: zhouhaibin2008.ok <zhouhaibin2008.ok.163.com>
Date: Fri, 5 Aug 2011 11:58:14 +0800 (CST)

Hi all,
Like amber tutorial B3, I wanted to simulate TC5b in linear state to see if it can reach the folded state . In this simulation , I used two different force fields , ff03.r1 and ff99SB, to built the linear states . When run sander , I used the same input file for the two simulations to optimise , heat the system and run MD production . Here is the MD production input file :
 &cntrl
 imin = 0, irest = 1,
 ntx = 7, ntb = 0,
 ntr = 0, igb = 5,
 cut = 999.9, rgbmax = 999.9, ntc = 2,
 ntf = 2, ntt = 3,
 gamma_ln = 4.0, tempi = 325.0,
  temp0 = 325.0, nstlim = 5000000,
 dt = 0.002, ntwr = 500,
 ntpr = 500, ntwx = 500
 /
When these simulations ended , I found the linear state using the ff03.r1 force field reached the folded state , but the linear state using ff99SB not . However when I check the out file , I found the EPTOT using ff03.r1 was about -200, but the EPTOT using ff99SB was about -400.I don't known why the difference is so big , and I think the folded state should have a lower EPTOT, but the it does not.
Does the difference result from the different charges used by different force fields .

Thank you in advance !



_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Aug 04 2011 - 21:00:02 PDT
Custom Search