field limit does not correspond to the limit of fast diffusion on networks if the interaction is *one-to-all* as is the case in our mafia model. Then, the dynamics of the network-structured society is actually better described within an extended mean-field theory which accounts for the fact that for each given node there is a finite set of neighborhood compositions, cf. Appendix A. In the following we discuss the standard mean-field approach in order to highlight the novel effects introduced by the network structure and the *one-to-all* interaction scheme.

We non-dimensionalize by choosing the inverse death rate, $d^{-1}$, as our basic time scale and define dimensionless time and all other parameters accordingly: $\tau = td$, $\beta = b/d$, $\sigma_i = s_i/d$, $\omega_{ij} = w_{ij}/d$, with $i,j \in \{c,m\}$. Then, the mean-field equations for the time evolution of the different population fractions read: $$\begin{aligned}
\dot{c}  & = & \phi\,\beta + m\,\omega_{mc}  - c\,\omega_{cm}  - c,\label{eq: dotc}\\
\dot{m} & = &  -m\,\omega_{mc}  + c\,\omega_{cm}  - m,\label{eq: dotm}\\
\dot{\phi} & = & -\phi\,\beta  +  (1-\phi ),\label{eq: dotphi}\end{aligned}$$ where a dot signifies a time derivative with respect to $\tau$. From Eq.  we immediately see that the stationary fraction of empty places is independent of the social dynamics encoded in the transition rates $\omega_{cm}$ and $\omega_{mc}$, but only depends on the dimensionless birth rate $\beta$: $$\phi=\frac{1}{1+\beta}\,.$$ In other words, the birth rates determines how densely populated the system is. Inserting this result together with the constraint $c+m+\phi=1$ into Eq.  yields an implicit equation, whose solutions are the fixed points or stationary states of the population dynamics: $$m\,\omega_{mc}(c,m,p) - \left(\frac{\beta}{1+\beta} - m\right)\,\omega_{cm}(c,m,p)  + m = 0.
\label{eq:master eq}$$


field limit does not correspond to the limit of firm dissemination on networks if the interaction is * one - to - all * as is the case in our mafia model. Then, the moral force of the network - structured society is actually well described within an extended average - field theory which account for the fact that for each given node there be a finite set of neighborhood compositions, cf. Appendix A. In the following we discuss the standard average - field approach in order to foreground the novel effects introduce by the net structure and the * one - to - all * interaction scheme. 

 We non - dimensionalize by choosing the inverse death rate, $ d^{-1}$, as our basic time plate and define dimensionless time and all other parameters accordingly: $ \tau = td$, $ \beta = b / d$, $ \sigma_i = s_i / d$, $ \omega_{ij } = w_{ij}/d$, with $ i, j \in \{c, m\}$. Then, the beggarly - plain equations for the time evolution of the different population fractions read: $ $ \begin{aligned } 
 \dot{c }   & = & \phi\,\beta + m\,\omega_{mc }   - c\,\omega_{cm }   - c,\label{eq: dotc}\\ 
 \dot{m } & = &   -m\,\omega_{mc }   + c\,\omega_{cm }   - m,\label{eq: dotm}\\ 
 \dot{\phi } & = & -\phi\,\beta   +   (1-\phi), \label{eq: dotphi}\end{aligned}$$ where a acid mean a time derivative with obedience to $ \tau$. From Eq.    we immediately visit that the stationary fraction of empty places is independent of the social moral force encoded in the passage rates $ \omega_{cm}$ and $ \omega_{mc}$, but only depend on the dimensionless birth rate $ \beta$: $ $ \phi=\frac{1}{1+\beta}\,.$$ In other words, the birth rates determine how densely populated the system is. Inserting this resultant role together with the restraint $ c+m+\phi=1 $ into Eq.    give an implicit equation, whose solutions are the fixed point or stationary states of the population dynamics: $ $ m\,\omega_{mc}(c, m, p) - \left(\frac{\beta}{1+\beta } - m\right)\,\omega_{cm}(c, m, p)   + m = 0. 
 \label{eq: master eq}$$
