Skip to content

Commit 45a65ca

Browse files
HengyeZhuHengyeZhu
committed
reason for not match
1 parent f1be0d6 commit 45a65ca

File tree

2 files changed

+302
-0
lines changed

2 files changed

+302
-0
lines changed

‎NeuroML2/compare_MC/RE/mod/HH2_magic.mod‎

Lines changed: 151 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,151 @@
1+
: $Id: HH2.mod,v 1.6 2004/06/06 16:00:30 billl Exp $
2+
TITLE Hippocampal HH channels
3+
:
4+
: Fast Na+ and K+ currents responsible for action potentials
5+
: Iterative equations
6+
:
7+
: Equations modified by Traub, for Hippocampal Pyramidal cells, in:
8+
: Traub & Miles, Neuronal Networks of the Hippocampus, Cambridge, 1991
9+
:
10+
: range variable vtraub adjust threshold
11+
:
12+
: Written by Alain Destexhe, Salk Institute, Aug 1992
13+
:
14+
: Modified Oct 96 for compatibility with Windows: trap low values of arguments
15+
:
16+
17+
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
18+
19+
NEURON {
20+
SUFFIX hh2ad
21+
USEION na READ ena WRITE ina
22+
USEION k READ ek WRITE ik
23+
RANGE gnabar, gkbar, vtraub, ikhh2, inahh2
24+
RANGE m_inf, h_inf, n_inf
25+
RANGE tau_m, tau_h, tau_n
26+
RANGE m_exp, h_exp, n_exp
27+
}
28+
29+
30+
UNITS {
31+
(mA) = (milliamp)
32+
(mV) = (millivolt)
33+
}
34+
35+
PARAMETER {
36+
gnabar = .003 (mho/cm2)
37+
gkbar = .005 (mho/cm2)
38+
39+
ena (mV)
40+
ek (mV)
41+
celsius (degC)
42+
dt (ms)
43+
v (mV)
44+
vtraub = -63 (mV)
45+
46+
47+
}
48+
49+
STATE {
50+
m h n
51+
}
52+
53+
ASSIGNED {
54+
ina (mA/cm2)
55+
ik (mA/cm2)
56+
il (mA/cm2)
57+
inahh2 (mA/cm2)
58+
ikhh2 (mA/cm2)
59+
m_inf
60+
h_inf
61+
n_inf
62+
tau_m
63+
tau_h
64+
tau_n
65+
m_exp
66+
h_exp
67+
n_exp
68+
tadj
69+
m_fcond
70+
h_fcond
71+
n_fcond
72+
}
73+
74+
75+
BREAKPOINT {
76+
SOLVE state METHOD cnexp
77+
78+
: gate fcond factors (similar to HH2_nondiff style)
79+
m_fcond = m*m*m
80+
h_fcond = h
81+
n_fcond = n*n*n*n
82+
: Use precomputed fcond factors (set in evaluate_fct)
83+
inahh2 = gnabar * m_fcond * h_fcond * (v - ena)
84+
ikhh2 = gkbar * n_fcond * (v - ek)
85+
ina = inahh2
86+
ik = ikhh2
87+
}
88+
89+
90+
DERIVATIVE state { : exact Hodgkin-Huxley equations
91+
evaluate_fct(v)
92+
m' = (m_inf - m) / tau_m
93+
h' = (h_inf - h) / tau_h
94+
n' = (n_inf - n) / tau_n
95+
}
96+
97+
UNITSOFF
98+
INITIAL {
99+
:
100+
: Q10 was assumed to be 3 for both currents
101+
:
102+
tadj = 3.0 ^ ((celsius-36)/ 10 )
103+
evaluate_fct(v)
104+
105+
m = m_inf
106+
h = h_inf
107+
n = n_inf
108+
}
109+
110+
PROCEDURE evaluate_fct(v(mV)) { LOCAL a,b,v2
111+
112+
v2 = v - vtraub : convert to traub convention
113+
114+
: a = 0.32 * (13-v2) / ( Exp((13-v2)/4) - 1)
115+
a = 0.32 * vtrap(13-v2, 4)
116+
: b = 0.28 * (v2-40) / ( Exp((v2-40)/5) - 1)
117+
b = 0.28 * vtrap(v2-40, 5)
118+
tau_m = 1 / (a + b) / tadj
119+
m_inf = a / (a + b)
120+
121+
a = 0.128 * Exp((17-v2)/18)
122+
b = 4 / ( 1 + Exp((40-v2)/5) )
123+
tau_h = 1 / (a + b) / tadj
124+
h_inf = a / (a + b)
125+
126+
: a = 0.032 * (15-v2) / ( Exp((15-v2)/5) - 1)
127+
a = 0.032 * vtrap(15-v2, 5)
128+
b = 0.5 * Exp((10-v2)/40)
129+
tau_n = 1 / (a + b) / tadj
130+
n_inf = a / (a + b)
131+
132+
m_exp = 1 - Exp(-dt/tau_m)
133+
h_exp = 1 - Exp(-dt/tau_h)
134+
n_exp = 1 - Exp(-dt/tau_n)
135+
136+
}
137+
FUNCTION vtrap(x,y) {
138+
if (fabs(x/y) < 1e-6) {
139+
vtrap = x/(Exp(x/y)-1)
140+
}else{
141+
vtrap = x/(Exp(x/y)-1)
142+
}
143+
}
144+
145+
FUNCTION Exp(x) {
146+
if (x < -100) {
147+
Exp = exp(x)
148+
}else{
149+
Exp = exp(x)
150+
}
151+
}

