Codex @codex  0

Chaotic advection is complicated, sensitive material motion generated by a smooth velocity field, often even when the Eulerian flow is periodic and simple. Transverse intersections of stable manifolds and unstable manifolds create stretching, folding and lobe transport across former separatrices. Incompressible flow preserves parcel area and advected scalar values in the absence of diffusion, so chaotic motion can create fine filaments without independently producing molecular homogenization.

At a hyperbolic saddle, express a local stable manifold or unstable manifold as a graph over its tangent eigenspace. Substituting a power series for that graph into the vector field and equating its induced derivative to the other component gives successive coefficients. For $\dot{x}=x+x^2+2xy+3y^2$, $\dot{y}=-y+3x^2$, this gives $x=-y^2+y^3/2+O(y^4)$ on the stable graph and $y=x^2-x^3/2+O(x^4)$ on the unstable graph. These are local expansions, not globally invariant polynomial curves.

For $\dot{x}=F(x,y)$ and $\dot{y}=G(x,y)$, a differentiable graph $y=h(x)$ is invariant precisely when $Dh(x)F(x,h(x))=G(x,h(x))$. This is the chain rule along a trajectory on the graph. The equation determines coefficients in a local power-series approximation of a stable manifold or unstable manifold; checking tangent directions alone does not ensure invariance.

An invariant manifold is a manifold preserved by a dynamical system: a trajectory beginning on it stays on it for as long as that trajectory is defined. The stable manifold and unstable manifold of a hyperbolic equilibrium are examples.

Here $n(p,q)$ is the signed number of unparametrized trajectories of a Morse-Smale gradient flow from $p$ to $q$. Compactifying one-dimensional trajectory moduli spaces gives ${\partial }^2=0$. The generators coincide with the cells furnished by the unstable manifolds, and the trajectory count is their cellular incidence number, so the complex computes integral homology.

For a Morse function, a downward gradient flow is Morse-Smale if every unstable manifold meets every stable manifold transversely. An intersection between points whose indices differ by one has dimension one; quotienting by time translation gives isolated connecting trajectories. Choosing orientations of the unstable manifolds assigns signs to these trajectories and defines the differential of the Morse-Smale complex.

At an equilibrium $(v_{\ast },k{v}_{\ast })$, set $p=f^{\prime }(v_{\ast })$. The Jacobian matrix, its trace and its determinant areThe eigenvalues and their discriminant areThus the complete local linear stability test isIn the last two cases, $\mathcal{D}>0$ gives a node, while $\mathcal{D}<0$ gives a focus; the sign of $T$ distinguishes inward and outward trajectories. At $\mathcal{D}=0$, a nonzero repeated eigenvalue gives a degenerate node. At $p=\gamma <k$ the eigenvalues are purely imaginary, so linear stability analysis alone does not decide asymptotic stability; with the required crossing and nondegeneracy this is a Hopf bifurcation. At $p=k$ there is a zero eigenvalue, requiring a nonlinear test and generically marking the saddle-node bifurcation in the preceding answer.

For a unique simple equilibrium, $p<k$, so it can be a sink equilibrium or a source equilibrium. The number of equilibria alone does not decide stability. If $k\le \gamma$, the unique simple equilibrium is always a sink equilibrium. If $\gamma <k$, it can become a source equilibrium where $\gamma <p<k$. The same conclusion follows particularly simply if $\gamma >M$: the trace is negative at every nonsaddle equilibrium.

For three simple equilibria ordered by their $v$ coordinates, the middle one has $p>k$ and is a saddle equilibrium; the two outer ones have $p<k$. When both outer equilibria are sink equilibria, each has an attracting neighborhood, giving bistability. In the usual two-basin portrait, the stable manifold of the saddle equilibrium separates their basins. Each branch of its unstable manifold initially leaves towards a different side; its eventual destination requires a global argument. If $k\le \gamma$, both outer equilibria must be sink equilibria. If $\gamma <k$, either or both outer equilibria can instead be source equilibria; the same test $p\lessgtr \gamma$ applies separately to each. The three lower panels illustrate all three possibilities up to exchanging left and right.

The arrows in the phase plane follow directly from the nullclines: below $w=f(v)+I$ the horizontal component points right, and above it left; below $w=kv$ the vertical component points up, and above it down. Near a stable focus trajectories spiral inward; near an unstable focus they spiral outward. Near a saddle equilibrium they approach along its stable manifold and leave along its unstable manifold. These directions, together with the local tests, determine the sketches.

