Engineering Gold Shelled Nanomagnets for Pre-Setting the Operating Temperature for Magnetic Hyperthermia

This study investigated the fabrication of spherical gold shelled maghemite nanoparticles for use in magnetic hyperthermia (MHT) assays. A maghemite core (14 ± 3 nm) was used to fabricate two samples with different gold thicknesses, which presented gold (g)/maghemite (m) content ratios of 0.0376 and 0.0752. The samples were tested in MHT assays (temperature versus time) with varying frequencies (100–650 kHz) and field amplitudes (9–25 mT). The asymptotic temperatures (T∞) of the aqueous suspensions (40 mg Fe/mL) were found to be in the range of 59–77 °C (naked maghemite), 44–58 °C (g/m=0.0376) and 33–51 °C (g/m=0.0752). The MHT data revealed that T∞ could be successful controlled using the gold thickness and cover the range for cell apoptosis, thereby providing a new strategy for the safe use of MHT in practice. The highest SAR (specific absorption rate) value was achieved (75 kW/kg) using the thinner gold shell layer (334 kHz, 17 mT) and was roughly twenty times bigger than the best SAR value that has been reported for similar structures. Moreover, the time that was required to achieve T∞ could be modeled by changing the thermal conductivity of the shell layer and/or the shape/size of the structure. The MHT assays were pioneeringly modeled using a derived equation that was analytically identical to the Box–Lucas method (which was reported as phenomenological).


Introduction
The principles that underlie magnetic hyperthermia (MHT) using magnetic nanoparticles (MNPs) have been extensively discussed over recent years [1][2][3]. In short, heat is dissipated mainly from the relaxation losses of the MNPs, which creates localized heat in the surrounding environment of the dispersed nanomagnets [4].
In cancer treatment, controlled localized heating with the temperature of the surrounding biological environment stabilized in the range of 41-48 • C represents a key strategy, which is challenging to achieve in practice but is nevertheless required to promote tumor cell death by apoptosis instead of necrosis. When the latter is triggered at higher temperatures, it has negative impacts on both the health of the patient and that of the local normal tissue [5]. Therefore, the use of this strategy within the medical field requires the surrounding environment to be at a safe temperature in order to promote apoptosis at the tumor site and avoid necrosis or any other adverse effects in the surrounding normal tissue, such as inflammation processes that are caused by untargeted overheating [6].
However, the lack of robust and accurate (preferably noninvasive) temperature probing methods for targeted biological sites has been a limiting factor in the application of MHT technology in practice [7]. Moreover, numerous MHT protocols have been employed in clinically unsafe conditions, which aimed to achieve therapeutic temperatures with minimal doses of the administrated MNPs [8].
Therefore, in MHT protocols, the strategy of pre-setting the asymptotic operating temperature (T ∞ ) for the surrounding environment via the modulation of the morphological and magnetic properties of the MNPs, as well as the alternating magnetic field characteristics, is a very promising direction to explore. This approach could lead to important breakthroughs in the development of MHT technology from its current restricted status within clinical study to widespread clinical application. Moreover, regardless of the MHT protocol conditions, safe operation clearly requires the use of optimized biocompatible nanomagnets to deliver controlled heating outputs under low alternating magnetic field amplitudes [4,9]. Regarding heating control, different approaches have been reported in the literature, including the modulation of MNP concentrations within liquid media with varying viscosities [10], dipolar interaction via MNP clustering [11], the nonmagnetic metal doping of typical MNPs [12] and the use of core-shell magnetic nanostructures [13].
Surface dressed iron oxide NPs, such as maghemite (γ-Fe 2 O 3 ), have been extensively employed for MHT purposes owing to their high colloidal stability when dispersed in aqueous media, superior biocompatibility and relatively well-known magnetic behavior for heat generation [14,15]. Additionally, metal shelled MNPs that use nonmagnetic metal shells (such as gold) to dress magnetic cores can improve both the magnetostatic and magnetodynamic properties of shelled nanomagnets [16]. Moreover, the metal shell layers can substantially improve the heat transfer between the magnetic iron oxide cores and their surrounding environment. Additionally, there is evidence that gold shells improve the biocompatibility of magnetic iron oxide nanoparticles [17][18][19]. Likewise, the metallic shells of gold shelled MNPs allow for improved photothermal responses [20], improved X-ray contrast in computer tomography (CT) imaging [21] and improved magnetic resonance imaging (MRI) contrast characteristics [22], which can be ultimately combined with magnetically responsive cores (MHT, CT and MRI) for applications in theranostics [21,23]. It has been reported that the MHT performance (i.e., steady state condition (SSC) with short transient times and asymptotic temperatures in the surrounding environment) of gold shelled MNPs decreases as the thickness of the gold shell increases, which has been claimed to be due to the interaction at the shell-core interface and the loss of cooperative behavior between neighboring MNPs [23]. Interestingly, as strongly suggested in the present report, the fine tuning of the shape and size of the magnetic cores and the thickness of the metallic shells and their composition seems to provide a way to modulate T ∞ in MHT protocols, i.e., the typical time that is needed for the system (a host medium that incorporates the metal shelled MNPs) to achieve T ∞ and the value of T ∞ itself. The impact of the size and shape of MNPs on MHT performance has been explored in recent years [24,25].
The present study investigated the successful modulation of T ∞ in MHT protocols via the modulation of the thickness of the gold shells that were deposited onto maghemite NPs (with a fixed particle size and particle concentration) under different external conditions for the AC (alternating current) applied magnetic fields. Asymptotic temperatures that were in the range that is required for MHT to trigger cellular apoptosis were achieved for magnetic field frequencies and field amplitudes that were in the range of 100-650 kHz and 9-25 mT, respectively. Moreover, the time (t) dependence of the system's temperature (T) was recorded up to the SSC and curve-fitting using a mathematical function that was developed in this study, which allowed for the extraction of τ (typical heat transfer time) and T ∞ from the achieved SSC. Importantly, the successful modulation of the SSC in the MHT protocols represented robust proof that the pre-setting of certain parameters for application in cellular apoptosis was actually possible. Additionally, the presented mathematical derivation of the function that was used to model the temperature increases during MHT was pioneering and extended the similar phenomenological mathematical expression that has been long employed in the literature, i.e., the Box-Lucas method [26][27][28][29][30]. At this stage, it is worth mentioning that the Box-Lucas method has also been used to describe temperature increases that are due to photothermal effects [31,32]. An important metric for the application of MHT technology is the so-called specific absorption rate (SAR), which is defined as the thermal power that is dissipated per unit mass (kW/kg). A key piece of information that is required to assess the SAR value is the initial slope of the T versus t curve. A huge variety of magnetic nanomaterials (with varying morphologies and phases) has been tested for use in MHT and the corresponding SAR values have been reported in the literature. SAR values in the range of 10 to 300 kW/kg have been reported for spherical iron oxide NPs [33][34][35][36]. Different shapes and sizes of magnetic NPs have also been tested in regard to their performance in MHT technology and have produced SAR values in the range of 50 to 1000 kW/kg [37]. Importantly, different configurations of core-shell magnetic NPs have also been tested for MHT performance and have produced SAR values in the range of 60 to 400 kW/kg [38][39][40]. Finally, gold shelled core-shell spherical magnetic NPs (CoFe 2 O 4 @Au) have been tested for MHT performance and have produced SAR values in the range of 10 to 55 kW/kg [41]. Although gold shelled core-shell magnetic nanoparticles produce smaller SAR values (typically one order of magnitude smaller) compared to their gold unshelled counterparts, there is an increasing interest in developing this morphology for MHT. As well as protecting the magnetic cores against leaching, gold shelled core-shell magnetic NPs provide very unique surface functionalization capabilities, which allow for the easy anchoring of a range of bioactive molecules [42,43].
We envisage that the present report could contribute to the opening of a new chapter on breakthroughs in MHT technology by providing a new strategy for safely promoting cell apoptosis using thermotherapy, which could help MHT to achieve the status of a safe technology for use in practice. Last but not least, photothermal therapy could also benefit from the material strategy and analysis that are reported in this paper.

