Temperature dependence of spherical electron transfer in a nanosized [Fe14] complex.

The study of transition metal clusters exhibiting fast electron hopping or delocalization remains challenging, because intermetallic communications mediated through bridging ligands are normally weak. Herein, we report the synthesis of a nanosized complex, [Fe(Tp)(CN)3]8[Fe(H2O)(DMSO)]6 (abbreviated as [Fe14], Tp-, hydrotris(pyrazolyl)borate; DMSO, dimethyl sulfoxide), which has a fluctuating valence due to two mobile d-electrons in its atomic layer shell. The rate of electron transfer of [Fe14] complex demonstrates the Arrhenius-type temperature dependence in the nanosized spheric surface, wherein high-spin centers are ferromagnetically coupled, producing an S = 14 ground state. The electron-hopping rate at room temperature is faster than the time scale of Mössbauer measurements (<~10-8 s). Partial reduction of N-terminal high spin FeIII sites and electron mediation ability of CN ligands lead to the observation of both an extensive electron transfer and magnetic coupling properties in a precisely atomic layered shell structure of a nanosized [Fe14] complex.

Introduction

To date, numerous nanosized transition metal complexes have been synthesized[1], which can be predetermined to have various fascinating structures via molecular self-assembly process[2,3]. One of the challenges in realizing their functions is the introduction of fast electron transfer over a long distance that ubiquitously occurs in chemical and biological systems[4,5]. The low-nuclearity compounds have been intensively investigated as intramolecular electron transfer systems[6-13]. The recent interesting work on the [Fe2]V mixed-valence compound exhibited changes in the electron transfer rate observable by Mössbauer spectroscopy with a ferromagnetic S = 9/2 ground state[14]. It is still challenging to isolate nanoarchitectures where electron transfer occurs over more metal centers, toward revealing extraordinary spectral and electronic behavior[15]. In the presence of intermetallic coupling in high-nuclearity complexes, the number of hopping d electrons and their interactions would be expected to increase. This would be accompanied by some intriguing processes, such as electron transfer along topology-specific pathways, many-electron transfer processes, and externally induced charge separation in a confined nanospace[16-18]. However, intermetallic electron transfer-mediated through bridging ligands are normally weak; thus, properties related to confined electron transfer processes in discrete nanosized complexes have remained hypothetical thus far[19,20].

Prussian blue analogs have received considerable attention as they can demonstrate the coexistence of electron transfer and exchange interaction through a cyanide bridge[21,22]. The realization of partially reducing or oxidizing the interacting metal centers in a nanoarchitecture is proposed as a route to achieve the extensive multielectron transfer with novel magnetic and electronic properties[23,24]. Here, we report a cyanide-bridged complex [FeII(Tp)(CN)3]8[Fe2.667+(H2O)(DMSO)]6 (abbreviated as [Fe14], Tp− = hydrotris(pyrazolyl)borate; DMSO = dimethyl sulfoxide), in which two extra 3d electrons hop at a rate faster than the time scale of Mössbauer measurements (<~10−8 s) at room temperature. The key feature of this complex is that only two of the six N-terminal Fe-hs sites in the FeII–ls–CN–Fehs structure are successfully reduced (hs = high spin and ls = low spin), thereby reducing the potential barrier for two-electron hopping at the six N-terminal Fe-hs sites. Furthermore, the rates of the intramolecular electron transfer exhibit a distinct temperature dependence that can be described with an Arrhenius law in the nanosized complex. Another important characteristic is that a ferromagnetic interaction operates in [Fe14] with a ground spin S = 14, which implies that extra 3d electrons hop around the exchange-coupled atomic layer shell.

Results

Crystal structure

