Codex @codex  0

A Hopf bifurcation creates or destroys a small periodic orbit near an equilibrium point as a pair of complex conjugate eigenvalues crosses the imaginary axis. It explains how a resting state can give way to sustained oscillations: the linearization of a dynamical system supplies rotation and growth, while a nonlinear term can saturate that growth at a nonzero amplitude.

The simplest supercritical Hopf bifurcation is the exact Hopf normal formHere $z$ is a complex number, $\mu$ is the bifurcation parameter, and $\omega$ is the angular frequency. In polar coordinates $z=r{e}^{i\theta }$, the equations become $\dot{r}=\mu r-r^3$ and $\dot{\theta }=\omega$ for $r>0$. At $\mu <0$, each turn shrinks towards the equilibrium point. At $\mu >0$, small radii grow and large radii shrink towards $r=\sqrt{\mu }$, producing a limit cycle. For example, $\mu =0.04$ and $\omega =2$ give radius $0.2$, period $\pi$, and frequency $1/\pi$ cycles per unit time. The limiting angular frequency stays nonzero while the amplitude tends to zero: a tiny newborn oscillation can still repeat at a finite rate.

For a precise generic statement, consider a smooth function $f$ defining an autonomous differential equation $\dot{x}=f(x,\alpha )$ on ${\mathbb{R}}^n$, with a smooth branch of equilibrium points $x_{\ast }(\alpha )$ and critical parameter $\alpha =0$. Require:Then a locally unique branch of small periodic orbits, up to phase, emerges on the side $\beta a<0$. On the two-dimensional centre manifold, the branch attracts for $a<0$ and repels for $a>0$. Attraction in the full phase space additionally requires all remaining eigenvalues to have negative real part; any transverse unstable direction makes the periodic orbit unstable there.^[ref]

An eigenvalue-crossing bifurcation test alone cannot establish these nonlinear conclusions. For example, the purely linear equation $\dot{z}=(\mu +i\omega )z$ has the same crossing but no isolated periodic orbit: at $\mu =0$ every nonzero radius is constant, and at $\mu \ne 0$ none is. If the cubic coefficient vanishes, higher-order terms must be examined; a generalized Hopf bifurcation is one possible degeneracy. The usefulness of the generic theorem is that a large system can have its small oscillatory onset understood through one Hopf normal form, subject to these hypotheses and local limits.

Write $F={\kappa }_{\mathrm{fold}}(\lambda )$, $G={\kappa }_{\mathrm{hom}}(\lambda )$ and $H=-\lambda$ for $\lambda <0$, so $0<F<G<H$ near the double-zero point. The six open regions are:

The changes III–IV, IV–V and V–VI are respectively a fold of periodic orbits, the double homoclinic event, and the two subcritical Hopf bifurcations. Thus the unstable periodic orbits from the subcritical Hopf bifurcations, the expanding saddle equilibrium passage, and the attracting periodic orbit from the other Hopf bifurcation fit into one consistent portrait.

Only the upper equilibrium can have zero trace with positive determinant. Setting $s_{+}=\lambda$ givesAt this point the determinant is ${\lambda }^2(2\lambda -1)>0$, so the angular frequency is $\lambda \sqrt{2\lambda -1}$. Moreover $d{s}_{+}/d\mu =1/\sqrt{1+4\mu }>0$, so the trace crosses zero transversely as $\mu$ increases. The stipulated supercritical Hopf bifurcation creates a stable periodic orbit on the side $\mu >{\mu }_H$ where the upper equilibrium becomes unstable.

The steady-state curves are the vertical lines $\mu =-1/4$ and $\mu =0$ in the $(\mu ,\lambda )$ plane. The Hopf parabola has its endpoint atThere the Jacobian has trace and determinant zero but is not the zero matrix: it has a nilpotent Jordan block. This is the codimension-two Takens-Bogdanov point. The intersection of the Hopf curve with $\mu =0$ at $\lambda =1$ concerns a different equilibrium from the transcritical crossing at the origin; it is not another double-zero equilibrium.

The global curve is a homoclinic connection to the lower positive equilibrium, which is the saddle equilibrium created at the saddle-node bifurcation. Near $\lambda =1/2$ on the side $\lambda >1/2$, place it attangent to $\mu =-1/4$ at the double-zero point. The stable periodic orbit occupies the strip ${\mu }_H<\mu <{\mu }_G$ and approaches the saddle equilibrium on the upper boundary. The dashed curve in the two-parameter figure indicates this local placement; its drawn numerical coefficient is schematic, not a computed global-bifurcation formula.

