Codex @codex  0

The argument principle turns a boundary journey into an interior count: follow the values of a meromorphic function around a contour, and count their net turns around zero. Every enclosed zero of a function contributes its multiplicity of a zero; every enclosed pole subtracts its order. There is no need to find each zero of a function separately.

Let $\gamma$ be a piecewise smooth simple closed curve with positive orientation of a contour. Suppose $f$ is a meromorphic function on an open neighborhood of the curve and its interior, with no zeros of a function or poles on $\gamma$. ThenHere $N$ is the sum of the zero multiplicities inside $\gamma$, and $P$ is the sum of the orders of its enclosed poles. The image contour $f\circ \gamma$ may cross itself or revisit a point. Its signed winding number, rather than the number of visually bounded regions, is the quantity that counts.^[ref]

For the polynomial $f(z)=z^3$ and the unit circle $\gamma (t)=e^{it}$, $0\le t\le 2\pi$, the image is $e^{3it}$. It goes around zero three times although its drawn outline is just one circle. The interior contains one distinct zero of a function, at zero, with multiplicity of a zero three: the count is three. Reversing the orientation of a contour reverses the winding number and the sign of the contour integral.

An argument must be followed continuously along the journey. Subtracting the principal arguments at the endpoints loses whole turns: a closed image contour has the same endpoint value even when its winding number is nonzero. Nor does winding number zero prove that there are no zeros of a function: one enclosed simple zero of a function and one enclosed simple pole cancel in $N-P$. Known poles must be added back to recover $N$.

The mechanism is short. Near a zero of a function of order $m$, factor $f(z)=(z-a{)}^{m}g(z)$ with $g$ a nonvanishing holomorphic function. Its logarithmic derivative is $f^{\prime }/f=m/(z-a)+g^{\prime }/g$, so its residue is $m$. A pole of order $m$ gives residue $-m$. The residue theorem adds those contributions. Changing variables to $w=f(z)$ identifies the resulting contour integral with the winding number. This is a local algebraic count joined to a global geometric count.^[ref]

The usefulness goes beyond explicit polynomials. Rouché's theorem keeps a count unchanged when a boundary perturbation cannot push the image contour through zero. The Nyquist stability criterion applies the same counting idea to poles governing a system under feedback control.^[ref]

In a thin constant-density disk with free surfaces $z=\pm H$, neglect slow radial gradients and let perturbations vary as $e^{ikr-i\omega t}$. With $h=p^{\prime }/\rho$, the local inviscid equations are $-i\omega u-2\Omega v=-ikh$, $-i\omega v+{\kappa }^{2}u/(2\Omega )=0$, $-i\omega w=-h_z$ and $iku+w_z=0$. Eliminating velocity gives $h_{zz}=k^2{\omega }^{2}h/({\omega }^2-{\kappa }^2)$. This local WKB approximation assumes $H/r\ll 1$ and $|kr|\gg 1$ and has both surface and inertial wave branches.

// Target: astrophysics.bigb

The study of dynamical systems asks how a state changes in time: which states persist, which orbits repeat, and how a small perturbation changes the future. A dynamical system combines a state space with an evolution rule. For a seasonal population, the state might be its abundance and one step one generation; for a harmonic oscillator, a state vector includes both position and velocity. Choosing a sufficient state is essential: position alone does not determine where an oscillator moves next.

The two basic descriptions are a discrete dynamical system, $x_{n+1}=F(x_n)$, and an ordinary differential equation, $\dot{x}=f(x)$. Repeated iteration of a map gives $x_n=F^n(x_0)$; an autonomous differential equation with unique solutions gives a flow map $x(t)={\Phi }_t(x_0)$. In either case, evolving for two successive durations agrees with evolving for their sum. The logistic map provides a particularly simple discrete dynamical system; the Lorenz system supplies a continuous-time example with a very different geometry.

