This paper studies the spatio-temporal dynamics of a diffusive plant-sulphide model with toxicity delay. More specifically, the effects of discrete delay and distributed delay on the dynamics are explored, respectively. The deep analysis of eigenvalues indicates that both diffusion and delay can induce Hopf bifurcations. The normal form theory is used to set up an exact formula that determines the properties of Hopf bifurcation in a diffusive plant-sulphide model. A sufficiently small discrete delay does not affect the stability and a sufficiently large discrete delay destabilizes the system. Nonetheless, a sufficiently small or large distributed delay does not affect the stability. Both delays cause instability by inducing Hopf bifurcation rather than Turing bifurcation.

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 bifurcation and 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.

In this paper, we consider a general single population model with delay and patch structure, which could model the population loss during the dispersal. It is shown that the model admits a unique positive equilibrium when the dispersal rate is smaller than a critical value. The stability of the positive equilibrium and associated Hopf bifurcation are investigated when the dispersal rate is small or near the critical value. Moreover, we show the effect of network topology on Hopf bifurcation values for a delayed logistic population model.

We study a fractional-order delayed predator-prey model with Holling–Tanner-type functional response. Mainly, by choosing the delay time $\tau $ as the bifurcation parameter, we show that Hopf bifurcation can occur as the delay time $\tau $ passes some critical values. The local stability of a positive equilibrium and the existence of the Hopf bifurcations are established, and numerical simulations for justifying the theoretical analysis are also presented.

We analyse oscillatory instabilities for a coupled partial-ordinary differential equation (PDE-ODE) system modelling the communication between localised spatially segregated dynamically active signalling compartments that are coupled through a passive extracellular bulk diffusion field in a bounded 2D domain. Each signalling compartment is assumed to secrete a chemical into the extracellular medium (bulk region), and it can also sense the concentration of this chemical in the region around its boundary. This feedback from the bulk region, resulting from the entire collection of cells, in turn modifies the intracellular dynamics within each cell. In the limit where the signalling compartments are circular discs with a small common radius ɛ ≪ 1 and where the bulk diffusivity is asymptotically large, a matched asymptotic analysis is used to reduce the dimensionless PDE–ODE system into a nonlinear ODE system with global coupling. For Sel’kov reaction kinetics, this ODE system for the intracellular dynamics and the spatial average of the bulk diffusion field are then used to investigate oscillatory instabilities in the dynamics of the cells that are triggered due to the global coupling. In particular, numerical bifurcation software on the ODEs is used to study the overall effect of coupling defective cells (cells that behave differently from the remaining cells) to a group of identical cells. Moreover, when the number of cells is large, the Kuramoto order parameter is computed to predict the degree of phase synchronisation of the intracellular dynamics. Quorum sensing behaviour, characterised by the onset of collective behaviour in the intracellular dynamics as the number of cells increases above a threshold, is also studied. Our analysis shows that the cell population density plays a dual role of triggering and then quenching synchronous oscillations in the intracellular dynamics.

In aquatic microbial systems, high-magnitude variations in abundance, such as sudden blooms alternating with comparatively long periods of very low abundance (“apparent disappearance”), are relatively common. We suggest that in order for this to occur, such variations in abundance in microbial systems and, in particular, the apparent disappearance of species do not require seasonal or periodic forcing of any kind or external factors of any other nature. Instead, such variations can be caused by internal factors and, in particular, by bacteria–phage interaction. Specifically, we suggest that the variations in abundance and the apparent disappearance phenomenon can be a result of phage infection and the lysis of infected bacteria. To illustrate this idea, we consider a reasonably simple mathematical model of bacteria–phage interaction based on the model suggested by Beretta and Kuang, which assumes neither periodic forcing nor action of other external factors. The model admits a loss of stability via Andronov–Hopf bifurcation and exhibits dynamics which explains the phenomenon. These properties of the model are especially distinctive for spatially nonhomogeneous biosystems as well as biosystems with some sort of cooperation or community effects.

