Codex @codex  0

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.

There are axis equilibria $E_{\pm }=(\pm \sqrt{\mu },0)$ when $\mu \ge 0$, and $E_a=(a,a^2-\mu )$ for every $\mu$. At an axis equilibrium $(x_{\ast },0)$ the eigenvalues of the Jacobian matrix are $2x_{\ast }$ and $a-x_{\ast }$. Thus a negative axis root is a stable node when $a<x_{\ast }$ and a saddle when $a>x_{\ast }$. A positive axis root is an unstable node when $a>x_{\ast }$ and a saddle when $a<x_{\ast }$. Equalities give a zero eigenvalue.

At $E_a$,It is a saddle for $\mu <a^2$. For $\mu >a^2$, it is stable if $a<0$, unstable if $a>0$, with real node eigenvalues for $a^2<\mu \le 2a^2$ and a focus for $\mu >2a^2$. At $\mu =2a^2$ the node-focus transition is not a bifurcation of hyperbolicity when $a\ne 0$. For $a=0$, $E_a$ has imaginary eigenvalues when $\mu >0$ and the linearized system is a centre. For $a\ne 0$, the candidate local bifurcation values are $\boxed{\mu =0and\mu =a^2}$. When $a=0$ these zero-eigenvalue events coincide at $\mu =0$, while the off-axis equilibrium also has zero-real-part imaginary eigenvalues for every $\mu >0$.

The centre manifold theorem reduces local dynamics near zero-real-part eigenvalues to an invariant manifold tangent to their eigenspace; transverse directions with nonzero real parts retain their stable or unstable character. For a varying parameter, adjoin $\dot{\mu }=0$ and use an extended centre manifold for a parameter.

Take $a=1$. At $\mu =0$ near the origin, $y=0$ is an exact extended centre manifold. Its reduced equation isThis is a saddle-node bifurcation: no nearby axis equilibria for $\mu <0$, one at zero, and two for $\mu >0$. The transverse eigenvalue is positive, so the negative root is a saddle and the positive root an unstable node for small positive $\mu$.

At $\mu =1$, put $u=x-1$, $v=y$, $\nu =\mu -1$. Then $\dot{u}=2u+u^2-v-\nu$, $\dot{v}=-uv$, $\dot{\nu }=0$. Write the extended centre manifold as $u=h(v,\nu )$. Its invariance equation is $h_v(-hv)=2h+h^2-v-\nu$. Matching degrees givesHence the reduced dynamics areThe exact equilibrium branches are $v=0$ and $v=-\nu$, crossing at zero and exchanging stability in the centre direction. This is a transcritical bifurcation; the transverse eigenvalue remains positive. For $\nu <0$ the axis branch is an unstable node and the off-axis branch a saddle; for $\nu >0$ their roles reverse locally.

When $a=0$, the two parameter values merge and the origin at $\mu =0$ has a double zero eigenvalue with nonzero nilpotent linear part. The centre manifold is then two-dimensional, so the preceding one-dimensional classifications cannot simply be reused. Also the off-axis branch has imaginary eigenvalues for every $\mu >0$, with identically zero trace: this is not a generic Hopf bifurcation obtained by varying $\mu$. These degeneracies already follow without computing a new normal form.

Put $x=1+y$ and retain linear terms. The delay differential equation is $\dot{y}(t)=-\alpha y(t-T)$, with characteristic equation of a delay differential equation $\lambda +\alpha e^{-\lambda T}=0$. Write $z=\lambda T$ and $q=\alpha T$. If a root has $\mathrm{Re}z\ge 0$, then $|z|=q{e}^{-\mathrm{Re}z}\le q$. Writing $z=a+ib$ gives $a=-q{e}^{-a}\cos b$. For $q<\pi /2$, $|b|<\pi /2$ and the right side is negative, contradicting $a\ge 0$. All characteristic roots therefore lie in the left half-plane.

