Emergence of cooperation in nonlinear higher-order public goods games
Pith reviewed 2026-05-10 17:23 UTC · model grok-4.3
The pith
Nonlinear public goods games on hypergraphs allow mixed-order groups to sustain cooperation alongside bistability.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In well-mixed populations, single-order nonlinear public goods games display a change in the nature of the transition from continuous to discontinuous depending on the exact form of nonlinearity. By contrast, mixed-order public goods games, where games with different numbers of players coexist, exhibit a richer dynamical regime wherein a state of active coexistence of bistability and cooperation can arise. On scale-free hypergraphs, cooperation is promoted, with the initial placement of cooperators and the presence of hyperdegree correlations playing crucial roles.
What carries the argument
Nonlinear collective reinforcement in public goods games on hypergraphs, which encodes synergistic or discounting effects for each additional cooperator across groups of varying sizes.
If this is right
- Mixed-order games produce a stable regime in which bistability and cooperation coexist dynamically, unlike the simpler transitions seen in single-order games.
- Scale-free hypergraphs increase overall cooperation levels compared with uniform random hypergraphs.
- The long-term level of cooperation depends on the precise starting distribution of cooperators in the hypergraph.
- Hyperdegree correlations between groups alter the spread of cooperation under nonlinear payoffs.
Where Pith is reading between the lines
- Real social systems whose groups vary in size and overlap may support cooperation more readily than uniform-size models predict.
- Evolutionary models of multiplayer interactions could gain accuracy by representing structure as hypergraphs rather than ordinary networks.
- Interventions to encourage cooperation might deliberately create networks with scale-free properties and diverse group sizes.
- Laboratory experiments that allow participants to form groups of different sizes could test whether the predicted coexistence of bistability and cooperation appears in practice.
Load-bearing premise
The specific mathematical form chosen for the nonlinear payoff reinforcement and the assumption that a static hypergraph plus fixed initial cooperator locations fully capture the relevant interactions without mechanisms such as strategy-dependent group formation.
What would settle it
Numerical simulations of the replicator dynamics on a scale-free hypergraph with mixed-order games that show whether removing hyperdegree correlations eliminates the reported increase in cooperation relative to random hypergraphs.
Figures
read the original abstract
Evolutionary game theory has provided substantial contributions to explain the emergence of cooperation under unfavourable conditions in ecology, economics, and the social sciences. Recently, inspired by newly available empirical evidence on group interactions, higher-order networks have emerged as a natural framework to properly encode multiplayer games in structured populations. Here, we study the emergence of cooperation in a nonlinear public goods game (PGG) on hypergraphs, where collective reinforcement captures the synergistic or discounting effect associated with each additional cooperator. In well-mixed populations, single-order PGGs, where all games have the same number of players, display a change in the nature of transition from continuous to discontinuous depending on the exact form of nonlinearity. By contrast, mixed-order PGGs, where games with different number of players coexist, exhibit a richer dynamical regime wherein a state of active coexistence of bistability and cooperation can arise. We further find that scale-free hypergraphs promote cooperation, highlighting the crucial role played by both the initial placement of cooperators and the presence of hyperdegree correlations. Overall, our results provide a comprehensive characterization of nonlinear PGGs on hypergraphs and open up new avenues for richer models of evolutionary dynamics of multiplayer interactions on structured populations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the emergence of cooperation in nonlinear public goods games (PGGs) played on hypergraphs. It reports that single-order PGGs (fixed group size) exhibit a transition from continuous to discontinuous phase transitions depending on the form of nonlinearity in collective reinforcement, whereas mixed-order PGGs (coexisting group sizes) produce a richer regime featuring active coexistence of bistability and cooperation. The work further claims that scale-free hypergraphs promote cooperation and that both the initial placement of cooperators and hyperdegree correlations play crucial roles.
Significance. If the reported simulation results hold under broader conditions, the manuscript offers a useful extension of evolutionary game theory to higher-order interactions by showing how mixing interaction orders can qualitatively alter dynamical regimes beyond what single-order models predict. The explicit treatment of synergistic/discounting effects and hypergraph topology (including correlations) is a constructive step toward more realistic models of multiplayer cooperation in structured populations.
major comments (2)
- [Abstract] Abstract and results on mixed-order dynamics: the central claim that mixed-order PGGs generically produce an 'active coexistence of bistability and cooperation' rests on simulations using one specific family of nonlinear reinforcement functions (synergistic or discounting). No tests with qualitatively different nonlinearities (e.g., threshold or saturating forms outside the reported parametrization) are shown, so it remains unclear whether the richer regime is a structural consequence of mixing orders or an artifact of the chosen functional family.
- [Methods] Methods and numerical results: the reported transitions and promotion effects on scale-free hypergraphs depend on specific hypergraph generation rules and initial cooperator placement. Without explicit verification of numerical stability across random seeds, hypergraph ensemble sizes, or alternative generation algorithms (e.g., configuration model vs. preferential attachment variants), it is difficult to assess whether the claimed distinctions between single- and mixed-order regimes are robust.
minor comments (2)
- [Model definition] Notation for the nonlinearity exponent and synergy/discount factor should be introduced once with a clear equation reference and then used consistently; occasional redefinition risks confusion when comparing single- and mixed-order cases.
- [Figures] Figure captions for the phase diagrams should explicitly state the hypergraph parameters (e.g., power-law exponent, correlation strength) and the number of independent realizations averaged, to allow direct reproduction of the reported bistability regions.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. We address each major point below, indicating the revisions we will make to strengthen the manuscript.
read point-by-point responses
-
Referee: [Abstract] Abstract and results on mixed-order dynamics: the central claim that mixed-order PGGs generically produce an 'active coexistence of bistability and cooperation' rests on simulations using one specific family of nonlinear reinforcement functions (synergistic or discounting). No tests with qualitatively different nonlinearities (e.g., threshold or saturating forms outside the reported parametrization) are shown, so it remains unclear whether the richer regime is a structural consequence of mixing orders or an artifact of the chosen functional family.
Authors: We acknowledge the referee's concern regarding the scope of the nonlinearity. Our power-law reinforcement function (with exponent r) is chosen because it provides a continuous parametrization that spans the linear case (r=1) and the synergistic (r>1) and discounting (r<1) regimes commonly studied in the PGG literature. The active coexistence in mixed-order settings emerges from the interplay between heterogeneous group sizes and this nonlinearity, rather than from a narrow functional choice. Nevertheless, we agree that the claim of genericity would benefit from explicit qualification. In the revised manuscript we will expand the discussion in the abstract and main text to clarify the scope of the functional family, note that qualitatively different forms (e.g., sharp thresholds) may produce different quantitative boundaries, and suggest this as an avenue for future work. revision: partial
-
Referee: [Methods] Methods and numerical results: the reported transitions and promotion effects on scale-free hypergraphs depend on specific hypergraph generation rules and initial cooperator placement. Without explicit verification of numerical stability across random seeds, hypergraph ensemble sizes, or alternative generation algorithms (e.g., configuration model vs. preferential attachment variants), it is difficult to assess whether the claimed distinctions between single- and mixed-order regimes are robust.
Authors: We agree that explicit robustness information improves the manuscript. Our simulations employed the configuration-model hypergraph generator to fix the hyperdegree distribution while controlling for correlations, and all reported quantities were averaged over 100 independent hypergraph realizations with at least 50 independent initial-condition runs per realization (both uniform random and hub-biased cooperator placement). We will revise the Methods section to state these ensemble sizes, the number of random seeds, and the convergence criteria explicitly. While we did not compare every alternative generator, the promotion of cooperation on scale-free structures is driven by the heavy-tailed hyperdegree distribution itself; we will add a short paragraph noting that the qualitative distinction between single- and mixed-order regimes is expected to persist under other scale-free constructions, but that a systematic comparison lies beyond the present scope. revision: yes
Circularity Check
No significant circularity; results from explicit forward simulation of a defined model.
full rationale
The paper's claims rest on numerical simulations of evolutionary game dynamics on hypergraphs using an explicitly stated nonlinear payoff function for collective reinforcement. No equations or results reduce by construction to fitted parameters, self-definitions, or load-bearing self-citations. The functional form of nonlinearity is presented as a modeling choice whose consequences are explored, not derived from the outcomes themselves. The derivation chain is self-contained against the stated assumptions and simulation protocol.
Axiom & Free-Parameter Ledger
free parameters (2)
- nonlinearity exponent or synergy/discount factor
- hypergraph generation parameters (e.g., degree distribution, correlation strength)
axioms (2)
- domain assumption Population evolves according to imitation or birth-death rules on the hypergraph structure.
- domain assumption Hypergraphs accurately encode higher-order group interactions without additional temporal or strategy-dependent rewiring.
Reference graph
Works this paper leans on
-
[1]
M. A. Nowak and R. Highfield,SuperCooperators: altru- ism, evolution, and why we need each other to succeed, 1st ed. (Free Press, New York, NY, 2012)
work page 2012
-
[2]
N. A. Christakis and J. H. Fowler, Friendship and nat- ural selection, Proceedings of the National Academy of Sciences111, 10796 (2014)
work page 2014
-
[3]
N. A. Christakis,Blueprint: The evolutionary origins of a good society(Little, Brown Spark, 2019)
work page 2019
-
[4]
J. Von Neumann and O. Morgenstern,Theory of games and economic behavior(Princeton University Press, 1944)
work page 1944
-
[5]
J. M. Smith and G. R. Price, The Logic of Animal Con- flict, Nature246, 15 (1973)
work page 1973
-
[6]
C. Hauert and G. Szabó, Game theory and physics, American Journal of Physics73, 405 (2005)
work page 2005
-
[7]
M. Perc, J. J. Jordan, D. G. Rand, Z. Wang, S. Boc- caletti, and A. Szolnoki, Statistical physics of human co- operation, Physics Reports687, 1 (2017)
work page 2017
-
[8]
J. F. Nash Jr, Equilibrium points in n-person games, Proceedings of the National Academy of Sciences36, 48 (1950)
work page 1950
-
[9]
G. Hardin, The Tragedy of the Commons: The popu- lation problem has no technical solution; it requires a fundamental extension in morality., Science162, 1243 (1968)
work page 1968
-
[10]
J. Maynard Smith,Evolution and the theory of games (Cambridge University Press, Cambridge ; New York, 1982)
work page 1982
-
[11]
R.Dawkins,The selfish gene,30thed.(OxfordUniversity Press, Oxford ; New York, 2006)
work page 2006
-
[12]
P. D. Taylor and L. B. Jonker, Evolutionary stable strate- gies and game dynamics, Mathematical Biosciences40, 145 (1978)
work page 1978
-
[13]
J. Hofbauer and K. Sigmund,Evolutionary Games and Population Dynamics, 1st ed. (Cambridge University Press, 1998)
work page 1998
-
[14]
R. Albert and A.-L. Barabási, Statistical mechanics of complex networks, Reviews of Modern Physics74, 47 (2002)
work page 2002
-
[15]
Newman,Networks(Oxford University Press, 2018)
M. Newman,Networks(Oxford University Press, 2018)
work page 2018
-
[16]
S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D.- U. Hwang, Complex networks: Structure and dynamics, Physics Reports424, 175 (2006)
work page 2006
-
[17]
Berge,Hypergraphs: combinatorics of finite sets, Vol
C. Berge,Hypergraphs: combinatorics of finite sets, Vol. 45 (Elsevier, 1984)
work page 1984
-
[18]
F. Battiston, G. Cencetti, I. Iacopini, V. Latora, M. Lu- cas, A. Patania, J.-G. Young, and G. Petri, Networks beyond pairwise interactions: Structure and dynamics, Physics Reports874, 1 (2020)
work page 2020
-
[19]
F. Battiston, E. Amico, A. Barrat, G. Bianconi, G. Fer- raz de Arruda, B. Franceschiello, I. Iacopini, S. Kéfi, V. Latora, Y. Moreno,et al., The physics of higher-order interactions in complex systems, Nature physics17, 1093 (2021)
work page 2021
-
[20]
F. Battiston, V. Capraro, F. Karimi, S. Lehmann, A. B. Migliano, O. Sadekar, A. Sánchez, and M. Perc, Higher- order interactions shape collective human behaviour, Na- ture Human Behaviour , 1 (2025)
work page 2025
-
[21]
F. Battiston, C. Bick, M. Lucas, A. P. Millán, P. S. Skardal, and Y. Zhang, Collective dynamics on higher- order networks, Nature Reviews Physics , 1 (2026)
work page 2026
-
[22]
R. Pastor-Satorras and A. Vespignani, Epidemic spread- ing in scale-free networks, Physical Review Letters86, 3200 (2001)
work page 2001
-
[23]
F. C. Santos and J. M. Pacheco, Scale-free networks pro- vide a unifying framework for the emergence of coopera- tion, Physical Review Letters95, 098104 (2005)
work page 2005
-
[24]
F. C. Santos, J. M. Pacheco, and T. Lenaerts, Evolution- ary dynamics of social dilemmas in structured heteroge- neous populations, Proceedings of the National Academy of Sciences103, 3490 (2006)
work page 2006
-
[25]
M. A. Nowak and R. M. May, Evolutionary games and spatial chaos, Nature359, 826 (1992)
work page 1992
-
[26]
G. Szabó and G. Fath, Evolutionary games on graphs, Physics Reports446, 97 (2007)
work page 2007
-
[27]
M. A. Nowak, Five rules for the evolution of cooperation, Science314, 1560 (2006)
work page 2006
-
[28]
E. Lieberman, C. Hauert, and M. A. Nowak, Evolution- ary dynamics on graphs, Nature433, 312 (2005)
work page 2005
-
[29]
H. Ohtsuki, C. Hauert, E. Lieberman, and M. A. Nowak, A simple rule for the evolution of cooperation on graphs and social networks, Nature441, 502 (2006)
work page 2006
- [30]
-
[31]
F. C. Santos, M. D. Santos, and J. M. Pacheco, Social diversitypromotestheemergenceofcooperationinpublic goods games, Nature454, 213 (2008)
work page 2008
-
[32]
M. Perc, J. Gómez-Gardenes, A. Szolnoki, L. M. Floría, and Y. Moreno, Evolutionary dynamics of group interac- tions on structured populations: a review, Journal of the Royal Society Interface10(2013)
work page 2013
- [33]
-
[34]
U. Alvarez-Rodriguez, F. Battiston, G. F. de Arruda, Y. Moreno, M. Perc, and V. Latora, Evolutionary dy- namics of higher-order interactions in social networks, Nature Human Behaviour5, 586 (2021)
work page 2021
-
[35]
A. Civilini, O. Sadekar, F. Battiston, J. Gómez- Gardeñes, and V. Latora, Explosive cooperation in so- cial dilemmas on higher-order networks, Physical Review Letters132, 167401 (2024). 9
work page 2024
- [36]
-
[37]
H. Guo, D. Jia, I. Sendiña-Nadal, M. Zhang, Z. Wang, X. Li, K. Alfaro-Bittner, Y. Moreno, and S. Boccaletti, Evolutionary games on simplicial complexes, Chaos, Soli- tons & Fractals150, 111103 (2021)
work page 2021
-
[38]
M. Gao, Z. Li, T. Wu, and L. Wang, Evolutionary dy- namics of multiplayer ultimatum games on hypergraphs, Physical Review E111, 054305 (2025)
work page 2025
- [39]
-
[40]
O. Sadekar, A. Civilini, V. Latora, and F. Battiston, Drivers of cooperation in social dilemmas on higher- order networks, Journal of the Royal Society Interface 22, 20250134 (2025)
work page 2025
- [41]
-
[42]
D. Wang, P. Yi, Y. Hong, J. Chen, and G. Yan, Emer- gence of cooperation promoted by higher-order strat- egy updates, PLOS Computational Biology21, e1012891 (2025)
work page 2025
-
[43]
Y. Wang and S. Gao, Evolution of cooperation on hy- pergraphs with heterogeneous update dynamics, Chaos, Solitons & Fractals204, 117776 (2026)
work page 2026
-
[44]
X. Wang, L. Zhou, A. McAvoy, Z. Tian, and A. Li, Strat- egy evolution on temporal hypergraphs, Proceedings of the National Academy of Sciences123, e2516380123 (2026)
work page 2026
-
[45]
I. Iacopini, G. Petri, A. Barrat, and V. Latora, Simplicial models of social contagion, Nature Communications10, 2485 (2019)
work page 2019
-
[46]
P. S. Skardal and A. Arenas, Abrupt desynchronization and extensive multistability in globally coupled oscillator simplexes, Physical Review Letters122, 248301 (2019)
work page 2019
-
[47]
L. Neuhäuser, A. Mellor, and R. Lambiotte, Multibody interactions and nonlinear consensus dynamics on net- worked systems, Physical Review E101, 032310 (2020)
work page 2020
-
[48]
T. Robiglio, L. Di Gaetano, A. Altieri, G. Petri, and F. Battiston, Higher-order ising model on hypergraphs, Physical Review E112, L022301 (2025)
work page 2025
-
[49]
H. Pérez-Martínez, S. Lamata-Otín, F. Malizia, L. M. Floría, J. Gómez-Gardeñes, and D. Soriano-Paños, So- cial polarization promoted by sparse higher-order inter- actions, Communications Physics (2025)
work page 2025
- [50]
-
[51]
N. W. Landry and J. G. Restrepo, The effect of het- erogeneity on hypergraph contagion models, Chaos: An Interdisciplinary Journal of Nonlinear Science30(2020)
work page 2020
-
[52]
G. St-Onge, I. Iacopini, V. Latora, A. Barrat, G. Petri, A. Allard, and L. Hébert-Dufresne, Influential groups for seeding and sustaining nonlinear contagion in het- erogeneous hypergraphs, Communications Physics5, 25 (2022)
work page 2022
- [53]
-
[54]
F. Malizia, A. Guzmán, I. Iacopini, and I. Z. Kiss, Disen- tangling the role of heterogeneity and hyperedge overlap in explosive contagion on higher-order networks, Physical Review Letters135, 207401 (2025)
work page 2025
- [55]
-
[56]
G. Marwell and R. E. Ames, Experiments on the pro- vision of public goods. i. resources, interest, group size, and the free-rider problem, American Journal of Sociol- ogy84, 1335 (1979)
work page 1979
-
[57]
A. Traulsen, M. A. Nowak, and J. M. Pacheco, Stochastic dynamics of invasion and fixation, Physical Review E74, 011909 (2006)
work page 2006
-
[58]
Q. F. Lotito, M. Contisciani, C. De Bacco, L. Di Gae- tano, L. Gallo, A. Montresor, F. Musciotto, N. Ruggeri, and F. Battiston, Hypergraphx: a library for higher- order network analysis, Journal of Complex Networks11, cnad019 (2023). Appendix A: Details of the stochastic simulations. Extraction of stationary branches from temporal trajectories The stochas...
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.