A quantitative, predictive mechanism of metallofullerene cage formation: kinetics, route selection, and transient intermediates of the transformation network anchored by the C1(51383)-C84 'missing link' cage

Whether fullerenes and endohedral metallofullerenes (EMFs) form 'bottom-up' from small carbon fragments or 'top-down' by shrinkage of graphenic or giant cages has been open since the discovery of $C_{60}$ in 1985. Zhang et al. (Nature Chemistry 5, 880–885; 2013) isolated the asymmetric, pentalene-containing cage $C_1(51383)$-$C_{84}$ in $M_2C_2@C_{84}$ (M = Y, Gd) and proposed it as a preserved 'missing link' from which most well-known, high-symmetry EMF cages are reachable through enumerated Stone–Wales rotations and $C_2$ extrusions, explicitly framing the result as a plausible, experimentally anchored pathway network rather than a quantitative theory, and declaring fullerene formation 'still an open question' requiring further experimental and theoretical effort.

Background

Fullerenes are closed carbon cages built from pentagons and hexagons, of which the 60-atom buckminsterfullerene ($I_h$-$C_{60}$) is the archetype. Endohedral metallofullerenes (EMFs) are fullerenes that encapsulate metal atoms or clusters — for example a metal carbide $M_2C_2$ unit — and are typically produced by arc-discharge vaporization of metal-doped graphite at temperatures of several thousand kelvin. Almost all isolable fullerenes obey the isolated pentagon rule (IPR): no two pentagons share an edge. Cages that violate it, such as those containing a pentalene motif (two pentagons fused along a bond) or a seven-membered (heptagon) ring — so-called nonclassical cages — are destabilized and only rarely isolated, usually because the encapsulated cluster stabilizes the strained region. How these cages assemble has been debated since 1985. 'Bottom-up' scenarios posit accretion of small carbon fragments (classically $C_2$ units) into chains, rings, bowls and finally closed cages; 'top-down' scenarios posit that graphenic sheets or giant cages curl and then shrink by successive fragment loss to the observed small fullerenes. Two elementary cage operations are central to both: the Stone–Wales transformation (SWT), rotation of a cage bond that reshuffles pentagon–hexagon adjacency at constant cage size, and extrusion or insertion of a $C_2$ unit, which changes the cage size by two carbons, in some motifs passing through heptagon-containing intermediates.

In 2013, Zhang et al. reported the single-crystal X-ray structure of $Gd_2C_2@C_1(51383)$-$C_{84}$ and the ${}^{13}$C NMR characterization of $Y_2C_2@C_1(51383)$-$C_{84}$: a chiral, asymmetric $C_{84}$ cage containing pentalene motifs, i.e., precisely the kind of destabilized, non-IPR geometry predicted to appear transiently in a top-down shrinkage cascade. They proposed this cage as a preserved 'missing link' and constructed a transformation map in which specific $C_2$ extrusions from the pentalene and from indene motifs (a pentagon fused to a hexagon), together with Stone–Wales steps, convert it into many of the best-known, high-symmetry EMF cages ($C_s(6)$-$C_{82}$, $C_{2v}(9)$-$C_{82}$, $C_{3v}(8)$-$C_{82}$, $I_h(7)$-$C_{80}$ and others), a network that accounts for the majority of solvent-extractable metallofullerenes across cluster types (metallic, nitride, carbide, sulfide, oxide, cyanide and methanide). They further argued that large endohedral groups require large starting cages and that, during shrinkage, different sizes and charges of the encapsulated units cause the cage to select different routes in the map. Crucially, the authors presented the map as structural evidence anchoring a plausible pathway network — no step had a measured or computed barrier, and no intermediate had been observed under formation conditions — and declared fullerene formation 'still an open question', calling for 'further experimental and theoretical efforts with consideration of specific reaction conditions ... to unambiguously unravel the formation mechanism of fullerenes and metallofullerenes' and for further study of other, less stable missing-link fullerenes and metallofullerenes.