Materials and Methods
The gold shelled maghemite NPs were prepared using a slightly modified procedure (upscaling) that was described by Coelho et al. [19] and is summarized in Figure 1. Firstly, maghemite (γ-Fe 2 O 3 ) NPs were obtained from the coprecipitation of magnetite (Fe 3 O 4 ) NPs via acidic oxidation. Briefly, the protocol for coprecipitating the maghemite NPs comprised the use of 50 mL of solution (0.10 mol/L Fe 3+ , 0.05 mol/L Fe 2+ and 0.10 mol/L HCl), which was stirred with 250 mL of 1.0 mol/L NH 4 OH for 30 min. After washing (five times with water and once with 0.5 mol/L HNO 3 ), the obtained precipitate was boiled for 30 min with 0.50 mol/L Fe (NO 3 ) 3 . The magnetic precipitate was then separated with a permanent magnet and stirred with trisodium citrate solution (1.0 mol/L) at 80 • C for 30 min (with a molar ratio of citrate to iron of 0.1) in order to fabricate citrate-capped maghemite NPs. The obtained magnetic precipitate was washed with acetone, redispersed in water (with the pH adjusted to ∼7.4) and labeled as sample "CNP". Secondly, to obtain the gold shelled maghemite NPs, 320 mg of the citrate-capped maghemite NPs (CNP) was dispersed in 500 mL of water under ultrasonication. Then, reagents were added in the following HAuCl 4 /NaBH 4 sequence: 1.8 mL of 1.0 wt% aqueous HAuCl 4 (10 min sonication); 1.5 mL of 0.3 mol/L ethanolic NaBH 4 (10 min sonication). Four cycles of the HAuCl 4 /NaBH 4 sequence were employed to prepare the first gold shelled maghemite NP sample (CNP@Au1) and eight cycles of the HAuCl 4 /NaBH 4 sequence were performed to prepare the second gold shelled maghemite NP sample (CNP@Au2). The final dispersions that contained the gold shelled maghemite NPs started with the previously prepared samples (CNP@Au1 and CNP@Au2), which were then concentrated in a rotary evaporator, washed and redispersed in water at pH ∼7.4 with an Fe concentration of ∼40 mg Fe/mL. It is worth mentioning that the three investigated samples comprised the same maghemite core: gold unshelled maghemite (CNP); gold shelled maghemite with a thin gold shell (CNP@Au1); and gold shelled maghemite with a thick gold shell (CNP@Au2).
grid. Zetametry was employed to determine the zeta potential of the maghemite NPs, which was estimated through electrophoretic experiments using a Nano ZS ZetaSizer with a DTS 1070 disposable cuvette (Malvern Panalytical, London, UK). The results of the electrophoretic mobility were converted into zeta potential values using Henry's equation [44]. The gold and maghemite contents of the samples were determined through inductively coupled plasma optical emission spectrometry (ICP-OES) using a Perkin Elmer Optima 8000 (Perkin Elmer, Waltham, MA, USA) [45]. The gold/maghemite (g/m) ratios were obtained for samples CNP@Au1 (0.0376) and CNP@Au2 (0.0752). The room temperature magnetization curves were recorded through vibrating sample magnetometry (VSM) using the ADE Technologies model EV7 instrument (Microsense LLC, Lowell, MA, USA). To evaluate the heating performance of the synthesized magnetic samples, 1-mL aliquots of the corresponding aqueous solutions (CNP, CNP@Au1 and CNP@Au2), which all contained 40 mg Fe/mL, were introduced into vials and placed in the coil center of the MHT apparatus, which was a magneTherm system (nanoTherics Ltd., Dig Lane, UK). The The crystalline structures and the average size of the synthesized maghemite NPs were analyzed in an X-ray Miniflex 600 diffractometer (Rigaku, Tokyo, Japan) using CuKα radiation (λ = 1.541 Å) at 40 kV and 30 mA. For the X-ray diffraction (XRD) evaluation, the samples were deposited onto the surface of a zero-background sample holder. The average crystallite size of the maghemite NPs was estimated using Scherrer's formula for the broadening of the most intense XRD peak (311). Additionally, the morphology and size distribution of the maghemite NPs were assessed by transmission electron microscopy (TEM) using a JEM-2100 (JEOL, Akishima, Japan) microscope, which operated at 200 kV and was equipped with energy dispersive spectroscopy (EDS) from Thermo Scientific (Thermo Scientific, Waltham, MA, USA). The samples were dispersed in ultrapure water under ultrasonication and were further deposited onto a lacey carbon-supported copper grid. Zetametry was employed to determine the zeta potential of the maghemite NPs, which was estimated through electrophoretic experiments using a Nano ZS ZetaSizer with a DTS 1070 disposable cuvette (Malvern Panalytical, London, UK). The results of the electrophoretic mobility were converted into zeta potential values using Henry's equation [44]. The gold and maghemite contents of the samples were determined through inductively coupled plasma optical emission spectrometry (ICP-OES) using a Perkin Elmer Optima 8000 (Perkin Elmer, Waltham, MA, USA) [45]. The gold/maghemite (g/m) ratios were obtained for samples CNP@Au1 (0.0376) and CNP@Au2 (0.0752). The room temperature magnetization curves were recorded through vibrating sample magnetometry (VSM) using the ADE Technologies model EV7 instrument (Microsense LLC, Lowell, MA, USA).
To evaluate the heating performance of the synthesized magnetic samples, 1-mL aliquots of the corresponding aqueous solutions (CNP, CNP@Au1 and CNP@Au2), which all contained 40 mg Fe/mL, were introduced into vials and placed in the coil center of the MHT apparatus, which was a magneTherm system (nanoTherics Ltd., Dig Lane, UK). The temperature profiles were monitored with a fiber optic thermometer (LumaSense Technologies Inc., Denver, CO, USA) fiber-optic thermometer during the MHT procedure. Five different frequencies of the MHT apparatus were tested (112, 167, 334, 473 and 631 kHz). Likewise, five alternate magnetic field amplitudes were tested in combination with the above-mentioned frequencies (9, 11, 17, 18 and 25 mT). Figure 2a shows the XRD patterns of the synthesized samples (CNP, CNP@Au1 and CNP@Au2) and indicates the main XRD peak positions of the maghemite (JCPDS 39-1346) and metallic gold (JCPDS 04-0784). Only XRD peaks for the maghemite phase were observed in the CNP sample (gold unshelled). In contrast, in the gold shelled maghemite NPs (CNP@Au1 and CNP@Au2), a sub-group of XRD peaks for metallic gold (g) could be observed as well as the maghemite XRD peaks (m). The average crystallite sizes of the magnetic cores (d XRD ), which were calculated from the 311 XRD peak (maghemite phase) using Scherrer's formula, are listed in Table 1. It is worth mentioning that Scherrer's formula is given by d XRD = Kλ/β cos θ, where d XRD is actually correlated with the average volumetric size (average diameter) of the crystallite and K, λ, β and θ represent the shape factor (which is equal to 0.9 for spherical NPs), the employed X-ray wavelength, the corrected XRD line broadening at half of the maximum intensity (FWHM) and the Bragg angle, respectively [46]. Using the 311 XRD peak, the extracted d XRD values (see the first row in Table 1) indicated that the average size of the magnetic cores was fairly preserved following gold shelling as there were no significant differences between the d XRD values of the gold unshelled (CNP: d XRD = 11.8 nm) and gold shelled (CNP@Au1: d XRD = 12.3 nm; CNP@Au2: d XRD = 12.1 nm) samples.

