A diuranium carbide cluster stabilized inside a C80 fullerene cage.

Unsupported non-bridged uranium-carbon double bonds have long been sought after in actinide chemistry as fundamental synthetic targets in the study of actinide-ligand multiple bonding. Here we report that, utilizing Ih(7)-C80 fullerenes as nanocontainers, a diuranium carbide cluster, U=C=U, has been encapsulated and stabilized in the form of UCU@Ih(7)-C80. This endohedral fullerene was prepared utilizing the Krätschmer-Huffman arc discharge method, and was then co-crystallized with nickel(II) octaethylporphyrin (NiII-OEP) to produce UCU@Ih(7)-C80·[NiII-OEP] as single crystals. X-ray diffraction analysis reveals a cage-stabilized, carbide-bridged, bent UCU cluster with unexpectedly short uranium-carbon distances (2.03 Å) indicative of covalent U=C double-bond character. The quantum-chemical results suggest that both U atoms in the UCU unit have formal oxidation state of +5. The structural features of UCU@Ih(7)-C80 and the covalent nature of the U(f1)=C double bonds were further affirmed through various spectroscopic and theoretical analyses.

Introduction

Understanding the nature of actinide-ligand multiple bonding remains a modern day challenge owing to the complex electronic structures of these elements, and as a consequence, our comprehension of their chemistries lags behind that of the more commonly studied transition metals[1-5]. Importantly, recent developments in the synthesis and study of molecular uranium complexes containing a variety of U=E(L) or U≡E(L) moieties (E = O, S, Se, Te, N, NR, P, As; L = other ligands) have provided valuable insights regarding the participation of the 5f and 6d valence orbitals in chemical bonding[1,2,6-27]. Conspicuously missing from this list are U=C bonds, specifically those which lack ancillary heteroatom or chelating support. While metal–carbon double bonds (i.e., Schrock carbenes) are common in transition metal chemistry and catalysis, unsupported U=C bonds have long remained a major and outstanding synthetic target. Some success towards the synthesis of An=C bonds has been realized through the use of heteroatom stabilizing chelating ligands such as U=[C(Ph2PNSiMe3)2](O)Cl2 or U=[C(Ph2PS)] (BH4)2 (THF)2, among others[4,27-33]. In a few rare instances, the isolation of non-chelated U=C bonds has been achieved in complexes such as (η5-C5H5)3U=CHPMe2Ph and [N(SiMe3)2]3U=CHPPh3, using formally dianionic ylidic ligands[31,34]. Electronically, the U=C units found in these methanediide ([κ3-CL2]2−) or ylidic ([C(H)(PR3)]2−) examples are best described as highly polarized, nucleophilic carbenes with σ and modest π bonding overlap between uranium and carbon. The α-carbon bonded heteroatom(s) (e.g., phosphorous) aids the delocalization of the carbon-centered charge accumulation. While compounds possessing unsupported U=C bonds are difficult to prepare under typical synthetic conditions, Andrews and co-workers have identified alkylidene and even alkylidyne species such as H2C=UX2 and HC≡UX3 (X= H, F, Cl or Br) in low-temperature noble-gas matrix isolation experiments, providing evidence for such bonding motifs[35-37]. However, these compounds are too reactive to be isolated under typical synthetic conditions. The knowledge gained from studying stable compounds featuring U=C multiple-bonds can be extended to uranium carbide ceramics, which offer improved thermal density and higher conductivity over current UO2 nuclear fuels[38-40].

Based on our recent success with the isolation of new actinide endohedral fullerenes, we turned our attention towards the synthesis and isolation of carbon cages encapsulating clusters that possess U=C bonds. The internal hollow cavity of C2-fullerenes can encapsulate and stabilize novel metallic clusters, especially some which are highly reactive and virtually impossible to prepare independently[41-43]. Recently, we reported the first crystallographically characterized examples of actinide endohedral metallofullerenes (Th@C82 and U@C2, 2n = 74, 82) and described their unique electronic properties[43,44]. We therefore hypothesized that fullerene cages would provide an ideal architecture and electronic environment to trap and stabilize unique actinide clusters with novel bonding motifs.

