# Theoretical and experimental study of the role of cell-cell dipole interaction in dielectrophoretic devices: application to polynomial electrodes

- Massimo Camarda
^{1}Email author, - Giuseppe Fisicaro
^{1}, - Ruggero Anzalone
^{1}, - Silvia Scalese
^{1}, - Alessandra Alberti
^{1}, - Francesco La Via
^{1}, - Antonino La Magna
^{1}, - Andrea Ballo
^{2}, - Gianluca Giustolisi
^{2}, - Luigi Minafra
^{3}, - Francesco P Cammarata
^{3}, - Valentina Bravatà
^{3}, - Giusi I Forte
^{3}, - Giorgio Russo
^{3}and - Maria C Gilardi
^{3}

**13**:71

https://doi.org/10.1186/1475-925X-13-71

© Camarda et al.; licensee BioMed Central Ltd. 2014

**Received: **17 November 2013

**Accepted: **23 May 2014

**Published: **5 June 2014

## Abstract

### Background

We aimed to investigate the effect of cell-cell dipole interactions in the equilibrium distributions in dielectrophoretic devices.

### Methods

We used a three dimensional coupled Monte Carlo-Poisson method to theoretically study the final distribution of a system of uncharged polarizable particles suspended in a static liquid medium under the action of an oscillating non-uniform electric field generated by polynomial electrodes. The simulated distributions have been compared with experimental ones observed in the case of MDA-MB-231 cells in the same operating conditions.

### Results

The real and simulated distributions are consistent. In both cases the cells distribution near the electrodes is dominated by cell-cell dipole interactions which generate long chains.

### Conclusions

The agreement between real and simulated cells’ distributions demonstrate the method’s reliability. The distribution are dominated by cell-cell dipole interactions even at low density regimes (10^{5} cell/ml). An improved estimate for the density threshold governing the interaction free regime is suggested.

## Background

As first defined by Pohl [1, 2], the term “*dielectrophoresis*” is used to describe the *“ponderomotive”* force exerted by a non-uniform electric field on polarizable neutral particles. Such force allows for the controlled manipulation of micro and nano-sized particles dispersed in colloidal solutions. Application fields include: cell partitioning and isolation [3, 4], bio-structure assembling [5], nanostructure (e.g. carbon nanotube) deposition [6], filtration systems for oils purification [7] etc. Among these, the separation of rare cells [8] is a specifically promising one as *dielectrophoresis* allows the capture/separations of cells without the use of biomarkers; relying, instead, in the strong selectivity of the dielectrophoretic (DEP) response [9] which depends on the particle mass, shape and composition. Indeed, recently, this selectivity has permitted to discriminate the tumor cell types of the NCI-60 panel from Peripheral Blood MonoNuclear cells (PBMNs) [10]. However, although many intriguing micro-structures have been fabricated in research laboratories, DEP devices have hardly gone beyond the proof-of-concept stage [11]. One of the problems that are hindering development and engineerization of the devices is the limited use of accurate numerical tools for their design which, in turn, is due to the computational complications arising by the particle-particle dipole interaction. Indeed, particle kinetics (i.e. the particle velocity field) in DEP devices can be easily calculated by mean of Poisson solvers and direct integration of the equation of motion only in the non-interacting particle approximation [12]. This approximation is not valid in the accumulation regions of the DEP devices where, due to the increased particle concentration, dipole-dipole interactions become important and can promote the formation of clusters and significant rearrangements of the particle space distribution [13–16]. These many-particle effects can be accurately simulated solving directly the equations of motion in the few-particles limit [17] i.e. this approach is not applicable for the simulation and design of realistic systems. Another possible approach is the use of reaction–diffusion models [18–20] but this approach needs an “ad hoc” parameter calibration to effectively consider the dipole-dipole interactions in compact models.