fiepd limit does not correskond to the limit of fasv diffuaion on vetworks if the interaction ms *obe-to-aol* as is the case in ojr mafia lodel. Thwn, tie dynamics of tis netwovh-strhgtureb wociety is actoally better described witviv cn extended mean-field theory which ascounts flr the fact thwt fpw eadh given node there is a finite sef of nepghborhood composotions, cf. Appendix A. In thf foplowing we discuss the standats mqqn-field appruach in order to highljght the novel effects introducdd by the netwotk shtucture and vhe *onv-to-all* interagnion scveme.

We mon-dimensionallze bb chiosing the inverse deeth rate, $d^{-1}$, as our bafic time vccle and define dimensuobless tima ana alu ouhec pzrametfrs accordingmy: $\tau = td$, $\veta = b/d$, $\sigma_i = s_i/c$, $\jnega_{ij} = w_{ij}/d$, sith $i,t \yn \{c,m\}$. Then, the mean-field equations for uhe tjme evolution of the didferent population frwctions rqad: $$\begin{aligned}
\dot{c}  & = & \phi\,\beta + m\,\omega_{mc}  - c\,\omeca_{cm}  - c,\lcncl{eq: eohc}\\
\dot{m} & = &  -m\,\omega_{mc}  + c\,\omega_{cm}  - m,\label{eq: doti}\\
\sou{\php} & = & -\phi\,\beta  +  (1-\pmi ),\label{eq: dotphi}\emd{wlodned}$$ where a aot siyhiries a time derivahive wijh respect to $\twu$. Ftom Eq.  we immediately see tyat the statponaey fraction of empcy places is indgpendemt of the social dynamies enckded in the transitikv rates $\omega_{cm}$ xnd $\okega_{mc}$, but only depends on the dimeisionkess bifth tate $\beea$: $$\phi=\frac{1}{1+\heta}\,.$$ Lt other words, the hirth rdtes deterlines how densely populated the system is. Insgrthng this refult bogether with tre constraint $e+m+\phi=1$ inco Eq.  hields an pmplicit xquation, whofe solutions dte the fixed 'oints or starionqry stagds of the popukation dyucmics: $$m\,\omwga_{mc}(c,m,p) - \left(\frac{\neta}{1+\bgtz} - m\right)\,\omega_{cn}(c,m,k)  + m = 0.
\label{eq:maxtef ez}$$


field limit does not correspond to the fast on networks the interaction is in mafia model. Then, dynamics of the society is actually better described within extended mean-field theory which accounts for the fact that for each given node is a finite set of neighborhood compositions, cf. Appendix A. In the following discuss standard approach order to highlight the novel effects introduced by the network structure and the *one-to-all* interaction scheme. non-dimensionalize by choosing the inverse death rate, $d^{-1}$, our basic time scale define dimensionless time and all parameters $\tau = $\beta b/d$, = s_i/d$, $\omega_{ij} w_{ij}/d$, with $i,j \in \{c,m\}$. Then, the mean-field equations for the time evolution of the different population read: $$\begin{aligned} = & + - - c,\label{eq: dotc}\\ = & -m\,\omega_{mc} + c\,\omega_{cm} - & = & -\phi\,\beta + (1-\phi ),\label{eq: dotphi}\end{aligned}$$ a dot a time derivative with respect to From Eq. we immediately see that the stationary of empty places is independent of the social dynamics encoded in the transition rates $\omega_{cm}$ but only depends on dimensionless birth rate $$\phi=\frac{1}{1+\beta}\,.$$ other the rates determines densely populated the system is. Inserting this result together with the $c+m+\phi=1$ into Eq. yields an implicit equation, whose solutions are points stationary states of population dynamics: $$m\,\omega_{mc}(c,m,p) - - + m = 0.
field limit does not correspoNd to the limIt of fAst DifFuSion On neTworks if the intERactIon is *one-to-all* as is the caSe in oUr MAfia MOdEl. TheN, the dynAMiCS Of tHe NeTwoRk-STrUcturEd sOciety iS actually bEttEr Described witHIn An extended MeaN-field theory WhiCh accoUnTs fOR the fAct That fOr each GIven noDe there is A fINite seT Of neighBORhOod cOmpositions, cf. AppeNDiX a. In the followinG we disCuSS tHE StaNdaRd mean-fielD aPproaCH in ordeR To HIGHliGHt the novel effEcts introduCEd bY the neTwOrk STructuRe and ThE *One-To-all* interaCtioN scheme.

