EIF voltage traces and a recurrent network raster

13th May 2026

The same setup as nb000, but using exponential integrate-and-fire neurons in place of LIF.

The exponential integrate-and-fire (EIF) neuron is the natural next step up from LIF: same one-variable spirit, but with an explicit exponential term that models the upswing of a spike rather than a hard threshold-crossing rule.

The subthreshold dynamics are

τmdVdt=(VVrest)+ΔTexp ⁣(VVTΔT)+RmI,\tau_m \frac{dV}{dt} = -(V - V_\text{rest}) + \Delta_T\, \exp\!\left(\frac{V - V_T}{\Delta_T}\right) + R_m\, I,

with a peak-and-reset rule

V(t)Vpeak    spike at t,VVreset.V(t) \geq V_\text{peak} \;\Longrightarrow\; \text{spike at } t,\quad V \leftarrow V_\text{reset}.

Compared to LIF, the threshold VTV_T is no longer a discontinuity: it is the soft inflection point at which the exponential term takes over and the membrane potential blows up on its own. The slope factor ΔT\Delta_T controls how sharply this transition happens — small ΔT\Delta_T approximates LIF, larger ΔT\Delta_T smooths the spike initiation.

The code integrates with forward Euler at Δt=0.1ms\Delta t = 0.1\,\text{ms},

VV+Δtτm ⁣[(VVrest)+ΔTexp ⁣(VVTΔT)+RmI],V \leftarrow V + \frac{\Delta t}{\tau_m}\!\left[-(V - V_\text{rest}) + \Delta_T\, \exp\!\left(\tfrac{V - V_T}{\Delta_T}\right) + R_m\, I\right],

using τm=10ms\tau_m = 10\,\text{ms}, Vrest=65mVV_\text{rest} = -65\,\text{mV}, Vreset=70mVV_\text{reset} = -70\,\text{mV}, VT=50mVV_T = -50\,\text{mV}, ΔT=2mV\Delta_T = 2\,\text{mV}, Vpeak=0mVV_\text{peak} = 0\,\text{mV}, Rm=10MΩR_m = 10\,\text{M}\Omega, and tonic input I=2.5nAI = 2.5\,\text{nA}.

With the same tonic input as nb000, the EIF neuron fires at a different rate than the LIF neuron — the exponential term accelerates spike initiation, but the upswing itself takes time, so the period between spikes shifts.

EIF voltage trace

The network is structurally identical to nb000: 200 all-excitatory neurons, sparse pconn=0.1p_\text{conn} = 0.1 random connectivity, weights w=0.1w = 0.1, exponentially-decaying synaptic input with τsyn=5ms\tau_\text{syn} = 5\,\text{ms}, noisy per-neuron bias IibiasN(2.2,0.42)nAI_i^\text{bias} \sim \mathcal{N}(2.2, 0.4^2)\,\text{nA}, and a 2ms2\,\text{ms} refractory period after each spike. The only difference is that every neuron now obeys the EIF dynamics above.

EIF network raster and input current

Same plumbing as nb000 — neuron eif and neuron enet produce the figures, the runner copies them and writes numbers.json.

EIF run parameters

ParameterValue
commandeif
current2.5
duration100.0
dt0.1
firing_rate_hz60.0

EIF network run parameters

ParameterValue
commandenet
n200
duration500.0
dt0.1
seed0
mean_firing_rate_hz69.11
min_firing_rate_hz32.0
max_firing_rate_hz98.0