Existence and stability analysis of spiky solutions for the Gierer-Meinhardt system with large reaction rates

We study the Gierer-Meinhardt system in one dimension in the limit of large reaction rates. First we construct three types of solutions: (i) an interior spike; (ii) a boundary spike and (iii) two boundary spikes. Second we prove results on their stability. It is found that an interior spike is always unstable; a boundary spike is always stable. The two boundary spike configuration can be either stable or unstable, depending on the parameters. We fully classify the stability in this case. We characterise the destabilizing eigenfunctions in all cases. Numerical simulations are shown which are in full agreement with the analytical results

A Tale of Two Distributions: From Few To Many Vortices In Quasi-Two-Dimensional Bose-Einstein Condensates

Motivated by the recent successes of particle models in capturing the precession and interactions of vortex structures in quasi-two-dimensional Bose-Einstein condensates, we revisit the relevant systems of ordinary differential equations. We consider the number of vortices $N$ as a parameter and explore the prototypical configurations ("ground states") that arise in the case of few or many vortices. In the case of few vortices, we modify the classical result of Havelock [Phil. Mag. ${\bf 11}$, 617 (1931)] illustrating that vortex polygons in the form of a ring are unstable for $N \geq7$. Additionally, we reconcile this modification with the recent identification of symmetry breaking bifurcations for the cases of $N=2,\dots,5$. We also briefly discuss the case of a ring of vortices surrounding a central vortex (so-called $N+1$ configuration). We finally examine the opposite limit of large $N$ and illustrate how a coarse-graining, continuum approach enables the accurate identification of the radial distribution of vortices in that limit.Comment: 15 pages, 2 figure

Swarm dynamics and equilibria for a nonlocal aggregation model

We consider the aggregation equation ρt − ∇ · (ρ∇K ∗ ρ) = 0 in Rn, where the interaction potential K models short-range repulsion and long-range attraction. We study a family of interaction potentials for which the equilibria are of finite density and compact support. We show global well-posedness of solutions and investigate analytically and numerically the equilibria and their global stability. In particular, we consider a potential for which the corresponding equilibrium solutions are of uniform density inside a ball of Rn and zero outside. For such a potential, various explicit calculations can be carried out in detail. In one dimension we fully solve the temporal dynamics, and in two or higher dimensions we show the global stability of this steady state within the class of radially symmetric solutions. Finally, we solve the following restricted inverse problem: given a radially symmetric density ρ ̄ that is zero outside some ball of radius R and is polynomial inside the ball, construct an interaction potential K for which ρ̄ is the steady state solution of the corresponding aggregation equation. Throughout the paper, numerical simulations are used to motivate and validate the analytical results

Weakly Nonlinear Analysis of Vortex Formation in a Dissipative Variant of the Gross-Pitaevskii Equation

For a dissipative variant of the two-dimensional Gross-Pitaevskii equation with a parabolic trap under rotation, we study a symmetry breaking process that leads to the formation of vortices. The first symmetry breaking leads to the formation of many small vortices distributed uniformly near the Thomas-Fermi radius. The instability occurs as a result of a linear instability of a vortex-free steady state as the rotation is increased above a critical threshold. We focus on the second subsequent symmetry breaking, which occurs in the weakly nonlinear regime. At slightly above threshold, we derive a one dimensional amplitude equation that describes the slow evolution of the envelope of the initial instability. We show that the mechanism responsible for initiating vortex formation is a modulational instability of the amplitude equation. We also illustrate the role of dissipation in the symmetry breaking process. All analyses are confirmed by detailed numerical computations

Mesa-type patterns in the one-dimensional Brusselator and their stability

