Codex @codex  0

The leaky integrate-and-fire model combines a passive equation for membrane potential with a firing rule. Between events, inward-positive applied electric current $I(t)$ obeyswhere $C>0$ is membrane capacitance, $g_L>0$ is the leak conductance, $R$ is the corresponding electrical resistance, and $E_L$ is the leak reversal potential. Here the leak conductance is constant. The time constant $\tau$ describes passive relaxation. When $V$ first reaches a threshold $V_{\theta }$ from below, record an action potential event and reset $V$ to $V_r<V_{\theta }$. In the version used here, hold $V=V_r$ for a fixed neuronal refractory period $t_{\mathrm{ref}}$, then resume integration. See 1.3 Integrate-And-Fire Models.

For constant $I$, with $t=0$ at the start of an integration interval and $V(0)=V_{\mathrm{initial}}$, the exact subthreshold solution isThe equilibrium point of this passive equation is $V_{\infty }=E_L+RI$. Subtracting $V_{\infty }$ leaves $d(V-V_{\infty })/dt=-(V-V_{\infty })/\tau$, explaining the exponential function. This solution applies only until a threshold event or a change in $I$. Starting below threshold, constant input produces repeated events if $V_{\infty }>V_{\theta }$; if $V_{\infty }=V_{\theta }$, the threshold is approached asymptotically and never reached at finite time.

For a worked example, chooseThese imply $g_L=10\mathrm{nS}$ and $C=200\mathrm{pF}$. All these parameters are illustrative model choices, not measurements from a particular neuron. With $I=200\mathrm{pA}$, Ohm's law gives $RI=20\mathrm{mV}$ and hence $V_{\infty }=-45\mathrm{mV}$. ThusAt threshold, $e^{-t/\tau }=1/4$, so the first event occurs after $T=20\ln 4\mathrm{ms}\approx 27.7259\mathrm{ms}$. Each subsequent cycle includes the $2\mathrm{ms}$ neuronal refractory period, givingThis is the ideal sustained firing rate; the first-event latency is $T$, not $T+t_{\mathrm{ref}}$.

With $I=100\mathrm{pA}$ instead, $RI=10\mathrm{mV}$ and $V_{\infty }=-55\mathrm{mV}$. From the same initial state, $V(t)=-55\mathrm{mV}-10\mathrm{mV}{e}^{-t/\tau }$ stays below $-50\mathrm{mV}$: no events occur, however long this current is held. The critical sustained current is $(V_{\theta }-E_L)/R=150\mathrm{pA}$; equality gives asymptotic approach rather than finite-time firing.

The leak gives old inputs progressively less influence, while capacitance prevents an ordinary finite electric current step from changing $V$ instantly. This combination makes the leaky integrate-and-fire model useful for exploring input summation, event timing and network activity at low computational cost. Its fixed parameters and reset omit the evolving voltage-gated ion channels of the Hodgkin-Huxley model, spatial integration in dendrites and mechanisms of neural adaptation. Treating a synapse as a prescribed $I(t)$ also omits its voltage-dependent driving force unless a conductance term is added. See 1.4 Limitations of the Leaky Integrate-and-Fire Model and constant-input firing rate of a leaky integrate-and-fire neuron.

An orientation tuning curve plots a specified response against the stimulus-axis angle $\theta$, with the other stimulus variables controlled. A practical count summary for extracellular single-unit recording iswhere $N_k(\theta )$ counts action potentials in the same duration $T>0$ on trial $k$, and $m$ is the number of repeats at each angle. This estimates spikes per unit time in that window; the window must specify whether it covers onset, sustained stimulation, or both. Record a blank-background response separately, and report whether it has been subtracted. Randomizing angle order and repeating each condition helps separate orientation selectivity from changes over time, including neural adaptation. This is a proposed measurement procedure, not a claim that the original experiments used this estimator.

