极化电场对可激发介质中螺旋波的控制*

潘军廷 何银杰 夏远勋 张宏†

1)(浙江大学海洋学院,舟山 316021)

2)(浙江大学物理学系,浙江近代物理中心,杭州 310027)

螺旋波在不同的物理、化学和生物系统中普遍存在.周期外场,比如极化电场,尤其是具有旋转对称性的圆极化电场可对螺旋波动力学产生重要影响.本文综述了极化电场对可激发介质中螺旋波的控制,包括共振漂移、同步、手征对称性破缺、多臂螺旋波的稳定、次激发介质中的螺旋波、三维回卷波湍流态的控制、心脏组织中螺旋波的去钉扎、心脏组织中螺旋波湍流态的控制等.

1 引 言

许多生物膜、生理组织、化学反应系统都具有“可激发”的性质[1].当此类系统的局部区域处于静息状态时,对微扰是稳定的;但对于较强的扰动将有一个快速的响应,呈现激发状态,最后回到静息状态.螺旋波广泛地存在于各类可激发介质,包括Belousov-Zhabotinsky(BZ)化学反应[2]、一氧化碳在铂金表面的氧化反应[3]、心脏组织[4]等系统中.心脏实验表明,心律失常与螺旋波的自组织及螺旋波湍流态有密切关系[4,5],螺旋波的深入研究对心脏病研究将会产生重要影响.另外,与人们所熟悉的行波表现为平行的波前和波背完全一致地传播的情形不同,螺旋波的波前和波背在顶点处融为一体,形成具有奇异性结构的螺旋波端点.此端点构成了螺旋波的组织中心,是一个时空点拓扑缺陷(相位奇点),研究螺旋波可为人们研究拓扑缺陷动力学提供一种有效的途径.

电场对螺旋波影响的实验和理论工作比较丰富,包括直流电场和交流电场.在BZ反应实验中,Steinbock等[6]发现直流电场可使螺旋波的漂移速度形成分别与电场方向平行和垂直的两个分量,且其垂直分量的方向会因螺旋波的手征性不同而各异.Muñuzuri等[7]指出,当交流电场频率成两倍于螺旋波频率时,螺旋波会产生共振漂移.Zhang等[8]研究了直流、交流电场作用下的螺旋波漂移行为,得到了螺旋波漂移速度的一个近似解析公式.

上述研究的外加电场,只具有一个方向指向.由于螺旋波绕其端点旋转,具有旋转对称性,人们希望外加的电场也具有旋转对称性.Chen等[9]通过把两个相互垂直且频率相同的交流电场叠加,理论上提出了一个极化电场(图1),实现了外电场具有旋转对称性的特性.通过调节两交流电场间的相位差,可以得到不同模式的极化电场,比如线极化、椭圆极化、圆极化电场等(图2).近来,Ji等[10]在BZ反应实验中实现了外加极化电场.

图1 极化电场示意图[11] Ex,Ey 表示两个相互垂直的交流电场,E0,ωe 分别是交流电场的振幅和频率,φe,φxy分别是初相位和相位差Fig.1.Realization sketch of a polarized electric field[11]: Ex,Ey are two ac electric fields perpendicular to each other,where E0,ωe are the amplitude and the frequency of the electric field,respectively,and φe,φxy are the initial phase and the phase difference,respectively.

描述螺旋波动力学的反应扩散方程一般可表述为[12]

其中 u,v 是系统变量;f,g 代表反应函数(f,g 的具体形式可参见文献[12]);ε 描述介质的可激发性;Du,Dv是扩散系数.系统(1)式可外加一极化电场 Ex=E0cos(ωet+φe),Ey=E0cos(ωet+φe+φxy),其中 Ex,Ey表示两个相互垂直的交流电场.

如果研究的系统是化学系统,比如BZ反应,极化电场对系统的影响表现为[13]: 方程(1)中,第一个方程的右边加入第二个方程的右边加入其中 Mu,Mv表示离子迁移率.本文第2—7节所研究的系统为化学系统.

如果研究的系统是心脏系统(Dv为零),极化电场只对系统中的缺陷边界产生影响,表现为在缺陷边界有一零流边界条件[14]: n ·∇(u+E ·r)=0,其中n是缺陷边界处的法向矢量,r 为缺陷边界处的点,u 为跨膜电压,E=(Ex,Ey)为极化电场.本文第8节和第9节所研究的系统为心脏系统.

