[gmx-users] Lennard-Jones Parameters in ffnonbonded.itp

Lara Bunte lara.bunte at yahoo.de
Mon Feb 20 18:36:07 CET 2012



Hello
I use a charmm force field in gromacs. I almost finished my parametrization except for Lennard-Jones parameters. For this parameters I want to adopt the values in the charmm force field for the basic atom types. 

In gromacs I have to put this parameters in the ffnonbonded.itp file in the section [ pairtypes ] and I have to give following informations: 

; i     j       func    sigma1-4        epsilon1-4 ; THESE ARE 1-4 INTERACTIONS

in example

CP1     CP1     1       0.338541512893  0.04184

Now here is a section of information given in the charmm documentation. 

!V(Lennard-Jones) = Eps,i,j[(Rmin,i,j/ri,j)**12 - 2(Rmin,i,j/ri,j)**6]
!
!epsilon: kcal/mole, Eps,i,j = sqrt(eps,i * eps,j)
!Rmin/2: A, Rmin,i,j = Rmin/2,i + Rmin/2,j
!
!atom  ignored    epsilon      Rmin/2   ignored   eps,1-4       Rmin/2,1-4
!
HN1      0.0       -0.0460    0.2245 
HN2      0.0       -0.0460    0.2245 


The question is: 


How to adopt the values in the charmm force field for the basic atom types in gromacs? I do not really know what I have to write in the [ pairtypes ] section in gromacs. 


Thank you for helping me
Greetings
Lara




More information about the gromacs.org_gmx-users mailing list