I am currently attempting to model an AMPA synapse as a point process. I used an existing model

(http://senselab.med.yale.edu/modeldb/Sh ... model=3808)

as an example and template for my code.

Code: Select all

```
TITLE Ampa2
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
NEURON {
POINT_PROCESS AmpaSyn
RANGE onset, gMax, e, i, g, tau
NONSPECIFIC_CURRENT i
}
UNITS {
(nA) = (nanoamp)
(mV) = (millivolt)
(nS) = (nanomho)
}
PARAMETER {
onset (ms)
gMax = 0 (nS)
tau = 18 (ms)
e = 0.0 (mV)
v (mV)
}
ASSIGNED {
i (nA)
g (nS)
}
UNITSOFF
BREAKPOINT { LOCAL tt
tt = (t-onset)
if (t>onset) {
g = gMax * exp(-tt/tau) * (1-exp(-tt/2.2))/0.68
}
else {g = 0}
i = g * (v - e)
}
UNITSON
```

*seems*okay until I started playing with some values of gMax, and found that the resulting current values were a bit off.

Now if I'm correct with my calculations, I should have a current with an amplitude of -73pA, which doesn't seem to be the case. The g values seem to check out well, so I'm not sure what my problem is. Any suggestions?

Thanks!