Studying the effect of Zn-substitution on the magnetic and hyperthermic properties of cobalt ferrite nanoparticles

The possibility to finely control nanostructured cubic ferrites (M(II)Fe2O4) paves the way to design materials with the desired magnetic properties for specific applications. However, the strict and complex interrelation among the chemical composition, size, polydispersity, shape and surface coating renders their correlation with the magnetic properties not trivial to predict. In this context, this work aims to discuss the magnetic properties and the heating abilities of Zn-substituted cobalt ferrite nanoparticles with different zinc contents (ZnxCo1-xFe2O4 with 0 < x < 0.6), specifically prepared with similar particle sizes (∼7 nm) and size distributions having the crystallite size (∼6 nm) and capping agent amount of 15%. All samples have high saturation magnetisation (Ms) values at 5 K (>100 emu g(-1)). The increase in the zinc content up to x = 0.46 in the structure has resulted in an increase of the saturation magnetisation (Ms) at 5 K. High Ms values have also been revealed at room temperature (∼90 emu g(-1)) for both CoFe2O4 and Zn0.30Co0.70Fe2O4 samples and their heating ability has been tested. Despite a similar saturation magnetisation, the specific absorption rate value for the cobalt ferrite is three times higher than the Zn-substituted one. DC magnetometry results were not sufficient to justify these data, the experimental conditions of SAR and static measurements being quite different. The synergic combination of DC with AC magnetometry and (57)Fe Mössbauer spectroscopy represents a powerful tool to get new insights into the design of suitable heat mediators for magnetic fluid hyperthermia.


