In this study, the interaction energies between aromatic organic molecules (AOMs) and cellulose are explored using density functional theory (DFT) through the software SIESTA and the exchange-correlational functional VDW-DRSLL. Three AOMs will be modelled: benzene, benzamide and benzoic acid. Firstly, the interaction energies of the dimers of the AOMs are determined. Then, the obtained interaction energies of the cellulose-AOM complexes are compared to the former in order to decide which interaction is stronger. It is found that the studied AOMs are more likely to interact with cellulose than with another identical monomer; benzamide has the highest propensity to interact, followed by benzoic acid and benzene. Furthermore, for all interaction energy calculations a counterpoise correction term will be introduced as an addition to the SIESTA optimisation and it will be shown that without this correction the acquired energy minima will deviate significantly from accepted values from previous studies.