From labyrinthine patterns to spiral turbulence

The nonlinear initiation of side-branching by activator-inhibitor-substrate (Turing) morphogenesis

An understanding of the underlying mechanism of side-branching is paramount in controlling and/or therapeutically treating mammalian organs, such as lungs, kidneys, and glands. Motivated by an activator-inhibitor-substrate approach that is conjectured to dominate the initiation of side-branching in a pulmonary vascular pattern, I demonstrate a distinct transverse front instability in which new fingers grow out of an oscillatory breakup dynamics at the front line without any typical length scale. These two features are attributed to unstable peak solutions in 1D that subcritically emanate from Turing bifurcation and that exhibit repulsive interactions. The results are based on a bifurcation analysis and numerical simulations and provide a potential strategy toward also developing a framework of side-branching for other biological systems, such as plant roots and cellular protrusions.

Interplay between exogenous and endogenous factors in seasonal vegetation oscillations

A fundamental question in ecology is whether vegetation oscillations are merely a result of periodic environmental variability, or rather driven by endogenous factors. We address this question using a mathematical model of dryland vegetation subjected to annual rainfall periodicity. We show that while spontaneous oscillations do not exist in realistic parameter ranges, resonant response to periodic precipitation is still possible due to the existence of damped oscillatory modes. Using multiple time-scale analysis, in a restricted parameter range, we find that these endogenous modes can be pumped by the exogenous precipitation forcing to form sustained oscillations. The oscillations amplitude shows a resonance peak that depends on model parameters representing species traits and mean annual precipitation. Extending the study to bistability ranges of uniform vegetation and bare soil, we investigate numerically the implications of resonant oscillations for ecosystem function. We consider trait parameters that represent species with damped oscillatory modes and species that lack such modes, and compare their behaviors. We find that the former are less resilient to droughts, suffer from larger declines in their biomass production as the precipitation amplitude is increased, and, in the presence of spatial disturbances, are likely to go through abrupt collapse to bare soil, rather than gradual, domino-like collapse.

Implications of tristability in pattern-forming ecosystems

Many ecosystems show both self-organized spatial patterns and multistability of possible states. The combination of these two phenomena in different forms has a significant impact on the behavior of ecosystems in changing environments. One notable case is connected to tristability of two distinct uniform states together with patterned states, which has recently been found in model studies of dryland ecosystems. Using a simple model, we determine the extent of tristability in parameter space, explore its effects on the system dynamics, and consider its implications for state transitions or regime shifts. We analyze the bifurcation structure of model solutions that describe uniform states, periodic patterns, and hybrid states between the former two. We map out the parameter space where these states exist, and note how the different states interact with each other. We further focus on two special implications with ecological significance, breakdown of the snaking range and complex fronts. We find that the organization of the hybrid states within a homoclinic snaking structure breaks down as it meets a Maxwell point where simple fronts are stationary. We also discover a new series of complex fronts between the uniform states, each with its own velocity. We conclude with a brief discussion of the significance of these findings for the dynamics of regime shifts and their potential control.

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that is triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media with applications to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact for cardiac arrhythmias.

A Theoretical Model of Jigsaw-Puzzle Pattern Formation by Plant Leaf Epidermal Cells

Plant leaf epidermal cells exhibit a jigsaw puzzle-like pattern that is generated by interdigitation of the cell wall during leaf development. The contribution of two ROP GTPases, ROP2 and ROP6, to the cytoskeletal dynamics that regulate epidermal cell wall interdigitation has already been examined; however, how interactions between these molecules result in pattern formation remains to be elucidated. Here, we propose a simple interface equation model that incorporates both the cell wall remodeling activity of ROP GTPases and the diffusible signaling molecules by which they are regulated. This model successfully reproduces pattern formation observed in vivo, and explains the counterintuitive experimental results of decreased cellulose production and increased thickness. Our model also reproduces the dynamics of three-way cell wall junctions. Therefore, this model provides a possible mechanism for cell wall interdigitation formation in vivo.

Multi-stability of circadian phase wave within early postnatal suprachiasmatic nucleus