Introduction
Nanostructured metal oxides of the formula M II Fe 2 O 4 (M II = Fe 2+ , Co 2+ , Mn 2+ , Ni 2+ , etc.), commonly referred to as cubic ferrites, are one of the most important materials for their ferrimagnetic properties closely related to their spinel structure. In the spinel structure, oxygen ions are close-packed in a cubic arrangement (ccp) in which metallic cations occupy one out of the eight tetrahedral interstices (indicated as round brackets) and half of the octahedral ones (squared brackets). Tetrahedral sites are filled by divalent cations in a direct or normal spinel (M II )[M III ] 2 O 4 , whereas trivalent cations replaced them in an inverse spinel (M III )[M II ;M III ] 2 O 4 . 1 Actually, real systems crystallise in a partially inverse structure with the following general formula: where δ is a parameter related to the inversion degree (i = 1−δ) and it is equal to 1 for a normal spinel and 0 for an inverse one. 2 The nature of M II strongly affects the magnetic properties of the material, such as magnetisation or anisotropy. 3 CoFe 2 O 4 shows a predominantly inverse structure with Co 2+ ions mainly on octahedral sites (square bracket) and Fe 3+ ions almost equally located between octahedral and tetrahedral sites (round bracket) but the observed inversion degree is often lower than 1. 4 ZnFe 2 O 4 is assumed to be a normal spinel with all Fe 3+ ions on octahedral sites and all Zn 2+ on tetrahedral sites. As a consequence of their structure, CoFe 2 O 4 is ferrimagnetic below 860 K, while ZnFe 2 O 4 is antiferromagnetic below 9 K. As zinc substitutes cobalt in the cobalt ferrite structure, Zn 2+ ions tend to preferentially occupy tetrahedral positions while Fe 3+ in tetrahedral positions should move toward octahedral positions. However, due to the competition between the different ions for the two sites the cation distributions have to be determined experimentally. Therefore, changing the cationic distribution together with the chemical composition represents one of the possible strategies for modulating the magnetic behaviour due to the strict relationship existing between the spinel structure and its magnetism. 3,[5][6][7] For a given composition, it is worth noting that the physical properties can be tuned also as a function of the size, 3,8,9 the shape 10,11 and the capping agent. 3,11 Other phenomena such as the surface and internal non-collinear spin structure (spin-canting) can further influence the magnetic behaviour. Spinel ferrites, being one of the most versatile systems, have found applications in a wide variety of fields that include data storage, catalysis, energy, environment, and in particular, biomedicine. Indeed, the use of magnetic nanoparticles for biomedical purposes has been proposed to a large extent in recent years. 12,13 Early diagnosis and targeted therapies are key challenges in the field of biomedicine. The development of new technologies designed to answer to these issues has been a strong driving force in this research field. Different biomedical applications of magnetic nanoparticles have been developed: magnetic separation, magnetic resonance imaging (MRI), 14 drug delivery 15 and magnetic fluid hyperthermia (MFH). 13,16 In particular, MFH is based on the ability of magnetic nanoparticles, when an alternate external magnetic field is applied, to convert the electromagnetic energy into heat. Different heat release mechanisms are involved depending on the magnetic behaviour: superparamagnetic or ferromagnetic. Néel and Brownian relaxation mechanisms are responsible for the heat delivery of superparamagnetic nanoparticles, both of them involving the reversal of the magnetisation vector inside the nanoparticle (Néel) or through the physical rotation of the particles themselves (Brownian). Therefore, the efficiency of magnetic nanoparticles as heat mediators is strictly related to their magnetic properties (saturation magnetisation and reversal energy barrier), besides the hydrodynamic ones depending on the environmental medium (viscosity, temperature). The most studied materials for biomedical applications are superparamagnetic iron oxides (SPIOs), i.e. maghemite (γ-Fe 2 O 3 ) and magnetite (Fe 3 O 4 ), nanoparticles, due to their biocompatibility (approved by U.S. Food and Drug Administration) and their low-cost. Nevertheless, much of the research on this topic is still devoted to the optimisation of the magnetic properties aimed to increase the heating ability and to reduce the nanoparticle dose to be injected into the human body. In particular, much effort points towards an increase of the saturation magnetisation 17 investigating as alternative metalbased systems, [18][19][20][21] and substituted-ferrite nanoparticles. 22,23 Other strategies based on the coupling of ferrimagnetic spinels by building appropriate core-shell heterostructures have been proposed. 24 Among cubic ferrites, CoFe 2 O 4 is the material with the highest magnetocrystalline anisotropy and reasonably high magnetisation. Moreover, it shows excellent chemical and thermal stabilities, good mechanical properties 1 and is also easily synthesised with different approaches. [25][26][27][28][29] Despite its potential toxicity, cobalt-containing materials have been proposed as promising heat mediators due to their high anisotropy. 23,[30][31][32][33][34][35][36][37] The partial substitution of Co 2+ with less toxic divalent ions such as Zn 2+ (Zn x Co 1−x Fe 2 O 4 ) has been proposed in order to increase the saturation magnetisation and to lower the toxicity. 38,39 Moreover, the presence of zinc should induce the decrease of the Curie temperature opening the possibility to build auto-tuning systems. 5,33,39 The magnetic properties of Zn x Co 1−x Fe 2 O 4 (0 < x < 1) nanoparticles have been investigated over the last ten years. 5,33,[40][41][42][43][44][45][46][47][48][49][50] However, because of the strict interrelation among the different properties of the materials (size and size distribution, shape, chemical composition, capping agent), adjusting the magnetic properties exclusively on the basis of one parameter is quite rare. In particular, the systematic study of the effect of the composition on ferrite nanoparticles, leaving other parameters unchanged, is still lacking, to the best of our knowledge. Furthermore, only a few studies are devoted to the study of their heating ability. 38,39 In this framework, an ad hoc set of Zn-substituted cobalt ferrite nanoparticles with different zinc amounts but with the same particle size, particle size distribution, crystallite size and capping agent amount has been selected and studied in order to shed light on the complex magnetic properties' dependence on the properties of the material. The similarities among the samples are ideal to focus the investigation exclusively on the effect of the chemical composition on the magnetic structure (i.e., cationic distribution and spin canting), while the crystallite and particle size, polydispersity, and type and amount of capping agents remain constant. Room temperature properties have been studied and correlated to SAR values by means of a multi-technique approach that combines DC/AC magnetometry and 57 Fe Mössbauer spectroscopy, characterised by different experimental time windows (10-100 s DC magnetometry; 0.001-1 s AC susceptibility; 10 −9 -10 −7 s Mössbauer spectroscopy). 51 To the best of our knowledge, a systematic and fundamental study of the effect on the magnetic properties and heat release based on the composition (Zn x Co 1−x Fe 2 O 4 ) and cation distribution while leaving morphological and structural parameters almost unchanged is still lacking.

Synthesis
The samples were synthesised by a surfactant-assisted thermal decomposition method previously described by Sun et al. 52 Iron(III) acetylacetonate, cobalt(II) acetylacetonate, zinc(II) acetylacetonate, 1,2-hexadecanediol (10 mmol), oleic acid (6 mmol), oleylamine (6 mmol), and dibenzylether (20 mL) were added in a 250 mL three-neck round bottom flask. The system was heated to 200°C for 2 h and to 280°C for 1 h, under magnetic stirring and a blanket of nitrogen. The system was left to cool to room temperature. Then, 40 mL of absolute ethanol were added to precipitate the nanoparticles from the supernatant overnight. The as-obtained nanoparticles were separated from the supernatant by centrifugation at 4500 rpm (for 15 min), washed with ethanol and collected by centrifugation at 4500 rpm (for 15 min). The washing procedure was repeated several times and, finally, the nanoparticles were dispersed in hexane. Specific attention was paid to the control of the temperature during the synthesis in order to make it highly repeatable. To ensure it, an autotuning proportional-integral-derivative (PID) controller was used.

