Quantum Chemical Insight into Metallofullerenes M2C98: M2C2@C96

Oct 28, 2013 - Quantum Chemical Insight into Metallofullerenes M2C98: M2C2@C96 or M2@C98, Which Will Survive? Hong Zheng†, Xiang Zhao*†, Wei-Wei ...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/JPCC

Quantum Chemical Insight into Metallofullerenes M2C98: M2C2@C96 or M2@C98, Which Will Survive? Hong Zheng,† Xiang Zhao,*,† Wei-Wei Wang,† Jing-Shuang Dang,† and Shigeru Nagase‡ †

Institute for Chemical Physics & Department of Chemistry, State Key Laboratory of Electrical Insulation and Power Equipment, Xi’an Jiaotong University, Xi’an 710049, China ‡ Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan S Supporting Information *

ABSTRACT: Systematic studies on M2C98 (M = Sc, Y, La, Gd, Lu) series by means of density functional theory (DFT) methods disclose that having large concentrations within the temperature interval of fullerene formation, the metal−carbide endohedral fullerenes M2C2@C96 are found to be more stable than the dimetallofullerenes M2@C98. When encapsulating M2C2 clusters, three types of C96 cages will present excellent thermodynamic stabilities and become the largest metal−carbide endohedral fullerenes so far. The dynamic mechanisms of the internal C2 moieties in the most important M2C2@C96 species have been discussed, which provide profound understanding of metal−carbide endohedral fullerenes. Moreover, 13C NMR and IR spectra have been simulated to aid further experimental identification and characterization.



atoms or metalloclusters.15−19 On the other hand, in the case of M2C2n fullerenes, there always are disputations on M2C2@C2n−2 vs. M2@C2n, which have been persisted for a long time.20−28 Many middle-size M2C2n EMFs tend to take M2C2@C2n−2 model,20,23−29 but there are also reports on the stable M2@C2n EMFs which have been discovered in experimental or theoretical works such as Sc2@C70,22 Sc2@C80,30 La2@C72,31 Td2@C76,32 Dy2@C100,18 etc. Consequently, cautious studies on M2C2n EMFs including the non-IPR isomers indicate that the competition between M2C2@C2n−2 and M2@C2n should be taken into consideration. As for Gd2C98 EMF which was discovered in 200833 by Balch et al., one recent theoretical study has pointed out that a nonIPR isomer (Gd2@C1(168785)-C98) with one pentagon adjacency (PA) should be the most proper product in experiment.19 Nevertheless, it should be noted that the authors did not discuss the relative stabilities of Gd2C2@C96 isomers as comparison. The real structures of such high-spin EMFs, which determine their properties and further applications, should be discussed based on more comprehensive theoretical computations. Besides, Sc2C98 and Y2C98 have been synthesized even earlier in 1994 without isolation.34 One common character among these three kinds of rare earth metal elements (Sc, Y, Gd) is that they have similar valence electron structures as nd1(n + 1)s2. Moreover, as another two types of popular endohedral metal atoms in EMFs, La and Lu also have the same valence electron structures as Gd. Which kind of M2C98

INTRODUCTION Ever since the first endohedral fullerene La@C82 was discovered in 1991,1 endohedral metallofullerenes (EMFs) have attracted wide interests during the past few decades for their unique structural, electronic, and magnetic properties and are expected to be applied in many fields such as electronics, photovoltaics, and biomedicine.2−8 The charge transfer from encapsulated metal atoms or metalloclusters to fullerene cages in the EMFs will stabilize many fullerene isomers which cannot exist in neutral hollow state. Particularly, some fullerene cages violating the well-known isolated-pentagon rule (IPR),9,10 designated non-IPR fullerenes, can become more stable than IPR fullerenes when accept some metal atoms or metalloclusters.8,11−15 Therefore, theoretical study on the stability of EMFs should take both IPR and non-IPR isomers into consideration. As for the large size EMFs, many dimetallofullerenes and trimetallic nitride endohedral fullerenes with carbon cages ≥C90 have been synthesized without structural qualitative analysis. Since many large-scale EMFs can hardly be isolated and characterized in experiments, deep insights into their geometry structures and other significant properties can only be supplied by theoretical methods. However, detailed theoretical studies about the large-scale EMFs is still lacking, and previous works always neglected some important aspects when discussing their relative stabilities. On one hand, recent theoretical reports have better stabilities than the IPR ones when encapsulating certain metal atoms or metalloclusters.15−19 Recent theoretical reports have revealed that the stabilities of some non-IPR structures may surpass most IPR structures, although not always can become the most stable ones, when encapsulating certain metal © 2013 American Chemical Society