The subsequent literature has substantially reshaped the qualitative picture while leaving that quantitative agenda untouched. Transformation maps of the same topological kind were extended to europium mono-EMFs ($Eu@C_{2n}$, 2n = 74–86, claimed to cover about 98% of metallofullerenes with established cage structures in that size range), to uranium mono-metallofullerenes (where both top-down shrinkage and bottom-up growth were concluded to proceed simultaneously in the arc reactor), and to giant cages beyond $C_{100}$. Heptagon-containing cages were trapped as kinetically stabilized products. Charge transfer from the encapsulated cluster and geometric cluster–cage matching were identified as factors directing which cages form. And several studies — precursor-size-resolved mass spectrometry, quantum molecular dynamics of carbon vapor, in situ laser-vaporization studies, and the 2023 isolation of $Sc_2C_2@C_1(51383)$-$C_{84}$ itself as a metastable direct-$C_2$-insertion product of smaller, more stable precursors — established that bottom-up growth, top-down shrinkage and post-formation annealing interconversions are superimposed. What no later work provides is the quantitative core: a parameterized, predictive account of which routes are taken, with what probability, under which conditions.

Problem Statement

Establish a quantitative, predictive mechanism of fullerene and metallofullerene cage formation, anchored to the experimentally preserved asymmetric-cage transformation network of Zhang et al. (2013): convert the topological map based on the $C_1(51383)$-$C_{84}$ 'missing link' cage — a topological scaffold whose steps are not yet energetically or kinetically characterized — into a physical theory of which formation routes are actually taken, with what probability, under which conditions. Such a quantitative account comprises the following facets; facets (1), (3) and (5) together determine the partition between top-down shrinkage, bottom-up growth and post-formation interconversion for the standard arc-discharge synthesis of solvent-extractable metallofullerenes:

(1) Measured or computed energetic and kinetic parameters for the elementary steps of the network — activation barriers and rate constants for the specific Stone–Wales rotations and for $C_2$ extrusions from pentalene versus indene motifs, and lifetimes of the proposed heptagon-containing intermediates — sufficient to predict route probabilities through the map under arc-synthesis conditions.

(2) The generality and prevalence of the cascade across synthesis conditions and carbon sources: whether multiple distinct 'missing-link' asymmetric cages form under other conditions, and their abundance distribution among the products.

(3) A predictive account of how the size, charge state and electronic structure of the encapsulated atom or cluster bias which transformation route is taken (the source notes that 'depending on the size and charge of the endohedral units, the fullerene cage will choose different routes').

(4) Direct observation or trapping of the proposed transient intermediates of the map — for example the heptagon-containing precursors — under formation-relevant conditions, distinguishing genuine cascade intermediates from kinetically stabilized end products or post-synthesis artifacts.

(5) A quantitative, predictive model that connects the statistics of the initial carbon population (carbon vapor, graphenic fragments or giant cages) to the measured fullerene and metallofullerene product distributions.

The scientific goal is a quantitative account capable of validation against measured product distributions: route probabilities that are predicted rather than enumerated, with the partition between top-down shrinkage, bottom-up growth and post-formation interconversion determined for the standard arc-discharge synthesis of solvent-extractable metallofullerenes.

References