2 螺旋波的共振漂移

极化电场下,Chen 等[9]研究得出,当电场频率是螺旋波频率的2倍时,螺旋波呈直线漂移,即共振漂移.当电场为圆极化电场,且与螺旋波同向旋转时,螺旋波的漂移速率最大;反之,当圆极化电场与螺旋波反向旋转时,螺旋波不发生漂移.另外,通过改变极化电场的相位差,可改变螺旋波的漂移速度.

近来,Li等[15]基于响应函数理论[16−18],考虑方程(1)受弱电场扰动作用,研究了弱极化电场下螺旋波的漂移行为.研究解析给出了螺旋波的漂移速度公式,且该公式适用于强激发及弱激发介质(图3).当圆极化电场频率是螺旋波频率的2倍时,研究得到顺时针旋转螺旋波的漂移速度公式为

图2 不同相位差的极化电场[9]Fig.2.Polarized electric fields at different phase differences[9].

图3 无电场作用下,顺时针旋转的螺旋波[15](a)强激发介质;(b)弱激发介质Fig.3.Clockwise(cw)rotating spiral waves without electric field[15]:(a)Highly excitable medium;(b)weakly excitable medium.

3 螺旋波的同步

若周期外力使系统的振动模式受迫变为与外力驱动同频,即出现了锁频现象.锁频状态中,外力振幅与频率差呈现出舌状关联,即存在Arnold tongue区域[19].Chen等[20]数值研究了圆极化电场对螺旋波频率的影响,发现当圆极化电场和螺旋波同向旋转,且二者频率较为接近时,原来处于严格旋转或漫游状态的螺旋波均可被圆极化电场同步(图5),并与圆极化电场保持着恒定的旋转方向偏角.锁频后,圆极化电场强度与频率之间,出现了 Arnold tongue区域.

图4 极化电场作用下顺时针旋转螺旋波的漂移[15](a),(b)顺时针(φxy=0.5π)、逆时针(φxy=1.5π)旋转的圆极化电场作用下螺旋 波的漂移,其中,ωe=2ω,φe=0,Φ =0,ω 是 螺旋波的频率;(c),(d)漂移 速率与相位 差 φxy 的关系,实 线为理论结 果,圆圈为数值结果;(e),(f)漂移角与相位差 φxy 的关系,实线为理论结果,圆圈为数值结果.当漂移速率为0(φxy=1.5π),无漂移角;(a),(c),(e)强激发介质;(b),(d),(f)弱激发介质Fig.4.Drifting behaviors of cw spirals under the influence of a polarized electric field[15]:(a),(b)Drifting behaviors of spirals under the influence of a cw(φxy=0.5π)and a counterclockwise(ccw)(φxy=1.5π)circularly polarized electric fields(CPEFs)with ωe=2ω,φe=0,Φ=0,and ω being the frequency of the spiral waves;(c),(d)dependence of theoretical(lines)and numerical(circles)drift speeds on the phase difference φxy;(e),(f)dependence of theoretical(lines)and numerical(circles)drift angles on the phase difference φxy .When the drift speed is 0(φxy=1.5π),the drift angle cannot be defined.(a),(c),(e)Highly excitable medium;(b),(d),(f)Weakly excitable medium.

图5 螺旋波端点轨迹[20](a)−(e)无圆极化电场作用;(f)−(j)有圆极化电场作用,同一列上下两个图的介质尺寸相同Fig.5.Trajectories of spiral tips without control(a)−(e)and under control(f)−(j)of CPEF[20].The size in the same column is identical.

近来,Li等[21]基于响应函数理论,对 Chen 等[20]的发现进行了理论研究,解析给出了圆极化电场作用下,螺旋波的稳定区域及同步条件,从理论上得到了圆极化电场的 Arnold tongue区域(图 6),理论和数值结果符合较好.

图6 圆极化电场强度与频率之间的 Arnold tongue 区域[21]实线为理论结果,圆圈为数值结果Fig.6.Arnold tongue of the(ωe,E0)-plane of CPEF[21]:Lines and circles denote the theoretical and the numerical results,respectively.

4 手征对称性破缺