The synthesis of the [Fe14] complex proceeds via the reaction of Fe(BF4)2·6H2O with Bu4N[Fe(Tp)(CN)3] in a mixture of water and DMSO. The crystalline sample of [Fe14] is thermally stable at room temperature (Supplementary Figs. 1 and 2). The crystal structure of [Fe14] features a tetradecanuclear iron molecular cluster, [Fe(Tp)(CN)3]8[Fe(H2O)(DMSO)]6, as shown in Fig. 1a, b (Supplementary Tables 1 and 2), wherein Fe site in [Fe(Tp)(CN)3] and [Fe(H2O)(DMSO)] is hereafter referred to as the A and B site, respectively (Fig. 1c, d). Each [Fe(Tp)(CN)3] unit at A site is connected to three B-site Fe ions via a facial triad of cyanide-bridged units, where the B-site Fe center is in a weakly distorted octahedron coordinated by axial water and DMSO molecules. The equatorial coordination of B-site Fe is completed by the four N-terminals of cyanide ligands with an average Fe…N bond distance of 2.068(6) Å. Thus, eight [Fe(Tp)(CN)3] building units are symmetrically located at the vertices (A site) of a cubic cage, and six [Fe(DMSO)(H2O)] units occupy six centroid areas (B site), resulting in the face-capped cubic structure. When we consider only the mode of the Fe–C≡N–Fe linkage in the framework, the simplified [Fe14] cage comprises 62 atoms (Fe14C24N24) with a shell having an atomic layer thickness. The diameter is ca. 8.8 Å, which is slightly larger than that of C60 (ca. 7.1 Å)[25].

Fig. 1

X-ray crystal structure determination. a Complete structure of the [Fe14] spherical complex in a ball-and-stick style. The valence state on different Fe sites is indicated by two-color systems, that is, A-site Fe (blue) and B-site Fe (red), and O (pink), S (yellow), N (purple), C (gray), and H (white). b View of the cubic packing of the simplified [Fe14] complex along the crystallographic c-axis. The inset denotes the scanning electron microscope (SEM) image of a single crystal of [Fe14], illustrating the cubic faces. c The simplified [Fe14] atomic layer framework with a diameter of 8.8 Å. d B-site-selective isotope substitution in the molecule.

Electronic behavior

To confirm the valences and spin states of the Fe ions in the complex, 57Fe Mössbauer spectroscopic measurements were performed at various temperatures under zero magnetic field (Fig. 2a). A 57Fe isotopically enriched sample, [Fe857Fe6], was also prepared, consisting of isotopically enriched B sites (57Fe, 96%) and naturally abundant A sites (57Fe, 2%), to mainly produce Mössbauer signals of almost only the B sites. (Fig. 2b) Compared with the enriched sample, the natural [Fe14] sample exhibits one additional quadrupole doublet in the spectrum at both high and low temperatures. Considering its isomer shift (IS) of 0.07–0.15 mm s−1 and quadrupole splitting (QS) of ~0.45 mm s−1, the additional doublet can be assigned to FeII-ls in the [Fe(Tp)(CN)3] unit[21]. The result indicates that the reactant, [FeIII(Tp)(CN)3]−, is completely reduced to [FeII(Tp)(CN)3]2− during the formation of [Fe14]. Furthermore, a comparison of the Mössbauer spectra of [Fe14] and [Fe857Fe6] indicates that the role of the FeII ions at the A sites remains static (FeII-ls) across the entire temperature range. (Supplementary Tables 3 and 4) In the spectra of the site-selective 57Fe isotopically enriched sample, [Fe857Fe6], exclusive information can be observed for the B sites, indicating that there is no electron transfer between the A and B sites in [Fe14]. The spectra recorded below 144 K can be split into two doublets with QS values of 2.20 and 1.18 mm s−1 (at 10.5 K) with a 1:2 ratio of their fractions, assignable to FeII-hs and FeIII-hs, respectively. This finding demonstrates the apparent mixed-valence nature of the B site with the expected fraction of 2:4 for FeII vs. FeIII. Above 164 K, the B-site signals show a pronounced temperature dependence, especially at around 200 K. The signals merge gradually into a unique doublet upon an increase in the temperature, indicating that valence fluctuation occurs due to electron hopping at the B sites. When two extra electrons are shared with the six B-site Fe ions, the averaged B-site Fe valence is expected to be 2.667 for [FeII8]A[FeIII6 + 2e]B. To evaluate such a fluctuation, an electron-hopping relaxation model was applied to the merged spectra between 164 and 297 K for [Fe857Fe6][26]. As a result, the well-fitted profiles were produced with the FeII vs. FeIII area ratio fixed at 2:4. A few additional parameter constraints were necessary to prevent divergence in the calculations for the spectra at 245−297 K because of strong parameter correlations, especially between the linewidth and relaxation rate. The results indicate that the electron-hopping rate at room temperature is comparable to or faster than the fast-exchange limit (~10−8 s) of the Mössbauer time window. On the other hand, the fluctuation rate below 144 K is slower than the lower limit (τ = ~3 × 10−7 s) of the time window, and the electron-hopping relaxation model no longer describes the system (Fig. 2c).

