Properties of assembly of superparamagnetic nanoparticles in viscous liquid

Detailed calculations of the specific absorption rate (SAR) of a dilute assembly of iron oxide nanoparticles with effective uniaxial anisotropy dispersed in a liquid are performed depending on the particle diameters, the alternating (ac) magnetic field amplitude H0 and the liquid viscosity. For small and moderate H0 values with respect to particle anisotropy field Hk the SAR of the assembly as a function of the particle diameter passes through a characteristic maximum and then reaches a plateau, whereas for sufficiently large amplitudes, H0 ~ Hk, the SAR increases monotonically as a function of diameter. The realization of viscous and magnetic oscillation modes for particle unit magnetization vector and director for moderate and sufficiently large H0 values, respectively, explains this behavior. It is found that the SAR of the assembly changes inversely with the viscosity only in a viscous mode, for nanoparticles of sufficiently large diameters. In the magnetic mode the SAR of the assembly is practically independent of the viscosity, since in this case the nanoparticle director only weakly oscillates around the ac magnetic field direction. The conditions for the validity of the linear response theory have been clarified by comparison with the numerical simulation data.

Detailed calculations of the specific absorption rate (SAR) of a dilute assembly of iron oxide nanoparticles with effective uniaxial anisotropy dispersed in a liquid are performed depending on the particle diameters, the alternating (ac) magnetic field amplitude H 0 and the liquid viscosity. For small and moderate H 0 values with respect to particle anisotropy field H k the SAR of the assembly as a function of the particle diameter passes through a characteristic maximum and then reaches a plateau, whereas for sufficiently large amplitudes, H 0 ~ H k , the SAR increases monotonically as a function of diameter. The realization of viscous and magnetic oscillation modes for particle unit magnetization vector and director for moderate and sufficiently large H 0 values, respectively, explains this behavior. It is found that the SAR of the assembly changes inversely with the viscosity only in a viscous mode, for nanoparticles of sufficiently large diameters. In the magnetic mode the SAR of the assembly is practically independent of the viscosity, since in this case the nanoparticle director only weakly oscillates around the ac magnetic field direction. The conditions for the validity of the linear response theory have been clarified by comparison with the numerical simulation data.
The behavior of an assembly of magnetic nanoparticles in a static or alternating (ac) magnetic field is of great interest due to the important applications of these assemblies in biomedicine 1 . The ability of magnetic nanoparticles to effectively absorb the energy of ac magnetic field can be used in magnetic hyperthermia, a new promising method for treating cancer [2][3][4][5] . Generally speaking, for applications in magnetic hyperthermia it is necessary to study the specific absorption rate (SAR) of an assembly in a biological environment. However, in a biological medium there is an intense interaction of nanoparticles with intracellular structures 6,7 . This leads to the formation of clusters of particles 8,9 and to a significant increase in the intensity of magnetic dipole interaction in the assembly 10,11 . In addition, the rheological properties of biological media have not yet been fully studied 12 . All these factors complicate the analysis of experimental results obtained for nanoparticles in a biological medium. Therefore, a common practice is to measure SAR of prepared assemblies in a liquid since it helps to quickly assess the prospects of a particular assembly for use in magnetic hyperthermia.
To date, a large number of SAR measurements have been carried out in liquids of various viscosities 6,7,[13][14][15][16][17][18][19][20][21] . In a number of cases [13][14][15] , very high SAR values, more than 1 kW/g, were obtained at high frequencies, f = 300-700 kHz, and large amplitudes of the ac magnetic field, H 0 = 300-600 Oe. However, for medical reasons one should strive to obtain a sufficient thermal effect at moderate values of these parameters because the product f H 0 must satisfy the Brezovich 22 , or at least Hergt et al. 23 criterion, f H 0 ≤ 6.28 × 10 4 kHz Oe. Besides, the influence of the liquid viscosity, magnetic parameters and characteristic sizes of nanoparticles 13-21 on the SAR of the assembly requires a more detailed study.
An early paper 13 stated that the linear response theory (LRT) developed by Rosensweig 24 describes the experimental SAR measurements in a viscous liquid very satisfactorily. Later, however, significant deviations from the LRT predictions were observed 15,16 with an increase in the ac magnetic field amplitude. It should be noted that theoretical calculations of SAR of dilute assemblies of magnetic nanoparticles in a viscous liquid are complicated by the fact that a magnetic nanoparticle in a liquid can rotate as a whole under the action of fluctuating thermal torques and external ac magnetic field. The spatial orientation of a single-domain nanoparticle with uniaxial magnetic anisotropy can be specified by a unit vector, the so-called director n, which indicates the direction of the easy anisotropy axis. Furthermore, the directions of the particle unit magnetization vector α and the director n are not independent, since the magnetic anisotropy energy of the particle depends on the angle between these vectors. In the LRT 24 a number of simplifying assumptions were introduced that allow one to gracefully bypass www.nature.com/scientificreports/ a detailed theoretical analysis of physical processes leading to the absorption of the energy of ac magnetic field by an assembly of nanoparticles in a liquid. However, the area of applicability of the LRT is still poorly studied. In subsequent works [25][26][27][28][29][30] , various generalizations of LRT were developed, based on the introduction of the expressions for the dynamic magnetic susceptibility or for the relaxation times of superparamagnetic nanoparticles in a viscous liquid. An alternative approach was also developed 31-36 based on the use of a system of stochastic equations to describe the dynamics of the unit magnetization vector and the particle director at a finite temperature. A useful overview of the current state of the problem can be found in Ref. 37 .
To study the behavior of a dilute assembly of magnetic nanoparticles in a viscous liquid, a system of stochastic equations that describes the coupled dynamics of vectors α and n under the action of ac magnetic field and thermal fluctuations in the medium was proposed 31 . Two vibration modes (the so called viscous and magnetic) were shown to be realized for vectors α and n depending on the amplitude H 0 of ac magnetic field compared to the particle anisotropy field H k . In the viscous mode, which is realized at low and moderate amplitudes, H 0 < H k /2, the vectors α and n move approximately in unison. A magnetic mode arises with an increase in the field amplitude, H 0 ~ H k . In magnetic mode the particle director experiences only small oscillations about the magnetic field direction, whereas the magnetic vector periodically jumps between equivalent potential wells with the frequency of the external ac magnetic field.
The orientation of nanoparticle directors in a sufficiently strong ac magnetic field was already observed in the early experiment 14 . The existence of the magnetic oscillation mode has also been confirmed in recent experiments 38,39 . Meanwhile, the results of SAR measurements of nanoparticle assemblies distributed in a liquid are usually compared 13,15,16,[28][29][30] with the LRT predictions, although formally this theory can be valid only in a relatively weak magnetic field. In this work, using the previously developed approach 31 , detailed calculations of the SAR of a dilute assembly of iron oxide nanoparticles are performed depending on the particle diameters, the alternating magnetic field amplitude and the liquid viscosity. The conditions for the applicability of the LRT have been clarified.

