Pulse-Coupled Adaptive Winfree Network

• Pulse-coupled adaptive Winfree networks are phase-oscillator systems that combine pulse-based interactions with adaptive mechanisms like synaptic depression and Hebbian updates.
• They model diverse regimes including balanced neural dynamics, distributed clock synchronization, and chimera states using event-triggered pulses and phase-response curves.
• Mathematical analyses such as mean-field methods and bifurcation studies reveal how adaptive rules and pulse timing shape collective oscillations and synchronization.

A pulse-coupled adaptive Winfree network is a class of phase-oscillator systems in which the state of each oscillator evolves according to a Winfree-type interaction, namely a phase-response curve multiplied by a pulse-derived stimulus, while at least one coupling component adapts on a slower or event-triggered timescale. In the recent literature, the term covers several closely related constructions: mean-field excitatory–inhibitory phase networks with short-term synaptic depression, globally coupled pulse systems with Hebbian weight adaptation, graph-based hybrid pulse networks with adaptive refractory or gain variables, and delay-bearing excitable networks whose topology self-adjusts to achieve frequency synchronization (Mato et al., 2024). Across these variants, the defining ingredients are pulse timing, a phase response curve, and adaptive coupling or adaptive internal state variables that reshape the effective interaction structure.

1. Model class and formal definition

In the most direct Winfree formulation, oscillator phases evolve as

$\dot{\theta}_i=\omega_i+Z(\theta_i)\,S(t),$

where $Z$ is the phase response curve and $S$ is a stimulus built from population activity. The pulse-coupled adaptive variants retain this factorization but replace smooth mean fields by event-driven pulses, and replace static coupling by adaptive weights, synaptic efficacies, or auxiliary throttling variables (Mato et al., 2024).

A canonical two-population formulation is the excitatory–inhibitory model

$$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$

with reset to $0$ after threshold crossing at $\phi=1$, pulse fields

$$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$

and a Type-I polynomial phase-response curve

$Z(\phi)=16\,\phi^2(1-\phi)^2.$

The aggregate currents are

$$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$

with

$$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$

This setup is explicitly described as Winfree because the stimulus enters as a common mean-field pulse signal multiplied by the PRC, with distinct receiver-specific stimuli $Z$0 and $Z$1 (Mato et al., 2024).

Adaptive structure can take several non-equivalent forms. In the excitatory–inhibitory model, adaptation appears as short-term depression on excitatory-to-excitatory synapses:

$Z$2

so that each presynaptic spike reduces efficacy by factor $Z$3 and recovery occurs on timescale $Z$4 (Mato et al., 2024). In a distinct globally coupled adaptive pulse network, the coupling weights themselves evolve according to

$Z$5

with the phase dynamics

$Z$6

where $Z$7 is a frustration parameter entering the PRC and $Z$8 for the case emphasized in that study (Anand et al., 11 Mar 2026). A third adaptive construction uses leak plus spike-triggered plasticity,

$Z$9

with $S$0 and $S$1, producing an adaptive all-to-all directed pulse network with identical native frequencies (Kasatkin et al., 2018).

These formulations share the same structural motif: pulse timing is the carrier of interaction, whereas adaptation determines how strongly pulses matter.

2. Adaptive mechanisms

The adaptive variable in a pulse-coupled adaptive Winfree network need not be a weight matrix. The literature uses at least four distinct mechanisms.

First, short-term depression implements synaptic resource depletion. In the excitatory–inhibitory mean-field model, only excitatory-to-excitatory synapses are adaptive, with $S$2 and $S$3. Between spikes,

$S$4

and the efficacy at spike emission is

$S$5

This mechanism is presented as compensation for energy dissipated by pulse emission and as a resource-limiting nonlinearity essential for the balanced state (Mato et al., 2024).

Second, Hebbian adaptation makes coupling strengths coevolve with phase relations. In the frustration-mediated globally coupled model,

$S$6

so weights relax toward $S$7. In fully entrained states this gives $S$8, while in antipodal arrangements it yields approximately $S$9 across clusters and approximately $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$0 within clusters. No separate hard bounds or normalization constraints are imposed; the target confines effective weights to $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$1 because $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$2 (Anand et al., 11 Mar 2026).

Third, event-triggered plasticity can be expressed as a leak plus pulse-induced increment. In the itinerant-chimera model, between pulses the weights decay exponentially, and at each spike of $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$3 one has

$$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$4