螺旋波湍流态下,系统中具有顺时针、逆时针旋转方向的螺旋波数目一般情况下一样多,系统总体上呈现出手征对称性,即零旋转手征对称性.早年,Nicolis和Prigogine[22]预言,若反应扩散系统受手征性外场作用,可能产生手征对称性破缺.近来,Li等[23]研究了圆极化电场对反应扩散系统中螺旋波湍流态手征性的影响,发现外加的圆极化电场,可使原先具有零旋转手征对称性的系统出现手征对称性破缺,且使系统具有与圆极化电场的旋转手征性一致的手征性(图7).研究利用漫游螺旋波,进一步发现其原因在于,外加的圆极化电场,可使系统中与其同向旋转的螺旋波同步,进而产生频移,从而压制住具有另外旋转方向的螺旋波,使得系统出现了手征对称性破缺(图8).对于严格旋转的螺旋波,Li等[24]理论分析和数值模拟发现,与圆极化电场同向、反向旋转的螺旋波,均可通过调节电场频率对其手征性进行控制.

5 多臂螺旋波的稳定

两个或多个相同手征性的螺旋波可以形成一个多臂螺旋波[25],其稳定性与它的周期密切相关[26].比如,当双臂螺旋波的周期小于一定临界值时,两个螺旋波之间的排斥作用会增强,使其衰退变成两个单臂螺旋波.Deng等[27]研究发现,与双臂螺旋波反向旋转的圆极化电场,可使双臂螺旋波的周期增大,进而减弱两个螺旋波之间的相互作用,从而可使其处于稳定状态(图9).进一步研究表明,外加的圆极化电场,需要具有合适的强度和频率,才可使双臂螺旋波稳定(图10给出了详细的相图).

图7 螺旋波湍流态中的手征对称性破缺[23](a)螺旋波湍流态中逆时针(黑点)、顺时针(白点)旋转的螺旋波;(b)在逆时针旋转的圆极化电场作用下,系统中仅存留逆时针旋转的螺旋波;(c)当(b)中的圆极化电场变成顺时针旋转后,系统中仅存留顺时针旋转的螺旋波.ωe=1.25,E0=0.20Fig.7.Coherent state out of defect-mediated turbulence accompanied by chiral symmetry breaking[23]:(a)An initial defect-mediated turbulence state consists of ccw spiral defects(black dots)and cw ones(white dots);(b)coherent state with only ccw spiral waves exists in the asymptotic state when the system is coupled to a ccw CPEF with ωe=1.25 and E0=0.20;(c)similar to(b)but with a cw CPEF,and in such a case,only cw spiral waves survive in the system.

6 次激发介质中的螺旋波

若介质的激发性只能维持平面波的传播,但不能让螺旋波旋转起来,则介质属于次激发介质[28].Jung和Mayer-Kress[29,30]提出,噪声作用可助次激发介质产生螺旋波.Cai等[11]发现,圆极化电场可使次激发介质中原先不断收缩的半平面波不再收缩,而是发生卷曲形成螺旋波(图11).所形成的螺旋波可处于严格旋转状态,且与圆极化电场同步.基于运动学关系[31]

图8 逆时针旋转圆极化电场作用下,漫游螺旋波对的手征对称性破缺[23](a)E0=0;(b)E0=0.10,ωe=1.350＞ω0ms(ω0ms≈ 1.309),其中 ω0ms 是无圆极化电 场作用下,漫游螺旋波的主频率;(c)E0=0.24,ωe=1.307＜ ω0ms;(d)圆极化电场作用下,逆时针旋转螺旋波频率 ωsccw(实心圆圈)及顺时针旋转螺旋波频率 ω scw(空心圆圈)与圆极化电场频率 ωe 的 关系,其 中 E0=0.10;(e)圆极化电场作用下,逆时针旋转螺旋波频率 ωsccw(实心圆圈)及顺时针旋转螺旋波频率 ωscw(空心圆圈)与圆极化电场强度 E 0 的关系,其中 ωe=1.307Fig.8.Symmetry breaking of a meandering spiral pair under a ccw CPEF[23]:(a)E0=0;(b)E0=0.10,ωe=1.350＞ω0ms(ω0ms≈ 1.309),where ω0ms is the principal frequency of the meandering spiral without the CPEF;(c)E0=0.24,ωe=1.307＜ ω0ms;(d)dependence of ωsccw(the frequency of the ccw spiral wave)(full circles)and ωscw(the frequency of the cw spiral wave)(open circles)on ωe with E0=0.10;(e)dependence of ωsccw(full circles)and ωscw(open circles)on E0 with ωe=1.307 .

