Structure, Stability, and Cluster-Cage Interactions in Nitride Clusterfullerenes M3N@C2n (M = Sc, Y; 2n = 68−98): a Density Functional Theory Study

Substituted fullerenes as a promising capping ligand towards stabilization of exohedral Dy(III) based single-ion magnets: a theoretical study

Organometallic dysprosocenium-based molecular magnets are the forefront runners in offering giant magnetic anisotropy and blocking temperatures close to the boiling point of liquid nitrogen. Attaining linearity in the organometallic dysprosocenium complexes is the key to generating giant magnetic anisotropy and blocking barriers. In the present study, we have unravelled the coordination ability of the substituted fullerene (C55X5)- (where X = CCH3, B, and N) generated by fencing around the five-membered ring of fullerene towards stabilizing a new family of exohedral dysprosium organometallic complexes showcasing giant magnetic anisotropy and blockade barriers. Eight exohedral mononuclear dysprosium organometallic complexes, namely [Dy(η5-C55X5)(η4-C4H4)] (1), [Dy(η5-C55X5)(η5-Cp)]+ (2), [Dy(η5-C55X5)(η5-Cp*)]+ (3), [Dy(η5-C55X5)(η6-C6H6)]2+ (4), [Dy(η5-C55X5)(η8-C8H8)] (5), [Dy(η5-C55X5)2]+ (6) (where X = CCH3), [Dy(η5-C55B5)2]+ (7) and [Dy(η5-C55N5)2]+ (8), were studied using scalar relativistic density functional theory (SR-DFT) and the complete active space self-consistent field (CASSCF) methodology to shed light on the structure, stability, bonding and single-ion magnetic properties. SR-DFT calculations predict complexes 1-8 to be highly stable, with a strictly linear geometry around the Dy(III) ion in complexes 6-8. Energy Decomposition Analysis (EDA) predicts the following order for interaction energy (ΔEint value): 5 > 1 > 2 ≈ 3 > 6 > 7 > 8 > 4, with sizable 4f-ligand covalency in all the complexes. CASSCF calculations on complexes 1-8 predict stabilization of mJ |±15/2〉 as the ground state for all the complexes except for 5, with the following trend in the Ucal values: 6 (1573 cm-1) ≈ 3 (1569 cm-1) > 1 (1538 cm-1) > 8 (1347 cm-1) > 2 (1305 cm-1) > 7 (1284 cm-1) > 4 (1125 cm-1) > 5 (108 cm-1). Ab initio ligand field theory (AILFT) analysis provides a rationale for Ucal ordering, where π-type 4f-ligand interactions in complexes 1-4 and 6-8 offer giant barrier height while the large (C8H8)2- rings generate δ-type interaction in 5, which diminishes the axiality in the ligand field. Our detailed finding suggests that the exohedral organometallic dysprosocenium complexes are more linear compared to bent [DyCp*2]+ cations and display a giant barrier height exceeding 1500 cm-1 with negligible quantum tunnelling of magnetization (QTM) - a new approach to design highly anisotropic dysprosium organometallic complexes.

Thirty Years of Hide-and-Seek: Capturing Abundant but Elusive MIII@C3v(8)-C82 Isomer, and the Study of Magnetic Anisotropy Induced in Dy3+ Ion by the Fullerene π-Ligand

Our knowledge about endohedral metallofullerenes (EMFs) is restricted to the structures with sufficient kinetic stability to be extracted from the arc-discharge soot and processed by chromatographic and structural techniques. For the most abundant rare-earth monometallofullerene MIII@C82, experimental studies repeatedly demonstrated C2v(9) and Cs(6) carbon cage isomers, while computations predicted equal stability of the "missing" C3v(8) isomer. Here we report that this isomer is indeed formed but has not been recovered from soot using standard protocols. Using a combination of redox extraction and subsequent benzylation and trifluoromethylation with single-crystal XRD analysis of CF3 adduct, we prove that Dy@C3v(8)-C82 is one of the most abundantly produced metallofullerenes, which was not identified in earlier studies because of the low kinetic stability. Further, using the Dy@C3v(8)-C82(CF3) and Dy@C3v(8)-C82(CH2Ph) monoadducts for the case study, we analyzed the role of metal-fullerene bonding on the single-ion magnetic anisotropy of Dy in EMFs. The multitechnique approach, combining ab initio calculations, EPR spectroscopy, and SQUID magnetometry, demonstrated that coordination of the Dy ion to the fullerene cage induces moderate, nonaxial, and very fluid magnetic anisotropy, which strongly varies with small alterations in the Dy-fullerene coordination geometry. As a result, Dy@C3v(8)-C82(CH2Ph) is a weak field-induced single-molecule magnet (SMM), whose signatures of magnetic relaxation are detectable only below 3 K. Our results demonstrate that metal-cage interactions should have a detrimental effect on the SMM performance of EMFs. At the same time, the strong variability of the magnetic anisotropy with metal position suggests tunability and offers strategies for future progress.

Dual Stabilization of a Tri-Metallofullerene Radical Er3@C80: Exohedral Derivatization and Endohedral Three-Center Bonding

