Feasibility and Single Parameter Scaling of Extinctions in Large Ecological Communities
Pith reviewed 2026-05-21 20:23 UTC · model grok-4.3
The pith
In large ecological communities, feasibility is broken before stability, with extinctions following a single-parameter scaling law.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Multispecies ecosystems modelled by generalized Lotka-Volterra equations exhibit stationary population abundances. Using standard methods of random matrix theory, distributions of species abundances are Gaussian at equilibrium in the weakly interacting regime. One consequence is that feasibility is generically broken before stability for large enough number of species. An analytical expression for the probability that n=0,1,2,... species go extinct is derived and a single-parameter scaling law is conjectured to govern species extinctions.
What carries the argument
Gaussian distributions of equilibrium species abundances obtained from random matrix theory, which determines when feasibility is lost relative to stability.
If this is right
- Feasibility breaks before stability as the number of species increases.
- The probability of exactly n species extinctions has a closed-form analytical expression.
- Extinctions are controlled by a single scaling parameter.
- These predictions hold across a wide range of system parameters as confirmed by simulations.
Where Pith is reading between the lines
- This could mean that larger communities are inherently more prone to extinctions due to feasibility constraints alone.
- The scaling law might be testable in empirical data from real ecosystems by measuring how extinction rates scale with community size.
- Similar approaches could be applied to other dynamical systems beyond ecology to predict when feasibility fails.
Load-bearing premise
The system stays in the weakly interacting regime so that random matrix methods produce Gaussian abundance distributions at equilibrium.
What would settle it
Numerical simulations or analytical checks where, for large species counts, stability is lost at lower interaction strengths than feasibility, contradicting the generic breaking of feasibility first.
Figures
read the original abstract
Multispecies ecosystems modelled by generalized Lotka-Volterra equations exhibit stationary population abundances, where large number of species often coexist. Understanding the precise conditions under which this is at all feasible and what triggers species extinctions is a key, outstanding problem in theoretical ecology. Using standard methods of random matrix theory, I show that distributions of species abundances are Gaussian at equilibrium, in the weakly interacting regime. One consequence is that feasibility is generically broken before stability, for large enough number of species. I further derive an analytical expression for the probability that $n=0,1,2,...$ species go extinct and conjecture that a single-parameter scaling law governs species extinctions. These results are corroborated by numerical simulations in a wide range of system parameters.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript applies standard random matrix theory methods to the generalized Lotka-Volterra equations modeling large ecological communities. It establishes that species abundance distributions at equilibrium are Gaussian in the weakly interacting regime. This implies that feasibility is generically broken before stability for large numbers of species. The paper derives an analytical expression for the probability that a given number n of species go extinct and conjectures a single-parameter scaling law for species extinctions, which is supported by numerical simulations over a wide range of system parameters.
Significance. If the results are confirmed, this work advances the understanding of extinction dynamics in diverse ecosystems by linking random matrix theory to feasibility and stability. The analytical probability expression and the conjectured scaling law offer new tools for predicting biodiversity patterns in large communities. The corroboration with simulations enhances the reliability of the findings in theoretical ecology.
major comments (2)
- [Abstract] The statement that Gaussian distributions follow from standard random-matrix methods and that numerical simulations corroborate the scaling conjecture lacks the full derivation or error analysis, making it challenging to evaluate if post-hoc parameter choices influence the central claim that feasibility is broken before stability.
- [paragraph on distributions of species abundances] The Gaussian claim relies on the weakly interacting regime where the equilibrium solution remains well inside the positive orthant. However, since the probability that n species have negative abundances grows with N, the effective interaction matrix changes, which may invalidate the RMT-derived Gaussian form and the single-parameter scaling conjecture near the feasibility boundary.
minor comments (1)
- Clarify the notation for the interaction matrix and growth rates to ensure consistency throughout the manuscript.
Simulated Author's Rebuttal
We thank the referee for their thoughtful and constructive comments on our manuscript. We address each major comment below and indicate the revisions we will make to strengthen the presentation and clarify the scope of our results.
read point-by-point responses
-
Referee: [Abstract] The statement that Gaussian distributions follow from standard random-matrix methods and that numerical simulations corroborate the scaling conjecture lacks the full derivation or error analysis, making it challenging to evaluate if post-hoc parameter choices influence the central claim that feasibility is broken before stability.
Authors: The derivation of the Gaussian form for species abundances is obtained by applying standard random-matrix techniques to the equilibrium condition of the generalized Lotka-Volterra system in the large-N limit; the explicit steps appear in the section titled 'Distributions of species abundances.' We have revised the abstract to include a concise outline of these steps and added a new paragraph in the main text that quantifies the approximation error via direct comparison of the analytic Gaussian to the empirical histograms over the full range of interaction strengths examined. The parameter sweeps reported in the figures were performed on a regular grid rather than selected post hoc; the revised manuscript now states this explicitly and includes the corresponding error bars. revision: yes
-
Referee: [paragraph on distributions of species abundances] The Gaussian claim relies on the weakly interacting regime where the equilibrium solution remains well inside the positive orthant. However, since the probability that n species have negative abundances grows with N, the effective interaction matrix changes, which may invalidate the RMT-derived Gaussian form and the single-parameter scaling conjecture near the feasibility boundary.
Authors: We agree that the Gaussian derivation is strictly valid only when the probability of negative abundances remains small. As N grows, this probability increases and extinctions modify the effective interaction matrix. In the revised manuscript we have added an explicit discussion of the regime of validity, together with additional simulations that approach the feasibility boundary from within the weakly interacting regime. These tests show that the single-parameter scaling continues to describe the data well until the fraction of extinctions becomes O(1), at which point the conjecture is no longer expected to hold; we now state this limitation clearly. revision: yes
Circularity Check
No circularity: standard RMT derivation with explicit conjecture and numerical checks
full rationale
The paper applies standard random-matrix methods to obtain Gaussian abundance distributions in the stated weakly interacting regime, then derives an analytical extinction probability expression directly from that distribution. The single-parameter scaling is explicitly labeled a conjecture and corroborated by independent numerical simulations across parameter ranges rather than obtained by fitting a subset and relabeling the output. No load-bearing self-citations, self-definitional steps, or ansatz smuggling are present in the provided derivation chain; the central claims remain independent of the target results and rest on externally standard RMT techniques plus direct simulation validation.
Axiom & Free-Parameter Ledger
free parameters (1)
- interaction strength threshold
axioms (1)
- domain assumption Species abundances at equilibrium follow a Gaussian distribution when interactions are weak.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
distributions of species abundances are Gaussian at equilibrium... kurtosis K=3+O(σ^8,S^{-1})... single-parameter scaling law governs species extinctions
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
feasibility is generically broken before stability, for large enough number of species
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
- [1]
- [2]
- [3]
-
[4]
V. Zelevinsky, B.A. Brown, N. Frazier, M. Horoi, Phys. Rep.276, 85 (1996)
work page 1996
- [5]
- [6]
-
[7]
G. Montambaux, D. Poilblanc, J. Bellissard, and C. Sire, Phys. Rev. Lett.70, 497 (1993)
work page 1993
- [8]
-
[9]
P. Kos, M. Ljubotina, and T. Prosen, Phys. Rev. X8, 021062 (2018)
work page 2018
- [10]
-
[11]
V.V. Flambaum, A.A. Gribakin, G.F. Gribakina, and M.G. Kozlov, Phys. Rev. A50, 267 (1994). [13]Theoretical Ecology: Principles and Applications, edited by R. May and A.R. McLean (Oxford University Press, 2007)
work page 1994
- [12]
- [13]
- [14]
-
[15]
G.F. Fussmann, S.P. Ellner, K.W. Shertzer, and N.G. Hairston, Science290, 1358 (2000)
work page 2000
-
[16]
E.A. McGehee, N. Schutt, D.A. Vasquez, and E.Peacock- L´ opez, Intl. J. Bif. and Chaos18, 2223 (2008)
work page 2008
- [17]
- [18]
- [19]
- [20]
-
[21]
Y.V. Fyodorov and B.A. Khoruzhenko, Proc. Natl. Acad. Sc.113, 6827 (2016)
work page 2016
- [22]
- [23]
- [25]
-
[26]
M. Barbier, J.-F. Arnoldi, G. Bunin, and M. Loreau, Proc. Proc. Natl. Acad. Sc.115, 2156 (2018)
work page 2018
- [27]
- [28]
-
[29]
V. Ros, F. Roy, G. Biroli, G. Bunin, and A.M. Turner, Phys. Rev. Lett.130, 257401 (2023)
work page 2023
- [30]
- [31]
- [32]
- [33]
- [34]
- [35]
- [36]
- [37]
-
[38]
I. A.Hatton, O.Mazzarisi, A. Altieri, and M. Smerlak, Science383, eadg8488 (2024)
work page 2024
-
[39]
G. Aguad´ e-Gorgori´ o, I. Lajaaiti, J.-F. Arnoldi, and S. K´ efi, Oikos e10980 (2025)
work page 2025
- [40]
- [41]
- [42]
- [43]
- [44]
- [45]
-
[46]
V. Volterra, Mem. Acad. Lincei Roma.2, 31 (1926); en- glish translation :Animal Ecology, edited by R.N. Chap- man, McGraw-Hill (1931)
work page 1926
-
[47]
For a simple model with a random distribution ofA ij = ±1, the ratio of predator-prey pairs withA ij =−A ji is (1−γ)/2
-
[48]
F. Roy, M. Barbier, G. Biroli, and G. Bunin, PLoS Com- put Biol16, e1007827 (2020)
work page 2020
-
[49]
H.J. Sommers, A. Crisanti, H. Sompolinsky, and Y. Stein, Phys. Rev. Lett.60, 1895 (1988)
work page 1988
-
[50]
Durrett,Probability: Theory and Examples, 5th Ed., Cambridge University Press (2019)
R. Durrett,Probability: Theory and Examples, 5th Ed., Cambridge University Press (2019). 6 Single Parameter Scaling of Species Extinction in Generalized Lotka-V olterra Systems Supplemental Material In this Supplemental Material section I briefly sketch the two Neumann expansions I use to solve Eq (3). The general method is the following. Inserting the Ne...
work page 2019
-
[51]
with⟨δA ij⟩= 0 and⟨δA ij δAkl⟩=σ 2(δikδjl +γδ ilδjk)/S
First expansion The first Neumann expansion takesX=1andY=A=M+δA, withM ij =µ/S, andδAcontaining the zero-average fluctuations of interspecies interactions, i.e. with⟨δA ij⟩= 0 and⟨δA ij δAkl⟩=σ 2(δikδjl +γδ ilδjk)/S. One has N=κ X p X j (−1)p D [(M+δA) p]ij E ,(S11) All terms with even powers ofδAcontribute to the RMT average, as before, and forµ̸= 0, one...
-
[52]
Second expansion The second Neumann expansion takesX=1+MandY=δAin Eq. (S10). I obtain, instead of Eqs. (S4) and (S11), N= κ 1 +µ ∞X p=0 X j (−1)p Dh δA 1− M 1 +µ ip ij E ,(S15) where I used (1+M) −1 = (1−M/(1 +µ). Taking into account terms up top= 4, I obtain N=κ (1 +µ) −1 + (1 +µ) −2γσ 2 + [(1 +µ) −2 + (1 +µ) −3]γ2σ4 +O(σ 6, S−1).(S16) Remarkably, I reco...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.