Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMulti-Spike Patterns in the Gierer-Meinhardt System with a Non-Zero Activator Boundary Flux
73
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The structure, linear stability, and dynamics of localized solutions to singularly perturbed reaction-diffusion equations has been the focus of numerous rigorous, asymptotic, and numerical studies in the last few decades. However, with a few exceptions, these studies have often assumed homogeneous boundary conditions. Motivated by the recent focus on the analysis of bulk-surface coupled problems we consider the effect of inhomogeneous Neumann boundary conditions for the activator in the singularly perturbed one-dimensional Gierer-Meinhardt reaction-diffusion system. We show that these boundary conditions necessitate the formation of spikes that concentrate in a boundary layer near the domain boundaries. Using the method of matched asymptotic expansions we construct boundary layer spikes and derive a new class of shifted Nonlocal Eigenvalue Problems for which we rigorously prove partial stability results. Moreover by using a combination of asymptotic, rigorous, and numerical methods we investigate the structure and linear stability of selected one- and two-spike patterns. In particular we find that inhomogeneous Neumann boundary conditions increase both the range of parameter values over which asymmetric two-spike patterns exist and are stable.
rate research
Read More
We analyze a coupled bulk-membrane PDE model in which a scalar linear 2-D bulk diffusion process is coupled through a linear Robin boundary condition to a two-component 1-D reaction-diffusion (RD) system with Gierer-Meinhardt (nonlinear) reaction kinetics defined on the domain boundary. For this coupled model, in the singularly perturbed limit of a long-range inhibition and short-range activation for the membrane-bound species, asymptotic methods are used to analyze the existence of localized steady-state multi-spike membrane-bound patterns, and to derive a nonlocal eigenvalue problem (NLEP) characterizing $mathcal{O}(1)$ time-scale instabilities of these patterns. A central, and novel, feature of this NLEP is that it involves a membrane Greens function that is coupled nonlocally to a bulk Greens function. When the domain is a disk, or in the well-mixed shadow-system limit corresponding to an infinite bulk diffusivity, this Greens function problem is analytically tractable, and as a result we will use a hybrid analytical-numerical approach to determine unstable spectra of this NLEP. This analysis characterizes how the 2-D bulk diffusion process and the bulk-membrane coupling modifies the well-known linear stability properties of steady-state spike patterns for the 1-D Gierer-Meinhardt model in the absence of coupling. In particular, phase diagrams in parameter space for our coupled model characterizing either oscillatory instabilities due to Hopf bifurcations, or competition instabilities due to zero-eigenvalue crossings are constructed. Finally, linear stability predictions from the NLEP analysis are confirmed with full numerical finite-element simulations of the coupled PDE system.
In this paper we consider the existence and stability of multi-spike solutions to the fractional Gierer-Meinhardt model with periodic boundary conditions. In particular we rigorously prove the existence of symmetric and asymmetric two-spike solutions using a Lyapunov-Schmidt reduction. The linear stability of these two-spike solutions is then rigorously analyzed and found to be determined by the eigenvalues of a certain $2times 2$ matrix. Our rigorous results are complemented by formal calculations of $N$-spike solutions using the method of matched asymptotic expansions. In addition, we explicitly consider examples of one- and two-spike solutions for which we numerically calculate their relevant existence and stability thresholds. By considering a one-spike solution we determine that the introduction of fractional diffusion for the activator or inhibitor will respectively destabilize or stabilize a single spike solution with respect to oscillatory instabilities. Furthermore, when considering two-spike solutions we find that the range of parameter values for which asymmetric two-spike solutions exist and for which symmetric two-spike solutions are stable with respect to competition instabilities is expanded with the introduction of fractional inhibitor diffusivity. However our calculations indicate that asymmetric two-spike solutions are always linearly unstable.
Localized spot patterns, where one or more solution components concentrates at certain points in the domain, are a common class of localized pattern for reaction-diffusion systems, and they arise in a wide range of modeling scenarios. In an arbitrary bounded 3-D domain, the existence, linear stability, and slow dynamics of localized multi-spot patterns is analyzed for the well-known singularly perturbed Gierer-Meinhardt (GM) activator-inhibitor system in the limit of a small activator diffusivity $varepsilon^2ll 1$. Our main focus is to classify the different types of multi-spot patterns, and predict their linear stability properties, for different asymptotic ranges of the inhibitor diffusivity $D$. For the range $D={mathcal O}(varepsilon^{-1})gg 1$, although both symmetric and asymmetric quasi-equilibrium spot patterns can be constructed, the asymmetric patterns are shown to be always unstable. On this range of $D$, it is shown that symmetric spot patterns can undergo either competition instabilities or a Hopf bifurcation, leading to spot annihilation or temporal spot amplitude oscillations, respectively. For $D={mathcal O}(1)$, only symmetric spot quasi-equilibria exist and they are linearly stable on ${mathcal O}(1)$ time intervals. On this range, it is shown that the spot locations evolve slowly on an ${mathcal O}(varepsilon^{-3})$ time scale towards their equilibrium locations according to an ODE gradient flow, which is determined by a discrete energy involving the reduced-wave Greens function. The central role of the far-field behavior of a certain core problem, which characterizes the profile of a localized spot, for the construction of quasi-equilibria in the $D={mathcal O}(1)$ and $D={mathcal O}(varepsilon^{-1})$ regimes, and in establishing some of their linear stability properties, is emphasized.
Gap solitons near a band edge of a spatially periodic nonlinear PDE can be formally approximated by solutions of Coupled Mode Equations (CMEs). Here we study this approximation for the case of the 2D Periodic Nonlinear Schr{o}dinger / Gross-Pitaevskii Equation with a non-separable potential of finite contrast. We show that unlike in the case of separable potentials [T. Dohnal, D. Pelinovsky, and G. Schneider, J. Nonlin. Sci. {bf 19}, 95--131 (2009)] the CME derivation has to be carried out in Bloch rather than physical coordinates. Using the Lyapunov-Schmidt reduction we then give a rigorous justification of the CMEs as an asymptotic model for reversible non-degenerate gap solitons and even potentials and provide $H^s$ estimates for this approximation. The results are confirmed by numerical examples including some new families of CMEs and gap solitons absent for separable potentials.
Several mechanisms have been proposed to explain the spontaneous generation of self-organized patterns, hypothesised to play a role in the formation of many of the magnificent patterns observed in Nature. In several cases of interest, the system under scrutiny displays a homogeneous equilibrium, which is destabilized via a symmetry breaking instability which reflects the specificity of the problem being inspected. The Turing instability is among the most celebrated paradigms for pattern formation. In its original form, the diffusion constants of the two mobile species need to be quite different from each other for the instability to develop. Unfortunately, this condition limits the applicability of the theory. To overcome this impediment, and with the ambitious long term goal to eventually reconcile theory and experiments, we here propose an alternative mechanism for promoting the onset of patterns. To this end a multi-species reaction-diffusion system is studied on a discrete, network-like support: the instability is triggered by the non-normality of the embedding network. The non-normal character of the dynamics instigates a short time amplification of the imposed perturbation, thus making the system unstable for a choice of parameters that would yield stability under the conventional scenario. Importantly, non-normal networks are pervasively found, as we shall here briefly review.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