The enclosed space within fullerene molecules, capable of trapping metal clusters, offers an opportunity to investigate the behavior of metal atoms in a highly confined sub-nanometer environment. However, the studies on trimetallofullerenes M3@C80 have been very limited due to their difficult obtainability. In this paper, we present a new method for obtaining a tri-metallofullerene Er3@C80 through exohedral modification of the fullerene cage. Our findings reveal that Er3@C80 exhibits a radical character and can react with the dichlorobenzene radical to generate a stable derivative Er3@C80PhCl2. Theoretical calculations demonstrate the presence of a three-center two-electron metal-metal bond in the center of the fullerene cage. This bond serves to counterbalance the Coulomb repulsion between the Er ions. Consequently, both exohedral derivatization and endohedral three-center bonding contribute to the substantial stability of Er3@C80PhCl2. Furthermore, molecular dynamics simulations indicate that the Er3 cluster within the molecule possesses a rigid triangle structure. The availability of M3@C80 derivatives opens avenues for future investigations into interactions among metal atoms, such as magnetic coupling, within fullerene cages.

Stability and Electronic Properties of Mixed Rare-Earth Tri-Metallofullerenes YxDy3-x@C80 (x = 1 or 2)

Tri-metallofullerenes, specifically M3@C80 where M denotes rare-earth metal elements, are molecules that possess intriguing magnetic properties. Typically, only one metal element is involved in a given tri-metallofullerene molecule. However, mixed tri-metallofullerenes, denoted as M1xM23-x@C80 (x = 1 or 2, M1 and M2 denote different metal elements), have not been previously discovered. The investigation of such mixed tri-metallofullerenes is of interest due to the potential introduction of distinct properties resulting from the interaction between different metal atoms. This paper presents the preparation and theoretical analysis of mixed rare-earth tri-metallofullerenes, specifically YxDy3-x@C80 (x = 1 or 2). Through chemical oxidation of the arc-discharge produced soot, the formation of tri-metallofullerene cations, namely Y2Dy@C80+ and YDy2@C80+, has been observed. Density functional theory (DFT) calculations have revealed that the tri-metallofullerenes YxDy3-x@C80 (x = 1 or 2) exhibit a low oxidation potential, significantly lower than other fullerenes such as C60 and C70. This low oxidation potential can be attributed to the relatively high energy level of a singly occupied orbital. Additionally, the oxidized species demonstrate a large HOMO-LUMO gap similar to that of YxDy3-xN@C80, underscoring their high chemical stability. Theoretical investigations have uncovered the presence of a three-center two-electron metal-metal bond at the center of Y2DY@C80+ and YDy2@C80+. This unique multi-center bond assists in alleviating the electrostatic repulsion between the metal ions, thereby contributing to the overall stability of the cations. These mixed rare-earth tri-metallofullerenes hold promise as potential candidates for single-molecule magnets.

Dynamic Metastable Characteristics of Carbon Cages Embedded with Er2C2

Novel endohedral metallofullerenes (EMFs), namely, Er2C2@C2v(5)-C80, Er2C2@Cs(6)-C82, Er2C2@Cs(15)-C84, Er2C2@C2v(9)-C86, Er2C2@Cs(15)-C86, and Er2C2@Cs(32)-C88, had been experimentally synthesized, and the unique structures and many fascinating properties had also been widely explored. Nevertheless, the position of the Er atoms inside the cage shows a severe disorder within the stable EMF monomer, which is difficult to understand and explain from the experimental point of view. In this work, based on the density functional theoretical calculations, the Er2C2@Cs(6)-C82 has 73 directional isomers and 2 Er atoms that are far beyond from Er-Er single bonding and tend to be close to the cage side (marked as "shell"), and the core (Er2C2 units) takes on a butterfly shape as generally revealed. The energy difference between any two of the isomers is in the range of 0.05 to 25.6 kcal/mol, indicating a relatively easy thermodynamic transition between the isomers. The other five Er carbide cluster EMFs (Er2C2@C2v(5)-C80, Er2C2@Cs(15)-C84, Er2C2@C2v(9)-C86, Er2C2@Cs(15)-C86, and Er2C2@Cs(32)-C88) are also studied in the same way, and 30, 37, 39, and 43 most stable Er-oriented sites inside the cage, respectively, are obtained. In addition, the shape of the Er2C2 gradually changed from butterfly to linear. Moreover, the electronic structure and molecular orbital analyses show that it is easy for Er2C2@C80-88 to form a charge transfer state of [Er2C2]4+@[C80-88]4- via the dynamic core-shell coordination equilibrium. Er2C2 with a steep drop in chemical stability is restricted to forming varying degrees of metastable states in the shell, determined by the shell size, to ensure the overall stability. The lowest unoccupied molecular orbital energy level of these EMFs is increased by 0.5-1.1 eV compared with fullerenes C80-88, potentially providing favorable conditions for suitable energy level matching with EMF as an electron acceptor used in organic solar cell devices.

Two Metastable Endohedral Metallofullerenes Sc2C2@C1(39656)-C82 and Sc2C2@C1(51383)-C84: Direct-C2-Insertion Products from Their Most Stable Precursors

