Response to sub-threshold stimulus is enhanced by spatially heterogeneous activity

Sub-threshold stimuli cannot initiate excitations in active media, but surprisingly as we show in this paper, they can alter the time-evolution of spatially heterogeneous activity by modifying the recovery dynamics. This results in significant reduction of waveback velocity which may lead to spatial coherence, terminating all activity in the medium including spatiotemporal chaos. We analytically derive model-independent conditions for which such behavior can be observed.

💡 Research Summary

The paper investigates how sub‑threshold electrical stimulation—currents too weak to elicit an action potential in a quiescent excitable medium—can nevertheless dramatically alter the dynamics of a medium that already contains spatially heterogeneous activity such as spiral waves or spatiotemporal chaos. Using detailed ionic models (Luo‑Rudy I and the TNNP model) the authors simulate one‑dimensional cables and two‑dimensional sheets of excitable cells, applying a globally uniform sub‑threshold current I for a duration τ. They find that while the wave‑front velocity (c_f) remains essentially unchanged, the wave‑back velocity (c_b) is strongly reduced during stimulation. This reduction lengthens the refractory (recovery) period of each cell, enlarges the inactive region behind a propagating wave, and consequently narrows the excitable gap between successive waves. When the widened refractory zone exceeds the spacing λ between two successive excitation fronts, the recovery front of the leading wave collides with the excitation front of the following wave, causing both to terminate.

The authors formalize this mechanism by defining an average wave‑back speed ˜c_b(I) that depends only on the stimulus strength. They derive a simple analytical condition for the minimal stimulus duration required to extinguish a source of recurrent activity:

τ_min = (λ – c_f·τ_r) / (c_f – ˜c_b(I))

where τ_r is the intrinsic excited‑state duration and λ is the distance between successive wave fronts in the unperturbed medium. This “strength‑duration” relationship matches numerical results for both single spirals and chaotic states across a wide range of I and τ values. A lower bound I_low exists below which c_b does not decrease enough to affect the dynamics, regardless of τ.

Physiologically, the sub‑threshold current primarily reduces the conductance of outward potassium channels, thereby suppressing the slow outward K⁺ current that normally drives repolarization. This effect prolongs the refractory period without altering the depolarizing mechanisms that set c_f. Consequently, the phenomenon cannot be captured by simplistic excitable‑media models that ignore stimulus‑induced changes in recovery dynamics.

From an application perspective, current clinical strategies for terminating arrhythmias rely on supra‑threshold, high‑energy shocks that simultaneously depolarize the entire tissue. Such approaches are energetically costly and can cause tissue damage. In contrast, the sub‑threshold method described here requires only microampere‑per‑square‑centimeter currents applied globally, yet can rapidly synchronize and extinguish pathological wave activity. This suggests a low‑energy, minimally invasive alternative for controlling spatiotemporal chaos in biological excitable media, including cardiac fibrillation.

In summary, the study reveals that sub‑threshold stimulation, though ineffective in a resting medium, becomes a powerful control tool in the presence of heterogeneous activity by selectively lengthening the recovery phase, reducing wave‑back speed, and collapsing the excitable gap. The derived analytical condition provides a practical guideline for designing stimulus protocols, and the findings open new avenues for low‑energy defibrillation and for the broader control of excitable systems in chemistry, neuroscience, and engineering.