Recent advances in single molecule magnetism of dysprosium-metallofullerenes

Encapsulation of dysprosium ions in fullerenes results in efficient air stable single molecule magnets, which can be used in preparation of various 1D, 2D, and 3D assemblies.


Introduction
In 1993 the discovery of quantum tunneling of magnetization in an Mn 12 -complex by Sessoli et al. started the ongoing quest to find better single molecule magnets (SMMs). 1 Better in this case means attaining slower relaxation of magnetization at ever increasing temperatures, since single molecule magnets show great promise toward applications such as high density data storage, quantum computing and spintronics.As the name suggests, the magnetic properties of SMMs are determined from the intramolecular spin structure and are scalable down to the single molecule level.To be able to exhibit SMM properties, a molecule should have a bistable magnetic ground state and a certain energy barrier preventing fast reorientation of the spins in the ground state doublet.† The energy barrier is strongly related to magnetic anisotropy, which therefore is a cornerstone of the SMM field.The first SMMs were multinuclear complexes of transition metals, such as Mn or Fe.Since the beginning of the 2000s, lanthanides have been recognized as viable building blocks of SMMs due to their strong singleion anisotropy, 2 and Dy has become the most popular metal for creating new SMMs.The highest temperature of magnetic hysteresis in SMMs exceeding the liquid nitrogen temperature has recently been achieved in Dy-metallocenium salts. 3ndohedral metallofullerenes (EMFs) with lanthanide ions entered the field in 2012 when single molecule magnetism was proven for DySc 2 N@C 80 . 4Fullerenes facilitate stabilization and protection from ambient conditions of otherwise impossible or unstable molecular configurations within the confines of their carbon cages.The magnetic anisotropy necessary for single molecule magnetism in lanthanide EMFs can be provided by negatively charged nonmetallic species inside the cage and the carbon atoms of the negatively charged fullerene cages themselves.Fullerenes are stable in air under ambient conditions.They feature high thermal stability, allowing evaporation under high-vacuum conditions and growth of thin films via sublimation.Besides, EMFs exhibit rich addition chemistry, allowing modification of the cage surface with various functional groups without disrupting the structure of the endohedral species. 5This combination of physical and chemical properties, the tunability of the structure of the endohedral magnetic species and a possibility to create functional materials makes EMF-SMMs attractive objects for research.In this frontier, we provide an overview on the recent advances in single molecule magnetism of Dy-containing EMFs.First, we will give a brief overview of the general aspects of synthesis and magnetic characterization of EMFs, and then proceed with the description of different EMF-SMM families, as well as the studies of 1D, 2D, and 3D assemblies of EMF-SMMs.

Synthesis and structures of EMFs
The preparation of endohedral metallofullerenes usually starts with the so-called Krätschmer-Huffmann synthesis, 6 modified to suit the requirements of the desired system.This means arc discharge evaporation of graphite electrodes which are filled with a precursor (usually metal or metal oxide) of the desired endohedral species.The evaporation takes place at currents around 100 A at a pressure around 100 mbar under a He atmosphere.Reactive gases (NH 3 , CH 4 , etc.) or addition of organic compounds may be employed depending on the specifically desired fullerene system. 7s fullerenes are soluble in various organic solvents, in the next step they are extracted from the soot produced from arc discharge evaporation, usually by Soxhlet extraction or boiling under reflux in a suitable solvent.The dissolved fullerenes can be separated by High Performance Liquid Chromatography (HPLC) in multiple steps.Separation is usually the most timeconsuming step as the arc discharge evaporation may yield hundreds of different species.Still, isomerically pure compounds are attainable through the use of specialized HPLC columns.Yields are the downside of the otherwise very remarkable fullerene systems, ranging somewhere in milligrams to tens of milligrams of pure compounds produced per year.Molecular structure elucidation of isolated EMFs is then accomplished with conventional approaches such as singlecrystal X-ray diffraction, or various spectroscopic techniques.
Based on the composition of the endohedral species, EMFs can be classified into two large groups.Conventional EMFs have only metal atoms inside the carbon cage; depending on the number of metal atoms they can be mono-, di-, or trimetallofullerenes.In clusterfullerenes the endohedral species also include non-metal atoms such as C, N, S, or O.The interactions within EMF molecules have a considerable ionic character as metal atoms transfer their valence electrons to the fullerene cage.In clusterfullerenes, nonmetal atoms also bear a substantial negative charge (formally, N 3− , S 2− , C 2 2− , etc.).Molecular structures of representative Dy-EMFs discussed in this review are shown in Fig. 1.Further details on the syntheses, structures and properties of EMFs can be found in a number of comprehensive reviews and monographs. 7,8

