TheoremComplete

The Bendixson-Dulac Criterion

The Bendixson-Dulac criterion provides a sufficient condition for the non-existence of periodic orbits, complementing the Poincare-Bendixson existence theorem.


Statement

Theorem6.6Bendixson-Dulac theorem

Let D⊆R2D \subseteq \mathbb{R}^2 be a simply connected domain and F=(f,g)∈C1(D)\mathbf{F} = (f, g) \in C^1(D). If there exists a C1C^1 function B:D→RB: D \to \mathbb{R} (called a Dulac function) such that

∂(Bf)∂x+∂(Bg)∂y\frac{\partial(Bf)}{\partial x} + \frac{\partial(Bg)}{\partial y}

does not change sign in DD and is not identically zero on any open subset, then the system x′=f(x,y)x' = f(x,y), y′=g(x,y)y' = g(x,y) has no periodic orbits entirely contained in DD.


Proof

Proof

Suppose for contradiction that γ\gamma is a periodic orbit in DD. Let RR be the region enclosed by γ\gamma (which exists since DD is simply connected). By Green's theorem:

∬R[∂(Bf)∂x+∂(Bg)∂y]dA=∮γ(Bf dy−Bg dx).\iint_R \left[\frac{\partial(Bf)}{\partial x} + \frac{\partial(Bg)}{\partial y}\right]dA = \oint_\gamma (Bf\,dy - Bg\,dx).

On γ\gamma: dx=f dtdx = f\,dt and dy=g dtdy = g\,dt, so

∮γ(Bf dy−Bg dx)=∮γB(fg−gf) dt=0.\oint_\gamma (Bf\,dy - Bg\,dx) = \oint_\gamma B(fg - gf)\,dt = 0.

But the double integral is nonzero (since the integrand has one sign and is not identically zero). Contradiction. ■\blacksquare

■

Examples

ExampleSimple application

For x′=−y+x(x2+y2−1)x' = -y + x(x^2+y^2-1), y′=x+y(x2+y2−1)y' = x + y(x^2+y^2-1): take B=1B = 1.

∂f∂x+∂g∂y=(3x2+y2−1)+(x2+3y2−1)=4(x2+y2)−2.\frac{\partial f}{\partial x} + \frac{\partial g}{\partial y} = (3x^2+y^2-1) + (x^2+3y^2-1) = 4(x^2+y^2) - 2.

This changes sign (positive for x2+y2>1/2x^2+y^2 > 1/2, negative for x2+y2<1/2x^2+y^2 < 1/2), so B=1B = 1 does not work directly. A different Dulac function is needed.

However, in polar coordinates: r′=r(r2−1)r' = r(r^2-1) has a limit cycle at r=1r = 1, so the criterion correctly cannot rule it out.

ExampleChoosing a Dulac function

For x′=yx' = y, y′=−x−y+x2y' = -x - y + x^2: try B(x,y)=e2yB(x,y) = e^{2y}.

∂(Be2y⋅y)∂x+∂(Be2y(−x−y+x2))∂y=0+e2y[2(−x−y+x2)+(−1)]=e2y(2x2−2x−2y−1).\frac{\partial(Be^{2y}\cdot y)}{\partial x} + \frac{\partial(Be^{2y}(-x-y+x^2))}{\partial y} = 0 + e^{2y}[2(-x-y+x^2) + (-1)] = e^{2y}(2x^2-2x-2y-1).

For the region where 2x2−2x−2y−1<02x^2-2x-2y-1 < 0, no periodic orbits exist.

RemarkGradient and Hamiltonian systems

Gradient systems x′=−∇V(x)\mathbf{x}' = -\nabla V(\mathbf{x}) never have periodic orbits (along any trajectory, V˙=−∣∇V∣2≤0\dot{V} = -|\nabla V|^2 \leq 0, so VV decreases and cannot return to its initial value). Hamiltonian systems x′=∂H/∂yx' = \partial H/\partial y, y′=−∂H/∂xy' = -\partial H/\partial x satisfy div(F)=0\text{div}(\mathbf{F}) = 0 (Bendixson with B=1B=1), but this does not rule out periodic orbits — instead, Hamiltonian systems typically have families of periodic orbits.