Herein, using the Ih(7)-C80 fullerene cage as a molecular nanocontainer, we report, to the best of our knowledge, the first structurally characterized example of unsupported uranium U=C bonds found in UCU@Ih(7)-C80, possessing the unprecedented dimetallic-carbide cluster U=C=U. X-ray crystallographic analysis reveals two very short U=C bonds of 2.033(5)/2.028(5) Å, with an unexpected nonlinear U=C=U bond angle of 142.8(3)°.

Results

Synthesis of U2C@C80

U2C@C80 was synthesized by the Krätschmer–Huffman arc discharge method[45]. Graphite rods, packed with graphite and finely dispersed U3O8, were vaporized in an arcing chamber under a He atmosphere. The compound was isolated and purified using a multistage high-performance liquid-chromatography protocol (HPLC). The composition and purity of the isolated U2C@C80 was confirmed by high-resolution matrix-assisted laser desorption-ionization time-of-flight positive-ion-mode mass spectrometry (MALDI-TOF/MS), which presents a prominent molecular ion peak with a mass-to-charge ratio of m/z = 1448.103 (Supplementary Fig. 1), corresponding to the [U2C81]+ empirical formula. The mass spectral isotopic distribution pattern matches the theoretically predicted one, thus confirming the molecular composition. In addition, an energy dispersive spectroscopic analysis of the purified sample was employed to determine the elementary composition of the compound. The spectrum shows characteristic peaks of uranium and carbon (Supplementary Fig. 2), which confirmed the assignment of this molecule to U2C81. Moreover, the powder X-ray diffraction pattern of U2C@C80 proves that the purified sample of U2C@C80 used for the experimental characterizations below, exists as a pure phase (Supplementary Fig. 5).

Molecular structure of UCU@Ih(7)-C80·[NiII-OEP]

The molecular structure of UCU@C80 was determined by single-crystal X-ray diffraction analysis. Slow diffusion of nickel(II) octaethylporphyrin (NiII-OEP) in benzene into a CS2 solution of UCU@C80 yielded black UCU@Ih(7)-C80·[NiII-OEP] cocrystals (1). In cocrystal 1, UCU@C80 is observed to adopt a slightly distorted icosahedral Ih(7)-C80 cage structure (Fig. 1a) that is highly ordered. Similarly, inside the fullerene cage, the central carbon atom C0 of the endohedral UC0U unit is fully ordered, while the U atoms are slightly disordered. Two major U positions (U1 and U2) have common dominant occupancy of 0.853(3), while the residual occupancies are located on both sides within ca. ½ to 1 Å distance from the two main positions, possibly indicating some large amplitude motions of the UC0U cluster (Supplementary Fig. 6). The distances between the major U1,2 sites and the carbons of the adjacent aromatic rings of the cage all lie within a narrow range of around 2.50 Å (2.471(5) to 2.543(5) Å, Fig. 1b), corresponding to distances of the ring centers Cti to Ui of 2.042(1) Å and angles Cti−Ui−C0≈159(1)°. The U–C0 bonds are 2.03 Å (U1–C0 = 2.033(5) Å, U2–C0 = 2.028(5) Å), forming a nearly linear Ct1–U1—U2–Ct2 chain (Fig. 1b) bridged by a non-linear –C0– carbide anion.

Fig. 1

ORTEP drawing of UCU@Ih(7)-C80·[NiII-OEP] with 40% probability ellipsoids. a UCU@Ih(7)-C80·[NiII-OEP] structure showing the relationship between the fullerene cage and the [NiII-OEP] ligands. The two U1/U2 sites have common occupancy of 0.853(3). Four minor U sites (Supplementary Fig. 6) and the solvent molecules are omitted here for clarity. b Fragment view showing the interaction of the major U1–C0–U2 cluster with the closest aromatic ring fragments of the cage with centers Ct1 and Ct2. The orange line connects Ct1–U1–U2–Ct2