The suprachiasmatic nucleus (SCN) is a group of cells that functions as a biological master clock. In different SCN cells, oscillations of biochemical markers such as the expression-level of clock genes, are not synchronized but instead form slow circadian phase waves propagating over the whole cell population spatio-temporal struc- ture is a fixed property set by the anatomy of a given SCN. Here, we show that this is not the case in early postnatal SCN. Earlier studies presumed that their Based on bioluminescence imaging experiments with Per2-Luciferase mice SCN cultures which guided computer simulations of a realistic model of the SCN, we demonstrate that the wave is not unique but can be in various modes including phase- coherent oscillation, crescent-shaped wave, and most notably, a rotating pinwheel wave that conceptually resembles a wall clock with a rotating hand. Furthermore, mode transitions can be induced by a pulse of 38.5 °C temperature perturbation. Importantly, the waves support a significantly different period, suggesting that neither a spatially-fixed phase ordering nor a specialized pacemaker having a fixed period exist in these studied SCNs. These results lead to new important questions of what the observed multi-stability means for the proper function of an SCN and its arrhythmia.

Unstable spiral waves and local Euclidean symmetry in a model of cardiac tissue

This paper investigates the properties of unstable single-spiral wave solutions arising in the Karma model of two-dimensional cardiac tissue. In particular, we discuss how such solutions can be computed numerically on domains of arbitrary shape and study how their stability, rotational frequency, and spatial drift depend on the size of the domain as well as the position of the spiral core with respect to the boundaries. We also discuss how the breaking of local Euclidean symmetry due to finite size effects as well as the spatial discretization of the model is reflected in the structure and dynamics of spiral waves. This analysis allows identification of a self-sustaining process responsible for maintaining the state of spiral chaos featuring multiple interacting spirals.

Influence of the medium's dimensionality on defect-mediated turbulence

Spatiotemporal chaos in oscillatory and excitable media is often characterized by the presence of phase singularities called defects. Understanding such defect-mediated turbulence and its dependence on the dimensionality of a given system is an important challenge in nonlinear dynamics. This is especially true in the context of ventricular fibrillation in the heart, where the importance of the thickness of the ventricular wall is contentious. Here, we study defect-mediated turbulence arising in two different regimes in a conceptual model of excitable media and investigate how the statistical character of the turbulence changes if the thickness of the medium is changed from (quasi-) two- dimensional to three dimensional. We find that the thickness of the medium does not have a significant influence in, far from onset, fully developed turbulence while there is a clear transition if the system is close to a spiral instability. We provide clear evidence that the observed transition and change in the mechanism that drives the turbulent behavior is purely a consequence of the dimensionality of the medium. Using filament tracking, we further show that the statistical properties in the three-dimensional medium are different from those in turbulent regimes arising from filament instabilities like the negative line tension instability. Simulations also show that the presence of this unique three-dimensional turbulent dynamics is not model specific.

A mathematical model for mechanotransduction at the early steps of suture formation

Growth and patterning of craniofacial sutures is subjected to the effects of mechanical stress. Mechanotransduction processes occurring at the margins of the sutures are not precisely understood. Here, we propose a simple theoretical model based on the orientation of collagen fibres within the suture in response to local stress. We demonstrate that fibre alignment generates an instability leading to the emergence of interdigitations. We confirm the appearance of this instability both analytically and numerically. To support our model, we use histology and synchrotron X-ray microtomography and reveal the fine structure of fibres within the sutural mesenchyme and their insertion into the bone. Furthermore, using a mouse model with impaired mechanotransduction, we show that the architecture of sutures is disturbed when forces are not interpreted properly. Finally, by studying the structure of sutures in the mouse, the rat, an actinopterygian (Polypterus bichir) and a placoderm (Compagopiscis croucheri), we show that bone deposition patterns during dermal bone growth are conserved within jawed vertebrates. In total, these results support the role of mechanical constraints in the growth and patterning of craniofacial sutures, a process that was probably effective at the emergence of gnathostomes, and provide new directions for the understanding of normal and pathological suture fusion.

Speed of traveling fronts in a sigmoidal reaction-diffusion system

