Codex @codex  0

For $\dot{x}=y^3-y$, $\dot{y}=x-y^2$, set $a_{\pm }=(1\pm \sqrt{3})/2$ and $F_{\pm }=x-1-a_{\pm }(y^2-1)$. Direct differentiation gives ${\dot{F}}_{\pm }=-2a_{\pm }y{F}_{\pm }$, proving the displayed first integral. The invariant parabolas meet at saddles $(1,\pm 1)$ and bound a heteroclinic cycle. Inside that lens, the logarithm of $I$ has a strict nondegenerate maximum at the origin, proving a nonlinear center equilibrium, rather than inferring one solely from imaginary Jacobian eigenvalues.

For small positive forcing-frame frequency $\omega$, use $\mu =\hat{\mu }{\omega }^2$, $C=\omega z$ and $\tau =\omega T$ in the three-to-one spatially forced amplitude equation. The leading system is $z_{\tau }+iz={\overline{z}}^2$. Writing $z=x+iy$ gives the Hamiltonian system $(x_{\tau },y_{\tau })=(H_y/2,-H_x/2)$ and the displayed first integral. The origin is a center equilibrium; three saddle equilibria at $(0,1)$ and $(\pm \sqrt{3}/2,-1/2)$ lie on $H=1/3$. The factorization $H-1/3=(1+2y)[x^2-(y-1{)}^2/3]$ reveals a triangular heteroclinic cycle, containing closed periodic orbits for every $0<H<1/3$.

The equilibria are $0$ and $-\epsilon \pm \sqrt{{\epsilon }^2+\lambda }$. Their linear trace is $\kappa -v_{\ast }^2$ and determinant $\lambda -4\epsilon v_{\ast }-3v_{\ast }^2$. For nonzero $\epsilon$ the symmetry-protected pitchfork bifurcation splits into a transcritical bifurcation at $\lambda =0$ and a saddle-node bifurcation at $\lambda =-{\epsilon }^2$. The energy $H={\dot{v}}^2/2+\lambda v^2/2-2\epsilon v^3/3-v^4/4$ obeys $\dot{H}=(\kappa -v^2){\dot{v}}^2$. For positive $\lambda$, the two outer saddle equilibria have unequal barrier heights, so the symmetric heteroclinic cycle unfolds into separated one-way connections and a homoclinic orbit termination of the central stable limit cycle.

The fixed-point conditions are $y(y^2-1)=0$ and $x=y^2$. They give $(0,0)$ and $(1,\pm 1)$. The Jacobian matrix isAt $(1,1)$ its eigenvalues are $-1\pm \sqrt{3}$; at $(1,-1)$ they are $1\pm \sqrt{3}$. Both equilibria are saddle equilibria. At the origin the eigenvalues are $\pm i$. This linear calculation alone does not distinguish a nonlinear center from a weak focus.

To settle the origin and the global connections, setUsing $2a_{\pm }^2-2a_{\pm }-1=0$, direct differentiation gives ${\dot{F}}_{\pm }=-2a_{\pm }y{F}_{\pm }$. Hence each parabola $F_{\pm }=0$ is invariant and, away from them, the following is a first integral:The two parabolas intersect precisely at the saddle equilibria. Their arcs between $y=-1$ and $y=1$ enclose the lensInside it $F_{+}>0$, $F_{-}<0$. The origin is a critical point of $\log I$ with negative-definite Hessian $-2\sqrt{3}$ times the identity. It is therefore a strict local maximum of this first integral, proving a nonlinear center equilibrium. The regular interior levels are closed periodic orbits, circulating anticlockwise. In fact the only critical point of $I$ inside the lens is the origin, so these levels fill the lens up to its boundary.

On $F_{+}=0$, $\dot{y}=(a_{+}-1)(y^2-1)$, so the interior arc runs from the upper saddle equilibrium to the lower. On $F_{-}=0$ it runs from the lower to the upper. Together they form a heteroclinic cycle. Their outer branches supply the remaining stable and unstable separatrices.

