1. Introduction: Magnetorelaxometry in Linear Approximation
Magnetic nanoparticles are successfully used in diverse bioengineering and medical applications both as themselves (e.g., as contrasting agents in Magnetic Resonance Imaging (MRI)) or as essential components of physicochemical complexes (magnetic polymerosomes [
1,
2], microferrogels [
3]). These techniques span from object-oriented drug delivery [
4] to active thermal action on malignant cells (magneto-inductive hyperthermia [
5]) or forced penetration through cell membranes [
6].
Nowadays, one of the most developed applications of magnetic nanoparticles is their introduction as sensors for diagnosing the state and content of complex media both of non-organic and biological origin. In the context of conventional microrheology approach, it is assumed that the particles do not chemically react with the medium under study and, thus, serve as the means of specific “nondestructive testing”. On the contrary, in biochemical analysis, the nanoparticles are functionalized, and the main interest is focused on the degree of reaction between the molecules of tested solution with the markers grafted to the particles or with their bare but chemically active surface. The macromolecular coating that a particle acquires as a result of this adsorption is often termed as a protein corona [
7].
A unique advantage of magnetic nanosensors is that, with the aid of an applied magnetic field, one is able to remotely excite their motion and analyze the generated response to that.
The chemical processes resulting in formation of the corona on the particle surface affect its dynamic behavior. Those changes could be sensed by either measuring the magnetic spectrum of the system (if the probing is done with an AC field) or—in pulse regime—by registering the signal of the magnetization decay. The latter experimental technique is known as magnetorelaxometry (MRX). In many cases, MRX as a laboratory test is more preferable in comparison with magnetic spectroscopy, in particular, due to its much easier technical implementation [
8,
9,
10].
Consider a ferromagnet nanoparticle that, due to its smallness, is single-domain and, thus, bears a magnetic moment of constant magnitude. This particle floats in a water solution of macromolecules some of which, possessing affinity to the particle surface, adsorb on it, so that the particle environment transforms from a Newtonian liquid into a viscoelastic polymer gel.
For simplicity of the following considerations, we assume that vector
is “frozen” in the particle body and treat the nanoparticle as a miniature permanent magnet. This means neglecting the superparamagnetic effect: thermal fluctuations of
inside the particle. Such an approach implies that the Néel (internal) relaxation time of
is much greater than any other reference time of the problem. Certainly, this approximation is not universally valid. However, as it follows from a number of works on the subject, see, for example, reviews [
11,
12] for general considerations and papers [
13,
14] for realistic examples. The particles of 15–20 nm in diameter made of a moderately magnetically hard substance (e.g., cobalt ferrite) matches this requirement fairly well in the frequency range below 1 MHz, which is most relevant for MRX.
For the above-described magnetically hard nanoparticle suspended in a fluid medium, the only fluctuational process affecting vector is its rotary Brownian motion together with the particle. Due to that, in the absence of field, the directions of magnetic moments are distributed at random, and the net magnetization of the ensemble is zero. An MRX measurement begins with application of a uniform magnetic field that orients the particle magnetic moments, thus inducing a non-zero equilibrium magnetization of the system. Then, the magnitude (or direction) of the field is abruptly changed, and the signal generated by the magnetization evolving to a new equilibrium is registered for further analysis. By that, MRX not only accomplishes its general purpose—to evaluate the amount of particles (that is proportional to the signal intensity)—but delivers information on the details of the particle rotary motion.
If the exerted field is turned off completely, the magnetization of the particle ensemble decays freely down to zero, and this process is governed solely by the rotary Brownian diffusion. There are no external factors. If the field is switched between two finite values—from to —the same occurs to the system magnetization. Under those conditions, the MRX signal becomes a function of both values ( and ), by variation of which one can extract more information from the same experiment than in the case where the field is just turned off. A convenient way to study this relaxation is to change the bias field in small (∼) steps to be able to use the linear response theory. In that case, the obtained MRX signal depends parametrically only on that value H of the field, around which the variation is done. The “linear” MRX of that kind is the subject of the present work.
2. Model
We consider a nanoparticle that, as a result of formation of a corona, dwells in a viscoelastic environment, i.e., a medium with retarded response. As a rheological model for the latter, the Jeffreys scheme (see [
15], for example) is used because, unlike the plain Maxwell one, it is robust when applied to Brownian motion and is free of artifacts [
16,
17,
18]. The viscoelastic properties of the Jeffreys model are fully rendered by three parameters (see the scheme outlined in
Figure 1). A single-element chain (narrow damper) associated with the solvent there—a Newtonian fluid with the viscosity coefficient
—is set in parallel with a two-element Maxwell chain associated with the corona. The latter is assumed to possess both elasticity with modulus
G (the spring) and intrinsic viscosity
(wide damper). Evidently, in a medium with a pronounced visoelasticty, the Jeffreys viscosity coefficients are substantially different: the intrinsic viscosity of the macromolecular gel is much greater than that of the low-molecular solvent,
.
To justify application of the Jeffreys scheme to the rotary motion of the particles (
Figure 1), we remind readers that this phenomenology works rather well for the translational Brownian motion in semi-dilute polymer solutions. In such models [
19,
20], it is assumed that each particle is surrounded by a depletion layer where concentration of macromolecules is much lower than that in the bulk. Inside this thin layer, the viscosity is effectively small (
), and the diffusion process is fast (short times, small distances). However, when larger displacements (of the order of the particle size or greater) are considered, the particle experiences its environment as a medium with a high viscosity
. Similar behavior is inherent to biofilms [
21], which, at short time intervals, respond as low viscous fluids but at longer scale react as weakly-fluent gels.
Following this line, we infer that a scheme with two rather different viscosity coefficients should be appropriate for the rotary diffusion of the particles furnished with a corona. Indeed, a loose macromolecular coating, just slightly changing the total mass of the system, at the same time strongly affects its rotary friction. Moreover, the presence of corona cannot be accounted for by just a simple renormalization of the particle hydrodynamic diameter since the corona brings in the particle dynamics a substantial retardation component. That is, if the particle turns over large angles, this motion entrains the whole corona. On the other hand, the corona is virtually insensitive to small angle displacements because such a motion concerns only the molecular fragments in close vicinity of the particle surface, and, due to that, the contribution of corona to resistance and retardation is insignificant. Therefore, at short time intervals, the particle may be considered as floating in a Newtonian fluid with low viscosity .
For nanoparticles embedded in any viscous environment, the effect of inertia is negligible, and the process of rotary relaxation is always monotonic (overdamped regime). Given that, we chose as a main indicator of the state of corona the relaxation time of its elastic stress: . In the theoretical consideration below, we analyze how the presence and magnitude of is reflected in the MRX spectra of the particle ensemble subjected to a bias DC field.
As in the adopted model, the magnetic moment
is “frozen” in the particle, and we take unit vector
as a marker of the particle orientation. The equations of rotary motion for a Brownian particle in the inertialess limit take the form [
17,
18]:
Here,
is angular velocity of vector
, and
U is the orientation-dependent part of the particle energy, while
is the operator of infinitesimal rotation with respect to
. Vector
in Equation (
1) has the meaning of a torque acting on the particle on the part of the corona. The response coefficients of the viscoelastic medium are defined in a standard way [
22] as
where
V is the particle volume (we assume a sphere), and the subscript indicates either Newtonian (
) or Maxwell (
) viscosity (see
Figure 1). With Notations (
2), the stress relaxation time may be equivalently written as
. We remark that there are no universal expressions to replace Formulas (
2) provided the particles are non-spherical (anisometric). However, quite reliable estimates could be obtained if to approximate a particle with an ellipsoid of revolution (spheroid). For that shape, the effect of non-sphericity is rendered by a formfactor
that is to be inserted in Formulas (
2) alongside
and
G, respectively. The dependence of
as a function of the particle aspect ratio is known (see [
23]) for an example.
Correlators of the random forces that model thermal noise in the system are expressed with the aid of the fluctuation–dissipation theorem:
note that, hereafter, we scale temperature in energy units.
The kinetic equation for the distribution function
that corresponds to the set of stochastic Equation (
1) is obtained via standard procedure (see [
18,
24,
25]) for example:
Here, the following notations are used:
for the Debye time of orientational diffusion of a spherical particle in a fluid of viscosity, and
for non-dimensional rheological parameters. The first one defines “maxwellity” of the Jeffreys medium, so that, at
, the model reduces to an ordinary viscous fluid; the second parameter is the dynamic elasticity coefficient scaled with thermal energy.
For a magnetically hard particle, the energy
U that enters Equation (
4) reduces to the Zeeman interaction with uniform field
:
In non-dimensional form, the reference magnitude of that energy is
It is easy to verify that the equilibrium solution of Equation (
4) is given by an extended Boltzmann distribution
From Equation (
8), it follows that the equilibrium state of the system in the absence of field
is isotropic. Note that, although a constant external field endows the system with uniaxial anisotropy, the phase variables
and
remain statistically independent in that case as well.
In a dilute (the interparticle interaction is neglected) statistical ensemble of such particles, the magnetization and susceptibility are rendered by expressions
where
n is the number concentration of particles, whereas angular brackets denote averaging with the distribution function from Equation (
4).
3. Dynamic Susceptibility
The MRX relaxation function, i.e., the dependence
, where
is the projection of magnetization on the direction of the probing field
, and
H is the bias field strength (see
Figure 1), could be derived by several ways. We chose the one where the relaxation function is obtained in terms of the dynamic magnetic susceptibility
. The reason is two-fold. First, in the linear response approximation, the relation between the relaxation function and
has a simple form. Second, in Refs. [
18,
26], we have developed a workable technique to obtain the pertinent dynamic susceptibility for the case of zero external field. Using those results, one may skip a considerable part of lengthy calculations and just add only the extension allowing for the bias field. Besides that, in [
18,
26], while analyzing the results obtained from solving the kinetic equation (the analogue of Equation (
4)) in an exact way, i.e., using a long multi-moment expansion, we have shown that, to get a plausibly accurate approximation, it suffices to truncate the infinite moment set to just the first two equations for the dynamic variables
and
; here, angular brackets with index 0 denote averaging over the equilibrium distribution (
8).
The aforementioned approximation is well known as the effective field model, and it has proven its usefulness for diverse problems of orientational kinetics of nanoparticles many times [
12,
27,
28]. For the case of zero bias, the dynamic magnetic susceptibility in the effective field approximation is [
28]:
where
is the static magnetic susceptibility of an ensemble of noninteracting particles bearing magnetic moments
.
As seen from Equation (
10), under zero bias the dynamic susceptibility is isotropic, and the type of its frequency dependence is determined by the elasticity
and “maxwellity”
q of the medium (corona). At high temperatures
, the viscoelastic properties of the corona do not manifest themselves, and the susceptibility reduces to a plain Debye formula with the reference relaxation time
. For a medium with high elasticity (
), the dependence (
10) has two maxima and may be with plausible accuracy presented as a superposition of two relaxation modes: slow (
s) and fast (
f) in the form
Equation (
11) shows that, under enhancement of elasticity, the slow component grows, and its peak moves further to the low-frequency domain, which means gradual increase of its reference time. Concurrently, the contribution of the fast relaxation mode goes down monotonically with
.
In the presence of bias field
that makes the system uniaxially anisotropic, the susceptibility becomes a second-rank tensor that is diagonal in the coordinate frame whose
Oz axis points along the bias field:
For the case of perturbation induced by a weak probing field
(in non-dimensional form,
), the set of moment equations in the effective field approximation takes the form
with the coefficients
defined in terms of functions
On solving Equations (
12), one obtains the sought for dynamic susceptibility
where index
assumes one of two possible values: ⊥ or ∥.
Presenting Function (
15) as a sum of two relaxation modes—this time, for arbitrary values of parameters
and
and—one gets
where the relaxation modes decrements
are the roots of characteristic equation
with
. The explicit expressions for the mode amplitudes and decrements are
Let us consider the frequency dependence of the transverse (⊥) component of susceptibility (
16). Under a weak constant field, it only slightly differs from the one rendered by Formula (
10). However, in a medium with a pronounced elasticity, the fast mode might become dominating. Indeed, under a growing bias field, the decrement of the slow mode virtually does not change while
undergoes rapid increase, i.e., the relaxation time goes down. As a result, the spectrum as a whole shifts to the higher frequency range.