图9 圆极化电场对双臂螺旋波的稳定作用[27](a)无外场作用;(b)有圆极化电场作用,E0=1.0,ωe=1.22Fig.9.Stabilization of two-armed spiral by CPEF[27]:(a)Without external fields;(b)in the presence of a CPEF with E0=1.0,ωe=1.22 .

研究给出了圆极化电场使次激发介质产生螺旋波的定量解释.方程(3)中,ct是半平面波端点的切向速度,c0表 示平面波的速度.cB=c0(B-Bc)/K,其中,表征介质的激发性,Δ=δ3-3δ由系统参数 δ 确 定;Bc=0.535 是区分次激发介质(B ＞Bc)和激发介质(B ＜Bc)的临界值;K≈0.63是常数.cE=γ//E//+γ⊥E⊥,其中 E//=E0cosθ,E⊥=E0sinθ分别是与ct平行和垂直的电场分量;γ//≈-0.850,γ⊥≈0.929是常数;θ 是 ct与电场E 之间的夹角(图 12(a)).Hakim和Karma[31]提出,方程(3)中,若 ct＜c0,则介质可产生螺旋波.Cai等[11]通过半解析分析,得出与图11中的螺旋波对应的方程(3)中,cB+cE＜0(图 12(b)),即ct＜c0,从而使得次激发介质中形成了螺旋波,且数值与半解析结果取得了一致的结果(图12(c)).值得提出的是,图11中的系统参数 B=0.5515,c0≈0.8887,从而可有 cB＞0 .因此,cB+cE＜0是圆极化电场作用的结果.

图10 圆极化电场对双臂螺旋波作用的相图[27] BU,TS分别表示破碎、稳定区域;SS表示电场过弱,不足以稳定双臂螺旋波而使其衰退为单臂螺旋波的区域;ω0=1.24是单臂螺旋波的频率Fig.10.The phase diagram for the effects of CPEF on twoarmed spiral[27]: BU,TS denote the breakup and the stabilization regions,respectively,and SS means the region where the electric field is not strong enough to stabilize the two-armed spiral and it decays into two single-armed spirals.The frequency of the single-armed spiral ω0=1.24 .

图11 次激−发介质中半平面波的演化[11](a)−(c)无外场作用;(d)(f)有圆极化电场作用.E0=0.2,ωe=0.2Fig.11.The evolution of a broken plane wave in the subexcitable system without(a)−(c)and with(d)−(f)CPEFs[11].E0=0.2,ωe=0.2 .

图12 圆极化电场使次激发介质产生螺旋波的机制分析[11],E0=0.1(a)圆极化电场对螺旋波端点的作用示意图;(b)cB+cE 随 ωe的变化关系;(c)ct的半解析解与数值解的比较Fig.12.The mechanism analyses for spiral waves sustained by CPEF in subexcitable media[11],E0=0.1 :(a)The sketch of a spiral wave tip submitted to a CPEF;(b)results of cB+cE v arying with ωe;(c)the comparison of the semi-analytical ct with the numerical ct .

7 三维回卷波湍流态的控制

回卷波会因其奇异线张力为负而失稳成湍流态[32],相反,若奇异线张力为正,则可使回卷波处于稳定状态.在系统全空间加周期信号[33],或在系统的某一局部加一个较强的周期信号[34],这种张力为负引起的回卷波湍流态可以被控制到空间均匀状态.最近,Li等[35]发现,圆极化电场可抑制因张力为负引起的回卷波湍流态,而使其从无序湍流态走向有序回卷波态(图13).其原因在于圆极化电场可使回卷波锁相,进而将奇异线张力由负变为正,以及具有较原先高些的频率.研究进一步发现,在圆极化电场频率与回卷波频率较为接近的范围内,如此的回卷波湍流被抑制现象存在一定的范围(图 14).另外,基于响应函数理论,研究发现圆极化电场下,奇异线张力满足以下关系[35]:

其中 Γ1,γ1分别表示回卷波在有、无外电场作用下的奇异线张力,E0为 电场强度,Δ ω 表征电场与回卷波之间的频差.M,g//,g⊥三个参数可由理论计算得到.数值计算与理论值取得一致的结果(图15).

