Four methods of computational quantum chemistry are used in a study of hyperconjugation in protonated aromatic molecules. Benzene, benzenium, toluene, and four isomeric forms of toluenium are examined using the self-consistent field level of theory followed by configuration interaction and coupled cluster calculations, as well as density functional theory. Results for proton affinities, geometric parameters, atomic populations, dipole moments, and polarizabilities are reported. The calculated results are in good agreement with previous computational studies and with experimental data. The presence of hyperconjugation is evident from the shortened carbon-carbon bond lengths in the aromatic ring and concomitant changes in dipole moments and polarizabilities. The proton affinities of benzene and toluene compare well with experimental values. The examination of all of the toluenium isomers reveals that the position of the methyl group has a minor impact on the strength of hyperconjugation, although the most stable isomer is found to be the para form. Mulliken population analyses indicate that the addition of a proton contributes to aromatic hyperconjugation and increases the strength of π-bonds at the expense of σ-bonds.