Magnetic characterization of EMFs as SMMs
As the SMM behaviour is rooted in the rate of spin relaxation, all parameters employed to characterize SMMs implicitly involve magneto-kinetic characteristics.The latter may be defined and measured in different ways, which results in a certain ambiguity through the literature and complicates the comparison of different SMMs among each other.‡ A few key characteristics may be employed to characterize EMF-SMMs.
The blocking temperature of magnetization, magnetic hysteresis and the relaxation time of magnetization are outlined in this section.For more details on characterization of SMMs the reader may refer to specialized books and reviews listed in ref.
1 and 2. The blocking temperature of magnetization (T B ) is measured via two temperature scans.For one scan, magnetization is measured when the sample is cooled down to the lowest possible temperature in a magnetic field of 0.1-0.2T. For the other scan the sample is cooled in zero field, then the field is turned on and the temperature dependence of magnetization is measured during the temperature increase.These measurements reveal the point where magnetic relaxation becomes fast on the timescale of the measurement: the two curves coincide above T B and deviate below T B .The curve measured for the zero-field cooled sample usually develops a peak with its maximum at T B .T B is slightly dependent on the rate of the temperature sweep; the values reported by our group and discussed hereafter were measured with a temperature sweep rate of 5 K min −1 in a field of 0.2 T. The relaxation time at T B defined this way is ca. 10 s.The blocking temperature should not be confused with another universal parameter, the 100 seconds blocking temperature (T B100 ), which, as the name suggests, marks the temperature at which magnetic relaxation takes 100 s.Hysteresis curves are measured at fixed temperatures by sweeping the magnetic field between negative and positive values.The form of the magnetization curves obtained can give hints towards the magnetic behavior of the sample.A sudden drop of magnetization at zero magnetic field for instance is a sign of tunneling of magnetization (QTM). 9The temperature at which the magnetic hysteresis loop closes may also be used to characterize SMMs, but since it strongly depends on the sweep rate, this definition creates a lot of confusion when the values determined by different groups are compared.When the magnetic sweep rate of 2.9 mT s −1 is used (as in many of our studies), the closing temperature of hysteresis is close to T B as defined above.
Finally, the relaxation times of magnetization τ m (or their inverse, relaxation rates) are measured by magnetizing the sample in an external magnetic field at a fixed temperature and then switching the field to zero or another value.Then the evolution of magnetization can be observed over time, and the decay curve is fitted with an exponential function.Very often, single or even double exponential functions cannot describe the measured decay curves, and stretched exponential fitting is commonly used.This method allows the measurement of τ m values longer than 10-100 s.Shorter relaxation times may be accessible via AC magnetometry.Unfortunately, the latter requires considerably larger sample amounts than DC magnetometry, and only a few EMF samples have been characterized by this technique.
Analysis of the temperature dependence of relaxation times yields a better understanding of the spin relaxation mechanism, which is important for the design of improved SMMs.The relaxation rate is treated as a sum of rates for different processes: The first term describes the temperature-independent QTM; the second term corresponds to the single-phonon direct process.A(H) is field-dependent because the phonon frequency corresponds to the Zeeman energy gap of opposite spins; n 1 = 1 but may deviate when a phonon bottleneck occurs.The third term describes the two-phonon Raman mechanism, and n 2 is typically in the range of 5-9.The last term describes the Arrhenius behavior, usually associated with the Orbach mechanism.U eff then corresponds to the energy of the excited spin state involved in the relaxation.Careful analysis of the temperature dependence of τ m facilitates the identification of the dominant relaxation mechanisms in different temperature ranges.This phenomenological approach to the relaxation of magnetization in SMMs goes back to the studies of spin-phonon relaxation in paramagnetic salts. 10Significant limitations of this theory have been recognized in the SMM community during the last few years, 11 but more refined approaches are yet under development.

