Structure analysis. Single-crystal X-ray diffraction (SCXRD) demonstrates that the structure of Cl12@Gd158Co38 is similatr to Br12@Gd158Co38, here only Cl12@Gd158Co38 as example to be discussed in detail. Cl12@Gd158Co38 crystallizes in the trigonal crystal system, R-3 space group. The cationic core of Cl12@Gd158Co38 constitutes of six Cl12@Gd27Co7(CO3)15(OAc)4(µ3-OH)40(CH3NH2CH2COO)(MIDA)7(H2O)14 (Cl12@Gd27Co7, in which Gd4, Co4, and O29 are distributed in the C3 axis) building units (Fig. 1a, 1b). Cl12@Gd27Co7 is made up of one highly symmetrical cationic unit Gd26Co5(CO3)15(OAc)3(µ3-OH)40(MIDA)7(H2O)14 (Gd26Co5) and two Cl− anionic templates, one Gd, two Co ions, one OAc−, and one CH3NH2CH2COO− ligand (decomposition of H2MIDA, Fig. S14, Scheme S2). In addition, Gd26Co5 (Fig. 1c) can be regarded as three different motifs: type I, formulated as Gd5Co2(MIDA)2(OAc)(CO3)3(µ3-OH)4 (Gd5Co2), is distributed in the top of Gd26Co5; type II, formulated as Gd16(MIDA)2(CO3)9(µ3-OH)22 (Gd16), is distributed in the middle of Gd26Co5; type III, formulated as Gd5Co3(MIDA)3(OAc)(CO3)2(µ3-OH)4 (Gd5Co3), is distributed in the bottom of Gd26Co5. Besides, two Gd7(MIDA)(CO3)3(µ3-OH)11 (Gd7, Fig. 1e) is connected by two Gd(CO3) to form Gd16. Based on four CO32− and four µ3-OH− groups, Gd5Co2, Gd16, and Gd5Co3 motifs are joined together to form Gd26Co5. Cl12@Gd158Co38 has a 3-fold symmetric (C3) axis and an inversion center (Fig. S11-12). Three asymmetric building units Cl12@Gd27Co7 obtained by rotation are connected together alternately via CO32− and µ3-OH− anions to form a trimer Gd79Co19(CO3)45(OAc)12(µ3-OH)120(CH3NH2CH2COO)3(MIDA)21 (Gd79Co19). Two trimers Gd79Co19 obtained by inversion are further joined into Gd158Co38.
It is worth noting that a large series of CO32− anions as important templates and linkers among metal ions, deriving from the decomposition of organic ligands, exhibit a rich variety of coordination modes (Fig. S3), reflecting the complexity of the nanocluster structure, and the adaptability of the anion templates. Meanwhile, the main ligand MIDA2− also shows unusual and diverse coordination modes (Fig. S4). For example, high chemical affinity of N atom from organic ligands tends to coordinate with 3d ions38, but N atoms in this work are also linked with 4f ions.
Interestingly, in Cl2@Gd27Co7, two crystallographically independent halide ions (Cl1 and Cl2) are found in the center of similar two Gd-CO3-OH cages [(Gd16(CO3)6(µ3-OH)7, Fig. 2a, 2c], hydrogen-bonded to 7 µ3-OH− groups, respectively (Fig. 2d, 2f, distance of Cl···O: from 3.219 Å to 3.419 Å, angle of Cl-H-O: from 155.15° to 176.16°)22. Additionally, six Cl1 and six Cl2 ions obtained by rotainversion form one icosahedron [Cl12]12− cage (Cl12, distance of Cl···Cl: from 7.216 Å to 8.445 Å, Fig. 2e, 3b). The 12 Gd-CO3-OH cages templated by Cl− ions form the main structure [Gd116(CO3)66(µ3-OH)84] (Fig. 2b) of Ln-core [Gd158(CO3)66(µ3-OH)236] by sharing 4f metals or anions (CO32− and OH−). Although Cl− ions as the templates in metal clusters have been obtained (such as {Gd36Ni12}22 templated by 2 Cl−; {Gd8Cr4}47 templated by one Cl− and one ClO4−), this 3d-4f nanocluster with more than ten Cl− templates (Cl12 with one icosahedron pattern) is firstly reported, and it has important guiding significance for the prediction and construction of high-nuclearity 4f-containing nanoclusters.
For the sake to facilitate the illustration and comprehending of the complex metal skeleton, Cl12@Gd158Co38 can be disassembled into three components (Fig. S15): (i) six Gd(MIDA)2(CO3)(µ3-OH)2 groups, exhibiting a near-octahedral geometry; (ii) six Co(MIDA)3 groups, showing a hexagonal arrangement; (iii) one Cl12@Gd152Co32(CO3)24(µ3-OH)116 group (Cl12@Gd152Co32), revealing a typical core-shell structure and modular features. From the inside to outside, Cl12@Gd152Co32 (Fig. 3 and S16) is treated as Gd20@Cl12@Gd48Co32@(Gd)12@(Gd12)6, and presents one dodecahedron Ln-core of 20 Gd3+ ions [Gd20(CO3)12] (Gd20, Fig. 3a), one icosahedron of 12 Cl− ions (Cl12, Fig. 3b), one truncated-cube-like cage of Gd132Co32 (Fig. 3i). Gd132Co32 comprised of eight propeller-like heterometallic building blocks [Gd6Co4(µ3-OH)9] (Gd6Co4, Fig. 3e) as truncated-cubic vertices forming main metal framework [Gd48Co32(µ3-OH)44] (Gd48Co32, Fig. 3d), 12 Gd(µ3-OH)3 groups (Gd, Fig. 3f) as edges, six saddle-shaped motifs ([Gd12(CO3)2(µ3-OH)6], Gd12, Fig. 3g) as faces. So, Gd132Co32 can be viewed as (Gd6Co4)8@(Gd)12@(Gd12)6. Here are 10 Co3+ and 28 Co2+ ions in final products, determined by XPS (Fig. S33) and charge balance.
Althrough the metal structures of Br12@Gd158Co38 and Cl12@Gd158Co38 are very similar, ther are some difference between two compound: (i) six CH3O ligands and 30 anions (12 Cl− and 18 Br−) are exisisted in Br12@Gd158Co38 to balace the charge balance; (ii) but 36 anions (30 Cl−) are exisisted in Cl12@Gd158Co38.
Magnetic properties. The large presence of metal ions inspires us to investigate the magnetic properties of Br12@Gd158Co38 and Cl12@Gd158Co38. The plot of temperature-dependent magnetic susceptibility (χMT-T) was studied under 1.0 kOe direct current (dc) field with the scope of 1.8–300 K (Fig. S38, S39). The χMT values under the room temperature were 1380.75 (Cl12@Gd158Co38) and 1384.91 (Br12@Gd158Co38) cm3 K mol− 1, which is bigger than the theoretical value of 1296.75 cm3 K mol− 1 for 158 uncorrelated Gd3+ (S = 7 /2, g = 2) and 28 high-spin Co2+ (S = 3/2, g = 2)45. The difference between the theoretical and test values is put down to the remarkable orbital contributions of the high-spin Co2+ ions5. As the temperature goes down, the value of χMT gradually decreases and achieves 726.25 (Cl12@Gd158Co38) and 696.30 (Br12@Gd158Co38) cm3 K mol− 1 at 1.8 K. This behavior is mainly ascribed to the depopulation of Kramers excited state levels of octahedral coordination environment for Co2+ ion because the magnetic interaction between 3d and 4f or 4f and 4f ions. Based on the Curie-Weiss Law, fitting the plot of χM−1 versus T shows parameters, [C = 1415 cm3 K mol− 1 and θ = −6.04 K (Cl12@Gd158Co38); C = 1462.60 cm3 K mol− 1 and θ = −3.20 K (Br12@Gd158Co38)], for the sum contribution of orbital of Co2+ ion and the coupling between metal ions (Fig. S40, S41). The field-dependent magnetization (M-H) of Cl12@Gd158Co38 was performed in the range of 1.8–20 K at 0–7 T (Fig. S42, S43). The curves of M vs H represent a steady rise in magnetization and attain 1045.29 NµB (Cl12@Gd158Co38) and 1058.57 NµB (Br12@Gd158Co38) NµB at 7 T at 1.8 K, which are slightly lesser than the expected value 1190 NµB for 158 uncorrelated Gd3+ (S = 7/2, g = 2) and 28 high-spin Co2+ (S = 3/2, g = 2)45. It also results from the orbital effect of Co2+, which gives an effective spin Seff = 1/2 and the magnetization is usually ~ 2.1 NµB per Co2+ ion48.A large class of Gd3+ ions in Cl12@Gd158Co38 and Br12@Gd158Co38 urge us to explore its magnetocaloric effect. The calculated magnetic entropy changes ∆Sm were evaluated via using the Maxwell relation (∆Sm(T) = ∫[∂M(T, H)/∂T ]H dH)29. As shown in Fig. 4, S44, the values of -∆Smmax are 46.95 (Cl12@Gd158Co38) and 46.30 (Br12@Gd158Co38) J kg− 1 K− 1 at 2.0 K at 7 T. The values are smaller than their theoretical values (64.33 and 62.32 J kg− 1 K− 1) by applying the formula of -∆Sm = nRln(2S + 1), which are attributable to the presence of possible antiferromagnetic interaction30. Nevertheless, their values were much larger than known 3d-4f cluster complexes (Table S6). And the value of Cl12@Gd158Co38 is the largest at present. In low magnetic field, Br12@Gd158Co38 and Cl12@Gd158Co38 also show prominent magnetocaloric effect with -∆Sm = 20.13, 20.81 J kg− 1 K− 1 at 2.0 K and 2 T, respectively (Table S7)49.