arxiv:physics/ v1 15 Nov 2006

Transcription

1 Magnetophoresis of nonmagnetic particles in ferrofluids Y. Gao, Y. C. Jian, L. F. Zhang and J. P. Huang Surface Physics Laboratory (National Key Laboratory) and Department of Physics, Fudan University, Shanghai , China Taking into account the structural transition and long-range interaction (lattice effect), we resort to the Ewald- Kornfeld formulation and the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions to calculate the effective permeability of the inverse ferrofluid. And we also consider the effects of the volume fraction and geometric shape of the nonmagnetic particles to the magnetophoretic force on the nonmagnetic particles. It is found that the absolute value of the magnetophoretic force decreases as the volume fraction increases. For nonspherical case, the absolute value of the magnetophoretic force decreases as the aspect ratio (k) increases. arxiv:physics/ v1 15 Nov 2006 I. INTRODUCTION Magnetizable particles experience a force in a nonuniform magnetic field, which comes from the interaction between the induced dipole moment inside the particles and the external nonuniform magnetic field. This phenomenon, sometimes called magnetophoresis (M AP), which has been exploited in a variety of industrial and commercial process for separation and beneficiation of solids suspended in liquids [1]. The nonmagnetic particle suspended in the ferrofluid possesses magnetic moment, so it also experiences a force in a nonuniform magnetic field. Magnetophoresis may open up an approach to characterize and separate cells for biochemical analysis based on intrinsic and extrinsic magnetic properties of biological macromolecules [2]. In recent years, inverse ferrofluids with nonmagnetic colloidal particles (or magnetic holes) suspended in a ferrofluid have drawn much attention for their potential industrial and biomedical applications [3, 4, 5, 6]. The size of the nonmagnetic particle (such as polystyrene particle) is about 1 100µm. Ferrofluid, also known as magnetic fluid, is a colloidal suspension of single-domain magnetic particles, with typical dimension around 10nm, dispersed in a liquid carrier [7, 8]. Since the nonmagnetic particles are much larger than the magnetic nanoparticles, the magnetic fluid can be treated as a one-component continuum with respect to the larger nonmagnetic particles [9, 10]. Recently, Jian, Huang and Tao first theoretically suggested that the ground state of inverse ferrofluid is body-centered tetragonal (bct) lattice [11], just as the ground states of electrorheological fluids and magnetorheological fluids [12, 13]. Here, we propose an alternative structural transition from the bct to the face-centered cubic (f cc), through the application of external magnetic field. Then, we apply a nonuniform magnetic field to investigate the magnetophoresis of the nonmagnetic particles submerged in the ferrofluid. In fact, one has investigated the force on the microparticles in electrorheological solids under the application of a nonuniform ac electric field by using a self-consistent method and the spectral representation approach [14]. In our present work, we use the Ewald-Kornfeld formulation and the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions to calculate the effective permeability of the inverse ferrofluid, considering the lattice effect, the volume fraction and shape of the nonmagnetic particle. Based on this, the magnetophoretic force on the nonmagnetic particle can be calculated. This paper is organized as follows. In Sec. II, we introduce the theory of magnetophoresis. By using the Ewald-Kornfeld formulation and the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions, we calculate the effective permeability of the inverse ferrofluids. And we investigate the magnetophoretic force on the nonmagnetic particles in ferrofluids. This paper ends with a conclusion in Sec. III. II. THEORY AND FORMALISM A. Theory of magnetophoresis Let us first consider the magnetically linear particles. For a magnetically linear particle, the effective magnetic dipole moment vector m eff takes a form very similar to that for the effective dipole moment of dielectric sphere, m eff = 4πR 3 K(µ 1, µ 2 ) H 0, (1) K(µ 1, µ 2 ) = µ 1 µ 2 µ 1 + 2µ 2, (2) where R is the radius of the magnetically linear particle, µ 1 is the permeability of the magnetically linear particle, and µ 2 is the permeability of the host medium. The force exerted by a nonuniform magnetic field on a magnetic dipole, F dipole = µ 2 m eff H 0, (3) Combining this expression with that for the effective magnetic dipole moment, the magnetophoretic force exerted on a magnetizable spherical particle in a nonuniform magnetic field can be written as [1], Electronic address: jphuang@fudan.edu.cn F MAP = 2πµ 2 R 3 K(µ 1, µ 2 ) H 2 0, (4)