1. Zhang, J.; Bowles, F. L.; Bearden, D. W.; Ray, W. K.; Fuhrer, T.; Ye, Y.; Dixon, C.; Harich, K.; Helm, R. F.; Olmstead, M. M.; Balch, A. L.; Dorn, H. C. A missing link in the transformation from asymmetric to symmetric metallofullerene cages implies a top-down fullerene formation mechanism. Nature Chemistry 5, 880–885 (2013). DOI: 10.1038/nchem.1748. https://doi.org/10.1038/nchem.1748
2. Dunk, P. W.; Mulet-Gas, M.; Nakanishi, Y.; Kaiser, N. K.; Rodríguez-Fortea, A.; Shinohara, H.; Poblet, J. M.; Marshall, A. G.; Kroto, H. W. Bottom-up formation of endohedral mono-metallofullerenes is directed by charge transfer. Nature Communications 5, 5844 (2014). DOI: 10.1038/ncomms6844. https://doi.org/10.1038/ncomms6844
3. Zhang, Y.; Ghiassi, K.; Deng, Q.; Samoylova, N.; Olmstead, M. M.; Balch, A. L.; Popov, A. A. Synthesis and Structure of LaSc2N@Cs(hept)-C80 with one heptagon and thirteen pentagons. Angewandte Chemie International Edition 54, 495–499 (2015). DOI: 10.1002/anie.201409094. https://doi.org/10.1002/anie.201409094
4. Cai, W.; Li, F.; Bao, L.; Xie, Y.; Lu, X. Isolation and Crystallographic Characterization of La2C2@Cs(574)-C102 and La2C2@C2(816)-C104: Evidence for the Top-Down Formation Mechanism of Fullerenes. Journal of the American Chemical Society 138, 6670–6675 (2016). DOI: 10.1021/jacs.6b03934. https://doi.org/10.1021/jacs.6b03934
5. Hao, Y.; Tang, Q.; Li, X.; Zhang, M.; Wan, Y.; Feng, L.; Chen, N.; Slanina, Z.; Adamowicz, L.; Uhlik, F. Isomeric Sc2O@C78 Related by a Single-Step Stone–Wales Transformation: Key Links in an Unprecedented Fullerene Formation Pathway. Inorganic Chemistry 55, 11354–11361 (2016). DOI: 10.1021/acs.inorgchem.6b01894. https://doi.org/10.1021/acs.inorgchem.6b01894
6. Chen, C.; Abella, L.; Cerón, M. R.; Guerrero-Ayala, M.; Rodríguez-Fortea, A.; Olmstead, M. M.; Powers, X.; Balch, A. L.; Poblet, J. M.; Echegoyen, L. Zigzag Sc2C2 Carbide Cluster inside a [88]Fullerene Cage with One Heptagon, Sc2C2@Cs(hept)-C88: A Kinetically Trapped Fullerene Formed by C2 Insertion? Journal of the American Chemical Society 138, 13030–13037 (2016). DOI: 10.1021/jacs.6b07912. https://doi.org/10.1021/jacs.6b07912
7. Han, J.; Choi, T.; Kim, S.; Lee, J.; Ha, Y.; Jeong, K.; Kim, H.; Choi, H.; Kim, H. Probing Distinct Fullerene Formation Processes from Carbon Precursors of Different Sizes and Structures. Analytical Chemistry 88, 8232–8238 (2016). DOI: 10.1021/acs.analchem.6b02076. https://doi.org/10.1021/acs.analchem.6b02076
8. Qian, H.; Wang, Y.; Morokuma, K. Quantum mechanical simulation reveals the role of cold helium atoms and the coexistence of bottom-up and top-down formation mechanisms of buckminsterfullerene from carbon vapor. Carbon 114, 635–641 (2017). DOI: 10.1016/j.carbon.2016.12.062. https://doi.org/10.1016/j.carbon.2016.12.062
9. Mulet-Gas, M.; Abella, L.; Cerón, M. R.; Castro, E.; Marshall, A. G.; Rodríguez-Fortea, A.; Echegoyen, L.; Poblet, J. M.; Dunk, P. W. Transformation of doped graphite into cluster-encapsulated fullerene cages. Nature Communications 8, 1222 (2017). DOI: 10.1038/s41467-017-01295-9. https://doi.org/10.1038/s41467-017-01295-9
10. Berné, O.; Montillaud, J.; Joblin, C. Top-down formation of fullerenes in the interstellar medium. Astronomy & Astrophysics 577, A133 (2015). DOI: 10.1051/0004-6361/201425338. https://doi.org/10.1051/0004-6361/201425338
11. Bao, L.; Yu, P.; Pan, C.; Shen, W.; Lu, X. Crystallographic identification of Eu@C2n (2n = 88, 86 and 84): completing a transformation map for existing metallofullerenes. Chemical Science 10, 2153–2158 (2019). DOI: 10.1039/c8sc04906h. https://doi.org/10.1039/c8sc04906h
12. Li, M.; Zhao, Y.; Han, Y.; Yuan, K.; Nagase, S.; Ehara, M.; Zhao, X. Theoretical Investigation of the Key Roles in Fullerene-Formation Mechanisms: Enantiomer and Enthalpy. ACS Applied Nano Materials 3, 547–554 (2020). DOI: 10.1021/acsanm.9b02110. https://doi.org/10.1021/acsanm.9b02110
13. Wu, B.; Jiang, L.; Luo, Y.; Wang, C. The Effect of the Polyaromatic Hydrocarbon in the Formation of Fullerenes. Angewandte Chemie International Edition 59, 3942–3947 (2020). DOI: 10.1002/anie.201915228. https://doi.org/10.1002/anie.201915228
14. Cai, W.; Alvarado, J.; Metta-Magaña, A.; Chen, N.; Echegoyen, L. Interconversions between Uranium Mono-metallofullerenes: Mechanistic Implications and Role of Asymmetric Cages. Journal of the American Chemical Society 142, 13112–13119 (2020). DOI: 10.1021/jacs.0c04888. https://doi.org/10.1021/jacs.0c04888
15. Yin, F.; Liu, X.; Tan, K. Key heptagon-containing metallofullerene intermediates in the LaSc2N@D(6)-C80 formation pathway: A theoretical study. Computational and Theoretical Chemistry 1138, 158–162 (2018). DOI: 10.1016/j.comptc.2018.06.014. https://doi.org/10.1016/j.comptc.2018.06.014
16. Li, M.; Zhao, R.; Dang, J.; Zhao, X. Theoretical study on the stabilities, electronic structures, reaction and formation mechanisms of fullerenes and endohedral metallofullerenes. Coordination Chemistry Reviews 471, 214762 (2022). DOI: 10.1016/j.ccr.2022.214762. https://doi.org/10.1016/j.ccr.2022.214762
17. He, J.; Li, M.; Zhang, W.; Zhao, X. Stabilities, Geometries, Electronic Structures, and Conversion Rules of Carbide Cluster Metallofullerenes. The Chemical Record 22, e202200148 (2022). DOI: 10.1002/tcr.202200148. https://doi.org/10.1002/tcr.202200148
18. Yao, Y.; Chen, Z.; Chen, L.; Zheng, S.; Yang, S.; Deng, S.; Echegoyen, L.; Tan, Y.; Xie, S.; Zheng, L. Two Metastable Endohedral Metallofullerenes Sc2C2@C1(39656)-C82 and Sc2C2@C1(51383)-C84: Direct-C2-Insertion Products from Their Most Stable Precursors. Journal of the American Chemical Society 145, 16778–16786 (2023). DOI: 10.1021/jacs.3c04840. https://doi.org/10.1021/jacs.3c04840
19. Hou, Y.; Kong, X. Endometallofullerenes in the Gas Phase: Progress and Prospect. Inorganics 12, 68 (2024). DOI: 10.3390/inorganics12030068. https://doi.org/10.3390/inorganics12030068
20. Jia, S.; Hu, Z.; Shao, X.; Huang, C.; Zhang, J.; Guo, K.; Pan, C.; Chen, W.; Bao, L.; Lu, X. Directional Breathing of Fullerene Cages: Crystallographic Evidence for Endohedral C2-Triggered Elongation of Cs(6)-C82 Cage. Inorganic Chemistry 64, 16572–16577 (2025). DOI: 10.1021/acs.inorgchem.5c02622. https://doi.org/10.1021/acs.inorgchem.5c02622
21. Wu, D.; Liu, X.; Li, Y.; Ying, H.; Zhang, H.; Huo, L.; Cui, C.; Shi, Z.; Lin, S.-C.; Xie, S.-Y.; Zheng, L.-S. Revealing the chlorine-mediated mechanism of fullerene pyrolysis synthesis by in situ mass spectrometry monitoring. Science China Chemistry 69, 1846–1853 (2026). DOI: 10.1007/s11426-025-2977-3. https://doi.org/10.1007/s11426-025-2977-3
22. Zhang, J. Doctoral dissertation, Virginia Tech (2013); Chapter 5 adopted from the Nature Chemistry 2013 manuscript (ref. 1). VTechWorks handle 10919/25290. https://vtechworks.lib.vt.edu/server/api/core/bitstreams/cbabd21b-8497-4e3a-965a-23af580c7a5a/content