Abstract
The lowfield (below 5 Oe) ac and dc magnetic response of a magnetic fluid [MF] sample in the range of 305 to 360 K and 410 to 455 K was experimentally and theoretically investigated. We found a systematic deviation of Curie's law, which predicts a linear temperature dependence of inverse initial susceptibility in the range of our investigation. This finding, as we hypothesized, is due to the onset of a secondorderlike clustertomonomer transition with a critical exponent which is equal to 0.50. The susceptibility data were well fitted by a modified Langevin function, in which cluster dissociation into monomers, at the critical temperature [T*], was included. In the ac experiments, we found that T* was reducing from 381.8 to 380.4 K as the frequency of the applied field increases from 123 to 173 Hz. In addition, our ac experiments confirm that only monomers respond for the magnetic behavior of the MF sample above T*. Furthermore, our Monte Carlo simulation and analytical results support the hypothesis of a thermalassisted dissociation of chainlike structures.
PACS: 75.75.C; 75.30.Kz; 75.30.Cr.
Keywords:
magnetic fluid; inverse initial susceptibility; cluster disruption; critical temperature; Langevin functionIntroduction
The interest in magnetic fluids [MFs] has increased enormously in the last decade, particularly due to the opportunities they provide for applications in the medical field [17]. Among others, MFs have been used as an excellent material platform for the development of magnetic immunoassay [1,2], contrast agents for magnetic resonance imaging [3,4], and material devices for magnetohyperthermia [57]. Deep understanding of magnetic susceptibility, however, is a key issue not only from the fundamental point of view, but also while tailoring nanosized magnetic materials for medical applications [17]. The design of nanosized magnetic particles, taking into account the maximization of the materials' response in terms of their use for diagnosis, imaging, and therapy, requires the knowledge of the temperature dependence of the magnetic susceptibility under the action of applied dc and ac fields.
In this context, the widely accepted concept is a linear relationship between the inverse initial magnetic susceptibility (1/χ) and the temperature (T), which is accounted for by the firstorder Langevin function. Nevertheless, unusual deviations of linearity at temperatures within the range of interest for the medical applications, with no conclusive explanation yet, have been reported [8]. Since an interaction among particles in a MF sample, either modulated or not by external fields, cannot be ignored and leads to expontaneous agglomeration in clusters or chainlike structures (dimers, trimers, etc.) [911], the superlinear deviation of the 1/χ versus T curve has been attributed to magnetic dipolar interaction among the nanosized particles [12]. The Langevin function, however, includes no interaction among the suspended particles, and therefore, clusters are ruled out from the classical description. Nevertheless, at high enough temperatures, the phenomenon of thermally assisted cluster disruption within MF samples has been reported [13]. Therefore, to understand the underlying physics of this nonlinearity in magnetization, a more complete physical model is highly demanded. The new model should take into account both monomers and clusters in MFs.
In this study, we report the unusual superlinear deviation in the temperature dependence of the inverse initial magnetic susceptibility in a magnetitebased (Fe_{3}O_{4}) MF sample in the temperature range of 305 to 360 K. The observed breakdown of Curie's law, which scales linearly the inverse susceptibility with temperature, indicates that besides the usual tendency of alignment of magnetic moments with the applied field, there exists an additional thermally assisted physical process connected to cluster disruption within MF samples. Therefore, we propose a model in which a chainlike disruption at a typical transition temperature (T*) is incorporated. We found that the extended model reproduces quite well the observed superlinear deviation. Additionally, the abovementioned superlinear deviation of the (1/χ) × T data for chains of particles is supported by Monte Carlo [MC] simulation and herein incorporated.
Experimental description
In order to verify the universality of the experiment regarding the relationship of inverse initial susceptibility and temperature, we explored the experiments under dc and ac magnetic fields, the latter at different frequencies. The MF sample used in our experiment was a commercial magnetic colloid (EFH1, Ferrotec Corporation, Santa Clara, CA, USA), consisting of Fe_{3}O_{4 }magnetite nanoparticles (mean particle diameter of 10 nm) suspended in light mineral oil. The applied magnetic fields were 5 and 2 Oe (amplitude) for dc and ac experiments, respectively. The ac experiments were performed at 123 and 173 Hz.
The chord magnetic susceptibility (χ = M/H) can be obtained by measuring the sample's magnetization (M) induced by a weak external dc/ac field (H). The homemade experimental system designed to measure the MF sample's initial susceptibility is schematically shown in Figure 1. The Helmholtz coils were driven by a dc/ac power supply to generate a uniform dc/ac excitation field at the sample's position; the direction of the dc applied field can be altered by a relay. The temperature of the MF sample was measured using a thermocouple (Pt100) and converted using a temperature transmitter (RS1852290). A weak magnetization signal collected from the MF sample was detected using the twoaxis giant magnetoresistive sensors (HMC1022) with an amplification circuit. Both temperature and magnetization data of the sample were acquired through a data acquisition card (PCI6251, National Instruments, Austin, TX, USA). LabVIEW was used as a humancomputer interface for observation and control.
Figure 1. Schematic diagram of the measurement system suitable for dc/ac initial magnetic susceptibility detection.
The experiment includes three steps: Firstly, the sample is placed inside a water bath, slowly heated up to the water's boiling point temperature, and kept there for a certain period of time. While keeping the sample holder stably hot, a weak uniform dc/ac field was applied to the MF sample. The second step is to record the inverse susceptibility (χ^{1 }= H/M) versus temperature curves. The heating circuit is shut off, allowing the sample of cooling down naturally while synchronously recording both the temperature and the chord susceptibility. To avoid stray magnetic field interference, when the dc field is applied, the direction was set to invert every 0.5 s. The magnetization was evaluated according to the difference between the two measurements, before and after switching the relay. For ac susceptibility measurement, the digital phasesensitive detection algorithm was used to calculate the amplitude of the magnetization. The last step is data processing. Regarding the data obtained from repeated experiments at different excitation fields, we took a temperature point every 0.5 K in the temperature range of the experiments, and the measured values at temperature points within a permissible measurement error range are considered to be achieved at the same temperature. Accordingly, the mean values of repeated experiments were obtained to draw the curves of inverse susceptibility versus temperature. Typical (1/χ) × T data are shown in open symbols in Figures 2 and 3. Note the superlinear behaviors of the (1/χ) × T data at the higher temperature end in both Figures 2 and 3.
Figure 2. Experimental data and the fitting curves using the Langevin function. Experimental data under a 5Oe dc external field (open squares) and the fitting curves using the firstorder Langevin function (dashed line) and the modified Langevin function (solid line). Inset (a) is a zooming in the temperature range of 350 to 360 K. Inset (b) is the temperature dependence of P_{2}.
Figure 3. Experimental data and the best fitting curves. Experimental data under an ac external field (2Oe amplitude) with frequencies set at 123 and 173 Hz (open squares and open circles, respectively) and the best fitting curves using Equation 5 (solid lines).
Model and discussion
The susceptibility model for pure monomers
Considering only isolated particles (monomers) with no interparticle interaction, the MF sample's magnetization is well described by the Langevin function:
where φ is the concentration of the nanoparticle (the number of particles per unit volume), m is the nanoparticle's magnetic moment, H is the applied field, k is the Boltzmann constant, and T is the absolute temperature. The nanoparticle's magnetic moment (m) is described in terms of the saturation magnetization (M_{S}) as m = M_{s}V, where V is the nanoparticle's volume. At weak enough applied fields, the sample's susceptibility χ = ∂M/∂H approaches the chord susceptibility χ = M/H. Then, the inverse initial susceptibility can be written as:
Equation 2 describes the (1/χ) × T data based on the firstorder Langevin function. Note that the dashed line in Figure 2 is the curve fitting of the data (open symbols) using Equation 2, showing the expected linearity in the temperature range of interest. However, our experimental data (open symbols in both Figures 2 and 3) revealed a superlinear trend at the high temperature end, making the fitting procedure using Equation 2 visibly poor (see inset (a) of Figure 2). We then hypothesized that the assumption of only monomers in the MF sample in the temperature range of 305 to 360 K no longer holds.
Discussion of the susceptibility model including dimers
The explanation of the superlinear behavior observed on the (1/χ) × T experimental data displayed in Figures 2 and 3 (open symbols), as we claim, is due to the presence of a fraction of chainlike structures within the MF sample, in addition to monomers, more likely dimers [13,14], and the temperature dependence of the dimer fraction in a critical way around T*, as discussed below.
The probability of agglomeration and disruption of suspended nanoparticles in MFs is assumed to be dependent upon the relative strength of magnetic, van de Waals, electrostatic, and steric interactions and thermal energy. The literature [1517] describes that magnetic dipolar interaction held particles together, whereas electrostatic interaction and thermal energy work together taking nearby particles apart. It can be inferred from the experimental observations that there is a dynamic balance between the relative content of monomers and dimers within the simplest model picture of a magnetically textured MF sample. This balance will be broken when there is a change on the parameters governing the energy terms involved, and the relative content of monomers and dimers will change accordingly. As the thermal energy increases, the dimers tend to disrupt into monomers, leading to a decrease of the dimers' relative content. Inversely, within this simplest model picture, when lowering the temperature of a MF sample, the monomers tend to agglomerate into dimers, leading to a decrease of the monomers' relative content.
The already described process of the thermalassisted dimer disruption within a MF sample [12] is related to a secondorder phase transition at a critical temperature T*. Below T*, the suspended nanoparticles are found as monomers and dimers, whereas above T*, nanoparticles are essentially isolated within the MF sample while the relative content of dimers drops down critically to nearly zero around T*. Below the critical temperature T*, monomers and dimers do coexist within the MF sample. In order to describe this mixed (magnetically textured) system and the influence of the temperature upon it, we have modified Equation 1 including two contributions, namely terms due to monomers and dimers. Below T*, the probability of finding monomers (P_{1}) and dimers (P_{2}) within the MF sample scales with the temperature according to:
Considering the classical Landau's theory of secondorder phase transition, β in Equation 3 represents the critical exponent (β = 0.5) as long as P_{2 }properly describes the order parameter of the thermally assisted dimer disruption process. Thus, the modified Langevin function describing the actual magnetization of the MF sample is:
where m_{1 }and m_{2 }are the magnetic moments of monomers and dimers, respectively. In Equation 4, φP_{1 }and φP_{2 }represent the content of monomers and dimers, respectively. For abbreviation, we write ξ_{1 }= m_{1}H/kT and ξ_{2 }= m_{1}H/kT. According to Equations 3 and 4, the inverse susceptibility describing the magnetically textured MF sample can be written as:
The lowfield dc susceptibility data are shown in Figure 2. The solid line represents (see Figure 2) the best curve fitting of the (1/χ) × T experimental data (open squares) using Equation 5. Included here for comparison, the dashed line in Figure 2 represents the fitting of the data using Equation 2. Note that ferromagnetic resonance data (resonance line splitting) have been used to describe dimer disruption in a nickel ferritebased ionic MF sample, providing values of T* (340 K) and β (0.42), the latter in reasonable agreement with the value we found in the present study [13]. Our data indicate that above T*, the suspended particles within the MF sample investigated are essentially isolated (monomers). The inset (a) of Figure 2 shows a detail of the (1/χ) × T data, emphasizing the superlinear behavior at temperatures close to the typical dimer disruption temperature (T*). The inset (b) of Figure 2 shows the temperature dependence of P_{2}, here, describing the order parameter associated to a secondorderlike phase transition.
Experiments on lowfield ac susceptibility at different frequencies (123 and 173 Hz) are shown (symbols) in Figure 3. The saturation magnetization of monomers (M_{S1}) and dimers (M_{S2}) in the same MF sample was considered to be constant. M_{S1}and M_{S2}were firstly acquired by the fitting of any applied field (dc or ac) and frequency, and then both M_{S1}and M_{S2}were used as known parameters to fit different data sets. We found different susceptibility responses while changing the frequency of the applied field. The fittings of the experimental data (solid lines) using Equation 5 show that the critical temperature decreases from T* = 381.8 down to T* = 380.4 as the frequency of the ac excitation field increases from 123 to 173 Hz, in agreement with recently reported results [18].
Monte Carlo simulation of secondorderlike dimermonomer transition
Because a dipolar interaction between neighboring magnetic nanoparticles favors parallel alignment of their magnetic moments, a strong interaction may induce coalescence of nanoparticles and formation of magnetically coherent clusters such as dimers. Then, monomers and dimers may coexist in the system with an intermediate concentration of nanoparticles. This model picture allows us to describe the system using the Ising model in which we consider N_{p }identical magnetic nanoparticles disposed in onedimensional chains. We assume that each nanoparticle can be represented by a magnetic monodomain, with a magnetic moment along the direction of an external magnetic field, where σ_{i }= ± 1, i = 1, 2,..., N_{p}. Then, the total energy of the nanoparticle assembly is given by
where the sum is extended to the nearest neighbors only, J is the nearest neighbor exchange integral, and h= M_{s}H. The equilibrium magnetic moment configuration is obtained by a MC simulation using the standard Metropolis algorithm. It is interesting to find out that for systems such as MFs in which suspended nanoparticles are found as monomers and dimers, the magnetization and susceptibility can be calculated analytically through the partition function of the system, , where E_{n }is the energy of the system in the nth configuration. For instance, we can use the following four configurations in order to describe dimers σ_{1}σ_{2}>=++>,+>,+>,>.
Temperaturedependent disruption process of dimers, i.e., variation on the number of dimers (N_{d}) in MFs, can be modeled by a parameter γ, which is defined as , where N_{m }is the number of monomers, N_{p }= N_{m}+2N_{d}, and T_{i }(T_{f}) is the initial (final) temperature of the simulation. Notice that γ describes the system in a wide range of configurations in regard to the disruption process, for instance, from a MF with the same number of monomers and dimers (N_{d }= N_{m}) at T = T_{i }to a configuration in which dimers are completely disrupted (N_{d }= 0) at T = T_{f}. Then, the system's susceptibility consisting of monomers and dimers (complex system) is given by
where χ_{m }(χ_{d}) is the susceptibility of the monomers (dimers). Figure 3 shows the temperature dependence of the inverse susceptibility obtained by MC simulation for monomers (dashed line), dimers (dashdotted line), and their complex (solid line) at μH/J = 0.1, kT_{i}/J = 0.5, and kT_{f}/J = 10. Notice that the inverse susceptibility of monomers is larger than that of dimers, and both of them show almost linear dependence on temperature. However, an unusual superlinear deviation is observed in the monomerdimer complex curve. From Equation 7, one sees that χ mainly depends upon χ_{m}, χ_{d}, and γ. Therefore, the variation of γ with temperature, i.e., the change of the number of dimers in a MF sample, results in the deviation of the susceptibility of MFs from linearity. Finally, the analytical solution of the monomerdimer complex system, as shown in inset of Figure 4, gives a further strong support for this model picture.
Figure 4. Inverse susceptibility of monomers, dimers, and their complex. Inverse susceptibility of monomers (dashed lines), dimers (dotdashed lines), and their complex (solid lines) as a function of temperature for μH/J = 0.1, obtained by MC simulation. Inset presents the correspondent analytical solution.
The pure monomer scenario above T*
Experimental investigation of the inverse susceptibility versus temperature was also performed in the temperature range above T*. In order to prevent volatilization of the MF solvent (mineral oil), the highest temperature achieved in our experiment was restricted to 455 K. According to the analysis presented above, all dimers are assumed to be disrupted into monomers when the temperature of the MF sample is set above the critical temperature T*. This means that above T*, the relative content of dimers tends to be 0% whereas the relative content of monomers tends to be 100%. The inverse susceptibilitytemperature curve in the temperature range of 410 to 455 K is plotted in Figure 5. The experimental data shown in Figure 5 were fitted by the monomer model, herein described by:
Figure 5. Experimental data under ac external field (2 Oe and 83 Hz) above T* (open squares). The solid line represents the best curve fitting using Equation 8.
Note that A in Equation 8 means just a scale factor for the experimental data. The solid line in Figure 5 has an excellent agreement with the experimental data and presents the best curve fitting achieved using Equation 8. Furthermore, the monomer's saturation magnetization (M_{S}) obtained from the data recorded in the temperature range above T* (from 410 to 455 K) is very much close to the value found from the data recorded in the temperature range below T* (from 305 to 360 K). Although fitting of the data presented in Figure 5 (open symbols) can be performed using the dimermonomer model, represented by Equation 5, there is no agreement between the M_{S1}obtained from below and above T*. This indicates that dimers can be assumed to be totally disrupted into monomers when the temperature is increased above T*, and the critical temperature T* can be actually used to describe a MF sample.
Conclusion
In conclusion, the usual linear temperature dependence of the inverse initial susceptibility (dc or ac) of MFs at lower temperatures is found whereas at higher temperatures, an upward deviation of linearity is observed. This superlinear behavior is attributed to a thermalassisted disruption of dimers into monomers which is described by the classical Landau's approach for secondorder phase transitions. The experimental observations in the temperature ranges below and above the critical temperature T* are well fitted by the model picture proposed here, in which an extended Langevin function including the thermal criticality of the MF system is adopted below T*. Our findings are strongly supported by both MC simulation and analytic analysis.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
WL conceived the project. PCM contributed the magnetization theory and model; WL and JZ designed the experimental system; and JZ and QX performed the experiments. LOM and FQ performed the Monte Carlo simulation. All authors wrote the manuscript, discussed the results, and commented on the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The work was supported by NSFC 61174008 and 11104089, MOST 0S2012GR0121, R&D HB 2010BFA013, and GBIE HF09062011184.
References