Results and discussion
Dilute nanoparticle assembly in viscous liquid. A behavior of a dilute assembly of magnetic nanoparticles in a viscous liquid can be studied 31 by solving a set of stochastic equations that describe the coupled dynamics of the unit magnetization vector α and the director n of a particle in ac magnetic field. A brief description of this approach is given in "Methods" section. Using this technique in this work we investigate in detail the dependence of the SAR of a dilute assembly of quasi-spherical single-domain iron oxide nanoparticles with a saturation magnetization M s = 350 emu/cm 3 and an effective uniaxial anisotropy constant K 1 = 10 5 erg/cm 3 on the particle size in the range of diameters D = 15-55 nm, on the amplitude of the ac magnetic field, H 0 = 50-500 Oe, and the kinematic viscosity of the liquid, η = 0.01-0.2 g/(cm s). The density of magnetic nanoparticles is taken to be ρ = 5.0 g/cm 3 , the temperature of the system T = 300 K. The parameters mentioned are of greatest interest for the application of magnetic nanoparticles in magnetic hyperthermia [2][3][4][5] . The calculations were carried out at magnetic field frequency f = 300 kHz, that is a typical value for SAR measurements in a viscous liquid [13][14][15][16][17][18][19] . The results of the calculations performed are shown in Fig. 1.
In Fig. 1a one notes a significant difference in the dependence of SAR on the particle diameter for the cases of moderate, H 0 = 100, 200 Oe, and sufficiently large amplitudes, H 0 = 300, 500 Oe, of the ac magnetic field. In the first case the SAR as a function of the particle diameter passes through a characteristic maximum and then reaches a plateau, while in the second case it increases monotonically with an increase in the particle diameter.
As noted earlier 31 , this difference is due to the fact that in the range of amplitudes H 0 < H k /2, where H k = 2K 1 /M s = 571 Oe is the particle anisotropy field, the dynamics of the unit magnetization vector α and director n occurs mainly in the so called viscous oscillation mode, where the movement of both vectors is nearly in unison. At the same time, for sufficiently large amplitudes, H 0 ~ H k , a magnetic oscillation mode is realized. In the magnetic mode, in the course of time the particle directors are oriented along the line of action of the ac magnetic field, while the magnetic vector periodically jumps between two equivalent magnetic potential wells. In some cases, there is also an intermediate regime in which the viscous and magnetic motion of the vectors alternate. Figure 2a shows the nature of the motion of vectors α and n in the viscous mode. Here, the black and red curves show the projections of the unit magnetization vector and the director of an individual, randomly selected nanoparticle of the assembly onto the applied magnetic field direction for several periods of field variation. This type of motion of vectors α and n is characteristic when the SAR of the assembly reaches a plateau. In particular, Fig. 2a shows the oscillations of the z-components of these vectors at H 0 = 200 Oe, f = 300 kHz, D = 30 nm, which corresponds to point 1 in Fig. 1a. However, with an increase in the ac field amplitude, H 0 = 300 Oe (point 2 in Fig. 1a), the magnetic mode occurs in the same assembly. As Fig. 2b shows, in the magnetic mode the longitudinal projection of the director is close to n z = + 1 (the value n z = − 1 is equally possible), whereas projection of the unit magnetization vector oscillates between the values α z = ± 1 with the frequency of the ac magnetic field. It should be noted that for particles of large diameters, and in liquid with appreciable viscosity, the establishment of stationary motion in the magnetic mode can occur rather slowly, over several tens of periods of the ac magnetic field.
Finally, Fig. 2c shows the intermediate mode of motion of vectors α and n at point 3 in Fig. 1a. In this regime, the longitudinal projection of the director also periodically jumps between the values n z = ± 1, but the frequency of this jump is much lower than the frequency f of the applied ac magnetic field. Figure 2d shows the low frequency hysteresis loops of a dilute assembly consisting of N p = 400 nanoparticles with a diameter D = 30 nm. These loops correspond to the viscous and magnetic modes of the assembly at points 1 and 2 in Fig. 1a, respectively. It is easy to see that the shape of the hysteresis loop for the magnetic mode is close to that of an ideal rectangular hysteresis loop. Note, that for an ideal rectangular hysteresis loop the SAR of the assembly has the maximum possible value, SAR = 4M s H 0 f/ρ 40  www.nature.com/scientificreports/ It is interesting to note that at point 3 in Fig. 1a the SAR of the assembly peaks as a function of the particle diameter. For the case of H 0 = 200 Oe the SAR at the maximum, at D = 24 nm, is 972.5 W/g, while the characteristic SAR values at the plateau, at D > 30 nm, are only 320-350 W/g. Similarly, for the case of H 0 = 100 Oe, the maximum SAR at D = 20 nm is 283 W/g, while the SAR at the plateau, D > 24 nm is approximately 70 W/g. Consequently, to generate heat in a viscous liquid at moderate amplitudes of ac magnetic field it is preferable to use assemblies of iron oxide nanoparticles with characteristic diameters D = 20-24 nm. Of course, with an increase in the field amplitude up to H 0 = 300 Oe much higher SAR values can be obtained at the same frequency f = 300 kHz, on the order of 1600-1700 W/g, for nanoparticle diameters D > 30 nm. But the use of ac magnetic field of large amplitude is energetically expensive and unsafe in a medical clinic. Besides, according to the Hergt et al. criterion 23 , the product f·H 0 should be limited for medical reasons.
Let us consider the influence of liquid viscosity on the SAR values. A common place is the statement 6,13,16,18 that the SAR of an assembly in a liquid decreases with an increase in the liquid viscosity. However, as Fig. 1b-d, show, at least for a dilute assembly, neglecting the clustering of the particles, this statement requires a significant correction. The SAR of the assembly does indeed changes inversely with the viscosity, but only in the developed viscous mode, when it reaches a plateau with an increase in the nanoparticle diameter. On the other hand, as Fig. 1b,c show, an increase in the viscosity has only a weak effect on the SAR value at the maximum, in the viscous and intermediate mode, at moderate field amplitudes, H 0 < H k /2. In addition, according to Fig. 1d, in the developed magnetic mode the SAR is practically independent of the viscosity, since in this case the nanoparticle director only weakly oscillates around the ac magnetic field direction.
Effective relaxation time of nanoparticle in liquid. It has recently become popular to apply LRT 24 to interpret the data of the SAR measurements for assemblies distributed in a viscous fluid. In 24 the analytical expression for the SAR of an assembly was formally obtained in the limit of a weak ac magnetic field Here χ 0 = πM 2 s V 3k B T is the initial magnetic susceptibility of an assembly of superparamagnetic nanoparticles 41 , V = πD 3 /6 is the nanoparticle volume, k B is the Boltzmann constant, and ω = 2πf is the circular frequency. Equation (1) contains also the characteristic time of magnetic relaxation of the assembly τ ef . This quantity was taken 24 to be equal to the phenomenological relaxation time suggested by Shliomis 42 (1)  43 . Note, that in the presence of a nonmagnetic layer on the particle surface the particle volume V in this formula should be replaced by the effective volume V η . However, such a replacement does not lead to qualitatively new results; therefore, in this work we assume for simplicity V η = V. Further, τ N is the Néel relaxation time of the average magnetic moment of an assembly of immobile superparamagnetic particles. For this quantity Brown 44 obtained an analytical approximation in the limit of a sufficiently large reduced energy barrier, Here κ is the magnetic damping constant, γ 1 = γ/(1 + κ 2 ), γ is the gyromagnetic ratio.
Obviously, the relaxation of the average magnetic moment of an assembly in a liquid occurs due to the simultaneous action of two processes, namely, due to the motion of the unit magnetization vectors of the particles with respect to the directions of the easy anisotropy axes, and due to the rotation of particles in the liquid as a whole. According to the Shliomis hypothesis 42 , both of these processes occur independently of each other.
Although the assumption of the independence of relaxation processes is not obvious, the Shliomis hypothesis is widely used in the analysis of experimental results 13,16,[27][28][29][30] . In 26 the Shliomis relation was derived by an approximate solution of the Fokker-Planck equation in the limit of a weak external magnetic field. Due to the importance of this hypothesis, in this work the process of relaxation of a dilute assembly of nanoparticles in a liquid is also studied using numerical simulation. www.nature.com/scientificreports/ For this, let us prepare the assembly at time t = 0 in the initial ordered state, when the directors and unit magnetization vectors of all particles are parallel to the z-axis of the Cartesian coordinates, n i = α i = (0,0,1), i = 1, 2, , N p . The time evolution of the assembly at t > 0 can be studied by solving the system of stochastic equations for the motion of vectors α and n given in "Methods" section. The average reduced magnetic moment of the assembly and the average moment of the particle directors during the evolution of the assembly are calculated as follows The last vector characterizes the change in the total spatial orientation of the particles over time. Obviously, as a result of the relaxation of a dilute assembly of superparamagnetic nanoparticles to the equilibrium state, the average values of both vectors should be equal to zero, < m > = < n > = 0.
The calculation of relaxation processes in a liquid was performed for a dilute assembly with the same magnetic parameters, M s = 350 emu/cm 3 , K 1 = 10 5 erg/cm 3 , for particles of various diameters. To obtain statistically reliable results the calculations were carried out for a sufficiently large assemblies of particles, N p = 800-1000. The numerical time step in these calculations was chosen to be 1/30 of the characteristic precession time of unit magnetization vectors T p , since it was shown earlier 45 that such a choice guarantees sufficient accuracy of the numerical calculations.
In Fig. 3a-d black circular dots show the time dependence of the z-component of the average reduced magnetic moment of the assembly, < m z (t) > , for nanoparticles of various diameters, D = 19, 20, 21 and 23 nm, respectively. The red triangular dots show the time evolution of the component < n z (t) > that characterizes the average spatial orientation of nanoparticles of the assembly. For comparison, the blue square dots show the rate As seen from Fig. 3, the curves of orientational relaxation of the assembly in all cases very closely follow the corresponding Debye relaxation curves. Thus, the presence of magnetic degrees of freedom of nanoparticles has practically no effect on the process of orientational relaxation of the assembly in a viscous liquid.
On the other hand, the presence of orientational degrees of freedom of magnetic nanoparticles significantly affects the relaxation process of the average magnetic moment of the assembly. As Fig. 3a-c show, the magnetic relaxation curves are very closely approximated by exponents with the Shliomis relaxation time, Eq. (2). These approximations are shown in Fig. 3a-c with solid lines. Similar results were obtained for dilute assemblies of iron oxide nanoparticles in liquids of different viscosities, η = 0.01-0.2 g/cm/s. Thus, the Shliomis hypothesis is fully justified in the absence of applied magnetic field. It is reasonable to assume that it remains valid in a sufficiently weak static or ac magnetic field.
As Fig. 3 shows, for given magnetic parameters and liquid viscosity the relaxation times τ N and τ B turn out to be nearly equal for particles with diameter D ≈ 21 nm. Since the relaxation time τ N grows exponentially as a function of particle size, the ratio τ N /τ B becomes very large already at D = 23 nm. Then it follows from Eq. (2) that the effective relaxation time τ ef = τ B . Thus, for particles with diameters D ≥ 23 nm, the magnetic relaxation occurs only due to the rotation of particles in the liquid as a whole. Indeed, as Fig. 3d shows, for the case D = 23 nm the relaxation curves < m z (t) > and < n z (t) > practically coincide with each other and with the Debye relaxation curve.
On the other hand, for particles with diameters D ≤ 19 nm, the Néel-Brown relaxation time (3) is much shorter than the Debye relaxation time, therefore τ ef = τ N . Accordingly, in this range of diameters the magnetic relaxation of the assembly occurs due to the Néel-Brown mechanism. The slow orientational relaxation for nanoparticles of such diameters is insignificant. As Fig. 3 shows, the phenomenological Eq. (2) also describes well the time dependence of magnetic relaxation in a narrow range of nanoparticle diameters, 19 nm ≤ D ≤ 23 nm, in which the relaxation regime changes.
Note that under the assumption that the relaxation processes under consideration are independent, a qualitative explanation of the validity of Eq. (2) can be obtained as follows. Let us assume that without taking into account orientational relaxation, the process of magnetic relaxation proceeds according to the law m z (t) = m z (0) exp(− t/τ N ). But due to orientational relaxation, the projection of the average magnetic moment of the assembly on the selected axis decreases according to the equation m z ~ exp(− t/τ B ). As a result of the simultaneous and independent action of both processes, one obtains the time dependence m z (t) ~ exp(− t/τ B )exp(− t/τ N ) ~ exp(− t/τ ef ).