Received: October 6, 2013 Revised: October 28, 2013 Published: October 28, 2013 25195

dx.doi.org/10.1021/jp409915t | J. Phys. Chem. C 2013, 117, 25195−25204

The Journal of Physical Chemistry C

Article

Figure 1. Relative concentrations of low-energy M2C98 (M = Sc, Y, La, Gd, Lu) isomers.

C96 cages has been probed to supply much deeper understanding of M2C2@C96 EMFs. In addition, the simulated NMR and IR spectra of most stable M2C98 isomers have been calculated to aid further experimental characterization.

structure should be the most stable isomer, and what are the similarities and differences among these five M2C98 (M = Sc, Y, La, Gd, Lu) series? Only more systematic studies could figure out these questions and disclose some important properties of M2C98 EMFs which would be necessary and interesting. In this paper, comprehensive analysis on M2C98 (M = Sc, Y, La, Gd, Lu) series have been carried out with DFT methods based on full screening of 17941 C98 isomers and 12840 C96 isomers (PA ≤ 2). In conjunction with statistical mechanics calculations, the studies on relative stabilities reveal that M2C2@C96 isomers have optimal stabilities in all M2C98 series. Size effects of cages and metal atoms have been discussed in detail. The nanoscale parameters of various M2C2 clusters and the interaction between metalloclusters and C96 cages which have a crucial influence on the macroscopic properties of materials have been presented. Further analysis on electronic structures discloses electron transfer between the metal atoms and C98 cages. The dynamic behavior of the M2C2 clusters in



