4Ikerbasque, Basque Foundation for Science, 48013, Bilbao, Spain
§ Tel: +98 41 33393334; Fax: +98 41 33341244
Associate Editor: S. A. Claridge
Beilstein J. Nanotechnol. 2019, 10, 1348–1359. doi:10.3762/bjnano.10.133
Received 26 Feb 2019, Accepted 18 Jun 2019, Published 03 Jul 2019
Recent advances in the field of magnetic materials emphasize that the development of new and useful magnetic nanoparticles (NPs) requires an accurate and fundamental understanding of their collective magnetic behavior. Studies show that the magnetic properties are strongly affected by the magnetic anisotropy of NPs and by interparticle interactions that are the result of the collective magnetic behavior of NPs. Here we study these effects in more detail. For this purpose, we prepared CoxFe3−xO4 NPs, with x = 0–1 in steps of 0.2, from soft magnetic (Fe3O4) to hard magnetic (CoFe2O4) ferrite, with a significant variation of the magnetic anisotropy. The phase purity and the formation of crystalline NPs with a spinel structure were confirmed through Rietveld refinement. The effect of Co doping on structure, morphology and magnetic properties of CoxFe3−xO4 samples was investigated. In particular, we examined the interparticle interactions in the samples by δm graphs and Henkel plots that have not been reported before in literature. Finally, we studied the hyperthermia properties and observed that the heat efficiency of soft Fe3O4 is about 4 times larger than that of hard CoFe2O4 ferrite, which was attributed to the high coercive field of samples compared with the external field amplitude.
Keywords: anisotropy; cobalt; ferrite; Henkel plots; hyperthermia therapy; nanoparticles; Rietveld refinement
Technological advances in various fields have motivated the design and the fabrication of nanostructures with tuned and improved properties. Among nanostructured materials, magnetic nanoparticles (NPs) are interesting from both fundamental and technological points of view [1,2]. In recent years, ferrite nanoparticles with the general formula of MFe2O4 (M = Fe, Co, Ni, Mn) have attracted great attention of researchers due to their potential applications in biomedicine and industry . Magnetic anisotropy and interparticle interactions are important parameters that affect the magnetic properties and application fields of ferrite nanoparticles [1,4]. For example, NPs to be applied for data storage or magnetic recording must have a high coercivity, which is directly related to their magnetic anisotropy (the high coercivity keeps the recorded bits from being demagnetized) [5,6]. Magnetic interactions (e.g., exchange and dipolar interactions) have a strong effect on the magnetic behavior of a NP system (e.g., coercivity and blocking temperature) [7,8] and its potential for different applications. For example, there may be unfavorable effects in biomedical applications, such as aggregation of nanoparticles in different parts of the body . Hence, the study of this kind of interactions is of particular importance, both from a practical and a fundamental point of view. Recently, Muscas et al.  studied the magnetic behavior of mixed cobalt–nickel and pure cobalt ferrite NPs by using a random anisotropy model. Their results showed that the overall magnetic properties are the equilibrium of the interplay between the interparticle interactions and the anisotropy of the single particles. The authors of this paper believe that this study is of fundamental importance to understand the physics of nanoparticle ensembles, which, in turn, is needed to develop technological applications of these systems. Among ferrites, CoFe2O4 NPs are of considerable interest because of their moderate saturation magnetization, good chemical stability and high intrinsic magnetocrystalline anisotropy at room temperature . The anisotropy constant of CoFe2O4 (K = 2 × 105 J·m−3) is nearly one order of magnitude larger than that of Fe3O4 [11-13]. Fe3O4 NPs have been studied extensively for bio-medical applications, such as drug delivery , magnetic resonance imaging (MRI) and especially magnetic hyperthermia therapy, which is one of the efficient and new approaches for cancer treatment [4,15]. When magnetic NPs concentrated in tumor tissue are exposed to an ac magnetic field, the electromagnetic energy is converted into thermal energy, and the generated heat is used to destroy cancer cells through the elevated temperatures [16,17]. The heating efficiency of the NPs as heat sources under ac magnetic fields is often denominated as specific absorption rate (SAR), which is directly related to the area of the magnetic hysteresis loop of the nanoparticles by the following formula [18,19]:
where f is the field frequency, c is the weight concentration of the material and A is the area of the hysteresis loop. Size and shape of the particles, saturation magnetization and magnetic anisotropy, as well as field amplitude and frequency strongly affect the hyperthermia output of a NP system [15,20,21]. Sathya et al. prepared CoxFe3−xO4 nanocubes by a thermal decomposition method and showed that nanoparticles of 18–20 nm in size and a Co fraction of x = 0.5–0.7 have the highest SAR value and are suitable for hyperthermia applications . Nemati et al. prepared iron oxide nanodiscs and compared their heating efficiency with spherical NPs of similar volume at different field strengths . Their results indicated that the heating efficiency obtained for spherical nanoparticles is smaller than that measured for nanodiscs of similar volumes, especially at low field strengths. Barrera et al. prepared Co1−xZnxFe2O4 NPs and studied the dynamic energy losses of nanoparticles under an extended range of applied magnetic field strengths. They show that NPs with a larger anisotropy reveal smaller energy losses . For a more systematic study of the effect of magnetic anisotropy and magnetic interactions on properties of magnetic nanoparticles, in this work, a series of CoxFe3−xO4 (0 ≤ x ≤ 1) NPs was synthesized using a co-precipitation method. The effect of Co doping on the structural, magnetic and hyperthermia properties of CoxFe3−xO4 nanoparticles has been studied. We report a detailed study of the magnetic interactions in the samples through field-dependent measurements of remanent magnetization. In order to investigate the magnetic interactions the Henkel plot method was used, which is an effective and powerful method.
The phase purity of the samples was confirmed by X-ray diffraction (XRD) analysis. Figure 1a shows XRD patterns of CoxFe3−xO4 (0 ≤ x ≤1) nanoparticles. No secondary phases are found. The peaks intensities indicate that the samples are highly crystalline. The peaks match well with JCPDS cards (no. 01-088-0866 for Fe ferrite; no. 01-088-2152 for Fe–Co ferrites and no. 01-079-1744 for Co ferrite) indicating the formation of a cubic spinel structure with the space group Fd−3m (no. 227). Figure 1b showas that the representative (440) reflection shifts towards lower angles with cobalt ions increasingly substituting iron ions in the magnetite structure. This indicates an increase of the interplane distances (d) in the spinel structure. Similar results have been reported for NixCo1−xFe2O4 NPs by Caetano and co-workers .
According to Bragg’s law, λ = 2d·sin θ, (λ is the wavelength of X-ray wavelength, here λ = 0.154 nm, and θ is the diffraction angle), a shift of θ to lower values indicates an increasing lattice spacing d that is directly related to the lattice constant a as:
where h, k, l are the Miller indices. The values of a and and the unit-cell volume Vuc for all samples were calculated by Equation 2 and are listed in Table 1. These values are in good agreement with the values reported before for CoxFe3−xO4 (0 ≤ x ≤1) nanoparticles .
Table 1: Lattice constant (a), volume of unit cell (Vuc), crystallite size ⟨D⟩XRD and mean particle size ⟨D⟩SEM.
|parameter||x = 0.0||x = 0.2||x = 0.4||x = 0.6||x = 0.8||x = 1.0|
|a (Å) (from Equation 2)||8.36(4)||8.40(1)||8.37(1)||8.38(3)||8.37(2)||8.40(1)|
|a (Å) (from Rietveld refinement)||8.346||8.406||8.379||8.382||8.373||8.399|
|Vuc (Å3) (from Equation 2)||584(8)||593(2)||586(2)||588(6)||586(4)||593(2)|
|Vuc (Å3)( Rietveld )||581.3||594.0||588.3||588.9||587.0||592.5|
|⟨D⟩XRD (nm)||7.5 ± 0.8||8.6 ± 1.0||9.5 ± 2.0||10.3 ± 0.8||9.7 ± 1.1||13.1 ± 1.7|
|⟨D⟩SEM (nm)||40.3 ± 8.5||28.7 ± 7.4||31.5 ± 6.1||25.8 ± 6.2||24.0 ± 5.0||27.1 ± 6.5|
The radii of Co2+ ions (A-site: 58 Å, B-site: 74 Å) are slightly different than those of either Fe2+ ions (A-site: 61 Å, B-site: 78 Å) or Fe3+ ions (A-site: 49 Å, B-site: 64 Å) . Therefore, the unsystematic variations and insignificant (in the error range) difference in the lattice constant could be attributed to the change of the cation distribution in the A-and B-sites.
The XRD patterns of the samples were analyzed using the Rietveld refinement method implemented in the the “FullProf Suite” software. As an example, the Rietveld refinement pattern of the Co0.2Fe2.8O4 sample is shown in Figure 2. Table 1 shows that the parameters a and Vuc obtained from Rietveld refinement are in good agreement with those calculated with Equation 2.
The effect of Co doping on the average crystallite size was studied using the Scherrer equation:
where ⟨D⟩XRD is the average crystallite size, K ≈ 0.9 is the Scherrer constant and β is the full width at half-maximum (FWHM) of the XRD peaks. Table 1 shows that the crystallite size increases with increasing cobalt content. The increase of the crystallite size is attributed to the bond energy of Co–O (397 kJ/mol), which is smaller than that of Fe–O (407 kJ/mol) . The smaller bond energy speeds up the crystallization process, thus increasing the crystallite size in the samples.
In order to determine the particle size distribution and morphology of the samples, field-emission scanning electron microscopy (FE-SEM) was carried out. Figure 3 shows FE-SEM images of all the samples. The images reveal that particles are in the nanometer range and roughly spherical in shape. The particles size distribution in the samples was determined by measuring the size of 100 particles from the FE-SEM images fitting the size histogram with a log-normal function:
where D0 is the median diameter and σ is the dispersion. The mean diameter ⟨D⟩ = D0·exp(σ2/2) and standard deviation σD = ⟨D⟩·[exp(σ2) – 1]1/2 were determined using the fit parameters D0 and σ and are given in Table 1. The obtained values of ⟨D⟩SEM are larger than those obtained from the XRD patterns, which could be attributed to aggregation of the nanoparticles due to the presence of magnetic interactions between nanoparticles. The magnetic interactions are discussed in more detail in the following sections.
The qualitative chemical composition of the samples was investigated by using energy-dispersive X-ray spectroscopy (EDX). Figure 4 shows the EDX spectra for the samples with x = 0.2, 0.6 and 1. The EDX spectra confirm the presence of Fe, Co and O in the samples. The atomic ratio Co/Fe obtained from EDX is in a good agreement with the theoretical stoichiometry for all samples (Figure 4d).
The formation of the spinel phase and its crystal structure were verified by Fourier-transform infrared (FTIR) spectra measured at 300 K in the wave number range of 400–4000 cm−1. Figure 5 shows the FTIR spectra of the samples.
The absorption band observed at around 3385 cm−1 is attributed to the vibration mode of the O–H groups in the H2O molecules. The peak observed at around 1556 cm−1 is ascribed to amide II (NH2 deformation, N–H bending)  and the absorption peak at around 1330 cm−1 is related to the stretching vibration bands of the carboxylate group (C=O) . The latter two peaks (1556 and 1330 cm−1) are observed in all samples and can be ascribed to the presence of some impurity in the KBr pellets, which is used for FTIR analysis. Two main absorption bands are observed at frequencies below 1000 cm−1. The band around 569 cm−1 and the band around 444 cm−1 are related to the vibration of metal–oxygen (Me–O) bonds at tetrahedral and octahedral sites, respectively [28,29]. The presence of these two bands confirms the formation of the spinel structure in all the samples.
The inset in Figure 5 shows that the absorption bands related to the tetrahedral site shift towards higher frequencies with increasing cobalt content. This can be explained by considering that the Co2+ ions, being smaller than the Fe2+ ions, tend to occupy both the B-sites and the smaller A-sites (see Figure 6). This mixed occupancy in cobalt-substituted magnetite nanoparticles has been confirmed by Mössbauer spectroscopy . Therefore, it is expected that when cobalt ions substitute iron ions at the A-sites, an increasing Me–O bond distance will result. This leads to a weakening of the bond strength and a shift of the peak position towards higher frequencies.
The magnetic properties depend on the percentage of the magnetic material (ferrite) in the sample. Hence, the presence of non-magnetic impurities in the CoFe2O4 samples was checked by thermogravimetric analysis (TGA). Figure 7 shows two weight-loss stages. The first weight loss (about 2.8%), observed in the temperature range of 30–200 °C, is attributed to the vaporization of water from the sample. Since the possible decomposition of spinel ferrite is excluded because of the absence of secondary phases in the XRD patterns, the second weight loss (about 1.2%) between 200 and 500 °C can be attributed to the escape of oxygen atoms from the surface of the NPs  and the removal of impurities present in the initial raw materials.
Figure 8 shows the room-temperature magnetic hysteresis curves of the samples. The inset of Figure 8 shows that the magnetic behavior of the samples changes from soft ferrite (Fe3O4) to hard ferrite (CoFe2O4) as the cobalt content increases.
The values of saturation magnetization (Ms), remanent magnetization (Mr) and coercivity (Hc) obtained from the magnetization curves are given in Table 2.
Table 2: Saturation magnetization (MS), coercivity (HC) and remanence (Mr) of the CoxFe3−xO4 samples at room temperature.
|parameter||x = 0||x = 0.2||x = 0.4||x = 0.6||x = 0.8||x = 1|
|Ms (emu/g)||51.64 ± 0.05||53.30 ± 0.02||58.18 ± 0.01||66.26 ± 0.03||59.42 ± 0.04||70.44 ± 0.11|
|Hc (Oe)||3.36 ± 0.65||72.50 ± 0.32||206.97 ± 5.35||258.19 ± 0.99||205.36 ± 2.67||273.24 ± 5.33|
|Mr (emu/g)||0.26 ± 0.05||2.91 ± 0.03||7.04 ± 0.05||10.75 ± 0.02||8.84 ± 0.02||11.70 ± 0.07|
The difference in magnetization of the ferrites NPs is mainly attributed to the difference in particle size [32,33]. The inset of Figure 9 shows that Ms increases with increasing cobalt content due to increasing the particles size. In fact, the high surface-to-volume ratio in the smaller nanoparticles leads to an increase of the surface effects such as spin disorder and dead layer on the surface, eventually resulting in a decrease of the magnetization. The thickness of the surface dead layer (t) equaling to t = 0.56 and 0.26 nm for the samples Fe3O4 and CoFe2O4, respectively, was obtained as follows [34,35]:
where d is the particle diameter and Mb is the bulk saturation magnetization (93 and 80 emu/g for the samples Fe3O4 and CoFe2O4, respectively) [3,36].
It is expected that Co increases the magnetic anisotropy in the cubic spinel structure. The effective anisotropy constant (Keff) of particles was estimated using the law of approach to saturation (LAS), which describes the dependence of the magnetization (M) on the applied magnetic field (H) at high field strengths (H ≫ Hc). According to the LAS, the magnetization near the saturation (Ms) can be expressed as [3,21]:
where the parameter b is associated with the effective anisotropy constant as :
To calculate Keff, the experimental curves of M as a function of 1/H2 were fitted by Equation 6 at high magnetic field strengths (Figure 9). The obtained values of b and Ms were used to calculate Keff from Equation 7. The calculated values of Keff are presented in Table 3. The result for CoFe2O4 NPs is in a good agreement with that reported for NPs (ca. 10 nm) of the same compound at room temperature (Keff = 3.7 × 105 erg/cm3) . Figure 10 shows an increasing magnetic anisotropy with increasing cobalt content. This is due to the gradual occupation of the octahedral sites by cobalt ions and the stronger LS coupling originating from their strong orbital angular momentum [37,38]. The drop of anisotropy in the x = 0.8 sample might be due to the decrease of coercivity (because of the smaller size of the single-domain NPs). It is known that the cobalt ions exhibit a strong anisotropy at the octahedral sites of the cubic spinel structure . Also, Mössbauer spectroscopy showed a relatively high number of A-sites occupied by Co2+ ions in the x = 0.8 sample . Hence, another reason for the sudden drop of the magnetic anisotropy may be attributed to the increasing number of A-sites occupied with Co2+ ions, which leads to a reduced anisotropy because of the less anisotropic environment of the A-sites. Deepak et al.  observed a sharp decrease in the anisotropy for x > 0.6 in CoxFe3−xO4 nanoparticles. They attributed this effect to Co–Co interactions at high Co concentrations leading to a reduction of anisotropy, while Fe–Co interactions in the lattice increase the magnetocrystalline anisotropy. Figure 10 shows Hc as a function of the cobalt content. The increase of coercivity is mainly related to the increase of anisotropy. An increase of Hc with increasing crystallite size has been reported for single-domain NPs [13,42].
Table 3: Efective anisotropy constant (Keff), reverse field (Hr) and mean interaction field (Hint) of the CoxFe3−xO4 samples.
|parameter||x = 0.0||x = 0.2||x = 0.4||x = 0.6||x = 0.8||x = 1.0|
|Keff (× 105 erg(cm3)||3.8 ± 0.1||4.2 ± 0.1||4.5 ± 0.2||4.7 ± 0.2||4.2 ± 0.1||5.1 ± 0.2|
The analysis of remanent magnetization curves (isothermal remanent magnetization (MIRM) and direct current demagnetization (MDCD)) measured at 290 K (Figure 11) allowed us to study the mechanism of interparticle interactions. DCD and IRM curves are given in Figure 11a for all samples.
The parameter Hr shown in the inset of Figure 11a is the remanence coercivity, defined as the reverse field at MDCD = 0 . For a system of non-interacting single-domain nanoparticles with uniaxial anisotropy, MIRM and MDCD curves are related via the Wohlfarth equation [9,13]:
where mDCD (H) and mIRM (H) are normalized to the remanence saturation values MsIRM and MsDCD of the DCD and the IRM curve, respectively. The interactions in the samples can be quantitatively investigated by the Henkel plot (mDCD as a function of mIRM). According to the Wohlfarth relationship (Equation 8), the Henkel plot of non-interacting nanoparticles should yield a linear function with a slope of −2. Hence, a deviation from linear behavior indicates the presence of interactions between nanoparticles. Kelly et al. showed that in interacting systems, the Henkel plot has a deviation from linearity by an amount of δm = mDCD(H) − (1 − 2mIRM(H)) [44-46]. In particular, a negative peak (a negative deviation of the Henkel plot) in the δm curve indicates the dominance of dipole–dipole interactions, while a positive peak (a positive deviation of the Henkel plot) can be attributed to the dominance of exchange interactions. This is because the dipole–dipole interactions tend to hinder the magnetization (i.e., they have the effect of stabilizing the demagnetized state), while the exchange interactions promote a magnetization. Also, the intensity of the dip of the δm curve depends on the strength of the interactions [9,47,48]. The strength of the interactions can be estimated by calculating the mean interaction field (Hint) defined as [46,47]:
where HDCD and HIRM correspond to the peak position of the χirr (DCD) and the χirr (IRM) curve, respectively (see Figure 11). In our case, the obtained negative values for Hint (Table 3) confirmed the presence of dipole–dipole interactions in all samples. The corresponding δm curves and Henkel plots are shown in Figure. 12. The interaction field increases with increasing cobalt content, which can be related to the particle size and the larger magnetic moment of bigger nanoparticles [13,47]. The particle aggregation visible in FE-SEM images shows that the particles are interacting.
The δm plot in Figure 12 indicates that the magnetic interactions between particles are weakest in in the x = 0.6 sample. This is in good agreement with the Hint value and FE-SEM observations. The origin of the low intensity of the δm plot of the x = 0.8 sample can be attributed to the smaller particle size, which is clearly visible in the FE-SEM images (Figure 3).
In order to study the heat generation of the nanoparticles for a potential use in magnetic hyperthermia therapy, the samples were dispersed into deionized water at the same concentration (111 mg/mL) and exposed to an ac magnetic field. The increasing temperature as a function of the time was measured. Figure 13 shows a remarkable result, the temperature rise in the x = 0 sample with the smallest anisotropy and particle size was much larger than that of the other samples.
Recently, a similar behavior was reported for NixCo1−xFe2O4 by Caetano and co-workers . They showed that the heat generation increases with Ni content because of an increase of the minor hysteresis loop area. Our results can be related to the magnetic anisotropy of the NPs. When the coercivity of the NPs is larger than the field amplitude (H), the magnetization does not reach complete saturation and exhibits minor loops.
In order to characterize the SAR value of the samples, the curves of the temperature as a function of the time were fitted by the Box–Lucas model, T(t) = A·(1 – e−Bt), where A is the saturation temperature and B is a fit parameter. Here, the product A·B is the initial rate of the temperature rise. It is equivalent to the ratio dT/dt in the following equation [49,50]:
where cp is the specific heat capacity of the solution (here cp = 4.18 J/(g·K) for water), ms is the mass of the solution, ms is the mass of the nanoparticles and dT/dt is the initial slope of the heating curves. Figure 13 shows the fit curve using the Box–Lucas model (solid line) for the x = 0 sample. The SAR value, or specific loss power (SLP), was then obtained by using Equation 10.
Table 4: The values of SAR and ILP of the samples important for the hyperthermia therapy.
|parameter||x = 0||x = 0.2||x = 0.4||x = 0.6||x = 0.8||x = 1|
The SAR value is commonly used to characterize the behavior in magnetic hyperthermia. However, it is not an intrinsic property of a given system. It depends on the field amplitude and frequency. Therefore, the intrinsic loss power (ILP) parameter also is useful to compare the heating behavior measured under different values of f and H [51,52]:
The values of SAR and ILP decrease with increasing cobalt content. This is because in the samples containing cobalt (Hc > H) the system exhibits minor loops with a slight hysteresis losses.
In the present paper, we studied the effect of Co doping on the structural, magnetic and hyperthermia properties of CoxFe3−xO4 nanoparticles. The substitution of Fe by Co leads to an increase of crystallite size, saturation magnetization, coercivity and especially of the magnetic anisotropy of the nanoparticles. Interparticle interactions were disclosed by Henkel plots and δm curves. The negative deviation of the Henkel plots from linearity as well as the negative δm curves indicate a predominance of dipole–dipole interactions in all samples. It was observed that Co doping strongly reduces the specific absorption rate values (to about a fourth) in the samples, despite increasing the magnetic anisotropy, saturation magnetization and particle size. Our results showed that the heat-generation efficiency is highly impacted by the magnetic anisotropy of the nanoparticles.
CoxFe3−xO4 nanoparticles were synthesized using a facile co-precipitation method at 80 °C in air. The chemical reaction can be written as follows:
In the first step, stoichiometric amounts of the starting material (see Table 5), including FeCl3·6H2O (Merck, 99%), FeCl2·4H2O (Carlo Erba, 99%) and CoCl2·4H2O (Merck, 99%) were mixed and dissolved in 50 mL of deionized water. Also, NaOH was separately weighed and dissolved in 50 mL deionized water to a concentration of 8.0 M. In the second step, the temperature of the metal-salt solution was raised 80 °C under stirring. Then, the aqueous solution of NaOH was added quickly to the solution at 80 °C. The obtained black slurry was constantly stirred at 80 °C for 30 min. The prepared nanoparticles were washed with deionized water for several times and, finally, dried at room temperature for 24 h.
Table 5: Amounts of the starting materials to prepare 1.5 g of CoxFe3−xO4 nanoparticles.
|material||x = 0||x = 0.2||x = 0.4||x = 0.6||x = 0.8||x = 1|
|m (FeCl3·6H2O) (g)||3.502||3.493||3.484||3.475||3.465||3.456|
|m (FeCl2·4H2O) (g)||1.288||1.027||0.768||0.511||0.255||0.000|
|m (CoCl2·4H2O) (g)||0.000||0.307||0.613||0.918||1.220||1.521|
Phase purity and crystalline structure of the samples were characterized by using a Philips X′Pert Pro MPDX-ray diffractometer (XRD) with Cu Kα (λ = 0.154 nm) radiation. The diffraction patterns were analyzed using the FullProf-Suite (Version 6.0) software. Thermogravimetric analysis (TGA) was carried out in the temperature range from 30 to 650 °C with a heating rate of 10 °C/min under N2 flow using a TGA/SDTA 851 Mettler Toledo thermogravimetric analyzer. Fourier transform infrared (FTIR) spectra of the samples were obtained in the range of 400–4000 cm−1 by pressing the powders in KBr pellets. The morphology and elemental chemical composition of the samples were investigated using a Tescan Mira 3 field-emission scanning electron microscope (FE-SEM) equipped with an energy-dispersive X-ray spectrometer. The magnetic properties were studied at room temperature by a custom-built vibrating sample magnetometer (VSM) with a maximum applied field of 18 kOe. The field-dependence of remanent magnetization was measured by following the isothermal remanent magnetization (MIRM) and direct current demagnetization (MDCD) protocols. For the MIRM measurement, an external field was applied to a demagnetized sample, then it was switched off and the remanent magnetization was measured. This process was repeated, increasing the field up to 18 kOe. In the MDCD measurement, the samples were magnetized at −18 kOe. After that, a small field in the opposite direction of magnetization was applied, then the field was switched off and the remanence MDCD was measured. This process was repeated increasing the field strength up to +18 kOe. Magnetic hyperthermia properties were studied by using a custom-built setup at a frequency of 92 kHz and a field amplitude of 27 Oe.