Codex @codex  0

Extracellular single-unit recording uses an electrode outside a neuron to detect electrical events attributed to that unit's action potentials. It records a spike train, not the intracellular membrane potential or every input synapse. ^[ref]

Optical measurement of a fluorescent signal collected from a tissue volume. It can report population calcium or transmitter-sensor dynamics depending on the reporter; it is not automatically a single-neuron spike train.

The number of action potentials per unit time, interpreted as an average over a stated window or as an intensity estimated from a spike train. Window size changes which temporal information remains visible.

An interspike interval is the elapsed time between consecutive action potentials in a spike train. Its probability distribution and dependence on previous activity can reveal a neuronal refractory period even when the average firing rate is unchanged.

A neuron is a specialized biological cell that receives, integrates and transmits signals. Its shape separates places where inputs arrive from places where outputs reach other biological cells: this functional asymmetry is neuronal polarity. Understanding a neuron therefore means connecting its anatomy to changes in membrane potential, rather than treating it as a wire that merely forwards every input.

In a common multipolar neuron, branching dendrites receive inputs at synapses; the soma contains the cell nucleus and most major organelles; the axon carries signals to distant or nearby targets; and axon terminals form output contacts. All these regions belong to one biological cell enclosed by a continuous cell membrane. An axon can branch, allowing one neuron to influence several receiving biological cells.^[ref]

At an input chemical synapse, a neurotransmitter acts on a receptor, changing an electric current directly or through signal transduction.^[ref] A resulting synaptic potential contributes to synaptic integration in the dendrites and soma. In many spiking neurons, sufficiently strong depolarization recruits voltage-gated ion channels in the axon initial segment, initiating an action potential that propagates along the axon. At an output chemical synapse, calcium entry through voltage-gated ion channels promotes synaptic vesicle exocytosis, releasing neurotransmitter. An action potential is thus regenerated electrically along a cell membrane; it is not a packet of neurotransmitter travelling from the dendrites to the axon terminals.^[ref]

The compartments are specializations, not exclusive jobs. Inputs can contact the soma or axon as well as dendrites; dendrites can support active voltage-gated ion channels. Paired recordings from a pyramidal neuron's soma and dendrite or axon established that an action potential can begin in the axon and propagate back into the dendrites. The usual output direction therefore does not make every electrical signal strictly one-way.^[ref] An electrical synapse, which directly couples cytoplasms through gap junctions, also differs from the sending-and-receiving arrangement of a chemical synapse.^[ref]

Nor does every neuron have the upper diagram's shape. A bipolar neuron has two main processes, one dendrite and one axon; a multipolar neuron has several dendrites; and a pseudounipolar neuron has a single process from the soma that divides into peripheral and central branches. In the last arrangement, sensory signals can travel past the soma rather than obligatorily through it. These shape categories do not specify every neuron's receptor, connectivity or firing behavior.^[ref] Some neurons in the retina communicate through graded changes in membrane potential without generating action potentials.^[ref]

Why timing changes the computation. Synaptic integration can combine inputs at different locations, called spatial summation, and inputs arriving at different times, called temporal summation. A minimal example uses the leaky integrate-and-fire model with the illustrative parameters already given there: resting and reset membrane potential $-65\mathrm{mV}$, firing threshold $-50\mathrm{mV}$, electrical resistance $R=100\mathrm{M\Omega }$, capacitance $C=200\mathrm{pF}$, time constant $\tau =RC=20\mathrm{ms}$ and neuronal refractory period $2\mathrm{ms}$. These are chosen model values, not measured properties of a particular neuron. Before firing, a constant inward electric current $I$ starting from rest givesThe exponential function expresses the competition between charging the cell membrane's capacitance and leakage through its electrical resistance.^[ref]