We study a sigmoidal version of the FitzHugh-Nagumo reaction-diffusion system based on an analytic description using piecewise linear approximations of the reaction kinetics. We completely describe the dynamics of wave fronts and discuss the properties of the speed equation. The speed diagrams show front bifurcations between branches with one, three, or five fronts that differ significantly from the classical FitzHugh-Nagumo model. We examine how the number of fronts and their speed vary with the model parameters. We also investigate numerically the stability of the front solutions in a case when five fronts exist.

Pattern formation controlled by time-delayed feedback in bistable media

Effects of time-delayed feedback on pattern formation are studied both numerically and theoretically in a bistable reaction-diffusion model. The time-delayed feedback applied to the activator and/or the inhibitor alters the behavior of the nonequilibrium Ising-Bloch (NIB) bifurcation. If the intensities of the feedbacks applied to the two species are identical, only the velocities of Bloch fronts are changed. If the intensities are different, both the critical point of the NIB bifurcation and the threshold of stability of front to transverse perturbations are changed. The effect of time-delayed feedback on the activator opposes the effect of time-delayed feedback on the inhibitor. When the time-delayed feedback is applied individually to one of the species, positive and negative feedbacks make the bifurcation point shift to different directions. The time-delayed feedback provides a flexible way to control the NIB bifurcation and the pattern formation.

Wave propagation in a FitzHugh-Nagumo-type model with modified excitability

We examine a generalized FitzHugh-Nagumo (FHN) type model with modified excitability derived from the diffusive Morris-Lecar equations for neuronal activity. We obtain exact analytic solutions in the form of traveling waves using a piecewise linear approximation for the activator and inhibitor reaction terms. We study the existence and stability of waves and find that the inhibitor species exhibits different types of wave forms (fronts and pulses), while the activator wave maintains the usual kink (front) shape. The nonequilibrium Ising-Bloch bifurcation for the wave speed that occurs in the FHN model, where the control parameter is the ratio of inhibitor to activator time scales, persists when the strength of the inhibitor nonlinearity is taken as the bifurcation parameter.

Evolving mazes from images

We propose a novel reaction diffusion (RD) simulator to evolve image-resembling mazes. The evolved mazes faithfully preserve the salient interior structures in the source images. Since it is difficult to control the generation of desired patterns with traditional reaction diffusion, we develop our RD simulator on a different computational platform, cellular neural networks. Based on the proposed simulator, we can generate the mazes that exhibit both regular and organic appearance, with uniform and/or spatially varying passage spacing. Our simulator also provides high controllability of maze appearance. Users can directly and intuitively "paint" to modify the appearance of mazes in a spatially varying manner via a set of brushes. In addition, the evolutionary nature of our method naturally generates maze without any obvious seam even though the input image is a composite of multiple sources. The final maze is obtained by determining a solution path that follows the user-specified guiding curve. We validate our method by evolving several interesting mazes from different source images.

Mechanism of skull suture maintenance and interdigitation

Skull sutures serve as growth centers whose function involves multiple molecular pathways. During periods of brain growth the sutures remain thin and straight, later developing complex fractal interdigitations that provide interlocking strength. The nature of the relationship between the molecular interactions and suture pattern formation is not understood. Here we show that by classifying the molecules involved into two groups, stabilizing factors and substrate molecules, complex molecular networks can be modeled by a simple two-species reaction-diffusion model that recapitulates all the known behavior of suture pattern formation. This model reproduces the maintenance of thin sutural tissue at early stages, the later modification of the straight suture to form osseous interdigitations, and the formation of fractal structures. Predictions from the model are in good agreement with experimental observations, indicating that the model captures the essential nature of the interdigitation process.

Complex-anisotropy-induced pattern formation in bistable media

A construct of anisotropy in bistable media is adopted to characterize the effects of anisotropy on pattern formation by means of anisotropic line tension. A velocity curvature relation is further derived to account for the anisotropic wave propagations. Stability analysis of transverse perturbations indicates that a sufficiently strong complex anisotropy can induce dynamical instabilities and even lead to a breakup of the wave patterns. Numerical simulations show that complex anisotropy can induce rich spatiotemporal behaviors in bistable media. The results of analysis and simulations demonstrate that this method successfully incorporates complex anisotropy into the reaction diffusion model and has general significance.