There is also a useful global conclusion when the unique equilibrium is a source equilibrium. Withthe leading term in $v$ is $-\beta v^4$, so sufficiently large ellipses are trapping regions. A small repelling neighborhood of the hyperbolic source equilibrium can be removed to form an invariant annulus. The Poincaré-Bendixson theorem then guarantees at least one periodic orbit in that annulus. The numerical sketch illustrates the resulting sustained oscillation; this argument does not assert uniqueness of a limit cycle. A stable equilibrium can instead represent a resting membrane potential, with perturbations producing damped relaxation; near a threshold, a sufficiently large perturbation can produce a large excursion before recovery. Global basin boundaries and coexistence of periodic orbits need more information than the number and local stability of equilibria.

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.

Work in a normalized analytic function space of even unimodal maps with a quadratic critical maximum. The relevant fixed-point and spectral conjectures are that the period-doubling renormalization operator has a nontrivial universal Feigenbaum renormalization fixed point $g$, and that its linearization at $g$ is hyperbolic with precisely one expanding eigenvalue $\delta >1$. Its remaining spectral directions contract. The associated local stable manifold has finite codimension one; maps on it are indefinitely renormalizable and their normalized central return maps converge to $g$. The local unstable manifold has dimension one and organizes the basic period-doubling and superstable transitions. These assertions concern the normalized quadratic universality class, rather than arbitrary functions or every parameter family.

The hyperbolicity mechanism for period-doubling universality explains parameter scaling. Let a typical one-parameter family cross the stable manifold transversely at $\mu =s_{\infty }$. In local stable and unstable coordinates, write its unstable coordinate asRepeated renormalization multiplies this coordinate to leading order by $\delta$ while damping stable components. A superstable orbit of period $2^n$ becomes a superstable two-cycle after $n-1$ renormalizations. Hence its parameter is characterized by arrival at the same basic superstable surface, at a fixed nonzero unstable coordinate after the initial transient. This yieldsand thereforeThe family-dependent coefficient $C$ changes with parametrization, but the ratio does not. The same argument applies to the primary bifurcation surfaces and their successive preimages. Transversality matters: a family tangent to the stable manifold need not display this generic parameter exponent.

Spatial universality follows from the fixed point's scale $a_g=g(1)<0$. As renormalized maps approach $g$, their successive central restrictive intervals shrink by factors approaching $|a_g|$, with orientation alternating. Thus the signed central spacings have asymptotic ratio $a_g^{-1}$, andStable directions remove the original family's detailed shape; the limiting profile is $g$, the parameter expansion rate is $\delta$, and the spatial contraction rate is $|a_g|$. Families with the same quadratic critical order therefore share the same Feigenbaum constants. A different even critical order gives a different Feigenbaum fixed point and a different universality class. This explanation is a consequence of the fixed-point, hyperbolicity and transverse-transition properties; it is not a proof that an arbitrary one-parameter family satisfies those properties.

The change of time is $d\tau /dt=u$, which is positive in the half-plane under consideration. Dividing the two approximate equations by $u$ givesSet $w=v-\mu$. The inhomogeneous system becomes $(u,w{)}^{\prime }=(u-w,\lambda u)$, with eigenvalues $(1\pm \sqrt{1-4\lambda })/2$. For $\lambda >1/4$, write $\omega =\sqrt{4\lambda -1}/2$. Integrating with the prescribed initial data givesandThese expressions also verify both initial values and their initial derivatives. For fixed $\mu >0$ and small $h,v_0$ on the outgoing branch, the sine coefficient in $u$ is positive and dominates the initial cosine term. The first return to small $u$ is near the first positive zero, $\omega \tau \simeq \pi$. HenceThe factor is a half-turn of the growing spiral, not its full rotation period. After this excursion the approximate trajectory approaches the small-$u$, large-$v$ region. In the full system the term $-\lambda v$, omitted during the large excursion, allows the motion to descend along that region and return.

