1. Introduction
The giant magnetic moment (the superspin) of a nanoparticle in an ideal, non-interacting ensemble relaxes according to the Néel–Brown model, which holds that the time (τ) it takes the superspin to flip along an easy magnetization axis depends on the temperature (T) according to an Arrhenius-type activation law [
1,
2]:
where τ
0 is a characteristic time (typically of the order of 10
−10 to 10
−9 s), k
B = 1.380649 × 10
−23 J/K is the Boltzmann constant, and E
B is the energy barrier to magnetization (superspin) reversal. At temperatures above a certain threshold, T
B (the blocking temperature), a so-called superparamagnetic transition occurs: τ becomes shorter than the observation time τ
obs used to perform a magnetic measurement, so the ensemble exhibits (super)paramagnetic properties although the nanoparticles’ material is ferro- or ferri-magnetic. τ
obs is determined by the measurement type, e.g., τ
obs~1 s for dc-magnetization, while τ
obs = 1/2πf for ac-susceptibility measurements (where f is the frequency of the driving magnetic field). Therefore, the blocking temperature T
B depends on the energy barrier to superspin reversal E
B, as well as on the measurement technique.
In recent years, a lot of effort has been invested in studying the behavior of superparamagnetic nanoparticle ensembles. Both experimental and computational investigations have been carried out, and significant progress has been made in developing new and enhancing existing applications of these systems in imaging [
3,
4,
5], high-density magnetic recording [
6,
7], and medicine [
8,
9,
10]. One important observation that has emerged from these studies is that the functionality of many of the magnetic nanoparticle applications critically depends on the ability to determine and control T
B and E
B. For example, magnetic recording devices based on nanoparticle ensembles catastrophically fail (i.e., all the stored information is lost) at temperatures above T
B as the system becomes superparamagnetic [
11]. For fine nanoparticles of average size smaller than 10 nm that can occur below room temperature. For medical applications, an insufficient E
B value is the main factor that prevents magnetic nanoparticles from effectively functioning as stand-alone hyperthermia cancer therapy agents [
12]. Heat dissipation in superparamagnetic nanoparticle ensembles is not hysteretic but largely depends on the superspin flip across the energy barrier to magnetization reversal. For the ideal nanoparticle ensemble described above, predicting E
B and T
B is straightforward, as E
B = KV and
, where K is the magneto-crystalline anisotropy constant of the material and V is the nanoparticle’s volume. Real cases, however, are significantly more complicated. Even for nanoparticle ensembles with sharp size distributions and negligible interparticle interactions, measurements have revealed barriers to superspin reversal that are significantly greater than KV. This was ascribed to an additional surface component to the anisotropy constant, a hypothesis confirmed by several studies [
13,
14,
15]. In commonly used nanoparticle ensembles (e.g., in nanopowders), interparticle dipolar interactions also play an important role in determining the barrier to magnetization reversal and the blocking temperature, together with other factors such as the particle size distribution, geometrical arrangements of the particles, and orientation of the easy axes [
16,
17].
The effect of interparticle dipolar interactions on the collective superspin relaxation has been characterized by several empirical and phenomenological models. For example, the relative peak temperature variation per frequency decade,
, is typically used to assess the strength of the interparticle interactions in a magnetic nanoparticle ensemble. ΔT = T
2 − T
1 is the shift of the peak temperature between two in-phase magnetic susceptibility χ′ vs. T|
f curves measured at different frequencies, f
1 and f
2, of the driving magnetic field, whereas Δlog(f) = log(f
2) − log(f
1). ϕ has been used to distinguish between weak interactions (0.1 < ϕ < 0.2) leading to blocking–unblocking superparamagnetic transitions [
18] and stronger interactions (ϕ < 0.05), which can lead to the collective freezing of the superspins in a spin-glass fashion [
19,
20]. Within the former ϕ range, the temperature dependence of the superspin relaxation time is well described by a Vogel–Fulcher law,
This is very similar to the Néel–Brown activation law in Equation (1), with the exception of the additional parameter T
0, which describes the strength of the interparticle interactions [
21,
22]. It is important to note the phenomenological nature of this description of the superparamagnetic relaxation in the presence of weak dipolar interactions. Indeed, T
0 has no physical meaning, and the exponent in the Vogel–Fulcher law diverges as T
0 increases and approaches T. The latter confirms that the law does not apply for strong interactions.
More recently, Dormann, Bessais, and Fiorani (DBF) developed a physical model for the magnetic relaxation of a nanoparticle ensemble based on a statistical calculation of the dipolar energy [
23]. The DBF equation:
introduces an additional component to the energy barrier to magnetization reversal, E
ad, and a variable pre-factor, τ
r. The model predicts that E
ad and, implicitly, T
B increase with the increase of the interparticle interactions [
24,
25], and, for negligible interactions (i.e., E
ad = 0), it expectedly yields the Néel–Brown equation. Interestingly, a previous model by Mørup and Tronc (MT) based on calculations of the average dipolar field 〈B
i〉 that acts on a nanoparticle’s magnetic moment μ predicts a blocking temperature T
B that actually decreases with the increasing strength of the interparticle interactions [
26]. It is important to mention, however, that the predictions of the two models are not necessarily incompatible with one another, as the MT theory was developed for very weak interparticle interactions, i.e., μ〈B
i〉 << KV. Experiments aimed at validating such models are difficult, as disentangling the contributions to the superspin relaxation by different factors that act in non-ideal nanoparticle ensembles—such as the bulk and surface anisotropy, interparticle dipolar interactions, geometrical arrangements, and size distributions—is a very complex task. Moreover, synthesizing solid samples (e.g., nanopowders) that allow the interparticle distance to be varied within a wide range while avoiding agglomeration and ensuring a sharp size distribution, is also particularly challenging.
Here, we used ac-susceptibility measurements to investigate the magnetic relaxation of Fe3O4 nanoparticles of average diameter 〈D〉 = 8 nm dispersed in different volumes of a carrier fluid (hexane). By using this method, we managed to synthesize samples with low volume concentrations, c, right above the value that renders the interparticle interactions negligible (c = 4 × 10−4% v/v). Our goal was to test the predictions of the DBF model in highly diluted nanoparticle ensembles with weak interparticle dipolar interactions, which cannot be carried out using nanopowders. We used Fe3O4 because it represents a model system for the study of fundamental processes in magnetic nanoparticle ensembles. We collected both in-phase and out-of-phase susceptibility frequency-resolved data (χ′ vs. T|f and χ″ vs. T|f, respectively) upon heating on five samples of different concentrations ranging from 4 × 10−4% to 4 × 10−2% v/v. Our first significant finding is that the blocking temperature, TB, obtained from the peak of the χ′ vs. T curves, increases with increasing c as predicted by the DBF model. We also found that the observed TB vs. c dependence is excellently described by a power law TB = A∙cγ, with A = 64 K and γ = 0.056. We then analyzed the temperature dependence of the relaxation time, τ(T), obtained from the shift with frequency of the χ′′ vs. T|f data, in the framework of the DBF model. Fits of Equation (2) to the observed τ(T) datasets allowed us to determine the total energy barrier to magnetization reversal, EB + Ead, for the five nanoparticle ensembles of different concentrations. We found EB/kB = 366 K and a reduced additional energy barrier Ead/kB that increases linearly with the common logarithm of the volume concentration, from c = 8.3 × 10−4% to 745 K at c = 4 × 10−2%. Our results are significant because they confirm the predictions of the DBF model in highly diluted ensembles of superparamagnetic nanoparticles and reveal quantitative information about the TB vs. c dependence within this low concentration regime. They also provide a means to determine the contributions to the barrier to superspin reversal from the bulk and surface magnetocrystalline anisotropy and the interparticle interactions at different values of c. This opens new paths for the rational design of nanoparticle ensembles for magnetic recording and biomedical applications.
3. Results and Discussion
Figure 1a shows the temperature dependence of the in-phase susceptibility measured on the S3 sample (1.7 × 10
−3%
v/
v) at five different frequencies: 100 Hz, 300 Hz, 1000 Hz, 3000 Hz, and 10,000 Hz.
The main feature of these χ′ vs. T|
f curves is that they exhibit robust peaks, whose temperature increases with the increase in the measurement frequency. This allows the calculation of the relative peak temperature variation per frequency decade,
. We found ϕ = 0.15 for the S3 sample, a value that shows the superparamagnetic nature of the transition undergone by the nanoparticle ensemble upon heating above a blocking temperature T
B that depends on the measurement frequency. Similar χ′ vs. T|
f data and analyses on samples S1, S2, S4, and S5 yielded ϕ values within the 0.1–0.2 range. This confirms that for all the concentrations used here, the magnetic relaxation occurs in the presence of weak interparticle interactions, and the ensembles exhibit superspin blocking–unblocking transitions. The out-of-phase susceptibility curves χ″ vs. T|
f measured on sample S3 are shown in
Figure 1b. The magnitudes of these curves are reversed compared to their in-phase counterparts (i.e., the χ″ vs. T|
10,000 Hz peak has the largest magnitude), but the peak temperature still increases with the increase in the frequency f or, equivalently, the observation time τ
obs = 1/2πf. For this reason, χ″ vs. T|
f curves will be used to determine the observed τ(T) dependence in each of the five nanoparticle ensembles used in this study.
To determine which model to use to analyze the observed temperature dependence of the superspin relaxation, we first studied the variation of the blocking temperature with the nanoparticle ensemble concentration, T
B vs c.
Figure 2 shows the χ′ vs. T|
100 Hz curves measured on the five samples of different concentrations c.
As indicated above, these curves peak at the blocking temperature corresponding to the transition from the blocked to the superparamagntic state of the nanoparticle ensemble. Therefore, the data demonstrate that T
B increases monotonically with the increase in c, from T
B = 39 K at c = 4 × 10
−4% (S5) to T
B = 53 K at c = 4 × 10
−2% (S1). This is the type of behavior predicted by the DBF model, which might seem somewhat unexpected given the high level of dilution of the samples used in our work, where the concentrations are below 4 × 10
−2%
v/v. It is worth noting, however, that denser Fe
3O
4/hexane ferrofluids of
w/
v concentration 20 mg/mL (i.e., 4 × 10
−1%
v/v) have been shown by previous studies [
27] to correspond to strong interactions that lead to a collective superspin-glass-like freezing upon cooling below a critical temperature. This is markedly different from the blocking-unblocking superparamagnetic transition observed by us and cannot be analyzed in the framework of the DFB model. On the other hand, we cannot exclude a certain level of particle agglomeration upon the freezing of the carrier fluid in our studies, which is the main reason for which we report our results in terms of the volume concentration at synthesis and not in terms of the interparticle distance calculated assuming a uniform distribution of the nanoparticles.
We then investigated the quantitative aspects of the T
B vs. c behavior.
Figure 3 shows the blocking temperature measured as a function of the ferrofluid’s volume concentration (solid symbols) obtained from the ac-susceptibility data in
Figure 2. We made several attempts to fit these data using basic functions and found that that a power law:
excellently describes the observed behavior of the blocking temperature. Indeed, a least-squares fit of this newly determined equation to the T
B vs. c data converges to low residuals upon the simultaneous variation of parameters A and γ. The solid line shows the best fit that yields A = 64 K and γ = 0.056.
Next, we focus our analysis on the temperature dependence of superspin relaxation time, τ(T), for samples of different concentrations. We will then use those data to gain quantitative information on the energy barrier to magnetization reversal and finally establish how the barrier changes with the concentration.
Figure 4a shows the shift with frequency of the χ″ vs. T|
f curve peak temperature for the S3 sample. As the system’s relaxation time, τ, is the same as the observation time, τ
obs, (i.e., τ = τ
obs = 1/2πf) at the peak temperature, these curves allow us to determine the τ(T) dependence for the S3 sample. This is shown by the solid symbols in
Figure 4b, where the natural logarithm of τ is plotted as a function of inverse temperature 1/T. We then analyzed theses data in the framework of the DBF theory. We used Equation (2), where the relaxation time is calculated based on three quantities: the pre-factor τ
r, the energy barrier for the individual nanoparticle relaxation according to the Néel and Brown model E
B, and the additional energy barrier, E
ad, which corresponds to the collective superspin relaxation due to interparticle interactions. Two parameters, τ
r and E
ad, were allowed to vary simultaneously in the fit of the DBF equation to the data, while E
B/k
B was kept at 366 K for reasons explained below. The best fit, shown by the solid line, converges to low residuals and yields τ
r = 1.1 × 10
−12 s and E
ad/k
B = 289 K.
We determined the individual energy barrier to magnetization reversal, EB, as follows. For the lowest concentration, sample S5 (c = 4 × 10
−4%), we managed to successfully fit the observed τ(T) behavior using the Néel–Brown equation
with parameters τ
0 = 5.6 × 10
−10 s and E
B/k
B = 366 K. This is significant because it indicates that the Fe
3O
4 nanoparticle ensemble in S5 is diluted enough, so each particle relaxes individually across an energy barrier, E
B, that only depends on the bulk and surface magnetic anisotropy. The key value that supports this conclusion is the vertical axis intercept of the lnτ vs. 1/T dependence, lnτ
0, which in this case yields a τ
0 between 10
−9 and 10
−11 s. Values of the inverse attempt frequency (τ
0 = 1/2πf
0) within this range are a known signature of non-interacting nanoparticle ensembles where each superspin relaxes independently [
20,
27]. In addition, we note that E
B/k
B has the same value in all five samples because they were prepared by progressively diluting the same nanopowder.
Figure 5 shows the temperature dependence of the relaxation time measured on samples S1 (upright triangles), S2 (diamonds), S3 (inverted triangles), S4 (circles), and S5 (squares), using χ″ vs. T|
f data, such as the ones in
Figure 4a. The main feature here is that the slope of the lnτ vs. 1/T plots increases with the increase in c, indicating that the superspin needs to overcome a higher energy barrier in denser samples. The solid lines are best fits of Equation (2), which allow the energy barrier to superspin reversal E
B + E
ad to be determined. As discussed above, we found E
ad = 0 K and E
B/k
B = 366 K for the S5 sample. For the other four samples, the best fits yield additional energy barriers that increase with the increase in concentration, from E
ad/k
B = 138 K for S4 (c = 8.3 × 10
−4%) to E
ad/k
B = 745 K for S1 (c = 4 × 10
−2%). These results are summarized in the inset table and shown by the solid symbols in
Figure 6. The solid line here is a guide to the eye showing that the additional barrier to superspin reversal, E
ad/k
B, increases linearly with the common logarithm of the concentration log
10c. We note that upon the increase in the volume concentration, the additional energy barrier reaches values comparable to its E
B/k
B counterpart at c~5 × 10
−3%. For sample S1, we found E
ad/k
B > 2E
B/k
B, which clearly demonstrates the strong effect of the interparticle interactions even at low nanoparticle volume concentrations, below 10
−1%. At the microscopic level, the observed dependence of the barrier to magnetization reversal on concentration indicates that, up to c = 4 × 10
−4%, the nanoparticles’ magnetic moments (the superspins) only need to overcome E
B (i.e., the barrier resulting from the bulk and surface magneto-crystalline anisotropy) in order to flip along an easy magnetization axis via thermal activation. This is the same as the case of an isolated magnetic nanoparticle, as described by the Néel–Brown model. Further increasing the concentration leads to non-negligible magnetic dipolar interactions among the nanoparticles in the ensemble. This collective relaxation requires more energy for the superspin flips. In the DBS model, this is equivalent to an additional energy barrier to superspin reversal that increases with the concentration and becomes dominant above c = 4 × 10
−4%.
Our findings are significant for several reasons. First, they demonstrate that despite the high dilution level (4 × 10−2% ≥ c ≥ 8.3 × 10−4%), samples S1–S4 still exhibit a collective relaxation of the superspins, which is typically indicative of interparticle interactions. This is unequivocally confirmed by the existence of the additional energy barrier to magnetization reversal Ead that increases with the increase in c. The ability to identify the concentration range where this behavior occurs is a new result that is particularly important for applications of magnetic nanoparticles in high-density magnetic recording, where each superspin acts as one recording bit that needs to flip individually. Second, our observation that the nanoparticles in sample S5 relax individually allowed us to determine EB and, more importantly, separate the EB and Ead components in the denser samples. This is important because it adds to our understanding of the contributions of different factors that influence the superspin relaxation, such as bulk and surface magnetic anisotropy, nanoparticle size and size distribution, and interparticle interactions. Finally, the quantitative relations that we established between the blocking temperature, additional energy barrier to magnetization reversal, and nanoparticle concentration are significant for the current efforts to design functional nanomaterials using a rational approach instead of the common trial and error method.
Further work will be aimed at investigating the interplay between the Néel and the Brown superspin relaxation mechanisms at temperatures above the melting point of the carrier fluid. Such studies are important because they enhance the knowledge needed to advance magnetic nanoparticle hyperthermia from its current adjunct-to-other-treatments status to a powerful stand-alone cancer therapy. Our experimental strategy will be to increase the average size of the nanoparticle ensembles so that the superparamagnetic unblocking occurs at temperatures above the freezing point of the carrier fluid. In addition, we will use ensembles of different concentrations to vary the interparticle dipolar interactions, and carrier fluids of different viscosities to vary the energy barrier to physical rotation that governs the Brown relaxation. These experiments are currently underway.