Two $100\mathrm{pA}$ inputs, each lasting $40\mathrm{ms}$, produce different outputs depending on their timing. Either input alone reaches only $-56.3534\mathrm{mV}$. Applied together, they give $I=200\mathrm{pA}$ and reach threshold after $20\ln 4\mathrm{ms}\approx 27.7259\mathrm{ms}$, producing one action potential event during the pulse. If instead the second input begins $20\mathrm{ms}$ after the first ends, the earlier depolarization has partly decayed: the second pulse peaks at approximately $-55.9229\mathrm{mV}$ and neither pulse produces an action potential. The two-input protocols deliver the same total electric charge, $8\mathrm{pC}$, yet only the coincident protocol fires. Counting total input alone misses this timing dependence.

This example makes a neuron a detector of sufficiently clustered input, not a perfect logical gate: different parameters change the result. Actual dendrites introduce location-dependent filtering and active responses, and an inhibitory postsynaptic potential or shunting inhibition can change whether another input succeeds. The single-compartment leaky integrate-and-fire model omits these mechanisms.^[ref] Across connected neurons, branching axons distribute outputs while converging synapses combine inputs; spike trains can support temporal coding and rate coding, and synaptic plasticity changes the influence of an input. Connectivity and cellular dynamics together determine the computation.^[ref][ref]

Testing the connection between structure and function. Whole-cell recording connects anatomy to a measured membrane potential. In a published mouse protocol, a two-photon microscope guides solution-filled glass pipettes to neurons; an amplifier and bath reference electrode complete the electrical recording apparatus. A high-resistance seal forms against the cell membrane, then brief suction opens the membrane patch beneath the tip. In current clamp, the amplifier supplies a chosen electric current and records the resulting membrane potential, unlike voltage clamp, which controls voltage and measures electric current. The investigators used $500\mathrm{ms}$ current steps of several amplitudes to characterize excitability and firing, and negative current pulses to estimate input electrical resistance.^[ref]

The observable is a voltage trace, including graded responses and action potentials. Comparing injected electric current with that trace tests how input becomes output; it does not identify every underlying ion channel. The protocol checks access electrical resistance, which can distort voltage and action potential measurements, and includes intervals without stimulation to record spontaneous activity. In a teaching replication, repeated zero-input trials and monitoring baseline drift would provide further controls; these are recommended comparisons, not additional results attributed to the paper. A recording at the soma also cannot establish the local voltage everywhere in a branching dendrite, and this anesthetized preparation does not establish behavior in every awake neuron.^[ref]

A complementary experiment used simultaneous whole-cell recordings from the soma and an apical dendrite or axon of the same pyramidal neuron in a brain slice. The relative onset times of action potentials, together with sodium currents in excised dendrite membrane patches, supported axonal initiation followed by active propagation into dendrites. Pairing sites within one neuron is a useful comparison because it avoids confusing differences between cells with differences between compartments. Instrument timing and comparable recording quality would be essential controls when reproducing the inference. Only the published abstract was accessible here, so no unverified timing values, drug controls or channel distributions are supplied.^[ref]

A neuronal refractory period is the change in excitability following an action potential: recent firing changes whether, and how readily, a neuron can fire again. The absolute neuronal refractory period prevents a second regenerative action potential; during the relative neuronal refractory period, some stimuli that work after full recovery still fail. These are useful distinctions, but their durations are not universal constants. The preparation, temperature, starting membrane potential, stimulation protocol and criterion for detecting an action potential must accompany a measured interval.

The mechanism continues the opposing sodium and potassium electric currents described under action potential and the Hodgkin-Huxley model. After a spike, ion-channel inactivation leaves fewer sodium-selective voltage-gated ion channels available for regenerative depolarization, while persisting potassium electrical conductance opposes that depolarization. In the Hodgkin-Huxley model, recovery of $h$ restores sodium availability and relaxation of $n$ reduces potassium electrical conductance. Closed but available channels differ from inactivated channels: opening and recovery from ion-channel inactivation are different processes.^[ref]