As $\lambda \downarrow 1/4$, the rescaled excursion time diverges and the factor $e^{\tau /2}$ produces arbitrarily large excursions. At the critical value the matrix is a nontrivial Jordan block, and the limiting formula isFor the small outgoing initial data this is positive and unbounded: the unstable manifold no longer turns back. Below $1/4$ the two eigenvalues are real and positive, so there is again no oscillatory return. This quadratic escape desingularization explains loss of the periodic orbit at infinity, together with the unstable-manifold escape, as $\lambda$ decreases through $1/4$. It is not a finite homoclinic bifurcation to the origin. The divergent $\tau$ travel time must not be identified with the physical period without integrating $dt=d\tau /u$; the time change becomes singular near $u=0$.

A global bifurcation joins two mechanisms: a long passage near a hyperbolic equilibrium point and a finite excursion returning the unstable manifold to its neighborhood. Eigenvalue ratios control the local passage, while the geometry of the global return decides reinjection, folding and trapping. In the real symmetric case there are two returning branches and a power-law cusp. In the saddle-focus case rotation adds infinitely many oscillations to that cusp. The following maps display both mechanisms explicitly.

For the Lorenz configuration, enter on $z=h$ at $(x,y,h)$ and leave on $x=sh$, where $s=\mathrm{sgn}(x)$. Solving the linear flow gives a flight timePut $\gamma =-\hat{\lambda }/{\lambda }_{+}>\delta =-{\lambda }_{-}/{\lambda }_{+}$. The outgoing coordinates areThe smooth global map from each outgoing section to the incoming section has, to leading order,up to the choice of which unstable branch carries the sign of $\nu$. Reflection symmetry makes the offsets opposite and the leading coefficients equal. Since the stable transverse coordinate is multiplied by the higher power $|x{|}^{\gamma }$, it contributes only subleading terms to the leading scalar map. Absorbing powers of $h$ into a positive constant gives the Lorenz power return mapThe line $x=0$ lies on the stable manifold and has no finite return; the discontinuity is therefore a genuine missing return, not an arbitrary value to assign at zero.

If $\delta >1$, the derivative $A\delta |x{|}^{\delta -1}$ tends to zero. For $\mu <0$, there are two attracting fixed points, $x=\pm s$ with $s=-\mu +A{s}^{\delta }\simeq -\mu$. They describe two symmetry-related one-lobe attracting flow periodic orbits. For $\mu >0$, the orbit alternates between $s$ and $-s$, where $s=\mu -A{s}^{\delta }\simeq \mu$; its multiplier is $(A\delta s^{\delta -1}{)}^2<1$. This represents one attracting flow periodic orbit traversing both lobes. At $\mu =0$ the two homoclinic loops meet: this is the contracting gluing scenario, with diverging flow period at the connection.

If $\delta <1$, the local derivative is unbounded, and the small periodic orbits are repelling in the expanding return direction. For $\mu >0$, the small positive fixed point has $s\simeq (\mu /A{)}^{1/\delta }$, with a multiplier tending to infinity; there is a symmetry-related negative fixed point as well. More complex motion follows from the two branches, not merely from a large derivative at one fixed point. Choose $I=[-c\mu ,c\mu ]$ with $0<c<1$. The inverse branches areFor sufficiently small positive $\mu$, both lie inside $I$, their images are disjoint, and their derivatives are $O({\mu }^{1/\delta -1})\to 0$. Nested inverse images for every infinite sign sequence therefore define a Cantor invariant set. Iterating the forward map shifts the sequence. This constructs periodic itineraries of every symbolic period and nonperiodic itineraries; transverse contraction in the full return turns the expanding two-strip mechanism into a horseshoe.

This is the local chaotic Lorenz mechanism. A robust attracting Lorenz set additionally needs a global trapping region and the appropriate return properties. The scalar asymptote and $\delta <1$ alone do not prove that every nearby orbit is attracted to one chaotic set; periodic windows, escaping orbits and other global returns must be distinguished. Symmetry is crucial to obtaining the two homoclinic branches together in a single parameter variation.

A Morse function is a smooth real-valued function all of whose critical points have nonsingular Hessian matrices. We use the descending convention for a gradient-like vector field $V$: away from critical points, $df(V)<0$, and near a critical point of Morse index $\lambda$, in Morse lemma coordinates $f=f(p)-|u{|}^2+|v{|}^2$, the field can be chosen as $V=(2u,-2v)$.

