For the adult-survival annual pulse-breeding predator-prey model, a positive annual fixed point exists when and :
At this point the Jacobian matrix is
If , the two eigenvalues are a complex conjugate pair of moduli one. The annual samples show linearly neutral seasonal oscillations. If , they form a reciprocal negative real pair and the fixed point is a saddle fixed point of a map. Equality at four is a degenerate boundary needing nonlinear analysis. For replacement rather than surviving adults the same calculation gives trace and determinant one.
The absence of an attracting annual cycle is stronger than a linear observation. Put , . The update becomes
It is the composition of two nonlinear shear maps, each of determinant one: an area-preserving seasonal predator-prey map. No isolated positive periodic orbit has an open attracting basin. Neutral linear oscillations alone also do not prove nonlinear stability at every resonance. In particular, this density-independent model does not select an asymptotically stable oscillation amplitude.
The Nicholson-Bailey model instead counts parasitoid offspring from attacked hosts:
Its positive equilibrium has and . Its Jacobian matrix has
Since the determinant exceeds one while the positive trace is below two, the eigenvalues are complex with moduli greater than one: instability of the Nicholson-Bailey equilibrium produces outward oscillations near it. The two maps share an exponential prey-survival factor but have different reproduction laws and stability properties.
Figure 1.
Neutral seasonal predator-prey oscillations compared with the unstable Nicholson-Bailey equilibrium
. Details.
At a fixed point , the quadratic equation is . Consequently
The Jacobian matrix is , and its Floquet multipliers solve . The Jury stability criterion requires and , . With , this becomes , . The lower branch never satisfies it; the upper branch is attracting precisely when
At the lower boundary the two fixed points merge with Floquet multipliers , a saddle-node bifurcation for . At the upper boundary the upper point has Floquet multipliers , giving a generic period-doubling bifurcation for .
Its direction can be established explicitly. A distinct two-cycle alternates values satisfying , . Subtracting gives , and substitution gives
It exists beyond the fixed-point period-doubling bifurcation. The second-iterate Jacobian matrix has determinant and trace . Applying the same stability inequalities shows attraction for and
Thus the first period-doubling bifurcation is supercritical; the next upper boundary has a Floquet multiplier for the two-cycle. These are the fixed-point and two-cycle thresholds of the Hénon map.
At the map preserves phase-space area, so attraction is impossible. The fixed-point fold is at , with double Floquet multiplier , and is a conservative saddle-centre bifurcation. The upper branch is linearly elliptic for ; the lower branch is a saddle fixed point of a map. At , the upper branch has double Floquet multiplier and becomes a saddle fixed point of a map for larger . The two-cycle is elliptic for , then reaches double Floquet multiplier for the second iterate at . The partial sketch distinguishes elliptic, neutrally stable Floquet multipliers from attracting states. Resonances may require additional nonlinear analysis; unit Floquet multiplier moduli alone do not prove nonlinear Lyapunov stability.
The segment is a conservative Hénon stability boundary, not an ordinary nondegenerate dissipative Neimark–Sacker bifurcation. Indeed for constant phase-space area scaling is . A simple invariant closed curve would enclose a bounded invariant set of positive phase-space area, forcing . Thus no such closed invariant curve can bifurcate into within this constant-determinant family.
Generic folds and flips away from these degeneracies persist under sufficiently small smooth perturbations: their simple unit Floquet multiplier, transverse parameter crossing and nonzero leading nonlinear coefficient persist, while the other Floquet multiplier stays away from the unit circle. The implicit function theorem then continues the bifurcation curves and their criticality. The unqualified perturbation claim needs this restriction. Arbitrary perturbations need not preserve the conservative behavior. For example, in local complex coordinates at a nonresonant elliptic fixed point, the area-preserving map property gives a cubic normal form with zero cubic radial damping. Composing with the arbitrarily small radial perturbation introduces radial coefficient . Varying the modulus through one can now give a genuine Neimark–Sacker bifurcation. Conservative degeneracies and resonances are therefore not covered by a blanket structural-stability assertion.
Figure 1.
Hénon fixed points and two-cycles at unit area determinant, and local branches of the triangular cubic map
. Details.