Magnetic Hyperthermia in the 400–1,100 kHz Frequency Range Using MIONs of Condensed Colloidal Nanocrystal Clusters

In the current study, we explored the magnetic hyperthermia performance of condensed–clustered magnetic iron oxide nanoparticles (MIONs) in the range of 400 kHz to 1.1 MHz at low field amplitudes. The strong interparticle interactions, present in such systems, can influence the hyperthermia power produced by MIONs. Herein, the heat dependence, as an increasing function of frequency, with a fixed magnetic field strength of 3 mT is recorded, revealing a direct relationship between the two physical quantities and a high heating efficiency for the condensed–clustered MIONs. In particular, the specific loss power (SLP) (or specific absorption rate [SAR]) parameter, which is the ratio of the heat power in watts produced per nanoparticle mass in grams, is linear to a good degree to the oscillating frequency with a step of roughly 30 W/g per 100 kHz increase. In addition, all the measurements were within the safety limits proposed by Hergt and Dutz criterion of H f ≤ 5 × 109 A/ms for clinical application of magnetic fluid hyperthermia (MFH). Finally, the measured data of temperature vs. time at each frequency were interpreted in terms of simple thermodynamic arguments, thus extracting useful thermodynamic parameters for the heat power generated by the condensed–clustered MIONs.