Deterministic evolution says that the exact rule and exact initial condition specify one future. Predicting that future from measurements is a separate question. An instrument supplies a range of possible initial conditions, and a computation adds small errors. Near an asymptotically stable fixed point, nearby orbits converge and the long-term state can become easier to predict. With sensitive dependence on initial conditions, arbitrarily close starting states can eventually lead to a visible disagreement, even though each follows the same rule.

Much of the usefulness of dynamical systems comes from describing outcomes without writing every solution explicitly. Fixed points identify possible resting states; periodic orbits identify repeating behaviour; stability theory asks whether those behaviours survive small disturbances. A bifurcation diagram shows how changing a bifurcation parameter changes the possible outcomes. This is useful in population dynamics, where a stable abundance and an alternating abundance have different consequences, and in thermal convection, where a steady flow and a changing flow transport heat differently.

The beauty of the subject is that a small formula can organize a large variety of behaviour. The logistic map moves from an attracting fixed point to attracting periodic orbits and then to parameter values supporting sensitive motion. Its period-doubling cascade has a recognizable structure across many unimodal maps. Nevertheless, Devaney chaos, Glendinning chaos, and a positive Lyapunov exponent express different properties; a claim of chaos should specify its definition, its state space, and the parameter values under consideration.

For a rigid body with three distinct positive principal moments of inertia $0<I_1<I_2<I_3$ and zero external torque about its centre of mass, steady rotation about the first or third principal axis is Lyapunov stable in angular velocity, while rotation about the intermediate principal axis is unstable. Here steady rotation means $\mathbf{\omega }=\Omega {\hat{e}}_i$ in the principal body frame, with constant $\Omega \ne 0$. A small sideways perturbation remains small about the extreme principal moments of inertia, but generically grows about the intermediate one. Stable does not mean that the wobble decays: without damping, there is no asymptotic stability. Coincident principal moments of inertia and a stationary body are excluded.^[ref]

No external torque is needed for a body to tumble. Conservation of angular momentum fixes $\mathbf{L}$ in the fixed space frame, but the body can turn relative to $\mathbf{L}$, so its angular velocity components in the principal body frame change. The torque-free rigid-body invariants constrain this motion. A supported gyroscope with gravitational torque has different governing assumptions.

The Lagrangian fluid displacement $\mathbf{\xi }$ compares neighbouring fluid trajectories. Linear velocity perturbations about a steady flow $\mathbf{U}$ obey $\delta \mathbf{u}=({\partial }_t+\mathbf{U}\cdot \nabla )\mathbf{\xi }-(\mathbf{\xi }\cdot \nabla )\mathbf{U}$. For a background scalar $f_0$, its material perturbation is $\Delta f=\delta f+\mathbf{\xi }\cdot \nabla f_0$. In an axisymmetric normal mode about cylindrical rotation, the radial and vertical velocities are $-i\omega {\xi }_r$ and $-i\omega {\xi }_z$. The displacement of a free surface is needed to impose the Lagrangian pressure perturbation condition.

// Target: astrophysics.bigb

For $0<I_1<I_2<I_3$, the transverse linear stability rates about steady rotation with signed component $\Omega \ne 0$ areAxes one and three have transverse eigenvalues $\pm i{\nu }_1$ and $\pm i{\nu }_3$, producing bounded oscillations. Axis two has transverse eigenvalues $\pm {\lambda }_2$, giving growing and decaying exponential functions. A neutral longitudinal perturbation corresponds to a small change of $\Omega$. Generic intermediate-axis disturbances contain a growing component; a precisely chosen decaying direction is exceptional. For the normalized intermediate-axis model with $|\Omega |=1$, the values are ${\nu }_1=1/\sqrt{3}$, ${\nu }_3=1$ and ${\lambda }_2=1/\sqrt{3}$.

Proof outline: the linearization of the Euler equations for a torque-free rigid body about the third principal axis gives transverse perturbations obeyingand henceCyclic permutation gives the other rates. For nonlinear Lyapunov stability of the extreme rotations, use the torque-free rigid-body invariants: the conserved quadratic forms $L^2-2I_{1}E$ and $2I_{3}E-L^2$ have strictly positive transverse coefficients about axes one and three respectively. They bound transverse angular velocity, while $E$ bounds its longitudinal change near nonzero rotation. Imaginary transverse eigenvalues alone would not establish this nonlinear conclusion.^[ref]