2 2 Eq. (4) is similar to the expression of the dielectrophoretic force on particles. The magnetophoretic force on particles are proportional to the particle volume and the permeability of the suspension medium. And the force is directed along the gradient of the magnetic field intensity. There are two kinds of magnetophoresis depending upon the relative magnitudes of µ 1 and µ 2. If K > 0 (µ 1 > µ 2 ), particles are attracted to magnetic field intensity maxima and repelled from minima. This phenomenon is called positive magnetophoresis. On the contrary, if K < 0 (µ 1 < µ 2 ), particles are attracted to magnetic field intensity minima and repelled from maxima. This case is called negative magnetophoresis. Now, we consider the nonlinear magnetophoresis. There are two kinds of nonlinear magnetophoresis. The first case is ferromagnetic particles in a linear magnetizable medium. The other one is nonmagnetic particles in a magnetizable liquid, such as a ferrofluid. In both cases, The magnetic field is assumed to be dc and hysteresis is ignored. In this work, the second case will be considered. For the reason that the magnetic nonlinearity of ferrofluid causes difficulties in calculating the effective magnetic moment for an immersed nonmagnetic particle. Therefore, we use a Maxwell stress tensor formulation to predict the magnetic force F ff on nonmagnetic immersed particles, F ff = Σ { u 0 M 2 n 2 + µ 0 H0 0 M(H )dh } ˆnds. (5) Where Σ is the surface that just encloses the object, H 0 is the magnitude of the nonuniform magnetic field, ˆn is the outwarddirected unit vector normal the surface, and M n = ˆn M( H) is the normal component of the magnetization in the ferrofluid adjacent to the surface Σ enclosing the particle. In a weak magnetic field, we can use Eq. (4) by setting µ 1 = µ 0 to calculate the magnetophorestic force on the nonmagnetic particles submerged in a ferrofluid. B. Effective permeability The effective permeability can be obtained by performing an Ewald-Kornfeld formulation so that the structural transition and long-range interaction can be taken into account explicitly. The ground state of inverse ferrofluid is a bodycentered-tetragonal (bct) lattice, which can be regarded as a tetragonal lattice, plus a basis of two particles each of which is fixed with an induced point magnetic dipole at its center. One of the two particles is located at a corner and the other one at the body center of the tetragonal unit cell. The tetragonal lattice has a lattice constant c 1 = ξq along the z axis and lattice constants a 1 = b 1 = ξq 1 2 along x and y axes respectively. In this paper we identified a different structure transformation from the bct to fcc by changing the uniaxial lattice constant c 1 under hard-sphere constraint. The volume of the tetragonal unit cell remains certain value (ξ 3 ) as q varies. In this way, the degree of anisotropy of tetragonal lattice is measured by how q is deviated from unity and the uniaxial anisotropic axis is along the z axis. In particular, q = , 1 and represent the bct, bcc and fcc lattice, respectively. When an external magnetic field H 0 is along the x axis, the induced dipole moment m are perpendicular to the uniaxial anisotropic axis. At the lattice point R=0, the local field H L can be determined by using the Ewald-Kornfeld formulation [15, 16, 17], H L = m 2 [ γ 1 (R j ) + x 2 jq 2 γ 2 (R j )] j=1 R 0 4π m v c G 0 Π( G) G2 x G 2 exp ( G 2 4η 2 where γ 1 and γ 2 are two coefficients, given by, ) + 4 mη3 3 π.(6) γ 1 (r) = erfc(ηr) r 3 + 2η πr 2 exp( η2 r 2 ), (7) γ 2 (r) = 3erfc(ηr) r 5 + ( 4η3 + 6η )exp πr 2 πr 4 (η2 r 2 ), (8) where erfc(ηr) is the complementary error function, and η is an adjustable parameter making the summations in real and reciprocal lattices converge rapidly. In Eq. (6), R and G denote the lattice vector and reciprocal lattice vector, respectively, R = ξ(q 1 2 lˆx + q 1 2 mŷ + qnẑ), (9) G = 2π ξ (q 1 2 uˆx + q 1 2 vŷ + q 1 wẑ), (10) where l, m, n, u, v, w are integers. In addition, z j and R j in Eq. (6) are given by, z j = l j 1 2, R j = R j 1 2 (a 1ˆx + a 1 ŷ + c 1 ẑ), (11) and Π( G) = 1 + exp[i(u + v + w)/π] is the structure factor. The local field can be obtained by summing over all integer indices in real lattice and the reciprocal lattice. Now, let us define a local field factor α, α = 3v c 4π H L m. (12) It is worth remaking that α is a function of single variable q. Also there is a sum rule for α T and α L, 2α T + α L = 3 [18, 19], where α L denotes the local field factor in longitudinal field case, and α T denotes the local field in transverse field case. Here the longitudinal (or transverse) field case corresponds to the fact the H 0 field is parallel (or perpendicular) to the uniaxial anisotropic z axis. We can use the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions [20, 21],

