In this work, the Bingel–Hirsch addition of diethylbromomalonate to all non-equivalent bonds of Sc3N@D3h-C78 was studied using density functional theory calculations. The regioselectivities observed computationally allowed the proposal of a set of rules, the predictive aromaticity criteria (PAC), to identify the most reactive bonds of a given endohedral metallofullerene based on a simple evaluation of