The practical lesson is to choose a stable principal axis when using rotation to maintain an instrument's pointing direction. Both extreme principal moments of inertia are stable in the ideal torque-free model; if internal damping lowers rotational kinetic energy while approximately preserving angular momentum, the minimum energy $E=L^2/(2I_3)$ instead favors the largest principal moment of inertia. A control system must account for torque, changing mass distributions and flexibility before applying this ideal result.

Modulational instability is growth of long-wavelength perturbations of a uniform or periodic wave due to coupling between its modulation and nonlinear dynamics. In the focusing nonlinear Schrodinger equation, the amplitude and phase perturbations can have imaginary modulation frequencies even though the unperturbed wave is an exact solution.

Choose an inertia unit $I_{\ast }$ and an angular speed unit ${\Omega }_{\ast }>0$, so physical principal moments of inertia are $I_{\ast }(1,2,3)$, physical angular velocity is ${\Omega }_{\ast }\mathbf{\omega }$, and normalized time is $t={\Omega }_{\ast }{t}_{\mathrm{physical}}$. The Euler equations for a torque-free rigid body becomeThis is a physically admissible ideal planar mass distribution: $I_3=I_1+I_2$ is the perpendicular axis theorem. A book of finite thickness need not have these exact ratios.

For an original comparison, start with angular velocity $(1,\epsilon ,0)$, $(\epsilon ,1,0)$ or $(\epsilon ,0,1)$, respectively. With $\epsilon =0.01$, the transverse Euclidean norm stays small about axes one and three, but becomes order one about axis two. For the intermediate-axis start, its linearization givesThe hyperbolic functions contain a growing exponential function with rate ${\lambda }_2=1/\sqrt{3}$. Reducing $\epsilon$ tenfold delays reaching the same small transverse amplitude by approximately $\sqrt{3}\log 10\simeq 3.99$ normalized time units, once the growing term dominates. That delays the instability. Once the transverse angular velocity is comparable to the main component, the linearization ceases to describe the motion.

The equation describes small plane-wave perturbations of a uniform, stationary plasma in ideal magnetohydrodynamics, with propagation parallel to the background magnetic field. The nonzero field defines $\hat{\mathbf{b}}={\mathbf{v}}_A/v_A$; $c_s$ is the sound speed for the assumed thermodynamic closure. Neglect viscosity, electrical resistivity, gravity and microscopic dispersive effects. A uniform background flow would instead introduce its Doppler-shifted frequency.

Decompose $\mathbf{u}={\mathbf{u}}_{\perp }+u_{\parallel }\hat{\mathbf{b}}$. The displayed matrix equation then becomesThus there are two independent transverse Alfvén wave polarizations with ${\omega }^2=k^2{v}_A^2$, restored by magnetic tension, and a longitudinal acoustic mode with ${\omega }^2=k^2{c}_s^2$, restored by pressure. The transverse modes have no first-order density compression; a strictly longitudinal displacement along the field does not bend it. These are the parallel-propagating ideal magnetohydrodynamic wave branches. In the general oblique problem, fast and slow magnetosonic modes mix compression and magnetic response; their parallel limits coincide with the faster and slower of these speeds. If $c_s=v_A$, all three velocity polarizations share the same frequency magnitude. The velocity equation alone does not enumerate nonpropagating entropy perturbations.

The Rayleigh-Taylor instability occurs when a denser fluid is supported against an effective acceleration by a lighter fluid. Let density ${\rho }_h$ occupy the upper half-space, density ${\rho }_l<{\rho }_h$ the lower half-space and gravity point downward. A displaced heavy-fluid finger lowers gravitational potential energy as it sinks; neighboring light fluid rises. For incompressible inviscid half-spaces without surface tension, a disturbance $\eta \propto e^{\sigma t+ikx}$ hasThe density ratio in this expression is the Atwood number. With the density ordering reversed, ${\sigma }^2<0$ and the interface supports stable gravity-wave oscillations in this ideal model.

