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.

We study the Cauchy problem on the real line for the nonlocal Fisher-KPP equation in one spatial dimension,

\begin{equation*} u_t = D u_{xx} + u(1-\phi *u), \end{equation*}

$\phi *u$

$\phi (y) \equiv H\left (\frac{1}{4}-y^2\right )$

$u=1$

$D$

$\Delta _1 \approx 0.00297$

$D \ll 1$

$O(1)$

$u=O(1)$

$u$

$D \to 0^+$

$D \geq \Delta _1$

$2 \sqrt{D}$

$0 \lt D \lt \Delta _1$

$u=0$

$D$

Motivated by the impact of worsening climate conditions on vegetation patches, we study dynamic instabilities in an idealised Ginzburg–Landau model. Our main results predict time instances of sudden drops in wavenumber and the resulting target states. The changes in wavenumber correspond to the annihilation of individual vegetation patches when resources are scarce and cannot support the original number of patches. Drops happen well after the primary pattern has destabilised at the Eckhaus boundary and key to distinguishing between the disappearance of 1,2 or more patches during the drop are complex spatio-temporal resonances in the linearisation at the unstable pattern. We support our results with numerical simulations and expect our results to be conceptually applicable universally near the Eckhaus boundary, in particular in more realistic models.

We investigate the equilibrium configurations of closed planar elastic curves of fixed length, whose stiffness, also known as the bending rigidity, depends on an additional density variable. The underlying variational model relies on the minimisation of a bending energy with respect to shape and density and can be considered as a one-dimensional analogue of the Canham–Helfrich model for heterogeneous biological membranes. We present a generalised Euler–Bernoulli elastica functional featuring a density-dependent stiffness coefficient. In order to treat the inherent nonconvexity of the problem, we introduce an additional length scale in the model by means of a density gradient term. We derive the system of Euler–Lagrange equations and study the bifurcation structure of solutions with respect to the model parameters. Both analytical and numerical results are presented.

In this paper, we prove existence results of a one-dimensional periodic solution to equations with the fractional Laplacian of order $s\in (1/2,1)$, singular nonlinearity and gradient term under various situations, including nonlocal contra-part of classical Lienard vector equations, as well other nonlocal versions of classical results know only in the context of second-order ODE. Our proofs are based on degree theory and Perron's method, so before that we need to establish a variety of priori estimates under different assumptions on the nonlinearities appearing in the equations. Besides, we obtain also multiplicity results in a regime where a priori bounds are lost and bifurcation from infinity occurs.

We consider patterns formed in a two-dimensional thin film on a planar substrate with a Derjaguin disjoining pressure and periodic wettability stripes. We rigorously clarify some of the results obtained numerically by Honisch et al. [Langmuir 31: 10618–10631, 2015] and embed them in the general theory of thin-film equations. For the case of constant wettability, we elucidate the change in the global structure of branches of steady-state solutions as the average film thickness and the surface tension are varied. Specifically we find, by using methods of local bifurcation theory and the continuation software package AUTO, both nucleation and metastable regimes. We discuss admissible forms of spatially non-homogeneous disjoining pressure, arguing for a form that differs from the one used by Honisch et al., and study the dependence of the steady-state solutions on the wettability contrast in that case.

This paper is concerned with a prey–predator model with population flux by attractive transition. Our previous paper (Oeda and Kuto, 2018, Nonlinear Anal. RWA, 44, 589–615) obtained a bifurcation branch (connected set) of coexistence steady states which connects two semitrivial solutions. In Oeda and Kuto (2018, Nonlinear Anal. RWA, 44, 589–615), we also showed that any positive steady-state approaches a positive solution of either of two limiting systems, and moreover, one of the limiting systems is an equal diffusive competition model. This paper obtains the bifurcation structure of positive solutions to the other limiting system. Moreover, this paper implies that the global bifurcation branch of coexistence states consists of two parts, one of which is a simple curve running in a tubular domain near the set of positive solutions to the equal diffusive competition model, the other of which is a connected set characterized by positive solutions to the other limiting system.