Fig. 2

Mössbauer spectra analysis of the valence and spin state of the Fe ions. a–b Zero-field 57Fe Mössbauer spectra of a natural isotopic [Fe14] and a 57Fe-enriched sample on the N-terminal B site, [Fe857Fe6], at selected temperatures with fits considering the electron hopping between FeII-hs and FeIII-hs (purple area), FeII-hs (green area), FeIII-hs (orange area), and FeII-ls (striped yellow area). c Schematic of the valence behavior on the crystallographically equivalent B site in the Mössbauer time window. d Quadrupole splitting (QS) and isomer shift (IS) as a function of the temperature in the Mössbauer spectra. The values of FeII-hs (green) and FeIII-hs (orange) are determined from the 57Fe-enriched sample, whereas those of FeII-ls (striped yellow area) are determined from the natural sample. e Temperature dependence of the relaxation time (τ) of the electron hopping between FeIII-hs and FeII-hs in the temperature window of 164–297 K. The error bars were estimated as the respective deviations having the deference of the reduced chi-square value (Δχ2) within 2, considering the parameter correlations in the fitting. The red line represents the Arrhenius fitting in the high-temperature region.

Figure 2d summarizes the QS and IS values from the spectral analysis. A sudden change in the temperature dependences of IS and QS for FeII-hs and FeIII-hs was observed at around 220 K, whereas there was no such abrupt change for FeII-ls. This change seems to correlate with the electron hopping at the B sites of the complex. On cooling, the IS values of FeIII-hs and FeII-ls increased almost linearly according to a second-order Doppler shift, except for the sudden change. In contrast, the IS values of FeII-hs anomalously rise upon an increase in the temperature and do not obey a second-order Doppler shift effect at high temperatures. However, the average FeII-hs and FeIII-hs IS values at the B sites exhibit a linear shift with temperature. A similar behavior has previously been reported in a mixed-valence trinuclear Fe complex, [FeIIFeIII2O(CH3COO)(H2O)3][27], where the slope of the second-order Doppler shift was calculated as 6.3 × 10−4 mm s−1 T−1, comparable to the value of 5.6(3) × 10−4 mm s−1 K−1 for [Fe14]. The electron-hopping relaxation model provides a good fit in the temperature range of 164–297 K, with a relaxation time of τ as shown in Fig. 2e. A typical Mössbauer time window was observed to range from ~3 × 10−7 to ~1 × 10−8 s for the system. Actually, in a previous study using the relaxation model, it was indicated that such dynamics are detectable as spectral changes in the range of ~6.5 < −log(τ/s) < ~8.0[28], which is comparable to the case of [Fe14]. The electron-hopping activation energy (Ueff) was successfully evaluated using the Arrhenius equation, ln(τ) = ln(τ0) + Ueff/kBT. The best fit of the experimental data yielded Ueff = 943(47) cm−1 with a pre-exponential factor of τ0 = 5.0(16) × 10−11 s. Hence, the temperature dependence of the rates of intramolecular electron transfer in the Mössbauer spectroscopic analysis clearly demonstrated that six B-site Fe ions form a class II mixed-valence system according to the Robin–Day classification[29]. Electronic spectroscopy in the domain of intervalence transitions exhibits a wide peak centered at 11,100 cm−1 (900 nm) with a pronounced absorption tail due to the superposition of adjacent IVCT (FeII-CN-FeIII) and remote IVCT (FeII-NC-FeII-CN-FeIII), as indicated in the reported cyanide-bridged systems with N-terminal mixed-valence state[24,30,31]. (Supplementary Fig. 3)