In 1958 started the study of the families of algebraic limit cycles in the class of planar quadratic polynomial differential systems. In the present we known one family of algebraic limit cycles of degree 2 and four families of algebraic limit cycles of degree 4, and that there are no limit cycles of degree 3. All the families of algebraic limit cycles of degree 2 and 4 are known, this is not the case for the families of degree higher than 4. We also know that there exist two families of algebraic limit cycles of degree 5 and one family of degree 6, but we do not know if these families are all the families of degree 5 and 6. Until today it is an open problem to know if there are algebraic limit cycles of degree higher than 6 inside the class of quadratic polynomial differential systems. Here we investigate the birth and death of all the known families of algebraic limit cycles of quadratic polynomial differential systems.

The Hopf bifurcation from spike solutions for the classical Gierer–Meinhardt system in a onedimensional interval is considered. The existence of time-periodic solution near the Hopf bifurcation parameter for a boundary spike is rigorously proved by the classical Crandall–Rabinowitz theory. The criteria for the stability of the limit cycle are determined, and it is shown that the limit cycle is unstable.

In a singularly perturbed limit, we analyse the existence and linear stability of steady-state hotspot solutions for an extension of the 1-D three-component reaction-diffusion (RD) system formulated and studied numerically in Jones et. al. [Math. Models. Meth. Appl. Sci., 20, Suppl., (2010)], which models urban crime with police intervention. In our extended RD model, the field variables are the attractiveness field for burglary, the criminal population density and the police population density. Our model includes a scalar parameter that determines the strength of the police drift towards maxima of the attractiveness field. For a special choice of this parameter, we recover the ‘cops-on-the-dots’ policing strategy of Jones et. al., where the police mimic the drift of the criminals towards maxima of the attractiveness field. For our extended model, the method of matched asymptotic expansions is used to construct 1-D steady-state hotspot patterns as well as to derive nonlocal eigenvalue problems (NLEPs) that characterise the linear stability of these hotspot steady states to ${\cal O}$(1) timescale instabilities. For a cops-on-the-dots policing strategy, we prove that a multi-hotspot steady state is linearly stable to synchronous perturbations of the hotspot amplitudes. Alternatively, for asynchronous perturbations of the hotspot amplitudes, a hybrid analytical–numerical method is used to construct linear stability phase diagrams in the police vs. criminal diffusivity parameter space. In one particular region of these phase diagrams, the hotspot steady states are shown to be unstable to asynchronous oscillatory instabilities in the hotspot amplitudes that arise from a Hopf bifurcation. Within the context of our model, this provides a parameter range where the effect of a cops-on-the-dots policing strategy is to only displace crime temporally between neighbouring spatial regions. Our hybrid approach to study the NLEPs combines rigorous spectral results with a numerical parameterisation of any Hopf bifurcation threshold. For the cops-on-the-dots policing strategy, our linear stability predictions for steady-state hotspot patterns are confirmed from full numerical PDE simulations of the three-component RD system.

In a recent paper by Cantrell et al. [9], two-component KPP systems with competition of Lotka–Volterra type were analyzed and their long-time behaviour largely settled. In particular, the authors established that any constant positive steady state, if unique, is necessarily globally attractive. In the present paper, we give an explicit and biologically very natural example of oscillatory three-component system. Using elementary techniques or pre-established theorems, we show that it has a unique constant positive steady state with two-dimensional unstable manifold, a stable limit cycle, a predator–prey structure near the steady state, periodic wave trains and point-to-periodic rapid travelling waves. Numerically, we also show the existence of pulsating fronts and propagating terraces.