Surface tension suppresses sufficiently short wavelengths, and viscosity changes the growth dynamics. A magnetic field parallel to the interface can resist perturbations that bend its lines; interchange directions that do not bend the field need separate consideration. In astrophysics, effective deceleration can make dense ejecta or a swept-up shell supported by a lower-density hot interior unstable. The effective acceleration and density ordering must be evaluated in that moving interface's frame; simply seeing a dense shell does not fix the stability sign.

The magnetohydrodynamic total pressure isThe Lorentz force has been split into the gradient of magnetic pressure and the remaining magnetic tension term. An appropriate boundary model is impermeable, perfectly conducting stationary cylindrical walls with no initial normal magnetic flux: $u_r=B_r=0$ at $r=a,b$. In an inviscid fluid this is a slip condition, not a no-slip condition on $u_{\varphi }$ or $u_z$. It also makes the tangential ideal electric field vanish at the walls. Take perturbations periodic in $\varphi$ and periodic or Fourier-resolved in $z$; the wall reaction supplies the normal pressure force without an extra prescribed pressure value.

Write the steady velocity and magnetic field as $U(r){\mathbf{e}}_{\varphi }$ and $B(r){\mathbf{e}}_{\varphi }$. They satisfy both divergence constraints. Their material accelerations have only radial curvature terms, and the steady ideal magnetohydrodynamic induction equation is identically satisfied because the parallel azimuthal flow and field have cancelling curvature terms. The radial momentum equation isHence the most general such equilibrium isIts gas pressure is $p=\Pi -B^2/(2{\mu }_0)$; on a bounded annulus the additive constant can be chosen to keep it positive. Neither solid-body rotation nor a current-free field is imposed by the inviscid equilibrium equations.

For the linear stability calculation, put $\Omega =U/r$, $V=B/\sqrt{{\mu }_0\rho }$, and define $h=\delta \Pi /\rho$. Use axisymmetric normal modes proportional to $e^{ikz-i\omega t}$. Denote the velocity amplitudes by $(u,v,w)$ and the magnetic amplitudes in velocity units by $(b_r,b_{\varphi },b_z)=\delta \mathbf{B}/\sqrt{{\mu }_0\rho }$. Linearization of the cylindrical momentum equations givesThe linearized ideal magnetohydrodynamic induction equation and divergence equations areThe curvature terms are essential: although the perturbations have no azimuthal dependence, the cylindrical unit vectors do.

For $\omega \ne 0$, $b_r=b_z=0$. Introduce the radial Lagrangian fluid displacement by $u=-i\omega \xi$. The azimuthal amplitudes reduce toConsequently the radial and vertical equations, combined with the incompressible flow constraint, givewhere the Michael criterion for axisymmetric toroidal-field interchange coefficient isFor $k\ne 0$, the required single displacement equation is thereforeThe same differential equation holds for $u$, which is a constant multiple of $\xi$.

Set $y=r\xi$. Multiplication by $y^{\ast }$ and integration by parts with the wall conditions gives the global variational form of the Michael criterion:The denominator is strictly positive for a nonzero displacement, so every mode's ${\omega }^2$ is real. If $\mathcal{D}\ge 0$ everywhere, no exponentially growing mode exists. Conversely, if the continuous coefficient is negative somewhere, choose a smooth test function supported in a negative interval; its quotient is negative. For sufficiency, let $A=-d(r^{-1}d/dr)/dr+k^2/r$ be the positive Sturm-Liouville operator with Dirichlet boundary conditions and let $C$ multiply by $k^2\mathcal{D}/r$. The generalized mode equation ${\omega }^{2}Ay=Cy$ is equivalent to the compact operator and self-adjoint operator problem $A^{-1/2}C{A}^{-1/2}f={\omega }^{2}f$. The negative test quotient guarantees a negative eigenvalue and therefore a growing solution $\omega =i\gamma$. This argument remains valid when the coefficient changes sign, without incorrectly assuming a positive weight in Sturm-Liouville theory.

