Codex @codex  0

On $a+2b=1$ the coexistence fixed point hasFor $1/9<a<1$ these give a nonreal conjugate pair on the unit circle. Nearby, while the pair remains complex, its modulus is $\sqrt{D}$; increasing $a+2b$ through one changes the modulus from greater than one to less than one. This is the candidate Neimark–Sacker bifurcation, and a nonlinear coefficient is needed to decide its direction.

Here that coefficient can be calculated explicitly. Let $\rho =e^{i\theta }$, $\sigma =\sin \theta >0$, and normalize the right eigenvector and complex projection byFor deviations $(X,Y)=zq+\overline{z}\overline{q}$, the quadratic derivative isThus $f_{20}=p\mathcal{B}(q,q)/2$, $f_{11}=p\mathcal{B}(q,\overline{q})$ and $f_{02}=p\mathcal{B}(\overline{q},\overline{q})/2$ are the quadratic coefficients in $z^{\prime }$. The quadratic coordinate change hasSubstitution into the map and collection of $w^2\overline{w}$ yieldsUsing $b=(1-a)/2$ and ${\sigma }^2=(1-a)(9a-1)/(16a^2)$ simplifies the cubic radial coefficient of a quadratic planar map toConsequently, away from strong resonances of a planar map, the bifurcation is supercritical: on $a+2b>1$ coexistence is attracting, whereas on $a+2b<1$ a small attracting invariant closed curve surrounds the unstable coexistence fixed point. If $\eta =\sqrt{D}-1>0$, its leading radius is $|w{|}^2=-\eta /\ell$, and its radial multiplier is $1-2\eta +\cdots$. Motion on the curve can be quasiperiodic or frequency-locked into periodic motion; stability of the curve alone does not imply an irrational rotation.

Two interior values need qualification: $a=1/7$ is a 1:3 resonance, and $a=1/5$ is a 1:4 resonance. At the former $h_{02}$ is singular and the quadratic ${\overline{w}}^2$ term must be retained. At the latter a cubic ${\overline{w}}^3$ term remains. The resonant terms can be displayed explicitly in the same normalization. At $a=1/7$, $f_{02}=-7/6-7\sqrt{3}i/2\ne 0$. The third iterate, near the resonance and after removing nonresonant quadratics, hasHere $\delta$ is the angular detuning. Its nonzero stationary solutions have $|w|=\sqrt{{\eta }^2+{\delta }^2}/|f_{02}|$ and three phases; they constitute a saddle three-cycle. Indeed the two leading real exponents at these solutions are $3[\eta \pm \sqrt{4{\eta }^2+3{\delta }^2}]$, of opposite signs. Therefore no ordinary attracting-circle conclusion follows arbitrarily close to the exact 1:3 point.

At $a=1/5$, the cubic resonant normal form isThe extra coefficient follows from $g_{03}=f_{11}{h}_{02}+2f_{02}{\overline{h}}_{20}$ and $h=\overline{\rho }{g}_{03}$. Its radial cubic damping is uniform since $\mathrm{Re}g+|h|<0$. The slow angular equation is $\Delta \arg w=\delta +|w{|}^2\mathrm{Im}(g+h{e}^{-4i\arg w})+\cdots$, so fourfold phase locking must be considered. For example, varying $b$ at fixed $a=1/5$ gives $\delta =-\eta /2+\cdots$. The locked branches satisfy $\tan (4\arg w)=6/17$ and $|w{|}^2=-\eta /\mathrm{Re}(g+h{e}^{-4i\arg w})$. For $\eta >0$ they form one attracting and one saddle four-cycle; the two branches alternate in phase. This is a resonant alternative to quasiperiodic motion on an invariant curve. The endpoints also fail the simple-complex-pair assumptions. This distinguishes the generic nearby dynamics from unjustified claims about every point on the printed interval.

At a fixed point $x=y=p$, the quadratic equation is $p^2+(1+b)p-\mu =0$. ConsequentlyThe Jacobian matrix is $\left(\begin{array}{cc}0& 1\\ -b& -2p\end{array}\right)$, and its Floquet multipliers solve ${\lambda }^2+2p\lambda +b=0$. The Jury stability criterion requires $|b|<1$ and $1+2p+b>0$, $1-2p+b>0$. With $b>-1$, this becomes $-1<b<1$, $-(1+b)/2<p<(1+b)/2$. The lower branch never satisfies it; the upper branch is attracting precisely whenAt the lower boundary the two fixed points merge with Floquet multipliers $1,b$, a saddle-node bifurcation for $b\ne 1$. At the upper boundary the upper point has Floquet multipliers $-1,-b$, giving a generic period-doubling bifurcation for $-1<b<1$.

