Abstract
Macroscopic models of brain networks typically incorporate assumptions regarding the characteristics of afferent noise, which is used to represent input from distal brain regions or ongoing fluctuations in non-modelled parts of the brain. Such inputs are often modelled by Gaussian white noise which has a flat power spectrum. In contrast, macroscopic fluctuations in the brain typically follow a 1/f(b) spectrum. It is therefore important to understand the effect on brain dynamics of deviations from the assumption of white noise. In particular, we wish to understand the role that noise might play in eliciting aberrant rhythms in the epileptic brain. To address this question we study the response of a neural mass model to driving by stochastic, temporally correlated input. We characterise the model in terms of whether it generates "healthy" or "epileptiform" dynamics and observe which of these dynamics predominate under different choices of temporal correlation and amplitude of an Ornstein-Uhlenbeck process. We find that certain temporal correlations are prone to eliciting epileptiform dynamics, and that these correlations produce noise with maximal power in the δ and θ bands. Crucially, these are rhythms that are found to be enhanced prior to seizures in humans and animal models of epilepsy. In order to understand why these rhythms can generate epileptiform dynamics, we analyse the response of the model to sinusoidal driving and explain how the bifurcation structure of the model gives rise to these findings. Our results provide insight into how ongoing fluctuations in brain dynamics can facilitate the onset and propagation of epileptiform rhythms in brain networks. Furthermore, we highlight the need to combine large-scale models with noise of a variety of different types in order to understand brain (dys-)function.