Dynamically critical systems and power-law distributions: Avalanches revisited

In this paper we show that a well-known model of genetic regulatory networks, namely that of Random Boolean Networks (RBNs), allows one to study in depth the relationship between two important properties of complex systems, i.e. dynamical criticality and power-law distributions. The study is based upon an analysis of the response of a RBN to permanent perturbations, that may lead to avalanches of changes in activation levels, whose statistical properties are determined by the same parameter that characterizes the dynamical state of the network (ordered, critical or disordered). Under suitable approximations, in the case of large sparse random networks an analytical expression for the probability density of avalanches of different sizes is proposed, and it is shown that for not-too-small avalanches of critical systems it may be approximated by a power law. In the case of small networks the above-mentioned formula does not maintain its validity, because of the phenomenon of self-interference of avalanches, which is also explored by numerical simulations.