We investigated molecular polarity by using theoretical means and comparing empirical solvent polarity. Our approach employed electrostatic potentials at the molecular surface calculated by density functional methods. A number of molecular descriptors related to molecular polarities were computed from molecular surface electrostatic potentials. Among computed molecular descriptors, the most positive electrostatic potential provided the best correlation with the empirical solvent polarities. A regression equation was developed in order to predict molecular polarities of molecules whose experimental solvent polarities were unknown. The new regression equations were utilized in estimating solvent polarities of cubane derivatives which are considered important precusors of high-energy density meterials.