Sample Characterization
All of the synthesized NPs, gold unshelled as well as gold shelled, presented nearly spherical morphologies (as shown in Figure 2b, which presents a typical TEM micrograph of CNP@Au2). Figure 2c shows a high-resolution transmission electron microscopy (HRTEM) image of a single gold shelled maghemite NP and clearly demonstrates the synthesis of the spherical core-shell nanostructures. The fast Fourier transform (FFT) analysis of the HRTEM image was carried out on a single particle and assessed the spot diffraction patterns of both the core and shell regions (see the white rectangles in Figure 2c). The data analysis indicated that the spots in the shell region corresponded to the (111) and (311) atomic planes of the cubic phase of the metallic gold, with interplanar distances of 0.24 nm and 0.12 nm, respectively. Likewise, the spots in the core region corresponded to the (111), (311) and (511) atomic planes of the spinel phase, with interplanar distances of 0.48 nm, 0.25 nm and 0.16 nm, respectively. The revealed crystal structures were also supported by the XRD data and confirmed the successful formation of the γ-Fe 2 O 3 /Au core-shell structures, as assessed using the HRTEM images.
The associated size histograms (vertical bars) of the citrate-capped (CNP) and gold shelled (CNP@Au1 and CNP@Au2) maghemite NPs, which were drawn from the recorded TEM images and the corresponding curve-fitting (solid lines) using a log-normal distribution function, are shown in Figure 2d. It is worth mentioning that the functional form of the log-normal distribution function that was used to fit the particle size histograms that were extracted from the TEM micrographs was: where d TEM and σ TEM are the mean diameter (logged data) and standard deviation (logged data), respectively [47]. The parameters that were extracted from the histogram curve fitting are summarized in Table 1. Although the mean diameter (d TEM ) and standard deviation (σ TEM ) of the gold shelled maghemite NPs in the TEM images were very much the same, they were slightly larger than the mean diameter (d TEM ) and standard deviation (σ TEM ) of the gold unshelled maghemite NPs (see Table 1). Importantly, the mode (i.e., the particle size that corresponded to the global maximum of the histogram curve) of the solid lines (blue, black and red) that are displayed in Figure 2d systematically shifted to higher values as the g/m ratio increased, which strongly suggested the growth of the gold shell layers at the increasing thickness between the CNP@Au1 and CNP@Au2 samples. The HRTEM/EDS spectrum of the CNP@Au2 sample, as shown in Figure 2e, revealed peaks that corresponded to Au and Fe, which also supported the claim that the MNPs in CNP@Au1 and CNP@Au2 presented γ-Fe 2 O 3 -Au biphasic structures.  (d) the size distributions (histograms) of the NPs in CNP (blue), CNP@Au1 (gray) and CNP@Au2 (red); (e) a typical HRTEM/EDS spectrum of a selected region in CNP@Au2; (f) the zeta potential as a function of the pH for CNP (blue squares), CNP@Au1 (black triangles) and CNP@Au2 (red triangles); (g) the room temperature hysteresis loops of CNP (blue), CNP@Au1 (black) and CNP@Au2 (red) and the corresponding magnetization at reduced field range (inset). The solid lines in (f,g) are only presented as a guide. In addition to the above assessment of particle size using both the XRD and TEM techniques, we further deepened the analysis by correlating the values that were extracted using the two employed techniques. In this regard, the average crystallite size was assessed using Scherrer's equation scales for the mean diameter (d TEM ) and standard deviation (σ TEM ) values that were extracted from the analysis of the TEM particle size histograms via  [48]. It is important to stress that d TEM XRD represents the average crystallite size of the magnetic cores, which were estimated using the parameters that were extracted from the curve fitting of the particle size histograms (see Figure 2d). In Table 1, the second row shows the values of d TEM XRD that were calculated using the values in the third and fourth rows. A more realistic comparison of the average sizes could then be performed by looking at the values in the first and second rows of Table 1. Still, the values of d TEM XRD were systematically slightly above the values of d XRD . It is worth mentioning that majority of the data that have been reported in the literature have shown d TEM values that were larger than the d XRD values [49]. In the literature, this discrepancy has usually been attributed to the presence of an amorphous dead layer (which leads to the well-known magnetic dead layer in magnetic nanomaterials) at the surface of nanomaterials, which is not picked up by XRD measurements but is captured in TEM micrographs [50,51]. Therefore, it was quite clear that our comparison of particle sizes (XRD versus TEM) should be carried out using d TEM XRD and d XRD data instead of d TEM and d XRD data. Then, in line with the above-presented arguments, the slightly larger values of d TEM XRD (d TEM XRD d XRD ), as observed in Table 1 (first and second rows), could be reasonably credited to the surface amorphous dead layer. Moreover, it is worth mentioning that the d TEM XRD /d XRD ratio increased monotonically between CNP (about 0.8%) and CNP@Au1 (about 6.5%) and CNP@Au2 (about 11.6%). It was quite obvious that this trend was due to the increasing thickness of the gold shell layer, which was not picked up while using Scherrer's formula to extract the d XRD values from the 311 XRD peak of the maghemite phase. This observation provided extra evidence for the successful fabrication of the reported core-shell structures. Figure 2f shows that the zeta potential of the synthesized MNPs was strongly pHsensitive and varied from 0 mV to −45 mV within the investigated range of pH (2-12). As expected for the CNP sample, this finding could be associated with the deprotonation of the -COOH groups (from the adsorbed citrate at the surface of the MNPs) as the pH increased, which reflected the enhancement of the negative zeta potential values due to the increasing number of negatively charged surface groups [52]. Additionally, the gold deposition onto the citrate-capped maghemite NPs yielded gold shelled MNPs with citrate molecules attached to their surface [19]. Thus, as shown in Figure 2f, the zeta potential of the gold shelled MNPs was also highly negative at physiological pH values and did not differ significantly from the value that was found for the CNP sample. Moreover, these high values of negative zeta potential indicated the increasing stability against NP aggregation that was due to particle-particle electrostatic repulsion. It is generally assumed that zeta potential values of higher than 30 mV provide stable conditions in charged colloids [53]. This finding corroborated the experimental observations that the solutions of the fabricated MNPs were very much stable with negligible zeta potential variation after months of storage.
As shown in Figure 2g, all of the maghemite-based synthesized samples (CNP, CNP@Au1 and CNP@Au2) displayed room temperature superparamagnetic behavior with negligible remanence and coercivity (see the inset). Moreover, the saturation magnetization of the samples decreased as the nonmagnetic gold shell thickness increased (see Table 1), which was in good agreement with the literature [15]. Impressively, the saturation magnetization of the CNP, CNP@Au1 and CNP@Au2 samples (last row in Table 1) could be analyzed against the relative g/m contents (ICP-OES data) that were obtained for CNP@Au1 (0.0376) and CNP@Au2 (0.0752). Assuming that all three samples had the same magnetic core (48.5 emu/g), the reductions in the saturation magnetization that was due to the increasing gold content was estimated to be about 47 and 45 emu/g for CNP@Au1 and CNP@Au2, respectively. These findings could then be compared to the values that were recorded from the magnetometry (VSM), as shown in Table 1, i.e., about 46 and 44 emu/g for CNP@Au1 and CNP@Au2, respectively. The two sets of data revealed an impressive consistency between the values that were independently recorded from two different experimental techniques (ICP-OES and VSM).