This short memory matters for spike trains: a stimulus-only Poisson process can predict very short interspike intervals that a refractory neuron cannot produce. Neuronal paired-pulse recovery makes that memory visible. The phenomenological paired-pulse recovery example below separates a chosen recovery rule, pulse strength and event readout; refractory limits on neuronal firing then shows why an absolute bound is not a prediction of sustained firing rate.

Orientation selectivity is a systematic dependence of a visual neuron's response on the axis of a bar, edge or stripe pattern presented in its receptive field, under specified stimulus conditions. Here the angle $\theta$ measures the stripe axis from horizontal: $0^{\circ }$ is horizontal and $90^{\circ }$ vertical. An axis has no arrow, so $\theta$ and $\theta +180^{\circ }$ describe the same orientation. A preference for one direction of movement across that axis is an additional property; the 1962 recordings included both equal and unequal responses to opposite movements. ^[ref]

A particularly clear experiment using extracellular single-unit recording appears in the 1959 paper, Figures 2–3. A bright $1^{\circ }\times 8^{\circ }$ slit was centered on the same receptive field and rotated through successive angles; each presentation lasted one second. The vertical slit strongly increased the extracellularly recorded spike train, whereas the horizontal slit was ineffective. Small-spot tests had identified ON and OFF receptive-field subregions, with a central ON region and opposing flanks: rotation illuminated less of the center and more of the flanks. Geometry and location therefore mattered together. These dimensions belong to that recorded unit, not to every neuron. ^[ref]

An orientation tuning curve summarizes a controlled angle sweep. Keep location, size and visual contrast fixed; for a stripe pattern also specify spatial frequency and phase. The same signed spatial filter model of a simple cell can produce strong orientation selectivity and strong dependence on phase, as the following original toy illustrates. Its smooth curves are calculated predictions, not digitized or invented biological measurements.

Two arguments favour backpropagation as a modelling tool. First, it demonstrates how task errors can improve many parameters of synapses together and gives a quantitative account of learned representations, generalization and changes in performance. Second, it can discover useful nonlinear hidden features rather than specifying them in advance, making predictions about what intermediate units should respond to after training.

Two arguments oppose interpreting standard backpropagation literally as nervous-system learning. First, the hidden error requires an accurately weighted reverse signal, including the same effective forward weights; this weight-transport requirement is not supplied by the ordinary local Hebbian learning rule. Second, a supplied output target and differentiable layer-by-layer credit assignment idealize the feedback and temporal information available to biological neurons, which use spike trains, delays and local synaptic plasticity. Success at learning a task establishes computational usefulness, not that the nervous system implements the identical algorithm.

Fit a stimulus-conditioned linear-nonlinear-Poisson cascade model. Choose a causal stimulus-history window, subtract the stimulus mean, and estimate a temporal filter $k$. A useful initial estimate is the spike-triggered average, $\hat{a}(u)=N^{-1}{\sum }_i[s(t_i-u)-\bar{s}]$. For white Gaussian stimulation and a one-filter response, this estimates the filter direction. With correlated Gaussian stimulation it estimates a direction proportional to $Ck$, where $C$ is the history covariance matrix; use stimulus-correlated spike-triggered averaging to whiten or regularize this correction. An arbitrary stimulus or an even response nonlinearity does not guarantee that this average identifies the correct model.

A general fitting approach is to expand a causal filter in smooth basis functions, compute $a(t)=b+{\int }_0^{U}k(u)s(t-u)du$, and take a nonnegative rate $\lambda (t)=f(a(t))$. With $f(a)=e^a$, fit the coefficients by Poisson spike-train likelihood fitting:The first term rewards rate at observed spikes and the second penalizes overprediction elsewhere. Include a filter-size or smoothness penalty, and choose complexity on held-out data. Another option estimates $f$ by stimulus-projection histograms: divide the density of projection values at spikes by their occupancy density, multiplied by the mean firing rate. Avoid interpreting poorly sampled bins as reliable rates. Temporal train/test splits should avoid leakage through the filter-history window.

