Research analysis · Spiking network training

Spikes that never vanish mid-trial: a trainability dividend

Wenig, Memmesheimer and Klos ran a controlled head-to-head between the two most common spiking neuron models on a standard spoken-digit benchmark, with exact spike-based gradient descent and a thorough hyperparameter search for each. Quadratic integrate-and-fire (QIF) networks reached 90.1 percent test accuracy; leaky integrate-and-fire (LIF) networks reached 79.2 percent. Landscape imaging then attributes the gap to a specific mechanism: in LIF, arbitrarily small weight changes can make a spike appear or disappear mid-trial, fragmenting the loss surface, while in QIF every spike time moves continuously.

Source: Quadratic integrate-and-fire neurons exhibit less fragmented loss landscapes and outperform leaky integrate-and-fire neurons in spike-based gradient descent, arXiv:2606.03935, 2 June 2026. Primary source. Read the full arXiv HTML version, including the hyperparameter search protocol, landscape analyses, and appendices.

What the work claims

This is a primary empirical comparison. The authors train two-layer feedforward spiking networks, 128 hidden neurons and 20 output neurons, on the Spiking Heidelberg Digits (SHD) dataset using time-to-first-spike coding and exact spike-based gradient descent, where spike times are treated as differentiable functions of the weights. Both models receive the same care: a rough grid search, a finer search around the best region, then retraining of the selected configurations for 500 epochs with 10 seeds and batches of 1,000 samples.1

The best QIF network reaches 90.1 plus or minus 9 percent test accuracy against 79.2 plus or minus 13 percent for the best LIF network. Neither is close to the roughly 96 percent state of the art on SHD, and the authors say so; the architectures are deliberately simple because the question is the model comparison, not the leaderboard. In exploratory larger networks the ordering holds: 91.3 plus or minus 8 percent (QIF) versus 85.2 plus or minus 13 percent (LIF).1

The second claim is mechanistic. Loss and gradient landscapes around the training trajectories are fragmented for LIF networks, with discontinuities that coincide with changes in the temporal order of spikes, and the batch-averaged gradients are erratic. QIF landscapes show sharp but continuous loss boundaries and far less fragmentation. Consistently, the gradient is a better predictor of what training actually does: for QIF, 90 percent of steps where a loss decrease is expected from the gradient in fact produce one, versus 75 percent for the best LIF network.

How it works

The LIF neuron has a finite firing threshold. An input spike can shove the membrane potential across that threshold with essentially no warning, so an infinitesimal change in a weight can create or annihilate a spike at an arbitrary time. When that happens, every downstream spike time rearranges, the loss jumps, and the local gradient computed just before the jump says nothing about what lies across it. The authors show that loss-landscape fragment boundaries line up with such spike-order changes, and that hidden-layer spikes abruptly disappear at one time and reappear earlier at the boundary.1

The QIF neuron has an effectively infinite threshold: its membrane potential diverges as the spike approaches, so no input spike can force an immediate crossing, and spike times deform continuously as parameters change. New spikes can still be added through a construction the authors call pseudodynamics, an extended tail of the trial in which neurons spike as often as the task requires, with pseudospike times computed in a single feedforward pass. The one discontinuity that remains in QIF is confined to gradient events where input and output spike orderings swap, and even there the loss itself stays continuous.1

A second thread concerns robustness of the code itself. Trained LIF networks tend to emit nearly synchronous output spikes: the best LIF model has an average gap of only 3 milliseconds between its first two output spikes, against 28 milliseconds for the best QIF model. Synchronous spikes are fragile under the timing jitter that real neuromorphic hardware injects. When the authors selected the one LIF configuration with comparably distinct spikes (29 milliseconds), its accuracy collapsed to 50.8 plus or minus 40 percent, and 62 percent of its gradient components vanished entirely, against about 25 percent for the well-behaved models. So the LIF shortcut, clustering spikes to avoid landscape cliffs, buys jitter intolerance at a hidden cost.

Where a skeptic should push

The most load-bearing assumption is that a result on one dataset, one coding scheme, and one training method generalizes to the way spiking models are actually used. Everything here is time-to-first-spike coding on SHD with exact spike-based gradient descent; the neurons are intrinsically oscillating and the input currents act as voltage jumps rather than the more common exponentially decaying currents. The authors state these restrictions plainly. Whether QIF keeps its advantage under surrogate-gradient training, rate coding, or recurrent architectures is explicitly left open.1

Sample sizes deserve weight too. Ten seeds with standard deviations of 9 to 13 percentage points means the 90.1 versus 79.2 headline is a distribution-level claim resting on modest replication, though the paired protocol (same seeds for both models) and the consistency of the landscape diagnostics make the ordering plausible. Prior work has trained LIF networks successfully on SHD-class tasks with modified losses, careful initialization, and ad-hoc anti-silencing tricks; this paper's counter is that those tricks are evidence of the disease, not refutation of it, but a reader should hold both facts.

