Spontaneous cooperation for public goods

TL;DR

This paper demonstrates that spontaneous cooperation in public goods can emerge from simple local conformity dynamics influenced by noise, without the need for rational incentives or monitoring.

Contribution

It introduces an Ising model-based explanation for how random noise can trigger cascades of cooperation in public goods dilemmas.

Findings

Cooperation can spontaneously emerge at a critical noise level.

An analytic relationship links phase transition to asymmetry in cooperation and defection.

Random noise alone can resolve the cooperation dilemma without rational decision-making.

Abstract

Cooperation for public goods poses a dilemma, where individuals are tempted to free ride on others' contributions. Classic solutions involve monitoring, reputation maintenance and costly incentives, but there are important collective actions based on simple and cheap cues only, for example unplanned protests and revolts. This can be explained by an Ising model with the assumption that individuals in uncertain situations tend to conform to the local majority in their network. Among initial defectors, noise such as rumors or opponents' provocations causes some of them to cooperate accidentally. At a critical level of noise, these cooperators trigger a cascade of cooperation. We find an analytic relationship between the phase transition and the asymmetry of the Ising model, which in turn reflects the asymmetry of cooperation and defection. This study thereby shows that in principle, the dilemma of cooperation can be solved by nothing more than a portion of random noise, without rational decision making.

Equations

$${\Pi }_C^g$$

$${\Pi }_D^g$$

$$P_C^g$$

$$P_D^g$$

$$H=-\sum _{i\ne j}^n{a}_{ij}{S}_i{S}_j.$$

$$S=\{C,-D\}\to \{S_0+\Delta ,S_0-\Delta \}$$

$$H=-\sum _{ij}{a}_{ij}(S_0+{\hat{S}}_i)(S_0+{\hat{S}}_j).$$

$$-\sum _{ij}{a}_{ij}{\hat{S}}_i{\hat{S}}_j-S_0\left(\sum _{ij}{a}_{ij}{\hat{S}}_i+\sum _{ij}{a}_{ij}{\hat{S}}_j\right)-S_0^2\sum _{ij}{a}_{ij}.$$

$$H_{loc}=-S_0\sum _i\left(\sum _j{a}_{ij}+\sum _j{a}_{ji}\right){\hat{S}}_i.$$

$$H_{loc}^{row}=-S_0\sum _i\left(\sum _j{a}_{ij}\right){\hat{S}}_i-S_0\sum _i{\hat{S}}_i,$$

$$H_{loc}^{row}=-2S_0\sum _i{\hat{S}}_i.$$

$$M=\frac{(S_0-\Delta )e^{\frac{(S_0-\Delta )M}{T}}+(S_0+\Delta )e^{\frac{(S_0+\Delta )M}{T}}}{e^{\frac{(S_0-\Delta )M}{T}}+e^{\frac{(S_0+\Delta )M}{T}}}.$$

$$\frac{M}{\Delta }=\frac{S_0}{\Delta }+\tanh \left(\frac{{\Delta }^2}{T}\frac{M}{\Delta }\right),$$

$$\frac{S_0}{\Delta }=\frac{\sqrt{\frac{\gamma -1}{\gamma +1}}(\gamma +1)\gamma -{\cosh }^{-1}(\gamma )}{{\gamma }^2},$$

$$M_c=p_c(S_0-\Delta )+(1-p_c)(S_0+\Delta ).$$

$$\frac{M_c}{\Delta }=\frac{{-\cosh }^{-1}(\gamma )}{{\gamma }^2}.$$

$$p_c=\frac{1}{2}-\frac{S_0}{2\Delta }+\frac{1}{2}\frac{{\cosh }^{-1}(\gamma )}{{\gamma }^2},$$

Full text

People may want to realize or preserve public goods, for example democracy and clean air, but because contributors are disadvantaged in the face of free riders, there is a dilemma [\citeauthoryearOlsonOlson1965, \citeauthoryearHardinHardin1968, \citeauthoryearOstromOstrom2009, \citeauthoryearGavriletsGavrilets2015]. Solutions typically require efforts of the participants to monitor one another [\citeauthoryearRustagi, Engel, and KosfeldRustagi et al.2010] and spread information (gossip) [\citeauthoryearNowak and SigmundNowak and Sigmund2005] through their network reliably that establishes reputations [\citeauthoryearPanchanathan and BoydPanchanathan and Boyd2004], upon which some of them have to deliver individual rewards or (threats of) punishments [\citeauthoryearFehr and FischbacherFehr and Fischbacher2003], under pro-social norms to preclude arbitrariness. These provisions are not always (sufficiently) available, though, whereas in certain situations, participants still manage to self-organize cooperation, even without leaders. Cases in point are impromptu help at disasters, non-organized revolts against political regimes [\citeauthoryearLohmannLohmann1994, \citeauthoryearTillyTilly2002, \citeauthoryearTufekciTufekci2017] and spontaneous street fights between groups of young men.

