We report the development of a new, fast, surface‐based method for numerical calculations of the solvation energy of biomolecules with a large number of charged groups. The procedure scales linearly with the system size, both in time and memory requirements, produces explicit values for the reaction field potential and stable values of the polar energy within only a few percent error margin practically for any molecular configurations. The method works well both for large and small molecules and thus gives stable energy differences for quantities such as the polar energy contributions to molecular complex formation energies. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011