Nitride clusterfullerenes
DySc 2 N@C 80 -I h § was the first endohedral fullerene proven to be a single molecule magnet. 4,13The compound shows hysteresis in SQUID magnetometry measurements up to 6 K and a blocking temperature of T B = 7 K (Fig. 2a and 3a).It could be shown that the Dy ion in the compound is responsible for its magnetic properties by comparing the magnetization curves attained by SQUID magnetometry and X-ray magnetic circular dichroism (XMCD) at the Dy M 5 edge.The peculiar "butterfly shape" of the magnetization curves is attributed to the quantum tunneling of magnetization, the relaxation mecha-Fig.2 (a) Magnetic hysteresis of non-diluted DySc 2 N@C 80 powder measured at 1.8 K compared to the sample diluted with the diamagnetic fullerene Lu 3 N@C 80 , absorbed in the metal-organic framework DUT-51 (Zr) (@MOF), and dispersed in polymer polystyrene (@PS).Strong variation of the QTM-induced drop of magnetization near zero-field with dilution can be seen.The inset shows that all samples have the same blocking temperature of 7 K.(b) Relaxation times of magnetization measured at 1.8 K in different magnetic fields for non-diluted powder and for diluted samples in MOF, polystyrene (PS) and in a single-crystal (non-diluted, SC, and diluted with Lu 3 N@C 80 , SC-dil).The inset zooms into the small field range.Reproduced from ref. 13. § Since fullerenes have many different isomers, the use of a certain nomenclature is necessary to distinguish different cages.A standard approach is to use Fowler-Manolopoulos spiral algorithm, 43 which yields a unique number for each cage.In Table 1 we denote fullerene isomers by the formal cage symmetry followed by the spiral number in parenthesis, in the text the spiral numbers are omitted for readability.nism common in SMMs with a single magnetic metal ion in the molecule (so-called single-ion magnets).Dilution of the sample with nonmagnetic C 60 was used to prove that the magnetic properties were indeed a single molecule phenomenon instead of a collective effect.It could also be demonstrated, that dilution increases the relaxation time in zero magnetic field.An in-depth investigation of the relaxation mechanisms in DySc 2 N@C 80 powders and single-crystals as well as diluted fullerenes in three different diamagnetic matrices was presented in 2018. 13Very careful measurements showed a strong influence of dilution on the field dependent relaxation mechanism.For instance, it was shown that strong dilution of the magnetic fullerenes in polystyrene decreases the QTM resonance from 150 mT in an undiluted sample to <1 mT, a feature easily missed in measurements with commercial equipment (Fig. 2).Additionally, measurements of zero field relaxation times revealed a slight temperature dependence of the QTM between 2 and 5 K, which was tentatively attributed to slow energy dissipation through the lattice.Finally, it is noteworthy that the expected linear dependence of log(τ m ) vs. T −1 (Arrhenius coordinates) could not be confirmed until 87 K, where a low signal to noise ratio of the AC magnetometry data ends the reliable measurement range (Fig. 3).Ab initio calculations showed that the Dy ion in the nitride cluster experiences a strong uniaxial ligand field with an overall splitting in the order of 1500 cm −1 . 19The quantization axis is aligned parallel to the Dy-N bond, and the ground state is a Kramers doublet with J z = ±15/2 separated from higher energy states by ca.400 cm −1 .Thus, equilibrium magnetic properties of Dy-nitride clusterfullerenes up to room temperature are essentially determined by the ground state doublet.Strong uniaxial anisotropy of lanthanide ions in the LnSc 2 N@C 80 -I h molecules was also confirmed by paramagnetic NMR spectroscopy and point-charge ligand-field calculations. 20Dy was found to impose the strongest paramagnetic shift on the 45 Sc nuclear spin in the whole lanthanide series.
The first investigation of the magnetic properties of dinuclear Dy 2 ScN@C 80 -I h and trinuclear Dy 3 N@C 80 -I h was published in 2014, 14a followed by a detailed study of temperature dependence in Dy 2 ScN@C 80 in 2017.14b Dy 2 ScN@C 80 exhibits magnetic hysteresis and blocking of magnetization at 8 K and does not show the QTM present in its single-ion counterpart, DySc 2 N@C 80 .This can be attributed to the ferromagnetic coupling of the two Dy spins in Dy 2 ScN@C 80 .Flipping one of the Dy spins brings the system into an antiferromagnetically coupled state, which is higher in energy than the ground state with ferromagnetic coupling by 10 K. 4.6 K, roughly half of this energy, is attributed to dipolar interactions, and the rest to exchange coupling.This barrier prevents zero-field QTM in Dy 2 ScN@C 80 .Indeed, relaxation times show Arrhenius behavior at low temperature with the U eff corresponding to the energy difference between the ferromagnetic and antiferromagnetic states (Fig. 3c), proving that relaxation proceeds via the latter state.AC magnetometry was used to investigate magnetic relaxation at higher temperatures and revealed an Orbach mechanism with an exceptionally high thermal barrier of 1735 ± 21 K (Fig. 3c).Ab initio calculations helped to assign this barrier to the relaxation via the 5 th Kramers doublet.
The third member of the series, Dy 3 N@C 80 -I h , does not show remanence (Fig. 3a), which can be attributed to a frustrated magnetic ground state.14a,19a Ferromagnetic coupling in a triangular Dy 3 N cluster cannot be realized for all three Dy spins at once since the single-ion quantization axes are linked to corresponding Dy-N bonds arranged at 120°with respect to each other.
The promising properties of the Dy x Sc 3−x N@C 80 -I h family led to increased interest in endohedral fullerene SMMs and subsequently many systems were synthesized and checked for their magnetic properties.Basically, three parameters which potentially affected the magnetic properties were identified: the magnetic species themselves, the size or the specific isomer of the encapsulating cage, and the nonmetallic species that might also be encapsulated in the fullerene.
The influence of the carbon cage size and its isomerism on the magnetic properties of encapsulated DySc 2 N and Dy 2 ScN clusters has been studied recently. 12DySc 2 N@C 68 -D 3 , DySc 2 N@C 80 -D 5h , DySc 2 N@C 80 -I h , D y 2 ScN@C 80 -D 5h , Dy 2 ScN@C 80 -I h and Dy 2 ScN@C 84 -C s were compared for their key characteristic markers.It could be demonstrated that the C 80 -I h cage isomer yields the SMM with the highest blocking temperature and slowest relaxation of magnetization.It was hypothesized that free movement of the clusters inside the cage and subsequent week spin-phonon coupling seem to be the strongest factor in enhancing the magnetic properties of EMFs.