EXPERIMENTAL SECTION Since each of the encapsulated metal atoms would donate three electrons in the M2@C98 isomers and the M2C2 clusters always keep a +4 valence state in the fullerene cages, the energetics for all those 17 941 C98 isomers and 12 840 C96 isomers with no more than two PA fragments were first screened at the AM135 level on hexa- and tetraanion states, respectively. Then certain most stable isomers (including IPR and non-IPR isomers) were reoptimized at the hybrid density functional theory B3LYP36−38 level with 6-31G* basis set. The results of C986− series are from the previous report (ref 19). Geometry optimizations of M2@ C98 and M2C2@C96 series with five kinds of rare earth atoms were performed at the B3LYP/3-21G-CEP level (3-21G basis 25196

dx.doi.org/10.1021/jp409915t | J. Phys. Chem. C 2013, 117, 25195−25204

The Journal of Physical Chemistry C

Article

development of relative concentrations of M2C98 series within a broad temperature region is presented in Figure 1. The relative concentrations of Gd2C98 series present outstanding stability of Gd2C2@C2(191809)-C96 in a wide temperature interval (Figure 1d). Obviously, Gd2C2@C2(191809)-C96 always shares the largest relative concentration until the temperature reaches around 4200 K, meaning it will possess the best thermodynamic stability among all the isomers in the temperature interval for metallofullerene formation (below 3000 K).46 b. Relative Stabilities of M2C98 (M = Sc, Y, La, Lu) Series. The energy sequences of the rest of M2C98 series are presented in Table S2b−e, also taking the comparison of two probable structural patterns into consideration. Specially, the relative energies of M2@C98 (M = Sc, Y, La, Lu) series are introduced in Table S3 determining the lowest-energy structures of each M2@C98 series. Sc2@C2(230924)-C98, Y2@C1(168785)-C98, La2@C1(168785)-C98, and Lu2@C2(230979)-C98 are the lowest-energy isomers of the four M2@C98 series. Being similar to the Gd2C98 series, Sc2C2@C2(191809)-C96 is predicted as the lowest-energy structure of Sc2C98 series, and Sc2@C2(230924)-C98, the most stable Sc2@C98 isomer as shown in Table S3, has a large relative energy of 19.70 kcal/mol (Table S2b). A non-IPR isomer Sc2C2@Cs(189891)-C96 with one PA fragment has the second lowest energy. Reoptimizations at the B3LYP/6-31G*-CEP level show good accordance to the B3LYP/3-21G-CEP results in both energy sequences and HOMO−LUMO gaps. C2(191809)-C96 is also the carrier cage of the lowest-energy isomer of the Y2C98 series when accepting a Y2C2 cluster (Table S2c), followed by Y2C2@ D2(191838)-C96 with a large relative energy of 10.65 kcal/mol. However, in La2C98 series, La2C2@C2(191809)-C96 loses the energy advantage since it is 4.76 kcal/mol higher than La2@ C1(168785)-C98. The energy sequence of Lu2C98 series is much similar to that of Y2C98 series, which supports Lu2C2@ C2(191809)-C96 as the lowest-energy isomer. Statistical thermodynamic analyses of Sc2C98 series (Figure 1a) show contribution interchanges at around 3000−3200 K, where isomer Sc2C2@C2(191809)-C96 is surpassed by Sc2C2@ D 2 (191838)-C 96 and Sc 2 C 2 @C 1 (191812)-C 96 . Sc 2 C 2 @ C2(191809)-C96 shares the largest contributions in the temperature interval for metallofullerene formation, and so it should be the most important thermodynamic stability isomer. That Sc2C2@D2(191838)-C96 and Sc2C2@C1(191812)-C96 show notable fractions at high temperature indicates that both isomers could also be synthesized experimentally with smaller yields. The lowest-energy Sc2@C98 isomer Sc2@ C2(230924)-C98 is extremely disfavored by the entropy effect. For Y2C98 series, Y2C2@C2(191809)-C96 has the largest contributions below 2600 K (Figure 1b), and then, very surprisingly, it is surpassed by another isomer Y2C2@ D2d(191815)-C96 which has a large relative energy of 14.43 kcal/mol. Thus, the composition of Y2C98 should be mainly composed of Y2C2@C2(191809)-C96 and Y2C2@D2d(191815)C96. Although La2@C1(168785)-C98 is predicted as the lowestenergy isomer of La2C98 series, it is obviously disfavored by the entropy effect since its contribution quickly descends after 200 K and is surpassed by La2C2@C2(191809)-C96 at 580 K. The La2C2@C2(191809)-C96 structure is completely dominant in the temperature interval for metallofullerene formation, while La2C2@D2(191838)-C96 shows notable fractions over 30% at high temperatures as subordinate isomer. Like the case of Y2C98 fullerenes, in the Lu2C98 series, the contribution of the lowestenergy isomer Lu2C2@C2(191809)-C96 is surpassed by Lu2C2@

set for C atoms and CEP-31G basis set with the corresponding pseudopotential for metal atoms), just as the reported Gd2@ C98 series in ref 19. Specially, for the first three lowest-energy isomers of Sc2C98 series, reoptimizations with a larger basis set 6-31G* for C atoms has been carried out, and the results are introduced in Table S2b. Rotational−vibrational partition functions were supplied by the computed structural and vibrational data at the B3LYP/3-21G-CEP level of theory. The simulations of NMR spectra for M2C98 isomers using the gauge-independent atomic orbital (GIAO) method39 are performed at the B3LYP/6-311G-CEP level. The GIAO method has been proven to be very sufficiently accurate at simulating 13C NMR chemical shifts of both fullerenes and metallofullerenes.24,40−42 All computational works above were carried out with the Gaussian03 program package.43



