In this paper the relation between the molecular electron density distribution and the crystal growth morphology is investigated. Accurate charge densities derived from ab initio quantum chemical calculations were partitioned into multipole moments, to calculate the electrostatic contribution to the intermolecular interaction energy. For urea and alpha-glycine the F-faces or connected nets were determined according to the Hartman-Perdok PBC theory. From attachment energy and critical Ising temperature calculations, theoretical growth forms were constructed using different atom-atom potential models. These were compared to the Donnay-Harker model, equilibrium form and experimental growth forms. In the case of alpha-glycine, the theoretical growth forms are in good agreement with crystals grown from aqueous solution. Crystals obtained by sublimation seem to show some faces which are not F-faces sensu stricto.