Population dynamics of the modified theta model: macroscopic phase reduction and bifurcation analysis link microscopic neuronal interactions to macroscopic gamma oscillation

Abstract

Gamma oscillations of the local field potential are organized by collective dynamics of numerous neurons and have many functional roles in cognition and/or attention. To mathematically and physiologically analyse relationships between individual inhibitory neurons and macroscopic oscillations, we derive a modification of the theta model, which possesses voltage-dependent dynamics with appropriate synaptic interactions. Bifurcation analysis of the corresponding Fokker-Planck equation (FPE) enables us to consider how synaptic interactions organize collective oscillations. We also develop the adjoint method (infinitesimal phase resetting curve) for simultaneous equations consisting of ordinary differential equations representing synaptic dynamics and a partial differential equation for determining the probability distribution of the membrane potential. This method provides a macroscopic phase response function (PRF), which gives insights into how it is modulated by external perturbation or internal changes of parameters. We investigate the effects of synaptic time constants and shunting inhibition on these gamma oscillations. The sensitivity of rising and decaying time constants is analysed in the oscillatory parameter regions; we find that these sensitivities are not largely dependent on rate of synaptic coupling but, rather, on current and noise intensity. Analyses of shunting inhibition reveal that it can affect both promotion and elimination of gamma oscillations. When the macroscopic oscillation is far from the bifurcation, shunting promotes the gamma oscillations and the PRF becomes flatter as the reversal potential of the synapse increases, indicating the insensitivity of gamma oscillations to perturbations. By contrast, when the macroscopic oscillation is near the bifurcation, shunting eliminates gamma oscillations and a stable firing state appears. More interestingly, under appropriate balance of parameters, two branches of bifurcation are found in our analysis of the FPE. In this case, shunting inhibition can effect both promotion and elimination of the gamma oscillation depending only on the reversal potential.