传导延迟对螺旋波动力学行为的影响

采用Greenberg-Hastings离散激发介质模型研究传导延迟对螺旋波动力学行为的影响.模拟结果表明:某些情况下,延迟对螺旋波动力学行为无影响;在适当参数下,可观察到螺旋波波长和波速的改变、呼吸螺旋波、多臂螺旋波以及螺旋波和反螺旋波共存现象,对产生这些现象的物理机制作了简要的解释.     

Lévy噪声与周期力共同驱动下的螺旋波动力学

研究加性Lévy噪声与周期外力对FitzHugh-Nagumo可激系统中螺旋波动力学行为的影响.螺旋波波头的运动随外力周期在一定范围内呈规则变化,该规则变化可用相应的傅立叶谱理解,维持该规则变化的是锁频行为.Lévy噪声序列中包含着小概率的大尺度噪声,螺旋波波头运动改变主要来自于它们的影响,本文指出Lévy噪声对波头运动的影响也依赖于外力周期的取值.在适当的参数取值下,Lévy噪声的存在也能导致螺旋波的消失,这为螺旋波的控制、消除提供了一种方法.分析了系统周期与外力周期的锁定行为,给出了不同噪声强度下的Arnold舌,指出随机共振行为的存在.     

环域反馈下复Ginzburg-Landau系统的螺旋波动力学行为

针对自激振荡系统的复金兹伯格-朗道(Complex Ginzbury-Landau, 简称CGL)方程, 研究圆形环域与方形环域两种反馈控制下的螺旋波动力学。结果表明: 圆形环域反馈控制下, 螺旋波波头通常经过一段过渡漂移后进入圆形吸引子, 圆形吸引子的半径以及反馈刚启动时波头的漂移方向随环域参数呈周期性变化, 过渡漂移阶段波头轨道的平缓程度与复反馈信号模的时间函数中钟形部分的陡度有关, 且反馈增益的正负与大小也会影响受控螺旋波的动力学行为。方形环域反馈控制下的螺旋波波头的吸引子更为丰富, 主要包括方形吸引子、小的极限环吸引子、菱形吸引子以及点吸引子, 点吸引子通常位于方形环域的两条对角线上, 且波头运动随环域控制参数呈现规律性变换。     

通过被动介质耦合的两螺旋波的同步

采用Br模型研究了通过被动介质耦合的两二维可激发系统中螺旋波的同步,被动介质由可激发元素组成,这些元素之间不存在耦合.数值模拟结果表明,被动介质对螺旋波的同步有很大影响,当两系统中的初态螺旋波相同时,被动介质可导致稳定螺旋波发生漫游,螺旋波转变为螺旋波对或反靶波;当两系统中的初态螺旋波不同步时,在适当的参数下,两螺旋波可以实现同步、相同步,此外还观察到两螺旋波波头相互排斥、多螺旋波共存、同步的时空周期斑图、系统演化到静息态等现象.在被动介质中,一般可观察到波斑图,但是在某些情况下,被动介质会出现同步振荡现象.这些结果有助于人们理解心脏系统中出现的时空斑图.     

两个延迟耦合FitzHugh-Nagumo系统的动力学行为

考虑了两个延迟耦合FitzHugh-Nagumo系统，首先分析两点延迟耦合的动力学行为，给出静息态的局部稳定和不稳定的参数区. 螺旋波同步、共同的静息态以及两个系统不同的状态出现在静息态稳定参数区内，在这些态的过渡耦合强度参数区内会出现很多演化图样，它们反映了螺旋波波头、波体和波尾的不同，也反映空间尺度对螺旋波的影响.讨论了静息态局部不稳定参数区内，两点延迟耦合的动力学行为和相应参数下两个延迟耦合反应-扩散系统的斑图动力学行为. 关键词： FitzHugh-Nagumo系统 螺旋波 同步     

时空调制对可激发介质螺旋波波头动力学行为影响及控制研究