Öisjöen F, Schneiderman JF, Astalan AP, Kalabukhov A, Johansson C, Winkler D: The need for stable, monodispersed, and biofunctional magnetic nanoparticles for onestep magnetic immunoassays.
J Phys: Conf Ser 2010, 200:122006. Publisher Full Text

Enpuku K, Soejima1 K, Nishimoto T, Kuma H, Hamasaki H, Tsukamoto A, Saitoh K, Kandori A: Quantitative evaluation of magnetic immunoassay with remanence measurement.
Supercond Sci Technol 2006, 19:S257S260. Publisher Full Text

Jeong SY, Kim HJ, Kwak BK, Lee HY, Seong H, Shin BC, Yuk SH, Hwang SJ, Cho SH: Biocompatible polyhydroxyethylaspartamidebased micelles with gadolinium for MRI contrast agents.
Nanoscale Res Lett 2010, 5:19701976. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Chen Z, Yu D, Wang S, Zhang N, Ma C, Lu Z: Biocompatible nanocomplexes for molecular targeted MRI contrast agent.
Nanoscale Res Lett 2009, 4:618626. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Nair BG, Nagaoka Y, Morimoto H, Yoshida Y, Maekawa T, Kumar DS: Aptamer conjugated magnetic nanoparticles as nanosurgeons.
Nanotechnology 2010, 21:455102. PubMed Abstract  Publisher Full Text