The crystal structure was additionally analyzed to probe whether a temperature-induced charge separation occurs in [Fe14] at lower temperature, that is, 25 K, a temperature at which the electrons should be completely frozen (τ ~ 1013 s) according to aforementioned Arrhenius equation. However, no crystallographic differences in the coordination environment of the B-site Fe ions were observed between the structures at low and room temperature. At 25 K, the B-site Fe ions stay at the crystallographically identical sites in the unit cell as in the case of the RT (room temperature) structure. Therefore, although two FeII are expected to be located in the trans positions of the [Fe6] octahedron to minimize Coulombic repulsion, the positions of the FeII atoms cannot be assigned by X-ray diffraction. This is mostly due to the random distribution of the localized electrons at six B-site irons in the lattice. Thus, the charges at each B-site Fe virtually hold site occupation factors of 0.333 for FeII and 0.667 for FeIII, respectively. The heat capacity of [Fe14] under zero field at 7–300 K was investigated (Fig. 3a). No distinct anomaly was observed in the heat capacity data, further indicating no first- or second-order phase transition. The observation is consistent with the thermal dependence of the rate of intramolecular electron transfer as indicated in the Mössbauer spectra.

Fig. 3

Thermal dependence of the intramolecular electron transfer rate. a Zero-field molar heat capacity of the [Fe14] sample at temperatures of 7–300 K. The series 1, 2 and 3 denote the repeatability of the data in three-testing in the different temperature ranges. b In situ temperature-dependent infrared spectra of [Fe14] in the solid state. The dashed rectangles are guided for view to include the changing bands with temperature. c Stacked L-edge XAS for the [Fe14] complex measured at 300 and 3.5 K alongside the reference spectra of [FeIII(Tp)(CN)3] and [FeII(Tp)(CN)3]. The [FeII(Tp)(CN)3] overlaps with all features in c within the [Fe14] complex, whereas the [FeIII(Tp)(CN)3] does not (in particular the peak at ~706.8 eV associated with transitions into the t2g hole).

To elucidate the electron-hopping behavior in the [Fe14] complex, temperature-dependent infrared (IR) spectroscopic measurements were carried out at 78–350 K. The IR spectrum of [Fe14] at 78 K exhibits several strong νCN stretches at 2101, 2078, 2068, 2053, and 2042 cm−1 (Fig. 3b), indicative of the FeII-ls species bound to the C terminal of a cyanide bridge and both FeII and FeIII to the terminal N atom. Notably, only two νCN stretching bands were observed in the same region for a previously reported [Fe42] cluster, [{Fe(Tp)(CN)3}24{Fe(H2O)2}6{Fe(dpp)(H2O)}12(CF3SO3)6]·18H2O (dpp = 1,3-di(4-pyridyl)propane), wherein the valence-trapped FeII-ls–CN–FeIII-hs structure was preserved across the whole temperature range[32] (Supplementary Fig. 4). Furthermore, the νCN stretching frequency in Prussian blue with a FeII-ls–CN–FeIII-hs structure has been previously observed to be slightly higher compared with that in Prussian white with a FeII-ls–CN–FeII-hs structure[33]. The additional νCN stretching bands for [Fe14], which lack in [Fe42], indicates the presence of a FeII-ls–CN–FeIII-hs and FeII-ls–CN–FeII-hs mixed linkage, by in situ Fourier-transform infrared spectroscopy. These results mean that the electron-hopping rate is slower than the time scale of the IR technique (10−12–10−13 s)[34].