Endohedral metallofullerenes (EMFs) are sub-nano carbon materials with diverse applications, yet their formation mechanism, particularly for metastable isomers, remains ambiguous. The current theoretical methods focus mainly on the most stable isomers, leading to limited predictability of metastable ones due to their low stabilities and yields. Herein, we report the successful isolation and characterization of two metastable EMFs, Sc₂C₂@C₁(39656)-C₈₂ and Sc₂C₂@C₁(51383)-C₈₄, which violate the isolated pentagon rule (IPR). These two non-IPR EMFs exhibit a rare case of planar and pennant-like Sc₂C₂ clusters, which can be considered hybrids of the common butterfly-shaped and linear configurations. More importantly, the theoretical results reveal that despite being metastable, these two non-IPR EMFs survived as the products from their most stable precursors, Sc₂C₂@C_2v(5)-C₈₀ and Sc₂C₂@C_s(6)-C₈₂, via a C₂ insertion during the post-formation annealing stages. We propose a systematic theoretical method for predicting metastable EMFs during the post-formation stages. The unambiguous molecular-level structural evidence, combined with the theoretical calculation results, provides valuable insights into the formation mechanisms of EMFs, shedding light on the potential of post-formation mechanisms as a promising approach for EMF synthesis.

Functionalized Fullerenes and Their Applications in Electrochemistry, Solar Cells, and Nanoelectronics

Carbon-based nanomaterials have rapidly advanced over the last few decades. Fullerenes, carbon nanotubes, graphene and its derivatives, graphene oxide, nanodiamonds, and carbon-based quantum dots have been developed and intensively studied. Among them, fullerenes have attracted increasing research attention due to their unique chemical and physical properties, which have great potential in a wide range of applications. In this article, we offer a comprehensive review of recent progress in the synthesis and the chemical and physical properties of fullerenes and related composites. The review begins with the introduction of various methods for the synthesis of functionalized fullerenes. A discussion then follows on their chemical and physical properties. Thereafter, various intriguing applications, such as using carbon nanotubes as nanoreactors for fullerene chemical reactions, are highlighted. Finally, this review concludes with a summary of future research, major challenges to be met, and possible solutions.

Self-driven carbon atom implantation into fullerene embedding metal-carbon cluster

Hundreds of members have been synthesized and versatile applications have been promised for endofullerenes (EFs) in the past 30 y. However, the formation mechanism of EFs is still a long-standing puzzle to chemists, especially the mechanism of embedding clusters into charged carbon cages. Here, based on synthesis and structures of two representative vanadium-scandium-carbido/carbide EFs, VSc₂C@I_h (7)-C₈₀ and VSc₂C₂@I_h (7)-C₈₀, a reasonable mechanism-C₁ implantation (a carbon atom is implanted into carbon cage)-is proposed to interpret the evolution from VSc₂C carbido to VSc₂C₂ carbide cluster. Supported by theoretical calculations together with crystallographic characterization, the single electron on vanadium (V) in VSc₂C@I_h (7)-C₈₀ is proved to facilitate the C₁ implantation. While the V=C double bond is identified for VSc₂C@I_h (7)-C₈₀, after C₁ implantation the distance between V and C atoms in VSc₂C₂@I_h (7)-C₈₀ falls into the range of single bond lengths as previously shown in typical V-based organometallic complexes. This work exemplifies in situ self-driven implantation of an outer carbon atom into a charged carbon cage, which is different from previous heterogeneous implantation of nonmetal atoms (Group-V or -VIII atoms) driven by high-energy ion bombardment or high-pressure offline, and the proposed C₁ implantation mechanism represents a heretofore unknown metal-carbon cluster encapsulation mechanism and can be the fundamental basis for EF family genesis.

Cluster-Geometry-Associated Metal-Metal Bonding in Trimetallic Carbide Clusterfullerenes

Geometry configurations of the metallic clusters play a significant role in the involved bonding nature. Herein, we report the crystallographic characterization of unprecedented erbium-based trimetallic clusterfullerenes, namely, Er₃C₂@I_h(7)-C₈₀, in which the inner Er₃C₂ cluster presents a lifted bat ray configuration with the C₂ unit elevated by ∼1.62 Å above the Er₃ plane. Within the plane, the Er···Er distances for Er1···Er2, Er1···Er2A, and Er2···Er2A are 3.4051(15), 3.4051(15), and 3.3178(15) Å, respectively, falling into the range of the metal-metal bonding. Density functional theory calculations unveil the three-center-one-electron Er-Er-Er bond in Er₃C₂@I_h(7)-C₈₀ with one electron shared by three metals, and thus, its exceptional electronic structure can be expressed as (Er₃)⁸⁺(C₂)^2-@C₈₀^6-. Interestingly, with the further observation on the geometry configurations of the encapsulated clusters in M₃C₂@C_2n (M = Sc, Y, and Lu) series, we find that the lifted bat ray configuration of the inner cluster is explicitly associated with the formation of the bonding interactions between the inner metals. This finding provides insights into the nature of metal-metal bonding and gives guidelines for the design of the single-molecule magnet.

Overlooked Effects of La-4f Orbitals in Endohedral Metallofullerenes