图13 圆极化电场作用下,回卷波湍流态从无序走向有序[35] E0=0.4,ωe=ω0=1.2455,其中 ω0 表 示螺旋波的频率;t=0,施加圆极化电场作用;黄线表示奇异线Fig.13.Ordering of scroll wave turbulence by switching on a ccw CPEF at t=0 with E0=0.4 and rotation frequency ωe=1.2455 equal to the natural spiral wave frequency ω0[35].Filaments are shown in yellow.

图14 圆极化电场抑制(实心圆圈)回卷波湍流态的参数区域中,电场强度 E 0 与 归一化频率 ωe/ω0 的关系[35],交叉表示不能抑制回卷波湍流态Fig.14.Parameter region of scroll wave turbulence suppression(full circles)as a function of external field amplitude E0 and normalized frequency ωe/ω0[35].Crosses denote failure of ordering turbulence.

图15 锁相回卷波的奇异线张力[35]Fig.15.Filament tension of phase-locked scroll waves[35].

8 心脏组织中螺旋波的去钉扎

心脏系统中,螺旋波会引起心动过速,若螺旋波被缺陷钉扎住而无法驱赶走,则可能导致更严重的心脏疾病.为了使螺旋波去钉扎,人们研究发现,匀强电场在缺陷边界处可产生激发波[36−39],进而可使螺旋波去钉扎.Feng等[40]发现,相比于匀强电场,圆极化电场作用下缺陷边界的膜电势分布呈动态的“太极图”状(匀强电场下,呈静态的对称偶极子状),且与圆极化电场同步旋转(图16).研究发现,圆极化电场能有效地使螺旋波去钉扎(图17).由于圆极化电场作用产生的“太极图”式的膜电势分布,即去极化及超极化区域可与电场动态同步旋转,使得螺旋波的去钉扎较匀强电场更有效、所需的外加电压更低,以及有更大的适用范围(图18).

9 心脏组织中螺旋波湍流态的控制

心脏组织中的螺旋波湍流态,会引起严重的心脏病,比如心颤[41].人们通常采用单个高压电击作用于心脏而使其恢复正常心电活动[42−44],但因其副作用太大[43,44]而局限过大.缺陷激发波法[36−39]利用低电压匀强电场在缺陷边界产生激发波,从而可以去颤.但该方法只能在缺陷邻近处有效而使其效率受限(需要较多缺陷,才能产生出较多的激发波去颤).Feng等[45]发现,圆极化电场在低电压下,可使缺陷组织激发形成不断向外传播的圆形波,且其频率可高于湍流波的频率(图19).这些具有较高频率的圆形波可抑制住心脏组织中的湍流态,从而可以去颤(图20).相对而言,匀强电场虽也可使缺陷组织激发出圆形波,但其频率几乎都是低于湍流波的频率(图21),从而很难实现对湍流态的有效控制.

图16 圆极化电场(a)(c)和匀强电场(b)(d)作用下的膜电势分布[40](a)Luo-Rudy 模型,E0=0.05V/cm,ωe=0.2 rad/ms;(b)Luo-Rudy 模型,E0=0.05V/cm;(c)Barkley模型,E0=0.05,ωe=4;(d)Barkley 模型,E0=0.05 .R 表示缺陷半径,Luo-Rudy 模型中,R=0.32cm;Barkley模型中,R=3 .红色点箭头表示电场方向,红色曲线箭头表示圆极化电场逆时针旋转.围绕着缺陷的红色、蓝色分别表示去极化、超极化区域Fig.16.Distribution of the membrane potential induced by CPEF and uniform electric field(UEF)[40]:(a)CPEF in Luo-Rudy model,E0=0.05V/cm,ωe=0.2rad/ms;(b)UEF in Luo-Rudy model,E0=0.05V/cm;(c)CPEF in Barkley model,E0=0.05,ωe=4;(d)UEF in Barkley model,E0=0.05 .In Luo-Rudy model,the obstacle size R=0.32cm,and in Barkley model,R=3 .The red dotted arrows represent the directions of electric fields.The red curved arrows mean CPEFs rotate counterclockwise.The red and the blue regions around obstacles demonstrate de-polarizations and hyper-polarizations,respectively.