3 3 µ e µ 2 αµ e + (3 α)µ 2 = f µ 1 µ 2 µ 1 + 2µ 2, (13) where µ e stands for the effective permeability of the inverse ferrofluid, µ 1 stands for the permeability of the nonmagnetic particles, µ 2 stands for the permeability of ferrofluid, and f is the volume fraction of the nonmagnetic particles submerged in a ferrofluid. It should be noted that Eq. (13) is fit for the case of the spherical nonmagnetic particle. For the cases of nonsphere (such as prolate and oblate spheroids), we resort to the demagnetizing factor g to describe the shape of spheroids. For a prolate spheroid, there is a = b < c, where a, b, c are three principal axes of the spheriods, and the longitudinal demagnetizing factor g L along the major axis c is, g L = 1 k 2 1 [ k k2 1 ln(k + ] k 2 1) 1, (14) where k = c/a, and g L < 1/3. For an oblate spheroid with a < b = c, the longitudinal demagnetizing factor g L along the major axis a is [ ] k 2 k2 1 g L = 1 arcsin 1 (k 2 1) 3 2 k k 2. (15) 1 where g L > 1/3. And there is a sum rule for g L and g T, 2g T + g L = 1, where g T denotes the transverse demagnetizing factor along the minor axis. It should be noted that the nonmagnetic particles are spherical ones at g L = g T = 1/3. The effective permeability of the inverse ferrofluid can be determined by the well-known developed Maxwell-Garnett equation [22], α L3 (µ e µ 2 ) µ 2 + αl 3 (µ e µ 2 ) = f g L (µ 1 µ 2 ) µ 2 + g L (µ 1 µ 2 ), (16) where α L is local field factor in longitudinal field case, and g L is the longitudinal demagnetizing factor. C. Results and Discussion F = F MAP /2πR 3 H 2 0 F 1 = F MAP /(2πa 2 c H 2 0 ), If the shape of the nonmagnetic particle is spherical, the magnetophoretic force on the nonmagnetic particle can be calculated by the following equations: F MAP = 2πµ e R 3 K(µ 1, µ e ) H 2 0, (17) K(µ 1, µ e ) = µ 1 µ e µ 1 + 2µ e. (18) Considering the lattice effect, µ 2 is substituted by µ e. During the numerical calculation, for convenience, F stands for F MAP /2πR 3 H0 2. It is worth mentioning the negative sign. In our calculation, we choose the parameters µ 1 = 1.0, µ 2 = 2.9 (The permeability of the nonmagnetic particle is µ 0, and the permeability of the typical ferrofluid is about 2.9 µ 0 ). In this case, the parameter K < 0 is corresponding to the negative magnetophoresis. As a matter of fact, we want to obtain the absolute value of the magnetophoretic force on the nonmagnetic particle. Figure 1 displays the dependence of F on q when the external magnetic field is parallel to the uniaxial anisotropic axis for three different volume fractions f=0.15, f=0.25 and f=0.35. While Figure 2 shows the dependence of F on q when the external magnetic field is perpendicular to the uniaxial anisotropic axis for three different volume fractions f=0.15, f=0.25 and f=0.35. The two figures show that the magnetophoretic force on the nonmagnetic particles are different according to different lattices. Because of the anisotropy, Fig.1 displays that F decreases as q increases. Inverse behavior appears in Fig2. Without the effect of the structural transition and the long-range interaction (lattice effect), F equals to With the effect, take the bcc lattice for example, F equals to when the volume fraction equals to The difference between the above two cases is large, which means that the lattice effect plays an important role in Magnetophoresis of nonmagnetic particles in ferrofluids. Now, let us focus our interest on the effect of the volume fraction to the magnetophoretic force on the nonmagnetic particles. Fig.?? represents the relation between F and q for different volume fractions of the nonmagnetic particles. It shows that the absolute value of the magnetophoretic force on the nonmagnetic particles decreases as the volume fraction increases. The relation between F and f for the external magnetic field parallel to the uniaxial anisotropic axis is similar to that of the external magnetic field perpendicular to the uniaxial anisotropic axis, (see Fig. 2). As the volume fraction of the nonmagnetic particles increases, the lattice effect becomes stronger. By using the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions, a smaller effective permeability of the inverse ferrofluid can be obtained. Consequently, the absolute value of the magnetophoretic force on the nonmagnetic particles becomes smaller. Shape-induced anisotropy makes the magnetophorestic force expression more complicated for an ellipsoidal particle than it is for a sphere. For the special case of a prolate particle (a = b < c) aligned with its long axis parallel to and positioned on the axis of an axially symmetric field, the magnetophorestic force expression resembles that for a sphere, [ ] F MAP = 2πa 2 1 (µ 1 µ 2 ) c H ( µ1 µ2 0 2 µ 2 )g, (19) L where µ 1 is the permeability of the particle, µ 2 is the permeability of the host medium, and g L is the longitudinal demagnetizing factor. From the above equation, we can find that the force remains proportional to the scalar product of particle volume and excess permeability. Ellipsoidal particles exhibit

4 4 the same positive and negative magnetophoresis phenomenology, depending on whether µ 1 is greater than or less than µ 2. This expression is also valid for oblate spheroids (a < b = c), though such a particle with its short axis parallel to the magnetic field will be in unstable alignment. For a spherical particle, k = 1 and g L = 1/3, in which case Eq.(19) reduces to Eq. (4). In order to calculate the magnetophoretic force on the nonmagnetic particles, the lattice effect must be considered. So µ 2 should be substituted by µ e. In this paper, for convenience, F 1 stands for F MAP /(2πa 2 c H0), 2 which means that we only calculate the coefficient of the magnetophoretic force. Here, we assume the volume of the nonmagnetic particle is invariable, and H0 2 is a constant across the nonmagnetic particles. For certain lattice, such as bct, bcc and fcc, the relation between F 1 and the aspect ratio k are shown in Fig.?? and Fig.??. The absolute value of the magnetophoretic force on the nonmagnetic particle decreases as k increases for both prolate and oblate cases. In addition, the relation between F 1 and the volume fraction of the nonmagnetic particles is also investigated. The results are shown in Fig.?? and Fig.??. The absolute value of the magnetophoretic force on the nonmagnetic particle decreases as the volume fraction increases, which are consistent with the result of the spherical case. We can calculate the transverse case by using the sum rules 2α T + α L = 3 and 2g T + g L = 1. There are no figures shown here for such case. force on the nonmagnetic particle in the ferrofluid due to the nonuniform magnetic field can be affected significantly by the lattice effect, the volume fraction and shape of the nonmagnetic particle. Thus, it is possible to detect the structure of the inverse ferrofluid, volume fraction and shape of the nonmagnetic particle by measuring such magnetophoretic force. Acknowledgements The author Y. Gao greatly appreciated Prof. X. Sun for his generously help in his study. The authors acknowledge the financial support by the Shanghai Education Committee and the Shanghai Education Development Foundation (Shu Guang project) under Grant No. KBH , by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry, China, by the National Natural Science Foundation of China under Grant No III. CONCLUSIONS In fact, Magnetophoresis is a complicated phenomenon. Such as magnetophoresis with eddy current induction, which is corresponding to the electrically conductive magnetizable particle. In our present work, we only study the magnetophoretic force on the nonmagnetic particles (or magnetic holes) suspended in a host ferrofluid, by taking into account the effect of structural transition and long-range interaction (lattice effect). So we resort to the Ewald-Kornfeld formulation and the well-known developed Maxwell-Garnett theory to calculate the effective permeability of the inverse ferrofluid. And we also consider the effects of the volume fraction and geometric shape of the nonmagnetic particles to the magnetophoretic force on the nonmagnetic particles. It is found that the absolute value of the magnetophoretic force on the nonmagnetic particle decreases as the volume fraction increases. For nonspherical case, the absolute value of the magnetophoretic force on the nonmagnetic particle decreases as the aspect ratio(k) increases. As the nonmagnetic particles are located closely, the multipolar interaction between them is expected to play a role [23, 24]. In this regard, it is of value to extend the present work to investigate the effect of the multipolar interaction. In addition, it is also interesting to see what happens if the external field is a nonuniform ac magnetic field and the nonmagnetic particle is inhomogeneous. To sum up, our results have shown that the magnetophoretic