We Non-dimENsionalIZe by choOsing tHe iNveRse dEAtH rAte, $D^{-1}$, aS Our BAsIc tIMe sCale and dEfInE dimeNsioNLESS timE anD all Other Parameters accOrdInglY: $\Tau = Td$, $\betA = b/d$, $\siGma_i = S_i/D$, $\omegA_{ij} = w_{ij}/D$, with $I,j \In \{c,m\}$. Then, the mean-FielD equationS foR tHe tImE evolUTion of The DifFerent pOpulatiON frAcTIONs Read: $$\begin{aligned}
\doT{c}  & = & \PHI\,\bEta + m\,\omegA_{mc}  - c\,\omEGa_{Cm}  - C,\Label{eq: dOtC}\\
\doT{m} & = &  -m\,\oMEGa_{mc}  + c\,\OmegA_{Cm}  - M,\label{eq: Dotm}\\
\doT{\PhI} & = & -\pHi\,\beta  +  (1-\pHi ),\Label{eQ: dOtpHi}\eNd{aliGNed}$$ wHere a dOt signifIes a tIMe derivative wiTH respect to $\tau$. fRoM eQ.  wE ImmeDiaTely see that The sTAtioNary FRaCtiON of emPty plAcES iS Independent of the socIaL dynamIcs enCoded in the traNsition ratES $\OMega_{cm}$ anD $\omeGA_{mC}$, But only depends On the DimensionlESs birth rAte $\beTa$: $$\phi=\fraC{1}{1+\beta}\,.$$ In otHER words, thE biRth RatEs dETErMines how denseLY PopuLaTed the sYstEm is. InsErtIng ThiS reSuLt togetheR with the CoNsTrAiNt $c+M+\phi=1$ iNTo Eq.  yielDs An iMpLicIt equATion, whOse soLutiOnS aRE thE fixed pOInTS Or stAtIoNary StaTeS of thE popULatIon dynaMics: $$m\,\omegA_{mc}(C,M,p) - \leFt(\FrAc{\beta}{1+\bEta} - m\right)\,\omegA_{cM}(c,m,p)  + m = 0.
\label{Eq:MasTer eq}$$


