pith. machine review for the scientific record. sign in

arxiv: 2604.22113 · v1 · submitted 2026-04-23 · 🌊 nlin.PS

Recognition: unknown

Dynamic solutions of next generation neural field models with delays

Carlo R. Laing, Oleh E. Omel'chenko

Pith reviewed 2026-05-08 12:36 UTC · model grok-4.3

classification 🌊 nlin.PS
keywords theta neuronsneural field modelsdelaysHopf bifurcationstraveling wavesbump solutionspattern formationself-consistency equations
0
0 comments X

The pith

Hopf bifurcations from delays in theta neuron rings generate traveling waves and breathing bumps

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies ring networks of theta neurons with delayed interactions in the continuum limit, modeled as next generation neural field equations. Stability of uniform states and localized bumps is analyzed, showing they destabilize via Hopf bifurcations as delay-related parameters change. This leads to the emergence of traveling waves and breathing bump solutions. These solutions fulfill self-consistency equations that the authors demonstrate can be solved efficiently. By continuing these solutions with varying parameters, the work provides a comprehensive view of how different delay types influence pattern formation in spatially extended theta neuron networks.

Core claim

Networks of theta neurons on a ring with distributed and conduction delays are considered. In the continuum limit these are described by next generation neural field models with delays. Uniform and bump states are shown to undergo Hopf bifurcations as delay parameters vary, creating traveling waves and breathing bump solutions. These dynamic solutions satisfy self-consistency equations, which are solved efficiently to follow the solutions across parameter space and map the effects of delays on pattern formation.

What carries the argument

Self-consistency equations for traveling waves and periodic breathing solutions derived from the delayed next generation neural field models

If this is right

  • Traveling waves arise from the Hopf bifurcations in the delayed models.
  • Breathing bump solutions with periodic time dependence are generated similarly.
  • Efficient numerical methods for the self-consistency equations enable continuation of these solutions with parameters.
  • The global picture shows distinct influences of finite-support distributed delays, infinite-support delays, and conduction delays on the patterns.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Extensions to two-dimensional domains or heterogeneous networks could reveal more complex delay-driven patterns.
  • The efficient solving technique might apply to other neural field models with delays for faster computation.
  • These results suggest delays as a control parameter for desired spatiotemporal activity in neural systems.

Load-bearing premise

The continuum limit provides an accurate description of the dynamics in the finite ring network of theta neurons with the modeled delays.

What would settle it

Numerical simulations of large but finite rings of theta neurons incorporating the delays should reproduce the predicted traveling waves and breathing bumps if the continuum approximation holds.

Figures

Figures reproduced from arXiv: 2604.22113 by Carlo R. Laing, Oleh E. Omel'chenko.

