Using Tersoff bond order potential, a vibrational analysis of the spherical fullerene family, namely C-60, C-80, C-180, C-240, C-260, C-320, C-500, and C-720 was performed. To evaluate the validity of our results, we have compared our simulation results with available experimental data and also with DFT B3LYP/6-31G( d) calculations. In general, molecular stiffness tends to decrease with increasing size, but its variation is limited in cases where mostly the tension-compression interaction sites are active such as the breathing mode. Furthermore, the bond length of each molecule is derived and compared with experimental and theoretical values calculated for graphene. Finally, vibrational frequencies are plotted in a histogram to reveal the common frequency gap and concentration points of the frequency distribution. (C) 2016 Elsevier B. V. All rights reserved.