5 5 [1] T. B. Jones, Electromechanics of Particles (Cambridge University Press, Cambridge, 1995), Chap.III. [2] Maciej Zborowski, Graciela R. Ostera, Lee R. Moore, Sarah Milliron, Jeffrey J. Chalmers and Alan N. Schechter, Biophysical Journal 84, 2638 (2003) [3] A. T. Skjeltorp, Phys. Rew. Lett, 51, 2306 (1983). [4] J. Ugelstad et al, Blood Purif. 11, 349 (1993). [5] J. B. Hayter, R. Pynn, S. Charles, A. T. Skjeltorp, J. Trewhella, G. Stubbs, and P. Timmins, Phys. Rev. Lett. 62, 1667 (1989). [6] A. Koenig, P. Hebraud, C. Gosse, R. Dreyfus, J. Baudry, E. Bertrand, and J. Bibette, Phys. Rev. Lett. 95, (2005). [7] C.Scherer and A. M. Figueiredo Neto, Brazilian Journal of Physics, 35 3A, (2005). [8] Rosenweig, R. E. Ferrohydrodynamics (Cambridge University Press, Cambridge,1985). [9] Efraim Feinstein and Mara Prentiss, J. Appl. Phys. 99, (2006). [10] B. J. de Gans, C. Blom, A. P. Philipse, and J. Mellema, Phys. Rev. E 60, 4518 (1999). [11] Y. C. Jian, J. P. Huang and R. Tao, Structure of polydisperse inverese ferrofluids: Theory and computer simulation, (to be published). [12] R. Tao and J. M. Sun, Phys. Rev. Lett. 67, 398 (1991). [13] Lei Zhou, Weijia Wen, and Ping Sheng, Phys. Rev. Lett. 81, 1509 (1998). [14] J. P. Huang, Chem. Phys. Lett. 390, 380 (2004). [15] J. P. Huang, Phys. Rev. E 70, (2004). [16] P. P. Ewald, Ann. Phys. (Leipzig) 64, 253 (1921); H. Kornfeld, Z. Phys. 22, 27 (1924). [17] J. P. Huang and K. W. Yu, Appl. Phys. Lett. 87, (2005). [18] C.K. Lo and K.W. Yu, Phys. Rev. E 64, (2001). [19] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electrodynamics of Continuous Media, 2nd ed. (Pergamon, New York, 1984), Chap. II. [20] J. C. M. Garnett, Philos. Trans. R. Soc. London, Ser. A 203,385 (1904); 205, 237 (1906). [21] J. P. Huang, J. T. K. Wan, C. K. Lo, and K. W. Yu, Phys. Rev. E 64, (2001). [22] G. Wang, W. J. Tian and J. P. Huang, J. Phys. Chem. B 110, (2006) [23] K.W. Yu, J.T.K. Wan, Comput. Phys. Commun. 129, 177 (2000). [24] J.P. Huang, M. Karttunen, K.W. Yu, L. Dong, Phys. Rev. E 67, (2003). Figure Captions Fig.1. (Color online) Cases of spherical particles. F as a function of q for different volume fractions. The external nonuniform magnetic field is parallel to the uniaxial anisotropic axis. Parameters: µ 1 = 1.0, µ 2 = 2.9. Fig.2. (Color online) Same as Fig. 1, but for transverse case. Fig.3. (Color online) Cases of prolate spheroidal particles. F 1 as a function of k for different latices. The external nonuniform magnetic field is parallel to the uniaxial anisotropic axis. Parameters: µ 1 = 1.0, µ 2 = 2.9, f = Fig.4. (Color online) Same as Fig. 3, but for oblate spheroidal case. Fig.5. (Color online) Cases of prolate spheroidal particles. F 1 as a function of f for different latices. The external nonuniform magnetic field is parallel to the uniaxial anisotropic axis. Parameters: µ 1 = 1.0, µ 2 = 2.9, k = 3.0. Fig.6. (Color online) Same as Fig. 5, but for oblate spheroidal case.

6 FIG. 1: Gao, Jian, Zhang and Huang 6

7 FIG. 2: Gao, Jian, Zhang and Huang 7

8 FIG. 3: Gao, Jian, Zhang and Huang 8

