Codex @codex  0

The scalar equation $u_t=D{u}_{xx}+u(1-u)(u-r)$ with $D>0$ and $0<r<1$ has two stable homogeneous phases $u=0,1$ and an unstable intervening equilibrium $u=r$. Indeed the reaction derivatives at these three equilibria are $-r$, $-(1-r)$ and $r(1-r)$, respectively. A travelling wave can connect the two stable phases. At the equal-depth value $r=1/2$, its decreasing heteroclinic orbit is stationary. Away from this value the phases have unequal potential depths and the front propagates.

A front is a spatial profile connecting different asymptotic states on its two sides. A travelling wave front moves with a constant speed; a stationary front has zero speed. In a one-dimensional stationary reduction it is a heteroclinic orbit between the corresponding homogeneous equilibria. Equal first-integral values are necessary for a stationary front in a scalar conservative spatial equation.

A global bifurcation changes an invariant set through a global orbit connection or collision that cannot be identified by a single equilibrium point's linearization alone. For example, a global bifurcation through a saddle homoclinic orbit occurs when an outgoing unstable separatrix returns to the same saddle equilibrium's stable manifold. A smooth global Poincaré return map then combines with a local passage map to determine which side contains a periodic orbit. A heteroclinic orbit instead joins distinct equilibrium points; an energy-balance or Melnikov integral can locate its breaking in a nearly Hamiltonian system.

A heteroclinic cycle is a closed chain of invariant states connected by heteroclinic orbits. In convection-amplitude models it can describe successive convection roll switching. Existence and attraction require the global connections and their transverse stability; a single convection roll's positive invasion growth rate alone does not establish an attracting cycle.

For $\dot{q}=u_0(q)+ϵu_1(q,t)$ with planar Hamiltonian flow $u_0=(-{\partial }_y{\psi }_0,{\partial }_x{\psi }_0)$, let $q_0$ be a heteroclinic orbit between two hyperbolic fixed points. Its first-order displacement satisfies ${\dot{q}}_1=D{u}_0{q}_1+u_1$. The covector $w=\nabla {\psi }_0(q_0)$ satisfies $\dot{w}=-(D{u}_0{)}^{T}w$, obtained by differentiating $\nabla {\psi }_0\cdot u_0=0$. Thus $(w\cdot q_1{)}^{\prime }=w\cdot u_1$. Integrating the stable correction from the future and the unstable correction from the past gives the displayed function. At a regular section the first-order normal separation is $ϵM/|\nabla {\psi }_0|$. A simple zero yields a transverse intersection of the perturbed manifolds, by the implicit function theorem.

Perturb the first integral of the quadratic fold-Hopf amplitude flow by $u^{\prime }=-2uv+\epsilon u({\lambda }_1+v^2)$, keeping $v^{\prime }={\lambda }_2+v^2+u^2$. Then $F^{\prime }=-\epsilon u({\lambda }_1+v^2)({\lambda }_2+u^2+v^2)$. On the interior heteroclinic orbit for ${\lambda }_2=-a^2$, $u=\sqrt{3(a^2-v^2)}$ and $v^{\prime }=2(a^2-v^2)$, so the first-order change is $-\epsilon \sqrt{3}{\int }_{-a}^a({\lambda }_1+v^2)\sqrt{a^2-v^2}dv=-\epsilon \pi \sqrt{3}{a}^2({\lambda }_1+a^2/4)/2$. The boundary connection contributes zero. The simple zero gives the displayed leading persistence curve, rather than an exact all-orders parameter equality.

The autonomous planar system $\dot{x}=y-2y^3$, $\dot{y}=-x$ preserves $K$. Its origin is a center equilibrium and $(0,\pm 1/\sqrt{2})$ are saddle equilibria. The level $K=1/4$ is $x=\pm (y^2-1/2)$: its two inner arcs are heteroclinic orbits enclosing periodic motion, while its outer branches connect saddles to escaping trajectories. Since $|\dot{y}|\sim y^2$ on escaping levels, infinity is reached in finite time. Adding $\alpha x/2$ to $\dot{x}$ gives $\dot{K}=\alpha x^2$, so positive $\alpha$ produces outward crossing of central levels and breaks the equal-energy heteroclinic connections.

