Prerpints.org logo
Preprint
Article

Modeling of Electric Field and Dielectrophoretic Force in a Parallel-plate Cell-separation Device with an Electrode Lid and Analytical Formulation Using Fourier Series

Altmetrics

Downloads

100

Views

51

Comments

0

This version is not peer-reviewed

Submitted:

11 October 2024

Posted:

11 October 2024

You are already at the latest version

Alerts
Abstract
Dielectrophoresis (DEP) cell-separation technology is an effective means of separating rare cells. To develop highly efficient cell-separation devices, detailed analysis of the nonuniform electric field’s intensity distribution within the device is needed, as it affects separation performance. Here we analytically expressed the distributions of the electric field and DEP force in a parallel-plate cell-separation DEP device by employing electrostatic analysis through the Fourier series method. The solution was approximated by extrapolating a novel approximate equation as a boundary condition for the potential between adjacent fingers of interdigitated electrodes and changed the underlying differential equation into a solvable form. The distributions of the potential and electric fields obtained by the analytical solution were compared with those from numerical simulations using finite element method software to verify their accuracy. The results showed excellent agreement. Three-dimensional fluorescence imaging analysis used live non-tumorigenic human mammary (MCF10A) cells. The distribution of cell clusters adsorbed on the interdigitated electrodes was compared with the analytically obtained distribution of the DEP force, and the mechanism underlying cell adsorption on the electrode surface was discussed. Furthermore, parametric analysis using the width and spacing of these electrodes as variables revealed that spacing is crucial for determining DEP force. The results showed that, by optimizing cell-separation devices using interdigitated electrodes, adjusting electrode spacing significantly enhances device performance.
Keywords: 
Subject: Engineering  -   Bioengineering

1. Introduction

In recent years, the development of a rapid and highly accurate method of separating and detecting specific rare cells from large numbers of diverse cell populations has become an urgent priority. For instance, the presence of circulating tumor cells (CTCs) in the blood increases as cancer progresses, underscoring the effectiveness of the rapid and accurate detection of CTCs for ultra-early cancer diagnosis. However, it is challenging to isolate and detect rare cells, including very early-stage CTCs, which exist in minute quantities in the blood. In light of these challenges, dielectrophoresis (DEP) technology has recently gained attention as an effective means of separating and detecting rare cells [1,2,3,4].
DEP is a phenomenon in which dielectric microparticles move along the gradient of a nonuniform electric field, propelled by the induction of electrical forces (DEP forces) acting on the microparticles. This results from the interaction between the dipole moment induced in the microparticles and the nonuniform electric field. In particular, using DEP with an AC electric field for cell separation offers significant advantages: it allows the separation of biological particles such as cells, viruses, and DNA with high accuracy, and is non-invasive to such particles because it does not require pretreatments such as fluorescent labeling [5,6]. However, planar electrodes, which are widely used in DEP cell-separation technology, have a limited effective volume for processing cell samples, posing a challenge in achieving highly efficient separation of large quantities of cells [7]. Along with the increasing size of cell-separation devices and the optimization of the parallel-plate microfluidic channel structure, which has been the basis of the cell-separation device structure, devices with various shapes of microfluidic channels and electrodes have been proposed to enable the highly efficient separation of large numbers of cells [8,9,10,11,12,13,14,15,16,17,18,19,20]. However, issues such as increased applied voltage due to larger and more complex device structures have raised concerns about the adverse effects stemming from the generation of Joule heat and strong AC electromagnetic fields on the physiological functions of cells [21]. Therefore, there is a need to develop devices capable of highly efficient cell separation even under low-voltage loading conditions. Since the magnitude and distribution of the DEP forces acting on cells depend heavily on the electric field distribution, the precise determination of the electric field distribution will enable the optimal design of the device’s internal structure, electrode arrangement, and electrode architecture. This, in turn, is expected to facilitate the development of highly efficient cell-separation devices that can operate under low-voltage conditions.
Numerous studies have aimed to rigorously determine the electric field in microfluidic devices and the DEP force induced by nonuniform electric fields. A wide range of methods exists for obtaining electric field distributions. For example, analytical approaches have utilized Fourier series [18,22,23,24,25], Green’s functions [26,27], and Green’s theorem [28,29], while other approaches have been based on the principle of conformal mapping, including the Schwarz‒Christoffel transform [30,31,32,33,34], or on numerical simulations using finite element methods (FEMs) [35,36,37,38,39].
In this study, we derived analytical solutions for the distribution of the AC electric field and DEP forces generated in a parallel-plate cell-separation device that we proposed previously [20,40] by applying electrostatic analysis using the Fourier series method. In modeling the proposed device, whose top plate (lid) also serves as an electrode, we need to solve a mixed boundary value problem comprising three Dirichlet boundary conditions and one Neumann boundary condition for the electric potential, excluding the periodic boundary condition [41]. However, this problem differs significantly from conventional devices that have either a glass top plate or open-top interdigitated electrodes, where exact solutions are attainable, highlighting the severe challenge of obtaining closed-form solutions [12,18,29,31,33,42]. Therefore, in this study we extrapolated a newly proposed approximate equation to the potential between adjacent fingers of interdigitated electrodes, originally necessitating a Neumann boundary condition. This was replaced with the corresponding Dirichlet boundary condition to reformulate the underlying differential equation into a solvable form, thus enabling the derivation of an analytical solution. The accuracy of the distributions of the potential and electric fields was verified by comparing those obtained by the analytical solutions to the results of numerical simulations using commercially available FEM software. Three-dimensional (3D) fluorescence imaging analysis was subsequently performed using live non-tumorigenic human mammary cells (MCF10A). The analytical model was validated by comparing the distribution of live cell clusters adsorbed on the interdigitated electrode surface to the analytically obtained distribution of the DEP force. The mechanisms underlying cell adsorption and cell cluster formation at the interdigitated electrode surface were also discussed. Furthermore, a parametric analysis was performed using the electrode width and spacing of the interdigitated electrodes as variables, and the optimal electrode arrangement and architecture for highly efficient cell separation were explored and proposed.

