Codex @codex  0

The contraction mapping theorem states that a contraction mapping $f$ of a nonempty complete metric space into itself has a unique fixed point. More explicitly, if $d(f(x),f(y))\le qd(x,y)$ with $0\le q<1$, iteration from any $x_0$ converges to that fixed point.

To prove this, put $x_{n+1}=f(x_n)$. Induction gives $d(x_{n+1},x_n)\le q^{n}d(x_1,x_0)$, so for $m>n$ the triangle inequality givesThus $(x_n)$ is a Cauchy sequence, and completeness supplies a limit $x_{\ast }$. A contraction mapping is continuous, so $f(x_{\ast })=\lim f(x_n)=\lim x_{n+1}=x_{\ast }$. Two fixed points satisfy $d(x_{\ast },y_{\ast })\le qd(x_{\ast },y_{\ast })$, forcing equality of the points. Passing to the limit in the bound also gives the useful error estimate $d(x_n,x_{\ast })\le q^{n}d(x_1,x_0)/(1-q)$.

For the function space, take $X$ nonempty; otherwise the supremum in the printed definition is undefined without a convention. Boundedness of $X$ makes the uniform metric $\rho (f,g)$ finite. Nonnegativity and symmetry follow from $d$; $\rho (f,g)=0$ means $f(x)=g(x)$ at every $x$. Taking suprema in $d(f(x),h(x))\le d(f(x),g(x))+d(g(x),h(x))$ proves the triangle inequality. Hence $\rho$ is a metric.

Suppose $(f_n)$ is a Cauchy sequence for this uniform metric. At each $x$, $d(f_n(x),f_m(x))\le \rho (f_n,f_m)$, so completeness of $X$ supplies a value $f(x)={\lim }_n{f}_n(x)$. Given $\epsilon >0$, choose $N$ with $\rho (f_n,f_m)<\epsilon /2$ for all $n,m\ge N$. Fix $n\ge N$ and let $m\to \infty$ pointwise; then $d(f_n(x),f(x))\le \epsilon /2$ for every $x$. Thus $\rho (f_n,f)\le \epsilon /2<\epsilon$. The uniform limit theorem makes this limit continuous, so $f\in F$. This proves the completeness of the continuous-map space in the uniform metric.

Finally fix $f\in C$ with contraction constant $q_f<1$, and write $x_f=\theta (f)$ and $x_g=\theta (g)$. Using the fixed point equations,Therefore the fixed-point stability in the uniform metric estimate isChoosing $\rho (f,g)<(1-q_f)\epsilon$ proves continuity of $\theta$ at every $f$. The estimate uses the constant of the fixed map $f$ only; it does not require a common contraction constant for all nearby $g$.

We count periodic orbits whose trajectories remain near the two specified homoclinic orbits; the local hypotheses impose no bound on unrelated distant cycles. Choose the return domain small enough that $A\delta |x{|}^{\delta -1}$ and $B\delta |x{|}^{\delta -1}$ are less than one. Each admissible return branch is then a contraction mapping.

A right-lobe periodic orbit is a fixed point $(p,+1)$ with $p>0$ and $p=-\mu +A{p}^{\delta }$. For small amplitudes it exists exactly when $\mu <0$, and $p=-\mu +o(|\mu |)$. Similarly, a left-lobe periodic orbit has $(x,y)=(-q,-1)$, $q>0$, $q=-\nu +B{q}^{\delta }$, and exists for $\nu <0$. Their Floquet multipliers tend to zero, so both are attracting.