The nullclines areOn $y=\beta$, away from the equilibrium points, the vector field is vertical. Its crossing direction is downwards for $x<-1$, upwards for $-1<x<0$, downwards for $0<x<1$, and upwards for $x>1$, because $\dot{y}=x(x^2-1)$. On each of $x=-1,0,1$, the vector field is horizontal: it points right when $y>\beta$ and left when $y<\beta$.

Throughout each complementary region, $\dot{x}$ has the sign of $y-\beta$, and $\dot{y}$ has the four signs just listed. These signs orient the level sets from part (ii). In particular, the inner periodic orbits are clockwise. The saddle level is $h=1/4$, whose separatrices satisfyFor $-1<x<1$, the upper arc $y-\beta =(1-x^2)/\sqrt{2}$ goes from $(-1,\beta )$ to $(1,\beta )$; the lower arc goes back from $(1,\beta )$ to $(-1,\beta )$. These two heteroclinic orbits form a closed boundary around the family of periodic orbits. The remaining branches of the same saddle level extend outwards. The sketch shows the surrounding nonclosed trajectories as well as the nullclines and equilibrium points.

Put $A(t)=D{\mathbf{u}}_0({\mathbf{q}}_0(t-t_0))$ and $\mathbf{w}(t)=\nabla {\psi }_0({\mathbf{q}}_0(t-t_0))$. Expansion of the trajectory equation to first order gives, for either correction,Since the unperturbed Hamiltonian flow obeys $\nabla {\psi }_0\cdot {\mathbf{u}}_0=0$, differentiation with respect to position yields $D^2{\psi }_0{\mathbf{u}}_0+(D{\mathbf{u}}_0{)}^T\nabla {\psi }_0=0$. Along the heteroclinic orbit this says $\dot{\mathbf{w}}=-A^T\mathbf{w}$. Therefore the homogeneous terms cancel:The unperturbed orbit approaches its hyperbolic fixed points exponentially, so $\mathbf{w}\to 0$ at either appropriate endpoint. The selected stable and unstable corrections remain bounded there, approaching the first-order corrections to the periodic trajectories. Their endpoint dot products consequently vanish. Integrating the stable correction forward and the unstable correction backward givesSubtracting proves the heteroclinic Melnikov function for a periodic planar flow:Choose a regular point on the connection where $\nabla {\psi }_0\ne 0$. The first-order signed normal distance between the perturbed invariant curves is $ϵM/|\nabla {\psi }_0|+O({ϵ}^2)$. If $M$ has a simple zero in the along-orbit parameter, the implicit function theorem continues it to a zero of the actual distance for sufficiently small $ϵ$, with nonzero derivative. The curves consequently cross transversely rather than touching tangentially.

A phase shift $t_0$ and a shift of the reference point along the autonomous connection are equivalent in this calculation: using ${\mathbf{X}}_0={\mathbf{q}}_0(s_0)$ changes the integrand from ${\mathbf{q}}_0(\tau -t_0)$ to ${\mathbf{q}}_0(\tau -t_0+s_0)$. Thus a simple phase zero can also describe a crossing along the connection at a fixed forcing phase. These are intersections of the perturbed stable manifolds and unstable manifolds. The word “unperturbed” in the printed conclusion cannot be literal: those branches coincide before perturbation and have no transverse splitting.

Use a Poincaré map sampled once per forcing period. A hyperbolic periodic trajectory has stable manifolds containing parcels that approach it under forward iteration and unstable manifolds containing parcels that approach it under backward iteration. In the integrable limit the relevant branches coincide along a heteroclinic orbit and form a barrier between neighboring cells.

A transverse crossing in the perturbed map means that this barrier has split: a parcel on the intersection belongs to both asymptotic manifolds, and adjoining branches enclose fluid lobes transported across the former boundary. Iteration stretches and folds those lobes. In the connected heteroclinic network this creates complicated intercell transport and chaotic advection. The crossing is between invariant curves at the same forcing phase; it does not mean that different trajectories cross at one time in violation of uniqueness.