For the full exterior portrait, the two invariant parabolas partition the rest of the plane into four regions. Where both $F_{\pm }>0$, $\dot{y}>0$ everywhere; where both are negative, $\dot{y}<0$. In the region above both saddle equilibria, $y>1$ and $\dot{x}>0$; in the corresponding lower region, $y<-1$ and $\dot{x}<0$. These exterior trajectories are unbounded and cannot cross the separatrices. The nullclines are $y=0,\pm 1$ for $\dot{x}=0$ and $x=y^2$ for $\dot{y}=0$.

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.

For the supplied normal form, the four local topological phase portraits are as follows. When ${\mu }_2<0$, only the origin is present and it is a saddle equilibrium. When ${\mu }_2>0,{\mu }_1<0$, the origin is attracting, between the two outer saddle equilibria; inward-facing separatrices enter its basin. When ${\mu }_2>0$ and $0<{\mu }_1<h({\mu }_2)$, the origin is repelling and a stable limit cycle surrounds it, inside the separatrix boundary formed by the outer saddle equilibria. Finally, when ${\mu }_2>0,{\mu }_1>h({\mu }_2)$, the origin remains repelling but the local limit cycle has disappeared, and trajectories leave the central region through the opened separatrix channels. Whether the origin is a node or a focus within one of these regions does not change this orbit topology.

On the connection curve, the two outer saddle equilibria are joined in both directions; the period of the disappearing limit cycle diverges. This is a heteroclinic cycle, not a pair of saddle homoclinic orbits. The sketches below show the four open regions of the supplied system; they are not presented as phase portraits of the inequivalent printed spherical equations.

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.

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.

At $(\sigma ,\mu )=(1,0)$ the origin's matrix is $\left(\begin{array}{cc}0& 0\\ 2& 0\end{array}\right)$: it has a double zero eigenvalue and one eigenvector. Trace and determinant vary independently with $\mu$ and $\sigma$ there, so two parameters are needed. The symmetry of an odd function removes quadratic terms, giving a reflection-symmetric cubic double-zero normal form. The quadratic nonlinearity required for a generic Bogdanov–Takens bifurcation vanishes.

Use a prime for differentiation with respect to the rescaled time. Direct substitution, with all rescaled variables written without tildes, givesAt $\epsilon =0$, the equations are $u^{\prime }=-sv+v^3/2$, $v^{\prime }=2u$. Fordifferentiation gives $H^{\prime }=2u(-sv+v^3/2)+(sv-v^3/2)2u=0$. This is the Hamiltonian blow-up of a reflection-symmetric double-zero point.

For $s>0$, the equilibrium points are a centre equilibrium at $(0,0)$ and saddles at $(0,\pm \sqrt{2s})$. The centre equilibrium has eigenvalues $\pm i\sqrt{2s}$; at either saddle equilibrium they are $\pm 2\sqrt{s}$. Levels $0<H<s^2/2$ contain closed inner ovals around the centre equilibrium. The saddle equilibrium level hasThese two branches form a heteroclinic cycle connecting the saddles. The branch with $u>0$ runs from the lower saddle equilibrium to the upper one because $v^{\prime }=2u>0$; the other runs back. Outside the potential well, trajectories are open. The right panel of the preceding figure uses the coordinates $(v,u)$ to display the well and these directed connections.

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.

Local equilibrium and Hopf bifurcations do not locate the larger cycles' endpoints. In the large-$s$ case, the stable origin cycle is confined between two symmetry-related saddles. A natural global completion is a heteroclinic cycle formed by the two separatrix connections between those saddles. The stable periodic orbit expands into that connecting cycle and disappears; the two connections occur together because $(q,p)\mapsto (-q,-p)$ exchanges them. This is a global transition in the positive-$\kappa$, positive-$\lambda$ region where the Bendixson-Dulac criterion permits cycles. Beyond it, the cubic inverted potential allows trajectories to leave the local region; the local reduction does not specify their ultimate fate in the original four-dimensional system.