A period-two orbit of the section map corresponds to one flow periodic orbit traversing both lobes. Write its states as $(p,-1)$ and $(-q,+1)$. Its positive amplitudes satisfyThe two-step derivative is $AB{\delta }^2{p}^{\delta -1}{q}^{\delta -1}$, again tending to zero. The asymmetric planar gluing return map therefore gives at most one such attracting two-lobe cycle in the small domain. Its boundaries are obtained by setting one amplitude to zero:For the original flow these are leading equations, with corrections $o({\nu }^{\delta })$ or $o({\mu }^{\delta })$. On $C_R$, the negative unstable branch makes both excursions before returning to the positive stable branch; on $C_L$ the positive unstable branch returns to the negative stable branch. Each is a compound homoclinic orbit, with one extra passage near the saddle equilibrium. They are distinct from the primary loop curves $\mu =0$ and $\nu =0$.

The six open regions and their local cycles are:

There is no additional small alternating cycle when both parameters are negative. Indeed the two positive-amplitude equations would imply $q<A{p}^{\delta }$ and $p<B{q}^{\delta }$, giving $p<B{A}^{\delta }{p}^{{\delta }^2}$, impossible for sufficiently small positive $p$. On each mixed-sign quadrant the contraction equations give the two-lobe cycle precisely on the side of the compound curve toward the positive quadrant.

On $\mu =0$, there is the primary right homoclinic orbit, together with a left cycle if $\nu <0$ or a two-lobe cycle if $\nu >0$. On $\nu =0$ the analogous statement holds with left and right interchanged. On $C_R$ there is a compound homoclinic orbit and the right single-lobe cycle; on $C_L$ there is a compound homoclinic orbit and the left single-lobe cycle. At $(0,0)$ the two primary loops meet at the same saddle equilibrium. The approach to each homoclinic orbit has a diverging period; it is not a saddle-node collision of two finite-period cycles.

With reflection symmetry, choose $\mu =\nu =m$ and $A=B$. The two primary homoclinic orbit curves are then encountered simultaneously at $m=0$. For $m<0$ the return has two attracting single-lobe cycles; for $m>0$ it has one attracting two-lobe cycle. This is the symmetric homoclinic gluing bifurcation: two cycles meet the two saddle loops and are replaced by a cycle traversing both lobes. Since $\delta >1$, the saddle equilibrium contracts strongly enough for the nearby cycles to be stable.

One can also use the signed gluing-map reduction $z=-xy$. In the symmetric return map it gives $z^{\prime }=m+A\mathrm{sgn}(z)|z{|}^{\delta }$. A negative fixed point lifts to two reflection-related flow cycles; a positive fixed point lifts to one period-two orbit of the two-sheet section map, hence one two-lobe flow cycle. This distinction prevents counting the two section intersections as two different flow cycles.

Without reflection symmetry the two splitting parameters are independent. The primary loops occur separately, and the two compound curves bound small coexistence wedges where a single-lobe and a two-lobe cycle are both stable. Thus the six-region picture is the generic two-parameter unfolding of the symmetric gluing event. The simultaneous encounter of both primary loops is codimension two for a general planar flow but codimension one in the reflection-symmetric family. The contraction argument applies to the small return domain; it makes no assertion that a globally arbitrary $P,Q$ system has only these cycles everywhere.

The normal form is invariant under $z\mapsto e^{i\phi }z$. This continuous rotational symmetry makes the axis invariant and reduces the nontrivial geometry to two amplitude variables. The axis connection from $E_{+}$ to $E_{-}$ persists without parameter tuning, while rotating the outer connection produces a family of connections from $E_{-}$ back to $E_{+}$. A two-dimensional autonomous amplitude flow cannot by itself produce the proposed saddle-focus return chaos: symmetry reduction hides the phase-dependent reinjection needed for that mechanism.

Generic higher-order terms need not retain this exact rotational symmetry. A normal form can preserve it to a finite truncation order while the full equations have small symmetry-breaking remainders. These unfold the special connection geometry into phase-dependent returns; nearby parameter curves can contain homoclinic orbits to a saddle-focus, and the Shilnikov bifurcation mechanism becomes available. Their exact placement depends on the omitted terms and is not fixed by the displayed truncation alone.

