Next Article in Journal
Flying Real-Time Network to Coordinate Disaster Relief Activities in Urban Areas
Next Article in Special Issue
A Multi-Region Magnetoimpedance-Based Bio-Analytical System for Ultrasensitive Simultaneous Determination of Cardiac Biomarkers Myoglobin and C-Reactive Protein
Previous Article in Journal
In Situ Determination of Bisphenol A in Beverage Using a Molybdenum Selenide/Reduced Graphene Oxide Nanoparticle Composite Modified Glassy Carbon Electrode
Previous Article in Special Issue
Aptamer-Modified Magnetic Beads in Biosensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetorelaxometry in the Presence of a DC Bias Field of Ferromagnetic Nanoparticles Bearing a Viscoelastic Corona

Institute of Continuous Media Mechanics, Russian Academy of Sciences, Ural Branch, Perm 614013, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2018, 18(5), 1661; https://doi.org/10.3390/s18051661
Submission received: 16 April 2018 / Revised: 16 May 2018 / Accepted: 18 May 2018 / Published: 22 May 2018
(This article belongs to the Special Issue Magnetic Materials Based Biosensors)

Abstract

:
With allowance for orientational Brownian motion, the magnetorelaxometry (MRX) signal, i.e., the decay of magnetization generated by an ensemble of ferromagnet nanoparticles, each of which bears a macromolecular corona (a loose layer of polymer gel) is studied. The rheology of corona is modelled by the Jeffreys scheme. The latter, although comprising only three phenomenological parameters, enables one to describe a wide spectrum of viscoelastic media: from linearly viscous liquids to weakly-fluent gels. The “transverse” configuration of MRX is considered where the system is subjected to a DC (constant bias) field, whereas the probing field is applied perpendicularly to the bias one. The analysis shows that the rate of magnetization decay strongly depends on the state of corona and slows down with enhancement of the corona elasticity. In addition, for the case of “transverse” MRX, we consider the integral time, i.e., the characteristic that is applicable to relaxation processes with an arbitrary number of decay modes. Expressions for the dependence of the integral time on the corona elasticity parameter and temperature are derived.

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 H to H + H 1 —the same occurs to the system magnetization. Under those conditions, the MRX signal becomes a function of both values ( H and H 1 ), 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 (∼ H 1 H ) 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 η N —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 η M (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, η M η N .
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 ( η N ), 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 η M . 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 η N .
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: τ M = η M / G . In the theoretical consideration below, we analyze how the presence and magnitude of τ M 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 e = μ / μ 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]:
e ˙ = Ω × e , Ω = 1 ζ N L ^ U + Q + y N ( t ) , L ^ = e × e ,
1 + τ M t Q = ζ M Ω + y M ( t ) .
Here, Ω is angular velocity of vector e , and U is the orientation-dependent part of the particle energy, while L ^ is the operator of infinitesimal rotation with respect to e . Vector Q 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
ζ α = 6 η α V , K = 6 G V , α = N , M ,
where V is the particle volume (we assume a sphere), and the subscript indicates either Newtonian ( α = N ) or Maxwell ( α = M ) viscosity (see Figure 1). With Notations (2), the stress relaxation time may be equivalently written as τ M = ζ M / K . 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 F that is to be inserted in Formulas (2) alongside η α and G, respectively. The dependence of F 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:
y i α ( t ) y j β ( t + τ ) = 2 T ζ α δ α β δ i j δ ( τ ) , α = N , M ;
note that, hereafter, we scale temperature in energy units.
The kinetic equation for the distribution function W ( e , Q , t ) that corresponds to the set of stochastic Equation (1) is obtained via standard procedure (see [18,24,25]) for example:
2 τ D t W = ( 1 + q 1 ) β T Q L ^ 1 T Q + β T Q W β T Q L ^ 1 T L ^ U + L ^ W .
Here, the following notations are used:
τ D = ζ N / 2 T = 3 η N V / T
for the Debye time of orientational diffusion of a spherical particle in a fluid of viscosity, and
q = ζ M / ζ N = η M / η N , β = K / T
for non-dimensional rheological parameters. The first one defines “maxwellity” of the Jeffreys medium, so that, at q = 0 , 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 H :
U = e H H = μ H e h h , h 2 = 1 .
In non-dimensional form, the reference magnitude of that energy is
ξ = μ H / T .
It is easy to verify that the equilibrium solution of Equation (4) is given by an extended Boltzmann distribution
W 0 ( e , Q ) exp 1 T U + Q 2 2 K .
From Equation (8), it follows that the equilibrium state of the system in the absence of field U = 0 is isotropic. Note that, although a constant external field endows the system with uniaxial anisotropy, the phase variables e and Q 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
M = n μ e , χ i j = M i / H j ,
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 Δ M ( H , t ) , where Δ M is the projection of magnetization on the direction of the probing field H 1 , 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 δ e H = e e 0 and P = Q × e H ; 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]:
χ ˜ ( ω , β ) = χ ( ω , β ) χ 0 = 1 1 i ω τ D 1 + 1 2 β 1 + 1 2 β / q i ω τ D ,
where χ 0 = n μ 2 / 3 T 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 β 1 , the viscoelastic properties of the corona do not manifest themselves, and the susceptibility reduces to a plain Debye formula with the reference relaxation time τ D . For a medium with high elasticity ( β 1 ), 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
χ ˜ ( ω , β ) 1 2 / β 1 i ω τ s + 2 / β 1 i ω τ f , τ s = τ D 1 + 1 2 β 1 + β / 2 q , τ f = 2 τ D β , β 1 .
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 H 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:
χ ˜ i j ( ω , β , ξ ) = χ ˜ δ i j + ( χ ˜ χ ˜ ) h i h j .
For the case of perturbation induced by a weak probing field H 1 (in non-dimensional form, ξ 1 ), the set of moment equations in the effective field approximation takes the form
t + γ α ( e ) δ e t α = 1 2 ξ 1 f α + 1 2 β P ,
t + 1 2 β ( 1 + q 1 ) + γ α ( M ) P = 1 2 ξ 1 f α + γ α ( e ) δ e t α ,
with the coefficients
γ ( e ) = γ ( M ) = ( 1 + c 2 ) / 2 ( 1 c 2 ) , γ ( e ) = ( 1 c 2 ) / 2 ( c 2 c 1 2 ) ,
γ ( M ) = ( 3 c 2 ) / 2 ( 1 + c 2 ) , f = 1 2 ( 1 + c 2 ) , f = 1 c 2 ,
defined in terms of functions
c k ( ξ ) = 1 1 x k exp ( ξ x ) d x / 1 1 exp ( ξ x ) d x .
On solving Equations (12), one obtains the sought for dynamic susceptibility
χ ˜ α ( ω , β , ξ ) = χ α ( ω , β , ξ ) χ 0 = 1 1 i ω τ D γ α ( e ) 1 + 1 2 β γ α ( M ) + β / 2 q i ω τ D ,
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
χ ˜ α ( ω , β , ξ ) = A s ( α ) ( β , ξ ) 1 i ω τ s ( α ) + A f ( α ) ( β , ξ ) 1 i ω τ f ( α ) , τ s , f = τ D λ s , f ( α ) ,
where the relaxation modes decrements λ s , f ( α ) are the roots of characteristic equation
λ 2 Γ α + γ α ( e ) λ + γ α ( e ) ( Γ α 1 2 β ) = 0 ,
with Γ α = γ α ( M ) + 1 2 β 1 + q 1 . The explicit expressions for the mode amplitudes and decrements are
A s , f ( α ) ( β , ξ ) = 1 2 1 ± Γ α γ α ( e ) Γ α γ α ( e ) 2 + 2 β γ α ( e ) , λ s , f ( α ) = 1 2 Γ α + γ α ( e ) Γ α γ α ( e ) 2 + 2 β γ α ( e ) .
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 λ f ( ) undergoes rapid increase, i.e., the relaxation time goes down. As a result, the spectrum as a whole shifts to the higher frequency range.

