File: examples-frompapers_Wang_Buszaki_1996.txt

package info (click to toggle)
brian 1.4.3-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, stretch
  • size: 23,436 kB
  • sloc: python: 68,707; cpp: 29,040; ansic: 5,182; sh: 111; makefile: 61
file content (63 lines) | stat: -rw-r--r-- 1,538 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
.. currentmodule:: brian

.. index::
   pair: example usage; NeuronGroup
   pair: example usage; plot
   pair: example usage; StateMonitor
   pair: example usage; run
   pair: example usage; show

.. _example-frompapers_Wang_Buszaki_1996:

Example: Wang_Buszaki_1996 (frompapers)
=======================================

Wang-Buszaki model
------------------

J Neurosci. 1996 Oct 15;16(20):6402-13.
Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network model.
Wang XJ, Buzsaki G.

Note that implicit integration (exponential Euler) cannot be used, and therefore
simulation is rather slow.

::

    from brian import *
    
    defaultclock.dt=0.01*ms
    
    Cm=1*uF # /cm**2
    Iapp=2*uA
    gL=0.1*msiemens
    EL=-65*mV
    ENa=55*mV
    EK=-90*mV
    gNa=35*msiemens
    gK=9*msiemens
    
    eqs='''
    dv/dt=(-gNa*m**3*h*(v-ENa)-gK*n**4*(v-EK)-gL*(v-EL)+Iapp)/Cm : volt
    m=alpham/(alpham+betam) : 1
    alpham=-0.1/mV*(v+35*mV)/(exp(-0.1/mV*(v+35*mV))-1)/ms : Hz
    betam=4*exp(-(v+60*mV)/(18*mV))/ms : Hz
    dh/dt=5*(alphah*(1-h)-betah*h) : 1
    alphah=0.07*exp(-(v+58*mV)/(20*mV))/ms : Hz
    betah=1./(exp(-0.1/mV*(v+28*mV))+1)/ms : Hz
    dn/dt=5*(alphan*(1-n)-betan*n) : 1
    alphan=-0.01/mV*(v+34*mV)/(exp(-0.1/mV*(v+34*mV))-1)/ms : Hz
    betan=0.125*exp(-(v+44*mV)/(80*mV))/ms : Hz
    '''
    
    neuron=NeuronGroup(1,eqs)
    neuron.v=-70*mV
    neuron.h=1
    M=StateMonitor(neuron,'v',record=0)
    
    run(100*ms,report='text')
    
    plot(M.times/ms,M[0]/mV)
    show()