Composite Catalyst Surfaces: Effect of Inert and Active Heterogeneities on Pattern Formation

Pattern localization to a domain edge

The formation of protein patterns inside cells is generically described by reaction-diffusion models. The study of such systems goes back to Turing, who showed how patterns can emerge from a homogenous steady state when two reactive components have different diffusivities (e.g., membrane-bound and cytosolic states). However, in nature, systems typically develop in a heterogeneous environment, where upstream protein patterns affect the formation of protein patterns downstream. Examples for this are the polarization of Cdc42 adjacent to the previous bud site in budding yeast and the formation of an actin-recruiter ring that forms around a PIP3 domain in macropinocytosis. This suggests that previously established protein patterns can serve as a template for downstream proteins and that these downstream proteins can "sense" the edge of the template. A mechanism for how this edge sensing may work remains elusive. Here we demonstrate and analyze a generic and robust edge-sensing mechanism, based on a two-component mass-conserving reaction-diffusion (McRD) model. Our analysis is rooted in a recently developed theoretical framework for McRD systems, termed local equilibria theory. We extend this framework to capture the spatially heterogeneous reaction kinetics due to the template. This enables us to graphically construct the stationary patterns in the phase space of the reaction kinetics. Furthermore, we show that the protein template can trigger a regional mass-redistribution instability near the template edge, leading to the accumulation of protein mass, which eventually results in a stationary peak at the template edge. We show that simple geometric criteria on the reactive nullcline's shape predict when this edge-sensing mechanism is operational. Thus, our results provide guidance for future studies of biological systems and for the design of synthetic pattern forming systems.

Spiral wave initiation in excitable media

Spiral waves represent an important example of dissipative structures observed in many distributed systems in chemistry, biology and physics. By definition, excitable media occupy a stationary resting state in the absence of external perturbations. However, a perturbation exceeding a threshold results in the initiation of an excitation wave propagating through the medium. These waves, in contrast to acoustic and optical ones, disappear at the medium's boundary or after a mutual collision, and the medium returns to the resting state. Nevertheless, an initiation of a rotating spiral wave results in a self-sustained activity. Such activity unexpectedly appearing in cardiac or neuronal tissues usually destroys their dynamics which results in life-threatening diseases. In this context, an understanding of possible scenarios of spiral wave initiation is of great theoretical importance with many practical applications.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 2)'.

Spiral wave initiation in excitable media

Spiral waves represent an important example of dissipative structures observed in many distributed systems in chemistry, biology and physics. By definition, excitable media occupy a stationary resting state in the absence of external perturbations. However, a perturbation exceeding a threshold results in the initiation of an excitation wave propagating through the medium. These waves, in contrast to acoustic and optical ones, disappear at the medium's boundary or after a mutual collision, and the medium returns to the resting state. Nevertheless, an initiation of a rotating spiral wave results in a self-sustained activity. Such activity unexpectedly appearing in cardiac or neuronal tissues usually destroys their dynamics which results in life-threatening diseases. In this context, an understanding of possible scenarios of spiral wave initiation is of great theoretical importance with many practical applications.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 2)'.

Signal Percolation within a Bacterial Community

Signal transmission among cells enables long-range coordination in biological systems. However, the scarcity of quantitative measurements hinders the development of theories that relate signal propagation to cellular heterogeneity and spatial organization. We address this problem in a bacterial community that employs electrochemical cell-to-cell communication. We developed a model based on percolation theory, which describes how signals propagate through a heterogeneous medium. Our model predicts that signal transmission becomes possible when the community is organized near a critical phase transition between a disconnected and a fully connected conduit of signaling cells. By measuring population-level signal transmission with single-cell resolution in wild-type and genetically modified communities, we confirm that the spatial distribution of signaling cells is organized at the predicted phase transition. Our findings suggest that at this critical point, the population-level benefit of signal transmission outweighs the single-cell level cost. The bacterial community thus appears to be organized according to a theoretically predicted spatial heterogeneity that promotes efficient signal transmission.

Achilles' heel of a traveling pulse subject to a local external stimulus

The response of a traveling pulse to a local external stimulus is considered numerically for a modified three-component Oregonator, which is a model system for the photosensitive Belousov-Zhabotinsky (BZ) reaction. The traveling pulse is traced and constantly stimulated, with the distance between the pulse and the stimulus being kept constant. We are interested in the minimal strength of the spatially localized stimulus in order to eliminate the pulse. The use of a stimulus of small width allows us to detect the point in the pulse most sensitive to the external stimulus, referred to as the "Achilles' heel" of the traveling pulse, at which minimal strength of stimulus causes a collapse of the pulse. Our findings are demonstrated experimentally as well with the photosensitive BZ reaction.

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.

