  Geometry optimization is efficient using generalized Gaunt coefficients,
which significantly limit the amount of cross differentiation for multi-center
integrals of high-angular-momentum solid-harmonic basis sets. We parameterize
the fully analytic formulation of density functional theory (ADFT), called the
Slater-Roothaan method, developed in our group to give the exact geometry of
C60 fullerene. The parametrized ADFT is subsequently used to optimize
geometries of most stable C240, C540, C960, C1500 and C2160 icosahedral
fullerenes. The calculations are all electron, the orbital basis set includes d
functions and the exchange-correlation-potential basis set includes f
functions. The calculation of C2160 fullerene employed about 39000 basis
functions and is the largest calculation reported on any isolated molecule
to-date. The evolution of interatomic distance and atomization energy from C60
to graphite has been investigated.