This is described as a phase-dependent Hebbian/STDP-like update in a phase-oscillator framework, with $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$5 corresponding to an STDP-like rule and $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$6 to a qualitatively Hebbian-like rule (Kasatkin et al., 2018).

Fourth, adaptation may be encoded in auxiliary node states rather than in synaptic weights. The adaptive 4-coupling on trees augments each node with

$$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$7

where $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$8 is the rested/refractory state and $$\dot{\phi}^{e/i}_n(t) = \omega^{e/i}_n + G\,C^{e/i}(t)\,Z\!\left(\phi^{e/i}_n(t)\right),\qquad 1\le n\le N,$$9 implements a pull counter. Excitation is defined by

$0$0

and the refractory variable evolves according to

$0$1

Pulses are ignored in refractory states, which the paper describes as throttling the input (Lyu, 2016).

A plausible implication is that “adaptation” in this literature is best understood functionally rather than structurally: it denotes any state-dependent mechanism that feeds past pulse activity back into future coupling efficacy.

3. Mean-field structure, asynchronous states, and balance

For the excitatory–inhibitory pulse-coupled adaptive Winfree network, the asynchronous regime is analyzed by a self-consistent mean-field method rather than by an Ott–Antonsen reduction. Writing

$0$2

and assuming constant $0$3, the interspike interval is

$0$4

The population-averaged fields satisfy

$0$5

together with the self-consistency conditions

$0$6

The frequency heterogeneity is introduced through compact-support distributions $0$7 with

$0$8

corresponding to means approximately $0$9 for excitatory and approximately $\phi=1$0 for inhibitory units (Mato et al., 2024).

The same work identifies a strong-coupling balanced asynchronous solution. In the limit $\phi=1$1, bounded fields require

$\phi=1$2

Below a threshold $\phi=1$3, the asynchronous regime has constant fields, the exact mean-field self-consistent solution matches simulations, and the fields remain finite and balanced (Mato et al., 2024).

The hybrid synchronization literature addresses a different balance question: under what conditions pulse-coupled Winfree-type networks converge to exact synchrony? For the rooted-graph hybrid model with identical frequencies, the phase transition curves $\phi=1$4 are required to satisfy the delay-advance property, and the diameter $\phi=1$5 of the shortest arc containing all phases is non-increasing when $\phi=1$6. Rootedness is sufficient, and also necessary, for synchronization in that setting (Proskurnikov et al., 2015).

These two strands are mathematically distinct. One studies mean-field balance with heterogeneity, excitation, inhibition, and synaptic depression; the other studies contraction of phase diameter under hybrid phase-jump maps. This suggests that “balance” in pulse-coupled adaptive Winfree networks may refer either to bounded collective fields or to contraction of phase spread, depending on the model class.

4. Collective regimes and bifurcation structure

The most detailed bifurcation picture currently available for the two-population excitatory–inhibitory model is organized by a Hopf destabilization of the asynchronous regime. Large $\phi=1$7 beyond a threshold $\phi=1$8 destabilizes the fluctuationless asynchronous state, with reported values

$\phi=1$9

Below $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$0 the system remains asynchronous; at $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$1 it undergoes a Hopf bifurcation; above $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$2 it enters an irregular oscillatory regime described as collective chaos (Mato et al., 2024).

The evidence for that regime is macroscopic and microscopic. The fields $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$3 and $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$4 exhibit large, irregular oscillations with broadband power spectra that are essentially independent of $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$5 for $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$6. Microscopic activity is strongly irregular, with inhibitory coefficients of variation approximately $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$7–$$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$8 and excitatory coefficients of variation approximately $$E^{e/i}(t)=\frac{1}{N}\sum_{k,m\,:\,t^{e}(k,m)<t}x^{e/i}_k(t)\,\delta\!\big(t-t^{e}(k,m)\big),\qquad I(t)=\frac{1}{N}\sum_{k,m\,:\,t^{i}(k,m)<t}\delta\!\big(t-t^{i}(k,m)\big),$$9–$Z(\phi)=16\,\phi^2(1-\phi)^2.$0, depending on quenched or annealed disorder. The largest Lyapunov exponent is not reported; the evidence consists instead of broadband spectra, irregular macroscopic fields, and irregular spiking statistics (Mato et al., 2024).