RESULTS AND DISCUSSION As depicted in Table S1, the D2(191838)-C964− cage is predicted as the most stable C964− anions, followed by another IPR isomer D2d(191815)-C964− with a small relative energy and a much larger HOMO−LUMO gap (1.60 eV). The smallest relative energy of non-IPR cages lies 16.6 kcal/mol above D2(191838)-C964−, indicating poor stabilities of non-IPR C964− anions. Relative Stabilities of M2C98 (M = Sc, Y, La, Gd, Lu) Series. a. Relative Stabilities of Gd2C98. As the most arrestive M2C98 EMFs which has been discussed before, the Gd2C98 series including both Gd2C2@C96 and Gd2@C98 structural patterns have shown rather surprising differences comparing to previous report according to our new results. The potential energy of Gd2C2@C2(191809)-C96 is 5.72 kcal/mol lower than that of the Gd2@C2(230924)-C98 which was predicted as the lowest-energy structure19 (Table S2a) and becomes the lowest energy isomer with a large HOMO−LUMO gap of 1.72 eV. The second lowest-energy isomer Gd2C2@D2(191838)-C96 is 9.26 kcal/mol higher than Gd2C2@C2(191809)-C96. Such a sequence reversal of EMF energy in comparison with the charged series of pristine fullerenes could be found in many theoretical studies on EMF systems.15−17,44 Usually discussion on the charged cages just deals with the possible charge transfer to empty cages. However, in an EMF system the charge transfer (from metallocluster to carbon cage) is not the only factor influencing its thermodynamic stability, steric strain of the whole EMF molecule (due to limited space of fullerene cage), and bonding between metal atom and carbon cage could also profoundly affect the thermodynamic stability of the EMF system.15,22 Therefore, it is reasonable that sometimes the energy sequence of EMF isomers may present some difference to that of the charged cages. Since stability interchanges induced by the enthalpy−entropy interplay always occur in EMFs systems, the relative stabilities of M2C98 series at high temperatures cannot be predicted by only the separation energy itself. Thus, to gain deep insight into the thermodynamic stability of the M2C98 fullerenes, we have investigated the entropy effect and evaluated the relative concentrations based on equilibrium statistical thermodynamic analysis.45 Such thermodynamic investigations have successfully depicted relative stability panoramas for various EMF series.14−16,18,19 Here, the isomers with lower energy gaps (1800 cm−1) in most of the M2C2@C96 isomers, which corresponds to C−C stretching of the internal C2 moieties.

concentrations. Having evident relative energies, the M2C2@ D2d(191815)-C96 (M = Y, Lu) structures are favored by the entropy effect with even larger concentrations than M2C2@ C2(191809)-C96 (M = Y, Lu) structures. As a consequence, the M2C2@C96 species present outstanding thermodynamical stabilities, and they are the largest metal−carbide endohedral fullerenes which are theoretically discovered heretofore. Observations on the geometry structures of M2C2@C96 isomers recover their evident dependence on the size of metal atoms. The analysis of electronic structures for M2C2@C96 isomers discloses that their electronic properties are mainly related to that of the carrier cages, but independent of the geometric structures of M2C2 clusters and the types of metal element. Interestingly, oscillation and rotation mechanisms of the internal C2 moieties in the most important M2C2@C96 species (excepting La2C2@D2(191838)-C96) with rather small barriers have been discovered, suggesting that these EMFs should maintain the same symmetries as their carrier cages. The predicted 13C NMR spectra of various M2C2@C96 species expected to aid future experiment reflect that the chemical shift of the internal C2 moiety may have a strong dependence on the compression and C−C distance of M2C2 clusters in the C96 cages. In addition, the IR spectra of the most significant M2C2@C96 isomers are theoretically predicted to recover the vibration characteristics which would be helpful for future experimental identification. To sum up, our work not only provides a profound understanding of thermodynamic stabilities and electronic/geometric peculiarities of M2C98 metallofullerenes but also assists further experimental research on these large size endohedral fullerenes containing rare earth metal atoms.