Dynamics and organization of laser-generated three-dimensional (3D) Richtmyer–Meshkov (RMI) and Rayleigh–Taylor instabilities (RMI and RTI) on metal target in the semiconfined configuration are different in the central region (CR) (Lugomer, 2016), near central region (NCR) (Lugomer, 2017) and the near periphery region (NPR) of the Gaussian-like spot. The RMI/RTI in the NPR evolve from the shock and series of reshocks associated with lateral expansion and increase of the vapor density, decrease of the Atwood number and momentum transfer. Scanning electron micrographs show irregular (chaotic) web of the base-plane walls, and mushroom spikes on its nodal points with disturbed two-dimensional (2D) lattice organization. Lattice disturbance is caused by the incoherent wavy motion of background fluid due to fast reshocks, which after series of reflections change their strength and direction. Reconstruction of the disturbed lattice reveals rectangular lattice of mushroom spikes with p2mm symmetry. The splitting (bifurcation) of mushroom diameter distribution on the large and small mushroom spikes increases with radial distance from the center of Gaussian-like spot. Dynamics of their evolution is represented by the orbits or stable periods in 2D phase space. The constant mushroom diameter – stable circulation or the stable periodic orbits – are the limit cycles between the unstable spiral orbits. Those with increasing periods represent supercritical Hopf bifurcation, while those leading to decrease and disappearance represent subcritical Hopf bifurcation. The empirical models of RMI, although predict dependence of the growth rate on radial distance (distance the reshocks travel to reach the interface), show many limitations. More appropriate interpretation of the simultaneous growth and lattice organization of small and large spikes give the fundamental model based on the interference of the perturbation modes depending on their amplitude, relative phase, and the symmetry. The late-time instability in the base-plane evolves into line solitons, vortex filaments and wave–vortex structures with chaotic rather than stochastic features.

A nonlinear mathematical model for innovation diffusion is proposed. The system of ordinary differential equations incorporates variable external influences (the cumulative density of marketing efforts), variable internal influences (the cumulative density of word of mouth) and a logistically growing human population (the variable potential consumers). The change in population density is due to various demographic processes such as intrinsic growth rate, emigration, death rate etc. Thus the problem involves two dynamic variables viz. a non-adopter population density and an adopter population density. The model is analysed qualitatively using the stability theory of differential equations, with the help of the corresponding characteristic equation of the system. The interior equilibrium point can be stable for all time delays to a critical value, beyond which the system becomes unstable and a Hopf bifurcation occurs at a second critical value. Employing normal form theory and a centre manifold theorem applicable to functional differential equations, we derive some explicit formulas determining the stability, the direction and other properties of the bifurcating periodic solutions. Our numerical simulations show that the system behaviour can become extremely complicated as the time delay increases, with a stable interior equilibrium point leading to a limit cycle with one local maximum and minimum per cycle (Hopf bifurcation), then limit cycles with more local maxima and minima per cycle, and finally chaotic solutions.

For certain singularly perturbed two-component reaction–diffusion systems, the bifurcation diagram of steady-state spike solutions is characterized by a saddle-node behaviour in terms of some parameter in the system. For some such systems, such as the Gray–Scott model, a spike self-replication behaviour is observed as the parameter varies across the saddle-node point. We demonstrate and analyse a qualitatively new type of transition as a parameter is slowly decreased below the saddle node value, which is characterized by a finite-time blow-up of the spike solution. More specifically, we use a blend of asymptotic analysis, linear stability theory, and full numerical computations to analyse a wide variety of dynamical instabilities, and ultimately finite-time blow-up behaviour, for localized spike solutions that occur as a parameter β is slowly ramped in time below various linear stability and existence thresholds associated with steady-state spike solutions. The transition or route to an ultimate finite-time blow-up can include spike nucleation, spike annihilation, or spike amplitude oscillation, depending on the specific parameter regime. Our detailed analysis of the existence and linear stability of multi-spike patterns, through the analysis of an explicitly solvable non-local eigenvalue problem, provides a theoretical guide for predicting which transition will be realized. Finally, we analyse the blow-up profile for a shadow limit of the reaction–diffusion system. For the resulting non-local scalar parabolic problem, we derive an explicit expression for the blow-up rate near the parameter range where blow-up is predicted. This blow-up rate is confirmed with full numerical simulations of the full PDE. Moreover, we analyse the linear stability of this solution that blows up in finite time.