Kok HP, de Greef M, Wiersma J, Bel A, Crezee J: The impact of the waveguide aperture size of the 3D 70 MHz AMC8 locoregional hyperthermia system on tumour coverage.
Phys Med Biol 2010, 55:48994916. PubMed Abstract  Publisher Full Text

Kita E, Oda T, Kayano T, Sato S, Minagawa M, Yanagihara H, Kishimoto M, Mitsumata C, Hashimoto S, Yamada K, Ohkohchi N: Ferromagnetic nanoparticles for magnetic hyperthermia and thermoablation therapy.
J Phys D: Appl Phys 2010, 43:474011. Publisher Full Text

Pan Y, Liu X, Du C, Li Z: An effective correction to the theoretical curve of magnetobirefringence of magnetic fluid.
J Appl Phys 1992, 71:39373939. Publisher Full Text

Wang A, Li J, Gao R: The structural force arising from magnetic interactions in polydisperse ferrofluids.
Appl Phys Lett 2009, 94:212501. Publisher Full Text

Liu Y, Liu X, Wang X: Biomimetic synthesis of gelatin polypeptideassisted noblemetal nanoparticles and their interaction study.

Fang WX, He ZH, Xu XQ, Mao ZQ, Shen H: Magneticfieldinduced chainlike assembly structures of Fe_{3}O_{4 }nanoparticles.
Europhys Lett 2007, 77:68004. Publisher Full Text