For hexa-coordinated uranium, Pyykkö found that two short and strong axial and four long and weak equatorial bonds are energetically near-equivalent to two long and weak axial and four short and strong equatorial bonds[46,47]. Therefore, particularly short axial bonds are expected when additional equatorial coordination is absent, as in the present case. Indeed, most U-Ct(arene) distances are reported between 2.5 and 2.9 Å, for instance 2.651(4) to 2.698(4) Å for the dinuclear U(V) inverse sandwich complex [{UV(Ts)}2(η6:η6-C6H5Me)][48]. However, by far the most salient features of 1 are the unusually short U–C0 bond lengths, (U1–C0 = 2.033(5) Å, U2–C0 = 2.028(5) Å) and the unusual bond angle of 142.8(3)° for U1–C0–U2 (Table 1). In fact, the shortest U–C bond distance in poly-coordinated U-complexes to date is 2.184(3) Å, reported by Liddle et al. for the U(VI) carbene complex U=[C(Ph2PNSiMe3)2](O)Cl2[27]. Interestingly, the U–C bond lengths of the UC0U unit are between the 2.067(7) Å and 1.948(7) Å, theoretically calculated distances for uranium alkylidene and alkylidyne complexes U=CH2(H)F and U≡CH(Cl)3[35-37], and are nearly identical to the 2.01 Å distance calculated from the sum of the covalent radii[49]. Furthermore, the U–C0 distances in 1 are longer than the triple bond of the U≡O unit (1.78 Å) in the uranyl dication UO22+, the U≡NR bonds (1.844(6) Å) of U(NtBu)2(I)2(THF)2), and the terminal U≡N bond (1.83(2) Å for the uranium nitride complex (TrenTIPS)UN), even when accounting for differences in the atomic radii of O, N, and C of ca. 0.05 Å each[18,23,50].

Table 1

Calculated and experimentally derived geometric parameters of UCU groups

Geometric parameter	Calc.a UCU@(3I)2 (molecule)	Calc.a UCU@(C7H7)2 (molecule)	Calc.a UCU@Ih(7)-C80 (molecule)	Exptl.b UCU@Ih(7)-C80 •[NiII-OEP] (crystal)
U–C /Å	2.024	2.074	2.022	2.033/2.028
U∙∙∙U /Å	3.918	4.118	3.917	3.849
Ligand→U /Å	2.33c	2.52–2.55	2.47–2.51	2.47–2.54
U–C–U / o	150.9	166.1	151.3	142.8
Fullerene C–C /Å	—	—	100 × 1.43 ± 0.01	110 × 1.43 ± 0.03
‘π-Donating’ C–C /Å	—	14 × 1.43 ± 0.01	20 × 1.46 ± 0.01	10 × 1.48 ± 0.01

aFrom scalar relativistic ZORA Kohn–Sham PBE VTZp approaches

bFrom crystal structure analysis by X-ray diffraction

cFor comparison of I−U with C−U, we have subtracted the I–C bond-radii difference of 0.55 Å from the actual I–U distance

Altogether, the bond lengths observed for 1 strongly support an Arene≡U=C=U≡Arene structure possessing two uranium–carbon double bonds formed by a bridging carbide atom, formally C4−. The bending of the U=C=U cluster of 1 was an unexpected exception to the common linear geometries of sp1–pπ2 hybridized carbon atoms found for main-group molecules of type E=C=E (E=O, NR, CR2), though a few cases with angles between 180° and 90° have been theoretically proposed and observed in recent years[51].

Computational studies of molecular structure and bonding in UCU@Ih(7)-C80