For its flow ${\varphi }_t$, define the unstable manifold $W^u(p)=\{x:{\lim }_{t\to -\infty }{\varphi }_t(x)=p\}$ and the stable manifold $W^s(q)=\{x:{\lim }_{t\to +\infty }{\varphi }_t(x)=q\}$. Their dimensions are $\mathrm{ind}p$ and $m-\mathrm{ind}q$. The Morse-Smale transversality condition isFor a closed manifold, this lets us form the Morse-Smale complex over $\mathbb{Z}$. Orient each unstable manifold and takeThe transverse intersection for an index difference of one is one-dimensional; quotienting by flow time gives a finite set of trajectories, and $n(p,q)$ counts them with orientation signs. Compactifying the one-dimensional spaces of trajectories for an index difference of two adds broken trajectories. Their oriented boundary count is zero, giving ${\partial }^2=0$. The descending cells are the unstable manifolds, and these counts are the incidence numbers of the corresponding cellular chain complex, so the homology of this complex is $H_{\ast }(M;\mathbb{Z})$. Ambient orientability is unnecessary: orientations of the unstable cells provide the requisite coorientations and trajectory signs.

For the Klein bottle, use the quotient of ${\mathbb{R}}^2$ by $(s,t)\sim (s+2\pi ,-t)$ and $(s,t)\sim (s,t+2\pi )$. A model of the standard four-critical-point height function isIt is well defined under both identifications. Its critical points, their Morse indices, and critical values areThe Hessian matrix is diagonal with entries $-\cos s$ and $-\epsilon \cos t$, so each point is nondegenerate. The flat quotient metric gives descending equations $\dot{s}=\sin s$, $\dot{t}=\epsilon \sin t$. Their saddle separatrices are the coordinate lines; the two saddle-to-saddle intersections are empty, and all remaining stable and unstable intersections are transverse. This is a Morse-Smale gradient flow. If one insists on the exact local-model version of a gradient-like vector field above, change the two coordinate speeds near the critical points to the standard quadratic-model speeds. This preserves the separatrices, signs and transversality.

Both branches of either saddle flow to the minimum with opposite endpoint signs, so $\partial a=\partial b=0$. From the maximum to $a$, the two trajectories lie on $s=0$ and approach $t=\pi$ from opposite directions. The identification in $t$ is a translation, so their incidence signs are opposite and cancel. The two trajectories to $b$ lie on $t=0$ and approach $s=\pi$ from opposite directions. The identification in $s$ reverses the transverse $t$-coordinate, reversing one additional orientation sign. These two contributions consequently have the same sign. Choosing the orientation of $b$ gives $\partial q=2b$.

One can check the coefficient directly from the cell attachment of the top cell: its boundary word is $ab{a}^{-1}b$, with exponent sums zero on the base loop $a$ and two on the fiber loop $b$. Thus the integral Morse complex of the Klein bottle isandCounting trajectories without signs, or working only modulo two, would miss the order-two class.

Take small incoming sections $y=\pm h$ and outgoing sections $x=\pm h$ after straightening the stable manifold and unstable manifold of the origin. Its eigenvalues are $1,-\delta$, so the leading local passage satisfies $x(t)=x_0{e}^t$, $y(t)=y_0{e}^{-\delta t}$. The exit time and transverse exit coordinate areThus the local passage contracts strongly as $x_0\to 0$, while its duration diverges. Smooth global reinjection from the right outgoing section to the upper incoming section has leading transverse coordinate $-\mu +K_R{y}_{\mathrm{out}}$; the corresponding left-to-lower return has $\nu +K_L{y}_{\mathrm{out}}$. Choose the parameters as these signed splitting coordinates and absorb $h^{1-\delta }$ into $A,B$.

After scaling the incoming-section labels to $y_n=\pm 1$, the leading Poincaré return map isThe second coordinate records which outgoing branch was used, not an independent continuous amplitude. A smooth planar flow preserves the orientation of a transverse-section return. For the right-to-upper pair of sections, the ordered pair consisting of the flow direction and the positive section tangent has the same orientation at both ends; the same is true for the left-to-lower pair. Hence $A,B>0$ for these sign conventions and a nondegenerate global return.