Recently a coupled Monte Carlo-Poisson (MC-P) method [21] has been implemented which allows simulating a large number of particles in large active zones (within the experimental range), explicitly including particle-particle interactions. The MC-P method has pointed out the relevance of this inclusion in the modeling predictions for the simplified condition of Two Dimensional (2D) electric field $\overrightarrow{\mathit{E}}\left(\overrightarrow{\mathit{r}}\right)$ distribution, where $\overrightarrow{\mathit{E}}\left(\overrightarrow{\mathit{r}}\right)$ explicitly depends on two space coordinates $\overrightarrow{\mathit{E}}\left(\overrightarrow{\mathit{r}}\right)\equiv \overrightarrow{\mathit{E}}\left(\mathit{x},\mathit{y}\right)$ as in the case of very long interdigitated electrodes. However, the possibility to apply this approach for the numerical design of devices exploiting more complex fully Three Dimensional (3D) electric field distributions has not been yet demonstrated. Moreover the MC-P predictions have never been compared with real cell distributions in dielectrophoretic devices, in order to confirm their reliability. Aiming to the two objectives of the model extension and validation, we have improved the application potentiality of the MC-P method to simulate the features of devices generating 3D electric field distributions. In addition we applied the simulation method to the case of polynomial electrodes which are known to produce well defined 3D non-uniform electric fields and are used for the study of negative dielectrophoresis [22] or for the determination of particle dielectrophoretic response through electrorotation analysis [23, 24]. We compare the simulated results with experimental distributions obtained in the same electrodes geometry to evaluate the role of p-p interactions and definitively demonstrate the predictive potential of this methodology.

## Methods

### Computational algorithm

A detailed description of the method can be found in ref. [21], here we summarize the key aspects of the simulations, specifically focusing on the 3D implementation.

*f*=

*ω*/2

*π*frequency, will be subjected to an effective averaged potential energy [25]:

*α*

_{ eff }is the average polarizability of the particle defined as:

*V*is the particle volume, ${\tilde{\mathit{\u03f5}}}_{1,\mathit{m}}={\mathit{\u03f5}}_{1,\mathit{m}}-\mathit{i}{\mathit{\sigma}}_{1,\mathit{m}}/\mathit{\omega}$ are the complex dielectric constants of the particle (1) and the media (m) and

*f*

_{ CM }(

*ω*) is the Clausius-Mossotti factor which fully characterizes the dielectric response of the particle in the given medium. The isolated particle approximation holds only in the diluted density limit, i.e. only if the average distance between two particles in the colloidal solution is always very large otherwise an effective particle-particle interaction has to be considered as a result of the local distortion of the electric field lines generated by itself and by the other particles. In this case, the DEP force acting on each suspended particle can be directly calculated by means of the Maxwell tensor [9]: $\overleftrightarrow{\mathit{T}}\equiv {\mathit{T}}_{\mathit{ij}}={\tilde{\mathit{\u03f5}}}_{\mathit{m}}\left({\mathit{E}}_{\mathit{i}}^{\mathit{tot}}{\mathit{E}}_{\mathit{j}}^{\mathit{tot}}-0.5{\mathit{\delta}}_{\mathit{ij}}{\mathit{E}}_{\mathit{k}}^{\mathit{tot}}{\mathit{E}}_{\mathit{k}}^{\mathit{tot}}\right)$ over the closing surface of the particles:

Note that whereas the $\overrightarrow{\mathit{E}}$ field in Eq. 1 is the field generated by the external electrodes only, ${\overrightarrow{\mathit{E}}}^{\mathit{tot}}$ in the Eq. 3 must be calculated *considering all the particle presence*. This direct calculation is not practically feasible in the kinetic simulation of large systems, since the particle distribution continuously changes in the space requiring an integration of the Poisson equation, $\nabla {\overrightarrow{\mathit{E}}}^{\mathit{tot}}=0$, at each simulation step. A more efficient approach, that requires only the evaluation of the external field, can be implemented approximating the total distorted electric field with the sum of the field generated by the external electrodes plus the contributions of the dipoles induced in all the particles: ${\mathbf{E}}^{\mathit{tot}}\left({\overrightarrow{\mathit{r}}}_{\mathit{i}}\right)\approx \mathbf{E}\left({\overrightarrow{\mathit{r}}}_{\mathit{i}}\right)+{\displaystyle \sum _{\mathit{j}}^{\mathit{all\; the\; particles}}{\mathbf{E}}_{\mathit{j}}\left({\overrightarrow{\mathit{r}}}_{\mathit{i}}\right)}$
[26].