2. Theory

2.1. Basic Principle of DEP

The DEP force induced in the cell subjected to a nonuniform AC electric field is given as follows,
F D E P = 2 π   ϵ l   d c   2 3 R e β   ∇ E 2
where ϵ l is the electric permittivity of the solution, d c is the cell diameter, ∇ is the Nabla operator, and E is the AC electric field vector. R e β denotes the real part of the Clausius–Mossotti (CM) factor
β = ϵ c * − ϵ l *   ϵ c * + 2   ϵ l *
where ϵ c is the electric permittivity of the cell. The symbol ‘ * ’ represents the complex electric permittivity and is expressed as
ϵ * = ϵ + σ   j ω   ,     j = − 1
where σ is the electric conductivity and ω is the angular frequency of the AC electric field ( = 2 π f ; f is the electric field frequency). For the application of the single-shell model to live mammalian cells, β can be expressed as [43]
β =   ω 2 τ l   τ c * − τ c   τ l * + j ω   τ l * − τ l −   τ c * − 1     ω 2 2 τ l   τ c * − τ c   τ l * − j ω   τ l * − 2 τ l −   τ c * − 2  
where
  τ c * =   c m d c       2 σ c ′     ,     τ c = ϵ c ′   σ c ′     ,       τ l * =   c m d c     2 σ l   ,     τ l =   ϵ l     σ l   .
Here, σ l is the electric conductivity of the solution, ϵ c ′ is the effective electric permittivity of the cell, σ c ′ is the effective electric conductivity of the cell, and c m is the capacitance of the cell membrane. As an example, Figure 1 shows the frequency spectra of the R e ( β ) of live MCF10A cells for various values of σ l . Values of the d c and other electrical properties were taken from Ref. [44].
To evaluate the magnitude of F D E P induced on the cell, it is necessary to know the values of R e ( β ) and ∇ E 2 . The value of R e ( β ) can be obtained from Equation (2) if the values of d c and f are known, and the value of ∇ E 2 can be determined if the distribution of E can be obtained by some means.

2.2. Electric Field Analysis

Figure 2a shows a schematic of the cell-separation device to be analyzed. The device has a parallel-plate microfluidic channel structure consisting of interdigitated electrodes on the bottom surface and a planar electrode on the top surface. By applying AC voltage to the device, a nonuniform electric field is generated across the whole volume of the microfluidic channel. Because the fingers of the interdigitated electrodes are sufficiently long and aligned along the flow direction of the microfluidic channel, the electric field analysis model in the device can be simplified into a two-dimensional (2D) model in the cross section of the microfluidic channel, as shown in Figure 2b.
The mathematical model of the electric field analysis is shown in Figure 3. In the present analysis, an electrostatic field analysis based on the AC effective value was adopted. The governing equation for the potential of the electrostatic field in the microfluidic channel, φ x , y , is the Laplace equation
  ∂ 2 φ x , y   ∂ x 2 +   ∂ 2 φ x , y   ∂ y 2 ＝ 0   .
The boundary conditions for the potential are
φ x , H = V r m s             0 ≤ x ≤ w
∂ φ x , H ∂ y = 0             w < x < L − w
φ x , H = 0               L − w ≤ x ≤ L
φ x , 0 = 0
∂ φ 0 , y ∂ x = 0
∂ φ L , y ∂ x = 0
where V r m s is the effective value of the applied voltage, L is the distance between the center lines of adjacent interdigitated electrodes, H is the height of the microfluidic channel, and w is half the width of the electrode. Here, the following dimensionless variables are introduced.
x * = x   L        
y * = y   H     .
Substituting Eqs. (10) and (11) into Eqs. (3) − (9) yields
∂ 2 φ x * , y * ∂ x * 2 + L   H   2 ∂ 2 φ x * , y * ∂ y * 2 = 0
φ x * , 1 = V r m s               0 ≤ x * ≤ w   L  
∂ φ x * , 1 ∂ y * = 0               w   L   < x * < 1 − w   L  
φ x * , 1 = 0               1 − w   L   ≤ x * ≤ 1
φ x * , 0 = 0
∂ φ 0 , y * ∂ x * = 0
∂ φ 1 , y * ∂ x * = 0
Equations (12) − (18) collectively form a mixed boundary value problem. It is particularly difficult to obtain a closed-form solution in the case of this analysis. Therefore, to obtain a solution by using the quadrature method, the Neumann boundary condition in Equation (14) was replaced by a Dirichlet boundary condition
φ x * , 1 = V x *               w   L   < x * < 1 − w   L     .
For an extrapolating function of V x * on the right-hand side of Equation (19), the expression
V x * =   V r m s   π cos − 1 ⁡ x *
was chosen. Equation (20) closely approximates the exact solution for the inter-electrode potential in the absence of a planar electrode on the top surface of the device ( y * = 0 ) [Appendix A]. To make Equation (20) integrable in the Fourier series method, the third-order terms of the McLaurin expansion of cos − 1 ⁡ x * were considered. This was done under the conditions where 1 − 2 w / L ≅ 0 , resulting in the reduction of Equation (20) as
V x * ⋍   V r m s   π π   2   − x * −   x * 3   6 + O x * 5   .
Equation (21), satisfying boundary conditions, V w / L = V r m s a n d V 1 − w / L = 0 , was obtained as
V x * =   V r m s   π π   2   − C x * − 1   2     1   2   − w   L     − π   2   − C x * − 1   2   3   6 1   2   − w   L   3     .
Here, C is an arbitrary constant, and its optimal value was determined by curve fitting Equation (22) to Equation (20) using the method of least squares, yielding C = 0.75480 . Figure 4 shows the distribution of Equation (22) normalized to the value V r m s for various values of C when H / L = 5 and w / L = 0.25 , together with the distribution of normalized Equation (20).

2.3. Analytical Solution by Fourier Series Method

