[gmx-users] ffamber03 top file problem

gromacs564 gromacs564 at 126.com
Tue Feb 15 09:57:29 CET 2011


Dear All:
 I want to obtain a top file for protein with ffamber03 force field,and correct some bonds and angles parameter.
 I added a new atom type HD(D2O) in ffamber03  [ atomtypes ], so I added some HD parameters in ffbond.itp、atomtypes.atp and ffnonbond.itp files,the output of the program indicate:


-------------------------------------------------------
Program grompp_m, VERSION 4.5.1
Source code file: toppush.c, line: 1631

Fatal error:
Incorrect number of parameters - found 1, expected 2 or 4 for Bond.
For more information and tips for troubleshooting, please check the GROM···
--------------------------------------------------------------------------------------------



I added some bonds and angles parameters like this:
······························
; residue  44 ASN rtp ASN  q  0.0
   664          N     44    ASN      N    664  -0.430106      14.01   ; qtot 4.57
   665         HD    44    ASN     HD    665   0.250466      2.016   ; qtot 4.824
   666         CT     44    ASN     CA    666   0.044609      12.01   ; qtot 4.869
···············
[ bonds ]
;  ai    aj funct            c0            c1            c2            c3
···
  676   678     1 
  678   679     1     0.10115 
  678   680     1 
  680   681     1 
·······················
[ angles ]
;  ai    aj    ak funct            c0            c1            c2            c3
························
  648   650   652     1 
  651   650   652     1     115.149 
  650   652   653     1 
  650   652   654     1 
····················

Any help will highly appreciated!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110215/34dac669/attachment.html>


More information about the gromacs.org_gmx-users mailing list