本文首先研究了时空调制对可激发介质中周期螺旋波波头动力学行为的影响. 随着时空调制的增大, 螺旋波经历了周期螺旋波、外滚螺旋波、旅行螺旋波和内滚螺旋波的显著变化. 通过定义序参量来定量的描述由时空调制引起的螺旋波在不同态之间非平衡跃迁的临界条件, 及漫游螺旋波波头圆滚圆半径随调制参数的变化情况. 当时空调制增大到某个临界值时, 螺旋波发生了破碎; 再增加时空调制, 螺旋波则发生了衰减, 系统最终演化为空间均匀静息态. 在文中给出了螺旋波发生破碎和衰减的机理和原因. 最后将时空调制方法运用于漫游螺旋波, 实现了将漫游螺旋波控制成周期螺旋波, 或将其控制为空间均匀静息态.     

动作电位动态变化对螺旋波演化行为的影响

在心脏病患者的心脏中观察到心室不同部位的心肌细胞动作电位时长(APD)恢复曲线陡峭程度差别很大,而陡峭的APD恢复曲线既可以在某些情况下导致螺旋波破碎和心室纤维性颤动,也可能在另一些情况下不导致螺旋波破碎,螺旋波动力学行为与陡峭的APD恢复曲线的关系仍未完全清楚,因此需要深入研究.本文采用二维可激发介质元胞自动机模型,研究陡峭程度不同的APD恢复曲线下螺旋波的动力学行为,数值模拟结果表明:陡峭的APD恢复曲线可以使漫游螺旋波稳定,也可以促进螺旋波漫游或引起破碎,甚至使螺旋波消失,观察到在APD恢复曲线总平均斜率大于1情况下螺旋波仍维持稳定或漫游,在APD恢复曲线总平均斜率比1小很多的情况下螺旋波出现破碎;在APD恢复曲线总平均斜率大于1的情况下观察到多普勒失稳、爱克豪斯失稳和APD交替变化三种螺旋波破碎方式,其中APD交替变化导致的螺旋波破碎可在APD恢复曲线总平均斜率比1小很多的情况下发生.观察到螺旋波通过漫游出系统边界和遇到传导障碍消失.此外,还发现通过增大元胞APD有利于防止螺旋波破碎.对产生这些现象的物理机制进行解释.     

非均匀可激发介质中的稀密螺旋波

在非均匀可激发介质中,采用Barkley模型数值模拟了稀螺旋波和密螺旋波, 并对二者的动力学行为随参数的变化进行了研究. 结果发现:稀螺旋波的旋转频率随参数b的增加迅速减小,之后趋于饱和, 显示出不同于密螺旋波的行为;两种螺旋波的周期和波长随参数ε 和非均匀区域尺寸R的增加而增加,相对稀螺旋波而言,密螺旋波的性质对R的依赖更为敏感; 稀螺旋波端点的波速随R的增加而减小,与密螺旋波波速变化趋势相反. 另外,由于非均匀区域的影响,当ε 或b 超过某一临界值时,螺旋波臂上出现缺陷点.     

时空调制引起的漫游螺旋波与旅行螺旋波共存现象

以Barkley模型为研究对象, 研究了在时空调制作用下螺旋波时空动力学行为特性. 发现在适当的调制参数的作用下, 能够在同一系统中同时观察到漫游螺旋波与旅行螺旋波. 通过数值模拟研究分析, 给出了能产生漫游螺旋波与旅行螺旋波共存现象的潜在机理, 并详细讨论了在Barkley 模型中要产生这种共存现象的两个必要条件. 关键词： 时空调制 漫游螺旋波 旅行螺旋波     

三层弱循环耦合可激发介质中螺旋波动力学