Sulfide clusterfullerenes
As seen in Dy 2 ScN@C 80 , coupling the two Dy atoms via a nitride ion leads to the suppression of quantum tunneling of magnetization.Checking different non-metal units in Dy-clusterfullerenes was therefore a logical next step.The sulfide clusterfullerenes Dy 2 S@C 82 -C 3v , Dy 2 S@C 82 -C s , and Dy 2 S@C 72 -C s were synthesized and purified to study this effect. 15The sulfide clusterfullerene with C 3v cage symmetry was found to be the best SMM among the three, showing hysteresis, which closes between 4 and 5 K, and a blocking temperature T B of 4 K (Fig. 4).The study revealed considerable cage dependence of the magnetization dynamics.Additionally, in the C 3v isomer, three Orbach processes with different thermal barriers could be identified (Fig. 4c), governing the relaxation behavior at different temperatures, respectively.Among the few SMMs with sulfur-ligated Dy reported so far, 21 Dy 2 S-clusterfullerenes have the longest relaxation times and the highest relaxation barriers.The reason is that in Dy 2 S@C 2n molecules, sulfur bears a substantially larger negative charge (the formal charge state is S 2− ) and Dy-S distances are at the same time much shorter, which altogether leads to a strong axial ligand field.

Carbide clusterfullerenes
Clusters featuring a TivC double bond and two more metal atoms enable the synthesis of SMMs with up to three different metal species in one cage. 16,17,22This is rare as chromatographic separation becomes more complex with every endohedral species added.Remarkably, it was found that DyYTiC@C 80 -I h shows a relatively high blocking temperature of magnetization of 7 K, which is comparable to that of DySc 2 N@C 80 -I h . 16The hysteresis of Dy 2 TiC@C 80 -I h on the other hand closes at 3 K, which compares very poorly to that of Dy 2 ScN@C 80 -I h . 17A comparison of dinuclear carbide and nitride clusterfullerenes shows that the anisotropy introduced by the nonmetallic unit plays a smaller role in magnetic behavior at low temperature than the exchange interaction.This becomes even more pronounced in the congener of Dy 2 TiC, Dy 2 TiC 2 , with one more carbon atom in the endohedral cluster, which only shows very narrow hysteresis at 1.8 K. 17 Another carbide clusterfullerene, Dy 2 C 2 @C 82 -C s , which is isostructural to the aforementioned sulfide clusterfullerene Dy 2 S@C 82 -C s , shows similar magnetic properties to the latter. 15Fitting of AC magnetometry data revealed barriers of 15.2 K and 17.4 K, respectively, between their ferromagnetic ground states and what is presumed to be an antiferromagnetic excited state.The main difference appears in their respective attempt times τ 01 , which are 2.9 ms for Dy 2 S@C 82 -C s and 0.5 ms for Dy 2 C 2 @C 82 -C s , yielding shorter relaxation times for the latter, therefore making it the weaker SMM.
In conclusion, the studies on clusterfullerenes yielded the following insights: (1) Nitride clusterfullerenes give the best SMMs, followed by sulfide, C 1 -carbide and C 2 -carbide clusterfullerenes in that order.
(2) The non-metal clusters facilitate the strong single-ion anisotropy needed to make SMMs.In dinuclear EMFs they also contribute to the coupling of the magnetic ions, suppressing QTM and thus giving SMMs with pronounced remanence.
(3) The fullerene cage is not just an inert container, but also plays a role in the relaxation of magnetization, as evidenced by the variation in SMM properties for different cage sizes and isomers.
Ab initio calculations predicted that oxide clusterfullerenes have the largest crystal field splitting among clusterfullerenes, making them an interesting subject for future investigation. 15,23till, the interaction between magnetic ions coupled by nonmetallic atoms is relatively weak, giving an energy barrier between ferromagnetic and antiferromagnetic states of less than 15 K.