This paper analyzes the existence of Hopf bifurcation and establishes the conditions under which the equilibrium path converges toward periodic solutions in a one-sector optimal growth model with delay. We establish the limits and the possibilities of nonlinear dynamics (i.e., cycles) vis-à-vis delays. In particular, we formulate a new method to further comprehend the root distribution of the characteristic equation of a standard optimal growth model with delayed investment structure. We show that nonmonotonic dynamics (limit cycles, persistent oscillations) occurs when the delayed investment causes permanent adjustment failures among the economic variables in the economy.

In this paper, a nonlinear stage-structured model for Lyme disease is considered. The model is a system of diòerential equations with two time delays. We derive the basic reproductive rate ${{R}_{0}}\left( {{\tau }_{1}},\,{{\tau }_{2}} \right)$. If ${{R}_{0}}\left( {{\tau }_{1}},\,{{\tau }_{2}} \right)\,<\,1$, then the boundary equilibrium is globally asymptotically stable. If ${{R}_{0}}\left( {{\tau }_{1}},\,{{\tau }_{2}} \right)\,>\,1$, then there exists a unique positive equilibrium whose local asymptotic stability and the existence of Hopf bifurcations are established by analyzing the distribution of the characteristic values. An explicit algorithm for determining the direction of Hopf bifurcations and the stability of the bifurcating periodic solutions is derived by using the normal form and the center manifold theory. Some numerical simulations are performed to confirm the correctness of theoretical analysis. Finally, some conclusions are given.

We consider a coupled, logistic predator–prey system with delay. Mainly, by choosing the delay time ${\it\tau}$ as a bifurcation parameter, we show that Hopf bifurcation can occur as the delay time ${\it\tau}$ passes some critical values. Based on the normal-form theory and the centre manifold theorem, we also derive formulae to obtain the direction, stability and the period of the bifurcating periodic solution at critical values of ${\it\tau}$. Finally, numerical simulations are investigated to support our theoretical results.

We examine the dynamics of fermentation process in a yeast cell. Our investigation focuses on the main branch pathway: pyruvate and acetaldehyde branch points. We formulate the kinetics for all enzymatic reactions as Michaelis–Menten models. Since the activity of an enzyme mainly depends on the conformational changes of the enzyme structure, the enzyme requires a certain period of time to reset its structure, until it is ready to bind substrates again. For this situation, a rate-limiting step exists, for which the catalytic process suffers a delay. Since all conversion processes are catalysed by enzymes, each reaction can experience a delay at a different time. To investigate how the delay affects the reaction processes, especially at the branch points, we propose that the rate-limiting step takes place at the first reaction. For this reason, a discrete time delay is introduced to the first kinetic model. We find a bifurcation diagram for the delay that depends on the rate of glucose supply and kinetic parameters of the first enzyme. By comparison, our analysis agrees with the numerical solution. Our numerical simulations also show that there is a certain glucose supply that may optimize ethanol production.

The aim of this paper is to give a detailed analysis of Hopf bifurcation of a ratio-dependent predator–prey system involving two discrete delays. A delay parameter is chosen as the bifurcation parameter for the analysis. Stability of the bifurcating periodic solutions is determined by using the centre manifold theorem and the normal form theory introduced by Hassard et al. Some of the bifurcation properties including the direction, stability and period are given. Finally, our theoretical results are supported by some numerical simulations.

A mathematical model, coupling the dynamics of short-term stem-like cells and matureleukocytes in leukemia with that of the immune system, is investigated. The model isdescribed by a system of seven delay differential equations with seven delays. Threeequilibrium points E0, E1,E2 are highlighted. The stability and the existence of theHopf bifurcation for the equilibrium points are investigated. In the analysis of themodel, the rate of asymmetric division and the rate of symmetric division are veryimportant.

A two dimensional two-delays differential system modeling the dynamics of stem-like cellsand white-blood cells in Chronic Myelogenous Leukemia is considered. All three types ofstem cell division (asymmetric division, symmetric renewal and symmetric differentiation)are present in the model. Stability of equilibria is investigated and emergence ofperiodic solutions of limit cycle type, as a result of a Hopf bifurcation, is eventuallyshown. The stability of these limit cycles is studied using the first Lyapunovcoefficient.