4. Magnetorelaxometry

4.1. Relaxation Functions

Consider an MRX experiment that is performed as follows. On the particle ensemble that dwells in equilibrium under a constant bias field H , a weak perturbing field H 1 is imposed transversely to H . This induces the magnetization component Δ M = n μ δ e t in the direction of H 1 . As soon as the magnetization attains its stationary value in the newly established equilibrium, the field H 1 is turned off, and the process of relaxation of Δ M is recorded. This “transverse” configuration seems to be more preferable than measuring the longitudinal relaxation by varying H in small steps because the greater the relative magnitude of the relaxing perturbations (transverse vs. longitudinal), the stronger the bias.
According to the linear response theory (see [24], for example), the reaction to a stepwise turning off the field, is obtained from the dynamic susceptibility (15) by way of the integral relation
χ ˜ ( ω , β , ξ ) = 1 + i ω 0 e i ω t R ( t ) d t ,
where R ( t ) = δ e / δ e 0 is the relaxation function, i.e., the transverse dynamic magnetization normalized by its initial value.
Subjecting Formula (19) to inverse Fourier transformation, one arrives at the expression for relaxation function
R ( t ) = 2 π 0 Im χ ˜ ( ω , β , ξ ) ω cos ω t d ω .
As mentioned above, we take the stress relaxation time τ M as the main characteristic of viscoelasticity of the macromolecular corona. If to describe the corona substance as a Jeffreys fluid, this time parameter is expressed in terms of the rheological coefficients as τ M = ζ M / K = η M / G . Assuming that formation of the corona does not affect the viscosity η N of the solvent, one finds that, for presenting the relaxation function R ( t ) , it is convenient to scale the time in units of τ D = 3 η N V / T . Moreover, the same quantity is appropriate as well for scaling the stress relaxation time τ M . On that basis, we introduce the non-dimensional parameter
κ = τ M / τ D = 2 ζ M T / ζ N K = 2 q / β .
As seen, κ grows with the increase of “maxwellity” of the corona and/or of temperature and goes down with enhancement of the corona elasticity.
The plots illustrating the behavior of relaxation function R ( t ) under variation of bias field H are presented in Figure 2. Their comparison reveals two main specific features. On the one hand, the more viscoelastic the corona, the slower the magnetic relaxation. On the other hand, the bias field, enhancing the effective restoring torque that acts on the particle, accelerates its orientational relaxation. Besides that, Figure 2 proves that magnetic relaxation, as soon as viscoelasticity (non-zero β ) is introduced in the particle environment, is multi-mode, and formally the number of modes makes an infinite countable set. Meanwhile, in our model, this set is replaced just by two modes, slow and fast (see Equation (11)). However, as the detailed multi-mode analysis shows (would be presented elsewhere), this truncation does not affect the above-obtained essential qualitative conclusion: the considered relaxation process has two main reference time scales whose ranges differ by orders of magnitude. In logarithmic representation of Figure 2, the plot of R ( t ) could be schematized as two quasi-straight lines—greater slope at short times and smaller slope at long times—connected by a crossover part. As seen from Figure 2—see curves 1 in all the panes—under zero bias, the crossover part is but weakly distinguishable. However, by using an appropriate bias field, the crossover part of R ( t ) function could be enhanced and exposed to observation (see curves 4 in all the panes).