Dimetallofullerenes
To enhance the coupling between magnetic ions a covalent metal-metal bond presents the most elegant solution.Dimetallofullerenes (di-EMFs) proved uniquely suited to this end.In EMFs, lanthanide ions tend to give their valence electrons away and fullerene cages tend to act as electron acceptors.The metal-metal bonding molecular orbital is one of the frontier orbitals in di-EMFs, and its population depends on its energy in relation to the cage MOs. 24M 2 @C 82 di-EMFs (M = Sc, Y, Er, Lu) were found to have occupied M-M bonding MOs with formal charges of +2 on both metal ions.Electrochemical manipulation of the Er-Er bonding orbital in Er 2 @C 82 was shown to effectively change coupling by creating a three spin system {Er 3+ -e-Er 3+ }. 25 In C 80 -I h cages the valence MOs of the La 2 dimer have a relatively high energy, which leads to all valence electrons being transferred to the fullerene cage and yielding La ions with a formal charge of +3.However, a single-electron reduction of La 2 @C 80 gives the monoanion with a single-electron La-La bond as evidenced by EPR spectroscopy. 26On the other hand, Y 2 and intermediate lanthanide dimers such as Gd 2 , Dy 2 , or Lu 2 give only five electrons to the cage, leaving each ion with a formal charge of +2.5 and the highly soughtafter single electron bond between them.The downside of this is the formation of radicals, which are hard to extract from the soot due to polymerization in neutral solvents. 27Stabilization of this unique electronic configuration can be achieved by the substitution of a carbon atom in the cage by nitrogen, giving azafullerenes M 2 @C 79 N, 28 or by extraction in a polar solvent such as dimethylformamide (DMF) and subsequent functionalization of the cage with a radical group to form monoadducts M 2 @C 80 -R (R = CF 3 , 27,29 benzyl 18 ).
In 2015, a computational study by Sing et al. predicted a strong magnetic exchange in Gd 2 @C 79 N and a large magnetization relaxation barrier in Dy 2 @C 79 N. 30 Indeed, the coupling constant between Gd and the unpaired electron residing on the Gd-Gd bond in Gd 2 @C 79 N has recently been found to be 170-175 cm −1 . 31ltimately, the EMF-SMM record set by Dy 2 ScN@C 80 was broken by a benzyl monoadduct of Dy 2 @C 80 -I h , Dy 2 @C 80 (CH 2 Ph).Dy 2 @C 80 could be extracted from soot with DMF, presumably in the form of a monoanion, whereas non-polar solvents such as toluene or carbon disulfide did not work in this regard.Functionalization of the cage with a benzyl group by treatment with benzyl bromide afforded toluene-soluble air-stable molecules which could be isolated by HPLC.An in-depth investigation into its synthesis and properties was published in 2017. 18The SMM properties displayed by this molecule are truly remarkable with a blocking temperature of magnetization of 21.9 K and hysteresis observable between 1.8 and 21 K (Fig. 5a).The 100 seconds blocking temperature T B100 was determined to be 18 K.Relaxation of magnetization in Dy 2 @C 80 (CH 2 Ph) in zero field between 1.8 and 5 K proceeds via QTM with a relaxation time of 3257 s.When QTM is quenched by the application of a finite magnetic field, an Orbach-like process, attributed to phonon assisted relaxation, with an effective barrier of 40 K starts at 3 K and becomes dominant between 10 and 18 K.Above 20 K another Orbach process with U eff = 613 K takes over (Fig. 5b).Dy 2 @C 80 (CH 2 Ph) can be treated as a three spin system {Dy 3+ -e-Dy 3+ }, with the Dy ions coupling ferromagnetically to the electron from both sides (Fig. 5c).Ab initio calculations were used to show that the negative charge between the positively charged ions enforces easy axis anisotropy along the Dy-Dy bond.The direct antiferromagnetic coupling between the Dy ions is very weak and can be neglected.Determination of the coupling constants between the Dy ions and the electron spin between them was achieved by simulating magnetization and χT curves to match experimental data, with a very large value of j Dy,e = 32 cm −1 or 46 K giving the best fit.Assuming this coupling constant, the energy of the exchange excited state, in which one Dy spin is flipped, was calculated to be 613 K.This energy matches the high-temperature Orbach barrier, determined from fitting magnetic relaxation data.The successful isolation of this elusive class of EMFs marks a breakthrough for the field.