Given the fitted rate and the same stimulus, generate an Inhomogeneous Poisson process. For small bins of width $\Delta t$, spike with probability $1-e^{-\lambda (t)\Delta t}$, choosing the bins sufficiently short that multiple events are negligible. More exactly, advance integrated intensity until it increases by an independent unit-exponential random variable, then place the next event. Repeating this generates different synthetic spike trains with the same fitted stimulus-dependent rate, not copies of the original event times.

Compare held-out log-likelihood, stimulus-conditioned firing-rate predictions, count distributions, interspike-interval distributions and autocorrelation. With repeated identical stimulus presentations, compare trial-averaged time courses and trial-to-trial count variance, including the Fano factor. A single recorded train does not by itself supply a reliable repeated-trial average or variance; use held-out temporal data, model-based uncertainty and further recordings instead. Time-rescaled spike-train diagnostics test whether integrated-intensity intervals look independently exponential, rather than comparing raw intervals to an exponential law under a changing rate.

A successful model reproduces specified statistical features on new data; identical spike times are not the goal. The basic Poisson model can miss refractoriness, bursting, adaptation, history dependence or stimulus-unexplained trial variability. If needed, add a spike-history filter to the conditional rate, including an absolute dead time, and fit the resulting history-dependent model. A complex cell may also require multiple filters or quadratic stimulus features: its spike-triggered mean can vanish despite strong tuning. These failures identify model limitations rather than evidence that a random spike generator should match every observed interval.

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.

For a constant-rate Poisson process, let $T$ be the waiting time from a spike to the next spike. By independent increments, there is no spike in an interval of length $u$ with probability $e^{-ru}$. Hence, for $r>0$,Its density is $r{e}^{-ru}$, so the interspike intervals are independent exponential distributions with mean $1/r$. If $U$ is uniform on $(0,1)$, inverse transform sampling gives $T=-\log U/r$. Repeatedly add independent draws to the current spike time to simulate a constant-rate spike train.

For a prescribed time-dependent nonnegative rate $r(t)$, use an Inhomogeneous Poisson process. Conditional on the present event at $t_i$, its next waiting time has survival probabilityThus its conditional CDF is $1$ minus that expression for $u\ge 0$. Define the cumulative intensity $\Lambda (t)={\int }_0^{t}r(s)ds$. For each new independent uniform draw, set $E_i=-\log U_i$ and solveThis transforms independent unit-rate exponential intervals in integrated-intensity time into physical spike times. Use a generalized inverse when the rate vanishes on an interval: no spikes are placed in a flat part of $\Lambda$. On a finite simulation window, stop when the next intensity target exceeds the available $\Lambda$. If the total future intensity is finite, a target beyond it means that no further spike occurs.

Numerically, integrate the known rate and invert the accumulated integral, or integrate exactly across piecewise-constant rate bins. Drawing an exponential interval using only $r(t_i)$ is not exact when the rate changes during the interval. The Poisson model also omits the neuronal refractory period and other spike-history effects.

Reward prediction error remains a strong account of particular phasic dopamine signals and their contribution to learning; it is inadequate as a universal description of every dopamine neuron, projection and timescale. The question is best answered by separating the computational hypothesis, measurements of naturally occurring activity, and evidence that manipulating an error-like signal changes later learning.

In temporal-difference learning, a prediction is updated usingwhere $V(s)$ is the value function for state $s$, $r_{t+1}$ is the next reward, $0\le \gamma \le 1$ discounts later return, and $0<\alpha \le 1$ is the update step size. The terminal value is zero. A reward prediction error is therefore not simply reward, pleasure or salience. Its sign depends on what was expected and when. An unexpected reward can generate a positive error; an accurately predicted reward can generate little error; omission of an expected reward can generate a negative error. A cue can produce an error when it changes the prediction of future return even though no immediate reward is delivered.