4.2. Integral Relaxation Time

For studying of multi-mode processes, a useful and experimentally easily determined characteristic is the so-called integral relaxation time that is defined as the area under the R ( t ) curve:
τ α int = 0 R α ( t ) d t .
In the linear response approximation that we use here, the integral time expressed in terms of the dynamic susceptibility is
τ α int = lim ω 0 Im χ ˜ α ( ω , β , ξ ) ω .
With Function (16), integral (23) is taken easily and yields
τ α int = τ D Γ α γ α ( e ) ( Γ α 1 2 β ) = τ D γ α ( M ) + 1 2 β ( 1 + q 1 ) γ α ( e ) γ α ( M ) + β / 2 q .
This formula (with accuracy of the effective field approximation) in a qualitatively correct way describes the dependence of magnetization relaxation on all the relevant material parameters of the system. For example, for an isotropic system ( ξ = 0 ) assuming that β > 1 and going back to dimensional time, one gets from Equation (23):
τ int = τ D 1 + 1 2 β ( 1 + q 1 ) 1 + β / 2 q τ D 1 + 1 2 β , β 2 q , 1 + q , β 2 q .
As the second line of Equation (25) shows, the asymptotic value of the integral time is defined by the slow diffusion. However, to really attain this regime, one needs a system with very strong elasticity: β 2 q . Even in dense polymer gels and “live polymers” where q 10 3 10 5 , one could hardly expect to match this condition. Under this limitation, one may always assume that the dependence of integral time (23) on the elasticity parameter is by and large linear.
In the presence of a bias field, full Expressions (13) for coefficients γ α ( e ) and γ α ( M ) should be used. Substitution of those in Equation (24) yields the formulas that enable one to estimate the effect of the bias field on the integral times of transverse and longitudinal relaxations of the dynamic magnetization:
τ int = 2 τ D 1 c 2 1 + c 2 1 + β 1 + c 2 3 c 2 , τ int = 2 τ D c 2 c 1 2 1 c 2 1 + β 1 c 2 1 + c 2 , β 2 q .
For a strong field ( ξ 1 ), Equation (26) expands to simple asymptotic expressions
τ int = τ H 1 + β 2 β ξ , τ int = 1 2 τ H 1 + 1 + β ξ , τ H ζ N μ H = 2 τ D ξ .
Therefore, as it follows from Equation (27), under a strong bias, the initial stage of relaxation is mostly due to the forced field-induced rotation of the particles with reference time τ int τ H 2 τ D / ξ .
The overall behavior of the integral time is illustrated in Figure 3, where its dependence on the bias magnitude is given for different values of parameter κ = τ M / τ D . The values of integral time under zero bias obey very well Formula (25) for β q : the integral time grows linearly with the elasticity parameter β = K / T . As already mentioned, this increase might halt only at β 2 q , where τ int ( q + 1 ) τ D , but that limit is practically unattainable. For the parameters used for plotting Figure 3, this upper limit would have ranged τ int / τ D 100 .