In the small-$s$ case, the two stable cycles born around the nonzero well equilibria can encounter the central saddle simultaneously. Their limiting orbits form a symmetric figure-eight pair of homoclinic orbits. A symmetric homoclinic gluing bifurcation can then replace the two one-well cycles by one larger cycle visiting both wells. There is also an outer unstable cycle associated with the subcritical origin Hopf branch, continued into the double-well region. A saddle-node bifurcation of periodic orbits between the large stable and unstable cycles provides a possible later termination. Thus the global curves needed for the usual completion are a symmetric heteroclinic connection in the large-$s$ case, and paired homoclinic gluing together with a fold of cycles in the small-$s$ case.

These connections, and any resulting basin changes, are not detected by the equilibrium Jacobians. Their precise parameter curves require a global return-map or separatrix calculation, which is not part of this requested qualitative description. Symmetry forces paired connections where indicated, but the local centre manifold and Hopf calculations alone should not be presented as numerical locations or a proof of the full large-amplitude diagram.

For $\omega \ne 0$, put $\mu =\hat{\mu }{\omega }^2$, $C=\omega z$ and $\tau =\omega T$, where $z=\hat{C}$. Then $C_T={\omega }^2{z}_{\tau }$ and division by ${\omega }^2$ givesWe may take $\omega >0$ without loss of the dynamics by conjugating the original equation if necessary; otherwise $\tau$ reverses the time orientation. The scaling is singular at $\omega =0$ and is not a transformation for that exactly zero-frequency case.

The Hamiltonian limit of three-to-one forcing drops the terms proportional to $\omega$. For $z=x+iy$ the remaining real system isThe proposed first integral isIndeed $x_{\tau }=H_y/2$ and $y_{\tau }=-H_x/2$, hence $H_{\tau }=0$. The unperturbed equation is a planar Hamiltonian system, with a center equilibrium at the origin and three saddle equilibria atAll saddle equilibria have $H=1/3$. The factorizationshows that their central separatrix consists of the three sides of an equilateral triangle. Each level $0<H<1/3$ inside this triangle is a closed periodic orbit. Indeed, inside the triangle $0<r<1$ and ${\partial }_{r}H=2r(1+r\sin 3\theta )>0$. Each ray from the origin therefore meets each such level once, giving a compact simple closed contour with no equilibrium point on it. The nonzero vector field traverses this contour periodically; the period grows without bound as the separatrix is approached. This supplies an infinite family, not a claim that every level outside the central region is closed.

Restore the small radial perturbation. Its exact effect on the first integral isConsequently the continuum of Hamiltonian system orbits generally does not persist. The origin becomes a weak attracting focus for $\hat{\mu }<0$ or a repelling focus for $\hat{\mu }>0$, and the three hyperbolic saddle equilibria persist with perturbed stable manifold and unstable manifold. For a small positive $\hat{\mu }$, outward drift on very small orbits balances cubic damping on somewhat larger ones, selecting a stable limit cycle rather than an arbitrary energy level. Near the center equilibrium $H\simeq r^2$, so its leading radius is $r\simeq \sqrt{\hat{\mu }}$ when $\hat{\mu }$ is also small.

For a more general closed unperturbed orbit ${\Gamma }_h$, the averaged area criterion for perturbed Hamiltonian cycles says that persistence requires the averaged energy drift ${\oint }_{{\Gamma }_h}{H}_{\tau }d\tau$ to vanish. Since the unperturbed speed is $|\nabla H|/2$, the planar divergence theorem converts this leading drift towhere ${\mathcal{D}}_h$ is the enclosed region. Isolated zeros select candidate periodic orbits; a drift changing from positive inside to negative outside gives an attracting limit cycle. The separatrix triangle has mean $r^2=1/4$, so its leading flux changes sign at $\hat{\mu }=1/2$. This marks the leading possible heteroclinic transition, with higher-order corrections needed to locate it precisely.

As a cycle approaches the saddle equilibria, long residence times and splitting of the heteroclinic cycle become important. Orbits can instead drift inward to the equilibrium point at the origin or leave the periodic island and approach one of the stable states with phase locking of the full canonical equation. Those upper-branch threefold phase-locked equilibria have $C=O(1)$, so they lie outside the local $C=O(\omega )$ scaling. Thus the small perturbation gives energy selection, attracting or repelling oscillations, and possible switching/locking transitions; it does not preserve a conserved $H$ or an infinite family of neutral periodic solutions. This qualitative picture does not assume all global parameter values have the same attractor.

