Number of attractors in the critical Kauffman model is exponential

We derive mostly analytically the scaling behavior of the number of nonfrozen and relevant nodes in critical Kauffman networks (with two inputs per node) in the thermodynamic limit. By defining and analyzing a stochastic process that determines the frozen core we can prove that the mean number of nonfrozen nodes scales with the network size N as N^{2/3}, with only N^{1/3} nonfrozen nodes having two nonfrozen inputs. We also show the probability distributions for the numbers o...

We determine the average number $ \vartheta (N, K) $, of \textit{NK}-Kauffman networks that give rise to the same binary function. We show that, for $ N \gg 1 $, there exists a connectivity critical value $ K_c $ such that $ \vartheta(N,K) \approx e^{\phi N} $ ($ \phi > 0 $) for $ K < K_c $ and $\vartheta(N,K) \approx 1 $ for $ K > K_c $. We find that $ K_c $ is not a constant, but scales very slowly with $ N $, as $ K_c \approx \log_2 \log_2 (2N / \ln 2) $. The problem of ge...

This is the first of two papers about the structure of Kauffman networks. In this paper we define the relevant elements of random networks of automata, following previous work by Flyvbjerg and Flyvbjerg and Kjaer, and we study numerically their probability distribution in the chaotic phase and on the critical line of the model. A simple approximate argument predicts that their number scales as sqrt(N) on the critical line, while it is linear with N in the chaotic phase and in...

This is the second paper of a series of two about the structural properties that influence the asymptotic dynamics of Random Boolean Networks. Here we study the functionally independent clusters in which the relevant elements, introduced and studied in our first paper, are subdivided. We show that the phase transition in Random Boolean Networks can also be described as a percolation transition. The statistical properties of the clusters of relevant elements (that we call modu...

We derive analytically the scaling behavior in the thermodynamic limit of the number of nonfrozen and relevant nodes in the most general class of critical Kauffman networks for any number of inputs per node, and for any choice of the probability distribution for the Boolean functions. By defining and analyzing a stochastic process that determines the frozen core we can prove that the mean number of nonfrozen nodes in any critical network with more than one input per node scal...

Random Boolean networks, originally invented as models of genetic regulatory networks, are simple models for a broad class of complex systems that show rich dynamical structures. From a biological perspective, the most interesting networks lie at or near a critical point in parameter space that divides ``ordered'' from ``chaotic'' attractor dynamics. In the ordered regime, we show rigorously that the average number of relevant nodes (the ones that determine the attractor dyna...

A model of cellular metabolism due to S. Kauffman is analyzed. It consists of a network of Boolean gates randomly assembled according to a probability distribution. It is shown that the behavior of the network depends very critically on certain simple algebraic parameters of the distribution. In some cases, the analytic results support conclusions based on simulations of random Boolean networks, but in other cases, they do not.

We investigate the influence of a deterministic but non-synchronous update on Random Boolean Networks, with a focus on critical networks. Knowing that ``relevant components'' determine the number and length of attractors, we focus on such relevant components and calculate how the length and number of attractors on these components are modified by delays at one or more nodes. The main findings are that attractors decrease in number when there are more delays, and that periods ...

Gene regulatory network (GRN)-based morphogenetic models have recently gained an increasing attention. However, the relationship between microscopic properties of intracellular GRNs and macroscopic properties of morphogenetic systems has not been fully understood yet. Here we propose a theoretical morphogenetic model representing an aggregation of cells, and reveal the relationship between criticality of GRNs and morphogenetic pattern formation. In our model, the positions of...

Complexity theory as practiced by physicists and computational complexity theory as practiced by computer scientists both characterize how difficult it is to solve complex problems. Here it is shown that the parameters of a specific model can be adjusted so that the problem of finding its global energy minimum is extremely sensitive to small changes in the problem statement. This result has implications not only for studies of the physics of random systems but may also lead t...