Breakup of propagating waves through the development of a transient unexcitable regime

Spontaneous breakup of propagating waves was investigated in this paper with a ferroin-catalyzed Belousov-Zhabotinsky reaction modified with the inclusion of a second substrate, 1,4-cyclohexanedione (CHD). The presence of CHD, which forms a separate chemical oscillator with acidic bromate, led to the development of a scattered unexcitable regime in the studied medium, where as waves passed through these unexcitable media, the propagation was disrupted, causing the creation of free ends. It thus presents a new avenue through which a chemical wave breaks up to form spirals. By manipulating the concentrations of CHD, sulfuric acid, and bromate, unexcitable regimes of different sizes with different survival times were obtained. Kinetic data on wave speed prior to and after the unexcitable window illustrates that the occurrence of an unexcitable regime is not due to depletion of components needed for wave propagation.

Nucleation, drift, and decay of phase bubbles in period-2 oscillatory wave trains in a reaction-diffusion system

This is a report on experimental observations of phase bubbles, simply closed boundaries between domains oscillating 2pi out of phase, associated with period-2 oscillatory traveling waves in a Belousov-Zhabotinsky reaction-diffusion system. The bubbles nucleate spontaneously through a fast localized phase slip, drift radially away from a neighboring spiral wave core in an oscillatory fashion, and gradually shrink to disappear. Their oscillatory drift along the radial direction is a consequence of "period adaptation," while their lateral shrinkage is an attribute of local curvature. Similar dynamic structures can be reproduced in a simple, three-species reaction-diffusion model that supports period-2 oscillatory wave trains.

Dynamics of domain walls in pattern formation with traveling-wave forcing

We study dynamics of domain walls in pattern forming systems that are externally forced by a moving space-periodic modulation close to 2:1 spatial resonance. The motion of the forcing induces nongradient dynamics, while the wave number mismatch breaks explicitly the chiral symmetry of the domain walls. The combination of both effects yields an imperfect nonequilibrium Ising-Bloch bifurcation, where all kinks (including the Ising-like one) drift. Kink velocities and interactions are studied within the generic amplitude equation. For nonzero mismatch, a transition to traveling bound kink-antikink pairs and chaotic wave trains occurs.

Noise-induced spatiotemporal patterns in a bistable reaction-diffusion system: photoelectron emission microscopy experiments and modeling of the oxidation reaction on Ir(111)

We use photoelectron emission microscopy (PEEM) measurements to study the spatiotemporal patterns obtained for the CO oxidation reaction on Ir(111) as a function of the noise strength we superpose on the CO and the oxygen fractions of the constant total reactant gas flux. The investigations are focused on the bistable regime this reaction displays including its monostable vicinity. Simultaneously we analyze numerically the underlying reaction-diffusion (RD) equations in two spatial dimensions. For intrinsic and/or small strength of the external noise we find transitions from the locally stable to the globally stable branch via slow nucleation and growth of islands of the globally stable state: oxygen or CO, respectively. With increasing noise strength the number of islands as well as their growth rate increases. These phenomena are very well reproduced by numerical calculations of the RD model. For sufficiently large noise strength we observe bursts from CO rich to oxygen rich and back as well as switching between the two states. While such phenomena are also obtained from the model calculations, their experimentally observed spatial scales were not satisfactorily reproduced using the same approach as for the lower noise strengths.

Anomalous nonequilibrium Ising-Bloch bifurcation in discrete systems

We investigate the self-organization of two-dimensional activator-inhibitor discrete bistable systems in the neighborhood of a nonequilibrium Ising-Bloch bifurcation. The system exhibits an anomalous transition--induced by discretization--whose signature is the coexistence of Ising and Bloch walls for selected values of the spatial coupling. After curvature reduction of Bloch walls, coexistence gives rise to a unique and striking spatiotemporal dynamics: Bloch walls drive the motion and Ising walls play the role of "extended defects" oriented along the background grid directions. Strong enough external noise asymptotically restores the scenario found in the continuum limit.