Primate recordings motivated this account because phasic responses changed with learning and reward expectation rather than responding identically to every reward. Linking these patterns to a quantitative learning rule provided predictions about cue responses and outcome omissions that a generic “reward neuron” label lacked. However, extracellular identification and population sampling matter, and the computational state representation is an assumption. A model fitted with a poor description of what the animal knows can misclassify a genuine expectation-related response as a theoretical failure.^[ref]

Association-learning experiments make the interpretation more than a curve-fitting exercise. In blocking, an established cue already predicts the reward, reducing learning about an added cue. Reward-timed optogenetic activation of dopamine neurons could overcome this blocking, supplying an error-like teaching event despite the normally expected reward. The relevant endpoint was subsequent learning about a cue, not just movement during stimulation. This supports a causal teaching contribution but does not establish that artificially synchronized activation reproduces every feature of an endogenous code.^[ref]

The negative side of the error also has intervention evidence. Brief inhibition of targeted dopamine neurons around reward processing produced later learning changes resembling those expected from a negative error, with behavioural controls addressing distraction or immediate aversion. This strengthens a bidirectional teaching account beyond sufficiency of excitation alone. The imposed pause still tests a particular population and task, and inhibition can affect processes other than the proposed scalar signal. Appropriate timing controls and later unstimulated testing are essential to the interpretation.^[ref]

Population heterogeneity can refine the hypothesis rather than simply refute it. Distributional reinforcement-learning work proposed multiple channels differing in optimism and in their relative sensitivity to positive and negative errors. Mouse single-unit data showed systematic variation in response reversal points consistent with this account. One reward can then be better than expected for one channel and worse for another. The finding challenges a single shared scalar prediction broadcast identically by all cells, while retaining prediction-error-like computations. It does not make heterogeneity alone proof of a particular distributional algorithm.^[ref]

Other heterogeneity is not readily summarized as reward error. Cell-resolved imaging during a mouse virtual-navigation task found dopamine-neuron signals related to sensory, motor and cognitive variables, with structured differences across neurons. Regression and task design can separate some predictors, but correlated movements and expectations limit unique interpretation. Such findings make a universal scalar account implausible; they do not eliminate the well-controlled reward-error results in other populations. Nor should every movement-related signal be re-labelled as error without an independently specified model that predicts it.^[ref]

The measurement compartment is another limitation. A comparison of identified VTA cell spiking with nucleus-accumbens dopamine release in a common decision task found release dynamics related to evolving reward expectations without corresponding changes in recorded cell-body spiking. Local release regulation therefore matters: a terminal sensor and a somatic spike train need not report the same quantity. This is especially important for interpreting slower ramps or motivational effects as direct readings of a broadcast phasic error.^[ref]

Further qualifications concern learning itself. An error depends on subjective reward, inferred state, uncertainty and task structure; it is not calculated from the experimenter’s reward schedule alone. Values informed by a model of the environment can also contribute, so identifying an error-like signal does not automatically prove exclusively model-free learning. Conversely, a very flexible state model can accommodate almost any observation unless its representation and parameters are constrained independently. Good tests manipulate expectation while matching outcome and movement, and compare explicit alternatives on held-out behaviour and neural data.

The account is thus still appropriate when stated at the right resolution: some temporally precise dopamine activity supplies expectation-dependent teaching signals, supported by both recordings and causal learning manipulations. A complete account of the dopamine system must additionally explain projection-specific coding, local transmitter dynamics, motivation and movement. The productive revision is a family of constrained computational and circuit hypotheses, rather than either one molecule–one number or abandonment of prediction error altogether.

Post-receptoral processing compares cone signals, organizes those comparisons across space and time, and uses them to represent surfaces and objects. Three receptor classes do not imply three neurons labelled with three color names. A complete answer must distinguish the information available at the receptors from its neural transformation and its perceptual interpretation.

The principle of univariance supplies the starting constraint. For one cone class, an absorption-dependent response does not identify wavelength independently of intensity. Comparing differently sensitive classes makes chromatic information available. Relative to an adapted background, cone contrasts can be written $c_j=(Q_j-Q_{j0})/Q_{j0}$, where $Q_{j0}>0$ is the background catch. This representation separates a change relative to the prevailing input from the raw number of absorbed photons. It is a useful coordinate system, not a claim that every neuron implements this exact formula.