For endohedral metallofullerenes (EMFs), a central issue is how to correctly describe the intracluster and metal-cage interactions, which are critical for understanding their structures, stabilities, and various properties. In this work, density functional theory calculations were carried out for 13 La-based EMFs covering all four reported types and a rather wide cage size range (C₃₂-C₁₀₄). The results reveal that the usually core-like lanthanide 4f subshell may play a critical role in the structural characteristics, energetic stabilities, frontier orbital energy levels, metal charges, and chemical reactivities of these endofullerenes. Regardless of the encapsulated forms, the La-4f contributions to the chemical bonding and structural stability increase with the reduced cage sizes because of the gradually enhanced cage confinement. The combination of metal-to-nonmetal charge transfer and compression of the cage cavity exposes and effectively activates the otherwise chemically inert 4f orbitals. By disclosing the important role of long-neglected metal orbitals inside fullerenes, the current work not only deepens our understanding of EMFs, but also provides new insights into the chemical bondings in general confined spaces at the subnanometer scale.

Orientational Effects on the Electronic Structure and Polarization in Sc3N@C80

Of particular interest in metal encapsulating fullerenes such as Sc₃N@C₈₀ is not just how the charge is transferred between atoms of the metal nitride core and the carbon cage but how the orientation of the core impacts the electronic structure of the entire molecule. A channel for the electron backdonation is identified which leads to a charge hole on the fullerene cage, just below the valence level, which is pinned to the orientation of the metal nitride. This electron hole is overcompensated by paired electron charge at deeper levels. Reorienting the metal nitride inside the cage results in a change in the surface charge distribution below the valence level and the lowest energy empty levels coupled to the metal ions. The charge separation between the metal nitride core and the carbon cage drives changes in the polarizability of the molecule depending on the orientation of the core. Furthermore, it is found that the electron affinity depends on the orientation of the metal nitride core. This is the result of the overlap between the scandium ions' d-orbitals with the fullerene cage 6:6:6 p orbitals. The electronic and geometric structures of the Sc₃N@C₈₀ metallofullerene are examined by using the density functional theory, and the findings are corroborated by an analysis of the density of states combined with charge density plots.

Ln3@C80+ (Ln = lanthanide): a new class of stable metallofullerene cations with multicenter metal-metal bonding in sub-nanometer confined space

Among the large number of members in the metallofullerene family, the nitride clusterfullerene M3N@C80 (M = trivalent metal) is a special one with unordinary high stability. It is generally thought...

Decisive role of non-rare earth metals in high-regioselectivity addition of μ3-carbido clusterfullerene

The reaction of μ3-CCF Dy2TiC@Ih-C80 with AdN2 affords only one [6,6]-open monoadduct along with the addition sites adjacent to the Ti4+ ion instead of the two Dy3+ ions, revealing the decisive role of the non-rare earth metal Ti(IV).

A spider hanging inside a carbon cage: off-center shift and pyramidalization of Sc3N clusters inside C84 and C86 fullerene cages

Structural analysis shows that, in Sc3N@Cs(51365)-C84 and Sc3N@D3(19)-C86, the Sc3N clusters are shifted to one side of the cages and unexpectedly pyramidalized inside the large cages of C84 and C86, which resembles a spider attached to a web.

Crystallographic Characterization of U@C2n (2n = 82-86): Insights about Metal-Cage Interactions for Mono-metallofullerenes

Endohedral mono-metallofullerenes are the prototypes to understand the fundamental nature and the unique interactions between the encapsulated metals and the fullerene cages. Herein, we report the crystallographic characterizations of four new U-based mono-metallofullerenes, namely, U@C_s(6)-C₈₂, U@C₂(8)-C₈₄, U@C_s(15)-C₈₄, and U@C₁(12)-C₈₆, among which the chiral cages C₂(8)-C₈₄ and C₁(12)-C₈₆ have never been previously reported for either endohedral or empty fullerenes. Symmetrical patterns, such as indacene, sumanene, and phenalene, and charge transfer are found to determine the metal positions inside the fullerene cages. In addition, a new finding concerning the metal positions inside the cages reveals that the encapsulated metal ions are always located on symmetry planes of the fullerene cages, as long as the fullerene cages possess mirror planes. DFT calculations show that the metal-fullerene motif interaction determines the stability of the metal position. In fullerenes containing symmetry planes, the metal prefers to occupy a symmetrical arrangement with respect to the interacting motifs, which share one of their symmetry planes with the fullerene. In all computationally analyzed fullerenes containing at least one symmetry plane, the actinide was found to be located on the mirror plane. This finding provides new insights into the nature of metal-cage interactions and gives new guidelines for structural determinations using crystallographic and theoretical methods.

Th@D5h(6)-C80: a highly symmetric fullerene cage stabilized by a single metal ion

A novel endohedral metallofullerene (mono-EMF), Th@D5h(6)-C80, has been successfully synthesized and fully characterized by mass spectrometry, single crystal X-ray diffraction, UV-vis-NIR and Raman spectroscopy and cyclic voltammetry. Single crystal XRD analysis unambiguously assigned the fullerene cage as D5h(6)-C80, the first example in which the highly symmetric cage is stabilized by a single metal ion. The combined experimental and theoretical studies further reveal that the D5h(6)-C80 cage, known only for its stabilization by 6-electron transfer, is stabilized by the 4-electron transfer from the encapsulated Th ion for the first time.

Structural, Electronic, and Nonlinear Optical Properties of C66H4 and C70Cl6 Encapsulating Li and F Atoms