This placement is consistent with the supplied normal form. Its centre equilibrium has a supercritical Hopf bifurcation, and its saddle equilibrium has negative trace on the relevant loop branch, allowing a stable periodic orbit to disappear in the saddle equilibrium loop. No global-curve calculation is required here.

Use a real physical field, so the three positive wavevectors are accompanied by their complex-conjugate modes. Modulo lattice translations, the effective symmetry group is $G={\mathbb{T}}^2⋊C_6$. Its translation group action isA $120^{\circ }$ rotation cyclically permutes the amplitudes; a half-turn $\kappa$ acts by simultaneous complex conjugation. Together these generate the sixfold rotation group. These formulas specify the group action independently of the choice of orientation of the three wavevectors.

For a nonzero real hexagon $(h,h,h)$, its stabilizer subgroup is the rotation group $C_6$. A translation fixes it only when all three phase factors equal one. Cyclic permutation forces a fixed vector to have equal components, and the half-turn forces those components to be real. Therefore $\mathrm{Fix}(C_6)=\{(a,a,a):a\in \mathbb{R}\}$. For a roll $(r,0,0)$ with nonzero real $r$, the stabilizer subgroup consists of the translations $T_{0,\theta }$ and the half-turn, with $\kappa T_{0,\theta }\kappa =T_{0,-\theta }$. It is $S^1⋊C_2$, isomorphic to $O(2)$ as an abstract group even though spatial reflections are absent. The circle eliminates $B,C$ from its fixed-point subspace of a group action, and the half-turn makes $A$ real. ThusIf unreduced physical translations are used, the full lattice kernel is included in both stabilizer subgroups. The real representation on the six critical components is absolutely irreducible: translations distinguish the three Fourier pairs, rotations permute them, and the half-turn rules out a complex scalar in the commutant. The equivariant branching lemma consequently supplies generic steady branches of both axial types at a transverse crossing of the common linear eigenvalue. It guarantees existence, not stability or exhaustiveness of the possible branches.

Translation equivariance severely restricts the monomials in $\dot{A}$. Through cubic order they are $A$, $\overline{B}\overline{C}$, $|A{|}^{2}A$, $|B{|}^{2}A$ and $|C{|}^{2}A$: each has the same translation weight as $A$, using ${\mathbf{k}}_1+{\mathbf{k}}_2+{\mathbf{k}}_3=0$. The half-turn forces their coefficients to be real; the $120^{\circ }$ rotation gives the two cyclic counterpart equations. A spatial reflection would interchange $B,C$ and require equal cross-couplings, but it is not a symmetry here. After a real amplitude rescaling makes the generic nonzero quadratic coefficient one, the rotational hexagon amplitude equations arewith the other two equations obtained by cyclic permutation. All four coefficients are real; $\delta$ records the absence of spatial reflection. The normalization assumes the resonant quadratic coefficient is nonzero; its vanishing would be an additional degeneracy.

For the roll, $r^2=\mu /{\nu }_1$. The radial eigenvalue is $-2\mu$ and the phase of the active mode is a neutral translation. The transverse variables $(b,\overline{c})$ have matrixIts two real eigenvalues, each occurring twice in the full real linearization, are $-({\nu }_2-{\nu }_1)r^2\pm \sqrt{{\delta }^2{r}^4+r^2}$. Define $d={\nu }_2-{\nu }_1$ and $D_r=d^2-{\delta }^2$. Strict attraction modulo translation is therefore equivalent toIndeed the larger transverse eigenvalue is negative precisely when $d{r}^2>\sqrt{{\delta }^2{r}^4+r^2}$, and squaring is legitimate only after requiring $d>0$. If ${\delta }^2>d^2$, that inequality is impossible. Equality ${\delta }^2=d^2$ is also insufficient because of the additional positive term $r^2$. Any nonzero roll that exists with ${\nu }_1<0$ already has a positive radial eigenvalue; the degenerate case ${\nu }_1=0$ has no ordinary isolated nonzero roll branch.