Figure 1
Figure 1. Figure 1: Blue: Hopf bifurcation (for eigenvalues satisfying (4.10)) of the spatially uniform state Z3. Red: saddle-node bifurcations of traveling wave solutions of (2.1)–(2.2). The spatially uniform state Z3 is stable in the blue shaded region. The dashed vertical lines refer to view at source ↗
Figure 2
Figure 2. Figure 2: Blue: Hopf bifurcations of the stable stationary bump. The bump is stable in the blue shaded region. Red curves show saddle-node bifurcations of periodic solutions. The dashed black line refers to view at source ↗
Figure 3
Figure 3. Figure 3: (a) Speed of traveling wave solution of (2.1)–(2.2) for two different values of τ (shown with dashed vertical lines in view at source ↗
Figure 4
Figure 4. Figure 4: Families of periodic solutions created in the Hopf bifurcations shown in view at source ↗
Figure 5
Figure 5. Figure 5: Blue: Hopf bifurcations of the spatially uniform steady state Z3 of (2.1) and (2.7). This state is stable in the blue shaded region. Red: saddle-node bifurcations of traveling waves. Parameters as in view at source ↗
Figure 6
Figure 6. Figure 6: Blue: Hopf bifurcations of a stable stationary bump. Red: saddle-node bifurcations of periodic solutions. Th bump is stable in the blue shaded region. Parameters as in view at source ↗
Figure 7
Figure 7. Figure 7: (a): a periodic solution of (2.1) with (2.7) for τ = 1.5. Firing frequency is shown in color, capped at 3. (b): continuation of the solution in panel (a) as τ is varied. a = 1 and other parameters as in view at source ↗
Figure 8
Figure 8. Figure 8: Hopf instabilities of a spatially uniform state of Eq. (2.1), (2.4), for perturbations with different wavenumbers, k. The state is stable when τ = 1/c = 0 and in the blue shaded region. The black dashed lines relate to later Figures. Parameters: η0 = 0.1 and other parameters as in view at source ↗
Figure 9
Figure 9. Figure 9: Traveling wave solutions of Eq. (2.1), (2.4). k is the angular wave number (2π divided by the wavelength). (a): wave speed v. (b): wave amplitude, defined as PM m=1 a 2 m + b 2 m 1/2 . The symbols are from numerical simulations of Eq. (2.1), (2.4). The waves with k = 4, 5 are numerically unstable. Parameters: τ = 0.5, η0 = 0.1, M = 20, and other parameters as in view at source ↗
Figure 10
Figure 10. Figure 10: Spatially-uniform periodic solutions of (2.1), (2.4). (a): Period T. (b): wave amplitude, defined as PM m=1 a 2 m + b 2 m 1/2 . Symbols are from numerical simulations of (2.1), (2.4). This type of solution loses stability to perturbations with spatial structure at a value of 1/c between 1.3 and 1.35, so is not stable upon its creation at the rightmost Hopf bifurcation. Parameters: τ = 1, η0 = 0.1, M = 5… view at source ↗
Figure 11
Figure 11. Figure 11: Solid: period, T, of periodic solutions of (2.1), (2.4) with spatial structure, for η0 = −0.15 (black) and η0 = −0.2 (magenta). Dashed: real part of the rightmost eigenvalues associated with the stability of a stationary bump solution, multiplied by 10 for the purpose of visualization. The dotted lines indicate where the eigenvalues cross the imaginary axis. Symbols show the result of numerically simulati… view at source ↗
Figure 12
Figure 12. Figure 12: A traveling wave solution of (B.1)–(B.3). Color shows (1 − cos θj ) 2 . Parameters as in view at source ↗
Figure 13
Figure 13. Figure 13: An approximately periodic solution of (B.1)–(B.3) arising from a Hopf bifurcation of a station￾ary bump. Color shows (1 − cos θj ) 2 . Parameters as in view at source ↗
Figure 14
Figure 14. Figure 14: Traveling waves with different angular wave numbers for (1.1), (1.6). (a): k = 3, 1/c = 0.5. (b): k = 4, 1/c = 1. (c): k = 5, 1/c = 1/5. Parameters as in view at source ↗
Figure 15
Figure 15. Figure 15: An approximately periodic solution for (1.1), (1.6). Parameters as in view at source ↗
read the original abstract

We study networks of theta neurons arranged on a ring with delayed interactions. In the continuum limit the systems are described by next generation neural field models with delays. We consider distributed delays with both finite and infinite support, and conduction delays. The stability of spatially uniform and localized bump states is determined, and we find that they undergo Hopf bifurcations as parameters related to the delays are varied. These bifurcations create traveling waves and ``breathing'' bump solutions. These dynamic solutions satisfy self-consistency equations and we show how to efficiently solve these equations. Following traveling waves and periodic solutions as parameters are varied provides a global picture of the influence of different delays on pattern formation processes in spatially extended networks of theta neurons.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The paper examines networks of theta neurons on a ring with delayed interactions. In the continuum limit, these yield next-generation neural field models incorporating distributed delays (finite and infinite support) and conduction delays. Stability of spatially uniform states and localized bumps is analyzed, revealing Hopf bifurcations as delay-related parameters vary. These bifurcations produce traveling waves and breathing bump solutions, which satisfy self-consistency equations that the authors solve efficiently. Parameter continuation of the dynamic solutions is used to map the global influence of different delay types on pattern formation.

Significance. If the derivations and reductions hold, the work provides a useful framework for tracking delay-induced spatiotemporal patterns in neural fields via self-consistency equations and continuation methods. This offers an efficient route to a global view of how distributed versus conduction delays shape traveling waves and oscillatory bumps, which is a clear methodological strength for exploring parameter dependence in spatially extended systems.

major comments (1)
  1. [model derivation and stability analysis] The central claims about Hopf bifurcations to traveling waves and breathing bumps rest on the N→∞ continuum limit commuting with the incorporation of distributed and conduction delays. The manuscript provides no finite-N simulations, error bounds, or analysis of finite-size corrections to the stability boundaries, which is load-bearing for the reported dynamic solutions (see the model derivation and stability sections).

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thoughtful summary and for recognizing the methodological value of using self-consistency equations and continuation to explore delay-induced patterns. We address the single major comment below.

read point-by-point responses

  1. Referee: [model derivation and stability analysis] The central claims about Hopf bifurcations to traveling waves and breathing bumps rest on the N→∞ continuum limit commuting with the incorporation of distributed and conduction delays. The manuscript provides no finite-N simulations, error bounds, or analysis of finite-size corrections to the stability boundaries, which is load-bearing for the reported dynamic solutions (see the model derivation and stability sections).

    Authors: We agree that direct numerical validation of the continuum limit in the presence of delays would strengthen the central claims. In the manuscript the delays are introduced at the level of the finite-N theta-neuron network on a ring; the Ott–Antonsen reduction is then applied to obtain the next-generation neural field equations, so that the N→∞ limit is taken after the delays have been incorporated. This ordering is standard for mean-field reductions of delayed pulse-coupled networks and ensures that the resulting continuum model inherits the delayed interactions. Nevertheless, the absence of finite-N benchmarks is a legitimate gap. In the revised version we will add a new subsection (or appendix) containing direct simulations of the finite-N system for representative values of the distributed-delay and conduction-delay parameters. These simulations will be used to compare the onset of Hopf bifurcations, the stability boundaries of uniform and bump states, and the existence of traveling-wave and breathing solutions against the continuum predictions. A short discussion of observed finite-size corrections will also be included. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper derives next-generation neural field equations from the continuum limit of a finite ring of theta neurons with distributed and conduction delays, performs linear stability analysis on uniform and bump states to locate Hopf bifurcations, and then solves the resulting self-consistency equations for traveling waves and breathing solutions. These steps are standard mathematical reductions from the model equations; the self-consistency relations are obtained directly from the field equations rather than being fitted to data or smuggled in via self-citation. No load-bearing premise reduces to a prior result by the same authors, and the continuum limit is treated as an independent modeling assumption whose validity is not claimed to be proven within the paper itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete; the work relies on standard continuum-limit assumptions and specific delay kernels whose independence from the target results cannot be checked.

axioms (1)
  • domain assumption Continuum limit of finite theta-neuron ring network yields a well-posed neural field equation
    Invoked when passing from discrete network to the next-generation field model with delays

pith-pipeline@v0.9.0 · 5414 in / 1226 out tokens · 29114 ms · 2026-05-08T12:36:45.945038+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

15 extracted references · 1 canonical work pages

  1. [1]

    Abrams, R

    [1]D. Abrams, R. Mirollo, S. Strogatz, and D. Wiley,Solvable model for chimera states of coupled oscillators, Phys. Rev. Lett., 101 (2008), p. 084103. [2]I. Al-Darabsah, S. A. Campbell, and B. Rahman,Distributed delay and desynchronization in a neural mass model, SIAM Journal on Applied Dynamical Systems, 23 (2024), pp. 3013–3051. [3]S.-i. Amari,Dynamics ...

    2008

  2. [2]

    [5]J. T. Ariaratnam and S. H. Strogatz,Phase diagram for the winfree model of coupled nonlinear oscillators, Physical Review Letters, 86 (2001), p

    2001

  3. [3]

    [6]F. M. Atay and A. Hutt,Stability and bifurcations in neural fields with finite propagation speed and general connectivity, SIAM Journal on Applied Mathematics, 65 (2005), pp. 644–666. [7]F. M. Atay and A. Hutt,Neural fields with distributed transmission speeds and long-range feedback delays, SIAM Journal on Applied Dynamical Systems, 5 (2006), pp. 670–...

    2005

  4. [4]

    Byrne, R

    ´A. Byrne, R. D. O’Dea, M. Forrester, J. Ross, and S. Coombes,Next-generation neural mass and field modeling, Journal of neurophysiology, (2020). [13]S. A. Campbell,Time delays in neural systems, in Handbook of brain connectivity, Springer, 2007, pp. 65–90. [14]S. Coombes,Waves, bumps, and patterns in neural field theories, Biol. Cybern., 93 (2005), pp. 9...

    2020

  5. [5]

    Coombes and C

    [18]S. Coombes and C. Laing,Delays in activity-based neural networks, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 367 (2009), pp. 1117–1129. [19]S. Coombes and H. Meijer,Synchrony in firing rate neural networks with multiple delays: A harmonic balance approach, arXiv preprint arXiv:2509.22495, (2025)...

    work page arXiv 2009

  6. [6]

    Ermentrout,Type I membranes, phase resetting curves, and synchrony, Neural computation, 8 (1996), pp

    [24]B. Ermentrout,Type I membranes, phase resetting curves, and synchrony, Neural computation, 8 (1996), pp. 979–1001. [25]B. Ermentrout,Neural networks as spatio-temporal pattern-forming systems, Rep. Prog. Phys., 61 (1998), pp. 353–430. [26]G. B. Ermentrout and N. Kopell,Parabolic bursting in an excitable system coupled with a slow oscillation, SIAM Jou...

    1996

  7. [7]

    Hartung, T

    [31]F. Hartung, T. Krisztin, H.-O. Walther, and J. Wu,Functional differential equations with state- dependent delays: theory and applications, Handbook of differential equations: ordinary differential equations, 3 (2006), pp. 435–545. [32]Z. Huang, Y. Zhang, P. Zou, X. Zong, and Q. Zhang,Myelin dysfunction in aging and brain disorders: mechanisms and ther...

    2006

  8. [8]

    [33]P. J. Hurtado and A. S. Kirosingh,Generalizations of the ‘linear chain trick’: incorporating more flexible dwell time distributions into mean field ode models, Journal of mathematical biology, 79 (2019), pp. 1831–1883. [34]J.-S. Kang, Y.-H. An, R.-S. Kim, and C.-U. Choe,Periodic external driving of transmission-delay- coupled phase oscillator system: ...

    2019

  9. [9]

    Laing and C

    [36]C. Laing and C. Chow,Stationary bumps in networks of spiking neurons, Neural Comput., 13 (2001), pp. 1473–1494. [37]C. R. Laing,The dynamics of chimera states in heterogeneous Kuramoto networks, Physica D, 238 (2009), pp. 1569–1588. [38]C. R. Laing,Derivation of a neural field model from a network of theta neurons, Physical Review E, 90 (2014), p. 010...

    2001

  10. [10]

    [40]C. R. Laing,Exact neural fields incorporating gap junctions, SIAM Journal on Applied Dynamical Systems, 14 (2015), pp. 1899–1929. [41]C. R. Laing and B. Krauskopf,Periodic solutions for a pair of delay-coupled excitable theta neurons, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 481 (2025), p. 20240897. [42]C. R...

    2015

  11. [11]

    [49]O. E. Omel’chenko,Periodic orbits in the Ott–Antonsen manifold, Nonlinearity, 36 (2023), pp. 845–861. [50]O. E. Omel’chenko and C. R. Laing,Activity patterns in ring networks of quadratic integrate-and-fire neurons with synaptic and gap junction coupling, Physical Review E, 110 (2024), p. 034411. [51]E. Ott and T. M. Antonsen,Low dimensional behavior ...

    2023

  12. [12]

    Paz ´o and E

    [54]D. Paz ´o and E. Montbri´o,Low-dimensional dynamics of populations of pulse-coupled oscillators, Phys- ical Review X, 4 (2014), p. 011009. [55]M. R. Qubbaj and V. K. Jirsa,Neural field dynamics under variation of local and global connectivity and finite transmission speed, Physica D, 238 (2009), pp. 2331–2346. [56]B. Rahman, K. B. Blyuss, and Y. N. Ky...

    2014

  13. [13]

    [69]S. A. van Gils, S. G. Janssens, Y. A. Kuznetsov, and S. Visser,On local bifurcations in neural field models with transmission delays, Journal of mathematical biology, 66 (2013), pp. 837–887. [70]R. Veltz,Interplay between synaptic delays and propagation delays in neural field equations, SIAM Journal on Applied Dynamical Systems, 12 (2013), pp. 1566–16...

    2013

  14. [14]

    [72]N. A. Venkov, S. Coombes, and P. C. Matthews,Dynamic instabilities in scalar neural field equations with space-dependent delays, Physica D, 232 (2007), pp. 1–15. [73]S. Visser, R. Nicks, O. Faugeras, and S. Coombes,Standing and travelling waves in a spherical brain model: the nunez model revisited, Physica D, 349 (2017), pp. 27–45. [74]H. R. Wilson an...

    2007

  15. [15]

    Zhang,Representation of spatial orientation by the intrinsic dynamics of the head-direction cell ensemble: a theory, The journal of neuroscience, 16 (1996), pp

    [78]K. Zhang,Representation of spatial orientation by the intrinsic dynamics of the head-direction cell ensemble: a theory, The journal of neuroscience, 16 (1996), pp. 2112–2126

    1996