The Brusselator is a generic reaction-diffusion model for a tri-molecular chemical reaction. We consider the case when the input and output reactions are slow. In this limit, we show the existence of $K$-periodic, spatially bi-stable structures, \emph{mesas}, and study their stability. Using singular perturbation techniques, we find a threshold for the stability of $K$ mesas. This threshold occurs in the regime where the exponentially small tails of the localized structures start to interact. By comparing our results with Turing analysis, we show that in the generic case, a Turing instability is followed by a slow coarsening process whereby logarithmically many mesas are annihilated before the system reaches a steady equilibrium state. We also study a ``breather''-type instability of a mesa, which occurs due to a Hopf bifurcation. Full numerical simulations are shown to confirm the analytical results.Comment: to appear, Physica

Mode transitions in a model reaction-diffusion system driven by domain growth and noise

Pattern formation in many biological systems takes place during growth of the underlying domain. We study a specific example of a reaction–diffusion (Turing) model in which peak splitting, driven by domain growth, generates a sequence of patterns. We have previously shown that the pattern sequences which are presented when the domain growth rate is sufficiently rapid exhibit a mode-doubling phenomenon. Such pattern sequences afford reliable selection of certain final patterns, thus addressing the robustness problem inherent of the Turing mechanism. At slower domain growth rates this regular mode doubling breaks down in the presence of small perturbations to the dynamics. In this paper we examine the breaking down of the mode doubling sequence and consider the implications of this behaviour in increasing the range of reliably selectable final patterns

Existence and Stability of a Spike in the Central Component for a Consumer Chain Model

We study a three-component consumer chain model which is based on Schnakenberg type kinetics. In this model there is one consumer feeding on the producer and a second consumer feeding on the first consumer. This means that the first consumer (central component) plays a hybrid role: it acts both as consumer and producer. The model is an extension of the Schnakenberg model suggested in \cite{gm,schn1} for which there is only one producer and one consumer. It is assumed that both the producer and second consumer diffuse much faster than the central component. We construct single spike solutions on an interval for which the profile of the first consumer is that of a spike. The profiles of the producer and the second consumer only vary on a much larger spatial scale due to faster diffusion of these components. It is shown that there exist two different single spike solutions if the feed rates are small enough: a large-amplitude and a small-amplitude spike. We study the stability properties of these solutions in terms of the system parameters. We use a rigorous analysis for the linearized operator around single spike solutions based on nonlocal eigenvalue problems. The following result is established: If the time-relaxation constants for both producer and second consumer vanish, the large-amplitude spike solution is stable and the small-amplitude spike solution is unstable. We also derive results on the stability of solutions when these two time-relaxation constants are small. We show a new effect: if the time-relaxation constant of the second consumer is very small, the large-amplitude spike solution becomes unstable. To the best of our knowledge this phenomenon has not been observed before for the stability of spike patterns. It seems that this behavior is not possible for two-component reaction-diffusion systems but that at least three components are required. Our main motivation to study this system is mathematical since the novel interaction of a spike in the central component with two other components results in new types of conditions for the existence and stability of a spike. This model is realistic if several assumptions are made: (i) cooperation of consumers is prevalent in the system, (ii) the producer and the second consumer diffuse much faster than the first consumer, and (iii) there is practically an unlimited pool of producer. The first assumption has been proven to be correct in many types of consumer groups or populations, the second assumption occurs if the central component has a much smaller mobility than the other two, the third assumption is realistic if the consumers do not feel the impact of the limited amount of producer due to its large quantity. This chain model plays a role in population biology, where consumer and producer are often called predator and prey. This system can also be used as a model for a sequence of irreversible autocatalytic reactions in a container which is in contact with a well-stirred reservoir

Parametric pattern selection in a reaction-diffusion model

We compare spot patterns generated by Turing mechanisms with those generated by replication cascades, in a model one-dimensional reaction-diffusion system. We determine the stability region of spot solutions in parameter space as a function of a natural control parameter (feed-rate) where degenerate patterns with different numbers of spots coexist for a fixed feed-rate. While it is possible to generate identical patterns via both mechanisms, we show that replication cascades lead to a wider choice of pattern profiles that can be selected through a tuning of the feed-rate, exploiting hysteresis and directionality effects of the different pattern pathways