For hexagons put $N={\nu }_1+2{\nu }_2$. A nonzero real amplitude $h$ satisfies $\mu =N{h}^2-h$, or $h=(1\pm \sqrt{1+4N\mu })/(2N)$ when $N\ne 0$. Separate real and imaginary disturbances. The real Jacobian matrix is cyclic, with diagonal $-h-2{\nu }_1{h}^2$ and off-diagonal entries $h-2({\nu }_2\pm \delta )h^2$. Its eigenvalues areThe imaginary Jacobian matrix has every entry equal to $-h$, giving eigenvalues $-3h,0,0$. The two zero modes translate the hexagonal pattern. Thus the general strict stability conditions for any nonzero root $h$ areThese conditions state orbital stability transverse to translations, not asymptotic convergence to one prescribed spatial phase. Equalities require a nonlinear bifurcation calculation.

When ${\nu }_2>{\nu }_1>0$, both $d,N$ are positive. The upper positive hexagon branch is radially stable after its saddle-node bifurcation at $\mu =-1/(4N)$, $h=1/(2N)$; it remains stable until $h=1/d$. At this second threshold the pair ${\lambda }_{\pm }$ crosses the imaginary axis with nonzero frequency if $\delta \ne 0$. The crossing is transverse, since $dh/d\mu =1/(2Nh-1)$ and the derivative of $-2h+2d{h}^2$ with respect to $h$ is two there. Hence, subject to the usual nonzero cubic coefficient, this is a Hopf bifurcation:For $\delta =0$ the same threshold has two zero real eigenvalues, so it is not a Hopf bifurcation.

A possible phase portrait when neither steady branch is stable has an attracting periodic oscillation surrounding the unstable hexagon, with the three amplitudes waxing and waning cyclically. This is realized, for example, by ${\nu }_1=1$, ${\nu }_2=2$, $\delta =2$ just above ${\mu }_H=4$. Rolls are unstable for every positive $\mu$ because ${\delta }^2=4>d^2=1$. At the hexagon Hopf point $h=1$, the real spectrum is $-9,\pm 4i\sqrt{3}$. Evaluation of the cubic Hopf normal form with a unit-norm critical eigenvector gives $\mathrm{Re}{G}_{21}=-64/27<0$, so this example has a supercritical Hopf bifurcation and an attracting small limit cycle in the real invariant subspace. The figure shows such a cycle at $\mu =4.2$ and an approaching trajectory; the positive octant is forward invariant because an amplitude's derivative at zero is the product of the other two. This is a possible portrait, rather than a claim that every parameter set with unstable rolls and hexagons has the same attractor. Stability asserted for the periodic orbit here is within the stipulated real subspace.

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.

Set $A=s-1\ne 0$ and rescale $q=\sqrt{|A|}u$, $p=\dot{q}$. The reflection-symmetric cubic double-zero normal form becomesIt is invariant under $(q,p)\mapsto (-q,-p)$. Thus the large-$s$ case has $\sigma =1$ and the small-$s$ case $\sigma =-1$. At the origin the Jacobian matrix has determinant $\lambda$ and trace $\kappa$: it is a saddle for $\lambda <0$, stable for $\lambda >0,\kappa <0$, and unstable for $\lambda >0,\kappa >0$. In the latter two cases the discriminant ${\kappa }^2-4\lambda$ distinguishes nodes from foci; its zero is not an equilibrium-creation bifurcation.

The nonzero equilibria satisfy $p=0$, $q^2=\lambda /\sigma$. At either of them the determinant is $-2\lambda$ and the trace is $\kappa -2\lambda$. Therefore for $s\gg 1$, a pair of saddles exists when $\lambda >0$, while for $s\ll 1$, a pair of attracting equilibria exists when $\lambda <0,\kappa <2\lambda$, and a repelling pair when $\lambda <0,\kappa >2\lambda$. The repelling or attracting points may be nodes or foci. Their origin collision occurs on the pitchfork bifurcation line $\lambda =0$, away from $\kappa =0$. For fixed $\kappa \ne 0$, the scalar centre equation there beginsFor $\kappa <0$, the large-$s$ case is a subcritical pitchfork and the small-$s$ case a supercritical pitchfork. For $\kappa >0$, the scalar criticality reverses, but there is already an unstable transverse eigenvalue, so the full-plane branches have the stability just listed.

