Complex patterns in reaction-diffusion systems: A tale of two front instabilities

Dispersal-induced resilience to stochastic environmental fluctuations in populations with Allee effect

Many species are unsustainable at small population densities (Allee effect); i.e., below the so-called Allee threshold, the population decreases instead of growing. In a closed local population, environmental fluctuations always lead to extinction. Here, we show how, in spatially extended habitats, dispersal can lead to a sustainable population in a region, provided the amplitude of environmental fluctuations is below an extinction threshold. We have identified two types of sustainable populations: high-density and low-density populations (through a mean-field approximation, valid in the limit of large dispersal length). Our results show that patches where population is high, low, or extinct coexist when the population is close to global extinction (even for homogeneous habitats). The extinction threshold is maximum for characteristic dispersal distances much larger than the spatial scale of synchrony of environmental fluctuations. The extinction threshold increases proportionally to the square root of the dispersal rate and decreases with the Allee threshold. The low-population-density solution can allow understanding of difficulties in recovery after harvesting. This theoretical framework provides a unique approach to address other factors, such as habitat fragmentation or harvesting, impacting population resilience to environmental fluctuations.

Belowground feedbacks as drivers of spatial self-organization and community assembly

Vegetation patterning in water-limited and other resource-limited ecosystems highlights spatial self-organization processes as potentially key drivers of community assembly. These processes provide insight into predictable landscape-level relationships between organisms and their abiotic environment in the form of regular and irregular patterns of biota and resources. However, two aspects have largely been overlooked; the roles played by plant - soil-biota feedbacks and allelopathy in spatial self-organization, and their potential contribution, along with plant-resource feedbacks, to community assembly through spatial self-organization. Here, we expand the drivers of spatial self-organization from a focus on plant-resource feedbacks to include plant - soil-biota feedbacks and allelopathy, and integrate concepts of nonlinear physics and community ecology to generate a new hypothesis. According to this hypothesis, below-ground processes can affect community assemblages through two types of spatial self-organization, global and local. The former occurs simultaneously across whole ecosystems, leading to self-organized patterns of biota, allelochemicals and resources, and niche partitioning. The latter occurs locally in ecotones, and determines ecotone structure and motion, invasion dynamics, and species coexistence. Studies of the two forms of spatial self-organization are important for understanding the organization of plant communities in drier climates which are likely to involve spatial patterning or re-patterning. Such studies are also important for developing new practices of ecosystem management, based on local manipulations at ecotones, to slow invasion dynamics or induce transitions from transitive to intransitive networks of interspecific interactions which increase species diversity.

Formation of localized states in dryland vegetation: Bifurcation structure and stability

We study theoretically the emergence of localized states of vegetation close to the onset of desertification. These states are formed through the locking of vegetation fronts, connecting a uniform vegetation state with a bare soil state, which occurs nearby the Maxwell point of the system. To study these structures we consider a universal model of vegetation dynamics in drylands, which has been obtained as the normal form for different vegetation models. Close to the Maxwell point localized gaps and spots of vegetation exist and undergo collapsed snaking. The presence of gaps strongly suggest that the ecosystem may undergo a recovering process. In contrast, the presence of spots may indicate that the ecosystem is close to desertification.

Continuum Modeling of Discrete Plant Communities: Why Does It Work and Why Is It Advantageous?

Understanding ecosystem response to drier climates calls for modeling the dynamics of dryland plant populations, which are crucial determinants of ecosystem function, as they constitute the basal level of whole food webs. Two modeling approaches are widely used in population dynamics, individual (agent)-based models and continuum partial-differential-equation (PDE) models. The latter are advantageous in lending themselves to powerful methodologies of mathematical analysis, but the question of whether they are suitable to describe small discrete plant populations, as is often found in dryland ecosystems, has remained largely unaddressed. In this paper, we first draw attention to two aspects of plants that distinguish them from most other organisms—high phenotypic plasticity and dispersal of stress-tolerant seeds—and argue in favor of PDE modeling, where the state variables that describe population sizes are not discrete number densities, but rather continuous biomass densities. We then discuss a few examples that demonstrate the utility of PDE models in providing deep insights into landscape-scale behaviors, such as the onset of pattern forming instabilities, multiplicity of stable ecosystem states, regular and irregular, and the possible roles of front instabilities in reversing desertification. We briefly mention a few additional examples, and conclude by outlining the nature of the information we should and should not expect to gain from PDE model studies.

Front Instabilities Can Reverse Desertification

