Abstract In this paper, we study a predator–prey mite model of Leslie type with generalized Holling IV functional response. The model is shown to have very rich bifurcation dynamics, including subcritical and supercritical Hopf bifurcations, degenerate Hopf bifurcation, focus‐type and cusp‐type degenerate Bogdanov–Takens bifurcations of codimension 3, originating from a nilpotent focus or cusp of codimension 3 that acts as the organizing center for the bifurcation set. Coexistence of multiple steady states, multiple limit cycles, and homoclinic cycles is also found. Interestingly, the coexistence of two limit cycles is guaranteed by investigating generalized Hopf bifurcation and degenerate homoclinic bifurcation, and we also find that two generalized Hopf bifurcation points are connected by a saddle‐node bifurcation curve of limit cycles, which indicates the existence of global regime for two limit cycles. Our work extends some results in the literature.
more »
« less
A degenerate bifurcation from simple eigenvalue theorem
A new bifurcation from simple eigenvalue theorem is proved for general nonlinear functional equations. It is shown that in this bifurcation scenario, the bifurcating solutions are on a curve which is tangent to the line of trivial solutions, while in typical bifurcations the curve of bifurcating solutions is transversal to the line of trivial ones. The stability of bifurcating solutions can be determined, and examples from partial differential equations are shown to demonstrate such bifurcations.
more »
« less
- Award ID(s):
- 1853598
- PAR ID:
- 10333069
- Date Published:
- Journal Name:
- Electronic Research Archive
- Volume:
- 30
- Issue:
- 1
- ISSN:
- 2688-1594
- Page Range / eLocation ID:
- 116 to 125
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
In this work we present a systematic review of novel and interesting behaviour we have observed in a simplified model of a MEMS oscillator. The model is third order and nonlinear, and we expressit as a single ODE for a displacement variable. We find that a single oscillator exhibits limitcycles whose amplitude is well approximated by perturbation methods. Two coupled identicaloscillators have in-phase and out-of-phase modes as well as more complicated motions.Bothof the simple modes are stable in some regions of the parameter space while the bifurcationstructure is quite complex in other regions. This structure is symmetric; the symmetry is brokenby the introduction of detuning between the two oscillators. Numerical integration of the fullsystem is used to check all bifurcation computations. Each individual oscillator is based on a MEMS structure which moves within a laser-driven interference pattern. As the structure vibrates, it changes the interference gap, causing the quantity of absorbed light to change, producing a feedback loop between the motion and the absorbed light and resulting in a limit cycle oscillation. A simplified model of this MEMS oscillator, omitting parametric feedback and structural damping, is investigated using Lindstedt's perturbation method. Conditions are derived on the parameters of the model for a limit cycle to exist. The original model of the MEMS oscillator consists of two equations: a second order ODE which describes the physical motion of a microbeam, and a first order ODE which describes the heat conduction due to the laser. Starting with these equations, we derive a single governing ODE which is of third order and which leads to the definition of a linear operator called the MEMS operator. The addition of nonlinear terms in the model is shown to produce limit cycle behavior. The differential equations of motion of the system of two coupled oscillators are numerically integrated for varying values of the coupling parameter. It is shown that the in-phase mode loses stability as the coupling parameter is reduced below a certain value, and is replaced by two new periodic motions which are born in a pitchfork bifurcation. Then as this parameter is further reduced, the form of the bifurcating periodic motions grows more complex, with yet additional bifurcations occurring. This sequence of bifurcations leads to a situation in which the only periodic motion is a stable out-of-phase mode. The complexity of the resulting sequence of bifurcations is illustrated through a series of diagrams based on numerical integration.more » « less
-
We study equilibrium configurations of double-stranded DNA in a cylindrical viral capsid. The state of the encapsidated DNA consists of a disordered inner core enclosed by an ordered outer region, next to the capsid wall. The DNA configuration is described by a unit helical vector field, tangent to an associated centre curve, passing through properly selected locations. We postulate an expression for the energy of the encapsulated DNA based on that of columnar chromonic liquid crystals. A thorough analysis of the Euler–Lagrange equations yields multiple solutions. We demonstrate that there is a trivial, non-helical solution, together with two solutions with non-zero helicity of opposite sign. Using bifurcation analysis, we derive the conditions for local stability and determine when the preferred coiling state is helical. The bifurcation parameters are the ratio of the twist versus the bend moduli of DNA and the ratio between the sizes of the ordered and the disordered regions.more » « less
An official website of the United States government

