It seems almost all the common implementations of Generalized Born (GB) implicit solvent models (in Amber, OpenMM etc) use the "pairwise descreening" approximation, which is described in Hawkins, Cramer & Truhlar 1995[1]: ... $H_k$ is eclipsed surface area of sphere of radius $r$ centered at ato...