SenseLab Home ModelDB Home

Large cortex model with map-based neurons (Rulkov et al 2004)
Accession: 45525
We develop a new computationally efficient approach for the analysis of complex large-scale neurobiological networks. Its key element is the use of a new phenomenological model of a neuron capable of replicating important spike pattern characteristics and designed in the form of a system of difference equations (a map). ... Interconnected with synaptic currents these model neurons demonstrated responses very similar to those found with Hodgkin-Huxley models and in experiments. We illustrate the efficacy of this approach in simulations of one- and two-dimensional cortical network models consisting of regular spiking neurons and fast spiking interneurons to model sleep and activated states of the thalamocortical system. See paper for more.
Reference: Rulkov NF, Timofeev I, Bazhenov M (2004) Oscillations in large-scale cortical networks: map-based model. J Comput Neurosci 17:203-23 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type:  Network;
Brain Region(s)/Organism:  Neocortex;
Cell Type(s):  Neocortical pyramidal neuron: deep;  Neocortical Fast Spiking (FS) interneuron;
Channel(s):   
Gap Junctions:  
Receptor(s):  GabaA; AMPA;
Gene(s):  
Transmitter(s):  
Simulation Environment:  C or C++ program;
Model Concept(s):  Activity Patterns; Oscillations; Spatio-temporal Activity Patterns; Simplified Models; Sleep;
Implementer(s):  Bazhenov, Maxim [Bazhenov at Salk.edu]; Rulkov, Nikolai [nrulkov at ucsd.edu];
Search NeuronDB for information about:  Neocortical pyramidal neuron: deep; GabaA; AMPA;
Model files   Download zip file             Help downloading and running models
\
largectx2004
index.html
import2D.m
colors.m
input2D.txt
movieMAP.m
network2D.cpp
                            

README

for the model associated with:

Rulkov NF, Timofeev I, Bazhenov M (2004)
Oscillations in large-scale cortical networks: map-based model.
J Comput Neurosci 17:203-23

This model is written in C++. The following instructions work have worked under redhat linux 9 although the code should work with little or no modification on other platforms with ansi compilers.

How to compile:

Change to the directory in which you have unziped the archieve file, then type:
g++ network2D.cpp -lm -O4 -ffast-math -o n.exe

How to run:

n.exe input2D.txt > tmp &
C++ code for simulating 2 layer 2D network (256x256 pyramidal cells - regular spiking type and 128x128 interneurons - fast spiking type) of map-based neurons. With some small modifications this code was used to generate Fig 9A (right plot) in our JCNS paper. Also there is a movie of spiral wave dynamics available (see pop-down window above model files in ModelDB.

Making the movie with matlab:

Matlab scripts which create images which are then used to create the movie:
import2D('graf_cx0')
movieMAP
To assemble the images into the movie and then save the movie:
for i=1:200; tmp=imread(['mov' num2str(i) '.png']); m(i)=im2frame(tmp); end;
movie2avi(m,'spiral2D.avi');

Additional note:

The simulations with 256x256 require > 1GB RAM (otherwise there is a huge time penalty for virtual memory use).


Correspondence:
Dr. Maxim Bazhenov, The Salk Institute,
10010 North Torrey Pines Road, La Jolla, CA 92037.
E-mail: bazhenov@salk.edu

ModelDB Home  SenseLab Home   Help
Questions, comments, problems? Email the ModelDB Administrator
How to cite ModelDB
This site is Copyright 2012 Shepherd Lab, Yale University