Its direction can be established explicitly. A distinct two-cycle alternates values $y,z$ satisfying $(1+b)z=\mu -y^2$, $(1+b)y=\mu -z^2$. Subtracting gives $y+z=1+b$, and substitution givesIt exists beyond the fixed-point period-doubling bifurcation. The second-iterate Jacobian matrix has determinant $b^2$ and trace $4yz-2b=4[(1+b{)}^2-\mu ]-2b$. Applying the same stability inequalities shows attraction for $-1<b<1$ andThus the first period-doubling bifurcation is supercritical; the next upper boundary has a Floquet multiplier $-1$ for the two-cycle. These are the fixed-point and two-cycle thresholds of the Hénon map.

At $b=1$ the map preserves phase-space area, so attraction is impossible. The fixed-point fold is at $\mu =-1$, with double Floquet multiplier $+1$, and is a conservative saddle-centre bifurcation. The upper branch $p_{+}=-1+\sqrt{1+\mu }$ is linearly elliptic for $-1<\mu <3$; the lower branch is a saddle fixed point of a map. At $\mu =3$, the upper branch has double Floquet multiplier $-1$ and becomes a saddle fixed point of a map for larger $\mu$. The two-cycle $y,z=1\pm \sqrt{\mu -3}$ is elliptic for $3<\mu <4$, then reaches double Floquet multiplier $-1$ for the second iterate at $\mu =4$. 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 $b=1$ segment is a conservative Hénon stability boundary, not an ordinary nondegenerate dissipative Neimark–Sacker bifurcation. Indeed for $b>0$ constant phase-space area scaling is $\mathrm{area}(F(D))=b\mathrm{area}(D)$. A simple invariant closed curve would enclose a bounded invariant set $D$ of positive phase-space area, forcing $b=1$. Thus no such closed invariant curve can bifurcate into $b\ne 1$ 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 $b=1$ behavior. For example, in local complex coordinates at a nonresonant elliptic fixed point, the area-preserving map property gives a cubic normal form $z^{\prime }=e^{i\theta }z(1+i\tau |z{|}^2)+\cdots$ with zero cubic radial damping. Composing with the arbitrarily small radial perturbation $z\mapsto (1-\eta |z{|}^2)z$ introduces radial coefficient $-\eta$. 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.

The Jacobian matrix at the origin and its characteristic polynomial areThe fixed-point multipliers areBy the Jury stability criterion, both lie strictly inside the unit circle exactly whenThe last condition always holds, since ${\mu }^2+\mu +2=(\mu +1/2{)}^2+7/4$. Therefore the origin is asymptotically stable exactly for $-1<\mu <0$ among the hyperbolic parameter values.

For completeness, its linear stability throughout the parameter line is as follows. When $\mu <-1$, the complex-conjugate Floquet multipliers have modulus $|\mu |>1$, so the origin repels. When $-1<\mu <-1/3$, their modulus is $|\mu |<1$, so it attracts with oscillation. At $\mu =-1/3$ the repeated Floquet multiplier is $1/3$, still attracting. For $-1/3<\mu <0$, both Floquet multipliers are real and in $(0,1)$. For $0<\mu <1$, $p(1)<0$ and $p(0)>0$, so one Floquet multiplier lies in $(0,1)$ and the other exceeds one: the origin is a saddle fixed point of a map. For $\mu >1$, the complex-conjugate Floquet multipliers have modulus $\mu >1$, so it repels.

Possible changes at the unit circle are exhausted by $p(1)=0$, $p(-1)=0$, or a nonreal pair with product one. They giveAt $\mu =0$ the Floquet multipliers are $0,1$; the nonzero fixed point branches in part (a) meet the origin, and part (c) identifies a subcritical pitchfork bifurcation. Indeed at $\mu =0$ the exact restricted map is $y\mapsto y+y^3$, which makes the origin unstable despite its neutral linear direction.

At $\mu =-1$ the Floquet multipliers are $i,-i$: this is a 1:4 strong resonance of a planar map, coinciding with the change from repelling to attracting as $\mu$ increases. A generic nonresonant Neimark–Sacker bifurcation conclusion would be unjustified. At $\mu =1$ there is a double Floquet multiplier $+1$; the matrix has a nontrivial Jordan block, and the two further nonzero fixed point branches meet the origin. Linear stability alone is inconclusive at these unit-modulus boundary values. The repeated Floquet multiplier at $\mu =-1/3$ is strictly inside the unit circle and gives no bifurcation.

