pith. sign in

arxiv: 2511.17499 · v2 · submitted 2025-11-21 · 🧬 q-bio.PE · cond-mat.dis-nn

Local equations for the generalized Lotka-Volterra model on sparse asymmetric graphs

David Machado , Pietro Valigi , Tommaso Tonolo , Maria Chiara Angelini This is my paper

Pith reviewed 2026-05-17 06:51 UTC · model grok-4.3

classification 🧬 q-bio.PE cond-mat.dis-nn
keywords generalized Lotka-Volterrasparse graphsasymmetric interactionsFokker-Planck equationmean-field closurestationary statesphase diagramecological stability
0
0 comments X

The pith

Local Fokker-Planck equations with mean-field closure compute stationary states for generalized Lotka-Volterra dynamics on sparse asymmetric graphs.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives local Fokker-Planck equations for the generalized Lotka-Volterra model on sparse graphs and closes them with a mean-field approximation to obtain stationary states. This approach works for both symmetric and asymmetric interaction matrices and avoids the computational cost of integrating the full stochastic dynamics. The resulting method reproduces known results for symmetric cases and produces the first phase diagram for asymmetric sparse networks, offering a practical route to stability analysis in ecological communities that match real-world sparsity and directionality.

Core claim

Deriving local Fokker-Planck equations for each node and applying a mean-field closure yields a closed set of ordinary differential equations whose fixed points give the stationary abundances; the same framework maps the phase diagram separating stable and unstable regimes on sparse asymmetric graphs.

What carries the argument

Local Fokker-Planck equations closed by a mean-field approximation that replaces higher-order moments with products of single-node averages.

If this is right

  • Stationary states and stability boundaries become computable for networks with thousands of species without full stochastic integration.
  • The phase diagram for asymmetric sparse interactions shows how the fraction of positive versus negative links controls the existence of a stable fixed point.
  • The same local-equation technique applies directly to models in economics and evolutionary game theory that share the same interaction structure.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could be extended to time-dependent driving or to networks whose topology itself evolves, by updating the local equations at each step.
  • One could test whether the predicted boundaries remain sharp when the graph is generated from empirical food-web data rather than random sparse ensembles.
  • Adding a small amount of noise to the interaction strengths in the closed equations might reveal how robust the phase diagram is to measurement error in real ecosystems.

Load-bearing premise

The mean-field closure remains accurate enough on sparse asymmetric graphs that higher-order correlations do not shift the location of the stationary states by a large amount.

What would settle it

A systematic comparison between the stationary states obtained from the closed equations and those measured in long direct simulations of the stochastic dynamics on the same graphs, for a range of asymmetry strengths and connectivities.

Figures

