Slow-Fast Duffing Neural Mass Model

Abstract

Epileptic seizures may be initiated by random neuronal fluctuations and/or by pathological slow regulatory dynamics of ion currents. This paper presents extensions to the Jansen and Rit neural mass model (JRNMM) to replicate paroxysmal transitions in intracranial electroencephalogram (iEEG) recordings. First, the Duffing NMM (DNMM) is introduced to emulate stochastic generators of seizures. The DNMM is constructed by applying perturbations to linear models of synaptic transmission in each neural population of the JRNMM. Then, the slow-fast DNMM is introduced by considering slow dynamics (relative to membrane potential and firing rate) of some internal parameters of the DNMM to replicate pathological evolution of ion currents. Through simulation, it is illustrated that the slow-fast DNMM exhibits transitions to and from seizures with etiologies that are linked either to random input fluctuations or pathological evolution of slow states. Estimation and optimization of a log likelihood function (LLF) using a continuous-discrete unscented Kalman filter (CD-UKF) and a genetic algorithm (GA) are performed to capture dynamics of iEEG data with paroxysmal transitions.