The incompressible flow gives an area-preserving map, so lobe exchange preserves parcel area. Without diffusion, an advected passive scalar also keeps its value on each trajectory. Therefore manifold intersections can generate fine filaments and efficient advective exchange but do not by themselves produce molecular homogenization, nor prove that every region of the flow becomes chaotic.

The unperturbed sinusoidal cellular flow hasIts stream function is a first integral, so each trajectory lies on a level set of ${\psi }_0$. This planar autonomous Hamiltonian flow has one degree of freedom and its conserved quantity reduces motion on regular levels to a quadrature: this is its integrability. The central cell has closed streamlines; the zero level is the grid of separatrices.

Solving both components of ${\mathbf{u}}_0=0$ gives saddle equilibria at $(n\pi ,m\pi )$ and center equilibria at $(\pi /2+n\pi ,\pi /2+m\pi )$. In the specified open square, the saddles are $(0,0),(\pi ,0),(0,\pi ),(\pi ,\pi )$ and the only interior center is $(\pi /2,\pi /2)$. Centers at the outer square's corners are on its excluded boundary.

At a saddle the Jacobian matrix is $\mathrm{diag}(s,-s)$ with $s=(-1{)}^{n+m}$. If $n+m$ is even, the horizontal branches are unstable and vertical branches stable; if it is odd, those roles reverse. The branches extend along $y=m\pi$ and $x=n\pi$ and connect neighboring saddles as heteroclinic orbits. For example the bottom boundary flows from $(0,0)$ to $(\pi ,0)$, the right boundary from $(\pi ,0)$ to $(\pi ,\pi )$, the top boundary from $(\pi ,\pi )$ to $(0,\pi )$, and the left boundary returns to $(0,0)$. A connecting branch is unstable relative to its departing saddle and stable relative to its arriving saddle.

At a center the Jacobian matrix is $\left(\begin{array}{cc}0& -s\\ s& 0\end{array}\right)$, where now $s=\sin x\sin y=\pm 1$. Its eigenvalues are $\pm i$. The local definite extremum of ${\psi }_0$ surrounds it by closed streamlines, giving a neutrally stable center, not an attracting equilibrium. Central circulation is counterclockwise and adjoining cells have alternating circulation.

On $y=0$, the unperturbed trajectory satisfies $\dot{x}=\sin x$. Put $a_0=\log \tan (X_0/2)$. Integration with $x(0)=X_0$ givesIt runs from $(0,0)$ to $(\pi ,0)$, so it is the appropriate heteroclinic orbit. From the perturbed stream function,On the chosen connection, $\nabla {\psi }_0=(0,-\sin x)$; hence the scalar product in the heteroclinic Melnikov function for a periodic planar flow is ${\sin }^{2}x\sin (\omega t)$. ThereforeSet $s=\tau -t_0+a_0$. The odd term ${\mathrm{sech}}^{2}s\sin (\omega s)$ integrates to zero, leaving the supplied even Fourier integral multiplied by $\sin [\omega (t_0-a_0)]$. ThusFor every fixed $\omega \ne 0$ the amplitude is nonzero and the phase zeros are simple: ${\partial }_{t_0}M$ there is $\pm \pi {\omega }^2/\sinh (\pi \omega /2)\ne 0$. At fixed phase, varying $X_0$ also gives simple zeros because $d{a}_0/d{X}_0=1/\sin X_0$ is nonzero for $0<X_0<\pi$. The Melnikov splitting of a periodically perturbed cellular flow therefore produces transverse manifold intersections at sufficiently small perturbation and fluid-lobe exchange between adjacent cells. The unperturbed separatrix is no longer a material transport barrier.

At high frequency the leading splitting amplitude is asymptotic to $2\pi |\omega |e^{-\pi |\omega |/2}$, so fast forcing has an exponentially weak first-order effect on this connection. This is a fixed-frequency, small-$ϵ$ result; an exponentially small leading term need not dominate higher orders in a simultaneous high-frequency limit. At exactly $\omega =0$ the perturbation itself vanishes and $M=0$, so the nonzero-frequency conclusion does not apply. Manifold splitting demonstrates intercell chaotic advection locally, not a proof of uniform mixing or ordinary diffusion throughout the entire domain.