An analysis is undertaken of the formation and stability of localised patterns in a 1D Schanckenberg model, with source terms in both the activator and inhibitor fields. The aim is to illustrate the connection between semi-strong asymptotic analysis and the theory of localised pattern formation within a pinning region created by a subcritical Turing bifurcation. A two-parameter bifurcation diagram of homogeneous, periodic and localised patterns is obtained numerically. A natural asymptotic scaling for semi-strong interaction theory is found where an activator source term \[a = O(\varepsilon )\] and the inhibitor source \[b = O({\varepsilon ^2})\], with ε2 being the diffusion ratio. The theory predicts a fold of spike solutions leading to onset of localised patterns upon increase of b from zero. Non-local eigenvalue arguments show that both branches emanating from the fold are unstable, with the higher intensity branch becoming stable through a Hopf bifurcation as b increases beyond the \[O(\varepsilon )\] regime. All analytical results are found to agree with numerics. In particular, the asymptotic expression for the fold is found to be accurate beyond its region of validity, and its extension into the pinning region is found to form the low b boundary of the so-called homoclinic snaking region. Further numerical results point to both sub and supercritical Hopf bifurcation and novel spikeinsertion dynamics.

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.

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.

This is the second part of our study on the spatially heterogeneous predator–prey model where the interaction is governed by a Crowley–Martin type functional response. In part I, we have proved that when the predator competition is strong (i.e. k is large), the model has at most one positive steady-state solution for any $\mu \in \mathbb {R}$, moreover it is globally asymptotically stable for any $\mu >0$. This part is denoted to study the effect of saturation. Our result shows that the large saturation coefficient (i.e. large m) can not only lead to the uniqueness of positive solutions, but also lead to the multiplicity of positive solutions, moreover the stability of the corresponding positive solutions is also completely obtained. This work can be regarded as a supplement of Ref. [10].

We consider the parabolic, initial-boundary value problem 1

$$\matrix{ {\displaystyle{{\partial v} \over {\partial t}} = \Delta _p(v) + f(x,v),} & {{\rm in}({\rm - 1},{\rm 1}) \times ({\rm 0},\infty ),} \cr {v( \pm 1,t) = 0,} \hfill \hfill \hfill & {{\rm t}\in [{\rm 0},\infty ),} \hfill \hfill \cr {v = v_0\in C_0^0 ([-1,1]),} & {{\rm in}[{\rm - 1},{\rm 1}] \times \{ {\rm 0}\} ,} \cr } $$

We also consider a bifurcation type problem similar to (1), having a line of trivial solutions. We characterize the stability or instability of the trivial solutions, and the bifurcating, non-trivial solutions, and show that there is an ‘exchange of stability’ at the bifurcation point, analogous to the well-known result when p = 2.

In the theory of spontaneous combustion, identifying the critical value of the Frank-Kamenetskii parameter corresponds to solving a bifurcation point problem. There are two different numerical methods used to solve this problem—the direct and indirect numerical methods. The latter finds the bifurcation point by solving a partial differential equation (PDE) problem. This is a better method to find the bifurcation point for complex geometries. This paper improves the indirect numerical method by combining the grid-domain extension method with the matrix equation computation method. We calculate the critical parameters of the Frank-Kamenetskii equation for some complex geometries using the indirect numerical method. Our results show that both the curve of the outer boundary and the height of the geometries have an effect on the values of the critical Frank-Kamenetskii parameter, however, they have little effect on the critical dimensionless temperature.

In this paper we study the problem

where Ω = ℝ N or Ω = B 1, N ⩾ 3, p > 1 and . Using a suitable map we transform problem (1) into another one without the singularity 1/|x|2. Then we obtain some bifurcation results from the radial solutions corresponding to some explicit values of λ.