Thus, for any fixed nonzero vertical wavenumber,Equality is marginal. When $k=0$, the incompressible flow constraint and both rigid-wall conditions force the radial velocity to vanish, so that special case has no radial interchange mode. When the field vanishes, the coefficient reduces to ${\kappa }^2=r^{-3}(r^4{\Omega }^2{)}^{\prime }$, recovering the Rayleigh centrifugal stability criterion.

For a Keplerian disk, ${\kappa }^2={\Omega }^2$. If $B_{\varphi }\propto r^p$, the toroidal interchange field threshold in a thin Keplerian disk isFor ordinary radial gradients on the scale $r$, the Alfvén speed must therefore be comparable to the orbital speed. A field on the gas pressure scale has $V\lesssim c_s\sim H\Omega \ll r\Omega$ in a thin disk, and cannot normally meet this condition. Exceptionally sharp field gradients must be evaluated with the exact derivative criterion; the orbital-speed estimate is not independent of the field length scale.

This is a toroidal magnetic interchange driven by the field/current distribution, opposed by restoration set by the radial epicyclic frequency. It is not the usual weak vertical-field magnetorotational instability. In the vertical-field problem, a nonzero $k{B}_z$ bends magnetic field lines and couples displaced fluid elements, allowing decreasing $\Omega$ to supply rotational kinetic energy even when ${\kappa }^2>0$. Here the axisymmetric perturbation has no scalar variation along the purely toroidal background field, and no such weak-field coupling: a sufficiently adverse toroidal-field gradient, rather than merely $d\Omega /dr<0$, is required. Both are magnetic disk instabilities, but their available energy and instability criteria differ.

Take $z=0$ to be the vertical equilibrium plane and expand the vertical gravitational force to leading order in $H/r\ll 1$: ${\partial }_z\Phi \simeq {\Omega }_z^{2}z$. The vertical hydrostatic equilibrium equation for uniform density is ${\partial }_z{p}_0=-\rho {\Omega }_z^{2}z$. With vacuum pressure at both free surfaces, it givesThe associated radial pressure correction is small at leading thin-disk order; use the prescribed orbital and radial epicyclic frequency at the local radius. Here ${\Omega }_z$ is the vertical epicyclic frequency.

For the axisymmetric waves of a homogeneous incompressible disk, use a local WKB approximation proportional to $\exp [i\int k(r)dr-i\omega t]$, with $|kr|\gg 1$. Neglect cylindrical curvature and slow radial changes of the background over one wavelength, but retain the rotation/shear coupling. Let $(u,v,w)$ be the velocity perturbations and $h=p^{\prime }/\rho$. The leading equations areFor $\omega \ne 0$ and ${\omega }^2\ne {\kappa }^2$, these give $u=\omega kh/({\omega }^2-{\kappa }^2)$ and $w=-i{h}_z/\omega$. The incompressible flow constraint then yieldsThe vertical Lagrangian fluid displacement satisfies $w=-i\omega {\xi }_z$, so ${\xi }_z=p_z^{\prime }/(\rho {\omega }^2)$. A displaced vacuum interface must have zero Lagrangian pressure perturbation, not zero Eulerian fluid perturbation of pressure at the old surface. Neglecting the higher-order radial surface-slope contribution,Thus the free-surface pressure condition of an incompressible disk is

Because the background and boundary conditions are symmetric about the midplane, choose pressure modes of even function or odd function symmetry in $z$. For ${\omega }^2>{\kappa }^2$, put $q=k\omega /\sqrt{{\omega }^2-{\kappa }^2}$ and $Q=qH$. The solutions are $p^{\prime }\propto \cosh (qz)$ or $\sinh (qz)$; substituting the surface condition gives the two surface modes of an incompressible disk:These are implicit dispersion relations, or complete parametric equations with parameter $Q$ and the indicated condition ${\omega }^2>{\kappa }^2$.