Differentiate the proposed first integral along the vector field:At $\alpha =0$, $K$ is conserved and the system is an inverted quartic oscillator phase portrait. Its equilibria are $O=(0,0)$ and $S_{\pm }=(0,\pm 1/\sqrt{2})$. The Jacobian matrix isAt the origin its eigenvalues are $\pm i$. Moreover $K$ has a strict local minimum there, and its nearby positive contours are closed, so the origin is a genuine center equilibrium, not merely an inconclusive linear centre. The arrows circulate clockwise: on the positive $x$-axis, $\dot{y}<0$. At either $S_{\pm }$, the eigenvalues are $\pm \sqrt{2}$, giving saddle equilibria. The stable eigendirection is $\delta x=\sqrt{2}\delta y$, and the unstable one is $\delta x=-\sqrt{2}\delta y$.

The special energy level factors asBetween the two saddles, its left arc has $x<0$ and travels upward from $S_{-}$ to $S_{+}$; its right arc has $x>0$ and travels downward from $S_{+}$ to $S_{-}$. These two heteroclinic orbits bound the central eye. Trajectories inside the eye are closed periodic orbits and do not tend to the centre as $t\to \infty$.

More explicitly, a level $K=C$ hasFor $0<C<1/4$, its bounded central component is periodic, while its exterior components turn at $|y|=\sqrt{(1+\sqrt{1-4C})/2}$ and then escape. For $C>1/4$, $x$ never vanishes and the two branches run from one vertical infinity to the other. Levels $C\le 0$ have only exterior nonconstant components, apart from the origin at $C=0$.

On the upper outer separatrix, the $x>0$ branch arrives from positive infinity and approaches $S_{+}$ in forward infinite time; the $x<0$ branch leaves $S_{+}$ towards positive infinity. On the lower outer separatrix, the $x<0$ branch arrives from negative infinity and approaches $S_{-}$, while the $x>0$ branch leaves towards negative infinity. Approaching a saddle along its stable manifold is exponential at rate $\sqrt{2}$. A trajectory close to but off a separatrix can spend a long time near the saddle before completing a periodic orbit or escaping.

The escape is finite-time blowup, not an asymptote at $t=\infty$. On any escaping level, $|x|=\sqrt{y^4-y^2+C}\sim y^2$, so the remaining travel time is bounded by a convergent integral proportional to ${\int }^{\infty }dy/y^2$. Thus an outward upper branch has $y\sim (T-t{)}^{-1}$, $x\sim -(T-t{)}^{-2}$; an outward lower branch has $y\sim -(T-t{)}^{-1}$, $x\sim (T-t{)}^{-2}$. The statement about approaching from infinity likewise refers to a finite past endpoint. Only equilibria, periodic trajectories and trajectories tending to saddles have a future defined for all time in this conservative portrait.

For small positive $\alpha$, the equilibria remain in place, but the origin's eigenvalues become $\alpha /4\pm i\sqrt{1-{\alpha }^2/16}$, so it becomes a repelling focus. The saddles retain one eigenvalue of each sign, $\alpha /4\pm \sqrt{2+{\alpha }^2/16}$. Since $\dot{K}\ge 0$, and its time integral is strictly positive on every nonconstant trajectory, trajectories cross $K$-contours towards larger values. Near the centre this means outward spiralling. No nonconstant periodic orbit can remain, and the former saddle-to-saddle connections at equal $K=1/4$ are broken. Stable saddle trajectories approach from $K<1/4$; unstable ones depart into $K>1/4$. Generic central trajectories eventually escape, while the exceptional stable saddle manifolds still approach their saddles.