Iron L-edge X-ray absorption spectra (XAS) were measured between 3.5 and 300 K to further characterize the electronic structure (Fig. 3c). Spin–orbit coupling of the 2p5 final state (2p6 3d → 2p5 3d) leads to the splitting of L2,3 edges into L3 (J = 3/2) and L2 (J = 1/2) absorption regions, which are separated in energy by ~12 eV. The L3 and L2 absorption edges of [Fe14] were found to have maxima at 708.5 and 722 eV, respectively. The [Fe14] L edge includes overlapping contributions from the A and B sites. The XAS spectra measured between 3.5 and 300 K are identical, confirming that there is no resolvable change in the local Fe coordination symmetry or crystal field splitting at the A or B sites between the low- and high-temperature states. Furthermore, no temperature-induced changes in the spectra were observed for the Kβ1,3 emission spectroscopy and high-resolution K-edge X-ray absorption near-edge structure measurements (Supplementary Fig. 5). The lack of evidence for the thermal transition by these techniques indicates that the electron-hopping rate is slower than the time scale of the X-ray spectroscopy measurements (10−14 – 10−15 s). A comparison of the [Fe14] L-edge XAS data with the monomeric references [FeIII(Tp)(CN)3] and [FeII(Tp)(CN)3] indicates that the valences of the eight A sites are diamagnetic FeII-ls[35-38], confirming the Mössbauer assignment.

Magnetic properties

To assess the magnetic properties of [Fe14], variable-temperature direct current magnetic susceptibility data were collected. Figure 4a shows that at 300 K the χMT value is 27.45 cm3 K mol−1. Upon a decrease in the temperature, the χMT value shows a steady increase, before beginning to abruptly increase at ca. 30 K with a maximum of 59.17 cm3 K mol−1 at 3.5 K. This magnetic behavior suggests an overall ferromagnetic intramolecular exchange for [Fe14]. The data in the range of 300–10 K can be fitted according to the Curie–Weiss law, yielding C = 27.17 cm3 K mol−1 and θ = 2.95 K. This C value is comparable to the contribution from two S = 2 FeII-hs, four S = 5/2 FeIII-hs, and eight S = 0 FeII-ls centers in total for an individual [Fe14] complex, which is further supported by the evaluated experimental magnetic entropy 87.1 J K−1 mol−1 ≈ R ln(52 × 64) (=86.4 J K−1 mol−1, expected magnetic entropy). The positive Weiss constant suggests an intramolecular ferromagnetic interaction[39]. Density functional theoretical (DFT) calculations revealed that the spin density is mainly concentrated on the six B-site Fe ions (avg. spin value of 3.535) rather than the A-site Fe ions (avg. spin value of 0.232) (Supplementary Fig. 6), which is consistent with the experimental observations. To determine the spin ground state of [Fe14], magnetization data were collected and were plotted as reduced magnetization (M/NμB) vs. applied field in the range of 0–5 T at 2 K, as shown in Fig. 4b. The onset of magnetization in the applied field is greater than that of six magnetically isolated centers (red line in Fig. 4b) and is closer to the Brillouin curve for one isotropic S = 14 center (blue line in Fig. 4b). Inclusion of weak ferromagnetic exchange coupling (green line) between B sites suitably reproduces the measured magnetization curve, indicating a ferromagnetic S = 14 ground state.

Fig. 4

Magnetic characterization. a Temperature dependence of χMT and χM−1 for [Fe14] (Hdc = 100 Oe). The black symbol is the experimental data and the red line is the fitting result. b Magnetization vs. external magnetic field curve for [Fe14] at 2 K. The blue line corresponds to the calculated Brillouin function for S = 14 (J → +∞). The red line represents the sum of the Brillouin functions that correspond to four magnetically isolated S = 5/2 and two S = 2 (J = 0). The green line represents the simulated curve that includes weak ferromagnetic coupling between B sites (J = 0.9 K, g = 2.0). c Top: the left and right circularly polarized L3-edge total fluorescence yield spectra of the [Fe14] complex (14 T, 3.5 K). Bottom: the field dependence of the [Fe14] L3-edge XMCD spectra (left–right) measured at 3.5 K. d Magnetic susceptibility under a 50 Oe field in the temperature region of 50–0.5 K. e The field-dependent heat capacity of [Fe14] from 0 to 9 T in the temperature region of 100–0.36 K.