采用Bär模型研究了具有循环反馈耦合的三层可激发介质中的螺旋波动力学行为,数值模拟结果显示: 在耦合强度较小时, 在各子系统中可观察到螺旋波漂移或漫游; 当耦合强度稍大时, 相互作用既可以使螺旋波漫游或漂移出系统边界而使子系统回到静息态,还可以使子系统的螺旋波态转变为靶波或湍流态, 并观察到子系统的渐近态依赖初值现象; 继续增大耦合强度, 三个子系统的螺旋波可达到近似广义同步; 当耦合强度更大时, 螺旋波演化为湍流态.     

White noise-induced spiral waves and multiple spatial coherence resonances in a neuronal network with type I excitability

White noise-induced pattern formation is studied in a network composed of Morris–Lecar neuronal models with type I excitability and with initial values higher than that of the resting potential. The appearance and disappearance of spiral waves, as well as the transitions between spiral wave patterns with different kinds of complexity characterized by the normalized spatial autocorrelation function, enable changes in the order of the network so as to exhibit a scenario with two or more locally maximal peaks, as can be clearly seen in the signal to noise ratio curves, as the noise intensity is adjusted from small to large in a wide range. A possible physical mechanism of the multiple resonances based on the dynamics of type I excitability and initial values is provided. The potential biological significance of the noise-induced spiral waves is discussed.     

Spatial patterns in a network composed of neurons with different excitabilities induced by autapse

It has been identified that autapse can modulate dynamics of single neurons and spatial patterns of neuronal networks. In the present paper, based on the results that autapse can induce type II excitability changed to type I excitability, spatial pattern transitions are simulated in a two-dimensional neuronal network composed of excitatory coupled neurons with autapse which can induce excitability transition. Different spatial patterns including random-like pattern, irregular wave, regular wave, and nearly synchronous behavior are simulated with increasing the percentage (σ) of neurons with type I excitability. When noise is introduced, spiral waves are induced. By calculating signal-to-noise ratio from the spatial structure function and the mean firing probability of neurons, regular waves and spiral waves exhibit optimal spatial correlation, implying the occurrence of spatial coherence resonance phenomenon. The changes of mean firing probability of neurons show that different firing frequency between type I excitability and type II excitability may be an important factor to modulate the spatial patterns. The results are helpful to understand the spatial patterns including spiral waves observed in the biological experiment on the rat cortex perfused with drugs which can induce single neurons changed from type II excitability to type I excitability and block the inhibitory couplings between neurons. The excitability transition, absence of inhibitory coupling, noise as well as the autapse are important factors to modulate the spatial patterns including spiral waves.     

External forcing of spiral waves

The effect of an external rhythm on rotating spiral waves in excitable media is investigated. Parameters of the unperturbed medium were chosen, such that the organizing spiral tip describes meandering (hypocyclic) trajectories, which are the most general shape for the experimentally observed systems. Periodical modulation of excitability in a model of the Belousov-Zhabotinsky (BZ) reaction forces meandering spiral tips to describe trajectories that are not found at corresponding stationary conditions. For different modulation periods, two types of resonance drift, phase-locked tip motion, a spectrum of hypocyclic trajectories, and complex multifrequency patterns were computed. The computational results are complemented by experimental data obtained for periodically changing illumination of the photosensitive BZ reaction. The observed drastic deformation of the tip trajectory is considered as an efficient means to study and to control wave processes in excitable media.     

无扩散功能的缺陷对螺旋波动力学行为的影响

采用Br 模型,在二维激发介质中引入无扩散功能的缺陷,研究了均匀分布的缺陷对螺旋波动力学行为的影响.研究发现,缺陷导致介质的激发性降低、波传播速度减少,在一定数量的缺陷均匀分布下,缺陷可以使原来稳定的螺旋波发生漫游或破碎,缺陷使原来不稳定的螺旋波稳定或漫游,首次在激发介质中观察到螺旋波因Doppler效应破碎形成小螺旋波和时空混沌共存现象.对产生这些现象的物理机理做了简要的讨论. 关键词： 激发介质 螺旋波 缺陷     

Hindmarsh-Rose神经元阵列自发产生螺旋波的研究