Characterization
Four samples of Zn x Co 1−x Fe 2 O 4 nanoparticles were synthesised as described above. The samples' chemical composition was studied by means of inductively coupled plasma-atomic emission spectrometry (ICP-AES). The analyses have been repeated two times on different portions of the samples. The samples have been labelled on the basis of the actual composition as CoFe_Znx for a Zn content (x) equal to 0, 0.30, 0.46, and 0.53, respectively.
X-Ray diffraction (XRD) patterns were collected by using two different instrument setups. The first one is a θ-θ Bragg-Brentano focalizing geometry Seifert X 3000 diffraction system equipped with a Cu Kα source (λ = 1.54056 Å), a graphite monochromator on the diffracted beam and a scintillation counter. Crystalline phases were identified by means of Analyze software. The second one is a PANalytical X'Pert PRO powder X-ray diffraction system equipped with a Co Kα source (λ = 1.78901 Å) and an X'Celerator detector. Phase identification was carried out by means of the X'Pert accompanying software program PANalytical High Score Plus. The mean crystallite size, <D XRD >, was obtained by Scherrer's equation: 53 , where, K is a constant related both to the crystallite shape and to the definition of both β and <D XRD >, λ is the wavelength of the X-rays, β is the full-width at half maximum of the peak occurring at 2θ. Here, K is assumed to be equal to 0.9, whereas β is defined by the Warren's correction, β ¼ where β exp is the experimental width of the peak and β std is the instrumental one estimated by means of the pattern of a standard obtained under the same experimental conditions. <D XRD > was calculated as a mean value by fitting, through Origin software, the most intense X-ray peaks (220), (311), (400), (422), (511) and (440) with the PseudoVoigt function, by using a 1 : 1 Gaussian : Lorentzian ratio (m u = 0.5): Refinement of the structural parameters was performed by the Rietveld method using the MAUD software 54 adopting recommended fitting procedures. 55 Structural models of the identified phases were obtained by an inorganic crystal structure database (ICSD, Karlsruhe, Germany).
The samples were analysed by transmission electron microscopy (TEM). The hexane colloidal dispersion was sonicated in an ultrasonic bath. Then, an aliquot was sampled and diluted with hexane. The diluted dispersion was then dropped on a carbon-coated copper grid and left to dry for TEM observations. The nanoparticles were observed in electron micrographs obtained with two different instruments. The first one is a TEM (JEOL JEM-1200 EX II) operating at 120 kV. The second microscope is a JEM 2010 UHR equipped with a Gatan imaging filter (GIF) with a 15 eV window and a 794 slow scan CCD camera. The mean particle size, <D TEM >, was obtained by measuring the average diameter of 350 particles using images collected in different parts of the grid. The polydispersity index, σ TEM (%), has been evaluated as the ratio between the standard deviation and the average particle size. The images were analysed by PEBBLES software 56 in a semi-automatic mode combined with a manual mode (in order to add, to the population, those particles not directly recognised by the software) and by adopting an ellipsoidal shape. The nanoparticle size distributions were fitted through Origin software using the log-normal function: Fourier transform-infrared spectroscopy (FT-IR) and thermogravimetric (TGA) analyses were carried out in order to study the capping agent and to estimate the organic phase content. FT-IR spectra were recorded in the region from 400 to 4000 cm −1 by using a Bruker Equinox 55 spectrophotometer on KBr-dispersed sample pellets. TGA curves were obtained on powders by using a Mettler-Toledo TGA/SDTA 851 in the 25-1000°C range, with a heating rate of 10°C min −1 under 50 mL min −1 argon flow.
Magnetic properties were studied by means of a Quantum Design MPMS SQUID VSM (H max = 70 kOe). Different kinds of magnetic measurements were carried out. Magnetisation vs. magnetic field curves were measured at 300 K and 5 K between −70 kOe and +70 kOe. The saturation magnetisation (M s ) was estimated by using the equation for H tending to ∞. 57 Field dependence of the remnant magnetisation was investigated through the DC demagnetisation (DCD) remanence protocol. The sample was saturated at −50 kOe and the DCD curve obtained by measuring the remanence M DCD (H) after applying reverse fields up to 50 kOe.
Calorimetric measurements of SAR were performed by means of a non-adiabatic experimental set-up built at the LAboratorio di Magnetismo Molecolare (LA.M.M) by means of a power supply CELESs MP6/400 (FIVES CELES), a watercooled heating station connected to the power supply and an induction coil. Heating curves were recorded under a magnetic field of 17 kA m −1 and 183 kHz for 300 s on water colloidal dispersions of the magnetic nanoparticles. Indeed, the hydrophobic nanoparticles have been converted to hydrophilic ones by the intercalation process with cetyltrimethylammonium bromide (CTAB, (C 16 H 33 )N(CH 3 ) 3 Br). The concentration of the colloidal dispersion was 3.5-4.0 mg ml −1 for all the samples. The temperature of the sample was monitored by an optical fiber probe (OPTOCON-FOTEMP) dipped into the solution. Samples were surrounded by polystyrene and hosted in a glass Dewar, equipped with an ethylene glycol thermostat, to ensure proper thermal isolation. The SAR values have been estimated by a linear curve fitting in the first 20 s of the heating curves (initial slope method). 57 Fe Mössbauer spectra were measured in the transmission mode with 57 Co diffuse into a Rh matrix as the source moving with constant acceleration. The spectrometer (Wissel) was calibrated by means of a standard α-Fe foil and the isomer shift was expressed with respect to this standard at 293 K. The samples were measured at 293 K and at 4 K under zero magnetic field. In order to get information on the distribution of the iron ions between the tetrahedral and octahedral sites of the spinel structure, field-spectra have been recorded at 4 K under 6 T in the perpendicular arrangement of the magnetic field vector with respect to the γ-beam. The fitting of the spectra was performed with the help of the NORMOS program using Lorentzian profiles. AC susceptibility measurements were performed on the water colloidal dispersion of the nanoparticles by a Quantum Design MPMS SQUID dedicated insert, at 5 log-spaced exciting frequencies (1-1000 Hz) and in the temperature range 10-300 K, in the absence of a static magnetic field.