Segmented waves from a spatiotemporal transverse wave instability

We observe traveling waves emitted from Turing spots in the chlorine dioxide-iodine-malonic acid reaction. The newborn waves are continuous, but they break into segments as they propagate, and the propagation of these segments ultimately gives rise to spatiotemporal chaos. We model the wave-breaking process and the motion of the chaotic segments. We find stable segmented spirals as well. We attribute the segmentation to an interaction between front rippling via a transverse instability and front symmetry breaking by a fast-diffusing inhibitor far from the codimension-2 Hopf-Turing bifurcation, and the chaos to a secondary instability of the periodic segmentation.

BLOCH-front turbulence in a periodically forced Belousov-Zhabotinsky reaction

Experiments on a periodically forced Belousov-Zhabotinsky chemical reaction show front breakup into a state of spatiotemporal disorder involving continual events of spiral-vortex nucleation and destruction. Using the amplitude equation for forced oscillatory systems and the normal form equations for a curved front line, we identify the mechanism of front breakup and explain the experimental observations.

Interaction of Ising-Bloch fronts with Dirichlet boundaries

We study the Ising-Bloch bifurcation in two systems, the complex Ginzburg Landau equation (CGLE) and a FitzHugh Nagumo (FN) model in the presence of spatial inhomogeneity introduced by Dirichlet boundary conditions. It is seen that the interaction of fronts with boundaries is similar in both systems, establishing the generality of the Ising-Bloch bifurcation. We derive reduced dynamical equations for the FN model that explain front dynamics close to the boundary. We find that front dynamics in a highly nonadiabatic (slow front) limit is controlled by fixed points of the reduced dynamical equations, that occur close to the boundary.

Transverse instability of line defects of period-2 spiral waves

Spiral waves that arise in period-2 oscillatory media extended in space generically bear "line defects" along which the local kinetics exhibits a period-1 oscillation. Locally, these defect structures can be viewed as a front separating two period-2 oscillatory domains oscillating 2pi out of phase. Here we show that their shape can become sinusoidal with a transverse instability as in bistable fronts. This instability eventually leads to a line-defect filled spatiotemporal chaotic state having erratic proliferations, annihilations, and regenerations of line defects. The same sequence of phenomena is observed in a model reaction-diffusion system as well as in an experimental system.

Vortex-pair dynamics in anisotropic bistable media: a kinematic approach

In isotropic bistable media, a vortex pair typically evolves into rotating spiral waves. In an anisotropic system, instead of spiral waves, the vortices can form wave fragments that propagate with a constant speed in a given direction determined by the system's anisotropy. The fragments may propagate invariably, shrink, or expand. We develop a kinematic approach for the study of vortex-pair dynamics in anisotropic bistable media and use it to capture the wave fragment dynamics.

Transition from spiral waves to defect-mediated turbulence induced by gradient effects in a reaction-diffusion system

The transition from spiral waves to defect-mediated turbulence was studied in a spatial open reactor using Belousov-Zhabotinsky reaction. The experimental results show a new mechanism of the transition from spirals to spatiotemporal chaos, in which the gradient effects in the three-dimensional system are essential. The transition scenario consists of two stages: first, the effects of gradients in the third dimension cause a splitting of the spiral tip and a deletion of certain wave segments, generating new wave sources; second, the waves sent by the new wave sources undergo a backfire instability, and the back waves are laterally unstable. As a result, defects are automatically generated and fill all over the system. The result of numerical simulation using the FitzHugh-Nagumo model essentially agrees with the experimental observation.

Spiral waves in a coupled network of sine-circle maps

A coupled two-dimensional lattice of sine-circle maps is investigated numerically as a simple model for coupled network of nonlinear oscillators under a spatially uniform, temporally periodic, external forcing. Various patterns, including quasiperiodic spiral waves, periodic, banded spiral waves in several different polygonal shapes, and domain patterns, are observed. The banded spiral waves and domain patterns match well with the results of earlier experimental studies. Several transitions are analyzed. Among others, the source-sink transition of a quasiperiodic spiral wave and the cascade of "side-doubling" bifurcations of polygonal spiral waves are of great interest.