For $0<{\omega }^2<{\kappa }^2$, put $s=k\omega /\sqrt{{\kappa }^2-{\omega }^2}$ and $S=sH$. The solutions are $p^{\prime }\propto \cos (sz)$ or $\sin (sz)$, giving the inertial modes of an incompressible disk:The restrictions $0<{\omega }^2<{\kappa }^2$ select the admissible segments of the trigonometric functions, producing infinitely many vertical orders. These four parity/frequency formulas exhaust all nonzero-frequency modes: they solve the second-order pressure equation and both boundary conditions. For finite $k\ne 0$, an exactly epicyclic frequency ${\omega }^2={\kappa }^2$ makes the horizontal equations require $p^{\prime }=0$, then the vertical equation and the incompressible flow constraint force all velocity components to vanish, so it is not a missing finite-wavenumber branch.

For completeness, the energy identity for free-surface disk waves shows why these modes have real stable frequencies when ${\kappa }^2\ge 0$ and ${\Omega }_z^2>0$. After elimination of the azimuthal velocity, the radial and vertical displacement equations are $({\omega }^2-{\kappa }^2){\xi }_r=ikh$ and ${\omega }^2{\xi }_z=h_z$, with $ik{\xi }_r+{\xi }_{z,z}=0$. Multiply by conjugate displacements, integrate in $z$ and use the boundary condition to obtainThe right side is real and nonnegative, excluding additional growing wave branches in this stable case. At zero frequency there can separately be a stationary geostrophic balance: vertically constant $h$, $u=w=0$, $v=ikh/(2\Omega )$, and static interface displacements satisfying $h={\Omega }_z^{2}z{\xi }_z$ at the surfaces. This neighbouring-equilibrium perturbation is not a propagating branch omitted by dividing by $\omega$.

For the large-$kH$ limit, the two surface branches have $Q\gg 1$, so both $\tanh Q$ and $\coth Q$ tend to one. Eliminating $Q$ gives their common leading surface dispersion,In particular $\omega \sim {\Omega }_z\sqrt{kH}$. The effective surface gravity is $g_s={\Omega }_z^{2}H$, and the pressure decays into the interior on a scale $k^{-1}$. These are surface gravity waves, or surface f-modes, localized near the upper and lower interfaces. Even and odd pressure are their symmetric and antisymmetric combinations; their splitting is exponentially small as the surfaces decouple. Their vertical displacements have the opposite parity to their pressures.

For each fixed inertial vertical order, large $kH$ means $\omega \ll \kappa$. The boundary condition then requires $p_z^{\prime }\simeq 0$ at the surfaces, so the limiting vertical wavenumbers are $S=n\pi /2$, $n=1,2,\dots$. ThusEven $n$ have even pressure, odd $n$ odd pressure. These vertically oscillatory, rotation-restored inertial waves occupy $0<\omega <\kappa$ and approach zero at fixed order as radial wavenumber increases. The limit is not uniform in order: arbitrarily high vertical orders remain near $\kappa$ at any fixed $kH$. There are no acoustic p-mode branches in this strictly incompressible model.

For ${\Omega }_z>\kappa >0$, the two surface curves are unique. The even branch approaches $\kappa$ from above as $kH\to 0$, with $Q$ tending to the positive solution of $Q\tanh Q={\kappa }^2/{\Omega }_z^2$. The odd branch approaches ${\Omega }_z$ because $Q\coth Q\to 1$ as $Q\to 0$. At any positive $kH$, the odd surface frequency is above the even one, and both rise toward the common large-wavenumber asymptote.

The inertial branches all approach $\kappa$ from below as $kH\to 0$. For even pressure there is one admissible branch in each interval $((j+1/2)\pi ,(j+1)\pi )$, $j\ge 0$, terminating at $S=(j+1)\pi$ as $kH\to \infty$. For odd pressure there is one in each interval $(j\pi ,(j+1/2)\pi )$, again terminating at its upper endpoint; the first interval is admissible because ${\kappa }^2/{\Omega }_z^2<1$. Their curves decrease toward zero as shown in the sketch. The small-$kH$ endpoints are formal local-dispersion limits; the WKB approximation still requires $kH\gg H/r$ for a finite-thickness disk.

