skip to main content

Explore Research Products in the NSF-PAR

Title: Localized patterns and semi-strong interaction, a unifying framework for reaction–diffusion systems
Abstract Systems of activator–inhibitor reaction–diffusion equations posed on an infinite line are studied using a variety of analytical and numerical methods. A canonical form is considered, which contains all known models with simple cubic autocatalytic nonlinearity and arbitrary constant and linear kinetics. Restricting attention to models that have a unique homogeneous equilibrium, this class includes the classical Schnakenberg and Brusselator models, as well as other systems proposed in the literature to model morphogenesis. Such models are known to feature Turing instability, when activator diffuses more slowly than inhibitor, leading to stable spatially periodic patterns. Conversely in the limit of small feed rates, semi-strong interaction asymptotic analysis shows existence of isolated spike-like patterns. This paper describes the broad bifurcation structures that connect these two regimes. A certain universal two-parameter state diagram is revealed in which the Turing bifurcation becomes sub-critical, leading to the onset of homoclinic snaking. This regime then morphs into the spike regime, with the outer-fold being predicted by the semi-strong asymptotics. A rescaling of parameters and field concentrations shows how this state diagram can be studied independently of the diffusion rates. Temporal dynamics is found to strongly depend on the diffusion ratio though. A Hopf bifurcation occurs along the branch of stable spikes, which is subcritical for small diffusion ratio, leading to collapse to the homogeneous state. As the diffusion ratio increases, this bifurcation typically becomes supercritical and interacts with the homoclinic snaking and also with a supercritical homogeneous Hopf bifurcation, leading to complex spatio-temporal dynamics. The details are worked out for a number of different models that fit the theory using a mixture of weakly nonlinear analysis, semi-strong asymptotics and different numerical continuation algorithms.  more » « less
Award ID(s):
1908891
NSF-PAR ID:
10332046
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
IMA Journal of Applied Mathematics
Volume:
86
Issue:
5
ISSN:
0272-4960
Page Range / eLocation ID:
1031 to 1065
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ( , IMA Journal of Applied Mathematics)
    null (Ed.)
    Abstract An activator–inhibitor–substrate model of side branching used in the context of pulmonary vascular and lung development is considered on the supposition that spatially localized concentrations of the activator trigger local side branching. The model consists of four coupled reaction–diffusion equations, and its steady localized solutions therefore obey an eight-dimensional spatial dynamical system in one spatial dimension (1D). Stationary localized structures within the model are found to be associated with a subcritical Turing instability and organized within a distinct type of foliated snaking bifurcation structure. This behavior is in turn associated with the presence of an exchange point in parameter space at which the complex leading spatial eigenvalues of the uniform concentration state are overtaken by a pair of real eigenvalues; this point plays the role of a Belyakov–Devaney point in this system. The primary foliated snaking structure consists of periodic spike or peak trains with $N$ identical equidistant peaks, $N=1,2,\dots \,$, together with cross-links consisting of nonidentical, nonequidistant peaks. The structure is complicated by a multitude of multipulse states, some of which are also computed, and spans the parameter range from the primary Turing bifurcation all the way to the fold of the $N=1$ state. These states form a complex template from which localized physical structures develop in the transverse direction in 2D. 
    more » « less
  2. ; ( , IMA Journal of Applied Mathematics)
    Abstract In this work, we revisit some general results on the dynamics of circular fronts between homogeneous states and the formation of localized structures in two dimensions (2D). We show how the bifurcation diagram of axisymmetric structures localized in radius fits within the framework of collapsed homoclinic snaking. In 2D, owing to curvature effects, the collapse of the snaking structure follows a different scaling that is determined by the so-called nucleation radius. Moreover, in the case of fronts between two symmetry-related states, the precise point in parameter space to which radial snaking collapses is not a ‘Maxwell’ point but is determined by the curvature-driven dynamics only. In this case, the snaking collapses to a ‘zero surface tension’ point. Near this point, the breaking of symmetry between the homogeneous states tilts the snaking diagram. A different scaling law is found for the collapse of the snaking curve in each case. Curvature effects on axisymmetric localized states with internal structure are also discussed, as are cellular structures separated from a homogeneous state by a circular front. While some of these results are well understood in terms of curvature-driven dynamics and front interactions, a proper mathematical description in terms of homoclinic trajectories in a radial spatial dynamics description is lacking. 
    more » « less
  3. ; ; ( , IMA Journal of Applied Mathematics)
    Abstract In this article, the recently discovered phenomenon of delayed Hopf bifurcations (DHB) in reaction–diffusion partial differential equations (PDEs) is analysed in the cubic Complex Ginzburg–Landau equation, as an equation in its own right, with a slowly varying parameter. We begin by using the classical asymptotic methods of stationary phase and steepest descents on the linearized PDE to show that solutions, which have approached the attracting quasi-steady state (QSS) before the Hopf bifurcation remain near that state for long times after the instantaneous Hopf bifurcation and the QSS has become repelling. In the complex time plane, the phase function of the linearized PDE has a saddle point, and the Stokes and anti-Stokes lines are central to the asymptotics. The non-linear terms are treated by applying an iterative method to the mild form of the PDE given by perturbations about the linear particular solution. This tracks the closeness of solutions near the attracting and repelling QSS in the full, non-linear PDE. Next, we show that beyond a key Stokes line through the saddle there is a curve in the space-time plane along which the particular solution of the linear PDE ceases to be exponentially small, causing the solution of the non-linear PDE to diverge from the repelling QSS and exhibit large-amplitude oscillations. This curve is called the space–time buffer curve. The homogeneous solution also stops being exponentially small in a spatially dependent manner, as determined also by the initial data and time. Hence, a competition arises between these two solutions, as to which one ceases to be exponentially small first, and this competition governs spatial dependence of the DHB. We find four different cases of DHB, depending on the outcomes of the competition, and we quantify to leading order how these depend on the main system parameters, including the Hopf frequency, initial time, initial data, source terms, and diffusivity. Examples are presented for each case, with source terms that are a uni-modal function, a smooth step function, a spatially periodic function and an algebraically growing function. Also, rich spatio-temporal dynamics are observed in the post-DHB oscillations. Finally, it is shown that large-amplitude source terms can be designed so that solutions spend substantially longer times near the repelling QSS, and hence, region-specific control over the delayed onset of oscillations can be achieved. 
    more » « less
  4. ; ; ; ( , Studies in Applied Mathematics)
    Abstract

    We investigate the effect of nonlocal intraspecific prey competition on the spatiotemporal dynamics of a Holling–Tanner predator–prey model with diffusion. We first establish the criteria for Hopf, Turing, double‐Hopf, and Turing–Hopf bifurcations, and determine the stable and unstable regions of the positive equilibrium. For Turing–Hopf bifurcation, by analyzing the normal form truncated to the third order, we derive that, with strong nonlocal interaction, the system exhibits the tristable phenomena, that is, the coexistence of a stable spatially nonhomogeneous periodic orbit and two nonconstant stable steady states, as well as the existence of periodic orbits with two spatial wave frequencies induced by the nonlocal interaction. The main analytical difficulty arises from the nonlocal interaction that prevents the direct application of formulas for the coefficients of the normal form. Biologically, the emerging spatiotemporal patterns suggest that the global intraspecific competition can promote the coexistence of the prey and predator by allowing the prey maintain a critical total population size, which may provide an alternative approach in explaining the group formation of some prey species under the risk of predation. Some coexistence patterns are destabilized by introducing another prey species with local intraspecific competition, leading to the coexistence of two preys and one predator, namely, the prey with nonlocal interaction is concentrated at a single spatial location, and the other prey is distributed uniformly in the rest of the habitat. Accordingly, the predator is forced to change its behavior as well.

     
    more » « less
  5. ; ; ; ( , IMA Journal of Applied Mathematics)
    Abstract Localized coherent structures can form in externally driven dispersive optical cavities with a Kerr-type non-linearity. Such systems are described by the Lugiato–Lefever (LL) equation, which supports a large variety of dynamical states. Here, we review our current knowledge of the formation, stability and bifurcation structure of localized structures in the one-dimensional LL equation. We do so by focusing on two main regimes of operation: anomalous and normal second-order dispersion. In the anomalous regime, localized patterns are organized in a homoclinic snaking scenario, which is eventually destroyed, leading to a foliated snaking bifurcation structure. In the normal regime, localized structures undergo a different type of bifurcation structure, known as collapsed snaking. The effects of third-order dispersion and various dynamical regimes are also described. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email