field limit does not corre spond to t he li mit  of  f astdiff usion on netwo r ks i f the interaction is * one-t o- a ll*a sis th e casei no u r m af ia  mo de l .Then,  th e dynam ics of the  ne tw ork-structur e dsociety is  ac tually bette r d escrib ed  wi t hin a n e xtend ed mea n -field  theory w hi c h acco u nts for t he  fac t that for each g i ve n  node there is  a fin it e  s e t  of  ne ighborhood  c ompos i tions,c f. A p pen d ix A. In thefollowing w e  di scussth e s t andard  mean -f i eld  approach i n or der to hi ghligh t  the no v el effe cts in tro duc ed b y  t he  ne tw o rks tr uct u reand the*o ne -to-a ll*i n t e ract ion  sch eme.
We non-dimen sio nali z e b y cho osing  the  i nvers e deat h rat e,  $d^{-1}$, as o ur b asic time  sc al e a nd  defi n e dime nsi onl ess tim e and a l l o th e r pa rameters according ly : $\ tau = td $, $\b e ta  = b/d$, $\ si gma _i = s _i/d$ , $\ o me ga_{ij}= w_{i j }/ d$ , with$i ,j \in  \ {c, m\} $. Th e n, t he mea n-fieldequat i ons for the ti m e evolution o f  t h e  d i ffer ent  population  fra c tion s re a d:  $$ \ begin {alig ne d }\ dot{c}  & = & \phi\ ,\ beta +  m\,\ omega_{mc}  -  c\,\omega _ { c m}  - c, \lab e l{ e q: dotc}\\
\do t{m}& = &  -m\ , \omega_{ mc} + c\,\om ega_{cm}-  m,\labe l{e q:dot m}\ \ \d ot{\phi} & =& -\ph i\ ,\beta + (1-\ph i ) ,\l abe l{e q:  dotphi}\ end{alig ne d} $$  w her e a d o t signif ie s a  t ime  deri v ativewithresp ec tt o $ \tau$.F ro m Eq. w eimme dia te ly se e th a t t he stat ionary fr act i on o fem pty pla ces is indepe nd ent of the  s oci al dyn a m ics enco ded in the transition r a tes $\o meg a_{cm }$ a nd $\omeg a_{ mc}$,but only d epends  on t he  di m e nsion l e ss  bi rt h rate $\b e t a$:  $$\p hi =\fr ac{1}{1 +\beta}\,.$$ In ot h erwords, the bi rth  rat e s  d ete r mi n esho w  de n s ely populated t he systemis .  I nserting t h isre sult to getherwitht he cons traint $c +m+\phi=1 $into E q. yields an  implici t equatio n , who s esolut ion s areth e f ixedpoints orstati onaryst ates o f the  p opulatio n dynamics: $$m\,\omega _{mc}( c,m,p ) -  \left(\f rac { \be ta}{1+\be ta}- m\right) \,\ ome ga_{c m}( c ,m,p)   +m  =  0. \labe l{eq : master eq } $$ 


field limit_does not_correspond to the limit_of fast_diffusion_on networks_if_the interaction is_*one-to-all* as is_the case in our_mafia model. Then,_the_dynamics of the network-structured society is actually better described within an extended mean-field theory_which_accounts for_the_fact_that for each given node_there is a finite set_of neighborhood_compositions, cf. Appendix A. In the following we_discuss_the standard mean-field_approach in order to highlight the novel effects introduced_by the network structure and the_*one-to-all* interaction scheme.

We_non-dimensionalize_by_choosing the inverse death_rate, $d^{-1}$, as our basic time_scale and define dimensionless time and_all other parameters accordingly: $\tau = td$,_$\beta = b/d$, $\sigma_i = s_i/d$,_$\omega_{ij} = w_{ij}/d$, with $i,j_\in \{c,m\}$._Then, the mean-field equations for_the time evolution_of the_different population fractions_read: $$\begin{aligned}
\dot{c}  & = &_\phi\,\beta + m\,\omega_{mc}_ - c\,\omega_{cm}  - c,\label{eq:_dotc}\\
\dot{m}_& = &__-m\,\omega_{mc}_ +_c\,\omega_{cm}  -_m,\label{eq:_dotm}\\
\dot{\phi} &_=_& -\phi\,\beta  +  (1-\phi_),\label{eq:_dotphi}\end{aligned}$$ where a dot signifies a time_derivative with respect to_$\tau$._From Eq.  we immediately_see that the stationary fraction_of empty places is independent of_the social_dynamics encoded_in the transition rates $\omega_{cm}$ and $\omega_{mc}$, but only depends on_the dimensionless birth rate $\beta$: $$\phi=\frac{1}{1+\beta}\,.$$_In other words, the_birth rates_determines_how densely populated_the_system is._Inserting this result together with the constraint_$c+m+\phi=1$ into_Eq.  yields an implicit equation, whose_solutions are the fixed_points_or stationary states of the population_dynamics: $$m\,\omega_{mc}(c,m,p) - \left(\frac{\beta}{1+\beta} - m\right)\,\omega_{cm}(c,m,p)_ + m = 0.
\label{eq:master_eq}$$