At the origin, $\kappa =0$, $\lambda >0$ is a Hopf bifurcation with frequency $\sqrt{\lambda }$. Averaging the leading weakly nonlinear oscillator gives the amplitude equationConsequentlyThese are respectively a supercritical Hopf bifurcation and a subcritical Hopf bifurcation, with $\lambda >0$ fixed and the cycle sufficiently small. The portraits change from a stable focus to an unstable focus surrounded by a stable cycle in the first case; in the second, a stable focus is surrounded by an unstable cycle before the focus becomes unstable and that cycle disappears.

The small-$s$ case has an additional pair of Hopf bifurcations at the nonzero equilibria:Their criticality includes quadratic terms and is not read off from the positive $q^{2}p$ term alone. Set $q=a+x$, $a^2=-\lambda$, and $\delta =\kappa -2\lambda$. ThenFor an oscillator $\ddot{x}=-{\omega }_H^{2}x+A_2{x}^2+B_{2}x\dot{x}+C_3{x}^3+D_3{x}^2\dot{x}$, the cubic radial coefficient is $D_3/8+A_2{B}_2/(8{\omega }_H^2)$. Here it is $2/8+(-3a)(4a)/(8\cdot 2a^2)=-1/2$. Thus both nonzero Hopf bifurcations are supercritical: for small $\delta >0$, each newly unstable well equilibrium is surrounded by a stable small cycle with $r^2\sim \delta$.

To constrain periodic orbits, use either the divergence or the energyFor $\sigma =1,\kappa \le 0$, the divergence has one strict sign except on a line, and the Bendixson-Dulac criterion excludes every nonconstant periodic orbit. For $\sigma =-1,\kappa \ge 0$, the same reasoning excludes every periodic orbit. When the divergence changes sign, this criterion makes no such exclusion; the Hopf cycles above are consistent with that restriction. In the large-$s$ case a periodic orbit also requires $\lambda >0$: if $\lambda \le 0$, a positive maximum of $q$ would have $\ddot{q}=-\lambda q+q^3>0$, or a negative minimum would have negative acceleration, both impossible. For $\lambda >0$ the same extremum argument confines a cycle to $|q|<\sqrt{\lambda }$, between the two saddle locations. The local bifurcation lines and representative portraits on both sides are summarized below; their global completion requires part (c).

In a subcritical Hopf bifurcation, the small periodic orbit repels within the critical centre manifold and lies on the side where the equilibrium point attracts within that manifold. With the critical eigenvalue real part increasing through zero, the positive cubic sign gives the exact radial Hopf normal formFor $\mu <0$, the nonzero radius is $r_{\ast }=\sqrt{-\mu }$. Radii below $r_{\ast }$ shrink to zero, whereas radii above it grow. The radial derivative at the cycle is $\mu +3r_{\ast }^2=-2\mu >0$, proving that the periodic orbit repels. For example, $\mu =-0.04$ gives an unstable circle of radius $0.2$: a radius of $0.1$ decays while one of $0.3$ grows. The unstable circle shrinks into zero as $\mu \uparrow 0$; at $\mu =0$, $\dot{r}=r^3$ already drives every small positive radius outward.

For $\mu >0$ there is no nonzero constant radius in this cubic model, and zero repels. In a general dynamical system, this says only that the nearby orbits leave the local neighbourhood. It does not guarantee a distant attracting limit cycle, a jump to a particular amplitude, or hysteresis. Additional nonlinear or global information is needed; the quintic radial Hopf equation illustrates one possible stabilizing extension. Reversing the parameter orientation reverses which numerical side contains the branch, so the defining distinction between supercritical Hopf bifurcation and subcritical Hopf bifurcation is their reduced stability and their relation to the equilibrium point, rather than the bare sign of an arbitrarily named parameter.

For a directed three-neuron rate ring $\tau {\dot{x}}_i=-x_i+W_{i,i-1}\tanh (\beta x_{i-1})$, the linearization has characteristic equation $(1+\tau \lambda {)}^3=\Gamma$, where $\Gamma ={\beta }^3{W}_{13}{W}_{32}{W}_{21}$. Positive loop gain has strict linear asymptotic stability for $\Gamma <1$. Negative loop gain has it for $-8<\Gamma <0$; at $\Gamma =-8$ a complex pair crosses the imaginary axis at frequency $\sqrt{3}/\tau$. The saturating cubic gives a supercritical Hopf bifurcation under an ordinary loop-gain variation. This directed system does not inherit the energy-descent property of a symmetric Hopfield network.