We report herein on the use of all-atom molecular dynamics simulations to investigate the solvation environment of C-60 and four C-60-derived fullerenes immersed in a variety of aromatic solvents. Utilizing a recently developed solvation shell analysis technique that quantifies the spatial relationships between fullerenes and solvent on a molecular level, we show that the number of fullerene substituents and solvent chemistry are crucial determinants of the solvation shell structure and thus fullerene solvation behavior. Specifically, it is shown for the derivatives investigated that the number of fullerene substituents is more critical to solvation behavior than the substituent chemistry. (C) 2017 Elsevier B.V. All rights reserved.