# Iterative Dipole Moment Method for the Dielectrophoretic Particle-Particle Interaction in a DC Electric Field

Authors: Qing Zhang; Kai Zhang
Journal: Journal of Nanotechnology
 (2018)

---

## Abstract

Electric force is the most popular technique for bioparticle transportation and manipulation in microfluidic systems. In this paper, the iterative dipole moment (IDM) method was used to calculate the dielectrophoretic (DEP) forces of particle-particle interactions in a two-dimensional DC electric field, and the Lagrangian method was used to solve the transportation of particles. It was found that the DEP properties and whether the connection line between initial positions of particles perpendicular or parallel to the electric field greatly affect the chain patterns. In addition, the dependence of the DEP particle interaction upon the particle diameters, initial particle positions, and the DEP properties have been studied in detail. The conclusions are advantageous in elelctrokinetic microfluidic systems where it may be desirable to control, manipulate, and assemble bioparticles.


---

## Body

## 1. Introduction

Microfluidic systems are widely used for biochemical analysis over the past decade, and electric forces are often used as an efficient and effective transport mechanism, which does not involve the intervention of moving parts and offers good manipulation over sample handling. Deposition of bioparticles and colloids onto the surface [1] is very important for many biochemical processes, such as assembly of carbon nanotubes [2, 3] and trapping of nanoparticles and nanoparticle synthesis [4–6]. In these processes, particles’ translation and rotation should be controlled accurately for counting and assembly of these particles [7, 8]. In microfluidic systems, electrical and hydrodynamic forces [9] dominate particles’ transportation. In the bulk fluid, particle transport is mainly affected by the hydrodynamic interactions and its neighboring particles.As one of the important forces of particle manipulation, dielectrophoretic (DEP) force has been drawing much attention in recent years. Equivalent dipole moment (EDM) has usually been used to calculate DEP force for single particle but cannot be used for the interaction between neighboring particles in an electric field. Maxwell stress tensor (MST) method [10] is complicated for the calculation of a large number of particles' interaction [11]. The iterative dipole moment (IDM) [12–14] method, which does not solve complicated differential equations of an electric field while the particle is moving, was used in this work for calculating the DEP forces of particle interactions.To control the assembling of bioparticles with high accuracy, an in-depth study is necessary to the underlying mechanisms of DEP particle-particle interaction in microfluidic systems. The impetus for the present study came from the importance in understanding the mechanics of particle assembling, which is often used to produce specialized material in a lab on a chip. In the present study, the DEP interaction of particle-particle [15] in a DC electric field will be studied, and a mathematical model based on the Lagrangian method incorporating DEP particle-particle interactions [16] will be presented to compute particle trajectories, and consequently, the assembling will be calculated. Finally, the dependence of particle assembling on the particle’s initial position and radius will be further concluded and discussed in detail.


## 2. Computational Model

As shown in Figure1, in a two-dimensional incompressible still media, DEP cylindrical particles are initially equispaced on a circle with radius R = 20 μm and an initial angle θ0 with respect to the applied electric field. Here, particle diameter d = 5 μm, for negative DEP, the permittivity of the media and particles are εm = 6.9e − 10 CV−1·m−1 and εp = 2.5 ∗8.8541878176e − 12 CV−1 · m−1; for positive dielectrophoresis, εm = 2.5 ∗8.8541878176e − 12 CV−1 · m−1 and εp = 6.9e − 10 CV−1· m−1.Figure 1
Cylindrical particles are equispaced on a circle in a two-dimensional DC uniform electric fieldE = 1 KV/m; here, red represents positive DEP and black represents negative DEP.The dielectrophoretic force,FDEP, acting on a spherical, homogeneous particle suspended in a local electric field gradient is given by the expression(1)FDEP=2πrp3εmReKw∇E2,where rp is the particle radius, εm is the permittivity of the suspending medium, ∇ is the Del vector (gradient) operator, E is the electric field incorporating additional fields due to particle interactions, and ReKw is the real part of the Clausius–Mossotti factor, which is given by(2)Kw=εp∗−εm∗εp∗+2εm∗,where εm∗ and εp∗ are the complex permittivities of the medium and particle, respectively, and ε∗=ε−jσ/w with σ is the conductivity, ε is the permittivity, and w is the angular frequency of the applied electric field. The limiting direct current (DC) case of the equation is(3)Kw=0=εp−εmεp+2εm.And the dielectrophoretic force can be modified as follows:(4)FDEP=4πrp3εmεp−εmεp+2εmE⋅∇E.Considering the influence of the dipole-induced field of other particles surrounding particle i, the modified electric field near the particle located atri is shown below:(5)Ei1r=E0r+∑j=1,j≠inEd0rj,r,i=1,2,3,…,N,where Ei1rdenotes the modified electric field and Ed0rj,r is the influence of the dipole-induced field of particle j; here,(6)Edx=−∂φd∂x=−Pi2πεmcosαiΔxi2+Δyi2−2ΔxiΔxicosαi+ΔyisinαiΔxi2+Δyi22,Edy=−∂φd∂y=−Pi2πε0sinαiΔxi2+Δyi2−2ΔyiΔxicosαi+ΔyisinαiΔxi2+Δyi22.The modified electric field induces a new dipole moment, which again induces a new electric field as follows:(7)Ei1ri→Pi1→Ed1ri,r,i=1,2,3,…,N,where Ei1ri, Pi1, and Ed1ri,r are the modified electric field, the modified dipole moment, and the modified dipole-induced electric field, respectively. These parameters can be iteratively modified until a converged value of the electric field Einri is obtained and the final dipole moment considering other particles’ influence can be achieved:(8)Pi=4πai3εmεi−εmεi+2εmEinri.These particles are assumed to be far from boundaries. The forces on particles are calculated by (1) while the dipole moment have been modified using IDM every time step, and particle transportation can be easily solved for particle trajectories and final particle chains.


## 3. Results and Discussion

Imagine a particle suspended in a dielectric fluid and subjected to a uniform electric field, which will polarize the dielectric particle and induce a dipole moment in it, then the dipole moment will induce an electrostatic potential as shown below:(9)φd=Picosθ4πεmr2.When particles are close to each other, one particle’s induced electrostatic potential will distort its neighbor’s electric field and make it nonuniform, and then the DEP force acting on its neighbor comes to a nonequilibrium state. The nonzero DEP forces on particles influencing each other suspended in a two-dimensional DC uniform electric field can be calculated by the IDM method.As shown in Figure2, if two heterogeneous DEP particles with same diameter are released at a small distance from each other and perpendicular to the electric field, they will attract each other and are likely to cluster at the center of line of their initial position. However, if all parameters except for the DEP property are the same, the results show that two homogeneous DEP particles will repel each other to move outwards symmetrically.Figure 2
Two particles (red or black dashed line) are initially equispaced on a circle (blue dashed line) withθ0 = π/2 in a uniform DC electric field with E = 1 kV/m along x direction, and then the final particle chains (solid particles) are, respectively, shown after DEP interaction under different statuses: (a) two negative (black) DEP particle, (b) one positive (red) and one negative DEP particle.
(a)
(b)From Figure3, it can be seen that if two homogeneous DEP particles with same diameter are released at a small distance from each other and parallel to the electric field, they will attract each other and cluster at the center of line of their initial position. However, if all parameters except for the DEP property are the same, the results show that two heterogeneous DEP particles will repel each other to move outwards symmetrically.Figure 3
The initial position of two particles are initially equispaced on a circle withθ0 = 0 in a uniform DC electric field with E = 1 kV/m along x direction, and the final particle chains after DEP interaction under different status: (a) two negative DEP particle, (b) two positive DEP particle, and (c) one positive and one negative DEP particle.
(a)
(b)
position: 227pt
(c)Figure4 shows the DEP interaction between two particles with different diameters. Figure 4(a) shows that when two heterogeneous DEP particles are released at a small distance from each other and perpendicular to the electric field, they will attract and cluster. However, as the difference of two particles’ diameters becomes larger, the position of final particle chains deviates from the center of line of their initial position, and the chains will move towards the side of the smaller particle. In addition, if all parameters except for the DEP property are the same, Figure 4(b) shows that the two homogeneous DEP particles with different diameters will repel each other to move outwards asymmetrically, and the smaller particle moves faster.Figure 4
The initial position of two particles is initially equispaced on a circle withθ0 = π/2 in a uniform DC electric field with E = 1 kV/m along x direction, and the final particle chains after DEP interaction under different status: (a) two heterogeneous particles with different diameters, (b) two negative DEP particles with different diameters.
(a)
(b)Figure5 shows that when two homogeneous DEP particles with different diameters are released at a small distance from each other and parallel to the electric field, they will attract and cluster, and the chains will move towards the side of the smaller particle.Figure 5
The initial position of two particles is initially equispaced on a circle withθ0 = 0 in a uniform DC electric field with E = 1 kV/m along x direction, and the final particle chains after DEP interaction under different status: (a) two negative DEP particle with different diameter and (b) two positive DEP particle with different diameters.
(a)
(b)From Figure6(a), it can be seen that if three homogeneous dielectrophoretic particles with same diameters are initially equispaced on a small circle, they will attract each other and finally cluster parallel to the electric field. However, from Figure 6(b), it can found that when one of the particles changed its DEP property, they still clustered but the chains were perpendicular to the electric field.Figure 6
The initial position of three particles is initially equispaced on a circle withθ0 = 0 in a uniform DC electric field with E = 1 kV/m along x direction, and the final particle chains after DEP interaction under different status: (a) three negative DEP particle and (b) two negative and one positive DEP particle.
(a)
(b)From Figure7(a), it can be seen that due to the symmetric configuration, homogeneous DEP particles 1 and 3 are attracted to move inwards because their connection line is parallel to the electric field; at the same time, positive DEP particles 2 and 4 are attracted by the chains of negative DEP particles 1 and 3 to move inward. While the connection line of the initial position of homogeneous particles 2 and 4 is perpendicular to the electric field, they cannot contact the chains of negative DEP particles 1 and 3. However, when the connection line between two homogeneous DEP particles 2 and 4 has a little shift θ0 = π/16 from the line perpendicular to the electric field, DEP interactions among particles cause four particles to cluster as shown in Figure 7(b).Figure 7
The initial position of two negative and two positive DEP particles are initially equispaced on a circle in a uniform DC electric field withE = 1 kV/m along x direction, and the final particle chains after DEP interaction under different status: (a) θ0 = 0, (b) θ0 = π/16.
(a)
(b)As shown in Figure8, three negative and two positive DEP particles are initially equispaced on a circle with θ0 = 0 in a uniform DC electric field. From the above conclusions, it can be found that the homogeneous adjacent particles whose connection line is not perpendicular to the electric field will cluster parallel to the electric field, and homogeneous adjacent heterogeneous particles will cluster perpendicular to the electric field, and eventually form the complex and asymmetric structure as shown in Figure 8.Figure 8
The initial position of five heterogeneous DEP particles is initially equispaced on a circle withθ0 = 0 in a uniform DC electric field with E = 1 kV/m along x direction, and the final particle chains after DEP interaction.According to these aforementioned conclusions, the regular aggregation patterns of large number of polarizable particles can also be well understood and can be used to control, manipulate, and assemble polarizable particles.


## 4. Conclusions

The IDM method was used to study multiple particle DEP interactions in a uniform DC electric field, and the DEP interactions among particles cause particles to cluster. It was found that the homogeneous adjacent particles whose connection line is not perpendicular to the electric field will cluster parallel to the electric field, and homogeneous adjacent heterogeneous particles will cluster perpendicular to the electric field; however, as the difference of two particles’ diameters becomes larger, the position of final particle chains deviates from the center of line of their initial position, and the chains will move towards the side of the smaller particle. In addition, the dependence of the DEP particle interaction upon the particle diameters, initial particle positions, and the DEP properties have been studied in detail. This can be advantageous in elelctrokinetic microfluidic systems where it may be desirable to control, manipulate, and assemble cylindrical bioparticles.



---

*Source: 3539075-2018-05-09.xml*