Mathematical foundations of Lotka-Volterra Phase Plane.
In mathematical biology and nonlinear dynamical systems, the Lotka-Volterra equations represent the foundational archetype of predator-prey coevolution. Formulated by Alfred Lotka in 1925 to model oscillating chemical reactions and independently by Vito Volterra in 1926 to explain Adriatic fishery population cycles, the non-linear coupled system dx/dt = α x - β x y, dy/dt = δ x y - γ y exhibits rich phase-plane geometry. Rather than converging to an asymptotic attractor, the trajectories form a continuous family of closed, periodic concentric orbits encircling a non-hyperbolic center equilibrium. Underpinning this stability is a hidden conservation law: the system possesses a time-invariant first integral V(x, y) = δ x - γ ln(x) + β y - α ln(y) that remains strictly constant along all physical trajectories. In this computational laboratory, we integrate the nonlinear phase plane, trace periodic orbits, analyze center stability, and verify the Hamiltonian invariant constant H = 2.000.
Viva voce questions & analytical theory
How did Alfred Lotka and Vito Volterra independently discover this system?
Alfred Lotka (1925) originally proposed the equations as a model for autocatalytic chemical reactions. Independently, Vito Volterra (1926) formulated the model to explain why predator fish (sharks, rays) in the Adriatic Sea increased in abundance relative to prey fish during World War I when commercial fishing ceased.
How is the conserved first integral V(x, y) derived mathematically?
Dividing the two differential equations eliminates time: dy/dx = (y (δ x - γ)) / (x (α - β y)). Separating variables yields ((α - β y) / y) dy = ((δ x - γ) / x) dx. Integrating both sides yields α ln(y) - β y = δ x - γ ln(x) - C, which rearranges to the conserved quantity V(x, y) = δ x - γ ln(x) + β y - α ln(y) = Constant.
Why are the closed orbits in the Lotka-Volterra model not limit cycles?
A limit cycle is an isolated closed trajectory that attracts or repels neighboring trajectories (as in the Van der Pol oscillator). In Lotka-Volterra, every non-equilibrium initial condition produces a distinct closed orbit depending continuously on the initial value of V(x_0, y_0), forming a continuous family of concentric nested curves like a harmonic oscillator.
What is the Volterra Principle regarding the effect of indiscriminate pesticide/fishing pressure?
Volterra proved that if both prey and predator populations are subjected to uniform proportional harvesting (replacing α by α - ε and γ by γ + ε), the average prey population increases from γ/δ to (γ + ε)/δ, while the average predator population decreases from α/β to (α - ε)/β. Indiscriminate insecticides paradoxically increase the average pest population.
Why does standard explicit Euler integration fail for Lotka-Volterra?
Because the eigenvalues of the Jacobian at (1, 1) are purely imaginary (λ = ± i), the stability amplification factor of Forward Euler is |1 + i h| = √(1 + h²) > 1 for all h > 0. Forward Euler artificially adds energy, causing trajectories to spiral outwards away from the closed orbit, whereas symplectic integrators preserve the invariant.
How can the Lotka-Volterra system be transformed into a canonical Hamiltonian system?
Using the logarithmic coordinate transformation u = ln(x) and v = ln(y), the system transforms into du/dt = -∂H/∂v and dv/dt = ∂H/∂u, where the Hamiltonian is H(u, v) = δ e^u - γ u + β e^v - α v. This proves the system is a true non-canonical Hamiltonian dynamical system.
What are the limitations of the classic Lotka-Volterra model in real-world ecology?
The model assumes unlimited prey growth in the absence of predators (exponential growth rather than logistic carrying capacity), linear predator functional response (predators eat at a rate proportional to prey density without saturation, ignoring handling time), and structural instability (any perturbation changes the orbit permanently).
Algorithm precautions and stability rules
- Use high-order integration (RK4) with small step sizes h ≤ 0.005 to prevent artificial numerical energy dissipation or growth.
- Ensure population states remain strictly in the positive quadrant x > 0, y > 0 where the logarithmic invariant ln(x) is mathematically defined.
- Verify that initial conditions are not chosen exactly at the stationary equilibrium (1, 1), which would trivially yield a single point.
- Monitor the conserved quantity V(t) at every integration step as a continuous sanity check on numerical integration fidelity.
- Ensure floating point precision remains 64-bit IEEE 754 to prevent catastrophic cancellation in the evaluation of logarithmic differences.
- Do not use explicit Euler integration, which artificially injects energy and spirals outward from the closed orbit.
Theoretical limitations
- The classic Lotka-Volterra model is structurally unstable: any small modification (such as adding logistic prey saturation) converts the center into a stable spiral or limit cycle.
- Standard non-symplectic numerical integrators will eventually exhibit secular drift in the invariant V over millions of orbital periods.
- The deterministic model ignores demographic stochasticity and extinction thresholds that occur in small real-world animal populations.