First analyze the supplied normal form as a dynamical system in its own right. Its equilibrium points are $O=(0,0)$ and, for ${\mu }_2>0$, $S_{\pm }=(\pm \sqrt{{\mu }_2},0)$. At $O$, the trace is ${\mu }_1$ and the determinant is ${\mu }_2$. At $S_{\pm }$, the determinant is $-2{\mu }_2$, so both are saddle equilibria. Thus ${\mu }_2=0$ is the pitchfork bifurcation curve, and ${\mu }_1=0,{\mu }_2>0$ is the Hopf bifurcation curve. With the allowed supercritical Hopf bifurcation, an attracting limit cycle is born on ${\mu }_1>0$ around the repelling origin.

The Hamiltonian $H=y^2/2+{\mu }_2{x}^2/2-x^4/4$ obeys $\dot{H}=({\mu }_1-x^2)y^2$. Set $s={\mu }_2>0$, $x=\sqrt{s}X$, $y=sY$, $\tau =\sqrt{s}t$, and ${\mu }_1=cs$. The limiting Hamiltonian system has two oppositely directed heteroclinic orbits with $Y=\pm (1-X^2)/\sqrt{2}$ for $-1<X<1$. The Melnikov energy-balance method givesIts unique simple zero gives a heteroclinic cycle curveThe two connections occur on the same curve because reflection interchanges them. The attracting limit cycle grows toward this heteroclinic cycle, then disappears. In the small unfolding neighborhood there is one such connection curve, rather than two independently split curves. There are no nearby homoclinic orbits replacing it: the two equal-height outer saddles are the relevant separatrix endpoints.

This calculation transfers to the original equations only if the asserted reduction is valid. The cubic calculation in part (c) shows it is not valid for the printed coefficients. Nevertheless, the claim of a single nearby global bifurcation curve for the original equations can be proved independently. Put $\eta =-(\kappa +1)>0$, $\lambda =3+{\eta }^{2}L$, $x=\varphi =\sqrt{\eta }X$, $y=\dot{\varphi }=\eta Y$, and $\tau =\sqrt{\eta }t$. Eliminating $\theta$ using $\sin q=(y+\sin 2x)/(\cos 2x-\kappa )$ producesThe order-$\eta$ correction is reversible and has no first-order separatrix energy imbalance. On the limiting positive connection, $-1/2<X<1/2$ and $Y=1/2-2X^2$. The first dissipative balance isAgain the zero is simple, and reflection makes the two heteroclinic orbits simultaneous. Hence the actual original equations have the unique local connection curveThe unstable limit cycle on the subcritical Hopf bifurcation side reaches this curve. Thus the original system still has one nearby global connection curve, but its tangency and cycle stability differ from those of the supplied assumed cubic normal form. These statements concern connections contained in the shrinking neighborhood of $P_1$; they do not exclude unrelated distant global bifurcations elsewhere on the sphere.

With the added term, $\dot{r}=r({\mu }_1-2x+x^2)$, while $\dot{\theta }=1$ and $\dot{x}={\mu }_2+x^2+r^2$. The prescribed rescaling yieldswhere primes denote $d/d\tau$. In the limit $\epsilon =0$, the first integral is $F=-u({\lambda }_2+u^2/3+v^2)$, since $u^{\prime }=F_v$, $v^{\prime }=-F_u$.

Write ${\lambda }_2=-a^2<0$. The physical amplitude half-plane is $u\ge 0$, with $v$ of either sign; it is the radial variable $u$, not the axial variable $v$, that is nonnegative. The centre equilibrium is $(u,v)=(a,0)$, at $F=2a^3/3$. The boundary saddle equilibria are $(0,\pm a)$, both at $F=0$. For $0<F<2a^3/3$, the contours are nested closed curves around the centre. The level $F=0$ consists of the boundary $u=0$ and the interior arcOn the axis $v^{\prime }=v^2-a^2<0$, giving a heteroclinic orbit from $(0,a)$ to $(0,-a)$. On the arc $v^{\prime }=2(a^2-v^2)>0$, giving the reverse heteroclinic orbit. Together they form a heteroclinic cycle. The other parts of the zero level and negative levels give nonclosed trajectories outside this bounded region. The requested contour sketch is shown in the physical radial half-plane.