Quantum-chemical calculations were performed to further verify the unique structural parameters of UCU@Ih(7)-C80 and gain insight about the bonding of the encapsulated bent U=C=U cluster. The molecular structure was modeled for free UCU@Ih(7)-C80 and for LUCUL (L = C7H7, 3I) model molecules, using quasi-relativistic density-functional (DF) approximations (Supplementary Figs. 13, 14) and ab-initio CASSCF(10e,12o)-PT2 approaches (Table 1). Additional information of the structure, orbitals and electronic states of LUCUL is available in Supplementary Figs. 15–19 and Supplementary Tables 3–5. For the three calculated cases, the formal oxidation state[52] is identified as 5+ for each U, derived from the formal [L3−≡>U5+<=C4−=>U5+<≡L3−] Lewis structure, where L3−≡>stands for the cage arene units that accept 3 electrons from each uranium atom and forms 3 dative pair bonds with the nearly empty U-5f6d valence shell. In the case of UCU@Ih(7)-C80, the two L3− units are the ligating arene units of the encapsulating C806– closed-shell cage. The C and I atoms of the model units (C7H7) and (3I) have similar electronegativity as C of (C80), and the model units also become closed-shell ligands upon accepting 3 electrons from each U. The effective physical partial charges in all three cases remain small, around −1 for formal C4−, around +½ to +1 for formal U5+, and small for the ligand atoms of C80, 2C7H7 or 6I (Supplementary Table 1).

The atomic free valence and the spin and orbital populations (Supplementary Table 1) indicate a single f-electron on each U atom. The experimental U–C0 bond lengths, though unprecedentedly short (2.03 Å), are well reproduced in all three computed models (2.02 Å for ligands C80 or (I3)2, 2.07 Å for (C7H7)2). Further, the observed unusual U–C0–U bond angle of 142.8° is comparable to those computed for UCU@C80 (151.3°) and (3I)UCU(3I) (150.9°), while the bulkier ligands of (C7H7)UCU(C7H7) led to 166.1°, which is still far from linear (Table 1). Neither the bond lengths nor the bond angle of the UCU unit appear to be significantly distorted by the encapsulation in C80.

We chose the localized equivalent molecular orbital picture of the two simpler model systems to sketch the correlated and spin–orbit coupled valence electronic structure of ligated UC0U. (C7H7)UCU(C7H7) (Fig. 2) and (3I)UCU(3I) (Supplementary Fig. 15) gave very similar results. Each U forms a strong σ-bond to the central C0 atom through covalent overlap of hybridized U(5 f6d) orbitals with C(2s2p) hybrids. The U–C bonds of dominant π-character can be described as three-center two-electron (3c2e) bonds that possess some σ-admixture. The large electronegativity difference of C and U (1.3 Pauling units) leads to bond-pair polarization toward C0, especially for the π-type pairs, which are strongly localized on C0 and lead to its effective negative charge. The U–C interactions are characterized by Mayer bond orders (BO) of 1.4 from DF calculations (comparable to 1.6 for the U=As double bond in [K(B15C5)2] {[N(CH2CH2NSiPr3)3] U=AsH})[26] and by a Roos effective bond order (EBO) of 1.9 (counting the formally bonding—antibonding NOs from our CASSCF calculation of (C7H7)UCU(C7H7), which substantiates the U=C double bond character).

Fig. 2

Natural Localized Molecular Orbitals (NLMOs) of sandwiched UCU. The UCU@(C7H7)2 molecule formally consists of a (UCU)6+ unit between two aromatic (C7H7)3− rings. Spin–orbit coupled multi-reference all-electron relativistic wavefunction calculations (state-averaged CASSCF) were analyzed with the NBO program [S27]. The contour values of the NLMOs generated thereby are ±0.042√(e/Å3). a The two 1-center U-5f1-type NLMOs, each occupied by one electron. b The U–C (left) and C–U(right) 2-center 2-electron (2c2e) σ-pair NLMOs, each consisting of 67% C-2sp1, 31% U-5f½ 6d and 2% tail-contribution from the other U. c The vertical π-type and the in-plane dominantly π-type U–C–U 3-center 2-electron (3c2e) NLMOs of 65% C-2ps and 17% of each of the two U-5f½ 6d

