“…These models were solvable by virtue of the specific nature of their architectures: one either chooses strictly symmetric dilution (so detailed balance and hence equilibrium analysis are preserved, e.g. [9,10,11]), or strictly asymmetric dilution, which ensures that neuron states are statistically independent on finite times [8] (now the local fields are described by Gaussian distributions, leading to simple dynamic order parameter equations). In the early models, the bond statistics were uniform over the entire network, leading to thin tails in its degree distribution, whereas the connectivity of a real neuron is known to vary strongly within the brain [12].…”