Retinal circuitry immediately transforms the signals. Light hyperpolarizes cone cells and changes their transmitter release. ON and OFF bipolar pathways respond differently to changes in cone output, while horizontal and amacrine circuitry contribute spatial and temporal interactions. Retinal ganglion cells transmit spike trains rather than simply forwarding a cone's voltage. Their center–surround organization compares regions, so chromatic and spatial structure are already intertwined. Near the central retina, the midget pathway can oppose an L- or M-dominated center to a surround containing other cone input. The same circuitry can carry information useful for fine spatial structure and for chromatic contrast; these are not necessarily separate wires.

A distinct S-cone route illustrates the anatomical specificity. Dennis M. Dacey and colleagues combined physiological characterization with identification of ganglion-cell morphology in primate retina, finding that the S-ON opponent response belonged to a small bistratified ganglion-cell type rather than simply another midget cell.^[ref] This supports parallel circuits with different cone combinations. Calling the pathway “blue-on” is shorthand for its stimulus response, not proof that the cell itself produces a unique blue experience.

At a simplified computational level, an achromatic-like combination and two chromatic combinations can be represented asThis invertible three-channel opponent transform reorganizes three coordinates rather than creating additional independent spectral information. Their matrix has determinant $-2$, so the transform is invertible: $c_L=(y_0+y_1)/2$, $c_M=(y_0-y_1)/2$ and $c_S=y_2+y_0/2$. Actual weights depend on cell, location, adaptation and stimulus conditions. This algebra is an illustrative opponent transform, not a fitted retinal wiring diagram or a complete perceptual color space.

The lateral geniculate nucleus preserves several parallel retinal pathways while also participating in cortical interactions. A. M. Derrington and colleagues studied macaque geniculate cells with gratings whose luminance and chromatic modulation could be varied independently on a calibrated color television display. Finding directions that changed cone stimulation without driving a cell constrained its linear cone weights.^[ref] For a nonzero linear weight vector, the zero-response contrasts obey $r={\mathbf{w}}^{\mathsf{T}}\mathbf{c}=0$, defining a chromatic null plane. A null response therefore tests a relation between inputs, rather than merely assigning the cell a preferred color name. The resulting opponent directions need not coincide with human unique hues.

Cortical processing adds more complex spatial comparisons. Bevil R. Conway mapped cone inputs to color-responsive neurons in alert macaque V1 using cone-isolating spots and eye-position-corrected reverse correlation. Some cells combined chromatic opponency with opposite tuning in adjacent spatial subregions.^[ref] Such double-opponent cells can respond strongly to a chromatic boundary while responding less to a spatially uniform change. This offers a mechanism useful for separating surface differences from shared illumination changes. It does not establish that one V1 population alone solves color constancy under arbitrary scenes.

Further processing must relate color to form, illumination, material and the current task. The same cone contrasts can be interpreted differently in a shadow, highlight or translucent object. Surface color judgments, discrimination thresholds and color naming therefore constrain different computations. Feedback and attention can influence which information is selected, while learning can alter how a signal is interpreted or normalized. A lesion or imaging response can constrain contributions to those tasks, but an area responsive to color is not automatically an isolated “color center.”

The strongest picture is consequently a hierarchy of interacting computations: cone-specific signals, opponent and spatial filtering, combinations tuned to richer image structure, and object-level judgments. Trichromatic sampling limits spectral information; later processing determines how that information supports boundaries, stable surface appearance and behavior.

The primary visual cortex, often abbreviated V1, is an early cortical stage of visual processing. Its neurons have receptive fields: changing a stimulus in a particular part of the visible scene can change their spike trains. Orientation selectivity adds a further distinction: at the same location, a properly aligned bar can evoke a strong response while a rotated bar evokes little. Simple cells make the dependence on local position and visual contrast especially clear. See the original experiments under orientation selectivity and simple cell.

