An empirical continuum solvation model, solvation free energy density (SFED), has been developed to calculate solvation free energies of a molecule in the most frequently used solvents. A generalized version of the SFED model, generalized-SFED (G-SFED), is proposed here to calculate molecular solvation free energies in virtually any solvent. G-SFED provides an accurate and fast generalized framework without a complicated description of a solution. In the model, the solvation free energy of a solute is represented as a linear combination of empirical functions of the solute properties representing the effects of solute on various solute-solvent interactions, and the complementary solvent effects on these interactions were reflected in the linear expansion coefficients with a few solvent properties. G-SFED works well for a wide range of sizes and polarities of solute molecules in various solvents as shown by a set of 5,753 solvation free energies of diverse combinations of 103 solvents and 890 solutes. Octanol-water partition coefficients of small organic compounds and peptides were calculated with G-SFED with accuracy within 0.4 log unit for each group. The G-SFED computation time depends linearly on the number of nonhydrogen atoms (n) in a molecule, O(n).implicit solvent | macromolecules | linear time A n accurate description of the effect of solvent on solvation is crucial for understanding chemical and biological phenomena. Because many peptide and protein drugs are currently being developed (1, 2), computational efficiency as well as accuracy are very important. Among the different approaches that have been used, continuum solvation models are appealing because of the simplified yet accurate description of the solvent effect (3-8).Continuum solvent models have been developed based on the early work of Born (9), Bell (10), Kirkwood (11), and Onsanger (12). They focus primarily on the description of the solute either at the quantum mechanical level or as a classical collection of point charges and try to represent the influence of the solvent based on an approximate treatment as a dielectric continuum medium. The charge distribution of the solute induces electric polarization of the surrounding solvent, and the electric field generated by the polarized solvent, the reaction field, in turn perturbs the solute, leading to a modification of the charge distribution of the solute. This electrostatic problem of the mutual polarization between solute and solvent can be recast by using boundary conditions at the cavity surface, leading to significant simplifications in the associated equations.In the classical approaches, most work is concerned with describing the electrostatic contribution to the solvation in terms of the Poisson-Boltzmann equation (PBE) (13). Several different computational techniques for solving the PBE have been developed, for example, finite difference methods in which the dielectric property of the solvent is described in terms of a 3D grid around the solute (14-17) and boundary element methods i...