The bonding motif is nicely illustrated by the electron localization function (ELF) map in Fig. 3, which shows the electron–density accumulation of the two two-center U–C σ-bonds and of the in-plane three-center U–C–U π-type bond, as well as the multiple maxima of the inner-valence non-bonding 5f electrons. Thus, the bending of the U=C=U unit can be interpreted in terms of the localized negative charge build-up on C0, which generates lone electron-pair density and gives rise to an sp1/sp2-hybridized type geometry at carbon with bond angles in the middle between 120° and 180°. This charge accumulation is partially offset by donation into the U(5f6d) orbitals. Thus the bonding model can be described as highly polarized U=C interactions with partial π-overlap strengthened by electrostatic attractive forces between the cationic U5+ atoms and the anionic C04− carbide bridge.

Fig. 3

ELF-plot of ligated UCU in the molecular plane. The left and right ligands are two (I3) groups. ELF is the ‘electron localization function’ of Becke and Edgecombe, showing electronic accumulations in the inner atomic shells, in the atomic lone-pair regions, and in covalent and dative bond regions

Spectroscopic properties of UCU@Ih(7)-C80

Nuclear Magnetic Resonance (NMR) spectroscopy was used to characterize UCU@Ih(7)-C80. The 13C NMR spectrum of UCU@Ih(7)–C80 measured at 298 K shows only two sharp signals at 138.53 and 125.14 ppm with a 3:1 intensity ratio, corresponding to the sets of 60 and 20 equivalent carbon atoms of the unperturbed Ih–C80 fullerene cage (Supplementary Fig. 7). The apparent symmetry is likely due to fast and large-amplitude librations of the cage, similar to those proposed for M2@Ih-C80 (M = La, Ce) and M3N@C80 (M = Sc, Y)[42,53,54]. Unfortunately, no signal was detected for the single bridging C0 carbon atom, likely due to limited sample amount, as well as to paramagnetic broadening effects from the unpaired 5f1 electrons on each uranium. Noteworthy is that the 13C NMR spectrum measured at lower temperature (i.e., 283 K) shows slightly shifted signals at 138.35 and 124.65 ppm, respectively, with a slightly larger chemical shift difference (Δδ = 13.70 ppm) than those observed at 298 K (Δδ = 13.39 ppm). Such a trend is very similar to that observed for Ce2@Ih–C80[55], consistent with the paramagnetic nature of UCU (see below) and with the presence of an unpaired f electron on the formal U5+ (5f16d07s0) ions (Supplementary Table 1)[54]. On the other hand, the slight temperature dependence of Δδ/ΔT = 0.02 ppm/K might also be accounted for by the libration of the UCU cluster that neutralizes the paramagnetic effects (mainly pseudocontact interactions) for all carbon atoms.

We also utilized Fourier transform infrared absorption (FTIR), Raman emission and photoluminescence (PL) spectroscopies together with quantum-chemical calculations to further characterize UCU@Ih(7)-C80. An experimental and theoretical spectral overview is presented in Fig. 4. The high-wavenumber range from 1600 to 1100 cm‒1 contains 104 carbon cage C–C stretching vibrations, resembling those of other Ih-C80 based endohedral fullerenes such as Ln3N@C80[41,56]. Typical characteristics are the major overlaid bands around 1380 cm−1 and the featureless gap between 1100 and 900 cm‒1. The remaining 130 collective deformations originating from bending and torsional motions of the carbon cage atoms show up between 900 and 200 cm−1.

Fig. 4

Observed and theoretically predicted vibrational features of UCU@Ih(7)-C80. The left upper curve (black) presents the observed Infrared absorption (IR) spectrum vs. wavenumber from 1600 to 600 cm‒1, with quantum-chemical density-functional simulation below (in blue). On the right is the observed Raman emission (in red) vs. wavenumber from 600 to 100 cm‒1, with calculated wavenumbers (‘Theor. Freqs.’) in the 600–0 cm‒1 range (below in magenta). Observed local UCU vibrations above 100 cm‒1 are indicated by heavy dots. Observed and calculated spectral C80-cage features represent either a single mode, or an overlay of several near-degenerate ones. indicates a contamination by CS2

