1.Ocean College, Zhejiang University, Zhoushan 316021, China 2.Zhejiang Institute of Modern Physics, Department of Physics, Zhejiang University, Hangzhou 310027, China
Abstract:Spiral waves are ubiquitous in diverse physical, chemical, and biological systems. Periodic external fields, such as polarized electric fields, especially circularly polarized electric fields which possess rotation symmetry may have significant effects on spiral wave dynamics. In this paper, control of spiral waves in excitable media under polarized electric fields is reviewed, including resonant drift, synchronization, chiral symmetry breaking, stabilization of multiarmed spiral waves, spiral waves in subexcitable media, control of scroll wave turbulence, unpinning of spiral waves in cardiac tissues, control of spiral wave turbulence in cardiac tissues, etc. Keywords:polarized electric fields/ excitable media/ spiral waves/ control
全文HTML
--> --> -->
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]${E_x}$, ${E_y}$表示两个相互垂直的交流电场, ${E_0}$, ${\omega _{\rm{e}}}$分别是交流电场的振幅和频率, ${\phi _{\rm{e}}}$, ${\phi _{xy}}$分别是初相位和相位差 Figure1. Realization sketch of a polarized electric field[11]: ${E_x}$, ${E_y}$ are two ac electric fields perpendicular to each other, where ${E_0}$, ${\omega _{\rm{e}}}$ are the amplitude and the frequency of the electric field, respectively, and ${\phi _{\rm{e}}}$, ${\phi _{xy}}$ are the initial phase and the phase difference, respectively.
图 2 不同相位差的极化电场[9] Figure2. Polarized electric fields at different phase differences[9].
其中$| {\dot R} |$表示螺旋波的漂移速率; $\varTheta $是漂移角; $\varPhi $是螺旋波的初相位; ${\mu _2}$, ${\upsilon _2}$是漂移系数, 可由响应函数理论计算得到. 从方程(2)可看出, 当圆极化电场与螺旋波同向旋转时(${\phi _{xy}} = 0.5{\text{π}}$), 螺旋波的漂移速率$| {\dot R} |$最大; 当圆极化电场与螺旋波反向旋转时(${\phi _{xy}} = 1.5{\text{π}}$), 螺旋波漂移速率$| {\dot R}| = 0$. 当$\varPhi $或${\phi _{\rm{e}}}$改变时, 漂移角$\varTheta $会均匀改变, 且$\Delta \varTheta = 2\Delta \varPhi $及$\Delta \varTheta = \Delta {\phi _{\rm{e}}}$. 当${\phi _{xy}}$改变时, 漂移的方向也会发生改变, 且$\Delta \varTheta = 0.5\Delta {\phi _{xy}}$. 研究表明, 理论结果与数值结果定量符合(图4). 图 4 极化电场作用下顺时针旋转螺旋波的漂移[15] (a),(b)顺时针(${\phi _{xy}} = 0.5{\text{π}}$)、逆时针(${\phi _{xy}} = 1.5{\text{π}}$)旋转的圆极化电场作用下螺旋波的漂移, 其中, ${\omega _{\rm{e}}} = 2\omega $, ${\phi _{\rm{e}}} = 0$, $\varPhi = 0$, $\omega $是螺旋波的频率; (c),(d)漂移速率与相位差${\phi _{xy}}$的关系, 实线为理论结果, 圆圈为数值结果; (e),(f)漂移角与相位差${\phi _{xy}}$的关系, 实线为理论结果, 圆圈为数值结果. 当漂移速率为0(${\phi _{xy}} = 1.5{\text{π}}$), 无漂移角; (a),(c),(e)强激发介质; (b),(d),(f)弱激发介质 Figure4. 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 (${\phi _{xy}} = 0.5{\text{π}}$) and a counterclockwise (ccw) (${\phi _{xy}} = 1.5{\text{π}}$) circularly polarized electric fields (CPEFs) with ${\omega _{\rm{e}}} = 2\omega $, ${\phi _{\rm{e}}} = 0$, $\varPhi = 0$, and $\omega $ being the frequency of the spiral waves; (c),(d) dependence of theoretical (lines) and numerical (circles) drift speeds on the phase difference ${\phi _{xy}}$; (e),(f) dependence of theoretical (lines) and numerical (circles) drift angles on the phase difference ${\phi _{xy}}$. When the drift speed is 0 (${\phi _{xy}} = 1.5{\text{π}}$), the drift angle cannot be defined. (a),(c),(e) Highly excitable medium; (b),(d),(f) Weakly excitable medium.
3.螺旋波的同步若周期外力使系统的振动模式受迫变为与外力驱动同频, 即出现了锁频现象. 锁频状态中, 外力振幅与频率差呈现出舌状关联, 即存在Arnold tongue区域[19]. Chen等[20]数值研究了圆极化电场对螺旋波频率的影响, 发现当圆极化电场和螺旋波同向旋转, 且二者频率较为接近时, 原来处于严格旋转或漫游状态的螺旋波均可被圆极化电场同步(图5), 并与圆极化电场保持着恒定的旋转方向偏角. 锁频后, 圆极化电场强度与频率之间, 出现了Arnold tongue区域. 图 5 螺旋波端点轨迹[20] (a)?(e)无圆极化电场作用; (f)?(j)有圆极化电场作用, 同一列上下两个图的介质尺寸相同 Figure5. 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]实线为理论结果, 圆圈为数值结果 Figure6. Arnold tongue of the (${\omega _{\rm{e}}}, \;{E_0}$)-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]理论分析和数值模拟发现, 与圆极化电场同向、反向旋转的螺旋波, 均可通过调节电场频率对其手征性进行控制. 图 7 螺旋波湍流态中的手征对称性破缺[23] (a)螺旋波湍流态中逆时针(黑点)、顺时针(白点)旋转的螺旋波; (b)在逆时针旋转的圆极化电场作用下, 系统中仅存留逆时针旋转的螺旋波; (c)当(b)中的圆极化电场变成顺时针旋转后, 系统中仅存留顺时针旋转的螺旋波. ${\omega _{\rm{e}}} = 1.25$, E0 = 0.20 Figure7. 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 ${\omega _{\rm{e}}} = 1.25$ and ${E_0} = 0.20$; (c) similar to (b) but with a cw CPEF, and in such a case, only cw spiral waves survive in the system.
图 8 逆时针旋转圆极化电场作用下, 漫游螺旋波对的手征对称性破缺[23] (a) ${E_0} = 0$; (b) ${E_0} = 0.10$, ${\omega _{\rm{e}}} = 1.350 > $${\omega _0}^{{\rm{ms}}}\;({\omega _0}^{{\rm{ms}}} \approx 1.309) $, 其中${\omega _0}^{{\rm{ms}}}$是无圆极化电场作用下, 漫游螺旋波的主频率; (c) ${E_0} = 0.24$, ${\omega _{\rm{e}}} = 1.307 < {\omega _0}^{{\rm{ms}}}$; (d)圆极化电场作用下, 逆时针旋转螺旋波频率${\omega _{\rm{s}}}^{{\rm{ccw}}}$(实心圆圈)及顺时针旋转螺旋波频率${\omega _{\rm{s}}}^{{\rm{cw}}}$(空心圆圈)与圆极化电场频率${\omega _{\rm{e}}}$的关系, 其中${E_0} = 0.10$; (e)圆极化电场作用下, 逆时针旋转螺旋波频率${\omega _{\rm{s}}}^{{\rm{ccw}}}$(实心圆圈)及顺时针旋转螺旋波频率${\omega _{\rm{s}}}^{{\rm{cw}}}$(空心圆圈)与圆极化电场强度${E_0}$的关系, 其中${\omega _{\rm{e}}} = 1.307$ Figure8. Symmetry breaking of a meandering spiral pair under a ccw CPEF[23]: (a) ${E_0} = 0$; (b) ${E_0} = 0.10$, ${\omega _{\rm{e}}} = 1.350 > $${\omega _0}^{{\rm{ms}}}\;({\omega _0}^{{\rm{ms}}} \approx 1.309) $, where ${\omega _0}^{{\rm{ms}}}$ is the principal frequency of the meandering spiral without the CPEF; (c) E0 = 0.24, ${\omega _{\rm{e}}} = 1.307 < {\omega _0}^{{\rm{ms}}}$; (d) dependence of ${\omega _{\rm{s}}}^{{\rm{ccw}}}$(the frequency of the ccw spiral wave) (full circles) and ${\omega _{\rm{s}}}^{{\rm{cw}}}$(the frequency of the cw spiral wave) (open circles) on ${\omega _{\rm{e}}}$ with ${E_0} = 0.10$;(e) dependence of ${\omega _{\rm{s}}}^{{\rm{ccw}}}$(full circles) and ${\omega _{\rm{s}}}^{{\rm{cw}}}$(open circles) on ${E_0}$ with ${\omega _{\rm{e}}} = 1.307$.
5.多臂螺旋波的稳定两个或多个相同手征性的螺旋波可以形成一个多臂螺旋波[25], 其稳定性与它的周期密切相关[26]. 比如, 当双臂螺旋波的周期小于一定临界值时, 两个螺旋波之间的排斥作用会增强, 使其衰退变成两个单臂螺旋波. Deng等[27]研究发现, 与双臂螺旋波反向旋转的圆极化电场, 可使双臂螺旋波的周期增大, 进而减弱两个螺旋波之间的相互作用, 从而可使其处于稳定状态(图9). 进一步研究表明, 外加的圆极化电场, 需要具有合适的强度和频率, 才可使双臂螺旋波稳定(图10给出了详细的相图). 图 9 圆极化电场对双臂螺旋波的稳定作用[27] (a)无外场作用; (b)有圆极化电场作用, ${E_0} = 1.0$, ${\omega _{\rm{e}}} = 1.22$ Figure9. Stabilization of two-armed spiral by CPEF[27]: (a) With-out external fields; (b) in the presence of a CPEF with ${E_0} = 1.0$, ${\omega _{\rm{e}}} = 1.22$.
图 10 圆极化电场对双臂螺旋波作用的相图[27] BU, TS分别表示破碎、稳定区域; SS表示电场过弱, 不足以稳定双臂螺旋波而使其衰退为单臂螺旋波的区域; ${\omega _0} = 1.24$是单臂螺旋波的频率 Figure10. The phase diagram for the effects of CPEF on two-armed spiral[27]: BU, TS denote the breakup and the stabi-lization 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 ${\omega _0} = 1.24$.
6.次激发介质中的螺旋波若介质的激发性只能维持平面波的传播, 但不能让螺旋波旋转起来, 则介质属于次激发介质[28]. Jung 和Mayer-Kress[29,30]提出, 噪声作用可助次激发介质产生螺旋波. Cai等[11]发现, 圆极化电场可使次激发介质中原先不断收缩的半平面波不再收缩, 而是发生卷曲形成螺旋波(图11). 所形成的螺旋波可处于严格旋转状态, 且与圆极化电场同步. 基于运动学关系[31] 图 11 次激发介质中半平面波的演化[11] (a)?(c)无外场作用; (d)?(f)有圆极化电场作用. ${E_0} = 0.2$, ${\omega _{\rm{e}}} = 0.2$ Figure11. The evolution of a broken plane wave in the subexcitable system without (a)?(c) and with (d)?(f) CPEFs[11]. ${E_0} = 0.2$, ${\omega _{\rm{e}}} = 0.2$.
${c_{\rm t}} = {c_0} + {c_B} + {c_E},$
研究给出了圆极化电场使次激发介质产生螺旋波的定量解释. 方程(3)中, ${c_{\rm{t}}}$是半平面波端点的切向速度, ${c_0}$表示平面波的速度. ${c_B} = {c_0}(B - {B_{\rm{c}}})/K$, 其中, $B = 4\sqrt 3 \varepsilon /{\varDelta ^3}$表征介质的激发性, $\varDelta = {\delta ^3} - 3\delta $由系统参数$\delta $确定; ${B_{\rm{c}}} = 0.535$是区分次激发介质($B > {B_{\rm{c}}}$)和激发介质($B < {B_{\rm{c}}}$)的临界值; $K \approx 0.63$是常数. ${c_E} = {\gamma _{//}}{E_{//}} + {\gamma _ \bot }{E_ \bot }$, 其中${E_{//}} = {E_0}\cos \theta $, ${E_ \bot } = {E_0}\sin \theta $分别是与${c_{\rm{t}}}$平行和垂直的电场分量; ${\gamma _{//}} \approx - 0.850$, ${\gamma _ \bot } \approx 0.929$是常数; $\theta $是${c_{\rm{t}}}$与电场E之间的夹角(图12(a)). Hakim 和Karma[31]提出, 方程(3)中, 若${c_{\rm{t}}} < {c_0}$, 则介质可产生螺旋波. Cai等[11]通过半解析分析, 得出与图11中的螺旋波对应的方程(3)中, ${c_B} + {c_E} < 0$(图12(b)), 即${c_{\rm{t}}} < {c_0}$, 从而使得次激发介质中形成了螺旋波, 且数值与半解析结果取得了一致的结果(图12(c)). 值得提出的是, 图11中的系统参数$B = 0.5515$, ${c_0} \approx 0.8887$, 从而可有${c_B} > 0$. 因此, ${c_B} + {c_E} < 0$是圆极化电场作用的结果. 图 12 圆极化电场使次激发介质产生螺旋波的机制分析[11], ${E_0} = 0.1$ (a)圆极化电场对螺旋波端点的作用示意图; (b)${c_B} + {c_E}$随${\omega _{\rm{e}}}$的变化关系; (c)${c_{\rm{t}}}$的半解析解与数值解的比较 Figure12. The mechanism analyses for spiral waves sustained by CPEF in subexcitable media[11], ${E_0} = 0.1$: (a) The sketch of a spiral wave tip submitted to a CPEF; (b) results of ${c_B} + {c_E}$varying with ${\omega _{\rm{e}}}$; (c) the comparison of the semi-analytical ${c_{\rm{t}}}$ with the numerical ${c_{\rm{t}}}$.
7.三维回卷波湍流态的控制回卷波会因其奇异线张力为负而失稳成湍流态[32], 相反, 若奇异线张力为正, 则可使回卷波处于稳定状态. 在系统全空间加周期信号[33], 或在系统的某一局部加一个较强的周期信号[34], 这种张力为负引起的回卷波湍流态可以被控制到空间均匀状态. 最近, Li等[35]发现, 圆极化电场可抑制因张力为负引起的回卷波湍流态, 而使其从无序湍流态走向有序回卷波态(图13). 其原因在于圆极化电场可使回卷波锁相, 进而将奇异线张力由负变为正, 以及具有较原先高些的频率. 研究进一步发现, 在圆极化电场频率与回卷波频率较为接近的范围内, 如此的回卷波湍流被抑制现象存在一定的范围(图14). 另外, 基于响应函数理论, 研究发现圆极化电场下, 奇异线张力满足以下关系[35]: 图 13 圆极化电场作用下, 回卷波湍流态从无序走向有序[35]${E_0} = 0.4$, ${\omega _{\rm{e}}} = {\omega _0} = 1.2455$, 其中${\omega _0}$表示螺旋波的频率; $t = 0$, 施加圆极化电场作用; 黄线表示奇异线 Figure13. Ordering of scroll wave turbulence by switching on a ccw CPEF at $t = 0$ with ${E_0} = 0.4$ and rotation frequency ${\omega _{\rm{e}}} = 1.2455$ equal to the natural spiral wave frequency ${\omega _0}$[35]. Filaments are shown in yellow.
图 14 圆极化电场抑制(实心圆圈)回卷波湍流态的参数区域中, 电场强度${E_0}$与归一化频率${\omega _{\rm{e}}}/{\omega _0}$的关系[35], 交叉表示不能抑制回卷波湍流态 Figure14. Parameter region of scroll wave turbulence suppression (full circles) as a function of external field amplitude ${E_0}$ and normalized frequency ${\omega _{\rm{e}}}/{\omega _0}$[35]. Crosses denote failure of ordering turbulence.
8.心脏组织中螺旋波的去钉扎心脏系统中, 螺旋波会引起心动过速, 若螺旋波被缺陷钉扎住而无法驱赶走, 则可能导致更严重的心脏疾病. 为了使螺旋波去钉扎, 人们研究发现, 匀强电场在缺陷边界处可产生激发波[36-39], 进而可使螺旋波去钉扎. Feng等[40]发现, 相比于匀强电场, 圆极化电场作用下缺陷边界的膜电势分布呈动态的“太极图”状(匀强电场下, 呈静态的对称偶极子状), 且与圆极化电场同步旋转(图16). 研究发现, 圆极化电场能有效地使螺旋波去钉扎(图17). 由于圆极化电场作用产生的“太极图”式的膜电势分布, 即去极化及超极化区域可与电场动态同步旋转, 使得螺旋波的去钉扎较匀强电场更有效、所需的外加电压更低, 以及有更大的适用范围(图18). 图 16 圆极化电场(a)(c)和匀强电场(b)(d)作用下的膜电势分布[40] (a) Luo-Rudy模型, ${E_0} = 0.05\;{\rm{V/cm}}$, ${\omega _{\rm{e}}} =$ 0.2 rad/ms; (b) Luo-Rudy模型, ${E_0} = 0.05\;{\rm{V/cm}}$; (c) Barkley模型, ${E_0} = 0.05$, ${\omega _{\rm{e}}} = 4$; (d) Barkley模型, ${E_0} = 0.05$. R表示缺陷半径, Luo-Rudy模型中, $R = 0.32\;{\rm{cm}}$; Barkley模型中, $R = 3$. 红色点箭头表示电场方向, 红色曲线箭头表示圆极化电场逆时针旋转. 围绕着缺陷的红色、蓝色分别表示去极化、超极化区域 Figure16. Distribution of the membrane potential induced by CPEF and uniform electric field (UEF)[40]: (a) CPEF in Luo-Rudy model, ${E_0} = 0.05\;{\rm{V/cm}}$, ${\omega _{\rm{e}}} = 0.2\;{\rm{rad/ms}}$; (b) UEF in Luo-Rudy model, ${E_0} = 0.05\;{\rm{V/cm}}$; (c) CPEF in Barkley model, ${E_0} = 0.05$, ${\omega _{\rm{e}}} = 4$; (d) UEF in Barkley model, ${E_0} = 0.05$.In Luo-Rudy model, the obstacle size $R = 0.32\;{\rm{cm}}$, 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模型, 螺旋波频率${\omega _{\rm{s}}} = 0.136\;{\rm{rad/ms}}$, 圆极化电场${E_0} = 0.7\;{\rm{V/cm}}$, ${\omega _{\rm{e}}} = 0.1\;{\rm{rad/ms}}$; 电场作用时间t = 0—46.2 ms, 圆极化电场相对于x轴的初相位为${\phi _{\rm{e}}}$; 螺旋波相对于x轴的初相位${\phi _{\rm{s}}} = 0$; (b) Barkley模型, 螺旋波频率${\omega _{\rm{s}}} = 1.024$, 圆极化电场${E_0} = 1.8$, ${\omega _{\rm{e}}} = 3.686$; 电场作用时间t = 0—6. N和N' 表示不同时刻圆极化电场激发产生的激发波, S和S' 分别表示初始的钉扎螺旋波、去钉扎后的螺旋波. 白色箭头表示波的传播方向 Figure17. Unpinning the cw rotating anchored spiral by CPEF[40]: (a) Luo-Rudy model, the frequency of spiral ${\omega _{\rm{s}}} = 0.136\;{\rm{rad/ms}}$; ${E_0} = 0.7\;{\rm{V/cm}}$, ${\omega _{\rm{e}}} = 0.1\;{\rm{rad/ms}}$; CPEF is applied from $t = 0$ to $t = 46.2\;{\rm{ms}}$. ${\phi _{\rm{e}}}$ is the initial phase of CPEF relative to x axis; ${\phi _{\rm{s}}}$ is the initial phase of the anchored spiral front relative to x axis and sets as zero; (b) Barkley model, ${\omega _{\rm{s}}} = 1.024$, ${E_0} = 1.8$, ${\omega _{\rm{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分别表示螺旋波、无波、收缩波和双稳区域; 圆极化电场${E_0} = 1.8$, ${\omega _{\rm{e}}} = 3.6\;{\omega _{\rm{s}}}$; 匀强电场${E_0} = 7$ Figure18. 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, ${E_0} = 1.8$, ${\omega _{\rm{e}}} = 3.6\;{\omega _{\rm{s}}}$ and for UEF, ${E_0} = 7$.
9.心脏组织中螺旋波湍流态的控制心脏组织中的螺旋波湍流态, 会引起严重的心脏病, 比如心颤[41]. 人们通常采用单个高压电击作用于心脏而使其恢复正常心电活动[42-44], 但因其副作用太大[43,44]而局限过大. 缺陷激发波法[36-39]利用低电压匀强电场在缺陷边界产生激发波, 从而可以去颤. 但该方法只能在缺陷邻近处有效而使其效率受限(需要较多缺陷, 才能产生出较多的激发波去颤). Feng等[45]发现, 圆极化电场在低电压下, 可使缺陷组织激发形成不断向外传播的圆形波, 且其频率可高于湍流波的频率(图19). 这些具有较高频率的圆形波可抑制住心脏组织中的湍流态, 从而可以去颤(图20). 相对而言, 匀强电场虽也可使缺陷组织激发出圆形波, 但其频率几乎都是低于湍流波的频率(图21), 从而很难实现对湍流态的有效控制. 图 19 静息态下圆极化电场及其激发产生的圆形波列之间的频率关系[45]${E_0} = 1.0\;{\rm{V/cm}}$, $0.065\;{\rm{rad/ms}} \leqslant {\omega _{{\rm{CPEF}}}} \leqslant $ 0.22 rad/ms; ${\omega _{{\rm{cir}}}}$是圆形波列的频率; ${\omega _{{\rm{tur}}}}$是螺旋波湍流态的主频率 Figure19. The frequency relations between the circular wave train and CPEF in a two-dimensional quiescent medium[45]: ${E_0} = 1.0\;{\rm{V/cm}}$, $0.065\;{\rm{rad/ms}} \leqslant {\omega _{{\rm{CPEF}}}} \leqslant $ 0.22 rad/ms; ${\omega _{{\rm{cir}}}}$ is the frequency of the circular wave trains, and ${\omega _{{\rm{tur}}}}$ is the dominant frequency of the spiral turbulence.
图 20 圆极化电场抑制螺旋波湍流态[45]${E_0} = 1.0\;{\rm{V/cm}}$, ${\omega _{{\rm{CPEF}}}} = 0.14\;{\rm{rad/ms}}$ (a) t = 0; (b) t = 1000 ms; (c) t = 1800 ms; (d) t = 2800 ms Figure20. Suppression of spiral turbulence by CPEF[45]: ${E_0} = 1.0\;{\rm{V/cm}}$, ${\omega _{{\rm{CPEF}}}} = 0.14\;{\rm{rad/ms}}$: (a) t = 0; (b) t = 1000 ms; (c) t = 1800 ms; (d) t = 2800 ms.
图 21 静息态下匀强电场及其激发产生的圆形波列之间的频率关系[45]${E_0} = 1.0\;{\rm{V/cm}}$, $0.065\;{\rm{rad/ms}} \leqslant {\omega _{{\rm{UEF}}}} \leqslant $ 0.22 rad/ms; 脉冲间隔为$10\;{\rm{ms}}$; ${\omega _{{\rm{cir}}}}$是圆形波列的频率; ${\omega _{{\rm{tur}}}}$是螺旋波湍流态的主频率 Figure21. The frequency relations between the circular wave train and UEF in a two-dimensional quiescent medium[45]: ${E_0} = 1.0\;{\rm{V/cm}}$, $0.065\;{\rm{rad/ms}} \leqslant {\omega _{{\rm{CPEF}}}} \leqslant 0.22\;{\rm{rad/ms}}$; the pulse duration is $10\;{\rm{ms}}$; ${\omega _{{\rm{cir}}}}$ is the frequency of the circular wave trains, and ${\omega _{{\rm{tur}}}}$ is the dominant frequency of the spiral turbulence.
图 1 极化电场示意图[11] 
图 2 不同相位差的极化电场[9]
图 3 无电场作用下, 顺时针旋转的螺旋波[15] (a)强激发介质; (b)弱激发介质
图 4 极化电场作用下顺时针旋转螺旋波的漂移[15] (a),(b)顺时针(
图 5 螺旋波端点轨迹[20] (a)?(e)无圆极化电场作用; (f)?(j)有圆极化电场作用, 同一列上下两个图的介质尺寸相同
图 6 圆极化电场强度与频率之间的Arnold tongue区域[21]实线为理论结果, 圆圈为数值结果
图 7 螺旋波湍流态中的手征对称性破缺[23] (a)螺旋波湍流态中逆时针(黑点)、顺时针(白点)旋转的螺旋波; (b)在逆时针旋转的圆极化电场作用下, 系统中仅存留逆时针旋转的螺旋波; (c)当(b)中的圆极化电场变成顺时针旋转后, 系统中仅存留顺时针旋转的螺旋波. 
图 8 逆时针旋转圆极化电场作用下, 漫游螺旋波对的手征对称性破缺[23] (a) 
图 9 圆极化电场对双臂螺旋波的稳定作用[27] (a)无外场作用; (b)有圆极化电场作用, 
图 10 圆极化电场对双臂螺旋波作用的相图[27] BU, TS分别表示破碎、稳定区域; SS表示电场过弱, 不足以稳定双臂螺旋波而使其衰退为单臂螺旋波的区域; 
图 11 次激发介质中半平面波的演化[11] (a)?(c)无外场作用; (d)?(f)有圆极化电场作用. 
图 12 圆极化电场使次激发介质产生螺旋波的机制分析[11], 
图 13 圆极化电场作用下, 回卷波湍流态从无序走向有序[35] 
图 14 圆极化电场抑制(实心圆圈)回卷波湍流态的参数区域中, 电场强度
图 15 锁相回卷波的奇异线张力[35]
图 16 圆极化电场(a)(c)和匀强电场(b)(d)作用下的膜电势分布[40] (a) Luo-Rudy模型, 
图 17 圆极化电场去除顺时针旋转的钉扎螺旋波[40] (a) Luo-Rudy模型, 螺旋波频率
图 18 Barkley模型中, 圆极化电场(灰色区域)和匀强电场(阴影区域)去除钉扎螺旋波的适用范围[40] SW, NW, RW和BI分别表示螺旋波、无波、收缩波和双稳区域; 圆极化电场
图 19 静息态下圆极化电场及其激发产生的圆形波列之间的频率关系[45] 
图 20 圆极化电场抑制螺旋波湍流态[45] 
图 21 静息态下匀强电场及其激发产生的圆形波列之间的频率关系[45] 