9 Magnetophoresis of nonmagnetic particles in ferrofluids Y. Gao, Y. C. Jian, L. F. Zhang and J. P. Huang Surface Physics Laboratory (National Key Laboratory) and Department of Physics, Fudan University, Shanghai , China Abstract arxiv:physics/ v1 15 Nov 2006 Taking into account the structural transition and long-range interaction (lattice effect), we resort to the Ewald-Kornfeld formulation and the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions to calculate the effective permeability of the inverse ferrofluid. And we also consider the effects of the volume fraction and geometric shape of the nonmagnetic particles to the magnetophoretic force on the nonmagnetic particles. It is found that the absolute value of the magnetophoretic force decreases as the volume fraction increases. For nonspherical case, the absolute value of the magnetophoretic force decreases as the aspect ratio (k) increases. PACS numbers: Electronic address: Typeset by REVTEX 1

10 I. INTRODUCTION Magnetizable particles experience a force in a nonuniform magnetic field, which comes from the interaction between the induced dipole moment inside the particles and the external nonuniform magnetic field. This phenomenon, sometimes called magnetophoresis (MAP), which has been exploited in a variety of industrial and commercial process for separation and beneficiation of solids suspended in liquids [1]. The nonmagnetic particle suspended in the ferrofluid possesses magnetic moment, so it also experiences a force in a nonuniform magnetic field. Magnetophoresis may open up an approach to characterize and separate cells for biochemical analysis based on intrinsic and extrinsic magnetic properties of biological macromolecules [2]. In recent years, inverse ferrofluids with nonmagnetic colloidal particles (or magnetic holes) suspended in a ferrofluid have drawn much attention for their potential industrial and biomedical applications [3, 4, 5, 6]. The size of the nonmagnetic particle (such as polystyrene particle) is about 1 100µm. Ferrofluid, also known as magnetic fluid, is a colloidal suspension of singledomain magnetic particles, with typical dimension around 10nm, dispersed in a liquid carrier [7, 8]. Since the nonmagnetic particles are much larger than the magnetic nanoparticles, the magnetic fluid can be treated as a one-component continuum with respect to the larger nonmagnetic particles [9, 10]. Recently, Jian, Huang and Tao first theoretically suggested that the ground state of inverse ferrofluid is body-centered tetragonal (bct) lattice [11], just as the ground states of electrorheological fluids and magnetorheological fluids [12, 13]. Here, we propose an alternative structural transition from the bct to the face-centered cubic (f cc), through the application of external magnetic field. Then, we apply a nonuniform magnetic field to investigate the magnetophoresis of the nonmagnetic particles submerged in the ferrofluid. In fact, one has investigated the force on the microparticles in electrorheological solids under the application of a nonuniform ac electric field by using a self-consistent method and the spectral representation approach [14]. In our present work, we use the Ewald-Kornfeld formulation and the well-known developed Maxwell- Garnett theory for uniaxially anisotropic suspensions to calculate the effective permeability of the inverse ferrofluid, considering the lattice effect, the volume fraction and shape of the nonmagnetic particle. Based on this, the magnetophoretic force on the nonmagnetic particle can be calculated. This paper is organized as follows. In Sec. II, we introduce the theory of magnetophoresis. By using the Ewald-Kornfeld formulation and the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions, we calculate the effective permeability of the inverse ferroflu- 2

11 ids. And we investigate the magnetophoretic force on the nonmagnetic particles in ferrofluids. This paper ends with a conclusion in Sec. III. II. THEORY AND FORMALISM A. Theory of magnetophoresis Let us first consider the magnetically linear particles. For a magnetically linear particle, the effective magnetic dipole moment vector m eff takes a form very similar to that for the effective dipole moment of dielectric sphere, m eff = 4πR 3 K(µ 1, µ 2 ) H 0, (1) K(µ 1, µ 2 ) = µ 1 µ 2 µ 1 + 2µ 2, (2) where R is the radius of the magnetically linear particle, µ 1 is the permeability of the magnetically linear particle, and µ 2 is the permeability of the host medium. The force exerted by a nonuniform magnetic field on a magnetic dipole, F dipole = µ 2 m eff H 0, (3) Combining this expression with that for the effective magnetic dipole moment, the magnetophoretic force exerted on a magnetizable spherical particle in a nonuniform magnetic field can be written as [1], F MAP = 2πµ 2 R 3 K(µ 1, µ 2 ) H0 2, (4) Eq. (4) is similar to the expression of the dielectrophoretic force on particles. The magnetophoretic force on particles are proportional to the particle volume and the permeability of the suspension medium. And the force is directed along the gradient of the magnetic field intensity. There are two kinds of magnetophoresis depending upon the relative magnitudes of µ 1 and µ 2. If K > 0 (µ 1 > µ 2 ), particles are attracted to magnetic field intensity maxima and repelled from minima. This phenomenon is called positive magnetophoresis. On the contrary, if K < 0 (µ 1 < µ 2 ), particles are attracted to magnetic field intensity minima and repelled from maxima. This case is called negative magnetophoresis. Now, we consider the nonlinear magnetophoresis. There are 3