Rotating Rayleigh-Bénard convection combines buoyancy-driven instability with the Coriolis force. Consider a plane layer of depth $d$ rotating uniformly about the vertical axis, heated from below. Adopt the Boussinesq approximation, fixed boundary temperatures and, for explicit formulas, impermeable stress-free boundary conditions. The conductive state of Rayleigh-Bénard convection is motionless with a linear temperature profile. The dimensionless controls are the Rayleigh number, Prandtl number and Taylor number:Here $\nu$ is kinematic viscosity and $\kappa$ thermal diffusivity. In thermal-diffusion time units, linear perturbations satisfyRotation does no direct mechanical work, since $\mathbf{u}\cdot ({\mathbf{e}}_z\times \mathbf{u})=0$, but couples vertical motion to vertical vorticity and changes the damping and oscillation balance.

For horizontal wavenumber $k$ and vertical mode $n$, put $m=n^2{\pi }^2$ and $a=k^2+m$. Use $w=W\sin (n\pi z)$, $\theta =\Theta \sin (n\pi z)$ and vertical vorticity $\zeta =Z\cos (n\pi z)$ with growth rate $\sigma$. Curling the momentum equation and eliminating pressure givesTheir determinant, without dividing by a possibly zero factor, is the rotating-convection growth-rate polynomialThis makes the linear stability analysis question precise: onset occurs when a root reaches zero real part and all other modes still decay.

A stationary neutral root has $\sigma =0$, givingRotation raises this stationary threshold. An oscillatory neutral root has $\sigma =i\varpi$ with $\varpi \ne 0$. Real-imaginary separation gives the oscillatory neutral curve of rotating convectionThis branch is admissible only when ${\varpi }^2>0$, requiring $P<1$ and sufficiently strong rotation. The restoring Coriolis force coupling permits an inertial/thermal oscillation whose phase-lagged buoyancy can overcome dissipation. At large $P$, temperature and momentum diffusion do not permit that overstability mechanism at primary onset, so the exchange of stabilities is stationary. The actual threshold is the minimum of the stationary and admissible oscillatory curves over all allowed modes, not an arbitrary formal value of $R_o$. Rigid plates require a different vertical eigenproblem and Ekman layers, so the explicit free-slip numbers are not universal.

At large Taylor number, the first vertical mode is selected in the ideal plane layer. Let $t=k^2,m={\pi }^2$. Minimizing the stationary curve givesThus the stationary neutral curve of rotating convection hasThe physical horizontal wavelength is $2\pi d/k_c$, hence decreases as ${\mathrm{Ta}}^{-1/6}$; its prefactor depends on the boundary convention. Thin nearly vertical cells reconcile the strong Coriolis force constraint with viscous and thermal diffusion. The oscillatory minimization replaces the right side of the wavenumber equation by $P^2\mathrm{Ta}m/(1+P{)}^2$, giving the same Taylor number exponent at fixed positive $P$. Where its frequency remains admissible,Equality is $8P^4-P-1=0$, whose positive root is approximately $P_{\ast }=0.6766$. Accordingly, for sufficiently rapid rotation in this free-slip problem, $P<P_{\ast }$ selects oscillatory onset and $P>P_{\ast }$ stationary onset. The weaker condition $P<1$ is only necessary for an oscillatory neutral mode; it does not by itself identify the first instability. Finite Taylor number, finite lateral geometry, allowed discrete wave numbers and plate conditions change the selection.