Recently, nonclassical fullerene derivatives C66H4 and C70Cl6, which both contain two negatively curved moieties of heptagons, have been successfully synthesized. Inspired by these experimental achievements, the structural and electronic properties of C66H4, C70Cl6, Li@C66H4, F@C66H4, Li@C70Cl6, and F@C70Cl6 were systematical studied through density functional theory calculations in this work. Our results show that the reduction of the front molecular orbital gap of fullerene derivatives occurs with the introduction of Li and F atoms. After quantitative analysis of back-donations of charge between an encapsulated atom and an external carbon cage, it is found that C66H4 and C70Cl6 prefer to act as electron acceptors. It is interesting to note that the strong covalent nature of the interactions between a F atom and a carbon cage is observed, whereas the weak covalent and strong ionic interactions occur between a Li atom and a carbon cage. On the other hand, according to the first hyperpolarizability results, the encapsulation of the Li atom enhances the nonlinear optical response of fullerene derivatives. This work provides a strategy to improve nonlinear optical properties of C66H4 and C70Cl6, reveals the internal mechanism of the contribution from Li and F atoms to endohedral fullerene derivatives, and will contribute to the designation of endohedral fullerene derivative devices.

A non-isolated pentagon rule C82 cage stabilized by a stretched Sc3N cluster

A novel cluster fullerene, Sc3N@Cs(39 663)-C82, has been synthesized and characterized. Crystallograpic charaterization unambiguously determines the non-IPR cage structure of Cs(39 663)-C82. Structural analyses further reveal that the Sc3N cluster is notably stretched to facilitate the interaction with the non-IPR fullerene cage.

A cyclic manipulation of cage isomers via anion exchange and thermal isomerism

A cyclic manipulation of peanut cage isomers has been achieved via anion exchange and unusual cage isomerism.

Significant Roles of a Particularly Stable Two-Center Two-Electron Lu-Lu σ Bond in Lu2@C86: Electronic Structure of Lu and Radius of Lu2

There is still dispute over the stability of endohedral metallofullerenes (EMFs) M₂C_2n, and recently, multiform lutetium-based dimetallofullerenes have been dropped in experiments. The thermodynamic stabilities of Lu₂C₈₆ EMFs are revealed by density functional theory (DFT) in conjunction with statistical thermodynamic analyses. Inevitably, besides the experimentally reported Lu₂@C_2v(63751)-C₈₆, Lu₂@C_s(63750)-C₈₆, and Lu₂@C_s(63757)-C₈₆, other three metal carbide clusterfullerenes, Lu₂C₂@D_2d(51591)-C₈₄, Lu₂C₂@C₁(51383)-C₈₄, and Lu₂C₂@C_s(id207430)-C₈₄, rather than Lu₂@C₈₆ are first characterized as thermodynamically stable isomers of Lu₂C₈₆. Specially, the C_s(id207430)-C₈₄ is a newly non-classical fullerene containing one heptagon, which is stabilized via encaging Lu₂C₂. Another interesting phenomenon is that the outer fullerene cages of thermodynamically stable Lu₂C_82-88 molecules are geometrically connected through C₂ addition/loss and Stone-Wales (SW) transformation, suggesting a special relationship between thermodynamic stabilities and geometries of Lu₂C_82-88 EMFs. Furthermore, the electronic configurations of (Lu₂)⁴⁺@C₈₆^4- and (Lu₂C₂)⁴⁺@C₈₄^4- were confirmed. A significantly stable two-center two-electron (2c-2e) Lu-Lu σ single bond is formed in Lu₂@C₈₆. By comparing M-M bonds in M₂@C_2v(63751)-C₈₆ (M = Sc, Y, La, and Lu), two significant factors, the valence atomic orbital (ns) of metal atoms and radius of M²⁺, are found to determine the stability of the M-M bond in the C_2v(63751)-C₈₆. Additionally, the simulated UV-vis-NIR spectra of thermodynamically stable Lu₂C₈₆ isomers were simulated, which further disclose their electronic features.

Crystallographic Characterization of Ti2C2@D3h(5)-C78, Ti2C2@C3v(8)-C82, and Ti2C2@Cs(6)-C82: Identification of Unsupported Ti2C2 Cluster with Cage-Dependent Configurations

Fullerene cages are ideal hosts to encapsulate otherwise unstable metallic clusters to form endohedral metallofullerenes (EMFs). Herein, a novel Ti₂C₂ cluster with two titanium atoms bridged by a C₂-unit has been stabilized by three different fullerene cages to form Ti₂C₂@D_3h(5)-C₇₈, Ti₂C₂@C_3v(8)-C₈₂, and Ti₂C₂@C_s(6)-C₈₂, representing the first examples of unsupported titanium carbide clusters. Crystallographic results show that the configuration of the Ti₂C₂ cluster changes upon cage variation. In detail, the Ti₂C₂ cluster adopts a butterfly shape in Ti₂C₂@C_3v(8)-C₈₂ and Ti₂C₂@C_s(6)-C₈₂ with Ti-C₂-Ti dihedral angles of 156.35 and 147.52° and Ti-Ti distances of 3.633 and 3.860 Å, respectively. In sharp contrast, a stretched planar geometry of Ti₂C₂ is observed in Ti₂C₂@D_3h(5)-C₇₈, where a Ti-C₂-Ti angle of 176.87° and a long Ti-Ti distance of 5.000 Å are presented. Consistently, theoretical calculations reveal that the cluster configuration is very sensitive to the cage shape which eventually determines the electronic structures of the hybrid EMF-molecules, thus adding new insights into modern coordination chemistry.

