ProofComplete

Bifurcation Theory - Key Proof

ProofProof Sketch of the Hopf Bifurcation Theorem

We outline the proof of the Hopf bifurcation theorem, emphasizing the key steps and techniques involved.

Setup: Consider x˙=Fμ(x)\dot{\mathbf{x}} = \mathbf{F}_\mu(\mathbf{x}) with F0(0)=0\mathbf{F}_0(0) = 0 and linearization having eigenvalues λ=±iω0\lambda = \pm i\omega_0 at μ=0\mu = 0. We assume the transversality condition α(0)>0\alpha'(0) > 0 (eigenvalues cross to the right half-plane as μ\mu increases).

Step 1: Complexification and change of variables

Convert to complex notation. Let z=x+iyz = x + iy represent the plane, so the system becomes z˙=λz+g(z,zˉ,μ)\dot{z} = \lambda z + g(z, \bar{z}, \mu) where gg contains nonlinear terms. The goal is to transform this into polar coordinates and analyze the radial dynamics.

Step 2: Center manifold reduction

For μ\mu near 0, the center manifold is two-dimensional (the plane itself in this case). The key is to eliminate non-resonant terms through near-identity changes of variables, leaving only the essential dynamics.

Step 3: Normal form computation

Using the method of normal forms, perform successive near-identity transformations to eliminate as many nonlinear terms as possible. The transformed system, in polar coordinates (r,θ)(r, \theta), takes the form:

r˙=α(μ)r+c(μ)r3+O(r5)\dot{r} = \alpha(\mu) r + c(\mu) r^3 + O(r^5) θ˙=ω(μ)+O(r2)\dot{\theta} = \omega(\mu) + O(r^2)

where α(0)=0\alpha(0) = 0, α(0)>0\alpha'(0) > 0, and c(0)=1c(0) = \ell_1 is the first Lyapunov coefficient.

Step 4: Analysis of the radial equation

For small μ|\mu|, the radial equation r˙=α(μ)r+1r3+O(r5)\dot{r} = \alpha(\mu) r + \ell_1 r^3 + O(r^5) determines stability:

  • If μ<0\mu < 0 (so α(μ)<0\alpha(\mu) < 0): The origin r=0r = 0 is stable; small perturbations decay exponentially.

  • If μ>0\mu > 0 (so α(μ)>0\alpha(\mu) > 0): The origin becomes unstable. If 1<0\ell_1 < 0 (supercritical), there exists a stable periodic orbit at:

r(μ)=α(μ)1=α(0)μ1+O(μ)r^*(\mu) = \sqrt{\frac{-\alpha(\mu)}{\ell_1}} = \sqrt{\frac{-\alpha'(0)\mu}{\ell_1}} + O(\mu)

The orbit is stable because ddr(αr+1r3)r=r=α+31(r)2=α2α=α<0\frac{d}{dr}\left(\alpha r + \ell_1 r^3\right)\bigg|_{r=r^*} = \alpha + 3\ell_1(r^*)^2 = \alpha - 2\alpha = -\alpha < 0 for μ>0\mu > 0.

Step 5: Existence and smoothness

The implicit function theorem guarantees that the periodic orbit r(μ)r^*(\mu) extends to a smooth family for small μ>0\mu > 0. The period is approximately:

T(μ)=2πω(μ)2πω0T(\mu) = \frac{2\pi}{\omega(\mu)} \approx \frac{2\pi}{\omega_0}

which varies smoothly with μ\mu.

Step 6: Return to original coordinates

The periodic orbit in polar coordinates corresponds to a limit cycle in the original (x,y)(x, y) coordinates. The inverse transformations (from normal form back to original system) preserve the existence and stability of the limit cycle, completing the proof.

Conclusion: Under the non-degeneracy and transversality conditions, a Hopf bifurcation creates a branch of periodic orbits bifurcating from the fixed point. The stability of the limit cycle is determined by the sign of the first Lyapunov coefficient 1\ell_1.

The proof combines several powerful techniques: complexification simplifies the eigenvalue problem, normal form theory eliminates inessential nonlinearities, and center manifold reduction focuses on the critical modes. This multi-step approach is paradigmatic for bifurcation analysis.

Remark

The computation of the first Lyapunov coefficient 1\ell_1 involves explicit calculation using third-order Taylor series coefficients of F\mathbf{F}. While algebraically involved, symbolic computation packages can automate this. The coefficient 1\ell_1 determines whether the bifurcation is supercritical (stable limit cycle, 1<0\ell_1 < 0) or subcritical (unstable limit cycle, 1>0\ell_1 > 0). In applications, this distinction is crucial: supercritical bifurcations lead to graceful onset of oscillations, while subcritical ones can cause sudden jumps to large-amplitude behavior.

ExampleBrusselator Lyapunov Coefficient

For the Brusselator x˙=a(b+1)x+x2y\dot{x} = a - (b+1)x + x^2y, y˙=bxx2y\dot{y} = bx - x^2y at the Hopf bifurcation b=1+a2b = 1 + a^2, explicit calculation yields:

1=a24(1+a2)<0\ell_1 = -\frac{a^2}{4}(1 + a^2) < 0

Since 1<0\ell_1 < 0, the Hopf bifurcation is supercritical: a stable limit cycle emerges smoothly as bb increases above 1+a21 + a^2. The amplitude grows as b(1+a2)\sqrt{b - (1 + a^2)}, and the oscillation represents periodic variation in chemical concentrations.

The Hopf bifurcation theorem exemplifies the power of modern dynamical systems theory: qualitative geometric insights (oscillations emerge from eigenvalue crossing) combined with quantitative analysis (normal forms, Lyapunov coefficients) provide complete understanding of local bifurcations. This methodology extends to other bifurcations (saddle-node, pitchfork, period-doubling) and forms the backbone of practical bifurcation analysis in engineering and science.