Efficient quantum-chemical geometry optimization and the structure of   large icosahedral fullerenes
  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. 
Brett I. Dunlap and Rajendra. R. Zope (US Naval Research Laboratory   and University of Texas at El Paso)