The reliability of this approximation has been demonstrated in Ref. [13] with the aid of the full calculation based on Eq. 3 for the case of two spherical particles immersed in a uniform external electric field: the magnitude, the angular dependency and the scaling with the distance of the calculated force are similar to those derived in the interacting dipoles approximation [26].

*r*

_{ j }is:

where ${\mathit{\theta}}_{\mathit{ij}}^{\mathit{i}},\phantom{\rule{0.5em}{0ex}}{\mathit{\theta}}_{\mathit{ij}}^{\mathit{j}}$ are the angles between the vectors $\overrightarrow{\mathit{E}}\left({\overrightarrow{\mathit{r}}}_{\mathit{i}}\right),\phantom{\rule{0.5em}{0ex}}\overrightarrow{\mathit{E}}\left({\overrightarrow{\mathit{r}}}_{\mathit{j}}\right)$, and $\stackrel{\rightharpoonup}{\mathit{n}}$ and ${\mathit{\alpha}}_{\mathit{eff}}^{\mathit{i}},\phantom{\rule{0.5em}{0ex}}{\mathit{\alpha}}_{\mathit{eff}}^{\mathit{j}}$ are the average polarizations for the *i* and *j* particles.

*r*

_{ i }and the configuration energy is

where ${\overline{\mathit{U}}}_{\mathit{eff}}$ and ${\overline{\mathit{U}}}_{\mathit{ij}}$ are calculated by means of the equations (1) and (6).

*Δt*between two consecutive displacement events as

where $\overline{\mathit{D}}\approx {\mathit{D}}_{\mathit{Brow}}$ is the effective diffusivity, *N* is the number of simulated particles, Δ*d* is the elementary particle displacement, *η* is the medium viscosity, *a* is the particle radius, *k*
_{
B
} is the Boltzmann constant and *T* is the system temperature.

### Experimental setup

The polynomial electrode design described in the previous section, has been fabricated by deposition of 10 nm of Titanium followed by 200 nm of Nickel on a standard microscope glass. The electrodes were delineated by lithographic methods followed by wet etching. The device has been energized using a Protek 9205C signal generator which applied, consistently with the simulated systems, a sinusoidal voltage signal of 8V_{pp} value at 1 MHz for 180 sec (long time allow for cells equilibration). The final distribution was observed with a standard 10× phase contrast inverted microscope. The human breast cancer cell line MDA-MB-231 were cultured according to American Type Culture Collection (ATCC) instructions. The cells, just before DEP tests, were suspended in a low conductive buffer (used as the elute in all our experiments) composed of 9.5% ultrapure sucrose (S7903, Sigma-Aldrich), 0.3% dextrose (Fisher D-16), and 0.1% Pluronic F68 (P1300, Sigma-Aldrich) titrated to a conductivity of 30 mS/m (consistent with Monte Carlo simulations) by adding KCl with the aid of a conductivity meter. The buffer had an osmolarity of 320 mOs/L and a pH of 7 and the experiments were conducted at room temperature (22°C). The cells, suspended in DEP buffer at concentration of 5×10^{5} cells/ml were pipetted into the chamber and occupied a total volume of about 100 μl when a cover slip was placed over the rubber o-ring.

## Results

*a = 6.2*μ

*m*and

*d = 10 nm*are the cell radius and membrane thickness

**.**The following values have been used: ${\mathit{\u03f5}}_{\mathit{me}}^{\text{'}}=24{\mathit{\u03f5}}_{0},\phantom{\rule{0.5em}{0ex}}{\mathit{\u03f5}}_{\mathit{in}}^{\text{'}}=50{\mathit{\u03f5}}_{0},\phantom{\rule{0.5em}{0ex}}{\mathit{\sigma}}_{\mathit{m}}=30\phantom{\rule{0.5em}{0ex}}\mathit{mS}/\mathit{m},\phantom{\rule{0.5em}{0ex}}{\mathit{\sigma}}_{\mathit{in}}=0.2\phantom{\rule{0.5em}{0ex}}\mathit{S}/\mathit{m},\phantom{\rule{0.5em}{0ex}}{\mathit{\sigma}}_{\mathit{me}}={10}^{-7}\mathit{S}/\mathit{m}$, [28], and the calculated real part of the factor

*f*

_{ CM }(

*ω*) is shown in Figure 1.

^{3}with four polynomial electrodes located at the bottom of the box (see Figure 2) whose shapes can be described by the following parametric system:

*D*represents half the distance of opposing electrodes whereas

*L*is related to the electrode width. Referring to Figure 2,

*right*and using Eq. 10, we have that ${\mathit{L}}_{1}=\sqrt{2}\mathit{D}$ and ${\mathit{L}}_{2}=\sqrt{2}\left(\mathit{L}-\sqrt{{\mathit{L}}^{2}-{\mathit{D}}^{2}}\right)$. To improve the capturing efficiency of the DEP system the electric fields in

*p*

_{ 1 }and

*p*

_{ 2 }must be of the same order (otherwise the trapping regions will be limited to the

*p*

_{ 2 }regions). Given that |

*E*

_{ rms }|(

*p*

_{1}) ∝

*V*

_{ pp }/

*L*

_{1}and |

*E*

_{ rms }|(

*p*

_{2}) ∝

*V*

_{ pp }/

*L*

_{2}

*D*must be of the same order of

*L*, in this study we sat

*D*= 390 μm and

*L*= 460 μm (note that in ref. [22]

*D =*64 μm and

*L*= 90 μm which favored negative-DEP only). Figure 2,

*left*shows the simulated box with the four electrodes situated at the bottom. The intensity map of $\left|{\mathit{E}}_{\mathit{rms}}\right|\left(\overrightarrow{\mathit{r}}\right)$ is also reported in grey scale for the oscillating four-electrodes configuration at a frequency of 1 MHz and a peak value of 8 V for

*V*

_{ pp }with 180° phase difference between neighbor electrodes (i.e.

*V*

_{1}(

*t*) = −

*V*

_{2}(

*t*) =

*V*

_{ pp }

*sin*(2

*πft*)). Figure 2,

*right*shows the electric field at 50 μm from the bottom surface, associated to the considered system. The regions of high electric field (to which the cells will move under positive-DEP) are located all around the edges of the electrodes. In the considered system the highest electrical fields are located at

*p*

_{ 1 }and

*p*

_{ 2 }equivalent points.

^{5}cells/ml) (

*left*) and the final (

*right*) equilibrium condition after

*2*×

*10*

^{ 8 }Monte Carlo iterations. The final distribution is the result of the minimization of Eq. 7, which induce a movement towards the

*p*

_{ 1 }and

*p*

_{ 2 }regions (minimization of the first term) and an aligned of the cells (causing cell-chains) along the electric field direction (minimization of the second term). Specifically particles tend to form larger chains in the regions of more intense electric field which are located mainly around

*p*

_{ 2 }equivalent points, similar distributions have been observed for densities as low as 1×10

^{5}cells/ml (not shown). It is important to note that this tendency does not depends on the DEP polarity because the ${\mathit{\alpha}}_{\mathit{eff}}^{\mathit{i}}{\mathit{\alpha}}_{\mathit{eff}}^{\mathit{j}}$ term in Eq. 7 is always positive for identical cells (whereas the ${\mathit{\alpha}}_{\mathit{eff}}^{\mathit{i}}$ term in Eq. 1 can be positive or negative, depending on the value of

*f*

_{ CM }(

*ω*)).

From these results we can infer that particle-particle interactions compete with the dielectrophoretic force-field, which would otherwise massively trap (in p-DEP conditions) the particles in the regions where the gradient of the electric field is larger. Note also the cells in the region far away from the electrodes which are not trapped by the DEP field, this allows for the definition of a depletion volume as the region where cells are effectively attracted to the electrodes. The connection between depletion region, electrodes geometry and particle-particle interaction is currently under investigation.

*d*falls below ~5. In a simple approximation we could assume that all cells will eventually settle to the bottom of the chamber, which gives:

*h*is the chamber height and

*C*is the cell concentration. Setting

*h*= 1500 μm and

*d*= 2 ×

*r*= 12.4

*μm*we obtain

*C*< 5 × 10

^{5}cells/ml. According to this rough estimate, the density used in the Monte Carlo simulation (

*C*= 10

^{5}cells/ml) should not lead to significant dipole-dipole interactions. The limit of this simple description is that it implicitly assumes a uniform distribution of cells at the bottom of the chamber, i.e. it neglects the fact that the cells, due to the DEP field, will concentrate at the electrodes edges (see Figure 4). In order to qualitatively take it into account this effect, it is possible to correct the formula using the following parameter:

*E*

^{2}|

_{ Max }and 〈∇

*E*

^{2}〉 are, the highest and average value of the electric field gradient calculated at the bottom of the chamber. Since the dielectrophoretic force is proportional to the electric field gradient,

*r*

_{ DEP }approximately represents the concentration factor generated by the DEP force. So that we improve the previous Eq. 11 as:

*all*the cells settle at the bottom of the chamber. This is not generally true (see Figure 3) and it depends on the allowed deposition time and on the electric fields. To take this into account we need to substitute, in Eq. 13 the chamber height

*h*with the effective capturing region (

*z*

_{ capture }):

Where 〈∇_{
z
}
*U*
_{
eff
}(*z*)〉 is the average DEP force in the vertical direction at distance z from the bottom of the chamber and μ_{
m
} is medium viscosity. In the case of the polynomial electrodes used and for a deposition time of 180 sec, *r*
_{
DEP
} ≅ 85 and *z*
_{
cap
} = 280*μm* ≅ (1/5)*h* so that the improved concentration threshold, to avoid dipole-dipole interaction, should be below 3 × 10^{4} cells/ml. We performed Monte Carlo simulations in the 10^{4} cells/ml range finding no significant evidence of cell chains formation (not shown), thus confirming that Eq. 14, together with Eq. 12, represent a better qualitative threshold to avoid cell-cell interaction requiring only a knowledge of the electric field in the DEP device.

## Conclusions

In conclusion, we have demonstrated that the effects of particle-particle interactions play a crucial role in the kinetic evolution of colloidal systems in DEP devices even at low density regimes (10^{5} cells/ml), lower than the ones currently used in DEP devices [30].

Monte Carlo methods allow for the simulation of sufficiently large systems in terms of size and number of particles (i.e. within the experimental scopes). The discrete approach (i.e. the particle resolution), as opposed to the fluid-flow methodologies, is the key ingredient of the method improvement. In the case of MDA-MB-231 tumor cells suspended in a static, low conductive, fluid under the action of a positive-DEP field generated by a polynomial schema we have elucidated the crucial role of particle-particle interactions on the trapping efficiency of the device, on the organization of cells in ordered chains and on the overall cell space distribution. We have also deduced a new qualitative concentration threshold to avoid cell-cell interaction which requires only a knowledge of the electric field in the DEP device. Clearly, to have an exact determination of the concentration threshold for the specific DEP device used, MC-P kinetic simulations varying the cells density, such as the one proposed in this paper, need to be performed.

Future works will be devoted to generalize the formalism here presented in order to include second-order effects such as cell-sedimentation and cell-stitching or including hydrodynamic forces to simulate cells distributions in dynamic separation systems.

## Declarations

### Acknowledgements

We thank Salvo Di Franco of CNR-IMM for the lithographic fabrication steps and Caterina Grillo for the English Grammar and Spelling review.

## Authors’ Affiliations

## References

- Pohl HA, Pethig R:
**Dielectric measurements using non-uniform electric field (dielectrophoretic) effects.***J Phys E Sci Instrum*1977,**10:**190. 10.1088/0022-3735/10/2/022View ArticleGoogle Scholar - Pohl H:
**The motion and precipitation of suspensoids in divergent electric fields.***J Appl Physics*1951,**22:**869. 10.1063/1.1700065View ArticleGoogle Scholar - Gascoyne PRC, Noshari J, Anderson TJ, Becker FF:
**Isolation of rare cells from cell mixtures by dielectrophoresis.***Electrophoresis*2009,**30:**1388. 10.1002/elps.200800373View ArticleGoogle Scholar - Huges MP:
**Strategies for dielectrophoretic separation in laboratory-on-a-chip systems.***Electrophoresis*2002,**23:**2569. 10.1002/1522-2683(200208)23:16<2569::AID-ELPS2569>3.0.CO;2-MView ArticleGoogle Scholar - Castillo J, Tanzi S, Dimaki M, Svendsen W:
**Manipulation of self-assembly amyloid peptide nanotubes by dielectrophoresis.***Electrophoresis*2008,**29:**5026. 10.1002/elps.200800260View ArticleGoogle Scholar - Duchamp M, Lee K, Dwir B, Seo JW, Kapon E, Forro L, Magrez A:
**Controlled positioning of carbon nanotubes by dielectrophoresis: insights into the solvent and substrate role.***ACS Nano*2010,**4:**279. 10.1021/nn901559qView ArticleGoogle Scholar - Shen Y, Elele E, Khusid B:
**A novel concept of dielectrophoretic engine oil filter.***Electrophoresis*2011,**32:**2559. 10.1002/elps.201100072View ArticleGoogle Scholar - Shin S, Stemke-Hale K, Tsimberidou AM, Noshari J, Anderson TE:
**Antibody-independent isolation of circulating tumor cells by continuous-flow dielectrophoresis.***Biomicrofluidics*2013,**7:**011807. 10.1063/1.4774304View ArticleGoogle Scholar - Cetin B, Li D:
**Dielectrophoresis in microfluidics technology.***Electrophoresis*2011,**32:**2410. 10.1002/elps.201100167View ArticleGoogle Scholar - Shim S, Stemke-Hale K, Noshari J, Becker F, Gascoyne PRC:
**Dielectrophoresis has broad applicability to marker-free isolation of tumor cells from blood by microfluidic systems.***Biomicrofluidics*2013,**7:**011808. 10.1063/1.4774307View ArticleGoogle Scholar - Nathan B:
**Microfluidics: the great divide.***Nat Methods*2009,**6:**683. 10.1038/nmeth0909-683View ArticleGoogle Scholar - Morgan H, Hughes MP, Green NG:
**Separation of submicron bioparticles by dielectrophoresis.***Biophys J*1999,**77:**516. 10.1016/S0006-3495(99)76908-0View ArticleGoogle Scholar - KH K, Li D:
**Dielectric force and relative motion between Two spherical particles in electrophoresis.***Langmuir*2006,**22:**1602. 10.1021/la052162kView ArticleGoogle Scholar - Lee D, Yu C, Papazoglou E, Farouk B, Noh HM:
**Dielectrophoretic particle–particle interaction under AC electrohydrodynamic flow conditions.***Electrophoresis*2011,**32:**2298.View ArticleGoogle Scholar - Lewpiriyawong N, Yang C:
**Dielectrophoresis field-flow fractionation for continuous-flow separation of particles and cells in microfluidic devices.***Advance in Transport Phenomena*2014,**3:**29. 10.1007/978-3-319-01793-8_2View ArticleGoogle Scholar - Mohtar MN, Hoettges KF, Hughes MP:
**Factors affecting particle collection by electro-osmosis in microfluidic systems.***Electrophoresis*2014,**35:**345. 10.1002/elps.201300420View ArticleGoogle Scholar - Arun TJ, Singh KP, Aubry N:
**Dielectrophoresis of nanoparticles.***Electrophoresis*2004,**25:**3625. 10.1002/elps.200406092View ArticleGoogle Scholar - Nicotra OE, La Magna A, Coffa S:
**A mean field approach to many-particles effects in dielectrophoresis.***Appl Phys Lett*2008,**93:**193902. 10.1063/1.3025840View ArticleGoogle Scholar - Nicotra OE, La Magna A, Coffa S:
**Particle-chain formation in a dc dielectrophoretic trap; a reaction–diffusion approach.***Appl Phys Lett*2009,**95:**073702. 10.1063/1.3177332View ArticleGoogle Scholar - Nicotra OE, La Magna A, Coffa S:
**Numerical study of particle-particle interactions in Lab-on-a-chip devices based on dielectrophoresis.***Sens Lett*2010,**8:**457. 10.1166/sl.2010.1294View ArticleGoogle Scholar - La Magna A, Camarda M, Deretzis I, Fisicaro G, Coffa S:
**Coupled Monte Carlo-Poisson method for the simulation of particle-particle effects in dielectrophoretic devices.***Appl Phys Lett*2012,**100:**134104. 10.1063/1.3697733View ArticleGoogle Scholar - Huang Y, Pethig R:
**Electrode design for negative dielectrophoresis.***Meas Sci Technol*1991,**2:**1142. 10.1088/0957-0233/2/12/005View ArticleGoogle Scholar - Gimsa J, Marszalek P, Loewe U, Tsong T:
**Dielectrophoresis and electrorotation of neurospora slime and murine myeloma cells.***J Biophys*1991,**60:**749. 10.1016/S0006-3495(91)82109-9View ArticleGoogle Scholar - Huang Y, Ralph H, Pethg R, Wang X:
**Differences in the AC electrodynamics of viable and non-viable yeast cells determined through combined dielectrophoresis and electrorotation studies.***Phys Med Biol*1992,**37:**1499. 10.1088/0031-9155/37/7/003View ArticleGoogle Scholar - Sancho M, Giner V, Martinez G:
**Monte Carlo simulation of dielectrophoretic particle chain formation.***Phys Rev E*1997,**55:**544.View ArticleGoogle Scholar - Simonova TS, Shilov VN, Shramko OA:
**Low-frequency dielectrophoresis and the polarization interaction of uncharged spherical particles with an induced Debye atmosphere of arbitrary thickness.***Colloid Jour*2001,**63:**108. 10.1023/A:1009458728584View ArticleGoogle Scholar -
**COMSOL 4.3b**http://www.comsol.com - Gascoyne PRC, Shim S, Noshari J, Becker FF, Stemke-Hale K:
**Correlations between the dielectric properties and exterior morphology of cells revealed by dielectrophoretic field-flow fractionation.***Electrophoresis*2013,**34:**1042. 10.1002/elps.201200496View ArticleGoogle Scholar - Jones TB:
*Electromechanics of Particles*. New York: Cambridge University Press; 1995. Chapters 6 and 7, pp.139–226View ArticleGoogle Scholar - Gupta V, Jafferji I, Garza M, Melnikova VO, Hasegawa DK, Pethig R, Davis DW:
**ApoStream™, a new dielectrophoretic device for antibody independent isolation and recovery of viable cancer cells from blood.***Biomicrofluidics*2012,**6:**024133. 10.1063/1.4731647View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.