These examples have in common a high uncertainty of outcomes, and unknown benefits and costs. Rational decision making is therefore not feasible. Participants who identify with their group or its goal [\citeauthoryearVan Stekelenburg and KlandermansVan Stekelenburg and Klandermans2013], and thereby feel group solidarity [\citeauthoryearDurkheimDurkheim1912], use the heuristic of conformism to the majority of their network neighbors \shortcitewu14, which can be based on no more than visual information. Human ancestors lived in groups for millions of years [\citeauthoryearShultz, Opie, and AtkinsonShultz et al.2011] and in all likelihood, solidarity and conformism are both cultural and genetic [\citeauthoryearBoydBoyd2018]. On an evolutionary time scale, conformism must have been beneficial on average when future benefits and costs were unknown [\citeauthoryearVan den Berg and WenseleersVan den Berg and Wenseleers2018]. To explain cooperation under conformism, we use an Ising model \shortciteweidlich71,galam82,stauffer08,dorogovtsev08,castellano09. This model recovers the critical mass that makes cooperation self-reinforcing but without the rationality assumptions of critical mass theory [\citeauthoryearMarwell and OliverMarwell and Oliver1993].

Model

A group $g$ of individuals who share an interest in a public good is modeled as a network with weighted and usually asymmetric ties $A_{ij}$ denoting $i$ paying attention to $j$. There is no assumption that $i$ and $j$ know each other before they meet at the site where collective action might take place. Consistent with other models of social influence [\citeauthoryearFriedkin and JohnsenFriedkin and Johnsen2011], the adjacency matrix is row-normalized, yielding cell values $a_{ij}=A_{ij}/\sum_{j}A_{ij}$, hence $\sum_{j}a_{ij}$ = 1.

Individuals have two behavioral options, defect ($D$) and cooperate ($C$), $C>D>0$, and all defect at the start. The average degree of cooperation among $n$ individuals is described by an order parameter $M=1/n\sum_{i=1}^{n}S_{i}$, where the behavioral variable $S_{i}$ can take the value $S_{i}=C$ or $S_{i}=-D$, for example $S=\{1,-1/2\}$. Everybody gets an equal share of the public good but cooperators incur a cost. A widely used definition of payoffs for cooperators $\Pi^{g}_{C}$ and defectors $\Pi^{g}_{D}$ in a group $g$ is the following \shortciteperc17,

[TABLE]

with $r>1$ an enhancement, or synergy, factor of cooperation, and $N_{C}$ the number of cooperators when the focal player decides. In our case,

[TABLE]

which is identical to Eq. 1 in standard game theory when $C=1$ and $D=0$. Our payoffs can be negative but that does not matter because they are used only comparatively. Other, for example non-linear, payoff functions [\citeauthoryearMarwell and OliverMarwell and Oliver1993] may also be used. The key point, however, is that under high uncertainty, participants do not maximize their payoff (directly) but align with others instead, thereby forming a collective lever that increases their payoff while avoiding exploitation. Behavior and network ties are expressed in the conventional, but here asymmetric, Ising model

[TABLE]

Solving the model boils down to minimizing $H$, where $H/n$ can be interpreted as average dissatisfaction. Minimizing can be done computationally with a Metropolis algorithm [\citeauthoryearBarrat, Barthelemy, and VespignaniBarrat et al.2008], where individuals decide sequentially as in many network models, or through a mean field analysis, elaborated below.

High-uncertainty situations to which the model applies are characterized by turmoil, $T$, or temperature in the original model. It causes arousal, measurable as heart rates \shortcitekonvalinka11, and produces noise [\citeauthoryearLewenstein, Nowak, and LatanéLewenstein et al.1992] in individuals’ information about the situation, which in turn becomes partly false, ambiguous, exaggerated or objectively irrelevant. Turmoil and its noise may consist of rumors, fire, provocations and violence. Some social movements produce turmoil by themselves, for instance an increasingly frequent posting of online messages \shortcitejohnson16. Arousal and noise entail “trembling hands” [\citeauthoryearDion and AxelrodDion and Axelrod1988] as game theorists say, which means a chance that some individuals accidentally change their behavior. The model is to show that few accidental cooperators entail a cascade of cooperation. An example of turmoil and its ramifications is the self-immolation of a street vendor in December 2010, which, in the given circumstances, set off the Tunisian revolution. Other examples are the revolts in East Germany [\citeauthoryearLohmannLohmann1994] and Romania in 1989 and in Egypt and Syria in 2011 [\citeauthoryearHussain and HowardHussain and Howard2013], where protesters were agitated by rumors about the events in neighboring countries. Autocratic rulers try to prevent revolts by suppressing turmoil, for instance by tightening media control.