Dynamic light scattering measurements
Dynamic light scattering (DLS) measurements were performed on water colloidal dispersions at the same concentration used for the calorimetric measurements by means of a Malvern Instrument Zeta Zetasizer Ver 7.03 equipped with a He-Ne laser (λ = 633 nm, max 5 mW) and operated at a scattering angle of 173°. All measurements were performed using a refraction index value of 2.42 for the material and of 1.330 for the dispersant. In all analyses, 1 mL of particle suspensions was placed in a 12 mm × 12 mm polystyrene cuvette.

Results and discussion
Composition, microstructure and morphology XRD patterns of four samples of Zn x Co 1−x Fe 2 O 4 , with the actual zinc content (x) obtained by ICP-AES equal to 0, 0.30, 0.46, and 0.53 suggest the presence of a unique spinel cubic phase (CoFe 2 O 4 , PDF card # 221086) (see Fig. 1a and ICP-AES analysis section in the ESI † for the experimental details). All the samples have a mean crystallite size of 6.2 ± 0.4 nm obtained with the X-ray copper source ( Fig. 1a and Table 1). These results have been confirmed by the use of a cobalt sourceequipped X-ray diffractometer (6.1 ± 0.4 nm) ( Table 1 and Fig. S1 †). The refinement of the structural parameters by the Rietveld method indicates that the lattice parameter linearly increases with increasing Zn content (Fig. 1c) in agreement with the Vegard rule 58 suggesting the insertion of the zinc in the spinel structure. Indeed, the lattice parameters of mixed ferrites are often a linear interpolation of the lattice constants of the extreme phases, i.e. CoFe 2 O 4 (8.38 Å) and ZnFe 2 O 4 (8.44 Å). In particular, we can better explain this trend by considering the ionic radii of the metal ions and the general preferences of the cations for a specific coordination. Co 2+ and Zn 2+ have the same ionic radius (0.82 Å 2 , 0.74 Å 59 ) whereas Fe 3+ is smaller (0.67 Å 2 , 0.65 Å 59 ). Taking into account that tetrahedral sites are smaller than the octahedral ones, a higher occupancy of the tetrahedral sites by bigger metal cations will lead to an expansion of the structure and, consequently, to an increase of the lattice parameter. This may suggest the substitution of the ferric ions by zinc ones in tetrahedral sites for the zinc-substituted samples, which is in agreement with the preference of zinc ions for the tetrahedral coordination. 2 A raw estimation of the distribution of the metallic cations in the tetrahedral and octahedral sites has been attempted (see Table S1 †). The data suggest a tendency to partially inverse structures. Moreover, the intensity ratios I(220)/I(400), I(220)/I(440) and I(422)/I(400) (Fig. S2 †) are considered to be sensitive to the cation occupancies. 46,60-62 A curve fitting confirms that these ratios linearly increase with increasing zinc content suggesting that the zinc ions occupy preferentially tetrahedral sites.
TEM images (Fig. 1b) show pseudo-spheroidal log-normaldistributed nanoparticles with a mean diameter of 7.5 ± 0.4 nm and a polydispersity index of 22 ± 2% ( Fig. 1d and Table 1). These values have been confirmed, as in the case of the crystallite size, by independent analyses with another instrument in another laboratory (Fig. S3 †). A mean diameter of 7.4 ± 0.4 nm and a polydispersity of 20 ± 2% (Table 1) have been obtained. The absence of aggregation suggests that the nanoparticles are well capped by organic molecules, as expected due to their dispersibility in organic solvents. The similar values obtained for crystallite (6 nm) and particle (7.5 nm) sizes indicate that the nanoparticles have high crystallinity. This is also confirmed by HR-TEM images showing continuous atomic lattice fringes across the particle and the absence of evidence for structural defects (Fig. 1e).
FT-IR analysis has been used mainly to verify the presence of organic molecules at the surface of the particles (capping agent) and to identify them. Indeed, the several washing steps used to purify the nanoparticles from the by-products could in principle completely remove the nanoparticle-bound organic molecules. All the spectra show as principal vibrational modes the ones typical for the hydrocarbon chain (asymmetric and symmetric CH stretching of CH 2 and CH 3 groups) and the carboxylate groups (asymmetric and symmetric COO − stretching) ( Fig. 2a and b; for further details on the complete assignments see Fig. S4 †) indicating that the nanoparticles are capped by oleate groups. 28    of a gradual substitution of cobalt ions by zinc ones within the spinel structure. TGA curves (Fig. 2c) under an inert atmosphere (50 ml min −1 of Ar) show for all the samples two weight losses occurring between 300 and 400°C and between 550 and 650°C. The first one can be associated with the decomposition of the capping molecules whereas the second one could be due to a reduction process of the inorganic core under an argon atmosphere, as suggested by other authors for oleic acid-capped zinc and cobalt ferrite nanoparticles. 64,65 In order to confirm this interpretation, a TGA curve has been recorded under an oxidative atmosphere (50 ml min −1 of O 2 ) on the CoFe_Zn0.53 sample (Fig. 2c). In this case, the decomposition of the capping agent shifts towards a lower temperature (between 200 and 300°C) but produces the same weight loss as the first weight loss of the Ar-measurement. No other weight losses have been observed, as expected for an oxidative atmosphere. Therefore, from the first weight loss in the Ar-TGA curves, the percentage of the capping agent has been calculated. Similar values have been obtained for all the samples with a mean per-centage of (15 ± 2)% (Table 1). It is easy to demonstrate that this amount of the capping agent corresponds to a monolayer of oleate molecules surrounding the surface of the nanoparticles ( Fig. 2d and S5 †). Therefore, the combined FT-IR and TGA data indicate that all the samples are composed by nanoparticles capped by a single layer of oleate molecules. Moreover, the first derivative of the thermogravimetric curve points out that the weight loss is made up by two steps in the case of the zinc-substituted samples (Fig. S5b †). This can be probably related to different interactions between the capping molecules and the inorganic cores that occur exclusively when the zinc is present in the structure. For instance, it can be related to different bond strengths between the oleate group and the different surface metal ions.
The similarities, highlighted by the careful and multi-technique characterisation, on the particle size, particle size distribution, crystallite size and capping agent amount, enable us to discuss the magnetic properties and the heating abilities exclusively as a function of the chemical composition effect caused by different zinc contents. Magnetic properties and magnetic structure M vs. H curves have been measured at 5 K (Fig. 3) and the values of the main magnetic quantities have been estimated ( Table 2). All the curves show a hysteretic behaviour without any anomalous shape associated with mixtures of hard/soft spinel phases, 66 consistently with the presence of a unique spinel cubic phase.
The hysteresis loop for the cobalt ferrite is characterised by a high coercive field, H c , (12 ± 1 kOe) near the values previously reported for nanoparticles of a similar size. 41,67-69 The coercive fields for the Zn-substituted samples were found to be equal to 4.6 ± 0.5 kOe (x = 0.30), 3.7 ± 0.4 kOe (x = 0.46) and 4.4 ± 0.4 kOe (x = 0.53). The reduced values of the coercive field for the Zn-containing samples with respect to the unsubstituted one is probably due to a decrease of the magnetocrystalline anisotropy occurring when Zn 2+ replaces Co 2+ . It is interesting to note that an increase of x between 0.3 and 0.53 does not change the coercivity significantly. The same trend has been recorded for the saturation field (H sat ), which represents the field necessary to reverse the moment of the particles with the highest anisotropy energy. H sat has been measured as the point in which the difference between the branches is under 3% of their maximum value. The saturation magnetisation (M s ) for the cobalt ferrite is about 110 emu g −1 , higher than the value usually reported for bulk CoFe 2 O 4 (80-90 emu g −1 ). 1,2 This may suggest a partially inverse structure. 70 As Zn 2+ enters the struc-ture, M s increases up to about 160 emu g −1 for a Zn content of 0.46 and decreases for a higher Zn concentration (∼140 emu g −1 ). The M s values are in agreement with the ones found for bulk Co-Zn ferrites 4,71 and higher than those reported for nanoparticles with similar compositions. 71 The differentiated remanence curve (Fig. 4), consisting of the derivative of M DCD with respect to H reverse (χ irr = dM DCD / dH reverse ), represents the irreversible component of the susceptibility. This quantity can be considered as a measure of the energy barrier distribution which, in a nanoparticle system, is associated with a distribution of particle coercivities, and it is generally called the switching field distribution (SFD). It is worth underlining that the cobalt ferrite nanoparticles exhibit a SFD centred at a higher magnetic field (15.5 kOe) with respect to the zinc-substituted ones (between 5 and 6 kOe) in agreement with the trends observed in the M vs. H curves at 5 K for the coercive and saturation fields. In addition it should be underlined that the single peak in SFD clearly indicates the presence of a unique magnetic phase, confirming the morphostructural characterization (i.e. XRD and TEM) and M vs. H measurements. It is worth noting that the presence of two magnetic phases, also exchange coupled, can be easily highlighted by DCD measurements. 7 The reduced remnant magnetisation (M r /M s ) has been found to be equal to 0.58 for the pure cobalt ferrite and lower than 0.5 for the Zn-doped samples. Although the value of 0.58 is far from that expected for the pure cubic anisotropy (the theoretical value is 0.83 72 ), it suggests that CoFe 2 O 4 nanoparticles have a mixed cubic/uniaxial anisotropy, whereas substitution by Zn 2+ leads to uniaxial anisotropic nanoparticles.
These data indicate that the insertion of Zn 2+ and the consequent decrease of Co 2+ produce at 5 K both an increase of the saturation magnetisation and a decrease of the anisotropy with respect to the non-substituted sample.
In order to interpret the magnetic behaviour, 57 Fe Mössbauer spectroscopy measurements have been carried out at 4.2 K in the absence (Fig. S6 and Table S2 †) and in the pres-  ence ( Fig. 5 and Table 3) of an intense magnetic field (6 T) in order to obtain information on the cationic distribution and spin-canting phenomena. The distribution of iron ions in the octahedral and tetrahedral sites of the spinel structure can be obtained by means of the in-field measurements carried out at low temperature ( Fig. 5 and Table 3). All the spectra can be fitted by two well separated sextets related to ferric ions located in the tetrahedral and octahedral sites of a spinel structure. The values of the isomer shift have been found to be higher than the values typical for ferric ions at room temperature (∼0. 33    . Taking into account that the exchange interactions between octahedral and tetrahedral atoms ( J AB ) are stronger than the ones between two tetrahedral atoms ( J AA ) or two octahedral ones ( J BB ), the effective field of the octahedral sub-lattice will be more affected than the effective field of the tetrahedral sublattice when the diamagnetic Zn atoms will be in tetrahedral positions. The consequence is that with the increasing of Zn substitution, the B eff values of the octahedral sub-lattice decrease while the B eff values of the tetrahedral sub-lattice are almost constant. This phenomenon corresponds well with the hypothesis that zinc occupies preferentially 4-fold coordination sites in the spinel structure. The Mössbauer measurement on the CoFe 2 O 4 provides the complete cationic distribution for this sample. The theoretical M s value for this cationic distribution can be calculated on the basis of the Néel model, which accounts for the existence of two magnetic sub-lattices in the spinel structure. This value has been found to be equal to ∼100 emu g −1 , which is similar to the experimental M s value (∼110 emu g −1 ), in agreement with the absence of spin canting. Indeed, the canting angle has been found to be around zero for all the samples (see Table S3 †). The Néel model can also be applied to calculate the theoretical M s values for the zinc-substituted samples by considering the cationic distribution obtained by the Mössbauer data and assuming that the zinc ions occupy tetrahedral sites. The theoretical values are ∼140 emu g −1 , ∼150 emu g −1 and ∼160 emu g −1 for a zinc content of 0.30, 0.46 and 0.53. The values are close to the experimental ones with the exception of the sample with the highest zinc content that probably has a part of the zinc in the octahedral sites, since no evident spin canting phenomena have been revealed by Mössbauer spectroscopy (see Table S3 †). Other authors report similar results, both for bulk and nanostructured Zn-Co ferrites, hypothesising cation distributions with a fraction of zinc ions in octahedral sites, 71 although this hypothesis is commonly excluded due to the strong preference of zinc ions for a tetrahedral coordination. 4,75 Thanks to the multi-technique approach used to characterise the sample, it is possible to reconstruct the cationic distribution on the CoFe_Zn0.53 sample by means of the experimental M s values obtained by DC magnetometry, the iron ion occupancies found by Mössbauer spectroscopy and the total content of cobalt measured by ICP-OES. The as-obtained cationic distribution is (Zn 0. 36

Magnetic properties and heating abilities
Since the increase of M s with zinc substitution (within a certain zinc content range) can be appealing for biomedical applications, we also studied the magnetic behaviour at 300 K. In this case, M vs. H curves (Fig. 6a) (Table 4). This confirms that the surface disorder phenomenon (i.e. spin-canting) is negligible in these samples in agreement with the Mössbauer results and as already observed in the literature. 77,78 A higher content of zinc produces a decrease of the saturation magnetisation at room temperature. The observed M s values for Zn-substituted cobalt ferrites are higher than the ones reported in the literature for nanoparticles 71 with similar compositions and similar to the ones found for bulk Co-Zn ferrites. 71,75 In the case of ordered magnetic systems, the thermal behaviour of the magnetization is related to the presence of low energy collective excitations (i.e. spin-waves or magnons), 79 depending on the particle size and on the chemical composition of the materials. 9,80 In the samples under investigation, the increasing of the diamagnetic ion content leads to a different temperature dependence of magnetization, justifying the different M s vs. zinc content at 300 K and 5 K.
Due to the high saturation magnetisation values, we tested the hyperthermal efficiency of all the samples by recording heating curves (Fig. 6b)   responsible for a heat release but not to the same extent. In order to quantify the amount of heat released, the SAR has been calculated by the initial slope method and values equal to 19 ± 3 and 6 ± 2 W g ox −1 , respectively, have been obtained (Table 4). A comparison with the literature data is not trivial because of the different procedures adopted to measure the heating curves, as, for instance, different frequencies and amplitudes of the external magnetic field. In this regard, in order to compare heating abilities obtained by different frequencies and amplitudes, Pankhurst et al. proposed to use a parameter named Intrinsic Loss Power (ILP) expressed in [nH m 2 kg −1 ], defined as follows: 81 ILP where SAR is expressed in W kg −1 , f in kHz and H 0 in kA m −1 .
The ILP value calculated for the CoFe_Zn0 sample is 0.36 ± 0.05 nH m 2 kg −1 . This value is comparable with those obtained by other authors for cobalt ferrite of a similar particle size 82 or similar crystallite size. 23 Higher ILP values have been reported in the literature for bigger particles, 36,82,83 for cobalt substituted iron oxide particles 23,32 or for cobalt doped magnetosome chains. 84 However, it is worth noting that the ILP is not a completely resolving tool for comparing with the literature data. Indeed, it can be applied only for superparamagnetic systems as it is based on linear response theory assumption. Furthermore, the difficulty in the comparison among the SAR values concerns other aspects that the ILP parameter does not deal with, such as differences in the experimental set-up adopted to measure the heating curve (adiabatic or non-adiabatic), in the solvents and in the analytical models adopted to determine the SAR values (initial slope method, Box-Lucas method, etc.). In addition, the comparison among literature data becomes even more difficult if we consider that often a complete chemicalphysical characterisation of the systems in terms of the crystallite size, crystallinity degree, particle size and size distribution, type and amount of capping agent, etc. is lacking.
The difference in the efficiency between CoFe_Zn0 and CoFe_Zn0.30 is considerable despite the similarities in terms of the M s , particle size, particle size distribution, crystallite size, and capping agent amount. It thus can be explained with the different anisotropies. Indeed, the magnetic measurements at 5 K have shown higher coercivity and saturation fields for the cobalt ferrite sample than the Zn-substituted one.
The results of DC magnetometry appear to be not sufficient to explain the different heating abilities of the samples. Indeed, we should consider that (i) the DC magnetometry is characterised by a time scale in the range 10-100 s; (ii) the relative extent of the measurement and relaxation time scales are critical to make the particles exhibit superparamagnetism rather than quasi-static properties. Therefore, 57 Fe Mössbauer spectroscopy and AC magnetometry have been used to study the dynamic properties of the two samples and to understand the physical reasons for the different hyperthermal efficiencies in different samples.
In order to confirm this idea, Mössbauer spectra at room temperature have been recorded ( Fig. 7 and S7 † for the spectra of all the samples). Both the spectra can be fitted by means of  (Table 5). The sextets account for the blocked spinel ferrite nanoparticles, the sharper singlet (e.g. subspectrum 1 in Fig. 7b) for superparamagnetic nanoparticles, whereas the broad singlet (subspectrum 3 in Fig. 7b) for the particles characterised by a relaxation time close to the measurement time scale. In agreement with the hypothesised scenario, the cobalt ferrite sample shows a higher area of the sextets than the Zn-substituted sample suggesting that it contains a higher percentage of blocked nanoparticles. In this view, AC magnetometry has been used to measure the temperature dependence of the in-phase (χ′) and out-ofphase (χ″) components of the magnetic susceptibility at different frequencies (1-1000 Hz) for both samples (Fig. 8). From these data, the Néel relaxation time, τ N , has been estimated at 300 K for both the samples by using the Vogel-Fulcher model 85 (for the fitting results see the Estimation of Néel relaxation time section in the ESI †): where τ 0 is the characteristic relaxation time, E b is the energy barrier against the magnetisation reversal, T is the absolute temperature and T 0 is the temperature value accounting for the strength of magnetic interactions. Values of τ N of 3.8 10 −6 s and 2.8 10 −9 s have been obtained for the CoFe_Zn0 and CoFe_Zn0.30 samples, respectively (Table 4). Besides the Néel relaxation, Brownian motion may also concur to the heat release of superparamagnetic nanoparticles. Therefore the Brownian relaxation time, τ B , was estimated using the hydrodynamic diameter obtained by dynamic light scattering (DLS) analyses. Similar distributions of the hydrodynamic diameter, centred at 29.2 and 31.5 nm, were obtained for the samples CoFe_Zn0 and CoFe_Zn0.30, resulting in τ B values of 8.5 × 10 −6 s and 1.1 × 10 −5 s, respectively. It arises that τ N and τ B for CoFe_Zn0 are of the same order of magnitude, suggesting that both the mechanisms contribute to the effective relaxation time. These estimates are nicely confirmed by the behaviour of the AC susceptibility at the melting point of the solutions (Fig. 8) where the Brownian mechanism begins. It should also be underlined that Hergt et al. 86 indicated 7 nm as the critical diameter at which τ N = τ B for cobalt ferrite nanoparticles dispersed in water, indeed a size similar to that of the sample CoFe_Zn0 (∼7.5 nm). In contrast, in the case of CoFe_Zn0.30, τ B is four orders of magnitude slower than τ N , and thus its contribution is negligible. The maximum heating efficiency will be reached when the time of the faster relaxation process matches the characteristic time of the hyperthermic measurement, τ SAR = 1/2πν = 8.7 × 10 −7 s. Therefore, cobalt ferrite, with a τ eff = (1/τ N + 1/τ B ) −1 = 2.6 × 10 −6 s, must be more efficient than the Zn-substituted sample (τ eff = of ∼2.8 × 10 −9 s).

Conclusions
Four Zn-substituted CoFe 2 O 4 samples with different zinc contents, the same crystallite size, particle size and particle size distribution, and capping agent weight percentage have been ad hoc prepared. This ideal premise has allowed the study of the magnetic properties and the heating abilities exclusively on the basis of the different chemical compositions, the other material features being unchanged. High saturation magnetisation values from 109 emu g −1 to 157 emu g −1 have been obtained at 5 K by increasing the Zn content up to 0.46. A further increase of the Zn content has resulted in a decrease of the saturation magnetisation to a value of 140 emu g −1 . The cationic distribution with the preferential 4-fold coordination of zinc ions has been found to justify the observed M s versus Zn content trend at 5 K. Two samples, CoFe 2 O 4 and Zn 0.30 Co 0.70 Fe 2 O 4 , have also shown at room temperature high and similar saturation magnetisation values (∼90 emu g −1 ). Despite the similarities between the samples, CoFe 2 O 4 has shown a triple SAR value of 19 W g −1 . This different behaviour has been justified by studying on the one hand the relaxation dynamics by combining AC magnetometry and 57 Fe Mössbauer spectroscopy and on the other hand by investigating the hydrodynamic properties of the samples. It has been found that the cobalt ferrite sample has both the Néel and Brownian relaxation times that match the characteristic time of the hyperthermic measurement whereas Zn 0.30 Co 0.70 Fe 2 O 4 has a faster effective relaxation time. This in the light of the linear response theory gives rise to the higher efficiency of the cobalt ferrite samples.
The study here presented demonstrates the power of a multitechnique approach in the comprehension of both the magnetic properties and the heating abilities. Moreover, the fundamental condition for such an in-depth comprehension remains a systematic study of the effect produced by only one parameter at a time and keeping the other parameters unchanged (in this case the chemical composition). This must be reached by the synthesis of the ad hoc prepared samples. Finally, this kind of study is also needed for the proper design of the material based on a critical choice of the features to be tuned.