A different adaptive pulse-coupled Winfree network, globally coupled and frustration-mediated, exhibits a broader catalog of collective states. For $Z(\phi)=16\,\phi^2(1-\phi)^2.$1, $Z(\phi)=16\,\phi^2(1-\phi)^2.$2, $Z(\phi)=16\,\phi^2(1-\phi)^2.$3, $Z(\phi)=16\,\phi^2(1-\phi)^2.$4, $Z(\phi)=16\,\phi^2(1-\phi)^2.$5, and $Z(\phi)=16\,\phi^2(1-\phi)^2.$6, the reported regimes are frequency-clustered states, entrainment, bump states, bump–frequency cluster states, antipodal and multi-antipodal cluster states, chimera states, and incoherent dynamics. The one- and two-parameter diagrams are organized in the $Z(\phi)=16\,\phi^2(1-\phi)^2.$7 plane. For $Z(\phi)=16\,\phi^2(1-\phi)^2.$8 the sequence as $Z(\phi)=16\,\phi^2(1-\phi)^2.$9 increases is FC $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$0 ENT $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$1 BFC $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$2 BS; for $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$3 it is AP $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$4 MAC $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$5 CHI; for $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$6 it is FC $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$7 CHI $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$8 INC (Anand et al., 11 Mar 2026).

Adaptive pulse networks on trees show yet another regime structure. The inhibitory 4-coupling synchronizes arbitrary initial data on finite trees with maximum degree $$C^{e}(t)\equiv g^{e}_e\,E^{e}(t)-g^{e}_i\,I(t),\qquad C^{i}(t)\equiv g^{i}_e\,E^{i}(t)-g^{i}_i\,I(t),$$9 by time $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$0, while the adaptive 4-coupling synchronizes arbitrary initial joint configurations on arbitrary trees by time $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$1. The non-adaptive rule admits non-synchronizing examples for $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$2, whereas the adaptive rule overcomes that obstruction through refractory throttling (Lyu, 2016).

Collectively, these results indicate that pulse-coupled adaptive Winfree networks do not possess a single universal phenomenology. Depending on whether adaptation acts as depression, Hebbian learning, or refractory throttling, the same Winfree-type coupling can support asynchronous balance, irregular collective oscillations, exact synchronization on trees, or frustration-induced clustered and chimera-like states.

5. Mechanisms of homeostasis, synchronization, and irregularity

In the two-population excitatory–inhibitory model, the central mechanism is a homeostatic interplay between PRC shape and adaptation. As synchrony increases with $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$3, pulses arrive at phases close to threshold or reset, where $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$4, so the effective drive $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$5 remains bounded even though $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$6 diverges. Quantitatively, conditioned order parameters at spike times satisfy

$$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$7

so the average PRC sampled at pulse arrivals scales as $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$8. Unconditioned order parameters remain strictly below $$g_e^e=1,\quad g_e^i=1,\quad g_i^e=\tfrac{1}{2},\quad g_i^i=2.$$9, and the PDF of $Z$00 develops an integrable power-law singularity at $Z$01 with exponent approximately $Z$02 (Mato et al., 2024).

Short-term depression supplies the second part of that mechanism. Each spike reduces the excitatory efficacy by factor $Z$03 and recovery occurs on timescale $Z$04; with the reported parameters $Z$05 and $Z$06, this acts as a resource-limiting nonlinearity. The paper explicitly states that STD is essential to support a balanced asynchronous regime and for the homeostatic reduction of effective excitatory drive (Mato et al., 2024).

In the frustration-mediated Hebbian model, the mechanism is different. The effective input is

$Z$07

so the phase lag does not shift a pairwise term $Z$08, but instead shifts the postsynaptic PRC. The reported interpretation is that varying $Z$09 moves the advance/delay window relative to the pulse times and postsynaptic phase, while the Hebbian rule strengthens couplings between near-in-phase oscillators and weakens or inverts couplings between antiphase oscillators. The paper attributes spontaneous entrainment, bump, and bump–frequency cluster states to this interplay of pulse coupling, a Type-II-like PRC with frustration, and Hebbian plasticity (Anand et al., 11 Mar 2026).

The synchronization proofs for rooted hybrid pulse networks rely on still another mechanism: the delay-advance property contracts phase diameter. When $Z$10 for $Z$11 and $Z$12 for $Z$13, pulses move leading oscillators backward and lagging oscillators forward without overshoot. Rootedness guarantees that this contraction propagates through the graph within finitely many firing rounds (Proskurnikov et al., 2015).

These mechanisms are not interchangeable. The excitatory–inhibitory model obtains strong yet imperfect synchronization together with irregularity; the rooted hybrid model proves exact asymptotic synchronization under identical frequencies; the frustration-mediated model stabilizes multiple clustered and chimera-like states. The shared Winfree structure therefore does not determine the collective state by itself; the adaptive law and the PRC geometry are decisive.