Magnetic Hyperthermia Data Analysis
The so-called Box-Lucas method is claimed to be a phenomenological model for fitting the time (t) dependence of the temperature (T) of a liquid medium that contains MNPs in suspension and is under the influence of an externally applied alternating (AC) magnetic field [26][27][28][29][30]. This data fitting approach has been systematically reported in the literature as being the most suitable method for modeling experimental T versus t data [29]. However, to the best of our knowledge, the present study demonstrated for the first time that the Box-Lucas method emerged naturally (instead of phenomenologically) from the definition of two well-established concepts. First, the thermal current (H S = dQ/dt), which flows across the interface between the MNPs and the outside liquid medium (i.e., the surrounding environment) and transfers heat (Q) at a given rate (dQ/dt) from the hot core of the MNPs to the colder liquid medium. Second, the thermal conductivity (k S ) of the MNP interface (gold shell) with the liquid medium, which is given by k S = −H S /A × (dT/dr), where A is the area that is perpendicular to the heat flow direction (r) at a given temperature gradient (dT/dr). Therefore, the thermal current that flows outward across the gold shell layer was given by: In order to move forward from Equation (1), we needed to take into account the morphology of the MNPs, i.e., their shapes and physical dimensions. Then, based on the recorded HRTEM micrographs (see Figure 2c, for instance), we took the simplest case of a spherical core-shell system with a magnetic core of radius R that was shelled by a nonmagnetic shell (interface) of thickness δ and thermal conductivity k S , as schematically By assuming the SSC of heat transfer with a constant thermal current ( ), Equation (2) could be integrated into the shell thickness ( ), i.e.: where is the average core temperature and ( < ) is the average outside temperature that extends to a radius of (see Figure 3). After integrating Equation (3), we obtained: The thermal current ( = d d ⁄ ) was now substituted into Equation (4) to provide: By integrating d = d into Equation (5), where and are the average values of the mass and specific heat of the suspension within the region of radius (surrounding medium), respectively, we obtained: which could be rearranged to provide: where the typical heat transfer time ( ) is given by = ( ]. The solution to Equation (7) was simply the so-called Box-Lucas method, which was obtained from very well-established concepts regarding thermal transfer across the interface between different media, i.e., the core of the MNPs and the liquid medium. Importantly. By taking this route to obtain Equation (7), we found a relationship between the typical heat transfer time ( ) and the characteristics of the interface ( ), which could be ultimately engineered during the synthesis protocol. It is worth mentioning that the relationship between and was inverse, meaning that in a typical MHT experiment, metal shelled MNPs would be expected to achieve the SSC later than unshelled MNPs. Additionally, the characteristics of the surrounding medium ( and ) were included in Equation (7) as well. Impressively, the solution to Equation (7) unveiled a relationship between the typical heat transfer time ( ) and the morphological (shape and size) aspects of the MNPs, as well as the thickness of the metal shell. The metal shell thickness and the shape and size of the MNPs were extremely important to consider while engineering the MNPs in order to meet specific requirements of the heat transfer time ( ). The claimed By assuming the SSC of heat transfer with a constant thermal current (H S ), Equation (2) could be integrated into the shell thickness (δ), i.e.: where T C is the average core temperature and T (T < T C ) is the average outside temperature that extends to a radius of R (see Figure 3). After integrating Equation (3), we obtained: The thermal current (H S = dQ/dt) was now substituted into Equation (4) to provide: By integrating dQ = m c dT into Equation (5), where m and c are the average values of the mass and specific heat of the suspension within the region of radius R (surrounding medium), respectively, we obtained: which could be rearranged to provide: where the typical heat transfer time (τ) is given by The solution to Equation (7) was simply the so-called Box-Lucas method, which was obtained from very well-established concepts regarding thermal transfer across the interface between different media, i.e., the core of the MNPs and the liquid medium. Importantly. By taking this route to obtain Equation (7), we found a relationship between the typical heat transfer time (τ) and the characteristics of the interface (k S ), which could be ultimately engineered during the synthesis protocol. It is worth mentioning that the relationship between τ and k S was inverse, meaning that in a typical MHT experiment, metal shelled MNPs would be expected to achieve the SSC later than unshelled MNPs. Additionally, the characteristics of the surrounding medium (m and c ) were included in Equation (7) as well. Impressively, the solution to Equation (7) unveiled a relationship between the typical heat transfer time (τ) and the morphological (shape and size) aspects of the MNPs, as well as the thickness of the metal shell. The metal shell thickness and the shape and size of the MNPs were extremely important to consider while engineering the MNPs in order to meet specific requirements of the heat transfer time (τ). The claimed phenomenological model of the Box-Lucas method failed to unveil all of the key aspects that were emphasized by the heat transfer time (τ). Finally, the solution to Equation (7) provided (analytically) the well-known Box-Lucas method: where T 0 is the average initial temperature of the liquid medium and T ∞ is the average asymptotic temperature (@ t → ∞ ) of the SSC for heat transfer. Notice that in Equation (8), It is worth mentioning that a very important piece of information was extracted from the fitting of the T versus t data (see Figure 4), namely the initial slope of the T versus t curve, i.e., (dT/dt)| t→0 . The initial slope was used to define the SAR, which could be estimated using SAR = c(m /m s ) (dT/dt)| t→0 , where c, m and m s are the solvent's specific heat capacity, mass of suspension and mass of solution (nanoparticles), respectively [54]. Although (dT/dt)| t→0 was extracted from the fitted data in practice, it should be kept in mind that Equation (8) could be used to find (dT/dt)| t→0 = (4πk s /m c )[R(R + δ)/δ](T ∞ − T 0 ). This was an interesting aspect when correlating the SAR values with the typical shapes and sizes of the NPs (e.g., R and δ in Figure 3). Within the context of the present discussion, we found that SAR ∼ = (4πk s /m s )(T ∞ − T 0 )[R(R + δ)/δ], as long as c ∼ = c . that were emphasized by the heat transfer time ( ). Finally, the solution to Equation (7) provided (analytically) the well-known Box-Lucas method: where 0 is the average initial temperature of the liquid medium and ∞ is the average asymptotic temperature (@ → ∞) of the SSC for heat transfer. Notice that in Equation (8), ∞ → as → ∞. It is worth mentioning that a very important piece of information was extracted from the fitting of the versus data (see Figure 4), namely the initial slope of the versus curve, i.e., (d d ⁄ )| →0 . The initial slope was used to define the SAR, which could be estimated using = ( ⁄ )(d d ⁄ )| →0 , where , and are the solvent's specific heat capacity, mass of suspension and mass of solution (nanoparticles), respectively [54]. Although (d d ⁄ )| →0 was extracted from the fitted data in practice, it should be kept in mind that Equation (8)