For the counterpropagating Hopf amplitudes in rotating convection near a simple oscillatory onset, the Hopf bifurcation produces slow complex amplitudes for counterpropagating roll waves. After separating the fast carrier oscillation, symmetry permits the cubic equationswith generally complex coefficients; an $i{\omega }_0{Z}_{\pm }$ term restores the fast frequency if desired. The real parts govern amplitude saturation and the imaginary parts give nonlinear frequency shifts. Write $a_s=\mathrm{Re}{g}_s$, $a_c=\mathrm{Re}{g}_c$. For a travelling wave from a supercritical bifurcation with only one amplitude nonzero, $|Z_{\pm }{|}^2=r/a_s$ with $r,a_s>0$, and the competing wave's growth rate is $r(1-a_c/a_s)$. It is amplitude-stable against that competitor when $a_c>a_s$. A standing wave has equal intensities $r/(a_s+a_c)$; provided this is positive, its intensity-difference mode is stable when $a_s>a_c$. Temporal and spatial phase symmetries leave neutral phase directions, so these are orbital/amplitude stability statements, not decay of every phase displacement.

These coefficients follow from nonlinear interactions and the Fredholm solvability condition obtained by projection onto an adjoint eigenfunction; symmetry alone cannot decide their signs. A negative saturating coefficient gives subcritical bifurcation behavior requiring higher-order terms. Spatial modulation leads to coupled complex Ginzburg–Landau equations with group velocities and diffusion; phase instabilities, mean-flow coupling and differently oriented rolls can destabilize a wave stable in the restricted two-amplitude system. A weakly nonlinear expansion of oscillations therefore predicts travelling waves or standing waves, frequency shifts, modulation and possible secondary mode competition, not a unique universal periodic state.

The Küppers–Lortz instability is a different route to time dependence: it destabilizes steady saturated rolls against oblique roll perturbations. For stationary-roll amplitudes of orientations ${\varphi }_j$, a leading competition system hasA pure roll has $|A_i{|}^2=r/g_0$ with $g_0>0$. An infinitesimal new roll at relative angle $\vartheta$ grows atFor sufficiently strong rotation in appropriate boundary and Prandtl number regimes, some finite oblique angle has $g(\vartheta )<g_0$, so a steady roll is unstable arbitrarily close above its stationary onset. Rotation is handed and allows $g(\vartheta )\ne g(-\vartheta )$, so replacement of one roll by another can favor a definite cyclic sense. Three or more competing orientations can form a heteroclinic cycle; whether it attracts depends on contraction/expansion rates and other modes. Noise, spatially varying domains and modulation can turn this competition into repeated orientation switching and irregular patterns.

The invading rolls are three-dimensional disturbances even when the original straight roll is described by a two-dimensional section. The finite-angle Küppers–Lortz instability mechanism should also be distinguished from the small-angle instability of rotating convection rolls mediated by large-scale mean flow at finite Prandtl number. Numerical thresholds and favored angles depend on mechanical boundaries and material parameters; the essential criterion is the cross-coupling relative to self-saturation. Rotation both changes primary onset and wavelength, and can prevent the resulting steady roll pattern from remaining a stable nonlinear state.

Consider a horizontal fluid layer under the Boussinesq approximation heated from below, rotating uniformly about its vertical axis. To make the linear formulas definite, take depth one, stress-free impermeable fixed-temperature plates, a horizontally infinite or sufficiently large periodic domain, and negligible centrifugal-buoyancy modifications. Rigid plates change the vertical eigenfunctions and thresholds and require their own boundary value problem; the explicit formulas below use the free-slip model. Let $P=\nu /{\kappa }_T$ be the Prandtl number, $\mathrm{Ta}=(2\Omega d^2/\nu {)}^2$ the Taylor number, and $R$ the Rayleigh number. Time is again measured on the thermal diffusion time.