Addition of CF2 group to endohedral fullerene Sc3N@Ih-C80

We report the first successful synthesis of a CF₂ derivative of the stable endohedral fullerene Sc₃N@I_h-C₈₀. Reaction with CF₂ClCOONa yields a single C_s-symmetric Sc₃N@C₈₀(CF₂) adduct where the CF₂ group is inserted into a [6,6]-bond and opens it to 2.3 Å between the bridgehead carbon atoms. As evidenced by absorption and fluorescence spectroscopy as well as cyclic voltammetry, both the HOMO and the LUMO level of Sc₃N@C₈₀(CF₂) are slightly (ca. 0.1 eV) downshifted with respect to the parent Sc₃N@I_h-C₈₀, so the HOMO-LUMO gap remains essentially unchanged. The DFT calculations suggest that the reaction mechanism is not the previously assumed [2 + 1]-cycloaddition of :CF₂ carbene but rather nucleophilic addition of CF₂Cl^- anion followed by elimination of Cl^- and closing of the CF₂ bridge via intramolecular nucleophilic substitution. Selective formation of the [6,6]-Sc₃N@C₈₀(CF₂) turns out to be kinetically controlled and promoted by a particular orientation of the endohedral Sc₃N cluster with respect to the CF₂Cl^- addition site. In its turn, the CF₂ addend partly hampers the rotation of Sc₃N the endohedral cluster compared to its quasi-free reorientations in the parent Sc₃N@I_h-C₈₀.

Intermolecular Coulombic Decay in Endohedral Fullerene at the 4d→4f Resonance

Intermolecular processes offer unique decay mechanisms for complex systems to internally relax. Here, we report the observation of an intermolecular Coulombic decay channel in an endohedral fullerene, a holmium nitride complex (Ho_{3}N) embedded within a C_{80} fullerene, between neighboring holmium ions, and between the holmium complex and the carbon cage. By measuring the ions and the electrons in coincidence after XUV photoabsorption, we can isolate the different decay channels, which are found to be more prevalent relative to intra-atomic Auger decay.

Th@C86, Th@C82, Th@C80, and Th@C76: role of thorium encapsulation in determining spherical aromatic and bonding properties on medium-sized endohedral metallofullerenes

Thorium encapsulated metallofullerenes (Th-EMFs) with external C76, C80, C82, and C86 cages have been synthesized, with the 13C-NMR spectrum recorded for Th@C82. Here, we explore computationally the chemical bonding, NMR and spherical aromaticity of Th@C82 and related thorium-encapsulated metallofullerenes. Our results show that these Th-EMFs are new examples of spherical aromatic structures, representing interesting low-symmetry exceptions to the Hirsch 2(N + 1)2 rule of spherical aromaticity. Their electronic structures are based on π-electron counts of 80, 84, 86, and 90, respectively, with a shell structure ranging from S2P6D10F14G18H22I8 to S2P6D10F14G18H22I18, where the partially filled I-shell remains as a frontier orbital. Their behavior is comparable to that of the spherical aromatic alkali-C606- phases, which in addition to the favorable endohedral Th-fullerene bonding account for their particular abundance exhibiting the ability to sustain a long-range shielding cone as a result of the favorable metal-cage bonding. This rationalization of such species as neutral spherical aromatic EMFs suggests the possibility of an extensive series of aromatic fullerenes with nuclearity larger than C60 buckminsterfullerene as stable building blocks towards nanostructured metal-organic materials.

Comparing Empty and Filled Fullerene Cages with High-Resolution Trapped Ion Mobility Spectrometry

We have used trapped ion mobility spectrometry (TIMS) to obtain highly accurate experimental collision cross sections (CCS) for the fullerene C₈₀^- and the endohedral metallofullerenes La₂@C₈₀^-, Sc₃N@C₈₀^-, and Er₃N@C₈₀^- in molecular nitrogen. The CCS values of the endohedral fullerenes are 0.2% larger than that of the empty cage. Using a combination of density functional theory and trajectory calculations, we were able to reproduce these experimental findings theoretically. Two effects are discussed that contribute to the CCS differences: (i) a small increase in fullerene cage size upon endohedral doping and (ii) charge transfer from the encapsulated moieties to the cage thus increasing the attractive charge-induced dipole interaction between the (endohedral) fullerene ion and the nitrogen bath gas molecules.

Isolation and Crystallographic Characterization of Two, Nonisolated Pentagon Endohedral Fullerenes: Ho3 N@C2 (22010)-C78 and Tb3 N@C2 (22010)-C78

Purified samples of Ho₃ N@C₂ (22010)-C₇₈ and Tb₃ N@C₂ (22010)-C₇₈ have been isolated by two distinct processes from the rich array of fullerenes and endohedral fullerenes present in carbon soot from graphite rods doped with Ho₂ O₃ or Tb₄ O₇ . Crystallographic analysis of the endohedral fullerenes as cocrystals with Ni(OEP) (in which OEP is the dianion of octaethylporphyrin) shows that both molecules contain the chiral C₂ (22010)-C₇₈ cage. This cage does not obey the isolated pentagon rule (IPR) but has two sites where two pentagons share a common C-C bond. These pentalene units bind two of the metal ions, whereas the third metal resides near a hexagon of the cage. Inside the cages, the Ho₃ N or Tb₃ N unit is planar. Ho₃ N@C₂ (22010)-C₇₈ and Tb₃ N@C₂ (22010)-C₇₈ use the same cage previously found for Gd₃ N@C₂ (22010)-C₇₈ rather than the IPR-obeying cage found in Sc₃ N@D_3h -C₇₈ .