Quantum-chemical calculations suggest that an endohedral cluster of n atoms (here n = 3) shows 3 frustrated translational rocking modes against the cage (at 33, 59, and 115 cm‒1), 3 frustrated torsional wagging modes (at 16, 50, and 142 cm‒1), and 3n – 6 = 3 internal vibrations. The latter comprise the UC0U bending mode (at 97 cm‒1), the symmetric stretch (coupled to near-degenerate cage deformations from 250 to 300 cm‒1), and the asymmetric stretch (at 780 cm‒1). The observed weak features in the Raman spectrum above 100 cm‒1 at 126, 148, and 277 cm‒1 agree reasonably well with the predicted UCU modes given above (in bold italics, estimated as weakly Raman active). Notably, the asymmetric UCU stretch appears as a pronounced feature in the IR spectrum at 785 cm‒1, which is consistent with the high negative charge on the C0 atom of the endohedral UCU cluster, and with the crystallographic observation of very short axial U=C bonds.

Further, an experimental violet-blue photo-luminescence progression of U2C@Ih(7)-C80 starting at a wavelength of 430 nm (~2.95 eV) in steps of ca. 1425 cm‒1 (Supplementary Fig. 8) may originate from a ligand-to-metal charge-transfer excitation of the highest occupied molecular orbital, which has dominant density on the fullerene cage near the U atoms. There are also Cring-breathing vibrations around 1400 cm−1 localized in the spatial Cring–U–C–U–Cring region.

The redox properties of UCU@Ih(7)-C80, as investigated by means of cyclic voltammetry, show a surprisingly small electrochemical gap of only 0.83 eV (Table 2, Supplementary Fig. 9). Most reported transition metal cluster fullerenes of Ih(7)-C80 have gaps larger than 1 eV[41]. Typical M3N@Ih(7)-C80 fullerenes have gaps of 1.8–2.2 V. In particular, the first reduction potential of UCU@Ih(7)-C80 (−0.41 V) is much more positive than that of any other reported M3N@Ih(7)-C80 fullerene[57], e.g., −1.26 V for Sc3N@Ih(7)-C80[58] and −0.94 V for TiLu2C@Ih(7)-C80[59]. UCU@Ih(7)-C80 obviously exhibits a much better electron-accepting ability than other I(7)-C80 clusterfullerenes. This significant difference between UCU@Ih(7)-C80 and other reported M3N@C80 fullerenes indicates a major difference of their electronic structures, which may be caused by the encapsulation of the U=C0=U cluster with an open d–f shell electronically interacting with the cage. Interestingly, when compared to most other Ih(7)-C80 cluster fullerenes, the redox potentials of UCU@Ih(7)-C80 are close to those of dilanthanide-fullerene, such as La2@Ih(7)-C80, which exhibits an electrochemical gap of 0.87 eV[60]. In general, the electrochemical gaps are related to the HOMO-LUMO gaps of the closed-shell fullerene subsystems. However, in the present case, the encapsulated U=C0=U cluster has a low-lying open U-5f shell, causing an abnormally small f–f gap of the UCU@Ih(7)-C80 molecule. The reversible reduction of UCU@Ih(7)-C80 is also quite remarkable, since M3N@Ih(7)-C80 fullerenes with group 3 metals exhibit irreversible reductions[57].

Table 2

Redox potentials and electrochemical gaps of UCU@Ih(7)–C80, La2@Ih(7)–C80, and Sc3N@Ih(7)–C80

Fullerenes	E 2+/+ a	E +/0 a	E 0/− a	E -/2− a	E 2−/3− a	E gap b	Ref.
UCU@Ih(7)-C80	+1.05d	+0.42c	−0.41c	−1.34c		0.83	This work
La2@Ih(7)-C80	+0.95c	+0.56c	−0.31c	−1.72c	−2.13d	0.87	[60]
Sc3N@Ih(7)-C80		+0.97c	−1.26d	−1.62d	−1.82d	2.23	[58]

aRedox potentials in V vs. ferrocene couple

bElectrochemical gaps in eV

cHalf-wave potential (reversible redox process)