INTRODUCTION
The application of nanotechnology in the field of medicine promises to revolutionize the treatment of malignant diseases such as cancer. This has led to the development of multifunctional nano-assemblies termed as theranostics (Xie et al., 2010;Hu et al., 2020;Misra et al., 2020;Shakeri et al., 2020) for their combined diagnostic and therapeutic action. Hybrid magnetic nanostructures (in particular, iron oxides) are among the most interesting classes of nano-theranostics and have been the topic of extensive research interest over the past decade (Yoo et al., 2011;Revia and Zhang, 2016;Lorkowski et al., 2021). Their excellent biocompatibility and inherent multifunctional nature of their magnetic core mark them suitable for a diverse set of biomedical applications, mainly including MRI (Zhou et al., 2014), physical tumor targeting through magnetic manipulation (Feng et al., 2017), and magnetic fluid hyperthermia (MHF) (Kozissnik et al., 2013;Deatsch and Evans, 2014;Das et al., 2019). The magnetic properties of magnetic iron oxide nanoparticles (MIONs) are strongly dependent on their physical characteristics. When the size of magnetic nanoparticles (MNPs) decreases below a certain limit (∼20 nm for MIONs), then the nanoparticles exhibit a single magnetic domain and the total magnetic moment of all unpaired electrons will have a preferred orientation (Krishnan et al., 2006). Further reduction in the size of nanoparticles results in a reduction in the energy barrier required to rotate this magnetic moment, until ultimately the thermal energy is sufficient to force the moment to random directions, a phenomenon known as "super-paramagnetism." In the absence of an external magnetic field, superparamagnetic particles do not exhibit magnetization phenomena, a property critical for their colloidal stability. However, under the application of an external magnetic field, the magnetic moments of the nanoparticles align along the direction of the field and the saturation magnetization can be achieved at relatively low fields. Another important physical parameter in work related to MNPs is the relaxation time τ which is an inherent characteristic time of the MNPs as a response to the alternating magnetic field (AMF), describing different rotational modes. Two such modes are usually present, namely the Néel rotation according to which the magnetization of the MNPs is rotating with the particle remaining fixed and the Brownian rotation where the magnetization of the MNPs is locked at some crystal axis and thus in an effort to follow the field, it forces the particle to rotate as well. MNPs have the potential to generate heat through the application of external AMF, property that can be exploited for in-vivo induction of hyperthermia (Fortin et al., 2007;Ortega and Pankhurst, 2012) aiming to the destruction of the cancerous tumors. The heating performance of different MNPs for the applications of magnetic fluid hyperthermia (MFH) is often quantified by the "specific loss power" (SLP) parameter (also known as "specific absorption rate" [SAR]), which is defined as the power diffused by a MNP per unit mass (Rosensweig, 2002). The unit of SLP is mW/g. Usually, the MNPs are suspended in aqueous solutions (ferrofluids) and this important property is determined experimentally by the simple formula (Sudame et al., 2020): where c is the ferrofluid-specific heat (of the carrier fluid in J/g · K), µ = m n /m is the mass fraction of the MNPs mass m n to the total solution mass m, T is the temperature rise of the ferrofluid, and t is the duration of the heating. As the T(t) curve is not linear but has an exponential rise-saturation behavior, the quantity T/ t is implied as the largest slope in Equation 1 which happens at the initial time t = 0 of the application of the AC field. The physical meaning of SLP is that the higher it is, the more heat is produced by the MNPs. Different experimental setups have been developed for the assessment of the heating efficiency of MNPs (Andreu and Natividad, 2013). These mainly include calorimetry experiments conducted in an adiabatic or non-adiabatic way (Rajan and Sahu, 2020;Boekelheide et al., 2021) as well as AC magnetometry experiments, where the dynamic magnetic hysteresis loops M(H) are measured with a pickup coil technology (Connord et al., 2014;Garaio et al., 2014;Salas et al., 2014). Nevertheless, the widely different experimental conditions (i.e., oscillating frequency, AMF amplitude, sample volume, etc.) employed by different groups, have led to the introduction of another parameter called "intrinsic loss power" (ILP) which is defined as where H is the magnitude of the applied external field and f is the excitation frequency. ILP has been proposed as a "normalized SLP parameter" as the magnetic excitation power is proportional to H 2 and the frequency f . ILP aims to be an independent parameter of AMF strength and frequency, eliminating the effect of different experimental conditions, thus enabling a direct comparison of the heating efficiency of MNPs (Kallumadil et al., 2009). However, the advantages of the concept of ILP are still under debate as SLP rarely appears as quadratic in field intensity (Périgo et al., 2015). In addition, the resonance phenomena due to the aforementioned relaxation times can also change this linear dependence of the measured signal to f . The theoretical models (Rosensweig, 2002) have also been employed for the prediction of power dissipation from the MNPs. These models are based on the intrinsic properties of MNPs, with the most commonly used being the linear response theory (LRT) model by Rosenweig (Rosensweig, 2002). In this study, it can be seen that the power loss [ (Rosensweig, 2002), Equation 14 therein] has a frequency dependence on the form f 2 /(1 + 2πf τ 2 ) which is quadratic at low f , leads to saturation at high f , and has roughly a linear increase at values somewhat lower than 2πf τ = 1. Nevertheless, this model refers to non-interacting particles, neglecting the interparticle interactions developed in concentrated samples (De La Presa et al., 2012), multicore systems (Yoon et al., 2011;Serantes et al., 2014), or condensed-clustered nanoparticles (Puntes et al., 2004;Zoppellaro et al., 2014) which can strongly affect the heating performance of MNPs (Deatsch and Evans, 2014;Salas et al., 2014). Condensed clustering of MIONs refers to a specific structural motif of colloidal nanocrystal clusters (CNCs) in which individual magnetic nanocrystallites are densely packed in such a fashion that their crystal planes adopt the same crystallographic orientation through the epitaxial aggregation during their crystal growth (Zoppellaro et al., 2014;Coral et al., 2018). This can result in large crystallites of ∼50-300 nm size to show a singlecrystal-like structure throughout the large areas of their volume and yet behave like the ensembles of smaller independent (5-15 nm) nanocrystallites (i.e., superparamagnetic behavior). This internal organization of MNPs can significantly enhance their magnetic properties due to the collective phenomena exhibiting high magnetic moment, high contrast with MRI applications (Xu et al., 2012), and improved hyperthermic performance (Lartigue et al., 2012). In these systems, the developed dipolar interparticle interactions are particularly strong due to the very close proximity of the magnetic nanocrystals. Such strong interactions may lead to increased magnetic anisotropy in the system which in turn could result in enhanced heating efficiency. Nonetheless, this is not always the case as more complex phenomena arising in such systems [e.g., surface spin disordering due to defected nanoparticle surface coordination and spin glass features (Kostopoulou et al., 2014)] could either favor or hinder the heating performance of clustered systems (Landi, 2014;Blanco-Andujar et al., 2015). In this context, the experimental evaluation of hyperthermia agents remains the most valid route for the extraction of SLP and related hyperthermia parameters, especially for strongly interacting systems with a high interest for the applications of magnetic hyperthermia, such as the condensedclustered MIONs.
This article aims to study the effects of frequency on the hyperthermia power, for a condensed-clustered MION-based system developed by Zoppellaro et al. (2014). The frequency range was 400 kHz to 1.1 MHz, (a) which covers almost all of the width of the fixed frequencies reported previously in the literature [see Table 1 that summarizes a quick outline of the parameters used in references (Bhayani et al., 2013;Hallali et al., 2016;Kumar et al., 2018;Bhardwaj et al., 2020;Brizi et al., 2020;Shaterabadi et al., 2020;Sudame et al., 2020;Khan et al., 2021;Rousseau et al., 2021;Yamaminami et al., 2021)] and (b) it is in agreement with the medically accepted range, as described in Khan et al. (2021), quoted here: "The choice of frequency must be appropriate for the site location of magnetic Hyperthermia applications in clinical trials like, preferably > 50 kHz to avoid skeletal muscles' stimulations and lower than 2 MHz (2, 000 kHz) for penetration depth inside the tissues." In addition, other reports state that frequencies <10 MHz are not attenuated by tissue and thus penetrate deep into the body (Attaluri et al., 2020). Furthermore, for in-vivo applications, safety limits regarding the strength of the magnetic field should also be considered to ensure that patients will not experience any discomfort during therapy sessions, caused by non-specific tissue heating due to eddy currents. To ensure that, a widespread limit known as "Brezovich criterion" (Atkinson et al., 1984) has been introduced, where the product of the applied magnetic field and oscillating frequency should be H f < 4.85 × 10 8 A/ms. Nevertheless, it is estimated that for exposure of adequate small body regions this strict limit can be exceeded up to 10 times, thus becoming H f ≤ 5 × 10 9 A/ms, as proposed by Hergt and Dutz (Hergt and Dutz, 2007). In the present study, all experiments were performed fulfilling the above Hf condition with the highest frequency tested to be Hf = 2.58 × 10 9 A/ms and therefore suitable for biomedical applications. Even though Néel relaxation times could result in resonance frequencies (Yamaminami et al., 2021) within our range, no such resonance was observed in our case and the measured temperature was proportional to the applied frequency. This could be explained by the increased interparticle interactions of the condensed magnetic clusters, as discussed later.