Degradation processes in living systems often take place gradually by front propagation. An important context of such processes is loss of biological productivity in drylands or desertification. Using a dryland-vegetation model, we analyze the stability and dynamics of desertification fronts, identify linear and nonlinear front instabilities, and highlight the significance of these instabilities in inducing self-recovery. The results are based on the derivation and analysis of a universal amplitude equation for pattern-forming living systems for which nonuniform instabilities cannot emerge from the nonviable (zero) state. The results may therefore be applicable to other contexts of animate matter where degradation processes occur by front propagation.

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.

Pattern formation in a reaction-diffusion system of Fitzhugh-Nagumo type before the onset of subcritical Turing bifurcation

We investigate numerically the behavior of a two-component reaction-diffusion system of Fitzhugh-Nagumo type before the onset of subcritical Turing bifurcation in response to local rigid perturbation. In a large region of parameters, the initial perturbation evolves into a localized structure. In a part of that region, closer to the bifurcation line, this structure turns out to be unstable and covers all the available space over the course of time in a process of self-completion. Depending on the parameter values in two-dimensional (2D) space, this process happens either through generation and evolution of new peaks on oscillatory tails of the initial pattern, or through the elongation, deformation, and rupture of initial structure, leading to space-filling nonbranching snakelike patterns. Transient regimes are also possible. Comparison of these results with 1D simulations shows that the prebifurcation region of parameters where the self-completion process is observed is much larger in the 2D case.

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.

Ternary eutectic dendrites: Pattern formation and scaling properties

Extending previous work [Pusztai et al., Phys. Rev. E 87, 032401 (2013)], we have studied the formation of eutectic dendrites in a model ternary system within the framework of the phase-field theory. We have mapped out the domain in which two-phase dendritic structures grow. With increasing pulling velocity, the following sequence of growth morphologies is observed: flat front lamellae → eutectic colonies → eutectic dendrites → dendrites with target pattern → partitionless dendrites → partitionless flat front. We confirm that the two-phase and one-phase dendrites have similar forms and display a similar scaling of the dendrite tip radius with the interface free energy. It is also found that the possible eutectic patterns include the target pattern, and single- and multiarm spirals, of which the thermal fluctuations choose. The most probable number of spiral arms increases with increasing tip radius and with decreasing kinetic anisotropy. Our numerical simulations confirm that in agreement with the assumptions of a recent analysis of two-phase dendrites [Akamatsu et al., Phys. Rev. Lett. 112, 105502 (2014)], the Jackson-Hunt scaling of the eutectic wavelength with pulling velocity is obeyed in the parameter domain explored, and that the natural eutectic wavelength is proportional to the tip radius of the two-phase dendrites. Finally, we find that it is very difficult/virtually impossible to form spiraling two-phase dendrites without anisotropy, an observation that seems to contradict the expectations of Akamatsu et al. Yet, it cannot be excluded that in isotropic systems, two-phase dendrites are rare events difficult to observe in simulations.

Wavy fronts in a hyperbolic FitzHugh-Nagumo system and the effects of cross diffusion

We study a hyperbolic version of the FitzHugh-Nagumo (also known as the Bonhoeffer-van der Pol) reaction-diffusion system. To be able to obtain analytical results, we employ a piecewise linear approximation of the nonlinear kinetic term. The hyperbolic version is compared with the standard parabolic FitzHugh-Nagumo system. We completely describe the dynamics of wavefronts and discuss the properties of the speed equation. The nonequilibrium Ising-Bloch bifurcation of traveling fronts is found to occur in the hyperbolic case as well as in the parabolic system. Waves in the hyperbolic case typically propagate with lower speeds, in absolute value, than waves in the parabolic one. We find the interesting feature that the hyperbolic and parabolic front trajectories coincide in the phase plane for the FitzHugh-Nagumo model with a diagonal diffusion matrix, which is the case of self-diffusion, and differ for the system with cross diffusion.

Interaction of multiarmed spirals in bistable media

We study the interaction of both dense and sparse multiarmed spirals in bistable media modeled by equations of the FitzHugh-Nagumo type. A dense one-armed spiral is characterized by its fixed tip. For dense multiarmed spirals, when the initial distance between tips is less than a critical value, the arms collide, connect, and disconnect continuously as the spirals rotate. The continuous reconstruction between the front and the back drives the tips to corotate along a rough circle and to meander zigzaggedly. The rotation frequency of tip, the frequency of zigzagged displacement, the frequency of spiral, the oscillation frequency of media, and the number of arms satisfy certain relations as long as the control parameters of the model are fixed. When the initial distance between tips is larger than the critical value, the behaviors of individual arms within either dense or sparse multiarmed spirals are identical to that of corresponding one-armed spirals.