Figures reproduced from arXiv: 2511.17499 by David Machado, Maria Chiara Angelini, Pietro Valigi, Tommaso Tonolo.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparing individual abundances from IBMF and [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Transitions obtained simulating the gLV model for [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Transitions of the gLV model for different system sizes [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Probability that IBMF does not converge ( [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Predictions of IBMF and BP in the presence of thermal nois [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Effects of damping in the results of IBMF for the [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Transitions of the Generalized Lotka-Volterra model [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Probability that IBMF does not converge ( [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Probability that IBMF does not converge ( [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Average runtime (in seconds) required to obtain [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
read the original abstract

Real ecosystems are characterized by sparse and asymmetric interactions, posing a major challenge to theoretical analysis. We introduce a new method to study the generalized Lotka-Volterra model with stochastic dynamics on sparse graphs. By deriving local Fokker-Planck equations and employing a mean-field closure, we can efficiently compute stationary states for both symmetric and asymmetric interactions. We validate our approach by comparing the results with the direct integration of the dynamical equations and by reproducing known results and, for the first time, we map the phase diagram for sparse asymmetric networks. Our framework provides a versatile tool for exploring stability in realistic ecological communities and can be generalized to applications in different contexts, such as economics and evolutionary game theory.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The manuscript derives local Fokker-Planck equations for the generalized Lotka-Volterra model with stochastic dynamics on sparse graphs and applies a mean-field closure to compute stationary states for both symmetric and asymmetric interactions. It reports validation against direct integration of the dynamical equations, reproduction of known results, and presents a phase diagram for sparse asymmetric networks for the first time.

Significance. If the mean-field closure is accurate, the work supplies an efficient computational framework for stationary-state analysis in realistic sparse and directed ecological networks, extending prior techniques beyond dense or symmetric cases and enabling phase-diagram mapping that could inform stability predictions in ecology and related fields.

major comments (1)
  1. [Validation and results sections] The mean-field closure for local statistics assumes neighbor independence, yet on sparse asymmetric graphs directed paths can generate correlations that alter effective growth rates and stability boundaries. This assumption is load-bearing for the claimed phase diagram; explicit quantitative error bounds (e.g., relative deviation in stationary abundances or boundary locations) versus direct integration must be shown specifically in the low-connectivity, high-asymmetry regime to rule out systematic bias.
minor comments (1)
  1. [Abstract] The abstract states that known results are reproduced but does not identify which results or parameter regimes are used for this check.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive feedback. We address the major comment on validation of the mean-field closure below.

read point-by-point responses

  1. Referee: The mean-field closure for local statistics assumes neighbor independence, yet on sparse asymmetric graphs directed paths can generate correlations that alter effective growth rates and stability boundaries. This assumption is load-bearing for the claimed phase diagram; explicit quantitative error bounds (e.g., relative deviation in stationary abundances or boundary locations) versus direct integration must be shown specifically in the low-connectivity, high-asymmetry regime to rule out systematic bias.

    Authors: We agree that the neighbor-independence assumption underlying the mean-field closure is an approximation and that directed paths on sparse asymmetric graphs could in principle induce correlations capable of shifting effective growth rates or stability boundaries. The manuscript already includes direct comparisons of the mean-field stationary states against numerical integration of the stochastic Lotka-Volterra dynamics for both symmetric and asymmetric sparse graphs, together with reproduction of known limiting cases. Nevertheless, the referee is correct that these comparisons do not yet provide explicit quantitative error bounds focused on the low-connectivity, high-asymmetry regime. In the revised manuscript we will add a dedicated subsection (or supplementary figure) that reports relative deviations in mean stationary abundances and in the locations of phase boundaries, obtained by comparing the closed local Fokker-Planck equations against direct integration at average degrees as low as 3 and across a range of asymmetry parameters. This will allow readers to assess the magnitude of any systematic bias in the reported phase diagram. revision: yes

Circularity Check

0 steps flagged

Derivation of local Fokker-Planck equations plus mean-field closure is independent of fitted inputs or self-citation chains

full rationale

The paper starts from the stochastic generalized Lotka-Volterra dynamics on sparse graphs, derives local Fokker-Planck equations, applies a standard mean-field closure to obtain stationary states, and validates the outputs against direct numerical integration plus reproduction of previously known results before mapping the new asymmetric phase diagram. None of the load-bearing steps reduce by construction to the target quantities (no parameter fitting to the same data being predicted, no self-definitional closure, no uniqueness theorem imported from the authors' prior work). The mean-field step is an explicit approximation whose accuracy is tested externally rather than assumed by definition, making the overall chain self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides insufficient technical detail to enumerate specific free parameters, axioms, or invented entities; the mean-field closure is the primary modeling assumption but its precise form and any fitted quantities are not stated.

pith-pipeline@v0.9.0 · 5421 in / 1114 out tokens · 46814 ms · 2026-05-17T06:51:23.407112+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

  1. [1]

    Epidemic spreading in real networks: an eigenvalue viewpoint

    Wang, Y., Chakrabarti, D., Wang, C., and Faloutsos, C. “Epidemic spreading in real networks: an eigenvalue viewpoint”. In: 22nd International Symposium on Reliable Dis- tributed Systems, 2003. Proceedings. 2003, pp. 25–34. 10.1109/RELDIS.2003.1238052

    work page doi:10.1109/reldis.2003.1238052 2003

  2. [2]

    Epidemic processes in complex networks

    Pastor-Satorras, R., Castellano, C., Van Mieghem, P., and Vespignani, A. “Epidemic processes in complex networks”. Reviews of Modern Physics 87.3 (2015), p. 925. 10. 1103/RevModPhys.87.925. S15

    work page 2015

  3. [3]

    Available: https://link.aps.org/doi/10.1103/PhysRevE

    Cator, E. and Van Mieghem, P. “Second-order mean-field susceptible-infected-susceptible epidemic threshold”. Physical review E 85.5 (2012), p. 056111. 10.1103/PhysRevE. 85.056111

    work page doi:10.1103/physreve 2012

  4. [4]

    Pair quenched mean-field theory for the susceptible- infected-susceptible model on complex networks

    Mata, A. S. and Ferreira, S. C. “Pair quenched mean-field theory for the susceptible- infected-susceptible model on complex networks”. EPL (Europhysics Letters) 103.4 (2013), p. 48003. 10.1209/0295-5075/103/48003

    work page doi:10.1209/0295-5075/103/48003 2013

  5. [5]

    Tonolo, M

    Tonolo, T., Angelini, M. C., Azaele, S., Maritan, A., and Gradenigo, G. General- ized Lotka-Volterra model with sparse interactions: non-Gaussian effects and topological multiple-equilibria phase. 2025. arXiv: 2503.20887 [cond-mat.stat-mech]

    work page arXiv 2025

  6. [6]

    Numerical implementation of dynamical mean field theory for disordered systems: application to the Lotka–Volterra model of ecosystems

    Roy, F., Biroli, G., Bunin, G., and Cammarota, C. “Numerical implementation of dy- namical mean field theory for disordered systems: application to the Lotka–Volterra model of ecosystems”. Journal of Physics A: Mathematical and Theoretical 52.48 (2019), p. 484001. 10.1088/1751-8121/ab1f32

    work page doi:10.1088/1751-8121/ab1f32 2019

  7. [7]

    and Ryzhik, I

    Gradshteyn, I. and Ryzhik, I. Table of Integrals, Series, and Products . Elsevier, 7th Edition, 2007. isbn: 9780123736376

    work page 2007

  8. [8]

    Machado, D., Valigi, P., Tonolo, T., and Angelini, M. C. GenLotkaVolterra_SparseGraphs: GitHub Repository for this article. 2025. 10.5281/zenodo.17612770. S16

    work page doi:10.5281/zenodo.17612770 2025