Dash waves in a reaction-diffusion system

Patterns in reaction-diffusion systems generally consist of smooth traveling waves or of stationary, discontinuous Turing structures. Hybrid patterns that blend the properties of waves and Turing structures have not previously been observed. We report observation of dash waves, which consist of wave segments regularly separated by gaps, moving coherently in the Belousov-Zhabotinsky system dispersed in water-in-oil microemulsion. Dash waves emerge from the interaction between excitable and pseudo-Turing-unstable steady states. We are able to generate dash waves in simulations with simple models.

Theory of domain patterns in systems with long-range interactions of Coulomb type

We develop a theory of the domain patterns in systems with competing short-range attractive interactions and long-range repulsive Coulomb interactions. We take an energetic approach, in which patterns are considered as critical points of a mean-field free energy functional. Close to the microphase separation transition, this functional takes on a universal form, allowing us to treat a number of diverse physical situations within a unified framework. We use asymptotic analysis to study domain patterns with sharp interfaces. We derive an interfacial representation of the pattern's free energy which remains valid in the fluctuating system, with a suitable renormalization of the Coulomb interaction's coupling constant. We also derive integro-differential equations describing stationary domain patterns of arbitrary shapes and their thermodynamic stability, coming from the first and second variations of the interfacial free energy. We show that the length scale of a stable domain pattern must obey a certain scaling law with the strength of the Coulomb interaction. We analyzed the existence and stability of localized (spots, stripes, annuli) and periodic (lamellar, hexagonal) patterns in two dimensions. We show that these patterns are metastable in certain ranges of the parameters and that they can undergo morphological instabilities leading to the formation of more complex patterns. We discuss nucleation of the domain patterns by thermal fluctuations and pattern formation scenarios for various thermal quenches. We argue that self-induced disorder is an intrinsic property of the domain patterns in the systems under consideration.

Pattern formation on anisotropic and heterogeneous catalytic surfaces

We review experimental and theoretical work addressing pattern formation on anisotropic and heterogeneous catalytic surfaces. These systems are typically modeled by reaction-diffusion equations reflecting the kinetics and transport of the involved chemical species. Here, we demonstrate the influence of anisotropy and heterogeneity in a simplified model, the FitzHugh-Nagumo equations. Anisotropy causes stratification of labyrinthine patterns and spiral defect chaos in bistable media. For heterogeneous media, we study the situation where the heterogeneity appears on a length scale shorter than the typical pattern length scale. Homogenization, i.e., computation of effective medium properties, is applied to an example and illustrated with simulations in one (fronts) and two dimensions (spirals). We conclude with a discussion of open questions and promising directions that comprise the coupling of the microscopic structure of the surface to the macroscopic concentration patterns and the fabrication of nanostructures with heterogeneous surfaces as templates. (c) 2002 American Institute of Physics.

Coexistence of large amplitude stationary structures in a model of reaction-diffusion system

The two-variable reaction-diffusion model of a chemical system describing the spatiotemporal evolution to large amplitude stationary periodical structures in a one-dimensional open, continuous-flow, unstirred reactor is investigated. Numerical solutions show that the structures are generated by divisions of the traveling impulse and its stopping at the boundary of the system. Analyses of projections of numerical solutions on the phase plane of two variables elaborated in the present paper allow qualitative explanation of the results. The coexistence of the large amplitude stationary periodical structures is shown. A number of coexisting structures grows strongly with increasing length of the reactor and may be as large as one wishes. The relationship of these results to biological systems is stressed.

Four-phase patterns in forced oscillatory systems