Chemical morphogenesis: recent experimental advances in reaction-diffusion system design and control

In his seminal 1952 paper, Alan Turing predicted that diffusion could spontaneously drive an initially uniform solution of reacting chemicals to develop stable spatially periodic concentration patterns. It took nearly 40 years before the first two unquestionable experimental demonstrations of such reaction-diffusion patterns could be made in isothermal single phase reaction systems. The number of these examples stagnated for nearly 20 years. We recently proposed a design method that made their number increase to six in less than 3 years. In this report, we formally justify our original semi-empirical method and support the approach with numerical simulations based on a simple but realistic kinetic model. To retain a number of basic properties of real spatial reactors but keep calculations to a minimal complexity, we introduce a new way to collapse the confined spatial direction of these reactors. Contrary to similar reduced descriptions, we take into account the effect of the geometric size in the confinement direction and the influence of the differences in the diffusion coefficient on exchange rates of species with their feed environment. We experimentally support the method by the observation of stationary patterns in red-ox reactions not based on oxihalogen chemistry. Emphasis is also brought on how one of these new systems can process different initial conditions and memorize them in the form of localized patterns of different geometries.

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.

Self-sustained collective oscillation generated in an array of nonoscillatory cells

Oscillations are ubiquitous phenomena in biological systems. Conventional models of biological periodic oscillations usually invoke interconnecting transcriptional feedback loops. Some specific proteins function as transcription factors, which in turn negatively regulate the expression of the genes that encode these "clock proteins." These loops may lead to rhythmic changes in gene expression in a cell. In the case of multicellular tissue, collective oscillation is often due to the synchronization of these cells, which manifest themselves as autonomous oscillators. In contrast, we propose here a different scenario for the occurrence of collective oscillation in a group of nonoscillatory cells. Neither periodic external stimulation nor pacemaker cells with intrinsically oscillator are included in the present system. By adopting a spatially inhomogeneous active factor, we observe and analyze a coupling-induced oscillation, inherent to the phenomenon of wave propagation due to intracellular communication.

Pattern formation in the ferrocyanide-iodate-sulfite reaction: the control of space scale separation

We revisit the conditions for the development of reaction-diffusion patterns in the ferrocyanide-iodate-sulfite bistable and oscillatory reaction. This hydrogen ion autoactivated reaction is the only example known to produce sustained stationary lamellar patterns and a wealth of other spatio-temporal phenomena including self-replication and localized oscillatory domain of spots, due to repulsive front interactions and to a parity-breaking front bifurcation (nonequilibrium Ising-Bloch bifurcation). We show experimentally that the space scale separation necessary for the observation of stationary patterns is mediated by the presence of low mobility weak acid functional groups. The presence of such groups was overlooked in the original observations made with hydrolyzable polyacrylamide gels. This missing information made the original observations difficult to reproduce and frustrated further experimental exploitation of the fantastic potentialities of this system. Using one-side-fed spatial reactors filled with agarose gel, we can reproduce all the previous pattern observations, in particular the stationary labyrinthine patterns, by introducing, above a critical concentration, well controlled amounts of polyacrylate chains in the gel network. We use two different geometries of spatial reactors (annular and disk shapes) to provide complementary information on the actual three-dimensional character of spatial patterns. We also reinvestigate the role of other feed parameters and show that the system exhibits both a domain of spatial bistability and of large-amplitude pH oscillations associated in a typical cross-shape diagram. The experimental method presented here can be adapted to produce patterns in the large number of oscillatory and bistable reactions, since the iodate-sulfite-ferrocynide reaction is a prototype of these systems.

Front motion and localized states in an asymmetric bistable activator-inhibitor system with saturation

We study the spatiotemporal properties of coherent states (peaks, holes, and fronts) in a bistable activator-inhibitor system that exhibits biochemical saturated autocatalysis, and in which fronts do not preserve spatial parity symmetry. Using the Gierer-Meinhardt prototype model, we find the conditions in which two distinct pinning regions are formed. The first pinning type is known in the context of variational systems while the second is structurally different due to the presence of a heteroclinic bifurcation between two uniform states. The bifurcation also separates the parameter regions of counterpropagating fronts, leading in turn to the growth or contraction of activator domains. These phenomena expand the range of pattern formation theory and its biomedical applications: activator domain retraction suggests potential therapeutic strategies for patterned pathologies, such as cardiovascular calcification.

Wavy fronts and speed bifurcation in excitable systems with cross diffusion