Noise is different from a stable bias, for instance the ideology of an autocratic regime, which entails revolts against it less often than a weakened regime, or stumbling opponents in street fights. Opponents’ weakness gives off noisy signals that they might be overcome, which readily entail collective actions against them [\citeauthoryearSkocpolSkocpol1979, \citeauthoryearGoldstoneGoldstone2001, \citeauthoryearCollinsCollins2008]. Whereas responses to noise are typically spontaneous, collective responses to stable signals, biased or not, tend to be mounted by organized groups with norms, incentives and all that [\citeauthoryearTillyTilly2002, \citeauthoryearTufekciTufekci2017, \citeauthoryearGoldstoneGoldstone2001]. Combinations of signal and noise also occur, of course, which can result in, for example, an organized peaceful demonstration to suddenly turn violent. Our focus is on spontaneous cooperation.

Results

Our general result is that within finite time and at low turmoil, cooperation does not get off the ground, but it does emerge at a critical level $T_{c}$. This pattern is illustrated by the red line in Fig. 1(b) along the direction of the arrows. The figure was obtained with a mean field approach, but numerical simulations with the Metropolis algorithm show up the same pattern, with lower $T_{c}$ for small networks (Fig 3a). The model thus shows that at a critical level of turmoil, few accidental cooperators can trigger a cascade of cooperation. If $T$ keeps increasing way beyond $T_{c}$, cooperators co-exist with increasing numbers of defectors, until the two behaviors become equally frequent. If in actuality cooperation then collapses completely is an issue for further study. Otherwise cooperation ends when the public good is achieved, the participants run out of steam, or others intervene.

Alternatively, if participants get to understand an enduring situation, their uncertainty will reduce and they may start acting strategically, which requires pro-social norms to prevent. If the participants then develop such norms prescribing rewards and punishments, these norms can be easily modeled as field(s) by adding term(s) $-h\sum S_{i}$ to the Hamiltonian (Eq.3). Consequently, cooperation emerges without phase transition. The actual maintenance of these norms, however, will entail additional costs over and above the contributions, whereas spontaneous cooperation is relatively cheap.

Comparison with the symmetric Ising model

Rewriting the asymmetric Ising model in a symmetric form enables a direct comparison with results in the literature for symmetric models, and a generalization to arbitrary values of $S=\{C,-D\}$. A model with asymmetric values can be reformulated as a symmetric model with an offset, or bias, $S_{0}=(C-D)/2$ and an increment $\Delta=(C+D)/2$ by the mapping

[TABLE]

Accordingly, the values chosen in Fig. 1(b), $S=\{1,-1/2\}$, imply $\Delta=0.75$ and $S_{0}=0.25$. For given $\Delta$, increasing $S_{0}$ means increasing interest in the public good, or, in line with the literature on protests, increasing grievances [\citeauthoryearVan Stekelenburg and KlandermansVan Stekelenburg and Klandermans2013].

Substitution of $S_{0}$ and $\hat{S}_{i}$ chosen from $\{\Delta,-\Delta\}$ in $H$ yields

[TABLE]

Expanding $H$ in orders of $S_{0}$ yields

[TABLE]

The first term in the expansion $H_{sym}=-\sum_{ij}a_{ij}\hat{S}_{i}\hat{S}_{j}$ is a symmetric model with the same adjacency matrix as the original asymmetric model. The second term $H_{loc}=-S_{0}(\sum_{ij}a_{ij}\hat{S}_{i}+\sum_{ij}a_{ij}\hat{S}_{j})$ is proportional to $S_{0}$ and can be interpreted as a local field that modifies $H_{sym}$. The contribution of this local field can be expressed in terms of row and column sums of $a_{ij}$ as

[TABLE]

For row-normalized adjacency matrices, with $\sum_{j}a_{ij}=1$ for all rows $i$, $H_{loc}$ becomes

[TABLE]

where the first term is a local field varying for each $\hat{S}_{i}$, and the second term is a homogeneous external field independent of $a_{ij}$. The third term in the expansion of $H$ is independent of the values of $\hat{S}$ and is a constant depending on $a_{ij}$ only. Hence it does not play a role in the minimization of $H$. For a connected network with row-normalization, the last expression can be further simplified to

[TABLE]

The asymmetry in $S$ is then equivalent to a symmetric system with an external field $2S_{0}$.

Mean field analysis