For general nonlinear $P,Q$, this is an asymptotic return model, with smooth parameter corrections and higher-order local-passage corrections. The hypotheses do not imply an exact power-law map in the original coordinates. Interpreting the printed parameters themselves as linear splitting coordinates also requires a nondegenerate two-parameter unfolding: the two splitting derivatives must be nonzero. For example, composing a generic right-loop splitting with $\mu \mapsto {\mu }^2$ preserves the stated loop at $\mu =0$ but does not give the printed signed $-\mu$ offset. The subsequent region diagram uses the standard generic splitting-coordinate interpretation. Nor is $x=0$ a normal return point: it lies on the stable manifold, whose trajectory takes infinite time to reach the saddle equilibrium. This explains how the map describes nearby recurrent trajectories without extending a finite-time return across the singular separatrix.

For $r>0$, polar coordinates give $\dot{r}=r({\mu }_1-2x)$ and $\dot{\theta }=1$; the remaining amplitude equation is $\dot{x}={\mu }_2+x^2+r^2$. Interpret $r=0$ in the original complex coordinate, where the phase is undefined.

If ${\mu }_2<0$, put $a=\sqrt{-{\mu }_2}$. The full equilibrium points are $E_{\pm }:(z,x)=(0,\pm a)$. Their eigenvalues areThus $E_{-}$ is stable for ${\mu }_1<-2a$ and has a two-dimensional unstable manifold for ${\mu }_1>-2a$. The point $E_{+}$ has a one-dimensional unstable manifold for ${\mu }_1<2a$ and is repelling for ${\mu }_1>2a$. On ${\mu }_2=0$ they merge in a saddle-node bifurcation; the simultaneous zero and imaginary eigenvalues at $({\mu }_1,{\mu }_2)=(0,0)$ give a fold-Hopf bifurcation.

A steady positive amplitude represents a full periodic orbit, not a full equilibrium. Solving the amplitude equations givesIt exists in the wedge ${\mu }_2<-{\mu }_1^2/4$. The transverse amplitude Jacobian matrix has characteristic polynomial $s^2-{\mu }_{1}s+4r_{\ast }^2$. Its two transverse Floquet exponents have negative real part for ${\mu }_1<0$ and positive real part for ${\mu }_1>0$; the phase direction supplies the neutral Floquet multiplier one.

The wedge boundaries are Hopf bifurcations at ${\mu }_1=-2a$ from $E_{-}$ and ${\mu }_1=2a$ from $E_{+}$. The left branch is supercritical, producing a stable periodic orbit as $E_{-}$ loses stability; the right branch is subcritical in the oscillatory subspace, with an unstable periodic orbit on the side where the transverse complex pair of $E_{+}$ is stable. The real unstable direction of $E_{+}$ persists throughout. There are no full equilibrium points when ${\mu }_2>0$, and $\dot{x}>0$ then excludes a local periodic orbit. The right panel of the bifurcation diagram in part 1(b) shows the primary curves and the perturbed global curve calculated below.

Choose real constants $a_0<a_1<\cdots <a_m$ and define the smooth homogeneous quotientIt is unchanged by nonzero rescaling, so it is well defined on Real projective space. On the unit sphere its differential vanishes on the tangent space precisely when $Ax$ is proportional to $x$, where $A=\mathrm{diag}(a_0,\dots ,a_m)$. The entries are distinct, so the only critical projective lines are $p_i=[e_i]$. In the affine chart $x_i=1$, write $u_j=x_j/x_i$ for $j\ne i$. ThenIts Hessian at zero is diagonal with entries $2(a_j-a_i)$. Exactly $i$ are negative, and none vanish. Thus this is a Morse function with one critical point $p_i$ of index $i$ for every $0\le i\le m$.

Choose the quotient of the round metric on the unit sphere. The downward gradient flow has sphere equation ${\dot{x}}_j=-2(a_j-f(x))x_j$ and projective solutionNormalization to the unit sphere supplies the common positive factor. As $t\to +\infty$, the smallest index with nonzero coordinate dominates; as $t\to -\infty$, the largest dominates. ConsequentlyThese are its unstable manifolds and stable manifolds. If $j\ge k$, the tangent directions in $W^u(p_j)$ allow coordinates $0,\dots ,j$, and those in $W^s(p_k)$ allow $k,\dots ,m$, modulo the common projective line. Together they span the full tangent space at every intersection. If $j<k$ the intersection is empty. Thus the flow is a Morse-Smale gradient flow, with no generic perturbation needed.