A bifurcation of excitation fronts induced by cross diffusion in two-component bistable reaction-diffusion systems of activator-inhibitor type is discovered. This bifurcation is similar to the nonequilibrium Ising-Bloch bifurcation. A different type of fronts, whose spatial profiles are characterized by oscillating tails, are associated with this bifurcation. These fronts are described using exact analytical solutions of piecewise linear reaction-diffusion equations.

Spatial bistability: a source of complex dynamics. From spatiotemporal reaction-diffusion patterns to chemomechanical structures

We show experimentally and theoretically that reaction systems characterized by a slow induction period followed by a fast evolution to equilibrium can readily generate "spatial bistability" when operated in thin gel reactors diffusively fed from one side. This phenomenon which corresponds to the coexistence of two different stable steady states, not breaking the symmetry of the boundary conditions, can be at the origin of diverse reaction-diffusion instabilities. Using different chemical reactions, we show how stationary pulses, labyrinthine patterns or spatiotemporal oscillations can be generated. Beyond simple reaction-diffusion instabilities, we also demonstrate that the cross coupling of spatial bistability with the size responsiveness of a chemosensitive gel can give rise to autonomous spatiotemporal shape patterns, referred to as chemomechanical structures.

Resonant and nonresonant patterns in forced oscillators

Uniform oscillations in spatially extended systems resonate with temporal periodic forcing within the Arnold tongues of single forced oscillators. The Arnold tongues are wedge-like domains in the parameter space spanned by the forcing amplitude and frequency, within which the oscillator's frequency is locked to a fraction of the forcing frequency. Spatial patterning can modify these domains. We describe here two pattern formation mechanisms affecting frequency locking at half the forcing frequency. The mechanisms are associated with phase-front instabilities and a Turing-like instability of the rest state. Our studies combine experiments on the ruthenium catalyzed light-sensitive Belousov-Zhabotinsky reaction forced by periodic illumination, and numerical and analytical studies of two model systems, the FitzHugh-Nagumo model and the complex Ginzburg-Landau equation, with additional terms describing periodic forcing.

Breaking of translational symmetry of a traveling planar impulse in a two-dimensional two-variable reaction-diffusion model

The stability of a planar impulse in rectangular spatial domains for a two-variable excitable reaction-diffusion system is numerically studied. The dependence of the stability on the size of the domain perpendicular to the direction of the propagation of the impulse is shown. The instability results in asymptotic stable curved impulses or an asymptotic spatiotemporal structure, which is generated similarly to the one-dimensional backfiring phenomenon.

Curvature effects on spiral spectra: generation of point eigenvalues near branch points

It is shown that spiral waves may possess many isolated point eigenvalues that appear near branch points of the linear dispersion relation. These eigenvalues are created by the same mechanism that leads to infinitely many bound states for selfadjoint Schrödinger operators with sufficiently weakly decaying long-range potentials. For spirals, the weak decay of the potential is due to the curvature effects on the profile of the spiral in an intermediate spatial range that separates the spiral core from the far field.

Flow-distributed oscillation, flow-velocity modulation, and resonance

We examine the effects of a periodically varying flow velocity on the standing- and traveling-wave patterns formed by the flow-distributed oscillation mechanism. In the kinematic (or diffusionless) limit, the phase fronts undergo a simple, spatiotemporally periodic longitudinal displacement. On the other hand, when the diffusion is significant, periodic modulation of the velocity can disrupt the wave pattern, giving rise in the downstream region to traveling waves whose frequency is a rational multiple of the velocity perturbation frequency. We observe frequency locking at ratios of 1:1, 2:1, and 3:1, depending on the amplitude and frequency of the velocity modulation. This phenomenon can be viewed as a novel, rather subtle type of resonant forcing.

Self-propelled running droplets on solid substrates driven by chemical reactions

We study chemically driven running droplets on a partially wetting solid substrate by means of coupled evolution equations for the thickness profile of the droplets and the density profile of an adsorbate layer. Two models are introduced corresponding to two qualitatively different types of experiments described in the literature. In both cases an adsorption or desorption reaction underneath the droplets induces a wettability gradient on the substrate and provides the driving force for droplet motion. The difference lies in the behavior of the substrate behind the droplet. In case I the substrate is irreversibly changed whereas in case II it recovers allowing for a periodic droplet movement (as long as the overall system stays far away from equilibrium). Both models allow for a non-saturated and a saturated regime of droplet movement depending on the ratio of the viscous and reactive time scales. In contrast to model I, model II allows for sitting drops at high reaction rate and zero diffusion along the substrate. The transition from running to sitting drops in model II occurs via a super- or subcritical drift-pitchfork bifurcation and may be strongly hysteretic implying a coexistence region of running and sitting drops.

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.