12 two kinds of nonlinear magnetophoresis. The first case is ferromagnetic particles in a linear magnetizable medium. The other one is nonmagnetic particles in a magnetizable liquid, such as a ferrofluid. In both cases, The magnetic field is assumed to be dc and hysteresis is ignored. In this work, the second case will be considered. For the reason that the magnetic nonlinearity of ferrofluid causes difficulties in calculating the effective magnetic moment for an immersed nonmagnetic particle. Therefore, we use a Maxwell stress tensor formulation to predict the magnetic force F ff on nonmagnetic immersed particles, { M F 2 H0 } n ff = u 0 Σ 2 + µ 0 M(H )dh ˆnds. (5) 0 Where Σ is the surface that just encloses the object, H 0 is the magnitude of the nonuniform magnetic field, ˆn is the outward-directed unit vector normal the surface, and M n = ˆn M( H) is the normal component of the magnetization in the ferrofluid adjacent to the surface Σ enclosing the particle. In a weak magnetic field, we can use Eq. (4) by setting µ 1 = µ 0 to calculate the magnetophorestic force on the nonmagnetic particles submerged in a ferrofluid. B. Effective permeability The effective permeability can be obtained by performing an Ewald-Kornfeld formulation so that the structural transition and long-range interaction can be taken into account explicitly. The ground state of inverse ferrofluid is a body-centered-tetragonal (bct) lattice, which can be regarded as a tetragonal lattice, plus a basis of two particles each of which is fixed with an induced point magnetic dipole at its center. One of the two particles is located at a corner and the other one at the body center of the tetragonal unit cell. The tetragonal lattice has a lattice constant c 1 = ξq along the z axis and lattice constants a 1 = b 1 = ξq 1 2 along x and y axes respectively. In this paper we identified a different structure transformation from the bct to f cc by changing the uniaxial lattice constant c 1 under hard-sphere constraint. The volume of the tetragonal unit cell remains certain value (ξ 3 ) as q varies. In this way, the degree of anisotropy of tetragonal lattice is measured by how q is deviated from unity and the uniaxial anisotropic axis is along the z axis. In particular, q = , 1 and represent the bct, bcc and fcc lattice, respectively. When an external magnetic field H 0 is along the x axis, the induced dipole moment m are perpendicular to the uniaxial anisotropic axis. At the lattice point R=0, the local field H L can be determined by using the Ewald-Kornfeld formulation [15, 16, 17], 4

13 H L = m 2 [ γ 1 (R j ) + x 2 j q2 γ 2 (R j )] j=1 R 0 4π m v c G 0 where γ 1 and γ 2 are two coefficients, given by, γ 1 (r) = erfc(ηr) r 3 Π( G) G2 x G 2 exp ( G 2 4η 2 ) + 4 mη3 3 π. (6) + 2η πr 2 exp( η2 r 2 ), (7) γ 2 (r) = 3erfc(ηr) r 5 + ( 4η3 πr 2 + 6η πr 4 ) exp (η2 r 2 ), (8) where erfc(ηr) is the complementary error function, and η is an adjustable parameter making the summations in real and reciprocal lattices converge rapidly. In Eq. (6), R and G denote the lattice vector and reciprocal lattice vector, respectively, R = ξ(q 1 2 lˆx + q 1 2 mŷ + qnẑ), (9) G = 2π ξ (q 1 2 uˆx + q 1 2 vŷ + q 1 wẑ), (10) where l, m, n, u, v, w are integers. In addition, z j and R j in Eq. (6) are given by, z j = l j 1 2, R j = R j 1 2 (a 1ˆx + a 1 ŷ + c 1 ẑ), (11) and Π( G) = 1 + exp[i(u + v + w)/π] is the structure factor. The local field can be obtained by summing over all integer indices in real lattice and the reciprocal lattice. Now, let us define a local field factor α, α = 3v c 4π H L m. (12) It is worth remaking that α is a function of single variable q. Also there is a sum rule for α T and α L, 2α T + α L = 3 [18, 19], where α L denotes the local field factor in longitudinal field case, and α T denotes the local field in transverse field case. Here the longitudinal (or transverse) field case corresponds to the fact the H 0 field is parallel (or perpendicular) to the uniaxial anisotropic 5

14 z axis. We can use the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions [20, 21], µ e µ 2 αµ e + (3 α)µ 2 = f µ 1 µ 2 µ 1 + 2µ 2, (13) where µ e stands for the effective permeability of the inverse ferrofluid, µ 1 stands for the permeability of the nonmagnetic particles, µ 2 stands for the permeability of ferrofluid, and f is the volume fraction of the nonmagnetic particles submerged in a ferrofluid. It should be noted that Eq. (13) is fit for the case of the spherical nonmagnetic particle. For the cases of nonsphere (such as prolate and oblate spheroids), we resort to the demagnetizing factor g to describe the shape of spheroids. For a prolate spheroid, there is a = b < c, where a, b, c are three principal axes of the spheriods, and the longitudinal demagnetizing factor g L along the major axis c is, g L = 1 [ k k 2 1 k2 1 ln(k + ] k 2 1) 1, (14) where k = c/a, and g L < 1/3. For an oblate spheroid with a < b = c, the longitudinal demagnetizing factor g L along the major axis a is [ g L = 1 k 2 (k 2 1) 3 2 arcsin k2 1 k ] 1 k 2 1. (15) where g L > 1/3. And there is a sum rule for g L and g T, 2g T + g L = 1, where g T denotes the transverse demagnetizing factor along the minor axis. It should be noted that the nonmagnetic particles are spherical ones at g L = g T = 1/3. The effective permeability of the inverse ferrofluid can be determined by the well-known developed Maxwell-Garnett equation [22], α L3 (µ e µ 2 ) µ 2 + α L 3 (µ e µ 2 ) = f g L (µ 1 µ 2 ) µ 2 + g L (µ 1 µ 2 ), (16) where α L is local field factor in longitudinal field case, and g L is the longitudinal demagnetizing factor. C. Results and Discussion F = F MAP /2πR 3 H 2 0 F 1 = F MAP /(2πa 2 c H 2 0 ), 6

