The hydroxylated surface structures of seven morphologically important crystallographic surfaces of hematite were calculated allowing the determination of the surface energies and, consequently,the hematite equilibrium morphology. The docking of two organic phosphonate additives (methylnitrilo-dimethylenephosphonicacid—MNDP, ethylenediaminetetraphosphonicacid—EDTP) was then simulated on these surfaces and there placement energies calculated. With these energies, the effect of MNDP and EDTP on the hematite equilibrium morphology could be predicted. Results without the additives how a bipyramidal morphology of (102) faces capped with (222)faces at each end. The interaction with the EDTP additive predicts the stabilization of the (101)face to be far more pronounced than for the MND Padditive. This leads to the appearance of (101) faces in the morphology predicted in presence of EDTP. Experimental validation by means of electron microscopy shows morphologies close to those calculated, confirming that a computational approach can be used for the prediction of morphologies of crystals grown in the presence of additives.