Synthesis of Condensed-Clustered MIONs
For the synthesis of the MIONs, we followed the protocol presented in our earlier study (Zoppellaro et al., 2014)

Physicochemical Characterization of MIONs
The morphology of the synthesized nanoparticles was investigated by transmission electron microscopy (TEM), wherein samples were prepared by casting a droplet of a dilute aqueous suspension of nanoparticles (0.01% w/v in Fe2O3) on copper grids coated by a Formvar carbon film. The images were obtained by a JEOL, JEM-2100 instrument operating at 200 kV. The determination of the hydrodynamic diameter (D h ) of nanoparticles dispersed in deionized H 2 O was performed with a ZetaSizer Nano series Nano-ZS (Malvern Instruments Ltd., Malvern, UK) equipped with a He-Ne laser beam at a wavelength of 633 nm and a fixed backscattering angle of 173 • . The concentration of the measured colloids was 0.0125% w/v (g/100 mL) in Fe 2 O 3 . The zeta potential of the nanoparticles was assessed with the same instrument as the average of 100  Hyperthermia Experimental Setup The Resonator Unit was connected to the RLC circuit which was formed by the resistance and the inductance of the excitation coil together with the capacitance of the capacitor. The circuit was brought to resonance to get the maximum current to the coil and thus a maximum magnetic field. To achieve a broad range of resonance frequencies, the matching capacitor had the following preset values as summarized in Table 2, which when they were combined with the inductance 3.27 ± 0.04 µH of the excitation coil, resulted in the resonance frequencies in the range of 400 to 1075 kHz. The excitation coil of length 3.9 cm was composed of 8 turns each with a diameter of 4.2 cm and it was made of copper wire of 1.05 mm diameter integrated with a cooling fin in order to be easily air-cooled by a nearby fan. In all experiments, the voltage of the resonator was adjusted so as to have a constant magnetic field of 3.0 ± 0.1 mT, which was measured indirectly by the induced voltage of a single-turn detection coil, since the high frequencies produce enough high voltages to be recorded on the oscilloscope.
Two different alcohol thermometers were used for the temperature measurements, one of which was directly immersed in the sample solution to record hyperthermia phenomena and the other to record the air temperature inside the detection coil, in order to ensure that temperature rises in the first coil were emerging from the solution itself and not the coil selfheating. The presence of high-frequency magnetic fields inside the detection coil makes it impossible to use conventional thermometers such as thermocouples or resistance temperature detectors, since they contain metals that will develop induced voltages with corresponding reading errors. Thus, the traditional alcohol thermometers were used and their values were read with the help of a magnifier in order to get more precise readings.