采用Hindmarsh-Rose(HR)神经元模型,研究了二维神经元阵列系统从一个具有随机相位分布的初态演化最终是否能自发产生螺旋波的问题.数值结果表明:系统是否出现螺旋波与单个HR神经元的状态、系统的初态和耦合强度有关,其中单个HR神经元的振荡状态起主要作用.当单个HR神经元处于一周期振荡态时,在一定的耦合强度范围内系统都会自发出现多个螺旋波和螺旋波对,出现螺旋波与系统初态无关,只要适当选择耦合强度,在系统中可以出现单个螺旋波.当耦合强度超过某一阈值后,继续增加耦合强度,系统会呈现三种不同的动力学行为,分别与三类初态有关.系统从第一类初态演化将偶尔出现单个螺旋波,系统从第二类和第三类初态演化将分别出现间歇性全局同步振荡和振荡死亡.当单个神经元处于二周期态时,只有当系统神经元的初相位比较均匀分布时,系统才能自发出现螺旋波,而且出现螺旋波的耦合强度范围大为减少.当神经元处于更高的周期态时,系统一般不容易自发出现螺旋波.这些结果有助于人们了解大脑皮层自发产生螺旋波的机制.     

Curvature-dependent excitation propagation in cultured cardiac tissue

The geometry of excitation wave front may play an important role on the propagation block and spiral wave formation. The wave front which is bent over the critical value due to interaction with the obstacles may partially cease to propagate and appearing wave breaks evolve into rotating waves or reentry. This scenario may explain how reentry spontaneously originates in a heart. We studied highly curved excitation wave fronts in the cardiac tissue culture and found that in the conditions of normal, non-inhibited excitability the curvature effects do not play essential role in the propagation. Neither narrow isthmuses nor sharp corners of the obstacles, being classical objects for production of extremely curved wave front, affect non-inhibited wave propagation. The curvature-related phenomena of the propagation block and wave detachment from the obstacle boundary were observed only after partial suppression of the sodium channels with Lidocaine. Computer simulations confirmed the experimental observations. The explanation of the observed phenomena refers to the fact that the heart tissue is made of finite size cells so that curvature radii smaller than the cardiomyocyte size loses sense, and in non-inhibited tissue the single cell is capable to transmit excitation to its neighbors.     

Dynamics of Meandering Spiral Waves Driven by Two-Point Feedback

We numerically study the dynamics of meandering spiral waves in theexcitable system subjected to a feedback signal coming from two measuring points located on a straight line together with the initial spiral core. The core location and size radius of the final attractors are computed, and they change with the position of the moving measuring point in a unique way. By the Fourier Spectral analysis, we find the frequency-locked behaviors different from the
driving scheme of the external periodic force. It is also found that the meandering spiral wave can be eliminated when the moving measuring point approaches closely the boundary and its feedback gain is large enough. This offers an effective and convenient method for eliminating meandering spiral waves.     

Spiral-wave dynamics in excitable medium with excitability modulated by rectangle wave

We numerically study the dynamics of spiral waves in the excitable system with the excitability modulated by a rectangle wave. The tip trajectories and their variations with the modulation period T are explained by the corresponding spectrum analysis. For a large T ,the external modulation leads to the occurrence of more frequency peaks and these frequencies change with the modulation period according to their specific rules,respectively. Some of the frequencies and a primary frequency f 1 determine the corresponding curvature periods,which are locked into rational multiplies of the modulation period. These frequency-locking behaviours and the limited life-span of the frequencies in their variations with the modulation period constitute many resonant entrainment bands in the T axis. In the main bands,which follow the relation T/T 12 = m/n,the size variable R x of the tip trajectory is a monotonic increasing function of T . The rest of the frequencies are linear combinations of the two ones. Due to the complex dynamics,many unique tip trajectories appear at some certain T . We find also that spiral waves are eliminated when T is chosen from the end of the main resonant bands. This offers a useful method of controling the spiral wave.