Information transfer of an Ising model on a brain network

We consider an anatomical connectivity matrix (66 nodes), obtained via diffusion spectrum imaging (DSI) and white matter tractography [1], describing the brain at a coarse scale, and implement on it an Ising model with Glauber dynamics, estimating the transfer of information between spins. Tuning the temperature to criticality, so as to render the system characterized by long range correlations, we find that the critical state is characterized by the maximal amount of total information transfer among variables and exhibits signature of the law of diminishing marginal returns, a fundamental principle of economics which states that when the amount of a variable resource is increased, while other resources are kept fixed, the resulting change in the output will eventually diminish [2]. The origin of such behavior resides in the structural constraint related to the fact that each node of the network may handle a limited amount of information. The modulation of this phenomenon is analyzed evaluating, at each node, the ratio R between the outgoing and the incoming information. It turns out that R is related to (but not fully explained by) the in-strength of the brain network: nodes with high R are those more prone to became bottlenecks as the information transfer increases. The regions which are recognized as potential bottlenecks by the present analysis are symmetrical in the two hemisphere, and are Superior Frontal Cortex, Precuneus, Superior Temporal Cortex, Medial and Lateral Orbitofrontal Cortex. Some of these regions are considered as hubs both for the structural and for the functional connectome. Furthermore the temperature of the Ising model has an influences on the ratio between intra-hemispheric and inter-hemispheric information transfer, which grows with the temperature, though not uniformly.