 One is that its survival neighbors will move to $i$ with probability $\sum_j \frac{(1-p_{\rm d}) a_{ij} n_j(t)}{k_j}$ and the other is that its survival neighbor $j$ selecting an occupied neighbor $l$ (with the probability $\frac{1}{k_j}\sum_l a_{jl}
n_l(t))$) reproduces a new individual on $i$, then we have the term $\sum_j \frac{(1-p_{\rm d}) a_{ij} n_j(t)}{k_j} \frac{p_{\rm
b}}{k_j}\sum_l a_{jl} n_l(t))$. Taking the average of Eq. (\[evolution\]), we obtain

$$\begin{aligned}
\langle n_i(t+1) | {\bf n}(t)\rangle &=& n_i(t)(1-p_{\rm
d})\frac{1}{k_i}\sum_j a_{ij} n_j(t)+(1-n_i(t))\left[\sum_j
\frac{(1-p_{\rm d}) a_{ij} n_j(t)}{k_j} (1+\frac{p_{\rm
b}}{k_j}\sum_l a_{jl} n_l(t))\right], \label{evolution_averaged}\end{aligned}$$

which describes the average evolution of the system, conditioned to the knowledge of its state at the previous time step. In the MF approximation, $\langle n_i(t)n_j(t)\rangle \equiv \langle n_i(t)
\rangle \langle n_j(t) \rangle$ and $\langle
n_i(t)n_j(t)n_l(t)\rangle \equiv \langle n_i(t) \rangle \langle
n_j(t) \rangle \langle n_l(t) \rangle$. Thus, after multiplying Eq. (\[evolution\_averaged\]) by the probability to find the system at state ${\bf n}$ at time $t$ and summing for all possible configurations, we obtain

$$\begin{aligned}
\rho_i(t+1) & \equiv & \langle n_i(t+1) \rangle =
\rho
One is that its survival neighbors will move to $ i$ with probability $ \sum_j \frac{(1 - p_{\rm d }) a_{ij } n_j(t)}{k_j}$ and the other is that its survival neighbor $ j$ selecting an occupy neighbor $ l$ (with the probability $ \frac{1}{k_j}\sum_l a_{jl } 
 n_l(t))$) reproduce a new individual on $ i$, then we have the terminus $ \sum_j \frac{(1 - p_{\rm d }) a_{ij } n_j(t)}{k_j } \frac{p_{\rm 
 b}}{k_j}\sum_l a_{jl } n_l(t))$. Taking the average of Eq.   (\[evolution\ ]), we obtain 

 $ $ \begin{aligned } 
 \langle n_i(t+1) | { \bf n}(t)\rangle & = & n_i(t)(1 - p_{\rm 
 d})\frac{1}{k_i}\sum_j a_{ij } n_j(t)+(1 - n_i(t))\left[\sum_j 
 \frac{(1 - p_{\rm d }) a_{ij } n_j(t)}{k_j } (1+\frac{p_{\rm 
 b}}{k_j}\sum_l a_{jl } n_l(t))\right ], \label{evolution_averaged}\end{aligned}$$ 

 which identify the median evolution of the arrangement, conditioned to the knowledge of its department of state at the previous time dance step. In the MF approximation, $ \langle n_i(t)n_j(t)\rangle \equiv \langle n_i(t) 
 \rangle \langle n_j(t) \rangle$ and $ \langle 
 n_i(t)n_j(t)n_l(t)\rangle \equiv \langle n_i(t) \rangle \langle 
 n_j(t) \rangle \langle n_l(t) \rangle$. Thus, after multiplying Eq.   (\[evolution\_averaged\ ]) by the probability to find the system at state of matter $ { \bf n}$ at time $ t$ and summing for all possible shape, we obtain 

 $ $ \begin{aligned } 
 \rho_i(t+1) & \equiv & \langle n_i(t+1) \rangle = 
 \rho
 Onf is that its survival ntighbors will movg ro $i$ wmth progability $\sum_j \frac{(1-p_{\rm d}) a_{ij} n_j(t)}{k_j}$ anv thw othtg is that its survivau neighbog $j$ selecring qn occupiev neighbor $l$ (wjbh thz 'robability $\frag{1}{k_j}\sum_l a_{jl}
t_l(t))$) reproduces a nzw individual on $i$, then we have the eerm $\suk_j \frac{(1-p_{\rm d}) a_{ij} n_j(t)}{l_t} \frzb{p_{\vm
b}}{k_j}\sum_l a_{jl} n_l(t))$. Taking the avedage of Eq. (\[evolution\]), ee obtain

