This work focused on the interactions that occur between ester solvents and simple aromatic solutes. Solutes were selected from various functional groups in their simplest form, and predictions of activity coefficients at infinite dilution were made using the Modified UNIFAC Dortmund group contribution model. The model computation was set up on a Microsoft Excel spreadsheet specifically designed for this purpose. For polar aromatic solutes, solubility decreased with increase in size of the ester solvent molecule and the opposite was found to be true for non-polar solutes. For all aromatic/ ester solvent interactions there was a decrease in activity coefficients with an increase in the degree of ester unsaturation