An official website of the United States government
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
Bifurcations and pattern formation in a host–parasitoid model with nonlocal effect
In this paper, we analyse Turing instability and bifurcations in a host–parasitoid model with nonlocal effect. For a ordinary differential equation model, we provide some preliminary analysis on Hopf bifurcation. For a reaction–diffusion model with local intraspecific prey competition, we first explore the Turing instability of spatially homogeneous steady states. Next, we show that the model can undergo Hopf bifurcation and Turing–Hopf bifurcation, and find that a pair of spatially nonhomogeneous periodic solutions is stable for a(8,0)-mode Turing–Hopf bifurcationand unstable for a(3,0)-mode Turing–Hopf bifurcation. For a reaction–diffusion model with nonlocal intraspecific prey competition, we study the existence of the Hopf bifurcation, double-Hopf bifurcation, Turing bifurcation, and Turing–Hopf bifurcation successively, and find that a spatially nonhomogeneous quasi-periodic solution is unstable for a(0,1)-mode double-Hopf bifurcation. Our results indicate that the model exhibits complex pattern formations, including transient states, monostability, bistability, and tristability. Finally, numerical simulations are provided to illustrate complex dynamics and verify our theoretical results.
more »
« less
- Award ID(s):
- 2052648
- PAR ID:
- 10520485
- Publisher / Repository:
- Proceedings of the Royal Society of Edinburgh
- Date Published:
- Journal Name:
- Proceedings of the Royal Society of Edinburgh: Section A Mathematics
- ISSN:
- 0308-2105
- Page Range / eLocation ID:
- 1 to 40
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract For multi-scale differential equations (or fast–slow equations), one often encounters problems in which a key system parameter slowly passes through a bifurcation. In this article, we show that a pair of prototypical reaction–diffusion equations in two space dimensions can exhibit delayed Hopf bifurcations. Solutions that approach attracting/stable states before the instantaneous Hopf point stay near these states for long, spatially dependent times after these states have become repelling/unstable. We use the complex Ginzburg–Landau equation and the Brusselator models as prototypes. We show that there exist two-dimensional spatio-temporal buffer surfaces and memory surfaces in the three-dimensional space-time. We derive asymptotic formulas for them for the complex Ginzburg–Landau equation and show numerically that they exist also for the Brusselator model. At each point in the domain, these surfaces determine how long the delay in the loss of stability lasts, that is, to leading order when the spatially dependent onset of the post-Hopf oscillations occurs. Also, the onset of the oscillations in these partial differential equations is a hard onset.more » « less
-
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
-
ABSTRACT Predator‐prey models, such as the Leslie‐Gower model, are essential for understanding population dynamics and stability within ecosystems. These models help explain the balance between species under natural conditions, but the inclusion of factors like the Allee effect and intraspecific competition adds complexity and realism to these interactions, enhancing our ability to predict system behavior under stress. To detect early indicators of population collapse, this study investigates the intricate dynamics of a modified Leslie‐Gower predator‐prey model with both Allee effect and intraspecific competition. We analyze the existence and stability of equilibria, as well as bifurcation phenomena, including saddle‐node bifurcations of codimension 2, Hopf bifurcations of codimension 2, and Bogdanov‐Takens bifurcations of codimension at least 4. Detailed transitions between bifurcation curves–specifically saddle‐node, Hopf, homoclinic, and limit cycle bifurcations–are also examined. We observe a novel transition phenomenon, where a system jumps from saddle‐node bifurcation to homoclinic and limit cycle bifurcations. This suggests that burst oscillations may serve as an early warning of system collapse rather than simply a tipping point. Our findings indicate that moderate levels of intraspecific competition or Allee effect support coexistence of both populations, while excessive levels may destabilize the entire biological system, leading to collapse. These insights offer valuable implications for ecological management and the early detection of risks in population dynamics.more » « less
-
Circular domains frequently appear in mathematical modeling in the fields of ecology, biology and chemistry. In this paper, we investigate the equivariant Hopf bifurcation of partial functional differential equations with Neumann boundary condition on a two-dimensional disk. The properties of these bifurcations at equilibriums are analyzed rigorously by studying the equivariant normal forms. Two reaction–diffusion systems with discrete time delays are selected as numerical examples to verify the theoretical results, in which spatially inhomogeneous periodic solutions including standing waves and rotating waves, and spatially homogeneous periodic solutions are found near the bifurcation points.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/prm.2024.24