The presence of a ferromagnetic interaction was also confirmed using X-ray magnetic circular dichroism (XMCD), a technique that is a probe of magnetization at the atomic level. Since they are diamagnetic, the A-site Fe ions have no XMCD intensity. Hence, the [Fe14] XMCD spectrum provides direct spectroscopic access to only the B-site Fe ions. The left and right circularly polarized total fluorescence yield (TFY) detected XAS at 3.5 K, and 14 T is shown in Fig. 4c, alongside the XMCD spectra at various applied magnetic fields. At 3.5 K, the B sites are in the valence-trapped state in a ratio of 4:2 (FeIII vs. FeII). The onset of the L3 edge for the octahedral FeII-hs component is known to be lower in energy than that of FeIII-hs by ~1 to 2 eV. Hence, the lower energy region of the L3 XMCD spectrum is expected to exhibit FeII-hs contributions with respect to the higher energy region, which is expected to show FeIII-hs contributions. The magnitude of the XMCD signal across the L3 edge increases monotonically with an applied field indicating that the ferrous and ferric spins at the B sites are ferromagnetically exchange coupled[40]. (Supplementary Fig. 7) The field dependence of the XMCD data is in accordance with that observed for the superconducting quantum interference device (SQUID) magnetization.

It should be noted that the maximum χMT at low temperatures is significantly reduced compared to the expected χMT value for an isolated S = 14 ground state (~105 cm3 K mol−1). (Supplementary Fig. 8) The susceptibility measurements conducted down to 0.5 K exhibits a sharp peak in χM at 0.8 K (Fig. 4d), showing a typical feature of antiferromagnetic long-range magnetic ordering. Heat capacity measurements were carried out and showed a λ-type anomaly centered at TN = 0.85 K in the zero field (Fig. 4e), revealing the onset of weak long-range magnetic ordering. The magnetic nature of this feature was confirmed by its disappearance upon the application of a strong magnetic field. The intermolecular interaction was also suggested by the high-frequency electron paramagnetic resonance spectra of [Fe14] (Supplementary Fig. 9). Since the high-spin sites are coordinated via diamagnetic [FeII-ls(Tp)(CN)3], the mechanism of intramolecular ferromagnetic exchange over the atomic layer of [Fe14] necessitates the delocalization of the t2g [FeII-ls(Tp)(CN)3] metal character to B-site t2g orbitals via bridging CN π* orbitals. A theoretical description of the intramolecular exchange on the atomic layer of [Fe14] requires consideration of the electron delocalization, as described by Mayoh and Day[41,42]. In practice, this requires the evaluation of electron transfer integrals and interelectronic Coulombic repulsions, in conjunction with the Heisenberg contribution to the exchange for charge transfer states[43]. However, the value of such an analysis is limited due to the onset of 3D magnetic ordering at low temperature.

In summary, we report on the synthesis and electronic behavior of a nanosize [Fe14] complex with the two-electron, four-hole mixed-valence state. Two extra 3d electrons in the complex are found to hop at N-terminal Fe centers on its atomically thin spheric surface. The valence of the A-site Fe is basically static, while the valence electrons on the B-site Fe exhibit thermal dependence of the rate of intramolecular electron transfer. Furthermore, [Fe14] has a high-spin ground state with S = 14, meaning that the two electrons are hopping around the exchange-coupled atom-layer thin surface. The nanoarchitecture in this work may be useful for application in future molecular electronic and chemical devices using the intermetallic electronic and magnetic interactions in the framework and precise nanospace.

Methods

Synthesis of the [Fe14] complex

A 2 mL aqueous solution containing 6 μmol of Fe(BF4)2·6H2O was layered over a 2 mL DMSO aqueous solution containing 8 μmol of Bu4N[Fe(Tp)(CN)3]. A 1 μmol aqueous solution of l-ascorbic acid was utilized as the middle buffer layer under aerobic conditions. Crystals suitable for X-ray diffraction analysis were obtained in a yield of 16% after a week. Anal. calcd. for C108H128B8Fe14N72O12S6: C, 35.18; H, 3.50; N, 27.35. Found: C, 35.12; H, 3.52; N, 27.40. The analysis of Fe content in the solid sample of [Fe14] was performed through inductively coupled plasma atomic emission spectroscopy. The experiment result is 21.27 wt.%, which is in good agreement with the calculated value of 21.20 wt.%.

