Entropy production of cyclic population dynamics

Entropy serves as a central observable in equilibrium thermodynamics. However, many biological and ecological systems operate far from thermal equilibrium. Here we show that entropy production can characterize the behavior of such nonequilibrium systems. To this end we calculate the entropy production for a population model that displays nonequilibrium behavior resulting from cyclic competition. At a critical point the dynamics exhibits a transition from large, limit-cycle like oscillations to small, erratic oscillations. We show that the entropy production peaks very close to the critical point and tends to zero upon deviating from it. We further provide analytical methods for computing the entropy production which agree excellently with numerical simulations.

💡 Research Summary

The paper investigates how entropy production, a central quantity in equilibrium thermodynamics, can be used to characterize far‑from‑equilibrium biological and ecological systems. The authors focus on a classic cyclic competition model involving three species that suppress each other in a rock‑paper‑scissors fashion. By formulating the dynamics as a continuous‑time Markov process, they write down the master equation for the probability distribution over population states and define the steady‑state probability currents (J_{ij}=P_iW_{ij}-P_jW_{ji}), where (W_{ij}) are the transition rates. Entropy production is then expressed as (\sigma=\sum_{i<j}J_{ij}\ln(W_{ij}/W_{ji})), which vanishes when detailed balance holds and becomes positive whenever a net cyclic flow persists.

Two control parameters are central: the system size (N) (total number of individuals) and a competition asymmetry (\varepsilon) that quantifies how much one direction of the cyclic interaction is favored over the opposite direction. In the deterministic limit ((N\to\infty)) the model reduces to ordinary rate equations that either settle to a fixed point (for symmetric interactions) or generate a stable limit cycle (for sufficiently strong asymmetry). When stochastic fluctuations are retained (finite (N)), the authors find a sharp transition at a critical asymmetry (\varepsilon_c\approx 1/N). For (\varepsilon<\varepsilon_c) the system exhibits small‑amplitude, erratic oscillations; for (\varepsilon>\varepsilon_c) it displays large, almost deterministic limit‑cycle‑like oscillations.

The main result is that the entropy production (\sigma) peaks very close to (\varepsilon_c). Near the critical point the probability currents become maximally asymmetric, leading to a pronounced increase in (\sigma). Using a combination of linear noise approximation and system‑size expansion, the authors derive an analytical scaling law (\sigma\propto(\varepsilon-\varepsilon_c)^2) in the vicinity of the transition. This prediction is validated by extensive Gillespie simulations, which show an excellent quantitative agreement with the theory.

The paper also discusses two limiting cases where (\sigma) vanishes. When the competition is completely one‑sided (one species always wins), the dynamics become effectively irreversible but end in an absorbing state with no circulating currents, so (\sigma=0). Conversely, when the interaction rates are perfectly symmetric, detailed balance holds, the steady state is an equilibrium distribution, and again (\sigma=0). These limits illustrate that entropy production directly measures the degree of nonequilibrium cyclic flow in the system.

Beyond the specific three‑species model, the authors argue that entropy production can serve as a universal diagnostic for nonequilibrium phase transitions in ecological networks, predator‑prey systems, and other complex biological assemblies. By quantifying how far a system deviates from detailed balance, (\sigma) provides insight into stability, susceptibility to perturbations, and the underlying mechanisms that sustain biodiversity in fluctuating environments. The work thus bridges concepts from stochastic thermodynamics with ecological dynamics, offering a new quantitative tool for the analysis of far‑from‑equilibrium biological systems.