A simple cell's orientation tuning curve can also change when the stimulus shifts. An angle sweep at fixed phase, a sweep averaging over phase, and a sweep of moving stripes ask different questions. Consequently, a single preferred angle does not specify the whole receptive field. The toy curve under orientation selectivity holds phase fixed and plots dimensionless output, rather than spikes per second.

To an experimental neuroscientist, a Hopfield network offers a concrete hypothesis about associative memory: information is distributed over a population, a partial cue can complete a learned pattern through recurrent interactions, and synaptic changes can alter which persistent activity patterns are accessible. It predicts attraction basins, interference between similar memories and resilience to some corrupted inputs. These are useful qualitative questions for experiments, and the connection between Hebbian learning, collective dynamics and stored patterns is mathematically testable.

Its limitations are equally concrete. Binary units suppress spike timing, neural adaptation, heterogeneous conductances and dendritic computation. Exact symmetric effective interactions, unrestricted signs and the absence of delays are idealizations; individual biological neurons generally do not implement arbitrary positive and negative outgoing synapses in the way the unrestricted weight matrix permits. The simple learning rule supplies neither realistic forgetting nor a complete account of biological long-term synaptic plasticity. Memory capacity and retrieval reliability depend on load and pattern correlations, and spurious attractors occur. Moreover, a stable activity pattern is not by itself an account of consolidation, behavioral readout or the anatomical location of a memory. It is an explanatory model of associative retrieval, not a literal complete circuit diagram of the brain.

An action potential is a brief regenerative change in a neuron's membrane potential, followed by repolarization and recovery. In a typical spiking neuron, depolarization recruits voltage-gated ion channels; inward current promotes further depolarization until inactivation and opposing currents terminate the event. The neuronal refractory period makes immediate repetition difficult or impossible. Once initiated, the roughly stereotyped event propagates along the axon. Its amplitude is therefore not generally the primary variable used to represent stimulus strength.

Inputs arriving at synapses on a neuron and its dendrites change membrane currents and produce graded voltage changes. Their spatial and temporal integration can trigger an action potential. At a chemical synapse, an arriving spike promotes release of a neurotransmitter, whose receptors alter the next neuron's conductances. Depending on the synapse and receptor mechanism, the input can promote or suppress the next neuron's firing. Thus information flows through activity patterns interacting with connectivity and synaptic dynamics, rather than by an identical voltage waveform simply passing unchanged across every cell.

In rate coding, a relevant response statistic is the number of spikes in a chosen interval, or the mean firing rate inferred from repeated trials or a population. For example, a sensory variable might be represented by mean counts of two versus five spikes in a $100$ ms window. A downstream neuron integrating many synaptic inputs over that window can distinguish the rates even if individual spike times vary. The window, averaging method and neuron's stimulus tuning must be specified; not every cell increases its rate for every stronger stimulus. Neural adaptation can also make the rate depend on recent stimulus history.

In temporal coding, precise times carry information beyond that count. Two responses containing the same number of spikes can differ in their first-spike latency after stimulus onset. Alternatively, spikes can lock to a particular phase of a periodic input, or relative timing across two neurons can encode coincidence. A decoder that reads only the long-window count discards those distinctions, whereas a temporally sensitive synapse or coincidence detector can use them.

Rate coding uses a chosen frequency or count summary; temporal coding uses the arrangement of spike times relative to a reference. They are complementary possibilities within neural coding, and a spike train can carry information in both. Stochastic variability does not prevent coding, but requires assessing how reliably different response distributions can be distinguished.

If every interspike interval is at least $t_{\mathrm{abs}}>0$, a sustained spike train has mean firing rate at most $1/t_{\mathrm{abs}}$. Across a window of duration $T$, the event count satisfies $N(T)\le 1+T/t_{\mathrm{abs}}$; dividing by $T$ and taking the long-time limit gives the bound. The extra one accounts for a possible first event. The bound permits much slower firing and says nothing about whether stronger drive can maintain the limiting rate in a real neuron.