ASSOCIATED CONTENT

* Supporting Information S

Relative energies and HOMO−LUMO gaps of partial C964− isomers and M2C98 series, energetics comparison between two kinds of M2C24+ cations, simulated IR spectra of M2C2@C98 isomers, complete author lists of refs 34, 43, and 47, and coordinates of main M2C2@C98 isomers. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*Fax +86 29 8266 8559; Tel +86 29 8266 5671; e-mail xzhao@ mail.xjtu.edu.cn (X.Z.).



Notes

CONCLUSIONS Detailed studies on M2C98 (M = Sc, Y, La, Gd, Lu) series by means of DFT calculations in conjunction with statistical thermodynamic analysis reveal that unlike the reported Gd2@ C98 series, the most stable isomers prefer to be metal−carbide endohedral fullerenes. The M2C2@C2(191809)-C96 structures become the lowest-energy isomers in most of the M2C98 series, except in the La2C98 series a La2@C98 structure is lower than isomer La2C2@C2(191809)-C96. Further studies by means of statistical thermodynamic analysis disclose that the most thermodynamically important isomers in all series should choose M2C2@C96 form rather than M2@C98 form. The M2C2@C2(191809)-C96 (M = La, Gd) structures should be the main product in electric arc since they share the largest concentrations within the temperature interval of fullerene formation. The M2C2@D2(191838)-C 96 (M = Sc, La) structures can also be synthesized due to their comparable

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work has been supported by the National Natural Science Foundation of China (21171138), National Key Basic Research Program of China (2012CB720904), and China Postdoctoral Science Foundation (2013M532038).



REFERENCES

(1) Chai, Y.; Cuo, T.; Jin, C.; Haufler, R. E.; Chibante, L. P. F.; Fure, J.; Wang, L.; Alford, J. M.; Smalley, R. E. Fullerenes with Metals Inside. J. Phys. Chem. 1991, 95, 7564−7568. (2) Wilson, L. J.; Cagle, D. W.; Thrash, T. P.; Kennel, S. J.; Mirzadeh, S.; Alford, J. M.; Ehrhardt, G. J. Metallofullerene Drug Design. Coord. Chem. Rev. 1999, 190−192, 199−207. (3) Akasaka, T.; Nagase, S. Endofullerenes: A New Family of Carbon Clusters; Kluwer Academic Publishers: Dordrecht, 2002.

25202

dx.doi.org/10.1021/jp409915t | J. Phys. Chem. C 2013, 117, 25195−25204

The Journal of Physical Chemistry C

Article