15 If the shape of the nonmagnetic particle is spherical, the magnetophoretic force on the nonmagnetic particle can be calculated by the following equations: F MAP = 2πµ e R 3 K(µ 1, µ e ) H 2 0, (17) K(µ 1, µ e ) = µ 1 µ e µ 1 + 2µ e. (18) Considering the lattice effect, µ 2 is substituted by µ e. During the numerical calculation, for convenience, F stands for F MAP /2πR 3 H0 2. It is worth mentioning the negative sign. In our calculation, we choose the parameters µ 1 = 1.0, µ 2 = 2.9 (The permeability of the nonmagnetic particle is µ 0, and the permeability of the typical ferrofluid is about 2.9 µ 0 ). In this case, the parameter K < 0 is corresponding to the negative magnetophoresis. As a matter of fact, we want to obtain the absolute value of the magnetophoretic force on the nonmagnetic particle. Figure 1 displays the dependence of F on q when the external magnetic field is parallel to the uniaxial anisotropic axis for three different volume fractions f=0.15, f=0.25 and f=0.35. While Figure 2 shows the dependence of F on q when the external magnetic field is perpendicular to the uniaxial anisotropic axis for three different volume fractions f=0.15, f=0.25 and f=0.35. The two figures show that the magnetophoretic force on the nonmagnetic particles are different according to different lattices. Because of the anisotropy, Fig.1 displays that F decreases as q increases. Inverse behavior appears in Fig2. Without the effect of the structural transition and the long-range interaction (lattice effect), F equals to With the effect, take the bcc lattice for example, F equals to when the volume fraction equals to The difference between the above two cases is large, which means that the lattice effect plays an important role in Magnetophoresis of nonmagnetic particles in ferrofluids. Now, let us focus our interest on the effect of the volume fraction to the magnetophoretic force on the nonmagnetic particles. Fig.?? represents the relation between F and q for different volume fractions of the nonmagnetic particles. It shows that the absolute value of the magnetophoretic force on the nonmagnetic particles decreases as the volume fraction increases. The relation between F and f for the external magnetic field parallel to the uniaxial anisotropic axis is similar to that of the external magnetic field perpendicular to the uniaxial anisotropic axis, (see Fig. 2). As the volume fraction of the nonmagnetic particles increases, the lattice effect becomes stronger. By using the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions, a smaller effective permeability of the inverse ferrofluid can be obtained. Consequently, the absolute value of the magnetophoretic force on the nonmagnetic particles becomes smaller. 7

16 Shape-induced anisotropy makes the magnetophorestic force expression more complicated for an ellipsoidal particle than it is for a sphere. For the special case of a prolate particle (a = b < c) aligned with its long axis parallel to and positioned on the axis of an axially symmetric field, the magnetophorestic force expression resembles that for a sphere, [ ] F MAP = 2πa 2 1 (µ 1 µ 2 ) c ( µ 1 µ 2 H0 2 µ 2 )g, (19) L where µ 1 is the permeability of the particle, µ 2 is the permeability of the host medium, and g L is the longitudinal demagnetizing factor. From the above equation, we can find that the force remains proportional to the scalar product of particle volume and excess permeability. Ellipsoidal particles exhibit the same positive and negative magnetophoresis phenomenology, depending on whether µ 1 is greater than or less than µ 2. This expression is also valid for oblate spheroids (a < b = c), though such a particle with its short axis parallel to the magnetic field will be in unstable alignment. For a spherical particle, k = 1 and g L = 1/3, in which case Eq.(19) reduces to Eq. (4). In order to calculate the magnetophoretic force on the nonmagnetic particles, the lattice effect must be considered. So µ 2 should be substituted by µ e. In this paper, for convenience, F 1 stands for F MAP /(2πa 2 c H 2 0), which means that we only calculate the coefficient of the magnetophoretic force. Here, we assume the volume of the nonmagnetic particle is invariable, and H 2 0 is a constant across the nonmagnetic particles. For certain lattice, such as bct, bcc and fcc, the relation between F 1 and the aspect ratio k are shown in Fig.?? and Fig.??. The absolute value of the magnetophoretic force on the nonmagnetic particle decreases as k increases for both prolate and oblate cases. In addition, the relation between F 1 and the volume fraction of the nonmagnetic particles is also investigated. The results are shown in Fig.?? and Fig.??. The absolute value of the magnetophoretic force on the nonmagnetic particle decreases as the volume fraction increases, which are consistent with the result of the spherical case. We can calculate the transverse case by using the sum rules 2α T + α L = 3 and 2g T + g L = 1. There are no figures shown here for such case. III. CONCLUSIONS In fact, Magnetophoresis is a complicated phenomenon. Such as magnetophoresis with eddy current induction, which is corresponding to the electrically conductive magnetizable particle. In 8