Beyond powder samples
The magnetic properties described in the previous section were obtained for bulk powder EMF samples.This is a first step in magnetic characterization of SMMs, but for the evaluation of their potential applications the studies of powder samples are insufficient.Playing to the strengths of SMMs, addressability of single molecules would be needed, which should be most easily achievable in 1D or 2D arrays.The welldefined positioning in 3D matrices is desirable as well to finetune the properties of future nanomaterials.

1D arrays: peapods
Single walled carbon nanotubes (SWCNTs) can provide a channel, in which endohedral fullerenes may line up in a 1D chain. 32These structures are commonly known as peapods and show great promise for applications in spintronic devices and quantum computation.
A study on chains of Dy 2 ScN@C 80 inside SWCNTs by XMCD was published recently by Avdoshenko et al. 33 A comparison of the encapsulated EMFs with a powder sample of the same composition revealed a reduction of the magnetic bistability caused by the encapsulation.Additionally, partial ordering of the clusters was observed (Fig. 6a and b).An explanation of this behavior was given with the help of calculations on isostructural Y 2 ScN@C 80 packed in SWCNTs.Depending on the relationship between the cage size and diameter of the used SWCNT, energetically preferable orientations of the clusters emerged.
DySc 2 N@C 80 was encapsulated in SWCNTs by Nakanishi et al. 34 (Fig. 6c).Here hysteresis was conserved and an increase in coercivity and a longer relaxation time compared to those of the powder sample could be observed (Fig. 6d and e).Thus, encapsulation within the SWCNT partially suppressed the QTM relaxation of DySc 2 N@C 80 , and the authors attributed this to a dilution effect which is also observable in bulk powder samples.

2D arrays: (sub)monolayers on substrates
Deposition of SMM molecules on conducting surfaces is an obvious route to their addressable 2D arrays. 35However, magnetic bistability in monolayers was observed for only a few of the many substances that show SMM behavior in bulk samples.The main difficulties in this route are caused by insufficient thermal or chemical stability of SMM molecules precluding the formation of monolayers and the detrimental effect of molecule-metal interactions on the SMM properties.
In 2014, Westerström et al. published a study of Dy 2 ScN@C 80 deposited onto a Rh(111) surface by evaporation under vacuum.36a In a submonolayer, ordering of the magnetic moments on the surface and hysteresis of magnetization could be observed at 4 K (Fig. 7).The relaxation time was estimated to be approximately 16 times faster than that for powder samples, although demagnetization by X-ray irradiation should be taken into account, 37 as the results were obtained by XMCD with synchrotron radiation.A recent XMCD study showed that deposition of Dy 2 ScN@C 80 onto h-BN/Rh(111) nanomesh resulted in a broader hysteresis than on a pure Rh(111) surface.36b Chemical functionalization of Dy 2 ScN@C 80 and DySc 2 N@C 80 with a thioether group was achieved via 1,3-dipolar cycloaddition by Chen et al. 38 The SMM behavior was observed in functionalized EMFs, but the magnetic properties changed noticeably in comparison with those of pristine EMFs (Fig. 8ac).The blocking temperature T B was increased by 1 K for DySc 2 N@C 80 , but decreased by 4 K for Dy 2 ScN@C 80 .The coercive field for functionalized Dy 2 ScN@C 80 was also visibly lower compared to that of the non-functionalized sample.The functionalized molecules were then able to attach to an Au(111) surface by physisorption.When deposited on gold both functionalized molecules showed hysteresis of magnetization at 2 K as proven by XMCD (Fig. 8d and e).DFT calculations showed that a horizontal configuration of the functionalized molecules, with the fullerenes touching the metal, is energetically favorable over a vertical configuration.Along with X-ray induced demagnetization this might be responsible for the clearly shortened relaxation times.Additionally, the structures are highly mobile at room temperature, leading to random orientations on the surface.Further analysis revealed a certain protective property of the cage π-system for the magnetic state of the cluster, even when the cage interacts strongly with the metal surface.