(25) Yumura, T.; Sato, Y.; Suenaga, K.; Iijima, S. Which Do Endohedral Ti2C80 Metallofullerenes Prefer Energetically: Ti2@C80 or Ti2C2@C78? A Theoretical Study. J. Phys. Chem. B 2005, 109, 20251− 20255. (26) Iiduka, Y.; Wakahara, T.; Nakajima, K.; Tsuchiya, T.; Nakahodo, T.; Maeda, Y.; Akasaka, T.; Mizorogi, N.; Nagase, S. 13C NMR Spectroscopic Study of Scandium Dimetallofullerene, Sc2@C84 vs. Sc2C2@C82. Chem. Commun. 2006, 2057−2059. (27) Lu, X.; Nakajima, K.; Iiduka, Y.; Nikawa, H.; Mizorogi, N.; Slanina, Z.; Tsuchiya, T.; Nagase, S.; Akasaka, T. Structural Elucidation and Regioselective Functionalization of An Unexplored Carbide Cluster Metallofullerene Sc2C2@Cs(6)-C82. J. Am. Chem. Soc. 2011, 133, 19553−19558. (28) Lu, X.; Nakajima, K.; Iiduka, Y.; Nikawa, H.; Tsuchiya, T.; Mizorogi, N.; Slanina, Z.; Nagase, S.; Akasaka, T. The Long-Believed Sc2@C2v(17)-C84 Is Actually Sc2C2@C2v(9)-C82: Unambiguous Structure Assignment and Chemical Functionalization. Angew. Chem., Int. Ed. 2012, 51, 5889−5892. (29) Feng, Y.; wang, T.; Wu, J.; Feng, L.; Xiang, J.; Ma, Y.; Zhang, Z.; Jiang, L.; Chu, C.; Wang, C. Structural and Electronic Studies of Metal Carbide Clusterfullerene Sc2C2@Cs-C72. Nanoscale 2013, 5, 6704− 6707. (30) Wu, J.; Wang, T.; Shu, C.; Lu, X.; Wang, C. Identification of the Most Stable Sc2C80 Isomers: Structure, Electronic Property, and Molecular Spectra Investigations. Chin. J. Chem. 2012, 30, 765−770. (31) Kato, H.; Taninaka, A.; Sugai, T.; Shinohara, H. Structure of a Missing-Caged Metallofullerene: La2@C72. J. Am. Chem. Soc. 2003, 125, 7782−7783. (32) Umemoto, H.; Ohashi, K.; Inoue, T.; Fukui, N.; Sugai, T.; Shinohara, H. Synthesis and UHV-STM Observation of the Tdsymmetric Lu Metallofullerene: Lu2@C76(Td). Chem. Commun. 2010, 46, 5653−5655. (33) Yang, H.; Lu, C.; Liu, Z.; Jin, H.; Che, Y.; Olmstead, M. M.; Balch, A. L. Detection of a Family of Gadolinium-Containing Endohedral Fullerenes and the Isolation and Crystallographic Characterization of One Member as a Metal-Carbide Encapsulated inside a Large Fullerene Cage. J. Am. Chem. Soc. 2008, 130, 17296− 17300. (34) Stevenson, S.; Dorn, H. C.; Burbank, P.; Harich, K.; Haynes, J.; Klang, C. H., Jr.; Salem, J. R.; DeVrles, M. S.; van Loosdrecht, P. H. M.; Johnson, R. D.; et al. Automated HPLC Separation of Endohedral Metallofullerene Sc@C2n and Y@C2n Fractions. Anal. Chem. 1994, 66, 2675−2679. (35) Dewar, M. J. S.; Zoebisch, E.; Healy, E. F.; Stewart, J. J. P. Development and Use of Quantum Mechanical Molecular Models. 76. AM1: A New General Purpose Quantum Mechanical Molecular Model. J. Am. Chem. Soc. 1985, 107, 3902−3909. (36) Becke, A. D. Density Functional Exchange Energy Approximation with Correct. Phys. Rev. A 1988, 38, 3098−3100. (37) Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (38) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785−789. (39) Wolinski, K.; Hinton, J.; Pulay, P. Efficient Implementation of the Gauge-Independent Atomic Orbital Method for NMR Chemical Shift Calculations. J. Am. Chem. Soc. 1990, 112, 8251−8260. (40) Sun, G.; Kertesz, M. Theoretical 13C NMR Spectra of IPR Isomers of Fullerenes C60, C70, C72, C74, C76, and C78 Studied by Density Functional Theory. J. Phys. Chem. A 2000, 104, 7398−7403. (41) Kobayashi, K.; Nagase, S. A Stable Unconventional Structure of Sc2@C66 Found by Density Functional Calculations. Chem. Phys. Lett. 2003, 362, 373−379. (42) Yang, T.; Zhao, X.; Osawa, E. Can a Metal−Metal Bond Hop in the Fullerene Cage? Chem.Eur. J. 2011, 17, 10230−10234. (43) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A.; Vreven, Jr., T.; Kudin, K. N.; Burant, J. C.; et al. Gaussian03 (Revision E.01); Gaussian, Inc.: Wallingford, CT, 2004.

