Masuda N, Goh K-I, Kahng B
Laboratory for Mathematical Neuroscience, RIKEN Brain Science Institute, 2-1, Hirosawa, Wako, Saitama 351-0198, Japan.
Phys Rev E Stat Nonlin Soft Matter Phys. 2005 Dec;72(6 Pt 2):066106. doi: 10.1103/PhysRevE.72.066106. Epub 2005 Dec 6.
The Bak-Sneppen model displaying punctuated equilibria in biological evolution is studied on random complex networks. By using the rate equation and the random walk approaches, we obtain the analytic solution of the fitness threshold xc to be 1/((k)f+1), where (k)f=(k2)/(k) (=(k)) in the quenched (annealed) updating case, where kn is the nth moment of the degree distribution. Thus, the threshold is zero (finite) for the degree exponent gamma<3 (gamma>3) for the quenched case in the thermodynamic limit. The theoretical value xc fits well to the numerical simulation data in the annealed case only. Avalanche size, defined as the duration of successive mutations below the threshold, exhibits a critical behavior as its distribution follows a power law, Pa(s) approximately s(-3/2).