U2C Unit in Fullerenes: Robust Multicenter Bonds with a Cluster Shape Controlled by Cage Size and Charge Transfer

Stimulated by the recent successful synthesis and crystallographic characterization of the first diuranium carbide endohedral metallofullerene (EMF) U₂C@I_h(7)-C₈₀ ( Zhang et al. Nat. Commun. ; 2018 ), density functional theory calculations were performed for a series of U₂C@C_2n (2n = 60, 68, 72, 78, 80, 88, 96, and 104) analogues. The internal UCU bond angle increases from 96.9° in I_h-C₆₀ to 180.0° in D_3d-C₁₀₄, exhibiting cage-size-dependent cluster configuration. However, further evidence suggests that the U₂C shape may be also affected by the amount of charge transferred from the cluster to the outer cage with 6e and 4e favoring bent and linear, respectively. The change of the bond angle closely correlates with the charge and hybrid state of the internal atom. Significantly, besides the covalent two-center two-electron (2c-2e) U-C bonds, the U₂C unit always features two 3c-2e bonds regardless of its size, shape, and charge state. Furthermore, for the cluster-cage interactions, besides the dominated electrostatic attractions, all these EMFs show an obvious covalent character with the substantial participation of U 5f valence orbitals.

Crystallographic characterization of Er3N@C2n (2n = 80, 82, 84, 88): the importance of a planar Er3N cluster

A series of Er-based nitride clusterfullerenes (NCFs), Er₃N@C_80-88, have been successfully synthesized and isolated. In particular, Er₃N@I_h(7)-C₈₀, Er₃N@D_5h(6)-C₈₀, Er₃N@C_2v(9)-C₈₂, Er₃N@C_s(51365)-C₈₄, and Er₃N@D₂(35)-C₈₈ have been characterized by single-crystal X-ray diffraction (XRD) for the first time. The planar configuration of the inserted Er₃N cluster is identified unambiguously and the Er-N distances increase in accordance with cage expansion to maintain strong metal-cage interactions. Additionally, the electrochemical properties of the Er₃N@C_80-88 series are studied by means of cyclic voltammetry. It is found that the first reduction potentials are roughly similar for all compounds under study, while the first oxidation potentials are cathodically shifted along with the increase of the cage size in the Er₃N@C_2n (2n = 80, 84, 86, 88) series, leading to a decrease in the corresponding electrochemical band gaps. Nevertheless, for Er₃N@C_2v(9)-C₈₂, a good electron donating ability is manifested by its relatively small first oxidation potential, which results from the relatively higher energy level of the highest occupied molecular orbital. The redox behaviors observed in such Er₃N-based NCFs may promise their great potential applications in donor-acceptor systems.

Mixed dysprosium-lanthanide nitride clusterfullerenes DyM2N@C80-Ih and Dy2MN@C80-Ih (M = Gd, Er, Tm, and Lu): synthesis, molecular structure, and quantum motion of the endohedral nitrogen atom

Systematic exploration of the synthesis of mixed-metal Dy-M nitride clusterfullerenes (NCFs, M = Gd, Er, Tm, Lu) is performed, and the impact of the second metal on the relative yield is evaluated. We demonstrate that the ionic radius of the metal appears to be the main factor allowing explanation of the relative yields in Dy-M mixed-metal systems with M = Sc, Lu, Er, and Gd. At the same time, Dy-Tm NCFs show anomalously low yields, which is not consistent with the relatively small ionic radius of Tm³⁺ but can be explained by the high third ionization potential of Tm. Complete separation of Dy-Gd and Dy-Er, as well as partial separation of Dy-Lu M₃N@C₈₀ nitride clusterfullerenes, is accomplished by recycling HPLC. The molecular structures of DyGd₂N@C₈₀ and DyEr₂N@C₈₀ are analyzed by means of single-crystal X-ray diffraction. A remarkable ordering of mixed-metal nitride clusters is found despite similar size and electronic properties of the metals. Possible pyramidalization of the nitride clusters in these and other nitride clusterfullerenes is critically analyzed with the help of DFT calculations and reconstruction of the nitrogen inversion barrier in M₃N@C₈₀ molecules is performed. Although a double-well potential with a pyramidal cluster structure is found to be common for most of them, the small size of the inversion barrier often leads to an apparent planar structure of the cluster. This situation is found for those M₃N@C₈₀ molecules in which the energy of the lowest vibrational level exceeds that of the inversion barrier, including Dy₃N@C₈₀ and DyEr₂N@C₈₀. The genuine pyramidal structure can be observed by X-ray diffraction only when the lowest vibrational level is below the inversion barrier, such as those found in Gd₃N@C₈₀ and DyGd₂N@C₈₀. The quantum nature of molecular vibrations becomes especially apparent when the size of the inversion barrier is comparable to the energy of the lowest vibrational levels.

Hybrid ADFT Study of the C104 and C106 IPR Isomers