3D arrays: MOFs
Metal organic frameworks (MOFs) provide highly ordered porous structures, which can accommodate SMM molecules. 39ang et al. introduced DySc 2 N@C 80 into the pores of MOF-177 and observed a suppression of QTM. 40Similar suppression of zero-field QTM was observed when DySc 2 N@C 80 was incorporated into the pores of DUT-51(Zr) 13,41 as mentioned above in  the discussion of the dilution effect on the QTM of DySc 2 N@C 80 .The distance between EMF molecules in the MOF is considerably longer than those in the powder EMF samples, which leads to much weaker dipolar magnetic fields and hence to a narrowing of the QTM resonance (i.e. to the decrease of the field range in which QTM can take place).Suppression of QTM in DySc 2 N@C 80 was also observed when it was encapsulated within the pores of an azobenzenefunctionalized MOF. 42The authors claimed that irradiation of the Azo MOF with light causing trans-cis isomerization of azobenzene moieties improved the SMM properties of absorbed DySc 2 N@C 80 .

Concluding remarks
Over the last few years Dy containing EMFs have been proven to comprise robust SMMs.Numerous recent studies have con-tributed to the understanding of the forces that govern their magnetic properties.Through this understanding, control over anisotropy and intramolecular interactions is attainable, by the choice of non-metal species as well as cage sizes and isomers.While among clusterfullerenes the nitrides show the strongest magnetic properties, followed by sulfides and carbides, oxides may be promising as well.The successful isolation of Dy 2 @C 80 -CH 2 Ph with a single electron bond gives a new direction to the field, providing molecules with the highest blocking temperatures measured for EMF SMMs and among the highest for SMMs in general.The high stability and protection provided by the fullerene cages make for perfect prospects towards use in future devices, as this enhances processability.Initial studies have already shown various routes toward 1D, 2D and 3D structures.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Fig. 1
Fig. 1 Molecular structures of selected Dy metallofullerenes showing single molecule magnetism.Dy is shown in green, Scmagenta, Ticyan, Yviolet, Nblue, Cgray, and Syellow.Only a part of the benzyl group of Dy 2 @C 80 (CH 2 Ph) can be seen.

Fig. 3
Fig. 3 (a) Hysteresis curves for Dy x Sc 3−x N@C 80 -I h (from left to right: x = 1, 2, and 3) recorded using SQUID magnetometry at 2 K with a field sweep rate of 0.8 mT s −1 .Reprinted with permission from Westerström et al., Phys.Rev. B: Condens.Matter Mater.Phys., 2014, 89, 060406.Copyright 2014 by the American Physical Society.(b) Relaxation times of magnetization of DySc 2 N@C 80 at temperatures of 2-87 K. Zero-field values are shown as full dots, and in-field (0.2 T) values are denoted as open dots.Relaxation times for non-diluted DySc 2 N@C 80 are shown in black, and the values for diluted samples are shown in blue (dilution with MOF) and green (diluted with polystyrene, PS).The times longer or shorter than 10 s were determined by DC and AC magnetometry, respectively.The blue line is the fit of the points in the 2-5 K range with the Orbach relaxation mechanism, and the black line represents the fit of the QTM-like zero-field relaxation with the power function of temperature.Reproduced from ref. 13.(c) Relaxation times of the magnetization of Dy 2 ScN@C 80 .Green dots denote the values from DC measurements in zero field; two in-field points (red crosses) are also shown.AC values are measured with MPMS XL (7-50 K; open, magenta, and blue dots) and with PPMS (brown dots, 52-76 K).Magenta and blue dots denote long and short times from double-τ fits of the AC data, respectively, and open dots denote single-τ fits.Reproduced from ref. 14b.