5. Conclusions

To summarise the results of the work, three essential points should be mentioned. First, the macromolecular corona covering the particles is considered as a viscoelastic entity. This approximation is qualitatively different from a customary approach where the effect of corona is reduced to just a change of the particle hydrodynamic diameter. Second, we study MRX that is performed in the presence of a bias field, so that the bias field strength is a controllable parameter of the experiment. Third, the main attention is focused on the MRX measured in the “transverse” geometry where the probing field is applied normally to the bias one. The advantage of that variant comes from a significant enhancement of the response in comparison with that obtained in “longitudinal” configuration. Therefore, the transverse configuration enables one to preserve the sensitivity of the method despite that in general the bias field diminishes the particle response to orientational perturbations. We are well aware that the considerations presented here are to a certain extent illustrative as they are obtained with a severely simplified kinetic description: the whole set of moment equations has been replaced by just a pair of them (effective field model). However, our goal is to show, in the first place, the qualitatively relevant features of the considered problem. The work on a complete solution of the kinetic equation is under way.

Author Contributions

V.V. and Y.R. have done all the work on the paper together: from conceiving the idea to writing of the final text.

Funding

The work was done under the auspices of RFBR grant 17-42-590504.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sanson, C.; Diou, O.; Thévenot, J.; Ibarboure, E.; Soum, A.; Brûlet, A.; Miraux, S.; Thiaudière, E.; Tan, S.; Brisson, A.; et al. Doxorubicin loaded magnetic polymersomes: Theranostic nanocarriers for MR imaging and magneto-chemotherapy. ACS Nano 2011, 5, 1122–1140. [Google Scholar] [CrossRef] [PubMed]
  2. Oliveira, H.; Pérez-Andrés, E.; Thevenot, J.; Sandre, O.; Berra, E.; Lecommandoux, S. Magnetic field triggered drug release from polymersomes for cancer therapeutics. J. Control. Release 2013, 169, 165–170. [Google Scholar] [CrossRef] [PubMed]
  3. Li, Y.; Huang, G.; Zhang, X.; Li, B.; Chen, Y.; Lu, T.; Lu, T.J.; Xu, F. Magnetic hydrogels and their potential biomedical applications. Adv. Funct. Mater. 2013, 23, 660–672. [Google Scholar] [CrossRef]
  4. Pankhurst, Q.A.; Thanh, N.K.T.; Jones, S.K.; Dobson, J. Progress in applications of magnetic nanoparticles in biomedicine. J. Phys. D Appl. Phys. 2009, 42, 224001. [Google Scholar] [CrossRef]
  5. Dutz, S.; Hergt, R. Magnetic nanoparticle heating and heat transfer on a microscale: Basic principles, realities and physical limitations of hyperthermia for tumour therapy. Int. J. Hyperth. 2013, 29, 790–800. [Google Scholar] [CrossRef] [PubMed]
  6. Sanchez, C.; El Hajj Diab, D.; Connord, V.; Clerc, P.; Meunier, E.; Pipy, B.; Payré, B.; Tan, R.P.; Gougeon, M.; Carrey, J.; et al. Targeting a G-protein-coupled receptor overexpressed in endocrine tumors by magnetic nanoparticles to induce cell death. ACS Nano 2014, 8, 1350–1363. [Google Scholar] [CrossRef] [PubMed]
  7. Foroozandeh, P.; Aziz, A.A. Merging worlds of nanomaterials and biological environment: Factors governing protein corona formation on nanoparticles and its biological consequences. Nanoscale Res. Lett. 2015, 10, 221. [Google Scholar] [CrossRef] [PubMed]
  8. Eberbeck, D.; Bergemann, C.; Wiekhorst, F.; Steinhoff, U.; Trahms, L. Quantification of specific bindings of biomolecules by magnetorelaxometry. J. Nanobiotechnol. 2008, 6. [Google Scholar] [CrossRef] [PubMed]
  9. Wiekhorst, F.; Steinhoff, U.; Eberbeck, D.; Trahms, L. Magnetorelaxometry assisting biomedical applications of magnetic nanoparticles. Pharm. Res. 2012, 29, 1189–1202. [Google Scholar] [CrossRef] [PubMed]
  10. Liebl, M.; Wiekhorst, F.; Eberbeck, D.; Radon, P.; Gutkelch, D.; Baumgartner, D.; Steinhoff, U.; Trahms, L. Magnetorelaxometry procedures for quantitative imaging and characterization of magnetic nanoparticles in biomedical applications. Biomed. Eng. 2015, 60, 427–443. [Google Scholar] [CrossRef] [PubMed]
  11. Shliomis, M.I.; Stepanov, V.I. Theory of the dynamic susceptibility of magnetic fluids. Adv. Chem. Phys. 1994, 78, 1–30. [Google Scholar]
  12. Raikher, Y.L.; Shliomis, M.I. Effective field method for orientational kinetics of magnetic fluids and liquid crystals. Adv. Chem. Phys. 1994, 78, 595–751. [Google Scholar]
  13. Fortin, J.-P.; Wilhelm, C.; Servais, J.; Ménager, C.; Bacri, J.-C.; Gazeau, F. Size-sorted anionic iron oxide nanomagnets as colloidal mediators for magnetic hyperthermia. J. Am. Chem. Soc. 2007, 129, 2628–2635. [Google Scholar] [CrossRef] [PubMed]
  14. Raikher, Y.L.; Stepanov, V.I. Physical aspects of magnetic hyperthermia: Low-frequency ac field absorption in a magnetic colloid. J. Magn. Magn. Mater. 2014, 368, 421–427. [Google Scholar] [CrossRef]
  15. Oswald, P. Rheophysics: The Deformation and Flow of Matter; Cambridge University Press: Cambridge, UK, 2009; ISBN 978-0-521-88362-7. [Google Scholar]
  16. Raikher, Y.L.; Rusakov, V.V. Brownian motion in a Jeffreys fluid. J. Exp. Theor. Phys. 2010, 111, 883–889. [Google Scholar] [CrossRef]
  17. Raikher, Y.L.; Rusakov, V.V.; Perzynski, R. Brownian motion in a viscoelastic medium modelled by a Jeffreys fluid. Soft Matter 2013, 9, 10857–10865. [Google Scholar] [CrossRef]
  18. Raikher, Y.L.; Rusakov, V.V.; Perzynski, R. Brownian motion in the fluids with complex rheology. Math. Model. Nat. Phenom. 2015, 10, 1–43. [Google Scholar] [CrossRef]
  19. Ochab-Marcinek, A.; Wieczorek, S.A.; Ziębacz, N.; Hołyst, R. The effect of depletion layer on diffusion of nanoparticles in solutions of flexible and polydisperse polymers. Soft Matter 2012, 8, 11173–11179. [Google Scholar] [CrossRef]
  20. Gomez-Solano, J.R.; Bechinger, C. Transient dynamics of a colloidal particle driven through a viscoelastic fluid. New J. Phys. 2015, 17, 103032. [Google Scholar] [CrossRef]
  21. Pavlovsky, L.; Younger, J.G.; Solomon, M.J. In situ rheology of Staphylococcus epidermidis bacterial biofilms. Soft Matter 2013, 9, 122–132. [Google Scholar] [CrossRef] [PubMed]
  22. Gardel, M.L.; Valentine, M.T.; Weitz, D.A. Microrheology. In Microscale Diagnostic Techniques; Breuer, K., Ed.; Springer: New York, NY, USA, 2005; ISBN 3-540-23099-8. [Google Scholar]
  23. Happel, J.; Brenner, H. Low Reynolds Number Hydrodynamics; Prentice-Hall: Boston, MA, USA, 1965. [Google Scholar]
  24. Coffey, W.T.; Kalmykov, Y.P. The Langevin Equation, 4rd ed.; World Scientific: Singapore, 2017; ISBN 978-981-3221-99-4. [Google Scholar]
  25. Ilg, P.; Evangelopoulos, E.A.S. Magnetic susceptibility, nanorheology, and magnetoviscosity of magnetic nanoparticles in viscoelastic environments. Phys. Rev. E 2018, 97, 032610. [Google Scholar] [CrossRef] [PubMed]
  26. Rusakov, V.V.; Raikher, Y.L. Magnetic response of a viscoelastic ferrodispertion: From a nearly Newtonian ferrofluid to a Jeffreys ferrogel. J. Chem. Phys. 2017, 147, 124903. [Google Scholar] [CrossRef] [PubMed]
  27. Raikher, Y.L.; Stepanov, V.I. Nonlinear dynamic susceptibilities and field-induced birefringence in magnetic particle assemblies. Adv. Chem. Phys. 2004, 129, 419–588. [Google Scholar]
  28. Raikher, Y.L.; Rusakov, V.V. Orientational Brownian motion in a viscoelastic medium. Colloid J. 2017, 79, 264–269. [Google Scholar] [CrossRef]
