Note
Go to the end to download the full example code.
Population of GIF neuron model with oscillatory behavior¶
Run this example as a Jupyter notebook:
See our guide for more information and troubleshooting.
This script simulates a population of generalized integrate-and-fire (GIF) model neurons driven by noise from a group of Poisson generators.
Due to spike-frequency adaptation, the GIF neurons tend to show oscillatory behavior on the time scale comparable with the time constant of adaptation elements (stc and sfa).
Population dynamics are visualized by raster plot and as average firing rate.
References¶
Import all necessary modules for simulation and plotting.
import matplotlib.pyplot as plt
import nest
import nest.raster_plot
nest.ResetKernel()
Assigning the simulation parameters to variables.
dt = 0.1
simtime = 2000.0
Definition of neural parameters for the GIF model. These parameters are extracted by fitting the model to experimental data [2].
neuron_params = {
"C_m": 83.1,
"g_L": 3.7,
"E_L": -67.0,
"Delta_V": 1.4,
"V_T_star": -39.6,
"t_ref": 4.0,
"V_reset": -36.7,
"lambda_0": 1.0,
"q_stc": [56.7, -6.9],
"tau_stc": [57.8, 218.2],
"q_sfa": [11.7, 1.8],
"tau_sfa": [53.8, 640.0],
"tau_syn_ex": 10.0,
}
Definition of the parameters for the population of GIF neurons.
N_ex = 100 # size of the population
p_ex = 0.3 # connection probability inside the population
w_ex = 30.0 # synaptic weights inside the population (pA)
Definition of the parameters for the Poisson group and its connection with GIF neurons population.
N_noise = 50 # size of Poisson group
rate_noise = 10.0 # firing rate of Poisson neurons (Hz)
w_noise = 20.0 # synaptic weights from Poisson to population neurons (pA)
Configuration of the simulation kernel with the previously defined time resolution.
nest.resolution = dt
Building a population of GIF neurons, a group of Poisson neurons and a spike recorder device for capturing spike times of the population.
population = nest.Create("gif_psc_exp", N_ex, params=neuron_params)
noise = nest.Create("poisson_generator", N_noise, params={"rate": rate_noise})
spike_det = nest.Create("spike_recorder")
Build connections inside the population of GIF neurons population, between Poisson group and the population, and also connecting spike recorder to the population.
nest.Connect(population, population, {"rule": "pairwise_bernoulli", "p": p_ex}, syn_spec={"weight": w_ex})
nest.Connect(noise, population, "all_to_all", syn_spec={"weight": w_noise})
nest.Connect(population, spike_det)
Simulation of the network.
nest.Simulate(simtime)
Plotting the results of simulation including raster plot and histogram of population activity.
nest.raster_plot.from_device(spike_det, hist=True)
plt.title("Population dynamics")
plt.show()