The 57Fe isotopically enriched sample, [Fe857Fe6], was prepared in a similar manner to that of [Fe14], except that the isotopically enriched salt, 57FeCl2·4H2O (57Fe, 96%), was utilized instead of Fe(BF4)2·6H2O. Through X-ray diffraction measurements, the crystal structure of the isotopically enriched sample, [Fe857Fe6], was observed to be identical to that of the natural isotopic [Fe14].

X-ray crystal structure determination

X-ray diffraction data at room temperature and 123 K for [Fe14] were collected on a BRUKER APEX-II CCD (Bruker Corp.) equipped with a graphite-monochromated Mo-Kα radiation source (λ = 0.71073 Å)[44-47]. Diffraction data at 25 K were collected under a cold helium gas stream on a Rigaku HPC X-ray diffractometer, using multi-layer mirror monochromated Mo-Kα radiation (λ = 0.71073 Å). Bragg spots were integrated using the CrysAlisPro program package, and empirical absorption correction (multi-scan) was applied using the SCALE3 ABSPACK program. The structures were solved by direct methods (SHELXT Version 2014/4) using full-matrix least-squares refinement (SHELXL Version 2018/1)[48]. The H atoms were geometrically placed on organic ligands in riding mode, and all of the non-H atoms were anisotropically refined by full-matrix least-squares refinement on F2 using the SHELXTL program[49]. A summary of the crystallographic data and refinement parameters is presented in Extended Data Table 1.

Computational details

Spin-polarized DFT calculations were carried out using the PWscf module in the Quantum Espresso 6.1 program package[50]. The exchange and correlation term in the Kohn−Sham equation was approximately treated using the Perdew–Burke–Ernzerhof method in terms of the gradient of the electronic density[51]. Projector-augmented wave potentials[52] were used along with a plane-wave basis set with a kinetic energy cut-off of 31.0 Ry. The simulation cell of [Fe(Tp)(CN)3]8[Fe(H2O)(DMSO)]6 contained 348 atoms in total. Atomic optimizations were carried out, starting from the experimental cell parameters. Integration in the first Brillouin zone for geometry optimizations was performed using 2 × 2 × 2 point Monkhorst–Pack sampling[53]. The SCF (self-consistent field) convergence and the total force convergence were set to be 1.0 × 10−6 (Ry) and 1.0 × 10−3 (Ry/au), respectively. The total magnetization was constrained to be 28.00 Bohr mag per cell. The spin-polarized DFT calculations further provided reliable information as to the charge and the spin state on each iron atom in [Fe14]. The singly occupied electrons at the B-site FeIII components partially delocalizes on the A-site FeII ions.

Magnetic analysis

The magnetic measurements of the samples were performed using a SQUID (MPMS-5S) magnetometer (Quantum Design Inc., USA). The magnetic susceptibility measurements shown in Fig. 4d were performed using MPMS-XL7AC (Quantum Design Inc., USA) apparatus with a 9 mm diameter 3He insert[54]. The data were corrected for diamagnetic contributions, calculated using Pascal’s constants[55].

Simulation of magnetization

The field dependence of magnetization for [Fe14] follows neither the Brillouin function for an S = 14 spin moment or the sum of uncoupled moments (four S = 5/2 and two S = 2). Inclusion of a weak ferromagnetic exchange interaction, acting between high spin sites was found to reproduce the measured magnetization curve. The magnetization curve was fit assuming one exchange constant J = 0.9 K and g = 2.0, based on the following Hamiltonian:where S1 and S6 represent the S = 2 sites and S2, S3, S4 and S5 represent the S = 5/2 sites, B is the applied magnetic field and µB is the Bohr magneton. This exchange model neglects the double-exchange component to the Hamiltonian present in Class II systems and effects due to dipolar fields from neighboring molecules. However, this simplified model indicates how the measured curve is consistent with weak ferromagnetic exchange coupling within [Fe14].

