Table 2.

Equations of the NL neuron model

Variable	Equation
Compartment	c = soma or node
Membrane potential
Na current	I_Na^c = ḡ_Na^c · m_c h_c · (E_Na − V_c)
KHVA current	I_KHVA^c = ḡ_KHVA^c · n_c · (E_K − V_c)
KLVA current	I_KLVA^c = ḡ_KLVA^c · d_c · (E_K − V_c)
Leak current	I_leak^c = g_leak^c · (E_leak − V_c)
Synaptic input	I_synaptic^c = g_synaptic^c · (E_synaptic − V_c)
Axonal current	I_axon^soma = g_axon · (V_node − V_soma), I_axon^node = −I_axon^soma
Channel variable (x = m, h, n, d)	dx_c/dt = φ{α_x(V_c) · (1 − x_c) − β_x(V_c) · x_c}
Temperature dependence	φ = Q₁₀^{(T−23.0)/10}(Q₁₀ = 2.5, T = 40°C)
Na channel activation	α_m(V) = 3.6exp((V + 34)/7.5), β_m(V) = 3.6exp(−(V + 34)/10.0)
Na channel inactivation	α_h(V) = 0.6exp(−(V + 57)/18.0), β_h(V) = 0.6exp((V + 57)/13.5)
KHVA channel activation	α_n(V) = 0.11exp((V + 19)/9.1), β_n(V) = 0.103exp(−(V + 19)/20)
KLVA channel activation	α_d(V) = 0.2exp((V + 60)/21.8), β_d(V) = 0.17exp(−(V + 60)/14)

This table summarizes the equations used in our simulation. In the single-compartment model, I_axon^node = I_axon^soma = 0 and only the somatic compartment was considered.