The classical and quantum structures of cationic polycyclic aromatic hydrocarbon molecules (benzene, pyrene, coronene, and circumcoronene) coated by helium atoms have been theoretically investigated using a variety of computational methods. Classical shell filling, as determined from global optimization, is generally found to proceed by epitaxial additions on the graphitic surfaces before peripheral closure. From the quantum mechanical perspective provided by path-integral molecular dynamics simulations, vibrational delocalization is found to generally decrease the size of this first solvation shell, but also to give rise to a variety of situations in which the helium atoms are more or less localized depending on their environment, with strong finite size effects depending both on the hydrocarbon cation and the number of coating helium atoms. While the graphitic planes tend to bind helium sufficiently to give rise to snowball precursors, the peripheral regions are less dense and more delocalized, not as liquid as the outer layers but within an intermediate slushy character.