Fe Mössbauer spectroscopy

The 57Fe Mössbauer spectra were measured using a conventional Mössbauer spectrometer (Topologic Systems, Kanagawa, Japan) in transmission mode with a 57Co/Rh γ-ray source. Low-temperature measurements were performed upon a CryoMini/CryoStat cryogenic refrigerator set (Iwatani Industrial Gases, Osaka, Japan). The samples were tightly sealed with silicon grease in an acrylic holder and the spectra were calibrated using α-Fe foil as a reference at room temperature. The spectral fitting was carried out using the MossWinn 4.0 program and the full zero-field 57Fe Mössbauer spectra at all investigated temperatures were provided in Supplementary Fig. 10. For the 57Fe-enriched [Fe857Fe6] sample, the spectra were analyzed by applying an electron-hopping relaxation model (164–297 K; see Supplementary Method and Supplementary Fig. 11) or two quadrupole doublets (10–144 K). The area ratio was fixed at the ideal value (FeII/FeIII = 1/2) according to the chemical formula in order to avoid overparameterization. Some parameter correlations were found, especially between the linewidth and relaxation rate, through curve fittings of the spectra at 245–297 K. Therefore, the linewidth and IS values for FeII-hs at the relevant temperatures were fixed using those at 224 K. For the analyses at 10–144 K, the FeII-hs doublet was regarded as an asymmetric doublet, rising as a result of paramagnetic relaxation because an alternative symmetric doublet was tried but presented no sufficient result. For the natural isotopic sample of [Fe14], the Mössbauer spectra were analyzed using an additional doublet (FeII-ls; site A) together with the corresponding FeII-hs and FeIII-hs signals (site B) reproduced by the fixed parameters obtained from [Fe857Fe6]. Across the entire temperature range, the FeII-ls/FeIII-hs/FeII-hs ratio of the sample was maintained at 8/4/2.

Heat capacity calorimetry

Heat capacity measurements were performed with a laboratory-made adiabatic microcalorimeter in the temperature range of 9–300 K (adiabatic method) and with a PPMS (Quantum Design Inc., USA) in the temperature range 0.36–100 K under magnetic fields of 0–9 T (relaxation method). In the adiabatic calorimetry, 0.06356 g of a polycrystalline sample, which was made a buoyancy correction, was loaded into a 0.09 cm3 gold-plated copper cell and sealed with an indium wire under helium gas atmosphere. Thermometry was performed using a rhodium–iron alloy resistance thermometer (nominal 27Ω, Oxford Instruments) calibrated on the basis of the international temperature scale of 1990 (ITS-90). In the relaxation calorimetry (PPMS), we used buoyancy-corrected 1.1031 mg of a polycrystalline sample formed into a pellet of 2.5 mm in diameter. For the measurements below 10 K, a 3He insert was employed[56].

X-ray absorption spectra and X-ray magnetic circular dichroism

XAS and XMCD measurements at the Fe L absorption edges (703–740 eV) were measured on beamline I10 at the synchrotron Diamond Light Source of the Harwell Science and Innovation Campus in Oxfordshire in the United Kingdom. The XMCD spectra were obtained by flipping the helicity of circularly polarized X-rays exhibiting a 100% degree of polarization in the case of fixed applied magnetic fields. The measurements were performed with the temperature of the sample holder being regulated between 3.5 and 300 K. The total electron yield was obtained by measuring the drain current of the sample, whereas the TFY was obtained using a photodiode. The powdered samples of [Fe14] were attached with indium, to a copper sample holder. Radiolysis was controlled through the attenuation of the incident X-ray flux to 7% of the optimized value. Multiple scans were performed at each sample location to maintain control of the radiolysis, which was indicated by an increase in intensity at the low-energy portion of the L3 edge. TFY-detected measurements were found to be less susceptible to radiolysis and were hence adopted for XMCD measurements.

High-frequency electron paramagnetic resonance

High-frequency electron paramagnetic resonance (HF-EPR) measurements were performed on a locally developed spectrometer at the Wuhan National High-magnetic Field Center with a pulsed magnetic field of up to 30 T.