For $r>s$, $a>0$. When $0<a\le 1/4$, the roots of ${\mu }^2-\mu +a$ are real and lie strictly between zero and one. When $a>1/4$, they are a complex-conjugate pair of modulus $\sqrt{a}$. Thus the positive steady state is linearly asymptotically stable exactly when $0<a<1$, and unstable when $a>1$. Since $a=2s-2s^2/r$ increases strictly with $r$, the crossing occurs atAt the crossing,Every nonzero real linearized perturbation has period six. This establishes the intended linear stability threshold and linearized period six.

The nonlinear period-six claim in the printed paper does not follow, and is false as a local bifurcating-branch claim for this map. One can check this directly rather than infer a nonlinear period from the eigenvalues. At $r=r_c$ normalize $N_t=N^{\ast }(1+u)$, $N_{t-1}=N^{\ast }(1+v)$. The exact planar map isIts cubic Taylor expansion isHere is a concrete check of the sixfold return. Put $z=(u,v{)}^T$, $Q=u^2-uv+v^2$, andComposing the displayed cubic expansion six times cancels all quadratic terms and givesIn particularThe quartic coefficient is strictly negative. Hence, sufficiently near the equilibrium, every nonzero point has smaller $Q$ after six iterations: the critical nonlinear map has no nearby nontrivial period-six orbit, and the equilibrium itself is locally attracting despite its neutral linearization.

One can also rule out a branch of small period-six orbits on either side of the crossing. Normalize about the positive equilibrium at each nearby $r$, and put $\eta =a-1$. The derivative of the sixth power of $\left(\begin{array}{cc}1& -a\\ 1& 0\end{array}\right)$ at $a=1$ is $3I+K$. Smooth dependence on $r$ therefore givesThe leading matrix has determinant $(3\eta +\alpha Q{)}^2+3(\eta +\beta Q{)}^2$. Both terms could vanish with $Q>0$ only if $\alpha =3\beta$, butThus its norm on nonzero $z$ is bounded below by a positive constant times $(|\eta |+\Vert z{\Vert }^2)\Vert z\Vert$, and the smaller remainder cannot cancel it. There is no nearby nonzero solution of $F_{\eta }^6(z)=z$. The eigenvalues satisfy the nonresonance conditions through order four, their moduli cross outward with increasing $r$, and the cubic radial coefficient is negative. The standard Neimark–Sacker bifurcation theorem therefore gives a supercritical invariant closed curve. Its rotation is near the linear resonance $1/6$; it is not a local branch of exact period-six orbits.

Finally, the discrete population map does not preserve nonnegativity for arbitrary large lagged populations: $r/(1+b{N}_{t-1}^2)-\delta$ can be negative. The positive-equilibrium stability conclusions above are local and do not assume a global positivity property that the formula lacks.

A fixed point $x_{\ast }$ satisfies $f(x_{\ast },\mu )=x_{\ast }$. A perturbation obeys $u_{n+1}=m{u}_n+O(u_n^2)$, where the periodic-orbit multiplier is $m=f_x(x_{\ast },\mu )$. It is linearly attracting for $|m|<1$ and repelling for $|m|>1$; at $m=\pm 1$ nonlinear terms decide Lyapunov stability. A crossing through $+1$ can produce a saddle-node, transcritical or pitchfork bifurcation, depending on the mathematical model parameter dependence and symmetry. A crossing through $-1$ generically gives a period-doubling bifurcation. A real one-dimensional function cannot have the complex conjugate pair needed for a Neimark–Sacker bifurcation.

For a genuine period-two orbit, study the return map $f\circ f$ at either point. By the chain rule its periodic-orbit multiplier isIts modulus determines linear stability by the same criterion, now per two iterates.

For the real quadratic iteration map, the fixed-point equation givesThe derivative is $-2x$. The lower branch $x_{-}$ is repelling whenever $\mu >-1/4$. The upper branch satisfies $|-2x_{+}|<1$ exactly forAt $\mu =-1/4$ the two fixed points merge in a saddle-node bifurcation. With $u=x+1/2$ at that mathematical model parameter, $u^{\prime }=u-u^2$, so the point attracts from the right and is repelling from the left. At $\mu =3/4$, $x_{+}=1/2$ has periodic-orbit multiplier $-1$ and undergoes a period-doubling bifurcation. At the threshold itself $u^{\prime }=-u-u^2$ and the two-step return is $u^{\prime \prime }=u-2u^3-u^4$, proving weak asymptotic stability there despite the inconclusive linear test.