6. Numerical methods, observables, and relation to adjacent model classes

The event structure of these networks strongly shapes their numerical treatment. The excitatory–inhibitory model is simulated by event-driven integration, exploiting exact evolution between spikes and nonlinear updates at spike times; the reported network sizes are $Z$14, $Z$15, and $Z$16, equally split between excitatory and inhibitory populations. Quenched disorder fixes $Z$17 per neuron, whereas annealed disorder redraws $Z$18 values from $Z$19 every $Z$20 spikes of the network. Initial transients of $Z$21 time units are discarded, and spectra are computed over $Z$22. For plotting only, pulse fields are filtered by the exponential kernel $Z$23 with $Z$24, while the dynamics themselves use $Z$25-pulses (Mato et al., 2024).

The principal macroscopic observables in that setting are the order parameters

$Z$26

population firing rates $Z$27, and the aggregate pulse fields entering the currents. In the frustration-mediated Hebbian model, regime classification is instead based on three incoherence measures: a frequency-based strength of incoherence $Z$28, an instantaneous-phase-based strength $Z$29, and a mean-frequency-per-bin strength $Z$30, using $Z$31 bins, $Z$32 oscillators per bin, and thresholds $Z$33, $Z$34, and $Z$35 (Anand et al., 11 Mar 2026).

Relative to adjacent model classes, pulse-coupled adaptive Winfree networks are distinguished by PRC-times-stimulus coupling. Classical Winfree coupling takes the form

$Z$36

whereas Kuramoto-type systems use phase-difference coupling,

$Z$37

The excitatory–inhibitory model explicitly contrasts these forms, noting that Kuramoto-type models do not include pulse timing or short-term depression (Mato et al., 2024). The broader Winfree literature also shows that non-adaptive pulse shape and PRC offset already organize distinct synchronization scenarios in the mean-field limit, including Bogdanov–Takens and mutated BT′ structures, but those studies treat non-adaptive baselines rather than adaptive pulse networks (Gallego et al., 2017).

Several limitations recur across the literature. The excitatory–inhibitory chaotic-synchronization model uses mean-field coupling and single-variable phase reductions; the rooted-hybrid theory assumes identical frequencies and an initial diameter smaller than $Z$38; the tree-synchronization results are rigorous for trees rather than arbitrary graphs; the frustration-mediated Hebbian study assumes identical oscillators, global coupling, no synaptic delays, no noise, and a specific PRC and pulse form (Mato et al., 2024). A plausible implication is that the term “pulse-coupled adaptive Winfree network” currently denotes a family of mathematically related but not yet unified models.

7. Biological relevance and broader significance

The biological interpretation is most explicit in the two-population excitatory–inhibitory model. Its architecture is described as mimicking a neural network composed of excitatory and inhibitory neurons, and the short-term depression is said to model finite synaptic resources and compensate for energy dissipated by pulse emission. The irregular single-neuron coefficients of variation, especially inhibitory values approaching approximately $Z$39, are presented as consistent with cortical spike-train statistics (Mato et al., 2024).

Adaptive pulse-coupled tree models and rooted-graph hybrid models emphasize a different application domain: distributed clock synchronization. In those works, pulses are short messages, PRCs are the clock update laws, and the hybrid event structure directly models threshold-reset clocks. The rootedness result establishes synchronization under minimal connectivity assumptions, while the adaptive 4-coupling yields a universal randomized distributed clock synchronization algorithm with $Z$40 memory per node and expected worst case running time

$Z$41

built from distance-$Z$42 coloring, a randomized spanning tree construction, and A4C/M synchronization on the resulting tree (Proskurnikov et al., 2015).

The adaptive delay-bearing excitable network extends the biological analogy further by using excitatory and inhibitory nodes, distance-related delays, and a local frequency-error adaptation rule on synaptic weights. That work reports sparse, anti-cluster, delay-structured solutions, a necessary minimum inhibitory fraction approximately $Z$43–$Z$44, and stronger long-range inhibitory projections than excitatory ones after adaptation (Gil, 2023).

Taken together, these developments establish pulse-coupled adaptive Winfree networks as a common language for several research programs: balanced neural population dynamics, exact synchronization in hybrid oscillator networks, adaptive clustered and chimera-like states, and pulse-based distributed coordination. What unifies them is not a single normal form, but the combination of three principles: pulse-mediated interaction, PRC-governed susceptibility, and adaptation driven by spiking or relative phase.