Five different aromaticity indexes are benchmarked for benzene, pyridine and the diazines in their ground states. A basis set study was performed using the Pople style, Karlsruhe and Dunning's correlation consistent basis sets. Ten different DFT functionals, including LSDA, PBE, PBE0, B3LYP, CAM-B3LYP, wB97XD, M06-2X, SOGGA11X, M11 and MN15 were benchmarked by comparison with CCSD, CASSCF and MP2. Large out-of-plane imaginary frequencies were observed for some of the optimized structures at the correlated wavefunction level of theory. It was found that the DFT functionals in general predict the para-delocalization index, multicenter index and aromatic fluctuation index to be approximately 70%, 50% and 45% larger, respectively, compared to the CCSD method. Comparisons of the DFT functionals showed that the wB97XD, CAM-B3LYP and M06-2X functionals performed the best. Furthermore, the basis set dependence of the DFT functionals was found to be large for the electron sharing indexes. Based on these findings, it is recommended to perform ground state calculations of aromaticity indexes at the wB97XD, CAM-B3LYP or M06-2X level of theory utilizing a simple basis set of triple-ζ quality.