For the standard vdW cavity construction, PEqS calculations reproduce analytic results to < 0.01 kcal/mol accuracy for the molecules FH and ClH, although the solvation energies themselves are |∆G elst | = 0.1-0.2 kcal/mol. For the Voronoi construction, the error for FH is large in relative terms but is only 0.1 kcal/mol in absolute terms, while the error for ClH is negligibly small. The NO - 3 ion has a much larger solvation energy and larger absolute errors for PEqS as compared to the analytic Kirkwood result, on the order of 0.3-0.4 kcal/mol and thus similar to discretization errors for the monatomic ions. As in that case, this error is only about 1% of |∆G elst |. This level of agreement suggests that we can use the PEqS code, in conjunction with the Voronoi construction for ε(r), to generate near-exact benchmarks that can be used as reference values for other methods. As such, only the Voronoi construction of ε(r) will be considered in subsequent PEqS calculations.