Applicability of linear response theory. As noted in the introduction, at the moment it is not clear what
is the range of validity of the Eqs. (1) and (2). In this section, the dependence of the SAR of the assembly on the amplitude of the ac magnetic field is investigated for magnetic nanoparticles of various diameters. As Fig. 4a,b show, for nanoparticles with diameters D ≤ 19 nm, the linear regime is limited by the value H 0 = 50 Oe. Note that for small particle diameters, when SAR is determined by the Néel mechanism, the relations τ ef ≈ τ N and ωτ ef ≪ 1 are fulfilled. Then one can neglect the factor (ωτ ef ) 2 in the denominator of Eq. (1) to obtain the equation Thus, in this limit the SAR of the assembly does not depend on the liquid viscosity. However, there is an exponential dependence on the reduced energy barrier and a quadratic dependence on the ac field frequency.
According to Fig. 4c,d, with an increase in the particle diameter the linear regime expands up to H 0 = 200 Oe for particles with D = 40 nm. Note that for particles with large diameters, when the condition ωτ ef ≫ 1 is satisfied, Eq. (1) takes the form Therefore, in this area the SAR of the assembly is inversely proportional to the liquid viscosity. In addition, it does not depend on either the particle diameter or the ac field frequency. Numerical data in Fig. 4d confirm the weak dependence of SAR on particle diameters in this limit.
It is interesting to note that according to Fig. 4a-c, with the development of the nonlinear regime the SAR values obtained by numerical simulation for particles of small diameters, D ≤ 20 nm, decrease in comparison with the LRT predictions. The opposite situation is observed in Fig. 4d for large particle diameters, D ≥ 23 nm. Figure 5 shows a comparison of the numerical simulation data with the LRT results depending on the particle diameters at different values of H 0 and η. One can see that numerical data are sufficiently close to LRT only for small particle diameters, D ≤ 16 nm, when the assembly SAR is determined by the Néel mechanism, and in the range of large diameters, D ≥ 24-30 nm, where the viscous friction in a liquid dominates. In the intermediate range of diameters the absorption peak in the LRT is significantly shifted to small particle diameters, so that the difference in SAR values obtained by the two methods is very large, especially for H 0 = 200 Oe. At even higher amplitudes of ac magnetic field the LRT is not applicable.
The difference between the numerical results and LRT shown in Fig. 5 for the SAR of a dilute assembly in a viscous fluid is expected, since LRT is applicable, strictly speaking, under the condition H ≪ H k only.
Indeed, Eq. (1) uses an expression for the initial magnetic susceptibility of the assembly, which is valid only in a weak magnetic field 41 . In addition, Eq. (2), which is essentially used in Eq. (1), characterizes the relaxation