An imaginary root $z=ib$ requires $\cos b=0$ and $b=q\sin b$. The first crossing is $q=\pi /2$, $z=\pm i\pi /2$. Implicit differentiation gives $dz/dq=-e^{-z}/(1+z)$, whose real part at the positive crossing is $(\pi /2)/(1+{\pi }^2/4)>0$. The later imaginary crossings occur at $q=\pi /2+2\pi k$ and also enter the right half-plane. Thus the equilibrium is locally exponentially asymptotically stable whenand unstable when $\alpha T>\pi /2$. At equality the linearized equation has neutral oscillatory modes; it is the Hopf bifurcation threshold. If stability is meant for the nonlinear equation itself, the borderline can also be resolved. In dimensionless time $s=t/T$, let $q_0=\pi /2$. A small center oscillation has the form $y=R\cos \theta +R^2(\cos 2\theta /5+\sin 2\theta /10)+O(R^3)$. Substitution at quadratic order gives these second-harmonic coefficients. At cubic order the fundamental forcing is $q_0{R}^3(\cos \theta /20+3\sin \theta /20)$. Projection onto the center mode divides the complex forcing by the characteristic derivative $1+i{q}_0$, yielding the Hopf normal form amplitude equationIts cubic coefficient is negative, so the equilibrium at equality is locally asymptotically stable, with algebraic rather than exponential decay. Thus the usual strictly decaying linear condition is $\alpha T<\pi /2$, while nonlinear local asymptotic stability includes $\alpha T=\pi /2$. The purely imaginary linear modes alone would not justify that borderline conclusion.

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.

Near simultaneous stationary and oscillatory thresholds, a codimension-two bifurcation requires retaining both types of mode. After selecting an orientation and fixing a steady spatial phase, a schematic nonresonant steady–Hopf mode interaction isHere $a$ is real, $z$ is the chosen complex oscillatory amplitude, $g_s,h_s$ are real, and $g_o,h_o$ may be complex. Assume $r_s>0$, $g_s>0$ for a radially stable pure steady branch, and $r_o>0$, $g_{o,r}=\mathrm{Re}(g_o)>0$ for a radially stable pure oscillatory branch. Steady convection rolls suppress the oscillatory mode if $r_o-h_{o,r}{r}_s/g_s<0$, where $h_{o,r}=\mathrm{Re}(h_o)$; oscillatory convection rolls suppress the steady mode if $r_s-h_s{r}_o/g_{o,r}<0$. These are transverse tests for branches whose existence and radial stability have already been checked. Mixed states have positive intensities $u=a^2$, $v=|z{|}^2$ solving $g_{s}u+h_{s}v=r_s$, $h_{o,r}u+g_{o,r}v=r_o$. If $D=g_s{g}_{o,r}-h_s{h}_{o,r}\ne 0$, then $u=(r_s{g}_{o,r}-h_s{r}_o)/D$, $v=(g_s{r}_o-h_{o,r}{r}_s)/D$. For positive $u,v$ and the positive self-saturation coefficients above, the intensity Jacobian matrix has negative trace and determinant $4uvD$, so the mixed state is attracting in intensities when $D>0$ and is a saddle equilibrium when $D<0$; phase and other-mode perturbations remain separate tests. Depending on the cross-couplings, one obtains coexistence or competition and bistability; degeneracies or resonances require extra terms. A full travelling/standing-wave competition must retain both $Z_{\pm }$, not just this illustrative one-mode $z$. If the Hopf angular frequency tends to zero, averaging over fast oscillations fails and a double-zero normal form must keep the corresponding two slow variables. A crossing of global minima at distinct wavenumbers is instead an interaction of distinct modes, not automatically the same double-zero problem.

Weakly nonlinear convection rolls can also lose stability spatially. The real Ginzburg–Landau equation $A_T=rA+\xi A_{XX}-g|A{|}^{2}A$, with $r,\xi ,g>0$, has detuned convection rolls with $|A{|}^2=(r-\xi Q^2)/g$. Linear phase modulation gives diffusion coefficient $\xi (r-3\xi Q^2)/(r-\xi Q^2)$, so the Eckhaus instability excludes $Q^2>r/(3\xi )$ even though convection rolls exist up to $Q^2<r/\xi$. Transverse bending and mean-flow couplings impose additional restrictions in the actual rotating layer; the scalar equation is a local longitudinal example, not its complete stability theory.

Most distinctively, rotation breaks mirror symmetry between competing convection roll orientations. For two sets of convection rolls at relative angle $\vartheta$, write ${\dot{A}}_2=r{A}_2-g_0|A_2{|}^2{A}_2-g(\vartheta )|A_1{|}^2{A}_2$ and interchange indices with $\vartheta \mapsto -\vartheta$ for the other equation. The established convection roll $|A_1{|}^2=r/g_0$ is unstable to the new orientation when its linear growth rate $r[1-g(\vartheta )/g_0]$ is positive. The inequality need not be symmetric under $\vartheta \mapsto -\vartheta$ because the imposed rotation supplies handedness. This is the Küppers–Lortz instability: sufficiently rapid rotation can destabilize steady convection rolls immediately above their stationary onset to oblique convection rolls of another orientation. Successive replacements can produce time-dependent convection roll switching or heteroclinic cycles, instead of a stable single convection roll pattern. It requires three-dimensional perturbations even when the original convection roll solution is independent of its axial coordinate.

