The role of additive and diffusive coupling on the dynamics of neural populations

Dynamical models consisting of networks of neural masses commonly assume that the interactions between neural populations are via additive or diffusive coupling. When using the additive coupling, a population's activity is affected by the sum of the activities of neighbouring populations. In contrast, when using the diffusive coupling a neural population is affected by the sum of the differences between its activity and the activity of its neighbours. These two coupling functions have been used interchangeably for similar applications. In this study, we show that the choice of coupling can lead to strikingly different brain network dynamics. We focus on a phenomenological model of seizure transitions that has been used both with additive and diffusive coupling in the literature. We consider small networks with two and three nodes, as well as large random and scale-free networks with 64 nodes. We further assess resting-state functional networks inferred from magnetoencephalography (MEG) from people with juvenile myoclonic epilepsy (JME) and healthy controls. To characterize the seizure dynamics on these networks, we use the escape time, the brain network ictogenicity (BNI) and the node ictogenicity (NI), which are measures of the network's global and local ability to generate seizure activity. Our main result is that the level of ictogenicity of a network is strongly dependent on the coupling function. Overall, we show that networks with additive coupling have a higher propensity to generate seizures than those with diffusive coupling. We find that people with JME have higher additive BNI than controls, which is the hypothesized BNI deviation between groups, while the diffusive BNI provides opposite results. Moreover, we find that the nodes that are more likely to drive seizures in the additive coupling case are more likely to prevent seizures in the diffusive coupling case, and that these features correlate to the node's number of connections. Consequently, previous results in the literature involving such models to interrogate functional or structural brain networks could be highly dependent on the choice of coupling. Our results on the MEG functional networks and evidence from the literature suggest that the additive coupling may be a better modeling choice than the diffusive coupling, at least for BNI and NI studies. Thus, we highlight the need to motivate and validate the choice of coupling in future studies involving network models of brain activity.