The solution of φ = φ x * , y * is expressed as
φ = A ′ + B ′ x * + C ′ y * + D ′ x * y * + f x * g y *
where A ′ , B ′ , C ′ , and D ′ are unknown constants, and f and g are functions of only x * and y * satisfying Eqs. (13) − (18), respectively. Substituting Equation (16) into Equation (23) yields
A ′ = B ′ = 0
Furthermore, substitution of Eqs. (17) and (18) into Equation (23) yields
D ′ = 0   .
Using Eqs. (13), (15), (19), and (22), the functional form of φ is determined as follows, assuming k n = n π ( n = 1 ,   2 ,   ⋯ ) [Appendix B]
φ = a 0 y * + 2 ∑ n = 1 ∞ a n cos ⁡ k n x * sinh ⁡   k n H   L y *
a 0 =   V r m s   2
a n =   V r m s   k n 1 −   2 C     π   − 4 C   π k n   L − 2 w     L     tan ⁡   k n   2   L − 2 w     L   − 24 1 −   2 C     π     k n 2   L − 2 w     L   3     L − 2 w     L     k n 2   24   L − 2 w     L   2 − 1 + 2   k n     k n 2   8   L − 2 w     L   2 − 1 cot ⁡   k n   2   sin ⁡   k n   2 cos ⁡   k n   2   L − 2 w     L     sinh ⁡   k n H     L   #

2.4. Analytical Form of DEP Force

Since the electric field E = E x * , y * is a vector quantity, it can be represented by the vector expression as
E = E x e x + E y e y = − g r a d φ
where e x and e y denote the unit vectors in the x and y directions, respectively, and E x = E x x * , y * and E y = E y x * , y * are the x and y components of the electric field E . Using Eqs. (26) − (28), E x and E y become
E x = 2 ∑ n = 1 ∞ a n k n sin ⁡ k n x * sinh ⁡   k n H   L y *                            
E y =   V r m s   2 −   2 H   L ∑ n = 1 ∞ a n k n cos ⁡ k n x * cos h ⁡   k n H   L y *   .
To find F D E P , we can utilize the mathematical relationship
∇ E 2 = ∇ E x 2 + E y 2                                                                                                                                                                                       = 2 E x   d E x   d x * + E y   d E y   d x * e x + 2 E x   d E x   d y * + E y   d E y   d y * e y
Equation (1) becomes
F D E P = 4 π   ϵ l   d c   2 3 R e β E x   d E x   d x * + E y   d E y   d x * e x + E x   d E x   d y * + E y   d E y   d y * e y   .
F D E P can be found by substituting Eqs. (30) and (31) into Equation (33).

3. Experimental

3.1. Cell Sample Preparation

Live non-tumorigenic human mammary (MCF10A) cells were used as cell samples. A culture medium consisting of 500 mL D M E M / F 12 (11320033, Gibco), 20 ng/mL Human epidermal growth factor (E9644, ThermoFisher), 0.5 μg/mL Hydrocortisone (H4001, Sigma Aldrich), 100 ng/mL Cholera toxin (C8052, Sigma Aldrich), 5% (v/v) Horse serum (16050122, ThermoFisher), and 10 μg/mL Human insulin (I9278, Sigma Aldrich) was used for cell cultivation. The cultured cells were incubated in a 5% CO2 incubator for 2–3 days. Experiments were conducted using cells at passages P5–10.

3.2. Preparation of Cell Sample Solution

A cell sample solution was prepared by suspending cells in a 300 mM isotonic mannitol solution to achieve a cell density of ϕ = 8.0 × 10 5 cells/mL. For fluorescence imaging observations, cells were incubated with Calcein AM (C396, Dojindo; λ e x = 490 nm, λ e m = 515 nm). Calcein AM was introduced into the cells by incubating them in a PBS( − ) solution with a Calcein AM concentration of 20 µM for 30 min. The conductivity of the cell sample solution was adjusted to σ f = 1.0 × 10 − 3 S/m using the culture medium.

3.3. Fabrication of Cell-Separation Device

A cell-separation device using the counter-interdigitated electrodes was fabricated, and the relationship between the distribution of cell clusters formed on the electrode surface and the DEP force distribution determined through analytical solutions was examined.
Three types of interdigitated electrodes with different electrode widths and spacings were fabricated to generate different electric field distributions under the condition of a constant applied voltage. Figure 5a shows the dimensions of the electrode substrate used in the experiment, and Figure 5b shows the dimensions of the three types of interdigitated electrode. For the electrode width or spacing, combinations of 50 μm and 75 μm were considered. The interdigitated electrodes with (width : spacing) ratios were designated as Type A (50 μm : 50 μm), Type B (50 μm : 75 μm), and Type C (75 μm : 50 μm). A standard photolithographic method was used to fabricate the interdigitated electrodes. The procedure is briefly described as follows. First, the surfaces of the glass planar plate with dimensions of 50 × 90 mm and a thickness of 0.4 mm were ultrasonically cleaned with acetone and isopropyl alcohol for 5 min each, followed by ultrapure water for 5 min. Next, a 300-nm-thick aluminum film was deposited on the glass plate surface by vacuum evaporation. Electrode patterning was achieved using a photoresist (S1805G). Pre-baking was performed at 90 ℃ for 3 min, followed by a 20 s exposure and 1 min of development. Post-baking was performed at 130 ℃ for 3 min to fix the electrode to the glass plate surface. Then, etching treatment with mixed acid was performed at 40 ℃ for 90 s, and the photoresist was removed using an AZ 100 remover (AZ Electronic Materials). The fingers of the interdigitated electrodes were arranged in alternating fashion with high-voltage electrodes and grounded electrodes featuring 75 pairs for Type A and 60 pairs for Types B and C.
The cell-separation device has a parallel-plate microfluidic channel structure with electrode substrates on the upper and lower surfaces, as shown in Figure 2. The upper electrode substrate is a 1.1-mm-thick glass planar plate coated with an indium‒tin oxide (ITO) film on its surface, while the lower electrode substrate is a glass planar plate printed with interdigitated electrodes. Holes with a diameter of ϕ 1.0 mm were drilled in the upper electrode substrate using a leutor (2307396, Sea Force) for the introduction and drainage of the cell sample solution. Polypropylene female luer fittings were bonded to the drilled holes, and a conductive epoxy resin adhesive (CW2400, Chemtronics) was used to attach lead wires to the terminals of each electrode substrate. Finally, a 0.5-mm-thick silicone rubber spacer was placed between the upper and lower electrode substrates and compressed to complete device assembly.