(4) Dunsch, L.; Yang, S. Metal Nitride Cluster Fullerenes: Their Current State and Future Prospects. Small 2007, 3, 1298−1320. (5) Akasaka, T.; Wudl, F.; Nagase, S. Chemistry of Nanocarbons; Wiley-Blackwell: London, 2010. (6) Lu, X.; Akasaka, T.; Nagase, S. Chemistry of Endohedral Metallofullerenes: The Role of Metals. Chem. Commun. 2011, 47, 5942−5957. (7) Lu, X.; Akasaka, T.; Nagase, S. Carbide Cluster Metallofullerenes: Structure, Properties, and Possible Origin. Acc. Chem. Res. 2013, 46, 1627−1635. (8) Popov, A. A.; Yang, S.; Dunsch, L. Endohedral Fullerenes. Chem. Rev. 2013, 113, 5989−6113. (9) Fowler, P. W.; Manolopoulos, D. E. An Atlas of Fullerenes; Clarendon Press: Oxford, UK, 1995. (10) Kroto, H. W. The Stability of the Fullerenes Cn, with n = 24, 28, 32, 36, 50, 60 and 70. Nature 1987, 329, 529−531. (11) Kobayashi, K.; Nagase, S.; Yoshida, M.; Osawa, E. Endohedral Metallofullerenes. Are the Isolated Pentagon Rule and Fullerene Structures Always Satisfied? J. Am. Chem. Soc. 1997, 119, 12693− 12694. (12) Wang, C. R.; Kai, T.; Tomiyama, T.; Yoshida, T.; Kobayashi, Y.; Nishibori, E.; Takata, M.; Sakata, M.; Shinohara, H. Materials Science: C66 Fullerene Encaging a Scandium Dimer. Nature 2000, 408, 426− 427. (13) Stevenson, S.; Fowler, P. W.; Heine, T.; Duchamp, J. C.; Rice, G.; Glass, T.; Harich, K.; Hajdu, E.; Bible, R.; Dorn, H. C. Materials Science: A Stable Non-classical Metallofullerene Family. Nature 2000, 408, 427−428. (14) Yang, T.; Zhao, X.; Xu, Q.; Zhou, C.; He, L.; Nagase, S. NonIPR Endohedral Fullerene Yb@C76: Density Functional Theory Characterization. J. Mater. Chem. 2011, 21, 12206−12209. (15) Zheng, H.; Zhao, X.; Ren, T.; Wang, W. W. C74 Endohedral Metallofullerenes Violating the Isolated Pentagon Rule: A Density Functional Theory Study. Nanoscale 2012, 4, 4530−4536. (16) Zheng, J.; Zhao, X.; Dang, J.; Chen, Y.; Xu, Q.; Wang, W. Density Functional Theory Characterization of Lanthanum Nitride Endohedral Fullerene: La3N@C92. Chem. Phys. Lett. 2011, 514, 104− 108. (17) Popov, A. A.; Dunsch, L. Structure, Stability, and Cluster-Cage Interactions in Nitride Clusterfullerenes M3N@C2n (M = Sc, Y; 2n = 68−98): A Density Functional Theory Study. J. Am. Chem. Soc. 2007, 129, 11835−11849. (18) Yang, T.; Zhao, X.; Nagase, S. Di-lanthanide Encapsulated into Large Fullerene C100: A DFT Survey. Phys. Chem. Chem. Phys. 2011, 13, 5034−5037. (19) Zhao, X.; Gao, W. Y.; Yang, T.; Zheng, J. J.; Li, L. S.; He, L.; Cao, R. J.; Nagase, S. Violating the Isolated Pentagon Rule (IPR): Endohedral Non-IPR C98 Cages of Gd2@C98. Inorg. Chem. 2012, 51, 2039−2045. (20) Kurihara, H.; Lu, X.; Iiduka, Y.; Mizorogi, N.; Slanina, Z.; Tsuchiya, T.; Akasaka, T.; Nagase, S. Sc2C2@C80 Rather Than Sc2@ C82: Templated Formation of Unexpected C2v(5)-C80 and Temperature-Dependent Dynamic Motion of Internal Sc2C2 Cluster. J. Am. Chem. Soc. 2010, 133, 2382−2385. (21) Sc2C2@C68 vs. Sc2@C70: Shi, Z. Q.; Wu, X.; Wang, C. R.; Lu, X.; Shinohara, H. Isolation and Characterization of Sc2C2@C68: A MetalCarbide Endofullerene with a Non-IPR Carbon Cage. Angew. Chem., Int. Ed. 2006, 45, 2107−2111. (22) Zheng, H.; Zhao, X.; Wang, W. W.; Yang, T.; Nagase, S. Sc2@ C70 Rather Than Sc2C2@C68: Density Functional Theory Characterization of Metallofullerene Sc2C70. J. Chem. Phys. 2012, 137, 014308. (23) Cao, B.; Hasegawa, M.; Okata, K.; Tomiyama, T.; Okazaki, T.; Suenaga, K.; Shinohara, H. EELS and 13C NMR Characterization of Pure Ti2@C80 Metallofullerene. J. Am. Chem. Soc. 2001, 123, 9679− 9680. (24) Tan, K.; Lu, X. Ti2C80 Is More Likely a Titanium Carbide Endohedral Metallofullerene (Ti2C2)@C78. Chem. Commun. 2005, 4444−4446. 25203