Discussion and conclusions
Based on the calculations performed we come to an interesting conclusion that the effect of viscosity on SAR of a dilute assembly of nanoparticles distributed in a liquid is rather limited. Namely, the SAR of the assembly changes inversely with the liquid viscosity only in the developed viscous mode, at small and moderate amplitudes of ac magnetic field, and for nanoparticles of sufficiently large diameters. At large field amplitudes, comparable to the particle anisotropy field, when the magnetic mode is realized, the SAR of assembly in a liquid significantly exceeds the SAR of the same randomly oriented assembly in a solid matrix. But this is a consequence of the orientation of the particle directors in liquid along the ac magnetic field direction. The viscosity of the liquid affects only the time necessary for establishment of the stationary magnetic mode. The latter can be as long as several tens or even hundreds of ac magnetic field periods for particles of large diameters. Obviously, the SAR of an assembly in a liquid should increase as a result of partial orientation of the particle directors due to the application of external static magnetic field H dc parallel to the ac field direction. Experimental confirmation of this fact was obtained in 46,47 . It is also known that the SAR of an assembly increases significantly upon the formation of one-dimensional dipole-bound chains of nanoparticles in a viscous liquid [48][49][50][51] . In this regard, it should be noted that the spontaneous orientation of the particle directors in a viscous fluid, predicted in 31 and confirmed experimentally 38,39 , occurs only under the action of ac magnetic field of sufficiently large amplitude, in the absence of external dc magnetic field, and in a sufficiently dilute assembly, neglecting the influence of magnetic dipole interactions between nanoparticles.
It is also worth noting that the calculated SAR values shown in Figs. 1, 4 and 5 were obtained in this work by numerically integrating low frequency hysteresis loops of assemblies of superparamagnetic nanoparticles according to the well-known 40,52,53 thermodynamic formula, SAR = M s f � αd � H ρ , without invoking any additional hypotheses about the dynamic magnetic susceptibility of the assembly 27 , or the relaxation times of nanoparticles 30,54 . Low-frequency hysteresis loops, examples of which are shown in Fig. 2d, were calculated according to the approach described in "Methods" section. As shown earlier 45 , this calculation procedure www.nature.com/scientificreports/ reproduces quite accurately the known analytical approximations for the relaxation times τ N of nanoparticles in a solid matrix. It also accurately reproduces the Debye time τ B of orientational relaxation of an assembly of nonmagnetic nanoparticles in a viscous liquid. It is worth noting that the low frequency hysteresis loops of the assembly allow direct experimental measurement 3,55,56 . All these facts constitute a solid basis for the approach outlined in this paper. At the same time, the recent attempts 30,54 to introduce the relaxation times, which depend on the ac magnetic field amplitude and frequency seem to us physically unreasonable at sufficiently large amplitudes and frequencies, which are essential for magnetic hyperthermia. Indeed, it is impossible to indicate a direct physical experiment in which such relaxation times could be directly measured. At the same time, the effective Schliomis relaxation time, Eq. (2), has a clear physical meaning, since this quantity, at least in principle, allows direct physical measurement. The Shliomis's hypothesis is not trivial. First, it refers only to the process of relaxation of the average magnetic moment of the assembly, and says nothing about the simultaneous orientational relaxation of the nanoparticles of the assembly. In this work it is shown that orientational relaxation in a liquid of finite viscosity occurs with the Debye relaxation time τ B . Besides, it was shown 57 that in a system with very low viscosity, close to zero, the strong interaction of the rotational and magnetic degrees of freedom of a nanoparticle can lead to new physical effects.

Methods: set of stochastic Landau-Lifshitz and director's equations
In Ref. 31 a system of stochastic equations was proposed to study the coupled dynamics of the unit magnetization vector of the particle α and the particle director n in a viscous liquid at a finite temperature. The dynamics of the unit magnetization vector is described by the stochastic Landau-Lifshitz equation 44,[58][59][60] Here the effective magnetic field is calculated as the derivative of the total particle energy where H is the vector of a static or variable external magnetic field, H(t) = H 0 cos(ωt). In turn, the dynamics of the nanoparticle director is determined by the stochastic equation Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.