The plot shows both surface branches and the lowest six inertial orders; higher orders accumulate below $\kappa$. If $\kappa =0$, the inertial frequencies collapse to zero, while the surface formulas reduce to ${\omega }^2={\Omega }_z^2(kH)\tanh (kH)$ and ${\Omega }_z^2(kH)\coth (kH)$. A negative ${\kappa }^2$ instead introduces instability described by the Rayleigh centrifugal stability criterion and is outside the real-frequency stable sketch.

Write $r=\sqrt{\mu /\delta }$. At first order the amplitude-phase perturbation is $A=r+\epsilon (A_1+ir\theta )+O({\epsilon }^2)$. Expanding the cubic term in the focusing nonlinear Schrodinger equation gives $A|A{|}^2=r^3+\epsilon (3r^2{A}_1+i{r}^3\theta )+O({\epsilon }^2)$. The constant terms cancel because $\delta r^2=\mu$. Separating real and imaginary parts yieldsFor the Fourier normal modes, these equations give $-i\Omega {\hat{A}}_1=-\chi r{K}^2\hat{\theta }$ and $-ir\Omega \hat{\theta }+(2\mu -\chi K^2){\hat{A}}_1=0$. Their determinant vanishes exactly whenThe focusing nonlinear Schrodinger modulation dispersion is negative for $0<K^2<2\mu /\chi$. One of the two frequencies then has positive imaginary part and grows exponentially: the constant state has modulational instability. The growth rate is $\sqrt{\chi K^2(2\mu -\chi K^2)}$, maximized at $K^2=\mu /\chi$ with value $\mu$. Modes beyond the unstable band have real frequencies. The zero-wavenumber phase mode is neutral, with possible linear phase drift from a uniform amplitude perturbation.

Multiply the homogeneous Squire equation by $\bar{\eta }$ and integrate over the channel. The zero boundary values remove the integration by parts boundary term, givingFor real $U$ and $k$, the middle term is purely imaginary. Taking real parts yields the dissipation identity for a homogeneous Squire mode:The strict inequality holds for a nonzero mode with zero endpoint values and positive Reynolds number. Hence all homogeneous Squire modes are temporally stable. In the complete shear-flow perturbation system, the Squire variable can also be forced by the wall-normal velocity; this homogeneous-mode result does not rule out instability of that coupled system or its transient growth.

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.

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.

Assume $\hat{s}>0$, as required for the subcritical case. The zero state exists for every $\mu$. For a nonzero uniform real state, put $B=A_0^2>0$. Its equilibrium point condition is $\mu +3\hat{s}B-10B^2=0$, givingBoth are positive for $-9{\hat{s}}^2/40<\mu <0$. At $\mu =0$, the lower branch reaches zero and the upper has $B=3\hat{s}/10$. For $\mu >0$ only the upper branch is positive. Each positive $B$ represents real states $A_0=\pm \sqrt{B}$, as well as their constant complex wave phase rotations.

The two nonzero branches merge atThe radial linearization is ${\lambda }_a=\mu +9\hat{s}B-50B^2=2B(3\hat{s}-20B)$. It is negative on the upper branch and positive on the lower. The zero state is stable for $\mu <0$ and unstable for $\mu >0$. Thus a subcritical pitchfork bifurcation branch leaves zero towards negative $\mu$, and its fold connects to the stable finite-amplitude branch. Constant wave phase shifts are neutral for nonzero complex states; the stated attraction is in amplitude, with spatial wave phase perturbations damped except for that uniform symmetry mode.

