We present a rigorous mathematical framework for analyzing dynamics of a broad class of Boolean network models. We use this framework to provide the first formal proof of many of the standard critical transition results in Boolean network analysis, and offer analogous characterizations for novel classes of random Boolean networks. We precisely connect the short-run dynamic behavior of a Boolean network to the average influence of the transfer functions. We show that some of the assumptions traditionally made in the more common mean-field analysis of Boolean networks do not hold in general. For example, we offer some evidence that imbalance, or expected internal inhomogeneity, of transfer functions is a crucial feature that tends to drive quiescent behavior far more strongly than previously observed.Introduction. Complex systems can usually be represented as a network of interdependent functional units. Boolean networks were proposed by Kauffman as models of genetic regulatory networks [1, 2] and have received considerable attention across several scientific disciplines. They model a variety of complex phenomena, particularly in theoretical biology and physics [3][4][5][6][7][8].A Boolean network N with n nodes can be described by a directed graph G = (V, E) and a set of transfer functions. We use V and E to denote the sets of nodes and edges respectively, and denote the indegree of node i by K i . Each node i is assigned a K i -ary Boolean function f i : {−1, +1} K i → {−1, +1}, termed transfer function. If the state of node i at time t is x i (t), its state at time t + 1 is described by