The workflow -CREST to B97-3c to DFT -is then applied to exploring the ensembles of a series of polycyclic aromatic and heteroaromatic homodimers. Studying these particular systems provides insight into the effects which are produced by substituting heteroatoms into polycyclic aromatic systems as they relate to non-covalent interactions, particularly π -π stacking.