3.4. Experimental Apparatus and Method

The experimental setup consisted primarily of the cell-separation device, a syringe pump, a waveform generator, and a confocal laser scanning microscope, as shown in Figure 6. To observe the behavior of cells in the device, the electric field frequency, f , that maximizes the value of R e β in Equation (1) was calculated using Equation (2) to maximize the DEP force acting on the cells. In the case of MCF10A cells, when the solution conductivity is   σ f = 1.0 × 10 − 3 S/m, the frequency that gives the maximum value of the R e β is   f ≅ 8.5 MHz. The value of R e β is R e β ~ 0.98 . From the above, the applied voltage was set as V = 10 Vpp with f = 8.5 MHz.
The experimental procedure was as follows. First, 2 mL of the cell sample solution was introduced into the device using the syringe pump, and the flow was stopped after the degassing was confirmed. Next, the waveform generator was turned on to produce a nonuniform electric field in the device. After confirming the stability of cell distribution in the device through bright-field observation under the microscope, the z -scan function of the microscope was used to capture a series of slice fluorescence 2D images of the cell clusters. A 3D image was then constructed using the accompanying image analysis program to examine the distribution of cells adsorbed on the surface of the interdigitated electrodes.

4. Results and Discussion

4.1. Accuracy Verification of the Analytical Solution

