[gmx-users] converting the normal hormonic potential bond to Fourth power potential bond?

[라지브간디]

Dear gmx,


I am getting confusion with converting the normal hormonic potential bond ( 
which is used in charmm27) to Fourth power potential bond ( gromos43a1). 


For Example. 
The value of C-O force constant, bond stretch and distance are in charmm is 
518816.0, 0.123 and 0.12. Hence I used this value for converting to gromos type 
by using Fourth power potential bond equation and end up with value of 
1.29189542054402. But the value of force constant in gromos for C-O is 
mentioned like this 


#define gb_4        0.1230  1.6600e+07
; C  - O        1200


Here i dont understand what i am missing here. 


Some guidance or suggestion or really appreciated. Thanks.

-- 
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists