2017, Vol. 38
介电泳如今已经成为微/纳米粒子控制的热门技术, 其广泛应用在细胞分析[1]、癌症检测[2]和药品研发[3]等诸多医药领域.介电泳力是一种电极化净力, 当电中性粒子处在非均匀电场中, 由于粒子两端的电极化力无法相互抵消而发生定向运动, 此力的方向与电场强度的梯度方向一致.因此, 非均匀电场是产生介电泳力实现粒子控制的前提条件.此外, 介电泳力的大小还直接与受控粒子的形状相关.介电泳力计算通常采用经典Pohl表达式[4-5], 该公式方法简单, 但是仅考虑二阶极子, 忽略了高阶极子的影响.针对分子动力学模拟非球形细胞承受介电泳力运动, 仿真结果难免与实验观测不符.导致误差的主要原因是高阶极子形成的高阶介电泳力会改变细胞运动方向.深入分析非球形细胞的高阶介电泳力能为实验上细胞的精确控制提供扎实的理论基础.
1 理论模型有效极矩法[6]将细胞在非均匀电场受到的介电泳力看作偶极子、四极子、八极子等高阶极子电场力的叠加.因此介电泳力FDEP如式(1)所示:
|
(1) |
式中:E0为电场强度；p1和p2分别是偶极子和四极子.经典Pohl介电泳力忽略了四极子、八极子、十六极子等高阶项的影响.然而, 当细胞为非规则球体, 或其尺寸与金属电极比较接近, 细胞高阶极子的电极化强度将不可忽视.因此, Jones根据电介质粒子的法向电场具有连续性, 且切向电场强度为零, 依靠泰勒级数推导出了任意形状细胞的高阶介电泳力表达式.在此, N阶极子主要由N-1阶极子组合而成, 如图 1所示.
|
图 1 单极子组成的高阶偶极子示意图 Fig.1 Generation of general dipole, quardrupole, and octupole starting from monopolar charges (a)—点电荷；(b)—偶极子；(c)—四极子；(d)—八极子. |
两个电荷相反的点电荷保持一定空间距离构成了偶极子, 然后两个偶极子又组合成四阶偶极子, 依次递增可形成高阶极子.根据式(1)可知, 精确介电泳力的计算与高阶极子的求解密不可分.偶极子或高阶偶极子的求解可通过半径Rt的球面作闭合曲面积分, 并乘以各阶数的权重系数得到.如图 2所示, 以球坐标系(R, φ, θ)为例, 各阶数偶极子的极化强度可根据式(2)得到.
|
图 2 球坐标系下细胞曲线积分 Fig.2 Cellular line integral based on spherical coordinate system |
|
(2) |
式中:εm为水溶液的介电常数; n为极矩阶数; Ln是以cosθ为自变量的Legendre多项式(L1=cosθ; L2= 
有关pn的积分运算先要求解φRt, 它代表Rt半径球表面的电势大小, 此电势是细胞受电场极化所产生的感应电势.因此它包含了各阶数电极化矢量的叠加, 如式(3)所示.
|
(3) |
在此基础上, 本文研究的细胞形状包括三种:规则球体、椭球和红细胞, 红细胞形状呈双面微凹之圆盘状, 三维形状如图 3所示.
|
图 3 细胞的三维形状示意图 Fig.3 Schematic of 3D shape of cells |
求解极化电势φRt, 本文采用有限元COMSOL软件对电场进行分析, 电场满足Laplace偏微分方程[7].为了简化计算量, 并且涉及的细胞形状具有中心对称性, 所以将细胞三维模型降至二维模型进行研究.首先假设细胞不在电场求解域中, εm溶液下可得到空间均匀电势φfree.然后介电常数为εp的细胞放入溶液中, 又可得到存在细胞条件下的电势分布φexist.在此, 细胞周围的极化电势可通过式(4)得到.
|
(4) |
细胞的电极化矢量pn可分解为三个方向(x, y, z), 本文主要研究细胞z方向的电极化强度, x和y方向可依此类推.
因为细胞与周围溶液的电属性不同, 极易导致垂直电场Ez发生扭曲.此时电场的变化也可看作是细胞发生电极化与外界电场相互作用的结果.对于准静态时谐电场的求解, 本文外界两端激励电压初始设为±2V.图 4是采用COMSOL软件得到三种形状细胞的极化电势分布图.细胞内部的电势基本保持与x轴平行, 电场E=-▽φRt.能够清晰发现三个细胞的电极化矢量pn与z轴平行.
|
图 4 基于COMSOL软件对三种形状细胞的电势仿真 Fig.4 Numerically simulated electrical potential for three type of cells based on COMSOL software (a)—球形细胞的极化电势分布；(b)—椭球细胞的极化电势分布; (c)—红细胞的极化电势分布. |
首先研究规则球体的高阶介电泳力, 如Pohl经典公式所述, 介电泳力等于电偶极子p1点乘▽E0, 电偶极子的大小如式(5)所示.
|
(5) |
式中:a是球体半径; fCM即俗称的Clausius-Mossotti因子, 其值为复数, 解析式如式(6)所示.
|
(6) |
式中, 
|
图 5 球体各阶数极子的系数曲线图 Fig.5 Factor of each order moment for sphere |
研究的球体半径a=1, Rt=2, 那么根据式(2)~(4)可求出各阶数极子的系数值, 图 5中以离散点表示.需要指出, 介电泳力主要作用于电中性粒子, 因此它的高阶作用力与奇数项有关, 对应偶数项为零, 如式(7)所示.
|
(7) |
观察图 5可发现, 有效极矩法得到的偶极子与CM因子误差很小.同时高阶极子系数为0.这说明对于规则球体的介电泳力而言, 它们仅承受偶极子作用力.此结果也反映出Pohl表达式求解规则球体的介电泳力具有很高的精度[8].
然而, 现实生物细胞中, 很少细胞保持完全规则的球体形结构.接下来研究椭球细胞的高阶极子影响.以大肠杆菌为例, 它具有细长的杆形结构非常类似于椭球体.与球体不同, 椭球体的偶极子如下:
|
(8) |
ax, ay和az是椭球在x, y, z方向的轴长.Ai是x, y, z各方向的去极化因子:
|
(9) |
本文定义偏心率e=az/ax, 首先研究e>1条件的高阶极子曲线, e取值分别是10:1, 5:1, 2:1, 结果如图 6所示.模型保持az=1不变, 对应变化ax的长度.闭合曲线积分半径Rt仍取2.分析发现在较低εp时候, 椭球的z方向极化系数非常小.但是随着εp逐渐增大(或者频率增加), 高阶极子系数逐渐增大并趋于稳定.与球体不同, 此时的高阶极子系数非零, 说明椭球存在高阶介电泳力作用.然而随着阶数的增大, 高阶电极化的影响也越来越弱.当椭球az和ax尺寸近似时, 在较低εp数值(或低频条件下), 各阶极化系数为负数, 如图 6c所示.表明介电泳力与电场方向相反, 会产生负介电泳力, 主要原因是细胞的介电常数εp低于溶液介电常数εm所致.同时, 模型研究e < 1条件下的各阶极子系数, 如图 7所示.与偏心率大于1不同, p1, p5, p9的趋势和p3, p7变化趋势正好截然相反.从椭球e的变化可发现, 形态特征对高阶介电泳力的大小和方向影响很大.
|
图 6 偏心率大于1的椭球各阶数极子系数曲线图 Fig.6 Factor of each order moment for ellipsoids with different eccentricities (e>1) (a)—偏心率10:1椭球；(b)—偏心率5:1椭球；(c)—偏心率2:1椭球. |
|
图 7 偏心率小于1的椭球各阶数极子系数曲线图 Fig.7 Factor of each order moment for ellipsoids with different eccentricities (e < 1) (a)—偏心率1:10椭球；(b)—偏心率1:5椭球；(c)—偏心率1:2椭球. |
最后研究红细胞各阶极子, 红细胞中间凹陷的形状参考了Ponder的实测离散点[9], 通过5次多项式拟合得到.红细胞的尺寸进行了等比例缩放.设定其正向x轴的最大长度为1, 如图 4c所示.积分半径Rt仍为2.红细胞的奇数阶电偶极子系数如图 8所示.本文提供偏心率为1:5的椭球与之进行对比.不难发现红细胞各阶电偶极子系数与偏心率为1:5的椭球趋势相同.该模拟结果与文献[6]的红细胞研究结果完全相符.
|
图 8 红细胞与偏心率1:5椭球的奇数阶偶极子曲线图 Fig.8 Factor of odd moment for erythrocyte and ellipsoid with eccentricity ratio 1:5 |
1) 基于有效极矩法研究了三种不同形状细胞的高阶极子, 并得到对应高阶极子的变化曲线, 通过高阶电极化力的叠加, 可精确计算出细胞在非均匀电场下的介电泳力大小.
2) 高阶电偶极子的分析证实了规则球体满足Pohl经典介电泳力表达式, 可忽略四极子、八极子、十六极子等高阶项的影响, 同时中间凹陷两端凸起的红细胞介电泳力近似于偏心率为1:5的椭球.
| [1] | Sabuncu A C, Beskok A. A separability parameter for dielectrophoretic cell separation[J]. Electrophoresis, 2013, 34(7): 1051–1058. DOI:10.1002/elps.v34.7 |
| [2] | Mulhall H J, Labeed F H, Kazmi B, et al. Cancer, pre-cancer and normal oral cells distinguished by dielectrophoresis[J]. Analytical and Bioanalytical Chemistry, 2011, 401(8): 2455–2463. DOI:10.1007/s00216-011-5337-0 |
| [3] | Shields C W, Reyes C D, Lopez G P. Microfluidic cell sorting:a review of the advances in the separation of cells from debulking to rare cell isolation[J]. Lab on a Chip, 2015, 15(5): 1230–1249. DOI:10.1039/C4LC01246A |
| [4] | Chen D F, Du H, Li W H. Bioparticle separation and manipulation using dielectrophoresis[J]. Sensors and Actuators A:Physical, 2007, 133(2): 329–334. DOI:10.1016/j.sna.2006.06.029 |
| [5] | Baek M K, Yu G J, Park I H. Numerical analysis and experiment for microparticle collector using dielectrophoretic force[J]. IEEE Transactions on Magnetics, 2014, 50(2): 1–4. DOI:10.1109/TMAG.2014.2301651 |
| [6] | Green N G, Jones T B. Numerical determination of the effective moments of non-spherical particles[J]. Journal of Physics D:Applied Physics, 2007, 40(1): 78–85. DOI:10.1088/0022-3727/40/1/S12 |
| [7] |
吴成东, 许可, 田孝军.
面向SWCNT-FET制造的介电泳装配技术[J]. 东北大学学报(自然科学版), 2011, 32(8): 1072–1075.
( Wu Cheng-dong, Xu ke, Tian Xiao-jun, et al. Assembly of single walled carbon nanotube field-effect transistor using dielectrophoresis[J]. Journal of Northeastern University(Natural Science), 2011, 32(8): 1072–1075. ) |
| [8] | Ogbi A, Nicolas L, Perrussel R, et al. Numerical identification of effective multipole moments of polarizable particles[J]. IEEE Transactions on Magnetics, 2012, 48(2): 675–678. DOI:10.1109/TMAG.2011.2173166 |
| [9] | Ponder E. The measurement of the diameter of erythrocytes.V.—the relation of the diameter to the thickness[J]. Quarterly Journal of Experimental Physiology, 1930, 20(1): 29–39. DOI:10.1113/expphysiol.1930.sp000484 |