Abstract:We carry out mathematical analyses, {\em à la} Helmholtz's and Boltzmann's 1884 studies of monocyclic Newtonian dynamics, for the Lotka-Volterra (LV) equation exhibiting predator-prey oscillations. In doing so a novel "thermodynamic theory" of ecology is introduced. An important feature, absent in the classical mechanics, of ecological systems is a natural stochastic population dynamic formulation of which the deterministic equation (e.g., the LV equation studied) is the infinite population limit. Invariant density for the stochastic dynamics plays a central role in the deterministic LV dynamics. We show how the conservation law along a single trajectory extends to incorporate both variations in a model parameter $\alpha$ and in initial conditions: Helmholtz's theorem establishes a broadly valid conservation law in a class of ecological dynamics. We analyze the relationships among mean ecological activeness $\theta$, quantities characterizing dynamic ranges of populations $\mathcal{A}$ and $\alpha$, and the ecological force $F_{\alpha}$. The analyses identify an entire orbit as a stationary ecology, and establish the notion of "equation of ecological states". Studies of the stochastic dynamics with finite populations show the LV equation as the robust, fast cyclic underlying behavior. The mathematical narrative provides a novel way of capturing long-term dynamical behaviors with an emergent {\em conservative ecology}.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to understand the dynamic behavior of ecosystems from the perspective of thermodynamics, especially the oscillation phenomenon in the predator - prey system described by the Lotka - Volterra (LV) equations. Specifically, the author aims to provide a new "thermodynamic theory" for ecology by introducing a generalized conservation law and stochastic dynamics analysis. The following are the core issues of the paper:
1. **Establishing the thermodynamic framework of the ecosystem**: The author attempts to apply the concepts of classical thermodynamics and statistical mechanics to ecosystems, especially by introducing Helmholtz's theorem to explain the dynamic behavior of LV equations. This includes defining concepts such as energy and entropy analogously and exploring the significance of these concepts in ecosystems.
2. **Analyzing the relationship between deterministic and stochastic dynamics**: The author studies the deterministic solutions of LV equations (i.e., the classical differential equation form) and their corresponding stochastic processes (such as birth - death processes). They attempt to reveal the connection between these two descriptions, especially how the deterministic solution becomes the average behavior of the stochastic process in the large - population limit.
3. **Exploring the influence of parameter changes on system behavior**: By introducing the parameter α (representing the ratio of the dynamic ranges of predator and prey populations), the author analyzes the influence of different parameter values on the long - term behavior of the system. They pay particular attention to how changing α affects the average ecological activity θ of the system and other related quantities.
4. **Proposing the equation of state for the ecosystem**: Similar to the equation of state in physics, the author proposes an equation of state that describes the relationships between different variables in the ecosystem. These equations describe the relationships between the geometric characteristics of the ecosystem (such as area A), the parameter α, and new ecological quantities (such as θ).
5. **Revealing the connection between conservative dynamics and stochastic dynamics**: The author shows that conservative dynamics (such as non - dissipative periodic motion) is an important part of stochastic population dynamics. In the large - population limit, a scale separation between fast periodic motion and slow stochastic traversal is observed.
### Formula Summary
- **Energy function \(H(x, y)\)**:
\[
H(x, y)=\alpha x + y-\ln(x^{\alpha}y)
\]
- **Average population range**:
\[
\frac{1}{\tau}\int_{0}^{\tau}(x(t) - 1)^{2}dt=\frac{\hat{A}}{\alpha\tau}, \quad \frac{1}{\tau}\int_{0}^{\tau}(y(t) - 1)^{2}dt=\frac{\alpha\hat{A}}{\tau}
\]
where \(\hat{A}\) is the area enclosed by \(D_{h}(\alpha)\).
- **Helmholtz's theorem**:
\[
dh = \theta(h, \alpha)d\ln A - F_{\alpha}(h, \alpha)d\alpha
\]
where
\[
\theta(h, \alpha)=A_{D_{h}(\alpha)}\left(\frac{\partial A}{\partial h}\right)_{\alpha}^{-1}=\langle\ln x\frac{\partial H}{\partial\ln x}\rangle_{t}=\langle\ln y\frac{\partial H}{\partial\ln y}\rangle_{t}
\]
and
\[
F_{\alpha}(h, \alpha)=\left(\frac{\partial A}{\partial\alpha}\right)_{h}\left(\frac{\partial A}{\partial h}\right)_{\alpha}^{-1}=-\langle\frac{\partial H(x, y, \alpha)}{\partial\alpha}\rangle_{t}