17 our present work, we only study the magnetophoretic force on the nonmagnetic particles (or magnetic holes) suspended in a host ferrofluid, by taking into account the effect of structural transition and long-range interaction (lattice effect). So we resort to the Ewald-Kornfeld formulation and the well-known developed Maxwell-Garnett theory to calculate the effective permeability of the inverse ferrofluid. And we also consider the effects of the volume fraction and geometric shape of the nonmagnetic particles to the magnetophoretic force on the nonmagnetic particles. It is found that the absolute value of the magnetophoretic force on the nonmagnetic particle decreases as the volume fraction increases. For nonspherical case, the absolute value of the magnetophoretic force on the nonmagnetic particle decreases as the aspect ratio(k) increases. As the nonmagnetic particles are located closely, the multipolar interaction between them is expected to play a role [23, 24]. In this regard, it is of value to extend the present work to investigate the effect of the multipolar interaction. In addition, it is also interesting to see what happens if the external field is a nonuniform ac magnetic field and the nonmagnetic particle is inhomogeneous. To sum up, our results have shown that the magnetophoretic force on the nonmagnetic particle in the ferrofluid due to the nonuniform magnetic field can be affected significantly by the lattice effect, the volume fraction and shape of the nonmagnetic particle. Thus, it is possible to detect the structure of the inverse ferrofluid, volume fraction and shape of the nonmagnetic particle by measuring such magnetophoretic force. Acknowledgements The author Y. Gao greatly appreciated Prof. X. Sun for his generously help in his study. The authors acknowledge the financial support by the Shanghai Education Committee and the Shanghai Education Development Foundation (Shu Guang project) under Grant No. KBH , by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry, China, by the National Natural Science Foundation of China under Grant No

18 [1] T. B. Jones, Electromechanics of Particles (Cambridge University Press, Cambridge, 1995), Chap.III. [2] Maciej Zborowski, Graciela R. Ostera, Lee R. Moore, Sarah Milliron, Jeffrey J. Chalmers and Alan N. Schechter, Biophysical Journal 84, 2638 (2003) [3] A. T. Skjeltorp, Phys. Rew. Lett, 51, 2306 (1983). [4] J. Ugelstad et al, Blood Purif. 11, 349 (1993). [5] J. B. Hayter, R. Pynn, S. Charles, A. T. Skjeltorp, J. Trewhella, G. Stubbs, and P. Timmins, Phys. Rev. Lett. 62, 1667 (1989). [6] A. Koenig, P. Hebraud, C. Gosse, R. Dreyfus, J. Baudry, E. Bertrand, and J. Bibette, Phys. Rev. Lett. 95, (2005). [7] C.Scherer and A. M. Figueiredo Neto, Brazilian Journal of Physics, 35 3A, (2005). [8] Rosenweig, R. E. Ferrohydrodynamics (Cambridge University Press, Cambridge,1985). [9] Efraim Feinstein and Mara Prentiss, J. Appl. Phys. 99, (2006). [10] B. J. de Gans, C. Blom, A. P. Philipse, and J. Mellema, Phys. Rev. E 60, 4518 (1999). [11] Y. C. Jian, J. P. Huang and R. Tao, Structure of polydisperse inverese ferrofluids: Theory and computer simulation, (to be published). [12] R. Tao and J. M. Sun, Phys. Rev. Lett. 67, 398 (1991). [13] Lei Zhou, Weijia Wen, and Ping Sheng, Phys. Rev. Lett. 81, 1509 (1998). [14] J. P. Huang, Chem. Phys. Lett. 390, 380 (2004). [15] J. P. Huang, Phys. Rev. E 70, (2004). [16] P. P. Ewald, Ann. Phys. (Leipzig) 64, 253 (1921); H. Kornfeld, Z. Phys. 22, 27 (1924). [17] J. P. Huang and K. W. Yu, Appl. Phys. Lett. 87, (2005). [18] C.K. Lo and K.W. Yu, Phys. Rev. E 64, (2001). [19] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electrodynamics of Continuous Media, 2nd ed. (Pergamon, New York, 1984), Chap. II. [20] J. C. M. Garnett, Philos. Trans. R. Soc. London, Ser. A 203,385 (1904); 205, 237 (1906). [21] J. P. Huang, J. T. K. Wan, C. K. Lo, and K. W. Yu, Phys. Rev. E 64, (2001). [22] G. Wang, W. J. Tian and J. P. Huang, J. Phys. Chem. B 110, (2006) [23] K.W. Yu, J.T.K. Wan, Comput. Phys. Commun. 129, 177 (2000). [24] J.P. Huang, M. Karttunen, K.W. Yu, L. Dong, Phys. Rev. E 67, (2003). 10

19 Figure Captions Fig.1. (Color online) Cases of spherical particles. F as a function of q for different volume fractions. The external nonuniform magnetic field is parallel to the uniaxial anisotropic axis. Parameters: µ 1 = 1.0, µ 2 = 2.9. Fig.2. (Color online) Same as Fig. 1, but for transverse case. Fig.3. (Color online) Cases of prolate spheroidal particles. F 1 as a function of k for different latices. The external nonuniform magnetic field is parallel to the uniaxial anisotropic axis. Parameters: µ 1 = 1.0, µ 2 = 2.9, f = Fig.4. (Color online) Same as Fig. 3, but for oblate spheroidal case. Fig.5. (Color online) Cases of prolate spheroidal particles. F 1 as a function of f for different latices. The external nonuniform magnetic field is parallel to the uniaxial anisotropic axis. Parameters: µ 1 = 1.0, µ 2 = 2.9, k = 3.0. Fig.6. (Color online) Same as Fig. 5, but for oblate spheroidal case. 11

20 FIG. 1: Gao, Jian, Zhang and Huang 12

21 FIG. 2: Gao, Jian, Zhang and Huang 13

22 FIG. 3: Gao, Jian, Zhang and Huang 14