图17 圆极化电场去除顺时针旋转的钉扎螺旋波[40](a)Luo-Rudy模型,螺旋波频率ωs=0.136rad/ms,圆极化电场E0=0.7V/cm,ωe=0.1rad/ms;电场作用时间 t=0—46.2 ms,圆极化电场相对于 x 轴的初相位为 φe;螺旋波相对于 x 轴的初相位 φs=0;(b)Barkley模型,螺旋波频率ωs=1.024,圆极化电场E0=1.8,ωe=3.686;电场作用时间t=0—6.N和N' 表示不同时刻圆极化电场激发产生的激发波,S和S' 分别表示初始的钉扎螺旋波、去钉扎后的螺旋波.白色箭头表示波的传播方向Fig.17.Unpinning the cw rotating anchored spiral by CPEF[40]:(a)Luo-Rudy model,the frequency of spiral ωs=0.136rad/ms;E0=0.7V/cm,ωe=0.1rad/ms;CPEF is applied from t=0 to t=46.2ms .φe is the initial phase of CPEF relative to x axis;φs is the initial phase of the anchored spiral front relative to x axis and sets as zero;(b)Barkley model,ωs=1.024,E0=1.8,ωe=3.686;CPEF is applied from t=0 to t=6 .N and N' represent different new waves nucleated by CPEF in different time.S and S' represent the initial anchord spiral and the new free spiral,respectively.White arrows are the propagation directions of waves.

图18 Barkley 模型中,圆极化电场(灰色区域)和匀强电场(阴影区域)去除钉扎螺旋波的适用范围[40] SW,NW,RW和BI分别表示螺旋波、无波、收缩波和双稳区域;圆极化电场 E0=1.8,ωe=3.6ωs;匀强电场E0=7Fig.18.Unpinning scope of CPEF(gray)and UEF(shaded)in Barkley model[40]: SW,NW,RW and BI regions represent spiral waves,no wave,retracting waves and bi-stability,respectively;for CPEF,E0=1.8,ωe=3.6ωs and for UEF,E0=7 .

图19 静息态下圆极化电场及其激发产生的圆形波列之间的频率关系[45] E0=1.0V/cm,0.065rad/ms≤ ωCPEF≤0.22 rad/ms;ωcir 是 圆形波列的频率;ωtur 是螺旋波湍流态的主频率Fig.19.The frequency relations between the circular wave train and CPEF in a two-dimensional quiescent medium[45]:E0=1.0V/cm,0 .065rad/ms≤ ωCPEF≤ 0.22 rad/ms;ωcir is the frequency of the circular wave trains,and ωtur is the dominant frequency of the spiral turbulence.

图20 圆极化电场抑制螺旋波湍流态[45]E0=1.0V/cm,ωCPEF=0.14rad/ms(a) t =0;(b)t=1000 ms;(c)t=1800 ms;(d)t=2800 msFig.20.Suppression of spiral turbulence by CPEF[45]:E0=1.0V/cm,ωCPEF=0.14rad/ms :(a)t=0;(b)t=1000 ms;(c)t=1800 ms;(d)t=2800 ms.

图21 静息态下匀强电场及其激发产生的圆形波列之间的频率关系[45] E0=1.0V/cm,0.065rad/ms≤ ωUEF≤0.22 rad/ms;脉冲间隔为10ms;ωcir是圆形波列的频率;ωtur是螺旋波湍流态的主频率Fig.21.The frequency relations between the circular wave train and UEF in a two-dimensional quiescent medium[45]:E0=1.0V/cm,0 .065rad/ms≤ ωCPEF≤ 0.22rad/ms;the pulse duration is 1 0ms;ωcir is the frequency of the circular wave trains,and ωtur is the dominant frequency of the spiral turbulence.

10 结 论

本文着重介绍了圆极化电场对可激发介质中螺旋波的控制.由于具有旋转对称性,圆极化电场可对螺旋波动力学产生独特而有效的影响,并可基于运动学、响应函数等理论,给出相应的理论解释.

近来,人们通过对两对相互垂直的电极板通以交流电,在BZ反应实验中实现了外加极化电场[10].心脏实验方面,人们一般把两对正负电极分别放置于心脏的四周,然后通以直流电来对心脏产生作用(文献[38]的图5(D)).很显然,只要把其中的直流电换成交流电,就可以类似于BZ反应实验一样,在心脏实验中实现外加极化电场.因此,极化电场在心脏病研究方面,也具有重要的意义.希望本文介绍的数值及理论结果可以在实验上实现.