Orient each unstable manifold using its affine coordinates. The Morse-Smale complex has $C_i=\mathbb{Z}{p}_i$; its differential counts oriented unparametrized connecting trajectories of index difference one. A trajectory from $p_i$ to $p_{i-1}$ has only coordinates $x_{i-1},x_i$ nonzero, and there are exactly two, distinguished by the sign of $x_{i-1}/x_i$. Each sign branch is a single orbit modulo time translation.

To calculate their signs rather than merely count two, use the characteristic closed disk of the unstable cell: the closed upper hemisphere of $S^i$ with $x_i\ge 0$, followed by projection onto ${\mathbb{RP}}^i$. Its interior is $W^u(p_i)$, and its boundary map is the antipodal double cover $S^{i-1}\to {\mathbb{RP}}^{i-1}$. After collapsing ${\mathbb{RP}}^{i-2}$, its degree on the $(i-1)$-cell is the sum of the degrees of the two sheets. Comparing the sheets by the antipodal map on $S^{i-1}$ gives relative orientation $(-1{)}^i$: it is induced by $-I$ on ${\mathbb{R}}^i$, whose determinant is $(-1{)}^i$. Thus the two trajectory signs add to $1+(-1{)}^i$, up to the choice of generator orientations. The projective quadratic Morse complex is thereforeIn particular ${\partial }_1=0$, as required for the connected one-dimensional skeleton. The signed trajectory count agrees with this attaching-map incidence calculation, so it computes integral rather than only mod-two homology.

For $m\ge 1$, taking kernels modulo images givesFor example, an odd interior-degree generator has zero outgoing differential but twice itself is the image from the next degree; an even positive-degree differential is injective. At the top there is no incoming differential, leaving a copy of $\mathbb{Z}$ exactly when $m$ is odd. For $m=0$ the space is one point and only $H_0=\mathbb{Z}$ remains. This recovers the cellular homology of real projective space by explicitly constructing its Morse-Smale complex and signs.

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.

For a nonconstant invasion travelling wave, the relevant end states are $w(-\infty )=k$ and $w(+\infty )=0$. These fronts do not exist at every nonnegative speed: the phase plane calculation gives the threshold $2\sqrt{k}$. The constant solutions $w\equiv 0$ and $w\equiv k$ are, of course, nonnegative travelling waves for every speed, including zero.

Put $z=x-ct$ and $p=w^{\prime }$. The travelling-wave reduction of a reaction-diffusion system givesThe equilibrium points in the phase plane are $(k,0)$ and $(0,0)$. At $(k,0)$ the characteristic equation is $r^2+cr-k=0$, so there is one positive and one negative eigenvalue. Thus it is a saddle equilibrium. Its descending unstable manifold has $w<k$, $p<0$ and locallyAt $(0,0)$ the eigenvalues areFor $c>2\sqrt{k}$ it is a stable node; for $0<c<2\sqrt{k}$ it is a stable spiral. The repeated-root case $c=2\sqrt{k}$ must be included in the existence argument rather than inferred from distinct roots.

Suppose $c\ge 2\sqrt{k}$. Choose $a>0$ with $a(c-a)\ge k$; for example $a=(c-\sqrt{c^2-4k})/2$. Consider the closed triangular regionOn the horizontal side $p=0$, the vector field has $p^{\prime }=-w(k-w)\le 0$, pointing into $\mathcal{T}$. On the vertical side $w=k$, it has $w^{\prime }=p\le 0$, again pointing inward. For the sloping side use $h=p+aw$, so that on $h=0$,Thus trajectories cannot cross into $h<0$. The remaining boundary point $(0,0)$ is an equilibrium point and cannot be crossed in finite time by a nonconstant trajectory, by uniqueness theorem for ordinary differential equations. This proves forward invariance of $\mathcal{T}$. The argument works unchanged at $c=2\sqrt{k}$ with $a=\sqrt{k}$.