dPeak potential (irreversible redox process)

The XPS spectrum of UCU@Ih(7)-C80 was also recorded (Supplementary Fig. 10). The U-4f7/2 ionization appears at 378.6 eV. In general, ionization energies increase with the formal oxidation state, e.g., from the U-4f7/2 value of ca. 3771/3 eV for the neutral uranium metal to ca. 3813/4 eV for some uranyl(VI) salts. The present value of 378.6 eV is consistent with the very low effective charges on U as computed for UC0U@C80 (see below, U ~ +1.0e; Co ~ −1.0e, Cfullerene ~ ± 0.0e), caused by the strong electron pair donation of formal C806− and C04− to the formal U5+ ions.

To gain further insight into the unique electronic characteristics of the U=C=U cluster in UCU@C80, the magnetic susceptibility was measured (see Supplementary Fig.11 and Methods). The magnetic curves exhibit two regimes, (i) one at “high temperatures” from 300 K down to ca. 60 K, and (ii) another one at “low temperatures” from ca. 40 K down to 2 K. The “high temperature” region exhibits a huge, basically temperature-independent paramagnetism (TIP) of yet unknown origin. It is unknown for uranium complexes with common ligands[21,22,48,61-65] and may possibly be due to the cooperative coupling of UCU and C80 in the solid and related to the low band-gap calculated for UCU@C80 (see Supplementary “Quantum Computational Methods”).

At “low temperatures”, we observed a Curie type paramagnetism ~μ2/T on top of the TIP, with net μ-values increasing with temperature from below to above a μBohr, which appears compatible with hardly coupled unique units of type η3L≡U(V)=C=U(V)≡Lη3, i.e., the endohedral di-uranium carbide fullerene molecules. The “low-temperature” literature values of μeff[21,22,48,61-65] for mono- and bridged di-metallic U(V) and U(IV) complexes, typically increase, too, from low values up to 3 μBohr. It shows the expected temperature dependent weak coupling of two spins and with the two partially quenched U-f±3 orbital angular momenta, the experimental curves being reasonably reproduced by numerical simulations (see Supplementary Method “Quasi-relativistic correlated ab initio approaches” on the (C7H7)UCU(C7H7) model).

Attempts to further resolve the electronic structure using EPR spectroscopy were unsuccessful, as no clearly defined signal was observed at 4 K (Supplementary Fig. 12). Factors such as near degenerate electronic states and cage shielding effects as well as the huge TIP complicate the analysis. Thus, additional studies of the electronic properties of this unique and unprecedented system are warranted and will be communicated in due time.

Discussion

In summary, the overall agreement between the crystallographic, luminescence, Raman, IR, core-electronic, magnetic and voltammetric results, and the quantum-computational findings, conclusively show that a pure UCU@Ih(7)-C80 compound was synthesized, which consists of a weakly perturbed Ih-C80 cage and a sandwiched, bent and strongly bound polar endohedral near-symmetric U=C=U unit. The quantum-chemical results suggest that both U atoms in the UCU unit have a formal +5 oxidation state. The encapsulation causes little bond length or angle deformation of the UCU fragment. The encapsulation protects the axially ligated U atoms from binding with harder or more electronegative ligands like O or N (as compared to carbon in the forms of C04− and C806−) so that C0 can participate in primarily axial L≡>U=C0 bonding at unusually short U–C0 distances of 2.03 Å, and is not pushed into the weaker bonding equatorial plane of an E=U=E building block. The discovery of unsupported U=C bonding in a molecular compound confirms the distinction between “diagonal” and “poly-coordinated” uranium, and between “axial” covalent and “equatorial” dative bonding mechanisms. The work reported here offers a deeper understanding of the fundamentals of uranium bonding properties. This study also demonstrates that fullerene cages can be utilized as effective nanocontainers to stabilize and study rare and reactive clusters which contain actinide metal–ligand bonds.

Methods

Synthesis, separation and purification of U2C@Ih(7)-C80