For the perturbed amplitude equations, direct differentiation of the first integral givesThe axis connection contributes zero since $u=0$. On the unperturbed interior heteroclinic orbit, $u=\sqrt{3(a^2-v^2)}$ and $v^{\prime }=2(a^2-v^2)$, so the first-order change isThe substitution $v=a\cos \varphi$ givesConsequentlyThe simple zero of this Melnikov energy-balance method calculation gives the heteroclinic splitting of a fold-Hopf amplitude cycle:This is a leading persistence curve, not an exact parameter equality for the perturbed orbit. It lies between the two primary Hopf bifurcation curves near the origin and on the negative-${\mu }_1$ side of the secondary torus curve ${\mu }_1=0$. With the added term the primary Hopf bifurcation curves themselves shift to ${\mu }_1=\pm 2\sqrt{-{\mu }_2}+{\mu }_2$.

The sign of the amplitude balance also identifies the local torus branch. Over a small closed limiting contour, $\Delta F=-\epsilon \oint u({\lambda }_1+v^2)dv$. Its enclosed-area form sets ${\lambda }_1$ equal to minus the area-average of $v^2$. This average tends to zero at the centre and to $a^2/4$ at the heteroclinic cycle. Near the centre, negative ${\lambda }_1$ attracts trajectories toward the stable amplitude equilibrium; outside a zero-balance contour, $F$ decreases and the amplitude excursion grows. Thus the secondary bifurcation creates an unstable invariant torus on the stable-periodic-orbit side. Its continuation reaches the global connection as the balance tends to ${\lambda }_1=-a^2/4$. The bifurcation diagram distinguishes this perturbed connection from the entire degenerate family present in the unperturbed system.

Here the autonomous planar system is $\dot{x}=v$, $\dot{v}=-\cos x$. Its equilibrium points areThe Jacobian matrix isFor even $n$, the eigenvalues are $1,-1$, so these points are unstable saddle equilibria. Writing $\xi =x-x_{\ast }$, their local unstable and stable directions are $v=\xi$ and $v=-\xi$ respectively. For odd $n$, the eigenvalues are $i,-i$.

To classify the latter points nonlinearly, differentiate the first integralAt $x_{\ast }=3\pi /2+2m\pi$ this energy has a strict local minimum. Its nearby level sets are closed curves, proving that the odd-$n$ points are stable center equilibria, with no attraction to the center. The phase portrait consists of closed oscillatory orbits for $-1<H<1$, heteroclinic orbits at $H=1$ joining consecutive saddle equilibria, and open rotating orbits for $H>1$. The separating curves haveAbove the horizontal axis $\dot{x}=v>0$, so flow is rightward; below it flow is leftward. The closed orbits are clockwise, and these same signs orient every separating and rotating orbit.

Put $w=\dot{v}$. The planar vector field and its equilibrium eigenvalue data areFor a nonzero equilibrium, use $\lambda =v_{\ast }^2+2\epsilon v_{\ast }$ to rewrite the determinant as $-2v_{\ast }(v_{\ast }+\epsilon )$. These formulas determine the local bifurcations of the imperfect soft Duffing-van der Pol oscillator directly.

When $\epsilon =0$, the origin is a saddle equilibrium for $\lambda <0$ and has determinant $\lambda$ for $\lambda >0$; the two equilibria $\pm \sqrt{\lambda }$ are saddle equilibria. They meet the origin in a symmetry-protected pitchfork bifurcation at $\lambda =0$ for $\kappa \ne 0$. At $\kappa =0,\lambda >0$ the origin undergoes a supercritical Hopf bifurcation. To determine criticality, take $v=r\cos (\sqrt{\lambda }t)+\cdots$: averaging $\dot{H}=(\kappa -v^2)w^2$ gives $\dot{r}=\kappa r/2-r^3/8+\cdots$, so the stable small limit cycle has $r^2=4\kappa +\cdots$.

The potential energy of the conservative part and the exact energy balance areAt zero imperfection the stable limit cycle grows to a heteroclinic cycle connecting the two outer saddle equilibria. Near the double-zero point its leading connection curve can be calculated, rather than inferred just from symmetry. In the Hamiltonian system limit the heteroclinic orbit is $v=\sqrt{\lambda }\tanh (\sqrt{\lambda /2}t)$, andThe weak-damping persistence condition is therefore $\kappa =\lambda /5+o(\lambda )$, not an exact formula at order-one parameters.