Jeun M, Bae S, Tomitaka A, Takemura Y, Park KH, Paek SH, Chung KW: Effects of particle dipole interaction on the ac magnetically induced heating characteristics of ferrite nanoparticles for hyperthermia.
Appl Phys Lett 2009, 95:082501. Publisher Full Text

Morais PC, Gonçalves GRR, Bakuzis AF, Skeff Neto K, Pelegrini F: Experimental evidence of dimer disruption in ionic ferrofluid: a ferromagnetic resonance investigation.
J Magn Magn Mater 2001, 225:8488. Publisher Full Text

Liu W, Zhong J, Xiang Q, Yang G, Zhou M: Discretization of magnetization curves and their application in size estimation of nanosized ferrofluid.

Skeff Neto K, Bakuzis AF, Morais PC, Pereira AR, Azevedo RB, Lacava LM, Lacava ZGM: The influence of aggregates and relative permeability on the magnetic birefringence in ionic magnetic fluids.
J Appl Phys 2001, 89:33623369. Publisher Full Text

Bakuzis AF, Skeff Neto K, Gravina PP, Figueiredo LC, Morais PC, Silva LP, Azevedo RB, Silva O: Magnetooptical properties of a highly transparent cadmium ferritebased magnetic fluid.
Appl Phys Lett 2004, 84:23552357. Publisher Full Text

Morais PC, Gravina PP, Skeff Neto K, Lacava BM: Birefringence of maghemitebased magnetic fluid in the highly concentrated regime.
Thin Sol Films 2006, 515:226230. Publisher Full Text

Xiang Q, Zhong J, Zhou M, Morais PC, Liu W: AC field dependence of cluster disruption in magnetic fluids.
J Appl Phys 2011, 109:07B317. Publisher Full Text