Experimental Procedure
The aim of this study was to record temperature rises in the MNP solutions at different frequencies and at different exposure times, in order to observe the effect of the exciting field frequency on the hyperthermia phenomenon. In order to achieve this task, we followed the steps below: Initial steps • The sample, which has the form of a solution in a sealed glass vial, was removed from the refrigerator where it was kept at a low temperature of 5 • C so as to be chemically inert. • The sealed vial was placed in an ultrasonic bath to homogenize and eliminate the agglomerations which were sometimes created during measurements. Magnetic particles tended to accumulate together due to their mutual magnetic attraction to each other. • Following the ultrasonication, the vial was placed at the center of the excitation coil and one of the thermometers was brought into contact with the solution through a snag hole on the vial lid. • The second thermometer was held by a clamp so as to have its measuring tip inside the excitation coil but not in contact with the glass vial.

Repeating steps
• The matching capacitance was set to an appropriate value so as to achieve the desired frequency. • Enough time was given to the sample (30 min) with no AC field present, to be brought to a thermal equilibrium with its environment. At this point, its initial temperature θ 0 was recorded. • The resonator was then turned on and its output voltage was adjusted so as to have a value of an AC magnetic field equal to 3.0 ± 0.1 mT, according to the oscilloscope reading. This happened when the timer was set to zero. • The temperatures θ and θ air of the two thermometers (sample and air correspondingly) were recorded every minute for a total period of 10 min. • The resonator was then turned off at t = 10 min and the vial was shaken by hand, without removing the thermometer, to avoid agglomeration and to have a homogeneous solution without precipitations. • The last five steps were repeated with a new frequency.

RESULTS AND DISCUSSION
The MNPs used herein (coded as MagAlg) have been studied and characterized in a previous study, where their physicochemical and magnetic properties are described in detail (Zoppellaro et al., 2014). The basic physicochemical characterization was conducted for the present study, on new and freshly synthesized MagAlg samples. As shown in Figures 2a,d, the TEM micrographs revealed the morphology of the densely packed MNPs, in which the individual magnetic nanocrystals, having dimensions in the range of 8-13 nm, aggregate together in a way to form larger clusters. Figure 2b shows the selected area electron diffraction (SAED) pattern recorded for the MagAlg assembly.  (Figures 2c,e) show the appearance of crystalline fringes in the condensed cluster nanoparticles (Figure 2c, the arrows indicate the spacing corresponding to the [311] plane, d = 2.52 Å). Furthermore, the large areas of the clustered magnetic material form crystallographically aligned nanocrystals. The results from Dynamic Light Scattering (DLS) and electrokinetic measurements are given in Figures 2f,g. The MagAlg nanoparticles exhibit an average size (D h ) of 93 nm and the zeta potential value of −40 mV, confirming the absence of much larger aggregates when dispersed in deionized H 2 O and ensuring colloidal stability of this system during the magnetic hyperthermia measurements. The magnetic behavior of MagAlg, in powder form, at 300 K is given in Figure 2h, which describes the variation of the magnetic susceptibility (M s ) vs. strength of the applied external field (B). The M s value reaches 71.11 Am 2 kg −1 at 7.5 T, fully consistent with the previous results (Zoppellaro et al., 2014). At a much lower field (B = 0.3 T), the M s value remains large (63.06 Am 2 kg −1 ) without the appearance of a coercive field. Taken together, these results translate into the possible implementation of the MagAlg system in hyperthermia treatments. It is recognized, however, that thorough in-vitro and in-vivo studies are required to establish the biocompatibility/safety of the MagAlg nanoparticles for any potential biomedical applications.
Concerning the calorimetric measurements of this study, Figure 3 shows the θ 10 = θ (10) − θ (0) data, where θ (10) is the final and θ (0) is the initial temperature for the 10-min interval, for the range of frequencies of 400 to 1075 kHz that our resonator was able to cover. The strength of the magnetic field was kept fixed at 3 mT (2.4 kA/m). The top data correspond to the MIONs solution temperature and the bottom to the air temperature inside the coil. It is noted from this graph that θ 10 for MNPs increases roughly linearly with the frequency, while the air temperature is practically constant. In fact, the θ 10 data are fitted better to a second-degree polynomial with an intercept at zero frequency equal to 1.4 • C which is within the error range in the graph so it can be assumed as zero, as expected for a zero frequency field (DC field). Upon judging the air data, it is concluded that the high-frequency AC current does not heat up the coil and subsequently the air inside it (remember that the coil is air-cooled externally by the help of its fins) and that any significant rise in the temperature in the graph is entirely due to the hyperthermia of the MIONs solution. Similar experiments performed in our lab with water replacing the MIONs solution (data not shown) confirm this conclusion as there were no temperature changes recorded by both thermometers.
To examine the time dependence of the hyperthermia effect, we plotted in Figure 4 the θ = θ (t) − θ (0) data vs. time for different frequencies in the range of 400 to 1075 kHz. It is obvious from this graph that θ increases with both time and frequency. In particular, θ seems to achieve saturation during the time interval of 10 min in our experiment. This is typical with heating experiments in which the system always reaches a timestanding condition where the temperature no longer changes with time. It will be shown in the next subsection that the θ (t) curve has an exponential dependence from which the important parameters of the solution, such as the heat power and the saturation temperature, can be extracted.