Fitting Procedure for Versus Data
Equation (8) was used to fit the collected versus data, as presented in Figure 4. Figure 4 also shows the representative data (filled symbols) and the corresponding fitting data (solid gray lines) that were obtained using Equation (8). Table 2 shows the parameters ( , ∞ and 0 ) that were obtained from the fitting data for all of the performed MHT experiments, with the typical average heat transfer time ( ) in seconds (s) and the average values of ∞ and 0 in degrees Celsius (°C). Moreover, the last column of Table 2 includes the estimated values of the specific absorption rate (SAR), which were extracted from the versus data. Notice that Table 2 also includes the characteristics of the externally applied alternating (AC) magnetic field, with field amplitude ( ) in mT and frequency ( ) in kHz. In this regard, it was important to compare the

Fitting Procedure for T Versus t Data
Equation (8) was used to fit the collected T versus t data, as presented in Figure 4. Figure 4 also shows the representative data (filled symbols) and the corresponding fitting data (solid gray lines) that were obtained using Equation (8). Table 2 shows the parameters (τ, T ∞ and T 0 ) that were obtained from the fitting data for all of the performed MHT experiments, with the typical average heat transfer time (τ) in seconds (s) and the average values of T ∞ and T 0 in degrees Celsius ( • C). Moreover, the last column of Table 2 includes the estimated values of the specific absorption rate (SAR), which were extracted from the T versus t data. Notice that Table 2 also includes the characteristics of the externally applied alternating (AC) magnetic field, with field amplitude (H) in mT and frequency ( f ) in kHz. In this regard, it was important to compare the SAR data of the gold shelled γ-Fe 2 O 3 NPs (see Table 2), which varied in the range of 2 to 75 kW/kg, to the SAR data that were reported by Sabale et al. [41] for gold shelled CoFe 2 O 4 NPs, which varied in the range of 10 to 55 kW/kg. The highest SAR value that was achieved for gold shelled CoFe 2 O 4 NPs was about 55 kW/kg (experimental conditions: 5 mg CoFe 2 O 4 /mL, 276 kHz and 34 mT). In contrast, the highest SAR value that was achieved for gold shelled γ-Fe 2 O 3 NPs (present report) was about 75 kW/kg (experimental conditions: 40 mg Fe/mL, 334 kHz and 17 mT). Actually, the 5 mg CoFe 2 O 4 /mL represented 2.4 mg Fe/mL. Therefore, by considering the inverse relationship of SAR ∝ (1/m s ) and neglecting any changes in the remaining parameters that were associated with SAR, we could estimate the SAR values for gold shelled CoFe 2 O 4 NPs with 40 mg Fe/mL, which were roughly equal to: SAR(40 mg Fe/mL) ∼ = SAR(2.4 mg Fe/mL) × (2.4/40) ∼ = 55 × (2.4/40) ∼ = 3.3kW/kg. So, roughly speaking, the gold shelled γ-Fe 2 O 3 NPs in this study presented the best SAR values (see Table 2: 75 kW/kg), which were about twenty times ((75/3.3) ∼ = 20. ) bigger than the best SAR values that have been reported for gold shelled CoFe 2 O 4 NPs [41]. Table 2. The parameters (τ, T ∞ and T 0 ) that were extracted from the fitting of the T versus t data (e.g., see Figure 4) at a fixed combination of the frequency ( f ) and amplitude (H) of the AC magnetic field. The estimated SAR values (in kW/kg-Fe) are included in last column. At this point, it was important to derive an expression for the ratio between the typical heat transfer time (τ = (m c /4πk S )[δ/R(R + δ)]) that was applied to CNP@Au2 and CNP@Au1, i.e., τ CNP@Au2 /τ CNP@Au1 = (R + e CNP@Au1 /R + e CNP@Au2 )(e CNP@Au2 /e CNP@Au1 ), and the specific mass of both gold (19.32 g/cm 3 ) and maghemite (4.88 g/cm 3 ) and the relative g/m contents that were obtained (ICP-OES data) for CNP@Au1 (0.0376) and CNP@Au2 (0.0752). These data allowed us to estimate the ratio τ CNP@Au2 /τ CNP@Au1 ∼ = 1.9 . Moreover, this finding (τ CNP@Au2 /τ CNP@Au1 ∼ = 1.9) could be compared to the ratio of the average τ-values that were included in Table 2 (τ CNP@Au2 /τ CNP@Au1 ∼ = 1.4), which were extracted from the fitting data of CNP@Au2 and CNP@Au1 under different values of AC magnetic field frequency and amplitude. The fair agreement between the two independent sets of data, i.e., the estimation of the τ CNP@Au2 /τ CNP@Au1 ratio from the assessed ICP-OES data (τ CNP@Au2 /τ CNP@Au1 ∼ = 1.9) and the fitting of the T versus t data using Equation (8) (τ CNP@Au2 /τ CNP@Au1 ∼ = 1.4), was quite impressive and pointed toward the impact of the thickness of the gold shells that were deposited on top of the maghemite cores when establishing the T versus t curves that were generated by MHT. In addition, the selected T ∞ and τ CNP@Au1 (τ CNP@Au2 ) data (dark gray and light gray boxes in Table 2) for similar operating frequencies revealed a linear relationship between T ∞ and τ, as shown in Figure 5.
∞ 0 comprehensive analysis of . At the end of the typical heat transfer time ( ), the surrounding medium reached the temperature , which was lower than ∞ by a fixed amount (∆ = ∞ − ≅ ⁄ ). The ∆ value scaled with the total temperature change (∆ 0 = ∞ − 0 ≅ ) via the following relationship: ∆ = ∆ 0 ⁄ . This means that we could fabricate MNPs to reach pre-set ∞ values at controllable times by modulating ( = ( 4 ⁄ )[ ( + ) ⁄ ]), e.g., 1 < 2 (see the inset of Figure 5), and changing the thermal conductivity of the metallic shells ( ) or the shapes and sizes of the MNPs. Once more, our findings confirmed the appropriateness of Equation (7) and its corresponding analytical solution that was represented by Equation (8) (which is known as the Box-Lucas method in the literature) for describing the versus curve for the time dependence of the temperature of a liquid medium during MHT.  Table 2) data (black circles) and ∞ versus CNP@Au2 (light gray boxes in Table 2) data (red circles). The inset shows a schematic diagram of the relationship between ∆ 0 and ∆ , i.e., ∆ = ∆ 0 ⁄ ≅ 37%∆ 0 .