$$\begin{aligned}
\langlf n_i(h+1) | {\bf n}(t)\rangle &=& n_i(h)(1-p_{\rm
d})\frac{1}{k_i}\wum_j q_{ij} n_j(t)+(1-n_i(t))\lewt[\sum_j
\frac{(1-k_{\rk d}) a_{ij} n_j(j)}{k_j} (1+\frac{p_{\rm
b}}{k_j}\sum_l a_{jl} n_l(t))\right], \uabel{zvolution_avgxqgef}\gnd{aligned}$$

whmch defcribes the average avolutipn of the systcm, coidituoned to the knowledgx of its state at thg previous tnme step. In the MF appriximajion, $\nanguw n_k(t)n_n(t)\canfle \eqkiv \langle n_i(f)
\rangle \lantle n_j(t) \rangle$ and $\kagtle
n_i(t)n_j(t)n_l(t)\rzngle \qqtiv \langle n_i(t) \rangle \langle
n_j(t) \rangle \nanfle n_l(t) \rangle$. Thus, aftwr multiplying Eq. (\[evolotion\_averaded\]) by the probability to find the system at stata ${\bf i}$ xt uine $t$ xbd summing for all possible configurations, we offaon

$$\negin{aligned}
\rho_i(b+1) & \equiv & \langle n_o(t+1) \rsggle =
\rho
One is that its survival neighbors will $i$ probability $\sum_j d}) a_{ij} n_j(t)}{k_j}$ its neighbor $j$ selecting occupied neighbor $l$ the probability $\frac{1}{k_j}\sum_l a_{jl} n_l(t))$) reproduces new individual on $i$, then we have the term $\sum_j \frac{(1-p_{\rm d}) a_{ij} \frac{p_{\rm b}}{k_j}\sum_l a_{jl} n_l(t))$. Taking the average of Eq. (\[evolution\]), we obtain $$\begin{aligned} n_i(t+1) {\bf &=& d})\frac{1}{k_i}\sum_j a_{ij} n_j(t)+(1-n_i(t))\left[\sum_j \frac{(1-p_{\rm d}) a_{ij} n_j(t)}{k_j} (1+\frac{p_{\rm b}}{k_j}\sum_l a_{jl} n_l(t))\right], \label{evolution_averaged}\end{aligned}$$ which describes the average of the system, conditioned to the knowledge of state at the previous step. In the MF approximation, n_i(t)n_j(t)\rangle \langle n_i(t) \langle \rangle$ $\langle n_i(t)n_j(t)n_l(t)\rangle \equiv n_i(t) \rangle \langle n_j(t) \rangle \langle n_l(t) \rangle$. Thus, after multiplying Eq. (\[evolution\_averaged\]) by the probability to the system ${\bf n}$ time and for all possible obtain $$\begin{aligned} \rho_i(t+1) & \equiv & = \rho
 One is that its survival neighBors will moVe to $i$ WitH prObAbilIty $\sUm_j \frac{(1-p_{\rm d}) a_{ij} N_J(t)}{k_j}$ And the other is that its suRvivaL nEIghbOR $j$ SelecTing an oCCuPIEd nEiGhBor $L$ (wITh The prObaBility $\fRac{1}{k_j}\sum_l a_{Jl}
n_L(t))$) Reproduces a nEW iNdividual oN $i$, tHen we have the TerM $\sum_j \fRaC{(1-p_{\rM D}) a_{ij} n_J(t)}{k_J} \frac{P_{\rm
b}}{k_j}\SUm_l a_{jl} N_l(t))$. Taking ThE AveragE Of Eq. (\[evoLUTiOn\]), we Obtain

$$\begin{aligneD}
\LaNGle n_i(t+1) | {\bf n}(t)\rangLe &=& n_i(t)(1-p_{\Rm
D})\FrAC{1}{K_i}\sUm_j A_{ij} n_j(t)+(1-n_i(t))\lEfT[\sum_j
\FRac{(1-p_{\rm d}) A_{Ij} N_J(T)}{K_j} (1+\fRAc{p_{\rm
b}}{k_j}\sum_l a_{Jl} n_l(t))\right], \lABel{EvolutIoN_avERaged}\eNd{aliGnED}$$

whIch describeS the Average evOlutioN Of the sySTem, condItioneD to The KnowLEdGe Of iTs STatE At The PRevIous time StEp. in the mF apPROXImatIon, $\LangLe n_i(t)N_j(t)\rangle \equiV \laNgle N_I(t)
\rAngle \LanglE n_j(t) \RaNgle$ aNd $\langLe
n_i(t)N_j(T)n_l(t)\rangle \equiv \LangLe n_i(t) \rangLe \lAnGle
N_j(T) \rangLE \langlE n_l(T) \raNgle$. ThuS, after mULtiPlYING EQ. (\[evolution\_averaged\]) By THE pRobabiliTy to fiND tHe SYstem at sTaTe ${\bF n}$ at TIMe $t$ anD sumMInG for all pOssiblE CoNfIguratiOnS, we obtAiN