Fig. 4
Fig. 4 Magnetization curves for (a) Dy 2 S@C 82 -C s and (b) Dy 2 S@C 82 -C 3v measured at T = 1.8-5K with a magnetic field sweep rate of 8.33 mT s −1 .The inset in each panel zooms into the region near zero-field.The inset in (b) shows the determination of T B of Dy 2 S@C 82 -C 3v from the peak in the susceptibility of the zero-field cooled sample (magnetic field: 0.2 T, temperature sweep rate: 5 K min −1 ).(c) Magnetization relaxation times of Dy 2 S@C 82 -C 3v ; dots are experimental points, red lines are results of a global fit with three Orbach processes; and green, magenta, and brown lines represent contributions of individual Orbach processes.The inset shows an enhancement of the high-temperature range.Reproduced from ref. 15.
This journal is © The Royal Society of Chemistry 2019 Open Access Article.Published on 13 February 2019.Downloaded on 12/19/2019 9:29:42 PM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Fig. 5
Fig. 5 (a) Magnetic hysteresis in Dy 2 @C 80 (CH 2 Ph) between 2 and 22 K, field sweep rate: 2.9 mT s −1 .(b) Magnetization relaxation times of Dy 2 @C 80 (CH 2 Ph) in zero-field and in a field of 0.4 T. The inset shows the out-of-phase dynamic susceptibility χ'' measured at different temperatures between 23 and 33 K. (c) Alignment of magnetic moments in the ground state of Dy 2 @C 80 (CH 2 Ph) and respective spin Hamiltonian (CF denotes the crystal field).Reproduced from ref. 18.

Fig. 6
Fig. 6 (a) X-ray absorption spectra of Dy 2 ScN@C 80 encapsulated in SWCNTs recorded using right (I + ) and left (I − ) circularly polarized X-rays.(b) A comparison of the normalized total absorption and XMCD spectra from bulk Dy 2 ScN@C 80 and Dy 2 ScN@C 80 encapsulated in SWCNTs.The temperature is 2 K, and an external magnetic field of 6.5 T is applied parallel to the X-ray beam and the surface normal to the samples.(c) TEM image and structural model of the [DySc 2 N@C 80 ]@SWCNT peapod.(d) Magnetization curves of [DySc 2 N@C 80 ]@SWCNT measured at different temperatures by SQUID magnetometry; (e) comparison of magnetic hysteresis curves for bulk DySc 2 N@C 80 and [DySc 2 N@C 80 ]@SWCNT peapod (T = 1.8 K).(a) and (b) reproduced from ref. 33.Reprinted with permission from Nakanishi et al., J. Am.Chem.Soc., 2018, 140, 10955.Copyright 2018 by the American Chemical Society.

Fig. 7
Fig. 7 (a) Sub-monolayer (ML) of Dy 2 ScN@C 80 /Rh(111), T = 4 K, μ 0 H = 6.5 T; measurement geometry is shown in the inset.The polarization dependent X-ray absorption spectra (left panel), and the corresponding XMCD spectra (right panel) measured at incidence angles of θ = 0°and θ = 60°.Strong angular dependence points to the preferential alignment of Dy spins parallel to the surface.(b, c) Hysteresis curves measured by XMCD from a multilayer (b) and a sub-ML (c) of Dy 2 ScN@C 80 /Rh(111) at a magnetic field sweep rate of 2 T min −1 and a sample temperature of ∼4 K.The drop in magnetization at zero field is a consequence of the time of 30 s it takes the magnet to switch polarity.Reprinted with permission from Westerström et al., Phys.Rev. Lett., 2015, 114, 087201.Copyright 2015 by the American Physical Society.

Fig. 8
Fig. 8 (a) Scheme of a Prato reaction to obtain EMF-R derivatives (EMF = DySc 2 N@C 80 (1), and Dy 2 ScN@C 80 (2), R denotes the functional group with a thioether linker).(b, c) Magnetization curves of (b) 1-R and 1, and (c) 2-R and 2 measured by SQUID magnetometry at T = 2 K (field sweep rate: 2.9 mT s −1 ); the insets show determination of the blocking temperatures of magnetization T B (temperature sweep rate: 5 K min −1 ).(d, e) magnetization curves of sub-monolayers of 1-R (d) and 2-R (e) on Au (111) measured by XMCD at 2 K with a sweep rate of 2 T min −1 (averaging over five measured curves, and error bars are standard deviations).Reproduced from ref. 38.

Table 1 T
B and T B100 parameters of Dy EMF-SMMs This journal is © The Royal Society of Chemistry 2019 Open Access Article.Published on 13 February 2019.Downloaded on 12/19/2019 9:29:42 PM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.