Conclusions
In conclusion, the present study investigated the successful modulation of the asymptotic temperature ( ∞ ) of an aqueous suspension that contained gold shelled maghemite nanoparticles, which were tested for applications in magnetic hyperthermia. Importantly, our data showed that ∞ was achieved in less than one hour. Moreover, the typical heat transfer time ( ) could be shortened by engineering the thermal conductivity of the shells and/or the shapes and sizes of the magnetic nanoparticles, which could encourage the exploration of new structures. Importantly, the present study showed that ∞ could be pre-set for increasing gold shell thicknesses and varying operating frequencies and amplitudes of the applied AC magnetic fields. Even more importantly, the achieved  Table 2) data (black circles) and T ∞ versus τ CNP@Au2 (light gray boxes in Table 2) data (red circles). The inset shows a schematic diagram of the relationship between ∆T 0 and ∆T τ , i.e., ∆T τ = ∆T 0 /δ ∼ = 37%∆T 0 .
Actually, the linear trend that is shown in Figure 5 emerged from the first derivative of Equation (8), i.e., dT/dt = [(T ∞ − T 0 )/τ]e −t/τ . In the linear regime of the T versus t curve (t τ), we had e −t/τ ∼ = 1 with the approximately constant dT/dt = η. Therefore, in this limit, (T ∞ − T 0 )/τ ∼ = η; thus, it was expected that T ∞ would scale linearly with τ, i.e., T ∞ ∼ = T 0 + ητ, as presented in Figure 5. Finally, the inset of Figure 5 provides a more comprehensive analysis of τ. At the end of the typical heat transfer time (τ), the surrounding medium reached the temperature T τ , which was lower than T ∞ by a fixed amount (∆T τ = T ∞ − T τ ∼ = ητ/δ). The ∆T τ value scaled with the total temperature change (∆T 0 = T ∞ − T 0 ∼ = ητ) via the following relationship: ∆T τ = ∆T 0 /δ. This means that we could fabricate MNPs to reach pre-set T ∞ values at controllable times by modulating τ (τ = (m c /4πk S )[δ/R(R + δ)]), e.g., τ 1 < τ 2 (see the inset of Figure 5), and changing the thermal conductivity of the metallic shells (k S ) or the shapes and sizes of the MNPs. Once more, our findings confirmed the appropriateness of Equation (7) and its corresponding analytical solution that was represented by Equation (8) (which is known as the Box-Lucas method in the literature) for describing the T versus t curve for the time dependence of the temperature of a liquid medium during MHT.