$$\beGin{AlignED}
\rho_I(t+1) & \equiV & \langle n_I(t+1) \ranGLe =
\rho
 One is that its survivalneighborswillmov e t o$i$with  probability $ \ sum_ j \frac{(1-p_{\rm d})a_{ij }n _j(t ) }{ k_j}$  and th e  o t h eris  t hat  i t ssurvi val  neighb or $j$ sel ect in g an occupie d  n eighbor $l $ ( with the pro bab ility$\ fra c {1}{k _j} \sum_ l a_{j l }
n_l( t))$) rep ro d uces a new ind i v id ualon $i$, then we h a ve the term $\sum _j \fr ac { (1 - p _{\ rmd}) a_{ij}  n _j(t) } {k_j} \ f ra c { p _{\ r m
b}}{k_j}\su m_l a_{jl}n _l( t))$.Ta kin g  the a verag eo f E q. (\[evolu tion \]), we o btain$$\begi n {aligne d}
\la ngl e n _i(t + 1)  |  {\ bf n}( t )\ ran g le&=& n_i( t) (1 -p_{\ rm
d } ) \ f rac{ 1}{ k_i} \sum_ j a_{ij} n_j( t)+ (1-n _ i(t ))\le ft[\s um_j 
\ frac{ (1-p_{ \rm d })  a_{ij} n_j(t)} {k_j } (1+\fra c{p _{ \rm 
b }}{k_ j }\sum_ l a _{j l} n_l( t))\rig h t],  \ l a b el {evolution_average d} \ e nd {aligned }$$

w h ic hd escribes  t heaver a g e evo luti o nof the s ystem, co nd itioned  t o thekn owl edg e ofi ts s tate a t the pr eviou s  time step. In the MF approx i ma t i on ,  $\l ang le n_i(t)n_ j(t) \ rang le \ e qu iv\ langl e n_i (t ) 
\ r angle \langle n_j(t )\rangl e$ an d $\langle
n_ i(t)n_j(t) n _ l (t)\rang le \ e qu i v \langle n_i( t) \r angle \lan g le
n_j(t ) \ra ngle \la ngle n_l( t )  \rangle $.Thu s,aft e r  m ultiplying Eq . (\[e vo lution\ _av eraged\ ])bythe  pr ob ability t o find t he  s ys te m a t sta t e ${\bfn} $ a ttim e $t$ and su mming  for  a ll pos sible c o nf i g urat io ns , we  ob ta in

$ $\be g in{ aligned }
\rho_i( t+1 )  & \ eq ui v & \la ngle n_i(t+1)  \ rangle =
\ rh o
 One_is that_its survival neighbors will_move to_$i$_with probability_$\sum_j_\frac{(1-p_{\rm d}) a_{ij}_n_j(t)}{k_j}$ and the_other is that its_survival neighbor $j$_selecting_an occupied neighbor $l$ (with the probability $\frac{1}{k_j}\sum_l a_{jl}
n_l(t))$) reproduces a new individual on_$i$,_then we_have_the_term $\sum_j \frac{(1-p_{\rm d}) a_{ij}_n_j(t)}{k_j} \frac{p_{\rm
b}}{k_j}\sum_l a_{jl} n_l(t))$. Taking_the average_of Eq. (\[evolution\]), we obtain

$$\begin{aligned}
\langle n_i(t+1) | {\bf n}(t)\rangle_&=&_n_i(t)(1-p_{\rm
d})\frac{1}{k_i}\sum_j a_{ij} n_j(t)+(1-n_i(t))\left[\sum_j
\frac{(1-p_{\rm_d}) a_{ij} n_j(t)}{k_j} (1+\frac{p_{\rm
b}}{k_j}\sum_l a_{jl} n_l(t))\right], \label{evolution_averaged}\end{aligned}$$

which describes the_average evolution of the system, conditioned_to the knowledge_of_its_state at the previous_time step. In the MF approximation,_$\langle n_i(t)n_j(t)\rangle \equiv \langle n_i(t)
\rangle \langle_n_j(t) \rangle$ and $\langle
n_i(t)n_j(t)n_l(t)\rangle \equiv \langle n_i(t)_\rangle \langle
n_j(t) \rangle \langle n_l(t) \rangle$._Thus, after multiplying Eq. (\[evolution\_averaged\]) by_the probability_to find the system at_state ${\bf n}$_at time_$t$ and summing_for all possible configurations, we obtain

$$\begin{aligned}
\rho_i(t+1)_& \equiv &_\langle n_i(t+1) \rangle =
\rho