Effective medium theory for reaction rates and diffusion coefficients of heterogeneous systems

A simple effective medium theory is derived for spatially heterogeneous nonlinear reaction-diffusion media. Its validity is tested through comparisons with simulations of front and pulse propagation in systems with spatially varying diffusion coefficients and reaction rates. The theory is able to predict wave speeds if the characteristic front width is much larger than the length scale of the heterogeneities. This condition is violated in media with isolated or weakly connected sites. However, the theory nevertheless provides good results in cases where it correctly predicts the percolation threshold of the medium.

Coherent wave patterns sustained by a localized inhomogeneity in an excitable medium

The influence of a localized inhomogeneity (oscillatory or stationary) on spatiotemporal chaotic state in an excitable reaction-diffusion system is investigated. We find that various coherent wave patterns, such as spiral waves (including multiarmed) and target wave patterns are able to be created by the inhomogeneity from the chaotic state. Due to the growth of these coherent wave patterns, the previously existing turbulent waves in the absence of inhomogeneity are suppressed. At last, the whole system is entrained by the coherent wave patterns. Closer investigations indicate that the possible mechanisms underlying the inhomogeneity sustained coherent wave patterns seem quite different for oscillatory and stationary inhomogeneities.

Guiding chemical pulses through geometry: Y junctions

We study computationally and experimentally the propagation of chemical pulses in complex geometries. The reaction of interest, CO oxidation, takes place on single crystal Pt(110) surfaces that are microlithographically patterned; they are also addressable through a focused laser beam, manipulated through galvanometer mirrors, capable of locally altering the crystal temperature and thus affecting pulse propagation. We focus on sudden changes in the domain shape (corners in a Y-junction geometry) that can affect the pulse dynamics; we also show how brief, localized temperature perturbations can be used to control reactive pulse propagation. The computational results are corroborated through experimental studies in which the pulses are visualized using reflection anisotropy microscopy.

Geometry-induced pulse instability in microdesigned catalysts: the effect of boundary curvature

We explore the effect of boundary curvature on the instability of reactive pulses in the catalytic oxidation of on microdesigned Pt catalysts. Using ring-shaped domains of various radii, we find that the pulses disappear (decollate from the inert boundary) at a turning point bifurcation, and we trace this boundary in both physical and geometrical parameter space. These computations corroborate experimental observations of pulse decollation.

Spatial synchronization of regular optical patterns

We investigate an extended, nonlinear optical experiment, exhibiting the spontaneous formation of hexagonal patterns out of a stationary bifurcation. The system is exposed to a two-dimensional spatially periodic forcing, namely, static hexagonal patterns, under variation of their spatial periodicity. Parameters are the strength of the forcing and the distance to pattern forming threshold. The system response is quantitatively characterized with different methods. We observe several locking regimes, where the system is entrained by the forcing. Most of the locking regimes can be related to resonances between the different critical wave numbers and the forcing wave number or its spatial harmonics. One particular locking appears to result from two of these simple resonances in a kind of generalized order m:n synchronization. The width of the locking regimes increases with forcing strength, apparently representing a spatial analog of Arnold tongues.

Global organization of dynamics in oscillatory heterogeneous excitable media

An oscillatory heterogeneous excitable medium undergoes a transition from periodic target patterns to a bursting rhythm driven by the spontaneous initiation and termination of spiral waves as coupling or density is reduced. We illustrate these phenomena in monolayers of chick embryonic heart cells using calcium-sensitive fluorescent dyes. These results are modeled in a heterogeneous cellular automaton in which the neighborhood of interaction and cell density is modified. Parameters that give rise to bursting rhythms are organized in distinct zones in parameter space, leading to a global organization that should be applicable to the dynamics in a large class of excitable media.

Damped oscillation in the NO+CO/Pt(100) reaction system

Considering the influence of the surface defects formed in the surface restructuring phase transition in the CO+NO/Pt(100) reaction system, we propose a lattice gas model to investigate the damped oscillation in the high-temperature oscillatory regime by means of a Monte Carlo simulation. The simulation results show that the persistent oscillation can change into a damped one when the fraction of the defects increases. The production rate of CO2 is near to the maximum value when the oscillation is damped to the end. Furthermore, it is found, in the early stage of the oscillation, that the NO decomposition mainly occurs in the 1 x 1 phase and the hex phase is inactive for the reaction. However, as the reaction proceeds, defects are gradually formed in the 1 x 1 <==> hex phase transition, the hex phase becomes active and dominative for the NO decomposition, and then the oscillation becomes damping. The simulation results give an explanation for some previous experimental phenomena.