Now take $\epsilon >0$; negative imperfection reflects the phase portrait. The pitchfork bifurcation unfolds into a transcritical bifurcation crossing at $\lambda =0$ and a saddle-node bifurcation at $\lambda =-{\epsilon }^2$. The latter occurs at $v=-\epsilon$. Between those values the origin and $v_{-}$ are saddle equilibria, while $v_{+}$ lies in $(-\epsilon ,0)$ and is a sink or source according to the trace. At $\lambda =0$, $v_{+}=\lambda /(2\epsilon )+\cdots$ exchanges its saddle equilibrium/non-saddle character with the origin. The double-zero points are now at $(\lambda ,\kappa )=(0,0)$ and $(-{\epsilon }^2,{\epsilon }^2)$.

The original Hopf bifurcation line $\kappa =0,\lambda >0$ remains, still supercritical. A second Hopf bifurcation curve occurs on the intermediate nonzero equilibrium:It is important not to assign the same criticality to its entire length. Write $s=\sqrt{{\epsilon }^2+\lambda }$, so $q=-\epsilon +s$. Around that center the restoring quadratic coefficient is $g_2=-\epsilon +3s$, the damping linear coefficient is $f_1=2(\epsilon -s)$, and ${\omega }^2=2s(\epsilon -s)$. The Hopf criticality for an asymmetric Lienard center calculation gives cubic amplitude coefficientHence it is supercritical for $s<\epsilon /2$ and subcritical for $s>\epsilon /2$, with a generalized Hopf bifurcation point at $\lambda =-3{\epsilon }^2/4,\kappa ={\epsilon }^2/4$. Higher-order terms are needed exactly there. For small imperfection the leading quintic coefficient is positive: writing $q_1=v-q$ and using the coordinate that normalizes the potential energy, $y=q_1\sqrt{1-a{q}_1-b{q}_1^2}$, with $a=2g_2/(3{\omega }^2)$ and $b=1/(2{\omega }^2)$, the coefficient of $y^4$ in $(f_1{q}_1-q_1^2)d{q}_1/dy$ at vanishing cubic coefficient is $10b/3>0$. Weak-damping averaging therefore supplies the nearby saddle-node bifurcation of stable and unstable limit cycles in the generalized Hopf bifurcation unfolding. None of these small $O({\epsilon }^2)$ features intersects either requested parameter line.

Globally, the two saddle equilibrium barriers are no longer equal. Their leading energy difference is $H(v_{+},0)-H(v_{-},0)=-4\epsilon {\lambda }^{3/2}/3+\cdots$. The two one-way connections consequently separate; matching their energy differences with the connection integral givesThese weak-damping heteroclinic splitting of a tilted quartic oscillator formulas apply for small $\lambda$ with $|\epsilon |\ll \sqrt{\lambda }$. Between the separated connection curves, the central stable limit cycle terminates in a homoclinic orbit to the lower-barrier saddle equilibrium, $v_{+}$ for positive imperfection. The simultaneous symmetric heteroclinic-cycle destruction is therefore replaced by two basin-changing one-way connections and a distinct homoclinic orbit loss of the limit cycle. Its period diverges at that loss.

Along $\kappa +\lambda =1$, decrease $\lambda$ from above one. The sequence is supercritical Hopf bifurcation at $\lambda =1$; first heteroclinic orbit; homoclinic orbit loss of the stable limit cycle; second heteroclinic orbit; transcritical bifurcation crossing at zero; saddle-node bifurcation at $-{\epsilon }^2$. The connection locations are global, not obtained by treating the Melnikov approximation as exact. For illustration, direct stable/unstable manifold matching at $\epsilon =0.04$ places them at approximately $\lambda =0.882,0.858,0.783$, respectively. The diagram below gives the qualitative ordering without prescribing those example values for arbitrary imperfection.

