A set of rules for determining the atomic radii of spheres used to build the molecular cavities in continuum solvation models are presented. The procedure is applied to compute the hydration free energy for molecules containing H, C, N, O, F, P, S, Cl, Br, and I at a computational level ͑Hartree-Fock with a medium size basis set͒ allowing the study of relatively large systems. The optimized radii reduce the mean error with respect to the experimental solvation energies below 0.20 kcal/mol for a set of 43 neutral solutes and around 1 kcal/mol for 27 ions. Moreover the correct trends are observed for the solvation energies of homolog series, like the series ammoniatrimethylamine, that are not correctly reproduced by usual solvation models. © 1997 American Institute of Physics. ͓S0021-9606͑97͒50631-3͔