This work presents a hybrid auxiliary density functional theory (ADFT) study of the neutral and hexaanionic C₁₀₄ and C₁₀₆ fullerenes with the aim to determine their ground state structures. To this end, all C₁₀₄ and C₁₀₆ fullerene structures that obey the isolated pentagon rule (IPR) were optimized with the Perdew-Burke-Ernzerhof generalized gradient approximation followed by a single-point energy calculation with the PBE0 hybrid functional. Our studies show that this composite approach yields relative energies of giant fullerenes that are accurate to around 1 kcal/mol. As a result, the ground states of C₁₀₄, C₁₀₄^6-, and C₁₀₆^6- can be assigned to the isomers 234:C_s, 821:D₂, and 891:C_s, respectively. On the other hand, the energetically lowest lying IPR isomers of C₁₀₆, 331:C_s, 1194:C₂, 534:C₁ are separated by less than 1 kcal/mol which makes an unequivocal ground state assignment by hybrid DFT methods impossible. To guide future experiments, we also report the simulated IR and Raman spectra of the most stable neutral and hexaanionic C₁₀₄ and C₁₀₆ fullerenes.

Aromaticity, Coulomb repulsion, π delocalization or strain: who is who in endohedral metallofullerene stability?

Three different models for endohedral metallofullerene structure prediction are compared, revealing the physical origin of the stability of these compounds.

Anion-encapsulating fullerenes behave as large anions: a DFT study

M06L/6-311++G(d,p)//M06L/6-31G(d,p) level density functional theory studies show that the endohedral reaction of C60 with X- (X = F, Cl, Br, OH, NH2, NO2, CN, and ClO) is exothermic by 37.8-65.2 kcal mol-1. The exothermic character of the reaction is drastically reduced in polar and nonpolar solvents due to the lack of direct solvation influence on the encapsulated anion. In all X-@C60, the occupied frontier molecular orbitals (FMOs) are located on X- while the energy levels of FMOs centered on C60 are very similar to those of the C60- radical anion. Molecular electrostatic potential (MESP) analysis of X-@C60 revealed that the negative character of the MESP minimum (Vmin) on the carbon cage increases by ∼72 fold compared to C60, which is very similar to the enhancement in the negative MESP observed on the C60- radical anion. The MESP data and quantum theory of atoms in molecules (QTAIM) analysis of charge, electron delocalization index, and Laplacian of bond critical point (bcp) support significant electron sharing from the anion to the carbon atoms of the fullerene cage, which makes the cage behave like a very large anion in a closed shell configuration. The data are also supportive of a multicenter charge-shift type of bonding interaction between the anion and the carbon cage. The anionic nature of the fullerene cage has been verified in the cases of larger systems such as Cl-@C70, Cl-@C84, and Cl-@C90. The binding of a counter cation K+ with X-@C60 is found to be highly exothermic (∼72 kcal mol-1) and very similar to the binding of K+ with the C60- radical anion (72.9 kcal mol-1), which suggests that C60 in X-@C60 behaves as a closed shell anion.

When metal clusters meet carbon cages: endohedral clusterfullerenes

Fullerenes have the characteristic of a hollow interior, and this unique feature triggers intuitive inspiration to entrap atoms, ions or clusters inside the carbon cage in the form of endohedral fullerenes. In particular, upon entrapping an otherwise unstable metal cluster into a carbon cage, the so-called endohedral clusterfullerenes fulfil the mutual stabilization of the inner metal cluster and the outer fullerene cage with a specific isomeric structure which is often unstable as an empty fullerene. A variety of metal clusters have been reported to form endohedral clusterfullerenes, including metal nitrides, carbides, oxides, sulfides, cyanides and so on, making endohedral clusterfullerenes the most variable and intriguing branch of endohedral fullerenes. In this review article, we present an exhaustive review on all types of endohedral clusterfullerenes reported to date, including their discoveries, syntheses, separations, molecular structures and properties as well as their potential applications in versatile fields such as biomedicine, energy conversion, and so on. At the end, we present an outlook on the prospect of endohedral clusterfullerenes.

Formation of Spherical Aromatic Endohedral Metallic Fullerenes. Evaluation of Magnetic Properties of M@C28 (M = Ti, Zr, and Hf) from DFT calculations

The small C₂₈ cage has been shown experimentally to encapsulate titanium, zirconium, and hafnium (M), among other elements. Here, we explore computationally its magnetic response properties accounting for both global and local shielding tensors. Our results exhibit a continuous shielding region for M@C₂₈ for an orientation-averaged applied field thereby differing from that observed for the hollow C₂₈ structure. Moreover, under a specific orientation of the applied field a long-ranged shielding cone is obtained supporting the spherical aromatic behavior expected by the 2(N + 1)² Hirsch rule for M@C₂₈, standing for its particular abundance. The comparison between the hollow and endohedral C₂₈ fullerenes exhibits a characteristic long-range behavior at the outside region of the structure. The particular shape of the local chemical shift anisotropy tensor at a representative carbon atom exhibits inherent patterns as a consequence of the spherical aromatic behavior. This shows the capabilities from NMR experiments to account for the nonaromatic → aromatic variation. We envisage that the current approach will be beneficial in comparative studies of aromatic and electronic structure properties, to gain a deeper understanding of the geometrical and electronic structure situation in other endohedral species beyond that available from the information provided by routine NMR measurements.