For example, at the perturbed positive-axis saddle, the real unstable eigenvalue is ${\lambda }_u=2a$ and the stable complex pair has real part ${\mu }_1-2a+a^2$. Near the calculated connection curve ${\mu }_1=-a^2/4+o(a^2)$, the contracting ratio isThis is the positive-saddle-value regime, explaining why the relevant saddle-focus return can be expanding.

In the Shilnikov return map, take a positive return coordinate $y$ and $0<\delta <1$, with nondegenerate $A\ne 0,B\ne 0$. At the connection parameter $\mu =0$, the fixed point condition isThe right side tends to zero, whereas the phase passes through infinitely many oscillations as $y\downarrow 0$. Hence there are infinitely many positive fixed points accumulating at zero, with logarithmically spaced amplitudes. For these points,whose magnitude tends to infinity along the accumulating fixed-point sequence. They represent long-period saddle periodic orbits of the full return, with period growing like ${\lambda }_u^{-1}\log (1/y)$, rather than infinitely many stable cycles.

Successive turning branches of this oscillatory return can stretch and fold across suitable intervals. The full saddle-focus Poincaré return map then admits invariant Smale horseshoes and chaotic itineraries near a generic homoclinic orbit at a Shilnikov bifurcation. Varying $\mu$ unfolds accumulating saddle-node and period-doubling events. Whether any resulting invariant set attracts trajectories depends on the other return direction and the omitted terms; the leading scalar map alone does not prove a chaotic attractor for every perturbation. Its domain also excludes negative iterates, which leave the chosen return branch. The degeneracies $A=0$ or $B=0$ would remove the oscillatory mechanism, so the conclusion uses generic, genuinely saddle-focus reinjection.

The contraction mapping theorem states: if $(X,d)$ is a nonempty complete metric space and $f:X\to X$ satisfies $d(fx,fy)\le qd(x,y)$ for a fixed $0\le q<1$, then it has exactly one fixed point and the iteration from every starting point converges to it.

To prove existence, choose $x_0$ and let $x_{n+1}=f(x_n)$. Induction gives $d(x_{n+1},x_n)\le q^{n}d(x_1,x_0)$. The triangle inequality and geometric sum then imply, for $m>n$,The iterates are a Cauchy sequence; completeness supplies their limit $x_{\ast }$. The contraction is continuous, so $f(x_{\ast })=\lim f(x_n)=\lim x_{n+1}=x_{\ast }$. If $y_{\ast }$ is another fixed point, $d(x_{\ast },y_{\ast })\le qd(x_{\ast },y_{\ast })$ forces equality of the points. The same estimate gives geometric convergence from every initial point.

For the particular map, the original PDF gives $X=[\sqrt{a/2},\infty )$, which is a closed, complete subset of the real line. For positive $x$, the arithmetic-geometric mean inequality gives $f(x)=(x+a/x)/2\ge \sqrt{a}$, proving $f(X)\subset X$. FurthermoreThe mean value theorem therefore gives $|f(x)-f(y)|\le |x-y|/2$, a genuine contraction. The fixed-point equation is $x^2=a$, and only the positive root belongs to $X$, soThe printed square-root endpoint is essential: the TeX transcription loses it.

An adjunction $L⊣U$ induces the monad $T=UL$ with unit $\eta$ and multiplication $U\epsilon L$. The Eilenberg-Moore comparison functor sends $B$ to $(UB,U{\epsilon }_B)$. A monadic adjunction is one for which that comparison is an equivalence of categories. Iterating comparisons, when their left adjoints exist, gives a tower; its monadic length is the number of steps until equivalence, with length zero for an equivalence and length one for a monadic adjunction that is not already an equivalence.

The precise monadicity theorem says that a right adjoint is monadic exactly when it reflects isomorphisms and creates coequalizers of pairs whose images admit split coequalizers.

Write $U_r:{\mathcal{C}}_r\to {\mathcal{C}}_{r-1}$ for one-step forgetting. For $r=1$, the free extension of nested partial unary operations is the familiar set $A\times \mathbb{N}$, with ${\alpha }_1(a,j)=(a,j+1)$ and unit $a\mapsto (a,0)$. A map into a set with endomorphism ${\beta }_1$ extends uniquely by $(a,j)\mapsto {\beta }_1^{j}f(a)$.