Conclusions
In conclusion, the present study investigated the successful modulation of the asymptotic temperature (T ∞ ) of an aqueous suspension that contained gold shelled maghemite nanoparticles, which were tested for applications in magnetic hyperthermia. Importantly, our data showed that T ∞ was achieved in less than one hour. Moreover, the typical heat transfer time (τ) could be shortened by engineering the thermal conductivity of the shells and/or the shapes and sizes of the magnetic nanoparticles, which could encourage the exploration of new structures. Importantly, the present study showed that T ∞ could be pre-set for increasing gold shell thicknesses and varying operating frequencies and amplitudes of the applied AC magnetic fields. Even more importantly, the achieved values of T ∞ could be pre-set to fall within the temperature range that is required to trigger cell apoptosis instead of cell necrosis with a precision of better than 0.3 • C. Impressively, the highest SAR (specific absorption rate) value that was achieved in our experiments (75 kW/kg) was roughly twenty times bigger than the best SAR value that has been reported in the literature for a similar structure. Our findings regarding the modulation of T ∞ are expected to provide a breakthrough for establishing clinically safe protocols for magnetic hyperthermia therapy in the near future. Finally, a mathematical expression that described the time dependence of the temperature of a magnetic suspension unveiled the characteristic heat transfer time and its dependence on the shapes and sizes of the core-shell structures. This mathematical derivation, which was based on the well-known concept of heat transfer, was presented for the first time and the derived equation analytically matched the phenomenological equation in the Box-Lucas method, which has been widely used in the literature.