The expected value of $M$ as function of $T$ can be obtained by assuming that the network is very large and by abstracting away from its topology; in the language of thermodynamics, by approximating the interaction energy by the energy of one spin (here, behavior) in the mean field of its neighbors [\citeauthoryearBarrat, Barthelemy, and VespignaniBarrat et al.2008], $M=\langle S\rangle$. The value of $M$ can now be expressed in closed form in terms of the probabilities given by the exponential of the Hamiltonian energy and $T$ as

[TABLE]

This reduces to an implicit relation

[TABLE]

where only dimensionless ratios of $M$, $S_{0}$ and $T$ with $\Delta$ remain in the expression. The mean degree $\langle k\rangle$, defined for binary ties, does not occur in it because the adjacency matrix is row-normalized and the mean weighted outdegree $\langle k_{w}\rangle=1$.

By analyzing the intersection of the line defined by $M/\Delta-S_{0}/\Delta$ and the tanh term on the right hand side of Eq. 11, the possible values for $M$ at a given $T$ can be found. For $T>T_{c}$ there is one stable high $T$ solution and for $T<T_{c}$ there is one stable solution of (nearly) full cooperation, another solution that is stable in finite time with (nearly) full defection, and one unstable solution. At $T=T_{c}$ the two stable solutions merge and the intersecting line coincides with the tangent line touching the tanh function; see Fig. 1(b). At that point a closed relation for $T_{c}$ in terms of $S_{0}$ and $\Delta$ can be found,

[TABLE]

where $\gamma=\sqrt{\frac{\Delta^{2}}{T_{c}}}.$ Eq. 12 is used in Fig. 3(a). It shows that if $S_{0}$ increases while keeping $\Delta$ constant, less agitation is required to motivate defectors to cooperate. When $\Delta$ decreases to $\Delta=S_{0}$, defection loses its appeal. The figure also shows that numerical simulations yield very similar results for large networks but diverge for small ones. This also holds true for $T_{c}$ in Fig 1, which is lower for smaller networks (not shown).

From the mean field approximation follows the proportion of defectors $p_{c}$ and cooperators $1-p_{c}$ at given $S_{0}$ and pertaining $T_{c}$, after a time long enough for the system to settle down. The proportion of cooperators $1-p_{c}$ at $T_{c}$ is the critical mass [\citeauthoryearMarwell and OliverMarwell and Oliver1993], and can be inferred from the value of $M_{c}$ at the phase transition,

[TABLE]

The mean field analysis of $T_{c}$ yields

[TABLE]

Note that the $\rm{cosh}^{-1}$ function in Eq. 14 only yields a result when $\gamma>1$, and sets a limit to $T_{c}$ for given $\Delta$. For the choice $\Delta=0.75$, the maximum value of $T_{c}=\Delta^{2}=0.565$. Solving for $p_{c}$ yields

[TABLE]

used for Fig. 3(b). It shows that the proportion of defectors $p_{c}$ at $T_{c}$ decreases with increasing $S_{0}$. In contrast to critical mass theory, however, the Ising model has no assumptions about initiative takers or leaders who win over the rest, rational decision making [\citeauthoryearMarwell and OliverMarwell and Oliver1993], or learning that would require fairly stable feedback [\citeauthoryearMacyMacy1991].

Discussion and Conclusion

We have shown that under high uncertainty, the dilemma of collective action can be solved by nothing more than a portion of random noise. The asymmetric Ising model does not require any knowledge or accurate expectations of the participants, and only depends on conformism, which can be empirically observed in synchronous motion, gestures or shouting [\citeauthoryearMcNeillMcNeill1995, \citeauthoryearJonesJones2013]. In particular, it has no assumptions about actors’ rationality, in contrast to critical mass theory, whereas it supports that theory’s key findings of the critical mass and the tipping point. Simulations add to our mean field result that turmoil-driven cooperation is most likely in small groups, where cooperation starts at relatively low levels of noise.

Shortly before we finished this manuscript, two other papers appeared where an Ising model was used to solve this dilemma [\citeauthoryearAdami and HintzeAdami and Hintze2018, \citeauthoryearSarkar and BenjaminSarkar and Benjamin2019], but their symmetric model requires complex quantum physics to define payoffs, in contrast to our simple definition.

In another paper, the effects of different networks and of multiple $S_{0}$ values on the tipping point are investigated [\citeauthoryearBruggeman and SprikBruggeman and Sprik2020]. Along with more empirical testing, perhaps also on other species, a future direction might be to explicitly model noisy information transmission on the group’s network [\citeauthoryearQuax, Apolloni, and SlootQuax et al.2013], too.

Acknowledgements

Thank you to José A. Cuesta, Raheel Dhattiwala, Don Weenink, and Alex van Venrooij for comments.

Author contributions

JB wrote the paper and made the asymmetric Ising model. RS related the asymmetric model to the symmetric one. RQ did the mean field analysis.