双层耦合介质中四边形图灵斑图的数值研究

采用双层线性耦合Lengyel-Epstein模型,在二维空间对简单正四边和超点阵四边形进行了数值分析.结果表明:当两子系统波数比N1时,随耦合强度的增大,基模的波矢空间共振形式发生改变,系统由简单六边形自发演化为结构复杂的新型斑图,除已报道的超六边形外,还获得了简单正四边和多种超点阵四边形,包括大小点、点线、白眼和环状超四边等斑图.当耦合系数α和β在一定范围内同步增大时,两子系统形成相同波长的Ⅰ型简单正四边;当α和β不同步增大时,由于两图灵模在短波子系统形成共振,系统斑图经相变发生Ⅰ型正四边→Ⅱ型正四边→超点阵四边形的转变;当系统失去耦合作用时,短波子系统波长为λ的Ⅰ型正四边斑图迅速失稳并形成波长为λ/N的Ⅰ型正四边,随模拟时间的延长,两子系统中不同波长的正四边均会经相变发生Ⅰ型正四边→Ⅱ型正四边→六边形的转变.     

Effect of a Local Source or Sink of Inhibitor on Turing Patterns

We investigate a reaction-diffusion model in which a Turing pattern develops and reproduces the formation of periodic segments behind a propagating chemical wave front. The chemical scheme involves two species known as activator and inhibitor. The model can be used to mimic the formation of prevertebrae during the early development of vertebrate embryo. Deterministic and stochastic analyses of the reaction-diffusion processes are performed for two typical sets of parameter values, far from and close to the Turing bifurcation. The effects of a local source or sink of inhibitor on the growing structure are studied and successfully compared with experiments performed on chick embryos. We show that fluctuations may lead to the formation of additional prevertebra.     

Turing pattern selection in a reactionben diffusion epidemic modelast

We present Turing pattern selection in a reaction—diffusion epidemic model under zero-flux boundary conditions. The value of this study is twofold. First, it establishes the amplitude equations for the excited modes, which determines the stability of amplitudes towards uniform and inhomogeneous perturbations. Second, it illustrates all five categories of Turing patterns close to the onset of Turing bifurcation via numerical simulations which indicates that the model dynamics exhibits complex pattern replication: on increasing the control parameter v, the sequence “H₀ hexagons → H₀-hexagon-stripe mixtures → stripes → H_π-hexagon-stripe mixtures → H_π hexagons” is observed. This may enrich the pattern dynamics in a diffusive epidemic model.     

Chaotic Turing pattern formation in spatiotemporal systems

The problem of Turing pattern formation has attracted much attention in nonlinear science as well as physics, chemistry and biology. So far spatially ordered Turing patterns have been observed in stationary and oscillatory media only. In this paper we find that spatially ordered Turing patterns exist in chaotic extended systems. And chaotic Turing patterns are strikingly rich and surprisingly beautiful with their space structures. These findings are in sharp contrast with the intuition of pseudo-randomness of chaos. The richness and beauty of the chaotic Turing patterns are attributed to a large variety of symmetry properties realized by various types of self-organizations of partial chaos synchronizations.     

反应扩散模型在图灵斑图中的应用及数值模拟

反应扩散方程模型常被用于描述生物学中斑图的形成.从反应扩散模型出发,理论推导得到GiererMeinhardt模型的斑图形成机理,解释了非线性常微分方程系统的稳定常数平衡态在加入扩散项后会发生失稳并产生图灵斑图的过程.通过计算该模型,得到图灵斑图产生的参数条件.数值方法中采用一类有效的高精度数值格式,即在空间离散条件下采用Chebyshev谱配置方法,在时间离散条件下采用紧致隐积分因子方法.该方法结合了谱方法和紧致隐积分因子方法的优点,具有精度高、稳定性好、存储量小等优点.数值模拟表明,在其他条件一定的情况下,系统控制参数κ取不同值对于斑图的产生具有重要的影响,数值结果验证了理论结果.     

Turing instability for a semi-discrete Gierer-Meinhardt system