The linearized complex equation couples a perturbation and its conjugate, so a single complex Fourier mode by itself is not a closed perturbation ansatz. Equivalently use real amplitude and wave phase perturbations, with Fourier modes paired to make each real field. Writing $A=A_0+\alpha +i\beta$ and $B=A_0^2$ givesA real amplitude mode of slow wavenumber $\ell >0$ therefore becomes neutral whenFor a nonzero uniform state the equilibrium point condition eliminates $\mu$, leaving amplitude growth $6\hat{s}B-40B^2-4{\ell }^2$. The wave phase mode growth is simply $-4{\ell }^2$, so it cannot give a nonzero-wavenumber bifurcation from these undetuned states.

The concave quadratic has maximumThus a nonzero uniform state needs $\ell \le 3\hat{s}/(4\sqrt{10})$ to have a modulation bifurcation. The slow domain length is ${\epsilon }^{2}L$, so ${\ell }_n=2\pi n/({\epsilon }^{2}L)$. If its smallest positive wavenumber already exceeds this bound, none can become neutral. Since $s={\epsilon }^2\hat{s}$, the conclusion isThis is the modulation cutoff of a cubic-quintic uniform pattern. If the bound is reversed, candidate squared amplitudes are $B=(3\hat{s}\pm \sqrt{9{\hat{s}}^2-160{\ell }^2})/40$. They lie on the lower branch; the upper branch is already damped in amplitude and cannot undergo this stationary modulation bifurcation.

The word “nonzero” is necessary in the final claim. At $A_0=0$, the equilibrium point exists independently of $\mu$, and the mode growth is $\mu -4{\ell }^2$. It becomes neutral at $\mu =4{\ell }_n^2$ for every finite domain. For example, take $\hat{s}=1$, ${\epsilon }^{2}L=20$: then $L=20/{\epsilon }^2<8\pi \sqrt{10}/(3s)$, but the zero state has a first-sideband threshold at $\mu ={\pi }^2/25$. Therefore the printed assertion, if read to include the trivial state, is false. The derived domain-size bound is the correct statement for nonzero pattern branches.

Put $g=\gamma P_0$. For a temporally uniform envelope, $A_z=igA$ and $A_{TT}=0$, so the propagation equation reduces to $-gA+\gamma P_{0}A=0$. Thus the power stays $P_0$ and the accumulated nonlinear wave phase is $gz$. The optical Kerr effect gives an intensity-dependent refractive index, and hence self-phase modulation. For the usual power normalization $\gamma \simeq n_2{\omega }_0/(c{A}_{\mathrm{eff}})$, with effective guided-mode area $A_{\mathrm{eff}}$. Losslessness preserves the background power; it does not remove the self-phase modulation. The carrier wave phase has already been factored out.

Write $A=(\sqrt{P_0}+ϵ)e^{igz}$. To first order,Cancelling the background and its phase derivative leavesThe last term couples a perturbation to its complex conjugate: the pump can transfer energy into a pair of sidebands. It is not an independent wave phase of each Fourier component. This sideband coupling is the linear origin of fibre modulational instability.

For real $K,\Omega ,c_1,c_2$, put $d={\beta }_2{\Omega }^2/2$. Matching the two exponential coefficients givesA nonzero pair requires a zero determinant, henceFor unstable complex $K$, the real-coefficient assumption is too restrictive. The consistent pair is $ϵ=u{e}^{iKz-i\Omega T}+v^{\ast }{e}^{-i{K}^{\ast }z+i\Omega T}$, with complex $u,v$. Its equations give the same determinant and allow the growing and decaying branches. Merely inserting imaginary $K$ while treating the two original coefficients as real would not respect complex conjugation.

For $P_0>0$ and nonzero ${\beta }_2,\gamma$, exponential growth requires ${\beta }_2\gamma <0$ andFor ordinary positive Kerr nonlinearity this means anomalous group velocity dispersion, ${\beta }_2<0$. Define ${\Omega }_c=\sqrt{4g/|{\beta }_2|}$. The growing branch is $K=-iG$, with amplitude growth rateSideband power grows at rate $2G$, rather than $G$. The cutoff and zero-frequency modes have no exponential growth; in a degenerate zero-frequency sector a secular phase response can remain. Positive dispersion with positive $\gamma$ gives real $K$ and oscillatory perturbations.