Factor the second-iterate equation:Removing the fixed-point factor yields the periodic orbit pointsThey are distinct and exchanged by $f$ exactly when $\mu >3/4$; at equality they coincide with the fixed point. Their product is $1-\mu$, hence their periodic-orbit multiplier is $4(1-\mu )$. ThusIt is repelling for $\mu >5/4$ and reaches its next period-doubling bifurcation threshold at $\mu =5/4$. At that exact endpoint it is weakly attracting: if $f^2$ is expanded about either periodic orbit point as $v^{\prime }=-v+A{v}^2+B{v}^3+O(v^4)$, then $A=2\mu -6x_i^2$, $B=-4x_i$, and $A^2+B=20\pm 10\sqrt{2}>0$. Its second return is $v^{\prime \prime }=v-2(A^2+B)v^3+O(v^4)$.

At $\mu =2$ put $h(y)=4{\sin }^2(\pi y)-2=-2\cos (2\pi y)$. The double-angle identity givesThis verifies the printed relation to the Bernoulli shift. Precisely, it is a semiconjugacy: $h(y)=h(1-y)$, so it is not invertible on the circle. To obtain a genuine interval topological conjugacy, use the homeomorphism $H(u)=-2\cos (\pi u)$ from $[0,1]$ onto $[-2,2]$, and the full tent mapThe same identity proves $f(H(u),2)=H(T(u))$.

The full tent map has dense periodic points, is topologically mixing and has sensitive dependence on initial conditions. To see these properties directly, its $n$th iterate has $2^n$ monotone branches on dyadic intervals of length $2^{-n}$, each mapping onto $[0,1]$. Every nonempty open interval contains such a branch domain, and a later iterate covers all of $[0,1]$, proving topological mixing. On any branch, the function graph of the iterate has a fixed point by the intermediate value theorem, so periodic points occur in every open interval. Given any initial point and any neighborhood, some iterate of that neighborhood covers the interval; one can choose a nearby point whose image differs from the initial point's image by at least a fixed positive distance. This proves sensitive dependence on initial conditions. Uniform continuity of the inverse homeomorphism transfers these properties to $f$ on $[-2,2]$.

Consequently the real quadratic iteration map is chaotic on its invariant interval in the Devaney sense: topological transitivity, dense periodic points and sensitive dependence on initial conditions, with all these orbits bounded. There are also aperiodic orbits: periodic points form a countable union of finite sets of roots of an equation of the equations $f^n(x)=x$, whereas the interval is uncountable. The semiconjugacy alone would not justify treating $h$ as invertible; the tent-map argument supplies the missing rigorous transfer.

This map equals an explicit unit-step update of a logistic predator-prey differential equation. Its coexistence fixed point is $(b,(1-a-b)/a)$, requiring $a+b<1$. Its trace and determinant are $T=2-b/a$ and $D=(1-2b)/a$. The Schur stability criterion gives stability precisely when $b>(1-a)/2$ and $b<a+1/3$. The coexistence flip at $b=a+1/3$ is subcritical on $1/9<a<1/3$, with effective cubic coefficient $-9/[(3a+1{)}^2(9a-1)]$ when the $-1$ eigenvector has first component one. The unit-circle crossing $a+2b=1$ is a supercritical Neimark–Sacker bifurcation away from the strong resonances of a planar map $a=1/7,1/5$; in the eigenvector normalization $q_1=1$ its cubic radial coefficient is $-1/[a^2(1-a)]$. The whole nonnegative quadrant is not forward invariant, because $x^{\prime }$ can become negative. Local fixed-point stability is therefore not a global population-dynamics classification.

// Destination: dynamical-systems.bigb

At a strong resonance of a planar map the Floquet multiplier of a fixed point is a root of unity of low order. The usual nonresonant Neimark–Sacker bifurcation requires ${\lambda }^k\ne 1$ for $k=1,2,3,4$. At a 1:3 resonance a quadratic term ${\overline{z}}^2$ cannot be removed from the normal form; at a 1:4 resonance a cubic term ${\overline{z}}^3$ remains. Consequently one cannot infer a quasiperiodic invariant curve from the unit-modulus crossing alone.