A soot containing uranium fullerenes was synthesized by a direct-current arc discharge method. Graphite rods of U3O8 and graphite powder were annealed in a tube furnace at 1000 °C for 20 h under an Ar atmosphere and then burned in the arcing chamber under a 300 Torr He atmosphere. In total 1.87 g of graphite powder and 1.83 g of U3O8 (24:1 atom number ratio) were packed in each rod (6.7 g, without filling). The collected raw soot was refluxed in chlorobenzene under Ar atmosphere for 12 h and the solution was filtered. After solvent removal, the fullerenes were extracted by dissolving in toluene. On average ca. 140 mg of crude fullerene mixture per rod was obtained. In total, 60 carbon rods were burned in this work. Separation and purification of U2C@Ih(7)-C80 was achieved by a three-stage HPLC procedure (Supplementary Fig. 3, using columns including a Buckyprep M column (25 × 250 mm2, Cosmosil, Nacalai Tesque, Japan), a Buckprep D column (10 × 250 mm2, Cosmosil, Nacalai Tesque, Japan) and a Buckprep column (10 × 250 mm2, Cosmosil, Nacalai Tesque, Japan). Toluene was used as the mobile phase and the UV detector was adjusted to 310 nm for fullerene detection. The HPLC procedures and corresponding MALDI-TOF spectra for the isolated fractions are shown in Supplementary Figs. 3, 4. In total, ca. 2 mg of highly purified U2C@C80 was obtained for characterization.

Single crystal and powder X-ray diffraction analysis

Black co-crystals of UCU@Ih-C80·[NiII-OEP] were obtained by allowing the benzene solution of [NiII-OEP](1.2 mg/mL, 0.8 mL) and the CS2 solution of UCU@Ih-C80(1 mg/mL, 0.5 mL) to diffuse together. X-ray diffraction data were collected at 120 K using a diffractometer (APEX II; Bruker Analytik GmbH) equipped with a CCD collector. The Multiscan method was used for the absorption correction. The structure was resolved using direct methods (SIR2004) and refined on F with the full-matrix least-squares approach using SHELXL2014[66,67] within the WinGX package[68].

Crystal Data for UCU@Ih-C80·[NiIIOEP]·1.5C6H6·CS2: C127H53N4NiS2U2, M = 2159.70, 0.05 × 0.05 × 0.03 mm3, monoclinic, P 21/c (No. 14), a = 17.678(4) Å, b = 16.970(3) Å, c = 26.695(5) Å, β = 106.65(3)°, V = 7673(3) Å3, Z = 4, ρ = 1.934 g cm−3, μ = 4.578 mm−1, θ = 2.688–27.092°, T = 120 K, R1 = 0.0610, wR2 = 0.1045 for all data; R1 = 0.0388, wR2 = 0.0610 for 16,771 reflections (I > 2.0σ(I)) with 1278 parameters and 784 restraints. Goodness of fit indicator 1.007. Maximum residual electron density 2.096 e Å−3. For further details see the Supplementary Data 1, 2.

In addition, a powder X-ray diffraction analysis of UCU@C80 was conducted at the Shanghai Synchrotron Radiation Facility (SSRF) with wavelength λ = 0.6199 Å, and the intensity vs. 2θ pattern was recorded, see Supplementary Fig. 5.

Computational methods

Open-shell Kohn–Sham single-reference calculations were carried out with the ADF-2016, Gaussian-09 and ORCA-4.0 software, applying scalar and spin–orbit-coupled quasi-relativistic approaches (RECP or ZORA) with triple-ζ double-polarized basis sets. Bonding analyses were then performed using the NBO5.0 and Multiwfn codes. Multi-reference CASSCF, CASPT2, and RAS-SI/SO wave-function calculations were performed using the MOLCAS-V8.1 software applying the quasi-relativistic DKH Hamiltonian. Further details are given in the Supplementary “Quantum Computational Methods”.

Data availability

The coordinates for the X-ray structure of 1 are available free of charge from the Cambridge Crystallographic Data Centre under deposition no. CCDC 1572891. All other data supporting the findings of this study are available from the corresponding authors on request.