To verify the accuracy of the analytical solution obtained, a numerical simulation was performed using commercially available FEM software, FEATOOLS (https:// www.featool.com), and the results were compared. The boundary conditions for the numerical simulation model were the same as those indicated in Figure 3. The numerical simulation was performed using approximately 2.6 × 10 6 grid points and with a relative error of 1 × 10 − 6 for the convergence criterion of the iterative calculation. If the number of grid points exceeded 2.6 × 10 6 , the average relative error of the electric field distribution was less than 0.15% compared to, for example, results obtained with 1 × 10 7 grid points, except for disparities in the electric field magnitude at the electrode edges ( x * = w / L , 1 − 2 w / L ) . To compute the analytical solution, FORTRAN code was written in-house, and the number of terms in the series expansion, n , was set to n = 10000 . The aspect ratios of the calculation domain were set to match those of the Type A device, H / L = 5 and w / L = 0.25 . For n > 10000 , the error of the solution for n = 20000 , for example, relative to that of n = 10000 , was ~ 3.0 × 10 − 10 %. Thus, n = 10000 gives sufficient accuracy for the solution. Figure 7 compares the analytical solution and numerical simulation for distributions of φ and E . The effective potential value, V r m s , was set to 5 / 2 V, matching the applied voltage adopted in the experiment, and φ and E were each normalized to φ / V r m s and E H / V r m s , respectively. For φ * and   E * at the interdigitated electrode surface y * = 1 , the average relative errors, ∆ φ * and ∆ E * , between the numerical simulation and the analytical solutions were 6.24 % and 6.29 %, respectively. These results showed good agreement, affirming that the analytical solutions provided reasonable outcomes. Figure 8 compares the analytical solution and numerical simulation results for the distribution of ∇ E * 2 along the bottom surface of the microfluidic channel. The analytical solution overestimated the magnitude of ∇ E * 2 in the inter-electrode region 0.25 < x * < 0.75 by about 10 % compared to the numerical simulation results. On the other hand, in the region from the electrode edge, where the value of ∇ E * 2 is at its maximum, extending across the entire electrode surface, the two were in relatively good agreement. From Equation (1), the magnitude of the DEP force is proportional to that of the gradient of the square of electric field, written as F D E P = F D E P ~ ∇ E 2 . Thus, the results demonstrated that the analytical solution of ∇ E * 2 gives a reasonable distribution of F D E P . While examining the inter-electrode distribution of ∇ E * 2 in the numerical simulation, despite utilizing a sufficient number of grid points ( 2.6 × 10 6 ), strong spatial oscillations occurred, preventing the correct distribution from being obtained. This occurs due to its
on the interdigitated electrode surface y * = 1 .
proximity to the two electrode edges, which align with the mathematical singular points, making the numerical integration extremely unstable. For example, the numerical integration generally gives a spatially continuous solution for F D E P , a function of the second-order derivative of the potential, but it does not always guarantee the smooth continuity of the solution. The effect is particularly pronounced in regions where local values change rapidly, such as positions near the singular point. The solution either has strong spatial oscillations or diverges without convergence. Spurious oscillations in the distribution of F D E P are inevitable even when other commercial FEM software is used [22,45]. On the other hand, the analytical solution is particularly useful for examining the force field in the DEP device, because it always gives a smooth and stable solution for physical quantities expressed in terms of the second- or higher orders derivative of φ , such as ∇ E 2 .

4.2. Distribution of DEP Force

Figure 9 shows the distributions of F D E P   = F D E P induced in cells by means of the three types of devices, Types A, B, and C. Note that F D E P is normalized as F D E P * = F D E P / F R . The F R represents the magnitude of the DEP force, defined as F R = 2 π   ϵ l d c / 2 3 R e β   φ 2 / H 3 . The analysis conditions are listed in Table 1. For all types of interdigitated electrodes, an extremely strong force field was created around the electrode edges. The strength of F D E P * was significant at the electrode surface and hardly decayed from the electrode surface to a height of ~ 20 µm, equivalent to the size of one or two cells. The value of F D E P * began to decay rapidly as the distance from the electrode surface increased, occurring at an approximate height of over 20 µm. At a height of ~ 30 µm, the value of F D E P * decayed by two to three orders of magnitude. The F D E P * reached a minimum value at the midpoint between two adjacent electrodes because the x * axis components of F D E P * generated by these two electrodes acting on the cell at this position canceled each other out. Therefore, F D E P * is an even function with respect to the vertical axis passing through this center position, leading to F D E P * satisfying ∂ 2 F D E P * / ∂ x 2 = 0 . Figure 8 shows that the analytical solution satisfies this condition ( ∂ 2 l o g ∇ E * 2 / ∂ x * 2 = 0 ) at x * = 0.5 , but the numerical simulation does not produce the correct distribution due to spatial oscillations.

4.3. Experimental Results

Figure 10 shows fluorescence images of the distribution of cells around the interdigitated electrodes of Types A, B, and C. For all three types, it can be seen that the cells were adsorbed along the electrode edges. This is because the F D E P induced in the cells at each electrode edge was extremely large, and the cells that approached the electrode were
attracted to this position first and then adsorbed. In addition, the cause of the low adsorption of cells to the central area of the electrode surface is that F D E P was weaker there because E x = 0 . Comparing the distribution of F D E P in Figure 9 with the distribution of cells in Figure 10, it can be seen that they are generally similar to each other. When Type A and Type C are compared, no significant difference in the cell distribution is found. This is likely due to the similarity in the F D E P distribution within the inter-electrode region between the two types, as illustrated in Figure 9. Next, when comparing Types A and B, Type A showed cell adsorption in the inter-electrode region, while Type B showed almost no cell adsorption in the inter-electrode region. This is because, as shown in Figure 9, F D E P between electrodes was weaker in Type B than in Type A. As indicated above, cells tend to adsorb in the inter-electrode region rather than on the electrode surface, highlighting that the amount of cell adsorption, which determines the device’s cell-separation performance, is affected more by adjusting the electrode spacing than by changing the electrode width.

4.4. Relationship between DEP Force and Electrode Width and Spacing

To investigate the effect of electrode architecture on DEP force, electrode width ( 2 w ) and electrode spacing
d = L − 2 w
were used as variables in the parametric analysis. The physical properties used in this analysis are shown in Table 1. Figure 11 shows a bird’s-eye view of the change in the average value of F D E P on the interdigitated electrode surface as d and w were changed independently from 5 µm to 60 µm in 200-point increments. As the figure shows, the value of F D E P changed more sensitively in response to changes in d than in response to changes in w . In addition, considering the connection between the distribution of F D E P and that of w , it is thought that optimizing a device for highly efficient cell separation entails electrode spacing ( d ) of 30–50 µm. This design accounts for the sizes of the cells, particularly for cells of a size similar to MCF10A cells ( d c ~ 15 μm) used in this study.

5. Conclusions

To investigate the distributions of the electric field intensity and DEP force within the newly proposed parallel-plate cell-separation device, we obtained analytical solutions using the Fourier series method. By extrapolating the inter-electrode potential with a newly proposed approximation, the boundary condition, initially given in Neumann type, was replaced by Dirichlet type. This adjustment allowed for the integration of the governing equations to derive an approximate solution. To verify the accuracy of the obtained solution, a comparison with numerical simulations was performed, demonstrating excellent agreement for both the potential and electric field intensity distributions. In addition, 3D fluorescence imaging analysis of the cell distribution in the cell-separation device was performed. The distribution of live cells near the interdigitated electrode surface was compared with that of the DEP force obtained by the analytical solution. This comparison revealed that the distribution of live cells was very similar to that of the DEP forces. In addition, a parametric analysis involving the electrode width and spacing as variables revealed for the first time that electrode spacing was the key parameter influencing the level of cell adsorption onto the electrode. In cell-separation processes utilizing devices with interdigitated electrodes, electrode width is not very important for highly efficient cell separation. For example, when separating human cells with a diameter of about 15 μm, electrode spacing of 30–50 μm is desirable.

Author Contributions

Conceptualization, D.N. and S.T.; methodology, D.N. and S.T; software, D.N. and S.T.; validation, D.N.,Y.S. and S.T.; formal analysis, D.N. and S.T.; investigation, D.N.,Y.S. and S.T.; resources, S.T.; data curation, D.N. and S.T.; writing—original draft preparation, D.N.; writing—review and editing, S.T.; visualization, D.N. and S.T; supervision, S.T.; project administration, S.T.; funding acquisition, S.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by JSPS KAKENHI Grant Number JP23K03669.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be provided on suitable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The exact solution for the inter-electrode potential in the absence of the planar electrode on the top surface   y * = 0 is expressed by an incomplete elliptic integral of the first kind with an elliptical modulus k ( = ( 1 − 2 w ) ⁄ L ) [33]
F λ , k = ∫ 0 λ d θ   1 − k 2 sin 2 ⁡ θ  
Since the integration of F ( λ , k ) by the quadrature method is difficult, we will seek an approximate solution in integrable form. The function F ( λ , k ) can be expressed by using the inverse function of Jacobi’s elliptic function, S n , as
F λ , k = S n − 1 ( s i n λ ,   k )
If λ = sin − 1 ⁡ x ′ , then for k → 0 , i.e., 1 − 2 d / L → 0 , we have
S n x ′ ,   0 ≅ s i n x ′
Therefore, Equation (A2) can be approximated as
F λ , k ≅ sin − 1 ⁡ x ′
Equation (20) is obtained by performing the variable transformation x ′ → x * in Equation (A4) and applying the appropriate similarity transformation to satisfy the boundary conditions Eqs. (13) and (15). To evaluate the accuracy of Equation (20), the inter-electrode potential obtained using it was compared with numerical simulation results obtained using commercial FEM software. The distribution of φ = φ x * , 1 obtained from the exact solution of Equation (A2) and that obtained from Equation (20) are shown in Fig. A1. For comparison, distributions of φ
Preprints 120926 g0a1
Figure A1. Inter-electrode potential distribution.
Figure A1. Inter-electrode potential distribution.
Preprints 120926 g0a1
obtained by numerical simulations with three different aspect ratios, H / L , in Figure 3 are also shown. The results of the exact solution and Equation (20) show remarkable agreement. Equation (20) gives a good approximation of the inter-electrode potential when H / L > 5 , exhibiting relative errors within 5% of the commercial FEM software results (not shown).

Appendix B

When the Laplace equation, Equation (3), is solved analytically through the separation-of-variables technique, the general form of f x * g y * on the right-hand side of Equation (23) is expressed as
f x * g y * = ∑ n = 1 ∞ a n cos ⁡ k n x * + b n sin ⁡ k n x * ( c n e k n y * + d n e − k n y * )
where a n , b n , c n , d n and k n are arbitrary coefficients. From Eqs. (17) and (18)
b n = 0
k n = n π
Furthermore, from Equation (16)
c n = − d n
and therefore, Equation (B1) becomes
f x * g y * = 2 ∑ n = 1 ∞ a n cos ⁡ k n x * sinh ⁡ k n H L y *  
Using Eqs. (13), (15), (19), and (22), the functional form of φ x * , y * is determined as
φ x * , y * = a 0 y * + 2 ∑ n = 1 ∞ a n cos ⁡ k n x * sinh ⁡ k n H L y *
where a 0 and a n are the Fourier coefficients. Because φ x * , 1 is a periodic function of period x * = 1 , the Fourier coefficient a 0 can be expressed as
a 0 = ∫ 0 1 φ x * , 1 d x *
Because the right-hand side of Equation (B7) is equal to the area of the figure bounded by the curve of potential φ x * , 1 and the x * and y * axes shown in Figure 4, a 0 can be readily obtained as
a 0 = V r m s w   L   + 1   2   V r m s 1 − w   L   − w   L   =   V r m s   2
For a n , if y * = 1 in Equation (B6), further multiplying both sides by cos ⁡ k m x * and integrating over the interval 0,1 yields
∫ 0 1 φ x * , 1 cos ⁡ k m x * d x * = ∫ 0 1 a 0 cos ⁡ k m x * d x * + 2 ∫ 0 1 ∑ n = 1 ∞ a n cos ⁡ k n x * sinh ⁡ k n H L cos ⁡ k m x * d x *
Due to the periodicity and orthogonality of trigonometric functions
∫ 0 1 a 0 cos ⁡ k n x * d x * = 0 2 ∫ 0 1 ∑ n = 1 ∞ a n cos ⁡ k n x * sinh ⁡ k n H L cos ⁡ k m x * d x * = 2 a n sinh ⁡ k n H L ∫ 0 1 cos 2 ⁡ k n x * d x * = a n sinh ⁡ k n H L   #
Since it is required that a n is
a n = 1 sinh ⁡ k n H L   ∫ 0 1 φ x * , 1 cos ⁡ k n x * d x *  
substituting Eqs. (13), (15), (19), and (22) into Equation (B11) yields Equation (28).

References

  1. Farasat, M.; Chavoshi, S.M.; Bakhshi, A.; Valipour, A.; Badieirostami, M. A dielectrophoresis-based microfluidic chip for trapping circulating tumor cells using a porous membrane. J Micromech Microeng 2021, 32, 015008. [Google Scholar] [CrossRef]
  2. Lv, B.; Cai, J. Simulation and analysis of geometric parameters based on Taguchi method in YY microfluidic device for circulating tumor cell separation by alternating current dielectrophoresis. J Chromatogr A 2023, 1693, 463894. [Google Scholar] [CrossRef] [PubMed]
  3. Valijam, S.; Salehi, A.; Andersson, M. Design of a low-voltage dielectrophoresis lab-on-the chip to separate tumor and blood cells. Microfluid Nanofluidics 2023, 27, 22. [Google Scholar] [CrossRef]
  4. Varmazyari, V.; Habibiyan, H.; Ghafoorifard, H.; Ebrahimi, M.; Ghafouri-Fard, S. A dielectrophoresis-based microfluidic system having double-sided optimized 3D electrodes for label-free cancer cell separation with preserving cell viability. Sci Rep 2022, 12, 12100. [Google Scholar] [CrossRef]
  5. Lin, S.H.; Su, T.C.; Huang, S.J.; Jen, C.P. Enhancing the efficiency of lung cancer cell capture using microfluidic dielectrophoresis and aptamer-based surface modification. Electrophoresis 2024, 45, 1088–1098. [Google Scholar] [CrossRef]
  6. Sarno, B.; Heineck, D.; Heller, M. J.; Ibsen, S.D. Dielectrophoresis: Developments and applications from 2010 to 2020. Electrophoresis 2021, 42, 539–564. [Google Scholar] [CrossRef]
  7. Shokouhmand, H.; Abdollahi, A. Detection of cell-free DNA nanoparticles in insulator based dielectrophoresis systems. J Chromatogr A 2020, 1626, 461262. [Google Scholar] [CrossRef]
  8. Huang, Y.; Wang, X.B.; Becker, F.F.; Gascoyne, P.R.C. Introducing dielectrophoresis as a new force field for field-flow fractionation. Biophys J 1997, 73, 1118–1129. [Google Scholar] [CrossRef]
  9. Fiedler, S.; Shirley, S.G.; Schnelle, T.; Fuhr, G. Dielectrophoretic sorting of particles and cells in a microsystem. Anal Chem 1998, 70, 1909–1915. [Google Scholar] [CrossRef] [PubMed]
  10. Morgan, H.; Izquierdo, A.G.; Bakewell, D.; Green, N.G.; Ramos, A. The dielectrophoretic and travelling wave forces generated by interdigitated electrode arrays: Analytical solution using Fourier series. J Phys D Appl Phys 2001, 34, 1553. [Google Scholar] [CrossRef]
  11. Choi, S.; Park, J.K. Microfluidic system for dielectrophoretic separation based on a trapezoidal electrode array. Lab Chip 2005, 5, 1161–1167. [Google Scholar] [CrossRef] [PubMed]
  12. Imasato, H.; Yamakawa, T. Measurement of dielectrophoretic force by employing controllable gravitational force. J Electrophoresis 2008, 52, 1–8. [Google Scholar] [CrossRef]
  13. Lewpiriyawong, N.; Yang, C.; Lam, Y.C. Continuous sorting and separation of microparticles by size using AC dielectrophoresis in a PDMS microfluidic device with 3-D conducting PDMS composite electrodes. Electrophoresis 2010, 31, 2622–2631. [Google Scholar] [CrossRef]
  14. So, J.H.; Dickey, M.D. Inherently aligned microfluidic electrodes composed of liquid metal. Lab Chip 2011, 11, 905–911. [Google Scholar] [CrossRef]
  15. Fathy, J.; Pourmand, A.; Ghavifekr, H. B. Design and simulation of a MEMS based cell separator utilizing 3D travelling-wave dielectrophoresis. Microsyst Technol 2017, 23, 1351–1360. [Google Scholar] [CrossRef]
  16. Jiang, T.; Ren, Y.; Liu, W.; Tang, D.; Tao, Y.; Xue, R.; Jiang, H. Dielectrophoretic separation with a floating-electrode array embedded in microfabricated fluidic networks. Phys Fluids 2018, 30, 112003. [Google Scholar] [CrossRef]
  17. Beech, J.P.; Keim, K.; Ho, B.D.; Guiducci, C.; Tegenfeldt, J.O. Active posts in deterministic lateral displacement devices. Adv Mater Technol 2019, 4, 1900339. [Google Scholar] [CrossRef]
  18. Alnaimat, F.; Mathew, B.; Hilal-Alnaqbi, A. Modeling a dielectrophoretic microfluidic device with vertical interdigitated transducer electrodes for separation of microparticles based on size. Micromachines (Basel) 2020, 11, 563. [Google Scholar] [CrossRef]
  19. Huang, X.; Torres-Castro, K.; Varhue, W.; Salahi, A.; Rasin, A.; Honrado, C.; Brown, A.; Guler, J.; Swami, N. S. Self-aligned sequential lateral field non-uniformities over channel depth for high throughput dielectrophoretic cell deflection. Lab Chip 2021, 21, 835–843. [Google Scholar] [CrossRef]
  20. Tada, S.; Seki, Y. Analysis of temperature field in the dielectrophoresis-based microfluidic cell separation device. Fluids 2022, 7, 263. [Google Scholar] [CrossRef]
  21. Zhang, J.; Song, Z.; Liu, Q.; Song, Y. Recent advances in dielectrophoresis-based cell viability assessment. Electrophoresis 2020, 41, 917–932. [Google Scholar] [CrossRef]
  22. Nerguizian, V.; Alazzam, A.; Roman, D.; Stiharu, I.; Burnier Jr, M. Analytical solutions and validation of electric field and dielectrophoretic force in a bio-microfluidic channel. Electrophoresis 2012, 33, 426–435. [Google Scholar] [CrossRef] [PubMed]
  23. Alazzam, A.; Roman, D.; Nerguizian, V.; Stiharu, I.; Bhat, R. Analytical formulation of electric field and dielectrophoretic force for moving dielectrophoresis using Fourier series. Microfluid Nanofluid 2010, 9, 1115–1124. [Google Scholar] [CrossRef]
  24. Chen, R.; Liu, R.; Shen, H. Modeling and analysis of electric field and electrostatic adhesion force generated by interdigital electrodes for wall climbing robots. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Tokyo, Japan, 3.11.2013. [Google Scholar]
  25. Alhammadi, F.; Waheed, W.; El-Khasawneh, B.; Alazzam, A. Mathematical model and verification of electric field and dielectrophoresis in a microfluidic device. In Proceedings of the 2018 Advances in Science and Engineering Technology International Conferences (ASET), Dubai, UAE, 4.4.2020. [Google Scholar]
  26. Clague, D.; Wheeler, E. Dielectrophoretic manipulation of macromolecules: The electric field. Phys Rev E 2001, 64, 026605. [Google Scholar] [CrossRef] [PubMed]
  27. Gurtner, M.; Hengster-Movric, K.; Hurák, Z. Green’s function-based control-oriented modeling of electric field for dielectrophoresis. J Appl Phys 2017, 122, 054903. [Google Scholar] [CrossRef]
  28. Garcia, M.; Clague, D. The 2D electric field above a planar sequence of independent strip electrodes. J Phys D Appl Phys 2000, 33, 1747. [Google Scholar] [CrossRef]
  29. Wang, X.; Wang, X.B.; Becker, F.; Gascoyne, P.R.C. A theoretical method of electrical field analysis for dielectrophoretic electrode arrays using Green’s theorem. J Phys D Appl Phys 1996, 29, 1649–1660. [Google Scholar] [CrossRef]
  30. Bai, S.; Tang, Y.; Ruan, L.; Song, R.; Chen, H.; Du, Y.; Lin, H.; Tang, Y.; Shan, Y. Investigation into the influence of interdigital parameters on electrochemical performance for in-plane supercapacitors via mathematical modeling and conformal mapping techniques. J Energy Storage 2023, 65, 107287. [Google Scholar] [CrossRef]
  31. Blume, S. O.; Ben-Mrad, R.; Sullivan, P.E. Modelling the capacitance of multi-layer conductor-facing interdigitated electrode structures. Sens Actuators B Chem 2015, 213, 423–433. [Google Scholar] [CrossRef]
  32. Pampin, R.S.; Raskin, J.; Huynen, I.; Flandre, D. Electrodes-oxide-semiconductor device for biosensing: Renewed conformal analysis and multilayer algorithm. J Electroanal Chem 2020, 856, 113651. [Google Scholar] [CrossRef]
  33. Sun, T.; Green, N.G.; Gawad, S.; Morgan, H. Analytical electric field and sensitivity analysis for two microfluidic impedance cytometer designs. IET Nanobiotechnol 2007, 1, 69–79. [Google Scholar] [CrossRef] [PubMed]
  34. Sun, T.; Morgan, H.; Green, N.G. Analytical solutions of ac electrokinetics in interdigitated electrode arrays: Electric field, dielectrophoretic and traveling-wave dielectrophoretic forces. Phys Rev E 2007, 76, 046610. [Google Scholar] [CrossRef] [PubMed]
  35. Bao, X.; Crupi, G.; Ocket, I.; Bao, J.; Ceyssens, F.; Kraft, M.; Nauwelaers, B.; Schreurs, D. Numerical modeling of two microwave sensors for biomedical applications. Int J Numer Model El 2021, 34, e2810. [Google Scholar] [CrossRef]
  36. Ghomian, T.; Hihath, J. Review of dielectrophoretic manipulation of micro and nanomaterials: Fundamentals, recent developments, and challenges. IEEE Trans Biomed Eng 2022, 70, 27–41. [Google Scholar] [CrossRef]
  37. Miura, T.; Uno, S. Computer simulation for electrochemical impedance of a living cell adhered on the inter-digitated electrode sensors. Jpn J Appl Phys 2019, 58, SBBG15. [Google Scholar]
  38. Santos-Neto, I.S.d.; Carvalho, C.D.; Filho, G.B.A.; Andrade, C.D.S.S.; Santos, G.C. d.O.; Barros, A.K.; Neto, J.V.d.F.; Casas, V.L.P.; Alencar, L.M.R.; Lopes, A.J.O. Interdigitated electrode for electrical characterization of commercial pseudo-binary biodiesel–diesel blends. Sensors 2021, 21(21), 7288. [Google Scholar] [CrossRef]
  39. Zaman, M.A.; Padhy, P.; Ren, W.; Wu, M.; Hesselink, L. Microparticle transport along a planar electrode array using moving dielectrophoresis. J Appl Phys 2021, 130(3), 034902. [Google Scholar] [CrossRef]
  40. Tada, S.; Hayashi, M.; Eguchi, M.; Tsukamoto, A. High-throughput separation of cells by dielectrophoresis enhanced with 3D gradient AC electric field. Biomicrofluidics 2017, 11, 064110. [Google Scholar] [CrossRef]
  41. Duffy, D.G. Mixed Boundary Value Problems, 1st ed.; Chapman and Hall/CRC: New York, USA, 2008; pp. 1–14. [Google Scholar]
  42. Feng, J.J.; Krishnamoorthy, S.; Chen, Z.J.; Makhijani, V.B. Numerical and analytical studies of AC electric field in dielectrophoretic electrode arrays. In Proceedings of the Technical Proceedings of the 2002 International Conference on Modeling and Simulation of Microsystems, Puerto Rico, USA, 22.4.2002. [Google Scholar]
  43. Jones, T.B. Electromechanics of Particles, 2nd ed.; Cambridge University Press: Cambridge, UK, 1995; pp. 34–48. [Google Scholar]
  44. Henslee, E.A.; Sano, M.B.; Rojas, A.D.; Schmelz, E.M.; Davalos, R.V. Selective concentration of human cancer cells using contactless dielectrophoresis. Electrophoresis 2011, 32, 2523–2529. [Google Scholar] [CrossRef] [PubMed]
  45. Gauthier, V.; Bolopion, A.; Gauthier, M. Analytical formulation of the electric field induced by electrode arrays: Towards automated dielectrophoretic cell sorting. Micromachines (Basel) 2017, 8, 253. [Google Scholar] [CrossRef]
Preprints 120926 g001
Figure 1. Frequency spectrum of R e ( β ) for various σ l .
Figure 1. Frequency spectrum of R e ( β ) for various σ l .
Preprints 120926 g001
Preprints 120926 g002
Figure 2. (a) Schematic of cell-separation device. (b) 2D model of cell-separation device.
Figure 2. (a) Schematic of cell-separation device. (b) 2D model of cell-separation device.
Preprints 120926 g002
Preprints 120926 g003
Figure 3. Mathematical model and boundary conditions.
Figure 3. Mathematical model and boundary conditions.
Preprints 120926 g003
Preprints 120926 g004
Figure 4. Distribution of φ ( = V x * ) for various values of C .
Figure 4. Distribution of φ ( = V x * ) for various values of C .
Preprints 120926 g004
Preprints 120926 g005
Figure 5. (a) Specifications of electrode substrate. (b) Configuration and dimensions of interdigitated electrodes.
Figure 5. (a) Specifications of electrode substrate. (b) Configuration and dimensions of interdigitated electrodes.
Preprints 120926 g005
Preprints 120926 g006
Figure 6. Schematic of the experimental setup.
Figure 6. Schematic of the experimental setup.
Preprints 120926 g006
Preprints 120926 g007
Figure 7. Comparison of analytical and numerical results for distributions of (left) electric potential and (right) electric field intensity.
Figure 7. Comparison of analytical and numerical results for distributions of (left) electric potential and (right) electric field intensity.
Preprints 120926 g007
Preprints 120926 g008
Figure 8. Distribution of ∇ E * 2
Figure 8. Distribution of ∇ E * 2
Preprints 120926 g008
Preprints 120926 g009
Figure 9. Distribution of DEP forces generated by three types of interdigitated electrodes.
Figure 9. Distribution of DEP forces generated by three types of interdigitated electrodes.
Preprints 120926 g009
Preprints 120926 g010
Figure 10. 3D fluorescent imaging of live cell distribution.
Figure 10. 3D fluorescent imaging of live cell distribution.
Preprints 120926 g010
Preprints 120926 g011
Figure 11. Variation of DEP force with variations in electrode width and spacing.
Figure 11. Variation of DEP force with variations in electrode width and spacing.
Preprints 120926 g011
Table 1. Analysis conditions.
Table 1. Analysis conditions.
Symbol Values
RMS voltage V r m s 5 / 2 V
electric permittivity of the solution ϵ l 6.9 × 10 − 10 F/m
electric conductivity of the solution σ l 1.0 × 10 − 3 S/m
cell diameter d c 15 μ m
height of the microchannel H 500 μ m
real part of the CM factor R e ( β ) 0.98
number of terms of the series expansion n 100,000
electrode width and spacing 　( 2 w ：d) Type A (50 μm：50 μm)
Type B (50 μm：75 μm)
Type C (75 μm：50 μm)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated