The dielectrophoretic gate and sorter system has been widely applied for preconcentrating and sorting of bioparticles for biodetection. In such systems, the dielectrophoretic force is generated by applying an AC electric field on the three dimensional electrode systems (containing a pair of electrodes on the top and bottom of the channel). Particles are held and sorted by balancing the DEP force with the hydrodynamic drag force. The holding capability is very important for such systems because it determines the preconcentration and sorting efficiency. In this paper, we investigate the holding capability of a simple dielectrophoretic gate system. Initially, a three dimensional numerical scheme was introduced to estimate the holding capability and was further validated by comparing with experimental results. Second, we systematically investigated the effects of the phase difference between the top and bottom electrodes; the height and width of the channel, and the relative position and size of top and bottom electrodes. The results demonstrated that the maximum holding capability is reached when the phase difference between the top and bottom electrodes is around 180 o. The results also show that the holding capability varied with the size and relative position of electrodes on the top and bottom, and the maximum holding capability is obtained when the top and bottom electrodes had the same size and the centers of both electrodes overlapped.
Rapid separation and detection of biochemical samples is necessary for most biological and chemical analysis. To achieve these processes, large laboratory equipment and big sample volumes are needed when using traditional techniques. In recent years, the rapid development of microfluidic systems has provided the technology to realize rapid separation and detection by integrating sample injection, mixing, separation, and detection onto a small chip [
In early DEP research, an array of interdigitated parallel electrodes was used as the typical design for dielectrophoretic separation. The sample suspension was quiescent over the electrodes on the bottom of the channel. When the AC electric potential was alternatively applied on the electrodes, a non-uniform electric field was established which induced a DEP force with polarizable particles. The DEP force caused the particles to move either toward or away from the electrodes depending upon the particle polarizability [9-11] . Subsequently, continuous flow DEP separation devices were developed for applications such as sorting particles [12-14] . In these devices, samples were injected into the channel with a constant flow rate, and an AC electric field was applied on the electrodes located on the bottom of the channel. Particles were separated by balancing the hydrodynamic stream force and the DEP force. Although planar electrodes on the bottom are simple and easily fabricated, some possible problems such as a much smaller DEP force and non-specific particle adhesion near the top of the channel caused low efficiency separation and trapping. To overcome these drawbacks, three dimensional (3D) electrode systems have been developed. The general idea of the 3D design is to fabricate a pair of electrodes which face each other in order to produce a stronger DEP trapping force. One electrode is on the top and the other electrode is on the bottom of the channel. Since the early 1990s, differrent shapes of 3D multi-electrodes to sort and manipulate particles have been explored [15-18] . Bennett et al. [
In such preconcentration and separation devices, particles were held due to the balance of the DEP force and the hydrodynamic force. However, particle would pass through the DEP gate if the hydrodynamic force is larger than the DEP force, leading a low efficiency on particle concentration and separation. Because the hydrodynamic force is linear with the sample flow rate, we introduced a threshold flow rate Qth, which is defined as the holding capability of the DEP force, to evaluate the holding and sorting efficiency. If the flow rate is bigger than Qth, the particle will escape and pass through the gate (bad holding and sorting efficiency). Conversely, if the flow rate is smaller than Qth, the particle will be held by the gate (good holding and sorting efficiency). Although the holding and sorting efficiency under 2D planar electrodes has been widely investigated [25-27] , the efficiency under a 3D electrode system has not been systematically studied and needs further investigation. Investigation of the critical voltage and critical velocity for sorting different size and types of particles has been performed using variable electrodes such as 3D parallel electrodes (alignment structure), 3D octopole electrodes (cages), and 3D curved electrodes [16-18] . Chen et al. [22,23] investigated the holding capability with different widths of the electrodes and different heights of the channel by using experimental and 2D simulation methods. However, important factors such as the relative position and size of the top and bottom electrodes, the width of the channel, and the phase difference between the top and bottom electrodes have not previously been investigated. In this paper, we present a systematical investigation of the holding capability of the DEP gate. First, a 3D numerical scheme is introduced to estimate the sholding capability and is validated by comparing to the experimental result. Second, the effects of the phase difference between the top and bottom electrodes, the height and width of the channel, and the relative position and size of electrodes is discussed in detail.
When a high frequency AC electric field is applied on the top electrode and the bottom electrode with a phase difference θ, a non-uniform electric field is generated within the channel. Particles flowing along the channel are polarized and a dipole moment is generated that interacts with the electric field. The time average DEP force for a single radius particle is determined by [
where a is the particle radius; is the vacuum permittivity; is the relative dielectric permittivity of the fluid, and ω is the angular frequency of applied electric field. is the phasor notation of the applied electric field; is the dipolar Clausius-Mossotti factor, and f is the frequency of the applied electric field. The first term on the right side depends on the spatially varying field magnitude which corresponds to the conventional dielectrophoretic (cDEP) force; the second term depends on the spatially varying phase, which corresponds to the traveling wave dielectrophoretic (twDEP) force. For a homoge-nous fluid medium, the electric potential satisfies the Laplace equation which is derived by the quasielec-trostatic form of Maxwell’s equations [
Using the Maxwell-Wagner expression for the complex dielectric permittivity, , where ε' is the permittivity, and σ is the conductivity. The polarizability, , can be expressed as follows [
The particle is attracted towards or repelled from the high electric-field regions depending on whether Re(λ) > 0 (positive dielectrophoretic force) or Re(λ) < 0 (negative dielectrophoretic force) with the range of Re(β) varying between 1 and −1/2. If the particle and the medium are not charged, λ is only a function of the frequency, f. There is a crossover frequency, fc , which indicates the change from positive to negative DEP. If the frequency f is below the crossover frequency, fc, the particle will undergo positive DEP (Re(λ) > 0), and if the frequency f is above the crossover frequency fc, the particle will undergo negative DEP (Re(λ) < 0). In this paper, DI water was used as the medium and 1 µm latex beads were used as the particles. The relative permittivity of the particle and the medium are 2.5 and 79 respectively. The conductivity of the DI water is 2 µS/cm. The crossover frequency is approximately 800 kHz for 1µm latex beads, and the calculated conductivity of the 1µm latex beads is 50 µS/cm.
If we assume the particles are neutrally buoyant and also neglect the gravitational effect, Brownian motion, and the electric-thermal effect, then the movement of the particles mainly depends on the hydrodynamic drag force and the DEP force by
where the hydrodynamic drag force in laminar flow is defined as, is the local flow velocity; is the particle velocity, and is the mass of particle. If the particle is held by the gate system, the velocity of the particle remains zero, and the hydrodynamic drag force and the DEP force are balanced. Because the flow direction is along the channel (z-direction), only the z-component of the DEP force needs to be considered. Equation (4) then becomes
To find the holding capability of the gate system, the maximum flow velocity needs to be found using equation (5). In the particular application, the flow rate Q is easier to control than the local fluid velocity in most applications. So, we introduced the maximum flow rate as the holding capability instead of the maximum flow velocity. In our investigation, the suspension of latex beads is driven through a rectangle channel using a syringe pump, and the stream flow along the channel (z-direction) is a fully developed pressure-driven Poiseuille flow. Thus, the distribution of the velocity is parabolic and the local velocity can be expressed as a classical Fourier series solution
where γn = (2n − 1) π/h and n is the series number. Uc corresponds to the fluid velocity in the center of the channel, and Rc(x, y) is defined as the position coefficient. The relationship between the flow rate Q and the flow velocity u(x, y) is calculated by
where is defined as the equivalent area of cross section.
A general model of the DEP gate is depicted in
The DEP gate consists of a pair of microelectrodes that span a fluid channel, perpendicular to the direction of fluid flow. Particles are injected using a syringe pump at a constant flow rate Q. When an AC voltage is applied to the microelectrodes, an electric field is generated within the fluid channel. The width of the channel is w, and the height of the channel is h.
tors were introduced to investigate the effect of electrode geometries and positions: one is the width/width ratio, , which is defined as the width of the top electrode divided the bottom electrode. And the other one is the translating rate, , which is defined as the distance between the centerline of two parallel electrodes divided by the width of the bottom electrode.
The 3D electric field was solved, and the DEP force was calculated by using the commercial software CFDACE + (ESI-CFD, Inc.). As an advanced technology, electric and electromagnetic simulation has been widely used in different research fields [30-32] .
For pressure driven flow, we calculated the maximum flow rate for particles held at position (x, y) by equation (7)
Here point A is called the critical point for particles located at the position (x, y). However, the maximum flow velocity um(x, y) or flow rate Qm(x, y) is related to the position (x, y) and only guarantees that particles will be held at the particular position. To find the holding capability of the DEP gate system, we need to obtain the maximum flow rate Qm(x, y) for each position (x, y). A critical surface is obtained by collecting the critical point for all x and y positions (a diagram of the edge of the critical surface is shown as the dashed line at
if the flow rate is not bigger than Qth, particles will be held at all the positions; otherwise, particles will pass through the channel at certain locations. We noticed that the shape of the contour surface varies case by case and is related to the distribution of DEP force and flow velocity. The zone where the minimum Qm(x, y) appears is known as the threshold zone (denoted with the ellipse at
Devices were fabricated in Sandia National Laboratories Microelectronics Development Laboratory using Sandia’s SwiFTTM (Surface micromachining With Integrated Fluid Technology) process (see [
Visualizations were made by using optical lenses (LMPlanFI 10×, 20×, 50×, Olympus), and movies were recorded using a CCD camera.
To validate the numerical method, we compared the predicted threshold flow rate Qth using the numerical method with the observed threshold flow rate Qth through experimental investigation. In the fabricated device, the height of the channel is 6 µm, the width of the electrodes is 10 µm, and the width of the channel is 200 µm. The numerical investigation has the same geometry and size. The applied electric field is 4 Volts peak to peak at 15 MHz, where Re(λ) = -0.47. The direction of the flow is from the left to the right. More and more particles are held and accumulated at the front of the DEP gate when the flow rate gradually increased from (20 to 40) µL/hour. The particles began to escape the DEP gate when the flow rate increased to 50 µL/hour.
However, the measured threshold flow rate determined by using this method (gradually increasing the flow until the particles began to escape the DEP gate) is much larger than the threshold flow rate predicted from the numerical simulation. A similar phenomenon was reported by [22,23]. The reason for the discrepancy is that more and more particles accumulating in front of the DEP gate cause strong inner-particle interactions. To avoid this phenomenon, another method is used for our experiments. By gradually increasing the flow rate and flushing out the trapped particles between each increase, we are able to determine the threshold flow velocity. The general idea is that we start from an initial flow rate which is so slow that the particles can be held by the DEP gate. Then, we turn off the electric field to flush away the accumulated particles. Next, we increase the flow rate by a constant small value to a higher flow rate and turn on the electric field again. If the particles are still held, we turn off the electric field and flush away the accumulated particles again, and increase the flow rate once more. This procedure is repeated until we observe that some particles escape from the DEP gate. The maximum flow rate at which particles can be held is measured as the threshold flow rate. In this method, the interaction between particles is almost eliminated because for every repeat most of the accumulated particles are flushed away before turning on the electric field.
Five different electric potentials (1, 2, 3, 4 and 5) Volts (corresponding to 2, 4, 6, 8 and 10 Volts peak to peak) were chosen to validate the numerical method.
As we know, the DEP force is directly affected by the applied electric field, which always includes the applied electric potential, electric frequency, and the phase difference between electrodes. Thus, the holding capability of the DEP gate system is sensitive to the electric field.
Although the effect of the applied electric potential and frequency for DEP gate and sorter systems has been investigated [16,18,22,23] , the effect of the phase difference between the top and bottom electrodes has not been studied. An AC electric field is always applied on one side of the electrode and the other side is grounded in most DEP gating systems. When we applied an AC electric potential with a phase difference on the other electrode, instead of the grounding the electrode, the electric strength varies with the phase difference. Meanwhile the traveling wave DEP force, FtwDEP, is generated and will also contribute to holding the particles. The holding capability may be enhanced due to these changes. In this section, we investigated the effect of the phase difference between the top and bottom electrodes (ranging from 0˚ to 360˚) on the holding capability of the DEP gate system. We applied the electric field with the amplitude 10 Volts peak-to-peak on both the top and bottom electrode while applying a phase difference θ between the top electrode and bottom electrodes. The frequency of the electric field is 2 MHz (Re(λ) = −0.3622 and Im(λ) = −0.3783), and the width of the channel is 200 µm. The numerical result shows that the traveling wave DEP force FtwDEP is much smaller than the conventional DEP force FcDEP (FtwDEP is only a ratio of 0.1% or less to FcDEP), which indicates the conventional DEP force makes the main contribution to the DEP gate. We also calculated and measured the threshold flow rate Qth, with phase differences between the top and bottom electrodes using numerical and experimental methods (the relationship is plotted in
The graph clearly shows the threshold velocity is a symmetric distribution around the phase difference θ = 180˚. When the phase difference θ is 0˚ (which means the electric field is the same on the top and bottom electrodes), particles will never be held because no DEP force is generated. When the phase difference increases up to 180˚, the threshold flow rate (holding capability)
also increased significantly. The threshold flow rate starts to decrease when the phase difference θ continuously increases from 180˚, and the value decreased back to 0 when the phase difference θ increased to 360˚. The result shows a maximum holding capability (which is approximately 180 µL/hour numerically and 190 µL/hour experimentally) was reached when the phase difference was 180˚. This is reasonable because the electric strength between the top and bottom increases when the phase difference changes from 0˚ to 180˚ and decreases when the phase difference changes from 180˚ to 360˚.
The height and width of the channel are also important factors for the DEP gating system because these parameters determine the distribution of the DEP force in the channel. In addition, the height and width of the channel need to be considered because the distribution of the velocity is also related to the height and the width of the channel (see Equation (6)). The flow rate Q is a function of the area of the cross section. To eliminate the effect of the area of the cross section on the flow rate, we investigate the holding capability by comparing the area-independent threshold flow rate qth, which is defined as the flow rate divided by the cross sectional area (qth = Qth/Wh). We keep the width of the bottom electrode d1 as 40 µm and only change the height and width of the channel. The applied electric field is 10 Volts peak-topeak at 15 MHz. Two other ratios γ = h/d1 and the aspect ratio κ = W/h were introduced.
We also investigated the effect of the aspect ratio. Four different values κ = 0.75, 1, 2 and 5 were considered.
In the above investigations, we used the same size for the top and bottom electrodes and also the centerlines of both electrodes are overlapped (α = 1 and β = 0). However, the size ratio of electrodes on the top and bottom and the relative position between the top and bottom electrodes will cause an asymmetric gradient of the electric field, which also affects the holding capability. To systematiccally investigate the effects of the size of electrodes and the relative positions between the top and bottom electrodes, the holding capability was investigated under different width/width ratio α and translating rate β using the validated numerical method. In the numerical simulation, the height of the channel h is 40 µm, and the width of the channel W is 200 µm. The width of the bottom electrode d1 is 40 µm. We considered different widths of the top electrode (α = 0.5, 1, 2, 3) and different translating distances of the top electrode away from the bottom electrode (β = 0, 0.5, 1, 2, 3). The applied electric field was 10 Volts peak-to-peak at 15MHz.
flow rate, Qth,, decreased with the increase of from 0.5 to 3. When β = 1, the largest threshold flow rate Qth corresponds to α = 2. The threshold flow rate decreased with the corresponding α = 2, 1, 3 and 0.5. When β = 0 was larger than 1 (β = 2, 3), the threshold velocity increased with the increase of α from 0.5 to 3. This means that the larger the top electrode; the larger the corresponding threshold velocity. Comparing the maximum threshold flow rate with different β indicates that the holding capability decreased from 1.28 µL/hour to 0.24 µL/hour when β increased from 0 to 3. Thus, the largest threshold velocity is 1.28 µL/hour and the corresponding parameters are α = 1 and β = 0, which indicates that electrodes equal in size with corresponding centerlines result in the largest holding capability.
It is also interesting to study the threshold zone because it is the weakest position though which the particle can pass. As shown in
positive value means the threshold zone is above the central plane and a negative value means the threshold zone is below the central plane. The figure shows the threshold zone is always below the center plane and near to the bottom electrode when α = 0.5 (the top electrode is smaller than the bottom electrode). For the case α = 1 (the top electrode is the same size as the bottom electrode), when β = 0, the threshold zone is located in the center plane (∆h = 0) because the gradient of the electric field is symmetric. However, the threshold zone is located above the central plane when β = 0.5, 1 and the threshold zone located below the central plane when β = 2, 3. For the case α = 2, the threshold zone is located above the central plane when β increases from 0 to 1 and is located below the central plane when β = 2, 3. For the case α = 3, the threshold zone is located above the central plane when β increases from 0 to 2 and is located below the central plane only for β = 3. Based on the analysis, we obtained a qualitative conclusion that the threshold zone is always located below the central plane when the size of the top electrode is smaller than that of the bottom electrode (α < 1), and the threshold zone is always located above the central plane when the size of the top electrode is larger than that of the bottom electrode (α > 1) for a low translating rate β. However the threshold zone will move under the central plane depending on the translating rate β.
In this paper, we studied the holding capability of a DEP gate or sorter system with a pair of parallel electrodes aligned on the top and bottom of the channel. To estimate the holding capability, a 3D numerical scheme was introduced and validated by comparing the experimental results at different applied electric potentials. This also indicates that the holding capability of DEP gate is linear with the square of applied electric potential and increases from approximately 4 µL/hour to 60 µL/hour when the applied electric potential increases from 2 Volt peak to peak to 10 Volt peak to peak. The effect of the phase difference between the top and bottom electrodes was investigated, and the result shows that the holding capability of the DEP gate increased when the phase difference between the top and bottom electrode increased from 0˚ to 180˚, and decreased when the phase difference changed from 180˚ to 360˚. A maximum holding capability (approximately 180 µL/hour in the numerical result and 190 µL/hour in the experimental result) was reached when the phase difference was 180˚. The effect of the height and width of the channel was also investigated. The result shows the holding capability was significantly increased when the height h decreased, especially smaller than the width of the electrode. However, the effect of the width of the channel is small and can be neglected when the height of the channel is larger than the width of the electrodes. Finally, the relative position and size of the top and bottom electrodes was studied. The results show the maximum holding capability is obtained when the top and bottom electrodes have the same width and the center lines of top and bottom electrodes are overlapped (α = 1 and β = 0). The results also indicate that the threshold zone will shift up or below the central plane when the gradient of electric field is asymmetric (α ≠ 1, β ≠ 0).
The authors would like to thank Conrad James at Sandia National Laboratories for his assistance regarding testing of the microfluidic devices. The authors would also like to thank Frances Ligler for her consultation in preparing this manuscript. We thank Sandia Laboratories for the use of the monolithic surface micromachined devices. D.B. thanks Sandia’s MESA Institute and the NSF MAGNET/SEM program for support. Sandia is a multiprogam laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy under Contract No. DE-ACo4-94-AL8500.