Figure 1. Sketch of a magnetic particle bearing a corona of Jeffreys medium. The system is set under constant field H , and the existing equilibrium is perturbed by the field H 1 switched on/off stepwise or alternating harmonically.
Figure 1. Sketch of a magnetic particle bearing a corona of Jeffreys medium. The system is set under constant field H , and the existing equilibrium is perturbed by the field H 1 switched on/off stepwise or alternating harmonically.
Sensors 18 01661 g001
Figure 2. Relaxation functions for the “transverse” MRX; parameter κ = 2 × 10 3 (a); 100 (b); 40 (c); 20 (d); the non-dimensional bias field is ξ = 0 . 1 (1), 2 (2), 5 (3), 10 (4).
Figure 2. Relaxation functions for the “transverse” MRX; parameter κ = 2 × 10 3 (a); 100 (b); 40 (c); 20 (d); the non-dimensional bias field is ξ = 0 . 1 (1), 2 (2), 5 (3), 10 (4).
Sensors 18 01661 g002
Figure 3. Dependence of the integral relaxation time on the strength of bias field for MRX in “transverse” configuration; the viscoelasticity parameter κ = 2 × 10 3 (1), 100 (2), 40 (3), 20 (4); “maxwellity” of the corona is q = 100.
Figure 3. Dependence of the integral relaxation time on the strength of bias field for MRX in “transverse” configuration; the viscoelasticity parameter κ = 2 × 10 3 (1), 100 (2), 40 (3), 20 (4); “maxwellity” of the corona is q = 100.
Sensors 18 01661 g003

Share and Cite

MDPI and ACS Style

Rusakov, V.; Raikher, Y. Magnetorelaxometry in the Presence of a DC Bias Field of Ferromagnetic Nanoparticles Bearing a Viscoelastic Corona. Sensors 2018, 18, 1661. https://doi.org/10.3390/s18051661

AMA Style

Rusakov V, Raikher Y. Magnetorelaxometry in the Presence of a DC Bias Field of Ferromagnetic Nanoparticles Bearing a Viscoelastic Corona. Sensors. 2018; 18(5):1661. https://doi.org/10.3390/s18051661

Chicago/Turabian Style

Rusakov, Victor, and Yuriy Raikher. 2018. "Magnetorelaxometry in the Presence of a DC Bias Field of Ferromagnetic Nanoparticles Bearing a Viscoelastic Corona" Sensors 18, no. 5: 1661. https://doi.org/10.3390/s18051661

APA Style

Rusakov, V., & Raikher, Y. (2018). Magnetorelaxometry in the Presence of a DC Bias Field of Ferromagnetic Nanoparticles Bearing a Viscoelastic Corona. Sensors, 18(5), 1661. https://doi.org/10.3390/s18051661

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop