ABSTRACT:Fullerenes represent unique species in many respects, one of them being one of the highest temperature regions ever used for preparation of chemical species. This very factor automatically substantially enhances the importance of the TS part of the Gibbs function. The essential feature of fullerenic systems has consistently been ignored in computational studies of fullerenes, almost always drawing their conclusions from potential energies or heats of formation. This report demonstrates that the simplified treatment cannot really describe all the features of complex mixtures of isomeric fullerenes. General formulas are surveyed and illustrated on C 96 -the biggest set of isomeric fullerenes characterized by 13 C NMR spectra to date. Computations are carried out for the complete set of 187 isolated-pentagon-rule (IPR) isomers of C 96 , using four semiempirical quantum chemical methods (MNDO, AM1, PM3, and SAM1). Entropy terms are also computed and the relative-stability problem is entirely treated in terms of the Gibbs function. In overall, the computed data match the observed data reasonably well, although no method can really reproduce the complete ten-membered observed set. There are however still several aspects to be checked and several possibilities how to improve quality of the computed terms. In principle, such developments are, however, only a question of available computing power. There are, nevertheless, some remaining questions, such as inclusion of anharmonicity effects, contribution from electronic partition function, or even interactions between individual, yet separated motions. This article also discusses a subsequent, more general task of relative stabilities of carbon cages of different dimensions, i.e., relative stabilities of nonisomeric fullerenes. The computational treatment of isomeric mixtures has several interesting features: the results depend on temperature, but not on pressure; only the relative, and not the absolute, values of the heats of formation are needed; and the form of the master equation allows for an ample cancellation of terms in the partition functions. In contrast, for the more general task-relative stabilities of fullerenes of different dimensions-the absolute values of the heats of formation are required, the mentioned cancellations do not operate there, and the scheme is essentially pressure dependent.