For $r\ge 2$, letKeep $A$ and all its old operations, and adjoin the disjoint set $S(A)\times \mathbb{N}$. On the new points define ${\alpha }_1(a,j)=(a,j+1)$; every old operation ${\alpha }_i$ with $i\ge 2$ is undefined there. Define the new top operation only at the old eligible points, by ${\alpha }_r(a)=(a,0)$ for $a\in S(A)$.

These domains satisfy the defining rules of the nested partial unary operation category: the new points have no ${\alpha }_1$-fixed points, so none has any higher operation defined. Given $f:A\to U_{r}B$, extend it byThe expression is defined because $f$ sends eligible points to eligible points and ${\beta }_1$ is total. It preserves every defined operation; preservation of ${\alpha }_r$ at old points determines $j=0$, and preservation of ${\alpha }_1$ determines the whole chain. Thus it is the unique extension, proving $L_r⊣U_r$.

For monadicity, $U_r$ reflects isomorphisms: an isomorphism of the lower structures takes $S(A)$ bijectively to $S(B)$, so the inverse of a map preserving the new operation preserves it too.

Next prove the split-coequalizer lifting for a fixed-point-domain operation. Let $f,g:A\rightrightarrows B$ be ${\mathcal{C}}_r$-morphisms and suppose their lower images have a split coequalizer $q:U_{r}B\to Q$, with $s:Q\to U_{r}B$, $t:U_{r}B\to U_{r}A$ satisfyingEquip $Q$ with ${\beta }_Q(z)=q{\beta }_B(sz)$ on its required eligible domain. The lower-structure maps $s,t$ preserve eligibility. For eligible $x\in B$,so $q$ preserves the top operation. If $h:B\to D$ coequalizes $f,g$, its lower factor $\overline{h}:Q\to U_{r}D$ satisfiesThus it is a morphism of the full structures. Uniqueness of the lifted operation follows from the eligible section $s$ and surjectivity of $q$. This creates the specified coequalizer, so the Beck monadicity theorem proves each one-step adjunction is monadic.

Finally, every free top operation just constructed has no fixed points: for $r=1$ it shifts the free chain, while for $r\ge 2$ it sends old eligible points to new distinct points and is undefined on the new points. Freely adjoining any further operation therefore changes nothing, since its required domain is empty. For $n>m$, the monad induced on ${\mathcal{C}}_m$ by the long free-and-forgetful adjunction is consequently the same monad, including unit and multiplication, as that for ${\mathcal{C}}_{m+1}\to {\mathcal{C}}_m$.

The first comparison category is thus equivalent to ${\mathcal{C}}_{m+1}$. Its algebra action evaluates the new formal value by ${\alpha }_{m+1}$, so the comparison from ${\mathcal{C}}_n$ retains precisely that next operation and is the remaining forgetful functor. Iteration recovers one additional operation at each step. These remaining comparisons are not equivalences while an operation is still forgotten: take a two-element set with all retained operations equal to the identity, and choose the next operation to be the identity in one extension and the transposition in another. Higher operations can be identities in the first case and undefined in the second. The identity on the retained structure is not a morphism between these extensions, so the forgetful functor is not full. This includes ${\mathcal{C}}_1\to \mathbf{Set}$.

The monadic tower for nested partial unary operations therefore has exactly the claimed length:

Let a scalar fixed point have magnitude $r>0$. If $z=-r<0$, then $x$ and $y$ have the same sign; the two lifts are the fixed points $(r,1)$ and $(-r,-1)$. If $z=r>0$, their signs are opposite, and the lift is the period-two orbit $(r,-1)\leftrightarrow (-r,1)$. The scalar fixed-point multiplier isThe two-return fixed-point multiplier of the period-two lift is $m^2$, so in either case dynamical stability requires $m<1$. The negative fixed point represents one of two mirror-related single-lobe periodic orbits in the flow, whereas the positive fixed point represents a double-lobe periodic orbit.

