[gmx-users] GB Parameters

Elizabeth Ploetz ploetz at ksu.edu
Mon Aug 23 19:27:03 CEST 2010

Dear Gromacs Users, 

I am interested in implementing implicit solvent using the OBC GBSA method.  In the Gromacs 4.5beta3 top directory, several force fields (oplsaa and multiple versions of amber) contain the gbsa.itp file, and one (charmm27), contains a gb.itp file.  They all have the same column titles (sar, st, pi, gbr, and hct).  I have read the six articles listed at http://www.gromacs.org/Documentation/Terminology/Implicit_Solvent, however I still have the following two questions regarding the contents of this file:

(1.)  What values are tabulated in these columns and what are the corresponding references for each of them?
(2.)  Which columns are used (versus ignored) by Gromacs in the implementation of each of the methods (Still/HCT/OBC)?

Here are my guesses for Question 1:
(a.)  sar: solvent accessible surface area?
(b.)  st: surface tension?
(c.)  pi: ?
(d.)  gbr: generalized born radius?  These values correspond to sigma/2 from Jorgensen and Tirado-Rives' 1988 JACS paper (v. 110, p. 1657).  However the generalized born radius is equal to Rvdw plus a dielectric offset distance, according to Still et al's 1997 J. Phys. Chem. A paper (v. 101, p. 3005).  This would mean the dielectric offset distance is set to zero (?).  
(e.)  hct: HCT parameter

With an eventual goal of implementing implicit solvent using my lab's force field, insight into the meaning of the columns and/or references for their origin would be greatly appreciated. 

Thanks in advance! 

Sincerely yours,

E. A. Ploetz

More information about the gromacs.org_gmx-users mailing list