In this paper, we investigate the spatial patterns of a Gierer-Meinhardt system where the space is discrete in two dimensions with the periodic boundary condition and time is continuous, in contrast to the continuum models. The conditions of Turing instability are obtained by linear analysis and a series of numerical simulations are performed. In the instability region, we have shown that this system can produce a number of different patterns such as stripes and snowflake pattern, other than ubiquitous fix-spotted patterns. As mentioned, the results are substantiated only by means of snapshots of the apatial grid. However, we also give some analysis by using the time series at three random grids and of the average value of states, that is, the stable state patterns can be observed. On the other hand, the effects of varying parameters on pattern formation are also discussed. Moreover, simulations for fixed parameters and special initial conditions indicate that the initial conditions can alter the structure of patterns. The patterns can form as a consequence of cellular interaction. So patterns arising from a semi-discrete model can present simulations on a geometrically accurate representation in biology. As a result, our work is interesting and important in ecology.     

Superlattice Patterns in Coupled Turing Systems

In this paper, superlattice patterns have been investigated by using a two linearly coupled Brusselator model. It is found that superlattice patterns can only be induced in the sub-system with the short wavelength. Three different coupling methods have been used in order to investigate the mode interaction between the two Turing modes. It is proved in the simulations that interaction between activators in the two sub-systems leads to spontaneous formation of black eye pattern and/or white eye patterns while interaction between inhibitors leads to spontaneous formation of super-hexagonal pattern. It is also demonstrated that the same symmetries of the two modes and suitable wavelength ratio of the two modes should also be satisfied to form superlattice patterns.     

Note on a universal quantum Turing machine

In this Letter, we construct a novel model of universal quantum Turing machine (QTM) by means of a property of Riemann zeta function, which is free from the specific time for an input data and efficiently simulates each step of a given QTM.     

双层非线性耦合反应扩散系统中复杂Turing斑图

采用双层耦合的Brusselator模型, 研究了两个子系统非线性耦合时Turing 模对斑图的影响, 发现两子系统Turing 模的波数比和耦合系数的大小对斑图的形成起着重要作用. 模拟结果表明: 斑图类型随波数比值的增加, 从简单斑图发展到复杂斑图; 非线性耦合项系数在0–0.1时, 系统1中短波模在系统2失稳模的影响下不仅可形成简单六边形、四边形和条纹斑图, 两模共振耦合还可以形成蜂窝六边形、超六边形和复杂的黑眼斑图等超点阵图形, 首次在一定范围内调整控制参量观察到由简单正四边形向超六边形斑图的转化过程; 耦合系数在0.1–1时, 系统1中短波模与系统2失稳模未发生共振耦合仅观察到与系统2相同形状的简单六边形、四边形和条纹斑图. 关键词： Brusselator模型 非线性耦合 Turing模     

Turing bifurcation in a reaction-diffusion system with density-dependent dispersal

Motivated by the recent finding [N. Kumar, G.M. Viswanathan, V.M. Kenkre, Physica A 388 (2009) 3687] that the dynamics of particles undergoing density-dependent nonlinear diffusion shows sub-diffusive behaviour, we study the Turing bifurcation in a two-variable system with this kind of dispersal. We perform a linear stability analysis of the uniform steady state to find the conditions for the Turing bifurcation and compare it with the standard Turing condition in a reaction-diffusion system, where dispersal is described by simple Fickian diffusion. While activator-inhibitor kinetics are a necessary condition for the Turing instability as in standard two-variable systems, the instability can occur even if the diffusion constant of the inhibitor is equal to or smaller than that of the activator. We apply these results to two model systems, the Brusselator and the Gierer-Meinhardt model.     

Pacemakers in a Reaction-Diffusion Mechanics System