For $\delta >1$, the nearby fixed point is $z=\mu +o(\mu )$ and its fixed-point multiplier tends to zero. Thus two stable single-lobe limit cycles for $\mu <0$ glue into one stable double-lobe limit cycle for $\mu >0$. Their periods diverge as the separatrix is approached; $\mu =0,z=0$ is not a finite-period limit cycle.

For $0<\delta <1$, the small branch instead satisfies $z\sim -\mathrm{sgn}(\mu )(|\mu |/A{)}^{1/\delta }$ and its fixed-point multiplier tends to infinity. Hence one unstable double-lobe limit cycle on $\mu <0$ splits into two unstable single-lobe limit cycles on $\mu >0$. Additional finite-size stable branches depend on the surrounding map; near the resonance $\delta =1$ with $A<1$ they are resolved by the saddle-node bifurcations in part (c).

For a fixed point write $\mu =h(z)=z-A\mathrm{sgn}(z)|z{|}^{\delta }$. A saddle-node bifurcation requires both $f(z)=z$ and $f^{\prime }(z)=1$. For $\delta \ne 1$ these conditions giveThe second derivative is nonzero at either nonzero saddle-node bifurcation, making it an ordinary saddle-node bifurcation. For $\delta <1$, defineThe saddle-node bifurcations occur at $(\mu ,z)=(-C,r)$ and $(C,-r)$, and three scalar fixed points coexist for $|\mu |<C$. The two outer branches have fixed-point multiplier below one and are stable, while the middle branch is unstable. Negative points lift to two symmetry-related limit cycles, so these scalar branch counts need not equal the number of flow limit cycles. The cusp records coexistence of distinct return-map branches; even outside it a negative branch represents a symmetry-related pair of flow cycles.

If $\delta =1-\eta$, $\eta \downarrow 0$, thenThis exponentially narrow region ends in a gluing-map cusp at $(\mu ,\delta )=(0,1)$. At its boundaries the graph of $f$ is tangent to the diagonal at $z=+r$ or $-r$, as shown below. At exactly $\delta =1$, $f=\mu +Az$ and $A<1$ gives one stable fixed point for nonzero $\mu$.

For $\delta =1+\eta$, the formal saddle-node bifurcation radius is $(A\delta {)}^{-1/\eta }\to \infty$, and the corresponding saddle-node bifurcation parameters also diverge. Thus the near-origin gluing-map cusp lies on the $\delta <1$ side; the formal large saddle-node bifurcations of the power-law formula on the other side are outside the local gluing-map approximation. They must not be drawn as a second small gluing-map cusp near the origin.

For $\delta =1+\epsilon >1$, the graph of $h(z)$ has a positive maximum at $z=r$ and a negative minimum at $z=-r$. The central branch $|z|<r$ is stable, while the formal outer branches $|z|>r$ are unstable. Near the gluing point only the stable central branch is relevant: it crosses from negative $z$ to positive $z$, representing the stable two-to-one gluing described in part (b). The formal full power-law diagram has three fixed points between its two large saddle-node bifurcations.

For $\delta =1+\epsilon <1$, the positive branch has its minimum at $(\mu ,z)=(-C,r)$, while the negative branch has its maximum at $(C,-r)$. The inner pieces $|z|<r$ are unstable and outer pieces stable. In the gluing-map cusp interval there is a stable negative branch, a stable positive branch, and an unstable central branch; the unstable branch passes through the singular gluing point with opposite sign to $\mu$. Outside the gluing-map cusp just one stable outer branch remains. The lower panels of the figure plot these diagrams in fold-scaled variables; dashed lines denote instability and solid lines dynamical stability. Neither diagram has a flip bifurcation, since the scalar derivative is positive wherever defined.