dx.doi.org/10.1021/jp409915t | J. Phys. Chem. C 2013, 117, 25195−25204

The Journal of Physical Chemistry C

Article

(44) Yang, T.; Zhao, X.; Nagase, S. Structural Determination on Yb@ C78 Reveals an Unexpected Relationship of Yb@C2n (2n = 74−80). J. Phys. Chem. C 2012, 116, 21640−21645. (45) Slanina, Z. Equilibrium Isomeric Mixtures: Potential Energy Hypersurfaces as the Origin of the Overall Thermodynamics and Kinetics. Int. Rev. Phys. Chem. 1987, 6, 251−267. (46) Slanina, Z.; Kobayashi, K.; Nagase, S. Ca@C74 Isomers: Relative Concentrations at Higher Temperatures. Chem. Phys. 2004, 301, 153− 157. (47) Zhang, J.; Fuhrer, T.; Fu, W.; Ge, J.; Bearden, D. W.; Dallas, J.; Duchamp, J.; Walker, K.; Champion, H.; Azurmendi, H.; et al. Nanoscale Fullerene Compression of an Yttrium Carbide Cluster. J. Am. Chem. Soc. 2012, 134, 8487−8493. (48) Nishimoto, Y.; Wang, Z.; Morokuma, K.; Irle, S. Molecular and Electronic Structures of Endohedral Fullerenes, Sc2C2@C3v-C82 and Sc2@C3v-C82: Benchmark for SCC-DFTB and Proposal of New Inner Cluster Structures. Phys. Status Solidi B 2012, 249, 324−334. (49) Inoue, T.; Tomiyama, T.; Sugai, T.; Shinohara, H. Spectroscopic and Structural Study of Y2C2 Carbide Encapsulating Endohedral Metallofullerene: (Y2C2)@C82. Chem. Phys. Lett. 2003, 382, 226−231. (50) Inoue, T.; Tomoyama, T.; Sugai, T.; Okazaki, T.; Suematsu, T.; Fujii, N.; Utsumi, H.; Nojima, K.; Shinohara, H. Trapping a C2 Radical in Endohedral Metallofullerenes: Synthesis and Structures of (Y2C2)@ C82 (Isomers I, II, and III). J. Phys. Chem. B 2004, 108, 7573−7579.

25204

dx.doi.org/10.1021/jp409915t | J. Phys. Chem. C 2013, 117, 25195−25204