Computational Study of the Properties and Stabilities of Endohedral Metallofullerenes

Abstract

The chemistry of fullerenes, which are a class of carbon allotropes that can be prepared by vaporization of graphite in an electric arc in a low pressure atmosphere,1 has become a topic of much experimental and theoretical study over the past 25 years. Herein we present a series of theoretical studies related to recently discovered or studied endohedral metallofullerenes (EMF) and a theory as to the selective stability of certain isomers of EMFs. Computational treatments of the anions of C80 and C94 are presented and compared in an effort to gain an understanding and predictive model for which isomers of each cage size EMF will be most stable. A model is proposed in which the pentagons of fullerene anions are seen as charge localization centers that repel one another, making the pyracyclene bonding motif much more unstable for fullerene anions than for fullerene neutral cages. Computational treatments are also presented for two newly discovered EMFs, Y2C2@C92 and Gd2@C79N. Y2C2@C92 is reported to exhibit a previously undiscovered mode of internal cluster rotation, while Gd2@C79N is shown to have unusual stability for an azofullerene with a large spin quantum number (15/2). Finally, computational techniques are employed to predict the thermodynamic feasibility of a chemical reaction replacing one metal atom in a trimetallic-nitride template (TNT) endohedral metallofullerene with different metal atom. At least two of these are predicted to be thermodynamically practical.