Analysis
Let the MNP sample solution have mass m, specific heat c, and a temperature θ which is assumed to be a function of time t. According to Newton's law of cooling, the solution loses heat at a rate of hA (θ − θ 0 ), where h is the heat transfer coefficient, A the surface area through which the transfer takes place, and θ 0 is the surrounding temperature. Additionally, let a heat source (the hyperthermia) supplying the solution with a heat rateQ H (heat per time). Even though this is an internal heat source, it can be assumed as an external source in order to better understand the thermal physics of the system. In this fashion, let the solution exchange an amount dQ of total heat with its surrounding (cooling plus warming) within a time interval dt. Then its heat rate dQ dt will be equal to the sum of the above two rates, one being negative as it describes heat loss and the other being positive as it describes heat gain: From the definition of the specific heat c, we have dQ = mcdθ, where dθ is an incremental change of θ upon an incremental exchange of heat dQ. Substituting this expression in Equation 3, it results in a first-order linear differential equation (DE) on θ (t):  This DE can be easily solved for θ (t) = θ (t) − θ 0 as where and k = hA/mc (7)  in the data of Figure 3. Equation 5 has a classic exponential increase-saturation behavior, like the curves formed by the data shown in Figure 4 and it is an excellent fit for them [the corresponding fit function in Mathematics is known as the empirical Box-Lucas mode, also used in a study by Kumar et al. (2018)]. The fits are shown in Figure 5, where only selected curves for specific frequencies are plotted for clarity. The results for the fitted parameters θ s and k are shown in Table 3 and a plot of k vs. the excitation frequency is presented in Figure 6. From Equation 7, we expected k to be a constant, independent of frequency. Indeed, it seems from Figure 6 that either this parameter is constant with a relatively large error, or it is partially constant, exhibiting a steplike behavior, with the lower frequency values being accumulated around the value k 1 = 0.20 min −1 while the higher frequency values being around k 2 = 0.25 min −1 . A possible explanation for this behavior is perhaps the loss of mass that we observed in our experiments, which is shown schematically in Figure 7. Even though our glass vial was sealed during the whole experiment with only a tight hole on its lid for the thermometer to be placed snag through it, we observed that the solution level dropped from the initial height of 11.6 mm to a final height of 9.8 mm. The origin of this loss is not understood, which could be possibly due to a partial evaporation. But even if this is true, it looks strange to have happened only at the narrow frequency range of 620-650 kHz. Anyhow, whatever the source of the mass loss is, it is quite clear from Equation 7 that it directly affects the value of k because it is inversely proportional to the mass m. In fact, the proportion 11.6/9.8 = 1.18 of the aforementioned levels is quite close to the ratio k 2 /k 1 = 1.25 of the two-step k values in Figure 6. In the rest of this study, it will be assumed that there are two different solution masses, the initial mass m 1 corresponding to k 1 and the final mass m 2 corresponding to k 2 . As our aqueous solution of 2.5 mg/mL is quite dilute, it can be safely assumed that its density is close to the water density of 1.0 g/mL and from it, we can calculate the two masses using the solution volumes as m 1 = 4.4 g and m 2 = 3.7 g. Equations 6 and 7 can be combined to get the hyperthermia heat produced per unit timeQ H as follows: Note that from Equation 5 it can be easily seen that the initial slope is equal to k θ s . Thus, Equation 8 is in agreement with Equation 1 above. The only thing missing to convertQ H to SLP is to convert the total solution mass m to the mass fraction µ = m n /m. This is easily done by using the water density ρ = 1.0 g/mL by denoting µ = m n /ρv = x/ρ where x is our MNP concentration of 2.5 mg/mL and v is the solution volume. Thus  For the specific heat c, the assumption is similar to that we did for the mass calculation, i.e., the dilute solution basically behaves like pure water, and we can take as c the value of the water c = 4.18 J/gK. Our SLP values were calculated with the help of Equation 9 and Table 3 at different frequencies and are shown in Figure 8 in Watts/g units (the minutes in our calculations have been converted to seconds). We concluded that the data of Figure 3 fit better to a quadratic function of the frequency rather than a linear one. However, the difference is rather small and for the purpose of simplicity and calibration, it can be assumed that the data points of Figure 8 are approximately lying on a straight line (the R 2 value is close to 0.91, indicating a fair linear fit), and get for the slope a value of 0.31 W/g · kHz, i.e., an SLP increase of 31 W/g per 100 kHz. Additionally, a power law fit leads to a power exponent of 1.1, and it is not a better fit (R 2 = 0.85) compared with the linear one. This observed linearity of SLP with frequency is in agreement with the previous studies of the clustered systems at low AMF field strength (Shubitidze et al., 2015) and could be possibly attributed to the inherent strong interparticle interactions in such ensembles. Such strong dipolar interparticle interactions can result in a coupling of Brownian and Néel processes leading to an increasing "effective" relaxation time of the system, as it has been proposed by computational studies (Ilg and Kröger, 2020), thus shifting the resonance of the condensed-clustered systems to lower frequencies (Shubitidze et al., 2015). As summarized in Table 1, the obtained SLP values are within the range of values found in the literature, wherever they are reported. In particular, our SLP increase shown in Figure 8 of about 200W/g in the range of 400 -1,000 kHz is in agreement with Garaio et al. (2014), Table 1 therein for the low field regime. However, in the study by Garaio et al., the SLP values were derived at higher magnetic fields (∼4 times higher) and with more concentrated samples (an order of magnitude higher), highlighting the superior performance of condensedclustered MIONs employed for the present study compared with non-clustered systems.

