Density functional theory BLYP/DNP was employed to optimize a series of fullerenes and their holmium endohedral compounds, including C(20), Ho@C(20), Ho(3+)@C(20), C(60), Ho@C(60), Ho(3+)@C(60),C(70), Ho@C(70), Ho(3+)@C(70) C(78), Ho@C(78), Ho(3+)@C(78), C(82),Ho@C(82) and Ho(3+)@C(82). DFT semi core pseudospot approximation was taken into consideration in the calculations of the element holmium because of its particular electronic structure. Fullerenes and their holmium endohedral compounds' aromaticity were studied in terms of structural criteria, energetic criteria, and reactivity criteria. The results indicate that the aromaticity of fullerenes was reduced when a holmium atom was introduced into the carbon cage, and the endohedral fullerenes' reactive activity enhance; but the aromaticity of the carbon cage increased when a Ho(3+) cation was encapsulated into a fullerene. Calculations of aromaticity and stability indicate that two paths can lead to the similar aim of preparing holmium endohedral fullerenes; that is, they can form from either a holmium atom or a holmium cation (Ho(3+)) reacting with fullerenes, respectively, and the latter is more favorable.