In their 1962 study, David Hubel and Torsten Wiesel recorded from forty cats under thiopental anesthesia, maintained light anesthesia while monitoring cortical electrical activity, and used succinylcholine to stabilize the eyes. White-light patterns were projected onto a screen $1.5\mathrm{m}$ away. Extracellular single-unit recording with tungsten electrodes sampled individual units; receptive fields were mapped separately for each eye. This connected controlled stimulus geometry to individual spike trains. ^[ref]

These acute, anesthetized-cat recordings do not directly measure human perception or responses during awake behavior. The useful idea is a local description of structure: a population can distinguish where a feature lies, which way it is aligned, and which side is brighter. The signed spatial filter model of a simple cell turns that idea into a computable example. One neuron's preferred stimulus does not establish recognition of an entire object. Likewise, orientation selectivity in one neuron and an orderly map of preferences across primary visual cortex are different findings: rodents can have strong orientation selectivity without the orderly arrangement typical of some other species. ^[ref] The existing correlation-based development of visual cortical maps and orientation-map self-organization sections describe hypotheses about how organization develops.

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.

A signed spatial filter model of a simple cell assigns a weight $K_{ij}$ to each sampled location and takes its dot product with the visual contrast image $C_{ij}$. Positive weights favor bright features; negative weights favor dark features. The signed driveis a linear function of the image, while the rectified linear unit output $R$ is nonlinear. $D$ and $R$ below are arbitrary model units. To obtain a predicted firing rate one would need a calibrated gain, background response and possibly a different output function. The model does not generate a spike train; a linear-nonlinear-Poisson cascade model adds a stochastic event-generation assumption.

For the first diagram, every one of five rows has weightsA bright central column, with $C_{i3}=1$ and zero elsewhere, gives $D=5/5=1$. A bright middle row gives $D=-1/20-1/20+1/5-1/20-1/20=0$; a uniformly bright patch also gives $D=0$, since every row sums to zero. Move the bright column one site right: $D=5(-1/20)=-1/4$, so $R=0$. Make that shifted column dark instead: $D=1/4$ and $R=1/4$. Thus an OFF-like weight can increase output for darkness. A negative weight is a descriptive sensitivity, not evidence for a particular inhibitory synapse.

For a smooth example use the real, even Gabor filterThe Gaussian function localizes the weights; the cosine supplies alternating polarity. Choose ${\sigma }_x=1/2$, ${\sigma }_y=1$ and $q=\pi$ in arbitrary spatial units. The stripes' spatial frequency is $q/(2\pi )=1/2$ cycles per unit. Their axis makes angle $\theta$ with horizontal, while brightness changes across the perpendicular axis. $\phi$ is their spatial phase. Normalize the full-plane integral of $K{C}_{\theta ,\phi }$ by its value at $\theta =\pi /2$, $\phi =0$. The Gaussian integral and the product identity for cosines givewhereAt $\phi =0$, horizontal stripes give $R\approx 0.004159$, $60^{\circ }$ stripes give $R\approx 0.286741$, and vertical stripes give $R=1$. At the same vertical orientation, $\phi =\pi /2$ gives $D=R=0$, while reversing the polarity with $\phi =\pi$ gives $D=-1$ and $R=0$. The preferred axis persists while the response changes dramatically with position. Shifting the vertical stripes by a quarter of their wavelength changes phase by $\pi /2$ up to sign, without rotating them.

These are two different toy filters, not discretizations of the same fitted receptive field. The discrete filter exactly rejects a uniform contrast patch; this Gabor filter need not reject uniform contrast because its weight integral need not vanish. Neither toy includes stimulus history, recurrent interactions, saturation or a measured response baseline. Their value is a transparent prediction: orientation selectivity can emerge from spatial summation, and it does not by itself imply tolerance to phase.

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.