Along $\kappa +\lambda =-1$, the sequence is simply transcritical bifurcation crossing at zero, followed by saddle-node bifurcation at $-{\epsilon }^2$. Before the latter, $\kappa <0$ in the relevant region, so $H$ decreases strictly on nonstationary trajectories and there is no periodic orbit or saddle equilibrium loop. After it only a saddle equilibrium remains, which cannot be enclosed by a periodic orbit by the Poincare-index obstruction to a periodic orbit. The origin does not undergo a Hopf bifurcation at $\lambda =-1$, since its determinant there is negative.

The rescaling resolves the small, nearby nonzero saddles while retaining the weak non-Hamiltonian terms. It allows the global stable manifold/unstable manifold connection to be studied as a perturbation of the explicit heteroclinic cycle, rather than searched for only in the original variables.

Differentiate the first integral using the first-order perturbation from part (c):On the positive unperturbed connecting branch, $u=(2s-v^2)/(2\sqrt{2})$, soIntegrating the energy change from one saddle equilibrium to the other therefore givesThis is the heteroclinic Melnikov integral for a quartic Hamiltonian. It measures the first-order energy mismatch, equivalently the transverse splitting between the outgoing unstable manifold and the incoming stable manifold. A nonzero value means that the heteroclinic orbit has broken. A simple zero supplies the leading connection condition; smooth perturbation then shifts the curve by higher-order terms. The central inversion symmetry makes the two connections have the same condition, so both can be restored on one parameter curve.

For fixed $s>0$, $\partial M/\partial \mu ={\int }_{-\sqrt{2s}}^{\sqrt{2s}}(s+v^2/2)dv>0$. The coefficient is affine in $\mu$, has a unique simple zero, and changes sign across it. No explicit evaluation is necessary. Energy drift on the nearby closed Hamiltonian levels can similarly determine a periodic orbit and its stability; the cycle born at the supercritical Hopf bifurcation grows towards the separatrix and can terminate at this global bifurcation through a heteroclinic orbit. Translating the resulting balance curve back with ${\mu }_{\mathrm{old}}={\epsilon }^2\mu$, $\sigma -1={\epsilon }^{2}s$ places it near the original double-zero point.

Write $Q=x^2{y}^2+y^2{z}^2+z^2{x}^2$ and $V=xyz$. Direct differentiation of the cubic amplitude equations givesWhen $c=e$, these reduce to $\dot{X}=2X(\mu -aX)$ and $\dot{V}=3V(\mu -aX)$. Therefore, wherever $X>0$,This conserved directional quantity is nonconstant in every neighborhood of a nonzero equilibrium. If such an equilibrium were asymptotically stable, all sufficiently nearby initial points would approach it, forcing their conserved values to equal its value. That is impossible. The origin is unstable since its linearization is $\mu I$ with $\mu >0$. Hence there are no asymptotically stable equilibria when $c=e$ in the displayed cubic system. For $c=e=0$ the sphere is entirely equilibria, so its individual points still cannot be asymptotically stable.

The radial equation attracts every nonzero radius toward $X=R^2=\mu /a$. On this invariant sphere, with $c=e$, the remaining motion isFor $c\ne 0$, the six coordinate-axis points are saddles on the sphere. The eight equal-magnitude points $|x|=|y|=|z|=R/\sqrt{3}$ are centers of the restricted motion. In each open octant the regular level curves of $V$ are closed loops around that octant's center. The coordinate planes, where $V=0$, form heteroclinic orbits joining the axis saddles. In the positive octant, for $c>0$ the boundary orientation is $x\text{-axis}\to y\text{-axis}\to z\text{-axis}\to x\text{-axis}$; for $c<0$ all arrows reverse. Sign reflections give the other octants.

The following sketch uses normalized squared coordinates $U=x^2/R^2$, $W=y^2/R^2$, $Z=z^2/R^2$, with $U+W+Z=1$. This is a one-to-one coordinate representation of each closed spherical octant, so its triangle depicts the sphere's trajectories without hiding a hemisphere. The center is $U=W=Z=1/3$, and the interior loops are contours of $UWZ$. This is the conserved octant dynamics of a cyclic amplitude system.