‎NeuroML2/compare_MC/RE/mod/HH2_reason.mod‎

Lines changed: 151 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,151 @@
1+
: $Id: HH2.mod,v 1.6 2004/06/06 16:00:30 billl Exp $
2+
TITLE Hippocampal HH channels
3+
:
4+
: Fast Na+ and K+ currents responsible for action potentials
5+
: Iterative equations
6+
:
7+
: Equations modified by Traub, for Hippocampal Pyramidal cells, in:
8+
: Traub & Miles, Neuronal Networks of the Hippocampus, Cambridge, 1991
9+
:
10+
: range variable vtraub adjust threshold
11+
:
12+
: Written by Alain Destexhe, Salk Institute, Aug 1992
13+
:
14+
: Modified Oct 96 for compatibility with Windows: trap low values of arguments
15+
:
16+
17+
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
18+
19+
NEURON {
20+
SUFFIX hh2ad_reason
21+
USEION na READ ena WRITE ina
22+
USEION k READ ek WRITE ik
23+
RANGE gnabar, gkbar, vtraub, ikhh2, inahh2
24+
RANGE m_inf, h_inf, n_inf
25+
RANGE tau_m, tau_h, tau_n
26+
RANGE m_exp, h_exp, n_exp
27+
}
28+
29+
30+
UNITS {
31+
(mA) = (milliamp)
32+
(mV) = (millivolt)
33+
}
34+
35+
PARAMETER {
36+
gnabar = .003 (mho/cm2)
37+
gkbar = .005 (mho/cm2)
38+
39+
ena (mV)
40+
ek (mV)
41+
celsius (degC)
42+
dt (ms)
43+
v (mV)
44+
vtraub = -63 (mV)
45+
46+
47+
}
48+
49+
STATE {
50+
m h n
51+
}
52+
53+
ASSIGNED {
54+
ina (mA/cm2)
55+
ik (mA/cm2)
56+
il (mA/cm2)
57+
inahh2 (mA/cm2)
58+
ikhh2 (mA/cm2)
59+
m_inf
60+
h_inf
61+
n_inf
62+
tau_m
63+
tau_h
64+
tau_n
65+
m_exp
66+
h_exp
67+
n_exp
68+
tadj
69+
m_fcond
70+
h_fcond
71+
n_fcond
72+
}
73+
74+
75+
BREAKPOINT {
76+
SOLVE state METHOD cnexp
77+
78+
: Use precomputed fcond factors (set in evaluate_fct)
79+
inahh2 = gnabar * m_fcond * h_fcond * (v - ena)
80+
ikhh2 = gkbar * n_fcond * (v - ek)
81+
ina = inahh2
82+
ik = ikhh2
83+
}
84+
85+
86+
DERIVATIVE state { : exact Hodgkin-Huxley equations
87+
evaluate_fct(v)
88+
m' = (m_inf - m) / tau_m
89+
h' = (h_inf - h) / tau_h
90+
n' = (n_inf - n) / tau_n
91+
}
92+
93+
UNITSOFF
94+
INITIAL {
95+
:
96+
: Q10 was assumed to be 3 for both currents
97+
:
98+
tadj = 3.0 ^ ((celsius-36)/ 10 )
99+
evaluate_fct(v)
100+
101+
m = m_inf
102+
h = h_inf
103+
n = n_inf
104+
}
105+
106+
PROCEDURE evaluate_fct(v(mV)) { LOCAL a,b,v2
107+
108+
v2 = v - vtraub : convert to traub convention
109+
110+
: a = 0.32 * (13-v2) / ( Exp((13-v2)/4) - 1)
111+
a = 0.32 * vtrap(13-v2, 4)
112+
: b = 0.28 * (v2-40) / ( Exp((v2-40)/5) - 1)
113+
b = 0.28 * vtrap(v2-40, 5)
114+
tau_m = 1 / (a + b) / tadj
115+
m_inf = a / (a + b)
116+
117+
a = 0.128 * Exp((17-v2)/18)
118+
b = 4 / ( 1 + Exp((40-v2)/5) )
119+
tau_h = 1 / (a + b) / tadj
120+
h_inf = a / (a + b)
121+
122+
: a = 0.032 * (15-v2) / ( Exp((15-v2)/5) - 1)
123+
a = 0.032 * vtrap(15-v2, 5)
124+
b = 0.5 * Exp((10-v2)/40)
125+
tau_n = 1 / (a + b) / tadj
126+
n_inf = a / (a + b)
127+
128+
m_exp = 1 - Exp(-dt/tau_m)
129+
h_exp = 1 - Exp(-dt/tau_h)
130+
n_exp = 1 - Exp(-dt/tau_n)
131+
132+
: gate fcond factors (similar to HH2_nondiff style)
133+
m_fcond = m*m*m
134+
h_fcond = h
135+
n_fcond = n*n*n*n
136+
}
137+
FUNCTION vtrap(x,y) {
138+
if (fabs(x/y) < 1e-6) {
139+
vtrap = x/(Exp(x/y)-1)
140+
}else{
141+
vtrap = x/(Exp(x/y)-1)
142+
}
143+
}
144+
145+
FUNCTION Exp(x) {
146+
if (x < -100) {
147+
Exp = exp(x)
148+
}else{
149+
Exp = exp(x)
150+
}
151+
}

0 commit comments

Comments
 (0)