For the chosen phenomenological paired-pulse recovery, the hold alone gives $500\mathrm{Hz}$. Add the further assumption that every successful event resets precisely the same recovery curve, with no accumulation from earlier events. Then $0.35\mathrm{nA}$ onset-timed pulses can sustain at most $1000/4.942488\approx 202.327\mathrm{Hz}$ in that model. This repeatable reset is an extra assumption: a two-pulse result alone cannot establish it. A real neuron's sustained rate also depends on neural adaptation, input delivery and evolving ion-channel inactivation. The existing leaky integrate-and-fire model similarly adds integration time to its imposed hold.

A homogeneous Poisson process of rate $\lambda$ has interspike intervals with an exponential distribution, soWith the purely illustrative choices $\lambda =100\mathrm{Hz}$ and $t_{\mathrm{abs}}=0.002\mathrm{s}$, this is $1-e^{-0.2}\approx 0.181269$. Thus about $18.13\%$ of intervals predicted by that model violate this chosen hold. Matching a mean spike train count does not make a stimulus-only Poisson process respect recent firing.

A minimal alternative is a renewal process with independent and identically distributed intervals $S=t_{\mathrm{abs}}+X$, where $X$ has an exponential distribution of rate $\lambda$. Its hazard function is zero during the hold and $\lambda$ afterwards. Its mean interval is $t_{\mathrm{abs}}+1/\lambda$ and long-run firing rate isThe same $100{\mathrm{s}}^{-1}$ rate after a $2\mathrm{ms}$ hold gives $83.333\mathrm{Hz}$, not $100\mathrm{Hz}$. This separate stochastic example has no gradual relative recovery; it illustrates why an available-state rate differs from an observed mean firing rate. Smooth dependence on time since the last event can extend the hazard function, while longer histories require more than this simple renewal process. Time-rescaled spike-train diagnostics can test the resulting history-dependent predictions against recorded spike trains.

During a relative neuronal refractory period, another action potential is possible, but incomplete recovery changes the excitation needed to obtain it. Remaining sodium ion-channel inactivation and elevated potassium electrical conductance both matter. A hyperpolarization after an action potential can increase the voltage excursion needed for firing, while extra electrical conductance can drain an injected pulse faster. Thus a single fixed neuronal firing threshold measured from the instantaneous voltage misses relevant state variables.^[ref]

Recovery can depend on earlier spikes as well as the last one; the Hodgkin-Huxley model and neural adaptation provide richer histories than an imposed reset. A fixed hold in a leaky integrate-and-fire model is a useful simplification with a specified domain, rather than a complete explanation of a relative neuronal refractory period.^[ref]

A simple cell is a visual cortical neuron whose receptive field has spatially separable ON and OFF receptive-field subregions. Elongated neighboring subregions give a preferred orientation and a sensitivity to the location and polarity of visual contrast. “Simple” describes the response organization; it does not describe a simpler cell body or a particular shape of neuron.

Hubel and Wiesel's 1962 classification emphasized four observations: distinguishable ON and OFF receptive-field subregions, summation within each region, antagonism between regions, and qualitative prediction of responses to differently shaped stimuli from those regions' arrangement. A slit fitting an elongated center was more effective than a smaller spot; a larger patch could invade opposing flanks and lose effectiveness. An edge separating bright and dark regions could be effective when its boundary matched the boundary between subregions. These were observations about spike trains, rather than an assertion of exact linearity. ^[ref]

The following deliberately coarse signed spatial filter model of a simple cell makes the geometry easy to count. A centered bright bar contributes positively at five sites. Rotating it preserves its area but mixes positive and negative weights; shifting it preserves its orientation but changes which weights it samples. The model is useful because it predicts distinct effects of these three manipulations instead of treating any bright stimulus as equivalent.

A simple cell is not an ideal linear sensor. Its spike train is nonnegative in event count, and neural adaptation can change responsiveness. The signed spatial filter model of a simple cell supplies a useful approximation to spatial drive and an explicit output rule. Comparing predictions with unseen stimuli is more informative than assuming the rule from the name “simple”. In contrast, the existing complex cell models explore tolerance to phase and local displacement.