Non-linear waves of excitation are found in various biological, physical and chemical systems and are often accompanied by deformations of the medium. In this paper, we numerically study wave propagation in a deforming excitable medium using a two-variable reaction-diffusion system coupled with equations of continuum mechanics. We study the appearance and dynamics of different excitation patterns organized by pacemakers that occur in the medium as a result of deformation. We also study the interaction of several pacemakers with each other and the characteristics of pacemakers in the presence of heterogeneities in the medium. We found that mechanical deformation not only induces pacemakers, but also has a pronounced effect on spatial organization of various excitation patterns. We show how these effects are modulated by the size of the medium, the location of the initial stimulus, and the properties of the reaction-diffusion-mechanics feedback.     

Boundary Effects on the Structural Stability of Stationary Patterns in a Bistable Reaction-Diffusion System

We study a piecewise linear version of a one-component, two-dimensional bistable reaction-diffusion system subjected to partially reflecting boundary conditions, with the aim of analyzing the structural stability of its stationary patterns. Dirichlet and Neumann boundary conditions are included as limiting cases. We find a critical line in the space of the parameters which divides different dynamical behaviors. That critical line merges as the locus of the coalescence of metastable and unstable nonuniform structures.     

Turing instability for a two-dimensional Logistic coupled map lattice

In this Letter, stability analysis is applied to a two-dimensional Logistic coupled map lattice with the periodic boundary conditions. The conditions of Turing instability are obtained, and various patterns can be exhibited by numerical simulations in the Turing instability region. For example, space-time periodic structures, periodic or quasiperiodic traveling wave solutions, stationary wave solutions, spiral waves, and spatiotemporal chaos, etc. have been observed. In particular, the different pattern structures have also been observed for same parameters and different initial values. That is, pattern structures also depend on the initial values. The similar patterns have also been seen in relevant references. However, the present Letter owes to pattern formation via diffusion-driven instabilities because the system is stable in the absence of diffusion.     

Turing pattern selection in a reaction-diffusion epidemic model

We present Turing pattern selection in a reaction-diffusion epidemic model under zero-flux boundary conditions.The value of this study is twofold.First,it establishes the amplitude equations for the excited modes,which determines the stability of amplitudes towards uniform and inhomogeneous perturbations.Second,it illustrates all five categories of Turing patterns close to the onset of Turing bifurcation via numerical simulations which indicates that the model dynamics exhibits complex pattern replication:on increasing the control parameter ν,the sequence "H 0 hexagons → H 0-hexagon-stripe mixtures → stripes → H π-hexagon-stripe mixtures → H π hexagons" is observed.This may enrich the pattern dynamics in a diffusive epidemic model.     

无限方势阱本征态的稳定性

我们用存活率(SP)的概念研究了无限方势阱本征态的稳定性,做了理论上的推导,并进行了数值计算,发现基态在外界的扰动下非常稳定,而激发态则非常不稳定,能级越高,越不稳定.     

Bifurcations in a Bistable Reaction-Diffusion System

The bifurcations of stationary solutions of bistable reaction-diffusion systems and especially the model due to Schlögl are investigated. The paper is based on the observation that these dynamical systems belong to the gradient type, i.e. methods of topological investigation of the corresponding potential may be applied. Three different types of models are discussed: 1. Homogeneous model; 2. Compartmental model with diffusion; 3. Continuous model with diffusion. Bifurcation maps and relations to catastrophe theory are given.     

Spiral and Antispiral Waves in Reaction-Diffusion Systems

Spiral waves are ubiquitous phenomena in nonlinear chemical, physical, and biological systems. But antispiral waves are infrequent to date. The transition between spiral and antispiral waves has been rarely explored. We have analyzed the extended Brusselator model and the extended Oregonator model by linear stability analysis. We have demonstrated that it is possible and plausible to realize the transition between them by control of diffusion coefficient of inactivator from theoretical analysis and numerical simulations.     

Spiral and Antispiral Waves in Reaction-Diffusion Systems

Spiral waves are ubiquitous phenomena in nonlinear chemical, physical, and biological systems. But antispiral waves are infrequent to date. The transition between spiral and antispiral waves has been rarely explored. We have analyzed the extended Brusselator model and the extended Oregonator model by linear stability analysis. We have demonstrated that it is possible and plausible to realize the transition between them by control of diffusion coefficient of inactivator from theoretical analysis and numerical simulations.