The linearized rotating Rayleigh-Bénard convection equations about conduction areThe Coriolis force does no work on the fluid, but it couples vertical motion to vertical vorticity and thereby changes the balance of buoyancy and viscous damping. On the plates require $w=\theta =0$ and ${\partial }_z{u}_x={\partial }_z{u}_y=0$. Write $\zeta =(\nabla \times \mathbf{u}{)}_z$. Taking a curl and a double curl removes pressure and givesA normal mode with horizontal wavenumber $k>0$, vertical index $n\ge 1$ and complex growth rate $s$ has $w=W\sin (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$, $\theta =\Theta \sin (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$ and $\zeta =Z\cos (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$. Put $a=k^2+n^2{\pi }^2$. The three equations become $a(s/P+a)W+\sqrt{\mathrm{Ta}}n\pi Z=R{k}^2\Theta$, $(s/P+a)Z=\sqrt{\mathrm{Ta}}n\pi W$, and $(s+a)\Theta =W$. Their determinant gives the rotating-convection growth-rate polynomialThis cubic includes viscous and thermal modes as well as the convective instability. It avoids excluding a root by division by $s/P+a$. Linear stability means every root has negative real part; stationary onset has $s=0$, while overstability has $s=\pm i\omega$ with nonzero angular frequency.

Setting $s=0$ gives the stationary neutral curve of rotating convectionFor oscillatory convection, set $s=i\omega$ and equate real and imaginary parts. Eliminating $R$ yieldsThe oscillatory neutral curve of rotating convection is physically admissible only when ${\omega }^2>0$. Thus $P<1$ and sufficiently rapid rotation are necessary; a formal minimum with ${\omega }^2\le 0$ is not a Hopf bifurcation. Convection begins at the smaller of the stationary minimum and the admissible oscillatory minimum. The angular frequency approaches zero where a fixed-wave-number stationary and oscillatory threshold meet; that is a double-zero limit of the cubic, requiring a different slow-time reduction.

For wavenumber selection in rotating convection, the lowest vertical index is $n=1$. To compare indices, set $k^2=n^2{\pi }^{2}y$ at fixed $y>0$. Both neutral thresholds have a positive term proportional to $n^4$ and a rotation term independent of $n$, while the admissible angular frequency squared decreases with $n$. Their minima therefore cannot improve on $n=1$. Set $t=k^2$, $m={\pi }^2$. Minimizing each neutral curve givesFor the latter equation also check ${\omega }^2>0$, or minimize over the admissible set instead. Without rotation, $k_c=\pi /\sqrt{2}$ and $R_c=27{\pi }^4/4$. At rapid rotation, $k_s\sim (\mathrm{Ta}{\pi }^2/2{)}^{1/6}$ and $R_{s,c}\sim 3(\mathrm{Ta}{\pi }^2/2{)}^{2/3}$: rotation narrows the horizontal convection rolls and raises their threshold. In a finite box only discrete wavenumbers are permitted, so the minimum must be taken over those allowed modes.

The rapid-rotation comparison has $R_{o,c}/R_{s,c}\sim 2P^{4/3}/(1+P{)}^{1/3}$. Its equality gives $8P^4-P-1=0$, with positive root about $0.67660$. Below this value, sufficiently rapid rotation can make oscillatory convection the primary instability; low Prandtl number allows inertial motions to interact with the thermal field before viscosity damps them. This estimate describes the selected neutral-curve comparison in the stated free-slip asymptotic model, rather than claiming every $P<1$ has oscillatory primary onset at every rotation rate.

Linear selection does not determine the nonlinear planform or its stability. A weakly nonlinear expansion with a solvability condition gives amplitude equations. A stationary convection roll amplitude has the Landau amplitude equation $A_T=rA-g|A{|}^{2}A$; a supercritical branch requires $g>0$ and has $|A{|}^2=r/g$, whereas $g<0$ requires higher-order saturation and permits subcritical behaviour. The coefficient depends on the physical parameters and boundary conditions, and should not be assumed positive merely because the linear threshold is known.

A nonzero-frequency Hopf bifurcation supplies oppositely travelling convection roll amplitudes $Z_{+},Z_{-}$. At a fixed orientation, a cubic normal form iswith complex $g,h$ and the understood symmetry exchanging propagation directions. For $r>0$ and $g_r=\mathrm{Re}g>0$, a travelling-roll branch has one nonzero amplitude, of squared modulus $r/g_r$; it is stable to the opposite travelling amplitude when $h_r=\mathrm{Re}h>g_r$. A standing-roll branch has equal squared moduli $r/(g_r+h_r)$, and is stable within this pair when $g_r>h_r$ and $g_r+h_r>0$. The imaginary parts shift angular frequencies. These conditions concern amplitude perturbations in that reduced subspace; orientation and modulation modes still have to be tested. This is the competition of travelling and standing convection rolls.