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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
// genesis
create neutral /library
// We don't want the library to try to calculate anything,
// so we disable it.
disable /library
//========================================================================
// CHANNELS
//========================================================================
create synchan /library/Na_channel
setfield /library/Na_channel \
Ek {ENA} \
tau1 { 2.0e-3 } \ // sec
tau2 { 2.0e-3 } \ // sec
gmax {DISTAL_GMAX_NA}
create synchan /library/K_channel
setfield /library/K_channel \
Ek { EK } \
tau1 { 100.0e-3 } \ // sec
tau2 { 100.0e-3 } \ // sec
gmax {GMAX_K} // S
create synchan /library/Cl_channel
setfield /library/Cl_channel \
Ek {ECL} \
tau1 { 4.0e-3 } \ // sec
tau2 { 4.0e-3 } \ // sec
gmax {GMAX_CL} // S
//========================================================================
// COMPARTMENT
//========================================================================
create compartment /library/compartment
setfield /library/compartment \
Cm {CM*SOMA_A} \ // F
Ra {RA*SOMA_L/SOMA_XA} \ // ohm
Em {EREST_ACT} \ // V
Rm {RM/SOMA_A} \ // ohm
inject 0.0
//========================================================================
// ACTIVE NA CHANNEL
//========================================================================
create hh_channel /library/HH_Na_channel
setfield /library/HH_Na_channel \
Ek {ENA_ACT} \ // V
Gbar { 1.2e3 * SOMA_A } \ // S
Xpower 3.0 \
Ypower 1.0 \
X_alpha_FORM {LINOID} \
X_alpha_A -0.1e6 \ // 1/V-sec
X_alpha_B -0.010 \ // V
X_alpha_V0 { 0.025 + EREST_ACT } \ // V
X_beta_FORM {EXPONENTIAL} \
X_beta_A 4.0e3 \ // 1/sec
X_beta_B -18.0e-3 \ // V
X_beta_V0 { 0.0 + EREST_ACT } \ // V
Y_alpha_FORM {EXPONENTIAL} \
Y_alpha_A 70.0 \ // 1/sec
Y_alpha_B -20.0e-3 \ // V
Y_alpha_V0 { 0.0 + EREST_ACT } \ // V
Y_beta_FORM {SIGMOID} \
Y_beta_A 1.0e3 \ // 1/sec
Y_beta_B -10.0e-3 \ // V
Y_beta_V0 { 30.0e-3 + EREST_ACT } // V
//========================================================================
// ACTIVE K CHANNEL
//========================================================================
create hh_channel /library/HH_K_channel
setfield /library/HH_K_channel \
Ek {EK} \ // V
Gbar {360.0*SOMA_A} \ // S
Xpower 4.0 \
Ypower 0.0 \
X_alpha_FORM {LINOID} \
X_alpha_A -10.0e3 \ // 1/V-sec
X_alpha_B -10.0e-3 \ // V
X_alpha_V0 {10.0e-3+EREST_ACT} \ // V
X_beta_FORM {EXPONENTIAL} \
X_beta_A 125.0 \ // 1/sec
X_beta_B -80.0e-3 \ // V
X_beta_V0 {0.0+EREST_ACT} // V
//========================================================================
// SPIKE DETECTOR
//========================================================================
create spikegen /library/spike
setfield /library/spike \
thresh -40e-3 \ // V
abs_refract { 10e-3 } \ // sec
output_amp 1
|