We investigate pattern formation in self-oscillating systems forced by an external periodic perturbation. Experimental observations and numerical studies of reaction-diffusion systems and an analysis of an amplitude equation are presented. The oscillations in each of these systems entrain to rational multiples of the perturbation frequency for certain values of the forcing frequency and amplitude. We focus on the subharmonic resonant case where the system locks at one-fourth the driving frequency, and four-phase rotating spiral patterns are observed at low forcing amplitudes. The spiral patterns are studied using an amplitude equation for periodically forced oscillating systems. The analysis predicts a bifurcation (with increasing forcing) from rotating four-phase spirals to standing two-phase patterns. This bifurcation is also found in periodically forced reaction-diffusion equations, the FitzHugh-Nagumo and Brusselator models, even far from the onset of oscillations where the amplitude equation analysis is not strictly valid. In a Belousov-Zhabotinsky chemical system periodically forced with light we also observe four-phase rotating spiral wave patterns. However, we have not observed the transition to standing two-phase patterns, possibly because with increasing light intensity the reaction kinetics become excitable rather than oscillatory.

Front propagation and pattern formation in anisotropic bistable media

The effects of diffusion anisotropy on pattern formation in bistable media are studied using a FitzHugh-Nagumo reaction-diffusion model. A relation between the normal velocity of a front and its curvature is derived and used to identify distinct spatiotemporal patterns induced by the diffusion anisotropy. In a wide parameter range anisotropy is found to have an ordering effect: initial patterns evolve into stationary or breathing periodic stripes parallel to one of the principal axes. In a different parameter range, anisotropy is found to induce spatiotemporal chaos confined to one space dimension, a state we term "stratified chaos."

Phase dynamics of nearly stationary patterns in activator-inhibitor systems

The slow dynamics of nearly stationary patterns in a FitzHugh-Nagumo model are studied using a phase dynamics approach. A Cross-Newell phase equation describing slow and weak modulations of periodic stationary solutions is derived. The derivation applies to the bistable, excitable, and Turing unstable regimes. In the bistable case stability thresholds are obtained for the Eckhaus and zigzag instabilities and for the transition to traveling waves. Neutral stability curves demonstrate the destabilization of stationary planar patterns at low wave numbers to zigzag and traveling modes. Numerical solutions of the model system support the theoretical findings.

Calcium, BOBs, QEDs, microdomains and a cellular decision: control of mitotic cell division in sand dollar blastomeres

The role of Ca2+ in controlling cell processes (e.g. mitosis) presents an enigma in its ubiquity and selectivity. Intracellular free Ca2+ (Ca2+i) is an essential regulator of specific biochemical and physiological aspects of mitosis (e.g. nuclear envelope breakdown (NEB)). Changes in Ca2+i concentrations during mitosis in second cell-cycle sand dollar (Echinaracnius parma) blastomeres were imaged as Ca(2+)-dependent luminescence of the photoprotein aequorin with multi-spectral analytical video microscopy. Photons of this luminescence were seen as bright observable blobs (BOBs). Spatiotemporal patterns of BOBs were followed through one or more cell cycles to detect directly changes in Ca2+i, and were seen to change in a characteristic fashion prior to NEB, the onset of anaphase chromosome movement, and during cytokinesis. These patterns were observed from one cell cycle to the next in a single cell, from cell to cell, and from egg batch to egg batch. In both mitosis and synaptic transmission increases in Ca2+i concentration occurs in discrete, short-lived, highly localized pulses we name quantum emission domains (QEDs) within regions we named microdomains. Signal and statistical optical analyses of spatiotemporal BOB patterns show that many BOBs are linked by constant displacements in space-time (velocity). Linked BOBs are thus nonrandom and are classified as QEDS. Analyses of QED patterns demonstrated that the calcium signals required for NEB are nonrandom, and are evoked by an agent(s) generated proximal to a Ca2+i-QED; models of waves, diffusible agonists and Ca(2+)-activated Ca2+ release do not fit pre-NEB cell data. Spatial and temporal resolution of this multispectral approach significantly exceeds that reported for other methods, and avoids the perturbations associated with many fluorescent Ca2+ reporters that interfere with cells being studied (Ca(2+)-buffering, UV toxicity, etc.). Spatiotemporal patterns of Ca2+i-QED can control so many different processes, i.e. specific frequencies used to control particular processes. Predictive and structured patterns of calcium signals (e.g. a language expressed in Ca2+) may selectively regulate specific Ca(2+)-dependent cellular processes.