CONCLUSIONS
In the current study, we evaluated the hyperthermia efficiency of condensed-clustered MIONs in a wide frequency range. According to this study, there is a direct relationship the power produced by the condensed-clustered MIONs on. f In particular, the SLP parameter, which is the ratio of the heat power in Watts produced per MNP grams, is linear to a good degree to the frequency with an increase of roughly 30W/g per 100 kHz. This linearity can be interpreted in the frame of the inherent strong interparticle interactions in such ensembles, resulting in the coupling of the Brownian and Néel processes, thus shifting the resonance of the condensed-clustered MIONs to lower frequencies away from the ones usually employed for in-vivo magnetic hyperthermia. The recorded SLP values were in the 100 to 300 W/g range of which is in agreement with the previous reports. Nevertheless, the lower magnetic field employed in the present study highlights the enhanced heating efficiency of condensed-clustered MIONs. The produced heat rateQ H was extracted by a simple thermodynamic analysis where the fit parameters are related to different physical quantities. The temperature vs. time curves show a simple exponential rise-saturation behavior. For the majority of studies reported in the literature, the AC frequency f is usually fixed, as it is hard to produce high magnetic fields at these frequencies.
In the present study, we were able to use a 3mT field (2.4kA/m) over a frequency scan of 400 kHz to 1.1 MHz being always within the safety limits as proposed by the Hergt and Dutz criterion of H.f ≤ 5 × 10 9 A/ms for the clinical application of magnetic hyperthermia. Evaluating the heating performance of condensed-clustered MIONs at higher frequencies and low magnetic field amplitudes can provide invaluable information for their potential use as magnetic hyperthermic agents in cases where low magnetic fields or low nanoparticle doses are required for safety reasons.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
DK: experiment design and theoretical model. GS: experiment design and hyperthermia measurements. AK-N: synthesis of nanoparticles and hyperthermia measurements. GZ: synthesis of nanoparticles and characterization of nanoparticles. KS: hyperthermia measurements. All authors contributed to the article and approved the submitted version.