Of concern is the stability and bifurcation analysis of a free boundary problem modelling the growth of multi-layer tumours. A remarkable feature of this problem lies in that the free boundary is imposed with nonlinear boundary conditions, where a Gibbs–Thomson relation is taken into account. By employing a functional approach, analytic semigroup theory and bifurcation theory, we prove that there exists a positive threshold value γ* of surface tension coefficient γ such that if γ > γ* then the unique flat stationary solution is asymptotically stable under non-flat perturbations, while for γ < γ* this unique flat stationary solution is unstable and there exists a series of non-flat stationary solutions bifurcating from it. The result indicates a significant phenomenon that a smaller value of surface tension coefficient γ may make tumours more aggressive.

Let Ω be a smooth bounded domain in RN, with N ≥ 5. We provide existence and bifurcation results for the elliptic fourth-order equation Δ2u − Δpu = f(λ, x, u) in Ω, under the Dirichlet boundary conditions u = 0 and ∇u = 0. Here λ is a positive real number, 1 < p ≤ 2# and f(.,., u) has a subcritical or a critical growth s, 1 < s ≤ 2*, where and . Our approach is variational, and it is based on the mountain-pass theorem, the Ekeland variational principle and the concentration-compactness principle.

We consider the stationary non-linear Schrödinger equation\begin{equation*}\Delta u + \{1 + \lambda g(x)\} u = f(u)\mbox{with}u \in H^{1} (\mathbb{R}^{N}), u \not\equiv 0,\end{equation*} where $\lambda >0$ and the functions $f$ and $g$ are such that\begin{equation*} \lim_{s \rightarrow 0}\frac{f(s)}{s} = 0 \mbox{and} 1 < \alpha + 1 = \lim _{|s| \rightarrow \infty}\frac{f(s)}{s} < \infty\end{equation*} and \begin{equation*} g(x)\equiv 0 \mbox{on} \bar{\Omega}, g(x)\in (0, 1] \mbox{on} {\mathbb{R}^{N}} \setminus {\overline{\Omega}} \mbox{and} \lim_{|x| \rightarrow + \infty} g(x) = 1 \end{equation*} for some bounded open set $\Omega \in \mathbb{R}^{N}$. We use topological methods to establish the existence of two connected sets $\mathcal{D}^{\pm}$ of positive/negative solutions in $\mathbb{R} \times W^{2, p} (\mathbb{R}^{N})$ where $p \in [2, \infty) \cap (\frac{N}{2},\infty)$ that cover the interval $(\alpha,\Lambda(\alpha))$ in the sense that \begin{align*} P \mathcal{D}^{\pm} & = (\alpha, \Lambda(\alpha)) \text{where}P(\lambda, u) = \lambda \text{and furthermore,} \\ \lim_{\lambda \rightarrow \Lambda(\alpha)-}\left\Vert u_{\lambda} \right\Vert _{L^{\infty} (\mathbb{R}^{N})} & = \lim_{\lambda \rightarrow \Lambda (\alpha )-} \left\Vert u_{\lambda} \right\Vert _{W^{2, p}(\mathbb{R}^{N})} = \infty \text{ for }(\lambda, u_{\lambda}) \in \mathcal{D}^{\pm}. \end{align*} The number $\Lambda(\alpha)$ is characterized as the unique value of $\lambda$ in the interval $(\alpha, \infty)$ for which the asymptotic linearization has a positive eigenfunction. Our work uses a degree for Fredholm maps of index zero.

An implicitly defined second order ODE is said to be singular if the second derivative cannot be smoothly written in terms of lower order variables. The standard existence and uniqueness theory cannot be applied to such ODE and the graphs of solutions may fail to be regular curves (i.e., the solutions may have isolated C0-points or may fail to be single valued). In this paper we describe a local analysis for a large class of implicit second order ODE whose singular points satisfy a regularity condition. Within this class of ODE there is a secondary notion of (contact) singularity which is analogous to rest points for regular ODE. Theorems 5, 6, 7 and 8 produce invariants for these singularities which control the existence, uniqueness and the level of regularity in solutions.