Pattern formation by boundary forcing in convectively unstable, oscillatory media with and without differential transport

Convectively unstable, open reactive flows of oscillatory media, whose phase is fixed or periodically modulated at the inflow boundary, are known to result in stationary and traveling waves, respectively. The latter are implicated in biological segmentation. The boundary-controlled pattern selection by this flow-distributed oscillator (FDO) mechanism has been generalized to include differential flow (DIFI) and differential diffusion (Turing) modes. Our present goal is to clarify the relationships among these mechanisms in the general case where there is differential flow as well as differential diffusion. To do so we analyze the dispersion relation for linear perturbations in the presence of periodic boundary forcing, and show how the solutions are affected by differential transport. We find that the DIFI and FDO modes are closely related and lie in the same frequency range, while the Turing mechanism gives rise to a distinct set of unstable modes in a separate frequency range. Finally, we substantiate the linear analysis by nonlinear simulations and touch upon the issue of competition of spatial modes.

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.

Front bifurcation in a tristable reaction-diffusion system under periodic forcing

A piecewise linear tristable reaction-diffusion equation under external forcing of periodic type is considered. A special feature of the forcing is that the force moves together with the traveling wave. Front velocity equations are obtained analytically using matching procedures for the front solutions. It is noted that there is a restriction in building of null-cline. For each choice of outer branches of null-cline the middle interfacial zone should not exceed some critical value. When this zone is larger the front does not exist. It is found that in the presence of forcing there exists a set of front solutions with different phases (matching point coordinates). The periodic forcing produces a change in the velocity-versus-phase diagram. For a specific choice of wave number, there is a bubble formation which corresponds to additional solutions when the velocity bifurcates to form three fronts.

General theory of nonlinear flow-distributed oscillations

We outline a general theory for the analysis of flow-distributed standing and traveling wave patterns in one-dimensional, open flows of oscillatory chemical media, emphasizing features that are generic to a variety of kinetic models. We draw particular attention to the cases far from a Hopf bifurcation and far from the so-called kinematic or zero-diffusion limit. We introduce a nonlinear formalism for both traveling and stationary waves and show that the wave forms and their amplitudes depend on a single reduced transport parameter that quantifies the departure from the kinematic limit. The nonlinear formalism can be applied to systems with more complex types of bifurcations (canards, period doublings, etc.).

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.

Traveling fronts in an array of coupled symmetric bistable units

Traveling fronts are shown to occur in an array of nearest-neighbor coupled symmetric bistable units. When the local dynamics is given by the Lorenz equations we observe the route: standing-->oscillating-->traveling front, as the coupling is increased. A key step in this route is a gluing bifurcation of two cycles in cylindrical coordinates. When this is mediated by a saddle with real leading eigenvalues, the asymptotic behavior of the front velocity is found straightforwardly. If the saddle is focus-type instead, the front's dynamics may become quite complex, displaying several oscillating and propagating regimes and including (Shil'nikov-type) chaotic front propagation. These results stand as well for other nearest-neighbor coupling schemes and local dynamics.

Modulated structures in electroconvection in nematic liquid crystals

Motivated by experiments in electroconvection in nematic liquid crystals with homeotropic alignment we study the coupled amplitude equations describing the formation of a stationary roll pattern in the presence of a weakly damped mode that breaks isotropy. The equations can be generalized to describe the planarly aligned case if the orienting effect of the boundaries is small, which can be achieved by a destabilizing magnetic field. The slow mode represents the in-plane director at the center of the cell. The simplest uniform states are normal rolls, which may undergo a pitchfork bifurcation to abnormal rolls with a misaligned in-plane director. We present a new class of defect-free solutions with spatial modulations perpendicular to the rolls. In a parameter range where the zigzag instability is not relevant these solutions are stable attractors, as observed in experiments. We also present two-dimensionally modulated states with and without defects which result from the destabilization of the one-dimensionally modulated structures. Finally, for no (or very small) damping, and away from the rotationally symmetric case, we find static chevrons made up of a periodic arrangement of defect chains (or bands of defects) separating homogeneous regions of oblique rolls with very small amplitude. These states may provide a model for a class of poorly understood stationary structures observed in various highly conducting materials ("prechevrons" or "broad domains").

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.

Absolute versus convective instability of spiral waves

Absolute and convective instabilities of spirals are investigated using the continuous and the so-called absolute spectrum. It is shown that the nature of transport, induced by an absolute instability, is determined by spectral data of the asymptotic wave trains. The results are applied to core and far-field breakup of spiral waves in excitable and oscillatory media.

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.