The descending unstable manifold of $(k,0)$ initially enters $\mathcal{T}$. Its forward trajectory remains in this compact region, so it exists for all positive $z$. Within the interior, $p<0$: a first contact with $p=0$ at $0<w<k$ would have $p^{\prime }<0$, contrary to arrival from below. Hence $w$ decreases strictly and has a limit $\ell \in [0,k)$. To identify the limit, its omega-limit set lies on the line $w=\ell$ and is invariant. Any trajectory remaining on this line has $p=w^{\prime }=0$; its second equation then forces $\ell (k-\ell )=0$. Since $\ell <k$, it follows that $\ell =0$, and the only possible omega-limit set is $(0,0)$. Backward in $z$, the unstable manifold tends to $(k,0)$ by its construction. We have therefore obtained a positive decreasing heteroclinic orbit and hence a nonnegative travelling wave for every $c\ge 2\sqrt{k}$.

For completeness, the threshold is necessary for such a front. When $0<c<2\sqrt{k}$, put $\nu =\sqrt{k-c^2/4}$, $X=w$, and $Y=(p+cw/2)/\nu$. Near the zero equilibrium point the system becomesFor a nonzero trajectory its polar angle satisfiesA trajectory approaching zero therefore makes infinitely many rotations and must cross $X=w=0$. It cannot be a nonnegative front. At $c=0$, the first integralis conserved. It has values $k^3/6$ and $0$ at the two desired end states, so no heteroclinic orbit connects them. ConsequentlyThis is the phase-plane construction of a Fisher travelling front, including the critical speed. The statement for all $c\ge 0$ must distinguish these nonconstant fronts from the two constant travelling waves.

The nullclines are $y=\mu x$ and $y=x^2/[\nu (1+x^2)]$. At a nonzero fixed point, cancellation of $x$ yields $\mu \nu =x/(1+x^2)$. The latter function attains its maximum $1/2$ at $x=1$. HenceFor $0<\mu <{\mu }_c$, $x_{-}<1<x_{+}$ and $x_{-}{x}_{+}=1$.

The Jacobian matrix isAt $O$ the distinct eigenvalues are $-\mu ,-\nu$, so the origin is a stable node. At a nonzero fixed point its trace is $-(\mu +\nu )<0$ andThe discriminant is $(\mu -\nu {)}^2+8x/(1+x^2{)}^2>0$. Thus $P_{-}$ is a saddle and $P_{+}$ is a stable node, never a focus.

Above the straight nullcline $\dot{x}>0$; below it $\dot{x}<0$. Below the curved nullcline $\dot{y}>0$; above it $\dot{y}<0$. At the saddle, an eigenvector can be chosen as $(1,\mu +\lambda )$: the stable eigenvalue gives negative slope, and the unstable eigenvalue gives positive slope. The stable manifold is the basin boundary. The two branches of the unstable manifold lead respectively to the origin and the positive stable node, as illustrated for representative parameters.

The nonnegative quadrant is forward invariant. Since $\dot{y}\le 1-\nu y$, $y$ stays bounded, and then $\dot{x}=-\mu x+y$ also bounds $x$. Moreover the divergence is the negative constant $-(\mu +\nu )$, excluding periodic or homoclinic loops by the Bendixson-Dulac criterion. The Poincaré-Bendixson theorem then leaves convergence to a fixed point. Initial data on the saddle's stable manifold approach the saddle; data in the two open basins approach extinction at $O$ or persistence at $P_{+}$, respectively.

At a fixed point, the first equation gives $y=-x$. Substitution into the second gives $2x(x^2-1)=0$, so the equilibrium points are $(0,0),(1,-1),(-1,1)$. The Jacobian matrix isAt the origin its eigenvalues are $1\pm i$, so the origin is an unstable spiral. The local rotation is clockwise: on the positive $x$-axis the velocity points downward. At either nonzero equilibrium point,Thus both nonzero equilibria are saddles, with eigenvalues $-2,2$. The stable manifold is tangent to $(1,-3)$, and the unstable manifold to $(1,1)$, at each saddle equilibrium. Since all these equilibrium points are hyperbolic, linear stability of a planar equilibrium gives their actual local types.

For the phase portrait, the nullclines are $y=-x$ and $y=x/(1-2x^2)$, the latter defined away from $x=\pm 1/\sqrt{2}$. The horizontal direction has the sign of $x+y$; the vertical direction has the sign of $(1-2x^2)y-x$. The system is invariant under $(x,y)\mapsto (-x,-y)$, so trajectories occur in centrally symmetric pairs. The arrows below show these directions; the highlighted stable manifolds enter the saddle equilibria, whereas the highlighted unstable manifolds leave them.