Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMulti-Spike Solutions to the Fractional Gierer-Meinhardt System in a One-Dimensional Domain
85
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
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.
It was recently demonstrated that two-dimensional Townes solitons (TSs) in two-component systems with cubic self-focusing, which are normally made unstable by the critical collapse, can be stabilized by linear spin-orbit coupling (SOC), in Bose-Einstein condensates and optics alike. We demonstrate that one-dimensional TSs, realized as optical spatial solitons in a planar dual-core waveguide with dominant quintic self-focusing, may be stabilized by SOC-like terms emulated by obliquity of the coupling between cores of the waveguide. Thus, SOC offers a universal mechanism for the stabilization of the TSs. A combination of systematic numerical considerations and analytical approximations identifies a vast stability area for skew-symmetric solitons in the systems main (semi-infinite) and annex (finite) bandgaps. Tilted (moving) solitons are unstable, spontaneously evolving into robust breathers. For broad solitons, diffraction, represented by second derivatives in the system, may be neglected, leading to a simplified model with a finite bandgap. It is populated by skew-antisymmetric gap solitons, which are nearly stable close to the gaps bottom.
A new class of solutions of three-dimensional equations from the Boussinesq paradigm are considered. The corresponding profiles are not localized functions in the sense of the integrability of the square over an infinite domain. For the new type of solutions, the gradient and the Hessian/Laplacian are square integrable. In the linear limiting case, analytical expressions for the profiles of the pseudolocalized solutions are found. The nonlinear case is treated numerically with a special approximation of the differential operators with spherical symmetry that allows for automatic acknowledgement of the behavioral conditions at the origin of the coordinate system. The asymptotic boundary conditions stem from the $1/r$ behavior at infinity of the pseudolocalized profile. A special approximation is devised that allows us to obtain the proper behavior for much smaller computational box. The pseudolocalized solutions are obtained for both quadratic and cubic nonlinearity.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