Propagation through heterogeneous substrates in simple excitable media models

The interaction of waves and obstacles is simulated by adding heterogeneities to a FitzHugh-Nagumo model and a cellular automata model. The cellular automata model is formulated to account for heterogeneities by modelling the interaction between current sources and current sinks. In both models, wave fronts propagate if the size of the heterogeneities is small, and block if the size of the heterogeneities is large. For intermediate values, wave fronts break up into numerous spiral waves. The theoretical models give insights concerning spiral wave formation in heterogeneous excitable media. (c) 2002 American Institute of Physics.

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.

Catalysis on microstructured bimetallic surfaces

Microstructured bimetallic Pt/Rh and Pt/Ti surfaces have been employed to study the dynamics of catalytic NO reduction and the O(2)+H(2) reaction at low pressure (p<10(-3) mbar). Photoelectron emission microscopy and scanning photoelectron microscopy were used as spatially resolved in situ methods to image the local work function changes and to identify chemical changes in the substrate and in the adsorbate layer. It is shown that diffusional coupling leads to dynamic effects which are dependent on the macroscopic size (&mgr;m range). With alkali metals on the surface, stationary patterns form whose mechanism of formation has been studied in detail. (c) 2002 American Institute of Physics.

Front initiation on microdesigned composite catalysts

We first briefly review the subject of spatiotemporal pattern formation on microdesigned composite catalysts. One of the most significant interaction mechanisms between different reacting domains (consisting of different metal catalysts such as Pt and Rh, coupled through surface diffusion) is the initiation of reaction fronts at the interface between them. We then explore in some detail the effect of two-dimensional composite geometry on this basic building block of composite catalyst dynamics. (c) 2002 American Institute of Physics.

Formation of two-dimensional concentration pulses on microdesigned composite catalyst surfaces

We study the effect of microdesigned composite geometries on pattern formation during the catalytic oxidation of CO on Pt-Ti, Pt-Rh, and Pt-Pd composite catalysts. In particular, we find experimentally (and rationalize through modeling) that adsorbate surface transport through the second (active) component drastically affects the shapes and interactions of concentration patterns (traveling pulses) observed on pure Pt.

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."

Two forms of spiral-wave reentry in an ionic model of ischemic ventricular myocardium

It is well known that there is considerable spatial inhomogeneity in the electrical properties of heart muscle, and that the many interventions that increase this initial degree of inhomogeneity all make it easier to induce certain cardiac arrhythmias. We consider here the specific example of myocardial ischemia, which greatly increases the electrical heterogeneity of ventricular tissue, and often triggers life-threatening cardiac arrhythmias such as ventricular tachycardia and ventricular fibrillation. There is growing evidence that spiral-wave activity underlies these reentrant arrhythmias. We thus investigate whether spiral waves might be induced in a realistic model of inhomogeneous ventricular myocardium. We first modify the Luo and Rudy [Circ. Res. 68, 1501-1526 (1991)] ionic model of cardiac ventricular muscle so as to obtain maintained spiral-wave activity in a two-dimensional homogeneous sheet of ventricular muscle. Regional ischemia is simulated by raising the external potassium concentration ([K(+)](o)) from its nominal value of 5.4 mM in a subsection of the sheet, thus creating a localized inhomogeneity. Spiral-wave activity is induced using a pacing protocol in which the pacing frequency is gradually increased. When [K(+)](o) is sufficiently high in the abnormal area (e.g., 20 mM), there is complete block of propagation of the action potential into that area, resulting in a free end or wave break as the activation wave front encounters the abnormal area. As pacing continues, the free end of the activation wave front traveling in the normal area increasingly separates or detaches from the border between normal and abnormal tissue, eventually resulting in the formation of a maintained spiral wave, whose core lies entirely within an area of normal tissue lying outside of the abnormal area ("type I" spiral wave). At lower [K(+)](o) (e.g., 10.5 mM) in the abnormal area, there is no longer complete block of propagation into the abnormal area; instead, there is partial entrance block into the abnormal area, as well as exit block out of that area. In this case, a different kind of spiral wave (transient "type II" spiral wave) can be evoked, whose induction involves retrograde propagation of the action potential through the abnormal area. The number of turns made by the type II spiral wave depends on several factors, including the level of [K(+)](o) within the abnormal area and its physical size. If the pacing protocol is changed by adding two additional stimuli, a type I spiral wave is instead produced at [K(+)](o)=10.5 mM. When pacing is continued beyond this point, apparently aperiodic multiple spiral-wave activity is seen during pacing. We discuss the relevance of our results for arrythmogenesis in both the ischemic and nonischemic heart. (c) 1998 American Institute of Physics.