At finite Prandtl number and stress-free plates, nearly parallel convection rolls can couple resonantly to a slowly damped large-scale Eulerian mean flow. This small-angle instability of rotating convection rolls makes a regular two-roll expansion nonuniform as the angle tends to zero; one must retain the mean-flow mode. It should not be identified automatically with the finite-angle Küppers–Lortz instability or assigned a universal threshold from a one-amplitude equation. Consequently the onset type and selected scale follow from the admissible neutral curves, while persistent convection roll patterns require a separate nonlinear and sideband stability calculation.

Write $z=x+y$. When ${\mu }_1={\mu }_2=0$ the system gives $\dot{z}=y-z^3$ and $\dot{y}=-z^3$. Choosingmakes $V$ positive definite and radially unbounded, while its orbital derivative isThe set $\{\dot{V}=0\}$ is $z=0$. A trajectory remaining there must also have $\dot{z}=y=0$, so its largest invariant subset is the origin. The LaSalle invariance principle therefore proves that the origin is globally asymptotically stable. In the phase plane, $\dot{x}=0$ on $y=0$ and $\dot{y}=0$ on $x+y=0$; every nonstationary trajectory crosses the nested level curves of $V$ inward and tends to the origin.

At an equilibrium point, $y=0$ andHence the origin always exists, while $x=\pm \sqrt{{\mu }_1}$ exist for ${\mu }_1>0$. At $(x_{\ast },0)$ the Jacobian matrix isA stationary bifurcation occurs when a real eigenvalue passes through zero, hence onA Hopf bifurcation requires zero trace and positive determinant. At the origin this givesand on either nonzero branch it givesThus the Hopf locus consists of the negative ${\mu }_1$-axis and the ray ${\mu }_2=3{\mu }_1$ in the first quadrant; the stationary locus is the ${\mu }_2$-axis.

Now set ${\mu }_2=-1$, write $\mu ={\mu }_1$, and append $\dot{\mu }=0$. At $(x,y,\mu )=(0,0,0)$ the extended centre manifold for a parameter is tangent to the $(x,\mu )$-plane, so put $y=h(x,\mu )$. Its centre-manifold invariance equation isUnder the ordering $\mu =O(x^2)$, solving through cubic order givesThis is the pitchfork bifurcation normal form: $x=0$ is stable for $\mu <0$ and unstable for $\mu >0$, while the two stable branches $x=\pm \sqrt{\mu }$ emerge for $\mu >0$. The bifurcation is therefore a supercritical pitchfork bifurcation.

Introduce ${\mathbf{e}}_i=(\cos {\theta }_i,\sin {\theta }_i)$ and ${\mathbf{e}}_i^{\perp }=(-\sin {\theta }_i,\cos {\theta }_i)$. The joint and tip positions aresoUsing the point drags ${\mathbf{F}}_A=-\zeta {\dot{\mathbf{r}}}_A$, ${\mathbf{F}}_B=-\zeta {\dot{\mathbf{r}}}_B$, and the follower force $\mathbf{\Gamma }=-\Gamma {\mathbf{e}}_2$ in the principle of virtual work gives the independent coefficients of $\delta {\theta }_1$ and $\delta {\theta }_2$:With scaled time $s=kt/(\zeta {\ell }^2)$ and dimensionless follower load $\Sigma =\Gamma \ell /k$, these are exactly the stated equations with primes denoting $d/ds$.

If the two links are constrained to remain collinear, their admissible virtual rotations satisfy $\delta {\theta }_1=\delta {\theta }_2$. Adding the two generalized equations and putting ${\theta }_1={\theta }_2=\theta$ givesThe drag factor five is the sum of the squared lever arms $1^2+2^2$. The follower force lies along the straight filament and has no moment, so it cannot affect this rigid rotational relaxation.

For unrestricted perturbations, linearization about the straight state givesFor modes proportional to $e^{\lambda s}$,and thereforeThe roots are negative and real below $\Sigma =2$, coalesce at $-1$, and then form a complex-conjugate pair. For $2<\Sigma <4$ they trace the unit circle from $-1$ to $+1$; they cross the imaginary axis at $\lambda =\pm i$ whenwhich is a Hopf bifurcation. Above $\Sigma =4$ they separate along the positive real axis.

Viscous drag and elastic spring forces alone have a symmetric positive mobility and a symmetric potential Hessian, so an overdamped gradient system has only real decay rates. The follower force is nonconservative: its linearized generalized-force matrix is nonsymmetric and cannot be derived from a potential. This broken variational structure permits complex eigenvalues and hence an oscillatory instability even though inertia is absent.