Finally, landscape imaging is suggestive, not causal proof. Continuous dynamics guarantee continuity of the loss; they do not by themselves guarantee a well-conditioned or easily optimizable surface, and the paper does not rule out that other continuous models could fragment in different ways. The claim that survives scrutiny is narrower and still useful: discontinuities of LIF type are absent in QIF, and measured trainability differences line up with that absence.

Picking the neuron model that stands in for tissue

Organoid intelligence runs on tissue, but it models tissue in software constantly: digital twins for closed-loop stimulation, surrogate networks for testing decoding strategies before electrodes go in, and trainable stand-ins whenever backpropagation through a biological loop is approximated rather than executed. Every one of those stand-ins picks a neuron model, and the field's default is LIF, because LIF is cheap and everywhere. This paper says that default has a hidden tax: the stand-in's loss surface is fractured in ways the real tissue is not, so training experiments fail for reasons that belong to the simulator, not the biology.1

The non-obvious implication cuts toward fidelity as well as convenience. Cortical neurons commonly show type-I excitability, meaning their firing rate grows continuously from zero as drive increases, the same continuous-onset regime the QIF model encodes, whereas the classic LIF with a fixed threshold is closer to a type-II caricature. A QIF-based digital twin may therefore be both easier to train and a more honest representation of the excitable membrane it stands in for. For a field that wants to run gradient-based optimization against organoid recordings, the twin's trainability is not an implementation detail; it decides whether the optimization pipeline converges at all.

The threat is subtler than a wrong number. If surrogate LIF models of organoid dynamics are systematically harder to train, the failures will be misread: teams will conclude the readout problem is hard, or the tissue is uncooperative, when part of the difficulty was discontinuities baked into the surrogate. That misdiagnosis has a cost in abandoned directions and over-built hardware. The opportunity, equally concrete, is a cheap experiment: repeat the paper's protocol against organoid MEA recordings instead of SHD, and ask whether QIF-family surrogates reach lower decoding error with less tuning. If they do, the neuron-model choice becomes a reportable experimental variable in every organoid computing paper, the way optimizer choice is in deep learning.

The bottom line

Established: under a matched protocol on SHD with exact spike-based gradient descent, QIF networks outperform LIF networks (90.1 versus 79.2 percent test accuracy), their loss and gradient landscapes are measurably less fragmented, their gradients predict actual loss changes more often (90 versus 75 percent of steps), and their output codes are more temporally distinct. Established: the LIF discontinuity mechanism is spike (dis)appearance at arbitrary times, and forcing LIF toward distinct spike timings destroys its accuracy. Not established: transfer to surrogate-gradient training, rate coding, recurrent networks, or any task beyond this benchmark and its variants.

For organoid intelligence, treat this as hygiene for the digital half of every wetware-digital loop. What would confirm relevance: surrogate-model studies of organoid recordings showing that continuous-onset models train more reliably than threshold models on the same data. What would break it: evidence that surrogate neuron-model choice stops mattering once models are scaled or once training moves off exact spike-based gradients. Either outcome is worth knowing before another training pipeline is blamed on the biology.

Frequently asked questions

What are LIF and QIF neurons?

Both are single-neuron models used in spiking networks. Leaky integrate-and-fire sums inputs with decay and emits a spike when a finite threshold is crossed. Quadratic integrate-and-fire uses a quadratic voltage equation whose potential diverges at the spike, giving continuous spike-time changes instead of abrupt threshold crossings.

What is exact spike-based gradient descent?

A training method that treats spike times as differentiable functions of the network parameters and backpropagates through them, rather than approximating the spike nonlinearity with a smooth surrogate function as surrogate-gradient methods do.

Why do LIF loss landscapes fragment?

Because an arbitrarily small parameter change can make a spike appear or vanish mid-trial, reshuffling all downstream spike times. The boundaries between landscape fragments coincide with these spike-order changes, and the loss jumps discontinuously across them.

How big is the measured accuracy gap?

On the Spiking Heidelberg Digits with time-to-first-spike coding, the best QIF network reached 90.1 plus or minus 9 percent test accuracy versus 79.2 plus or minus 13 percent for the best LIF network, with 10 seeds each. In larger exploratory networks the figures were 91.3 versus 85.2 percent.

Does QIF fix all discontinuities?

No. Gradient discontinuities remain where the temporal order of input and output spikes swaps, but the loss stays continuous there, new spikes enter only at trial end through pseudodynamics, and measured fragmentation is far lower than for LIF.

Why should organoid researchers care?

Every digital twin or surrogate of living tissue picks a neuron model, and LIF is the common default. If the default's discontinuities make training fail, those failures get blamed on the biology. Continuous-onset models like QIF are both more trainable and arguably closer to cortical type-I excitability, so the choice should be a reported experimental variable.

References

  1. Wenig C, Memmesheimer RM, Klos C. Quadratic integrate-and-fire neurons exhibit less fragmented loss landscapes and outperform leaky integrate-and-fire neurons in spike-based gradient descent. arXiv:2606.03935 [cs.NE]. 2026. https://arxiv.org/abs/2606.03935. Accessed 2026-09-16.