Pith. sign in

REVIEW 3 major objections 5 minor 40 references

Spectral theory for population density dynamics of spiking neurons with refractoriness

T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

Pith's one-line read This paper establishes a complete spectral theory for the Fokker-Planck operator governing spiking neural populations with a finite absolute refractory period, showing that the spectrum is discrete and generated by a single characteristic e

desk verdict Serious spectral-theory paper with a real result on refractoriness, but the central theorem leans on an unverified regularity assumption and the abstract overclaims the limit-cycle result. read the letter →

arxiv 2607.20699 v1 pith:TWLUJRWK submitted 2026-07-22 q-bio.NC math-phmath.MPmath.SP

classification q-bio.NCmath-phmath.MPmath.SP MSC 47A1035P0547D0692C20
keywords populationdensitydynamicsspikingneuronsabsoluterefractoryperiodFokker-Planckoperatornon-self-adjointspectraltheoryboundaryeigenvalueproblemlinearresponsetransferfunctionlimitcycles
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper proves that population density dynamics of spiking neurons with an absolute refractory period can be studied rigorously as a non-self-adjoint boundary eigenvalue problem. By adding a refractory compartment that stores the past firing rate, the dynamics become Markovian and the generator is shown to have a purely discrete spectrum: every eigenvalue solves an explicit characteristic equation, zero is always a simple eigenvalue, and every eigenspace is one-dimensional. This justifies spectral decomposition methods that were previously used heuristically, and yields an exact linear-response transfer function including new boundary-modulated terms. As an application, the paper shows that increasing refractoriness can drive a stable population into self-sustained limit-cycle oscillations by moving diffusion poles across the imaginary axis.

What carries the argument

The central object is the evolution operator Tγ acting on L^p of three compartments: (α,H), (H,θ) and the refractory interval (0,τ0), with five linear boundary conditions coupling the membrane potential density to the stored firing-rate history. The argument runs through the theory of non-self-adjoint boundary eigenvalue problems: a holomorphic operator function eTγ(λ) is shown to be Fredholm and globally equivalent to a 5×5 characteristic matrix Mγ(λ); the determinant of that matrix divided by Wronskians yields the scalar characteristic equation Δγ(λ)=0 that completely determines the spectrum. Holomorphic root functions built from the adjoint fundamental solutions supply the eigenfunctions

What would settle it

Discretize Tγ for the LIF model with τ0>0 (finite-element or spectral method on the three compartments) and compute the largest eigenvalues; compare with zeros of Δγ(λ). Finding a spectral eigenvalue that is not a zero of Δγ, or a zero of Δγ with no matching eigenvalue, would refute Theorem 3.3. Also check whether Δγ'(0) can vanish for some admissible γ, which would contradict the asserted simplicity of the zero eigenvalue.

Watch

Extended reading notes

Core claim

The central result is Theorem 3.3: for fixed input moments, the evolution operator Tγ has a discrete spectrum given by the zeros of Δγ(λ)=f1(θ,λ)/wγ(θ)−e^{−τ0λ}f1(H,λ)/wγ(H), with 0 always a simple eigenvalue and all eigenvalues of geometric multiplicity one. The characteristic equation emerges from the determinant of a 5×5 characteristic matrix built from fundamental solutions of the Fokker-Planck operator and an exponential factor encoding the refractory delay. Because the boundary eigenvalue operator function is globally equivalent to this matrix, eigenvalues, algebraic multiplicities, and generalized eigenfunctions of the operator are read off from the matrix. This converts previously he

Load-bearing premise

The proof leans on an abstract regularity condition for the boundary conditions (from the theory of non-self-adjoint boundary eigenvalue problems) that the paper invokes without verifying explicitly.

Editorial extensions

If this is right

  • Spectral mode expansions of population activity for neurons with τ0>0 are justified, placing low-dimensional firing-rate models on a rigorous footing.
  • The exact transfer function reveals two previously missed contributions: noise-modulation at threshold and the presence of the reflecting barrier at α.
  • Defective eigenvalues are identified as exceptional points where relaxational modes coalesce into oscillatory modes; the full Jordan-chain description is required near them.
  • In recurrent networks with delay, increasing the refractory period can destabilize the fixed point via diffusion poles, producing limit-cycle firing rates at realistic parameters without transmission delays.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the boundary regularity hypothesis used to invoke the abstract boundary-eigenvalue machinery fails for some parameter range, the characteristic equation would not necessarily capture the full spectrum; this is a concrete target for numerical verification.
  • The new boundary-modulated terms in the transfer function should also appear in conductance-based neuron models where noise is state-dependent; the paper suggests but does not demonstrate this.
  • A natural next step is to prove completeness of the biorthogonal eigenbasis via resolvent asymptotics; if achieved, the spectral decomposition would become a fully rigorous replacement for time-stepping simulations in these networks.
  • The Hopf-bifurcation scenario for excitatory populations suggests that the refractory period is a singular perturbation that can change network state; testing whether this persists in finite-size networks would check robustness.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper develops a spectral theory for the population-density Fokker-Planck operator of integrate-and-fire neurons with an absolute refractory period. The state space is augmented by a refractory-age density, and the evolution operator is formulated as a non-self-adjoint boundary eigenvalue problem on three compartments. The main claim (Theorem 3.3) is that the spectrum is purely discrete and coincides with the zeros of a characteristic function Δγ(λ) (Eq. 3.14), that 0 is always a simple eigenvalue, and that every eigenvalue has geometric multiplicity one. A transfer function is derived from the resolvent, including new boundary-modulation terms, and the theory is applied to linear stability of recurrent networks. The paper also identifies defective eigenvalues as exceptional points and presents numerical evidence that refractoriness can promote limit cycles.

Significance. If the central spectral theorem holds, this is a substantial contribution: it would put spectral decomposition methods for population density dynamics on a rigorous footing for arbitrary finite refractory periods, and it yields an explicit characteristic equation that reduces correctly to known results in the τ0→0 limit. The paper is self-contained from the FP model, introduces no free parameters, and makes falsifiable predictions about the location of eigenvalues and the structure of the transfer function. The treatment of defective eigenvalues and the corrected boundary contributions to the transfer function are also valuable. The numerical checks are consistent with the analytic expressions and the authors are candid about the non-rigorous nature of the nonlinear bifurcation claim.

major comments (3)
  1. [Lemma 3.1 and Theorem 3.3 (Eq. 3.8, Eq. 3.13)] The proof of the central spectral characterization is delegated to [28] without verifying the hypotheses. The boundary operator Bγ in Eq. (2.11) contains nonlocal transmission conditions (BC3, BC4) coupling p(θ) to p_r(0) and p(H±) to p_r(τ0). Regularity in the sense of [28] is not automatic for such non-self-adjoint boundary eigenvalue problems. Please verify the regularity assumptions of [28, Thm. 1.11.1 and Lemma 1.11.2] for Bγ, or state and prove them. In particular, show that the global equivalence (3.8) holds with invertible Cγ(λ), Dγ(λ). Without this, Theorem 3.3's identification of σ(Tγ) with the zeros of Δγ, and the algebraic/geometric multiplicity statements used in Lemma 3.4 and in Eq. (4.15), are not established.
  2. [Theorem 3.3, geometric multiplicity argument] The proof that every eigenvalue has geometric multiplicity one is compressed: 'if some λ_n had geometric multiplicity all minors of Mγ(λn) would vanish, forcing f1(θ,λn)=f1(H,λn)=0; the boundary conditions then give f2(θ,λn)=f2(H,λn)=0'. This step is not immediate for a 5×5 characteristic matrix with the nonlocal structure of (3.12). Please provide the explicit minor computation, since this uniqueness is used in the root-function construction (Lemma 3.4) and in the normalization of the eigenfunctions.
  3. [Section 4.3 and Abstract] The abstract states that refractoriness 'can facilitate the onset of limit cycles, that is, stable oscillations in the firing rate'. However, Section 4.3 explicitly says: 'the spectral analysis proves only the loss of linear stability of the fixed point; the emergence of a stable limit cycle is verified numerically and not proven here, as this would require a nonlinear bifurcation analysis beyond the present scope.' The abstract should be rephrased to distinguish the rigorous Hopf linear-stability result from the numerical observation of stable oscillations.
minor comments (5)
  1. [Eq. (2.3)] The second equation is written as ∂_t p^r_t(τ) = −∂_τ p^r_t(τ,t), but the argument (τ,t) is inconsistent with p^r_t(τ).
  2. [Section 3.4] The reference to the Lumer-Phillips theorem is misspelled as '[13, Therem 3.15]'. Also, the authors should mention explicitly that they are using the spectrum result of Proposition 3.9 to obtain the range condition required by Lumer-Phillips.
  3. [Proposition 3.9 proof (Appendix B)] The displayed form of an eigenfunction for eigenvalue iω writes the refractory component as exp(iτ ω), while the construction in Eq. (3.17) gives exp(−iτ s). The subsequent integral uses cos(n2πτ/τ0), which is consistent with the latter convention; please fix the sign for clarity.
  4. [Discussion, completeness] The paper uses the phrase 'complete spectral characterization' in the abstract, but the discussion lists 'the rigorous proof of the completeness of the eigenbasis' as an open problem. Consider softening 'complete' to 'explicit' or 'full point-spectrum' to avoid overstating what is proven.
  5. [Throughout] There are some typos, e.g., 'redunction' in Section 3.2 and 'spiking neuron networks' in reference [27]. A careful proofreading pass is recommended.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the spectral derivation is self-contained, with self-citations used only as baselines and one external-theorem applicability caveat.

full rationale

The central derivation chain is: define the augmented state space and evolution operator T_gamma (Section 2), recast the eigenvalue problem as an abstract boundary eigenvalue problem, introduce the characteristic matrix M_gamma, derive the characteristic equation Delta_gamma(lambda)=0 (Theorem 3.3), construct root functions from M_gamma and its adjoint (Lemma 3.4, Appendix A), and then obtain the transfer function from the resolvent (Lemma 4.1) and the coupled-system characteristic equation via a rank-one perturbation (Theorem 4.3). No parameter is fitted and later called a prediction; no 'uniqueness theorem' from the authors' own prior work is invoked to force a choice; no ansatz is smuggled in via self-citation. The characteristic equation is computed from the boundary operator B_gamma and the fundamental solutions of L_gamma f = lambda f, not assumed. Self-citations [25,40] are used as baselines to generalize ('this equation reduces to the known expression [20,25]', 'This formula generalizes the expression obtained in [40, Eq. (28)]', 'We refer to [40] for details in the derivation'), and as prior phenomenological context; they do not already contain the finite-refractoriness spectral characterization that the paper derives. The tau0 -> 0 limit is used as a consistency check, not as an input. The main caveat is Lemma 3.1, which is justified as 'a restatement of Theorem 1.11.1 and Lemma 1.11.2 in [28]' without explicitly verifying that the boundary operator B_gamma is regular in the sense of [28]. This is a correctness / applicability risk: if regularity fails, the equivalence of sigma(T_gamma) with zeros of det M_gamma could break. But it is not circularity, because [28] is an independent external mathematical source, not the paper's own conclusion. The paper itself also flags as open 'the rigorous proof of the completeness of the eigenbasis and a detailed characterization of the regularity properties of the generated semigroup', which again are rigor gaps, not circular reduction. For those reasons, the derivation is not circular and receives a low score, though not zero because of the minor non-load-bearing self-citations and the external-theorem applicability caveat.

Assumptions & free parameters 0 free parameters · 5 assumptions · 1 invented entities

The core derivation introduces no fitted free parameters; the results are consequences of the FP model with the stated boundary conditions. The main input assumptions are (i) the validity of the FP/diffusion approximation, (ii) the fixed deterministic refractory period, (iii) the applicability of the Mennicken–Möller boundary-eigenvalue framework (Lemma 3.1), and (iv) the extended mean-field approximation for coupled networks. The augmented refractory density is a bookkeeping variable, not a new physical entity.

assumptions (5)
  • domain assumption The Fokker–Planck equation (1.1) with boundary conditions (1.3)–(1.4) is a valid thermodynamic-limit description of an ensemble of integrate-and-fire neurons with absolute refractoriness.
    This is the object of study; the microscopic derivation is cited to [24] and not addressed here.
  • domain assumption The absolute refractory period τ0 is a fixed deterministic delay after each spike, after which neurons reset to H; the past firing rate is stored via pr_t(τ)=ν(t−τ) (2.1).
    Basis for the augmented state space and the coupled transport equation (2.2)–(2.4).
  • standard math The boundary eigenvalue problem (3.1)–(3.7) satisfies the hypotheses of [28, Thm. 1.11.1 and Lemma 1.11.2], so eTγ is Fredholm and globally equivalent to Mγ via invertible Cγ,Dγ (Lemma 3.1).
    The paper states this is a restatement of the cited theorem but does not verify the regularity hypotheses; the entire spectral characterization (Theorem 3.3) depends on this equivalence.
  • domain assumption For coupled populations, the extended mean-field approximation γ(t)=γ·ν(t−δ)+γ0 (4.11) gives the input moments.
    Used in Secs. 4.2–4.3 to derive the closed-loop characteristic equation H(s)=e^{δs} and the limit-cycle criterion.
  • domain assumption The synaptic input is approximated by white noise with infinitesimal moments (µ(t),D(t)) (diffusion approximation).
    Standard in the population density approach; introduced in Sec. 1 and used throughout.
invented entities (1)
  • Refractory-state density pr_t(τ)
    purpose: Adds a history variable so the system becomes Markovian; it stores the firing rate τ0 intervals in the past (Eq. 2.1) and obeys a transport equation (2.2).
    This is a mathematical bookkeeping device, not a new physical quantity; it is fully determined by the firing rate history and has no independent falsifiable signature.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Spectral theory for population density dynamics of spiking neurons with refractoriness." pith.science (2026). https://pith.science/paper/TWLUJRWK

@misc{pith2026260720699,
  author       = {Pith},
  title        = {Pith review of: Spectral theory for population density dynamics of spiking neurons with refractoriness},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TWLUJRWK}},
  note         = {Machine review of arXiv:2607.20699}
}
read the original abstract

Incorporating an absolute refractory period into the population density approach for spiking neurons remains an open problem, despite evidence that refractoriness can strongly affect nonlinear transfer functions and network stability. We develop a rigorous operator-theoretic framework for neuronal population dynamics with a finite refractory time by augmenting the state space to include refractory history and formulating the problem as a non-self-adjoint boundary eigenvalue problem for the Fokker-Planck operator. This yields a complete spectral characterization of the generator, proves dissipativity and the existence of a contraction semigroup, and identifies defective eigenvalues as exceptional points where oscillatory modes emerge from coalescing relaxational modes. Within the framework of linear response theory, we also derive an exact transfer function that accounts for boundary conditions modulated by external input, correcting previous heuristic derivations and revealing additional threshold-noise contributions. Using this transfer function under a mean-field approximation, we further show that refractoriness in populations of interacting neurons can facilitate the onset of limit cycles, that is, stable oscillations in the firing rate. These results provide a rigorous foundation for spectral decomposition methods in computational neuroscience, opening the way to their further rigorous mathematical analysis.

Figures

Figures reproduced from arXiv: 2607.20699 by the authors.

Figure 1
Figure 1. Eigenvalues λn of the evolution operator for different values of refractory period τ0. Current mean µ is varied at 0.1 mV steps keeping D = 8mV2 . (a) Real part of λn; (b) imaginary part. Red and blue branches are mixed and real λn, respectively. Different shadings represent τ0 from 0 τm (light) to 0.4 τm (dark). Parameters: τm = 1, H = 0 mV, θ = 20 mV, α = −5θ = −100 mV, 10 [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Defective eigenvalues at the real-to-complex transition in a mixed branch. (a) Real part [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Module (left) and phase (right) of the transfer function in response to a noise ( [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: (a) Locus of the leading eigenvalues (poles) of the coupled system as the refractory [PITH_FULL_IMAGE:figures/full_fig_p020_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references

  1. [28]

    Gulf Professional Publishing, 2003

    Reinhard Mennicken and Manfred Möller.Non-self-adjoint boundary eigenvalue problems, vol- ume 192. Gulf Professional Publishing, 2003

  2. [1]

    Asynchronous states in networks of pulse-coupled oscillators.Physical Review E, 48(2):1483–1490, 1993

    Larry F Abbott and Carl van Vreeswijk. Asynchronous states in networks of pulse-coupled oscillators.Physical Review E, 48(2):1483–1490, 1993

  3. [2]

    Model of global spontaneous activity and local structured activity during delay periods in the cerebral cortex.Cereb

    Daniel J Amit and Nicolas Brunel. Model of global spontaneous activity and local structured activity during delay periods in the cerebral cortex.Cereb. Cortex, 7(3):237–52, 1997

  4. [3]

    Quantitative study of attractor neural network retrieving at low spike rates: I

    Daniel J Amit and Misha Tsodyks. Quantitative study of attractor neural network retrieving at low spike rates: I. substrate—spikes, rates and neuronal gain.Netw. Comput. Neural Syst., 2(3):259–273, 1991

  5. [4]

    Low-dimensional spike rate models derived from networks of adaptive integrate-and-fire neurons: comparison and implementation.PLoS Comput

    Moritz Augustin, Josef Ladenbauer, Fabian Baumann, and Klaus Obermayer. Low-dimensional spike rate models derived from networks of adaptive integrate-and-fire neurons: comparison and implementation.PLoS Comput. Biol., 13(6):e1005545, 2017

  6. [5]

    A K Peters, Wellesley, MA, 2005

    András Bátkai and Susanna Piazzera.Semigroups for Delay Equations, volume 10 ofResearch Notes in Mathematics. A K Peters, Wellesley, MA, 2005. 25

  7. [6]

    Ergodic behavior of diffusions with random jumps from the boundary.Stochastic processes and their applications, 119(3):864–881, 2009

    Iddo Ben-Ari and Ross G Pinsky. Ergodic behavior of diffusions with random jumps from the boundary.Stochastic processes and their applications, 119(3):864–881, 2009

  8. [7]

    Refractoriness and neural precision.J

    Michael J Berry II and Markus Meister. Refractoriness and neural precision.J. Neurosci., 18(6):2200–2211, 1998

Show all 40 references
  1. [8]

    Metastable dynamics of neural circuits and networks.Applied Physics Reviews, 9(1):011313, 2022

    Braden A W Brinkman, Han Yan, Arianna Maffei, Il Memming Park, Alfredo Fontanini, Jin Wang, and Giancarlo La Camera. Metastable dynamics of neural circuits and networks.Applied Physics Reviews, 9(1):011313, 2022

  2. [9]

    Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons.Journal of computational neuroscience, 8(3):183–208, 2000

    Nicolas Brunel. Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons.Journal of computational neuroscience, 8(3):183–208, 2000

  3. [10]

    Fast global oscillations in networks of integrate-and-fire neurons with low firing rates.Neural computation, 11(7):1621–1671, 1999

    Nicolas Brunel and Vincent Hakim. Fast global oscillations in networks of integrate-and-fire neurons with low firing rates.Neural computation, 11(7):1621–1671, 1999

  4. [11]

    Analysis of nonlinear noisy integrate & fire neuron models: blow-up and steady states.The Journal of Mathematical Neuroscience, 1:1–33, 2011

    María J Cáceres, José A Carrillo, and Benoît Perthame. Analysis of nonlinear noisy integrate & fire neuron models: blow-up and steady states.The Journal of Mathematical Neuroscience, 1:1–33, 2011

  5. [12]

    Solving the two-dimensional fokker-planck equation for strongly correlated neurons.Physical Review E, 95(1):012412, 2017

    Taşkın Deniz and Stefan Rotter. Solving the two-dimensional fokker-planck equation for strongly correlated neurons.Physical Review E, 95(1):012412, 2017

  6. [13]

    One-parameter semigroups for linear evolution equa- tions.Semigroup Forum, 63:278–280, 1999

    Klaus-Jochen Engel and Rainer Nagel. One-parameter semigroups for linear evolution equa- tions.Semigroup Forum, 63:278–280, 1999

  7. [14]

    Collective behavior of networks with linear (vlsi) integrate- and-fire neurons.Neural Computation, 11(3):633–652, 1999

    Stefano Fusi and Maurizio Mattia. Collective behavior of networks with linear (vlsi) integrate- and-fire neurons.Neural Computation, 11(3):633–652, 1999

  8. [15]

    Collective behavior of networks with linear (VLSI) integrate- and-fire neurons.Neural Computation, 11(3):633–52, 1999

    Stefano Fusi and Maurizio Mattia. Collective behavior of networks with linear (VLSI) integrate- and-fire neurons.Neural Computation, 11(3):633–52, 1999

  9. [16]

    Time structure of the activity in neural network models.Phys

    Wulfram Gerstner. Time structure of the activity in neural network models.Phys. Rev. E, 51:738–758, Jan 1995

  10. [17]

    Kistler, Richard Naud, and Liam Paninski.Neuronal dynamics: From single neurons to networks and models of cognition

    Wulfram Gerstner, Werner M. Kistler, Richard Naud, and Liam Paninski.Neuronal dynamics: From single neurons to networks and models of cognition. Cambridge University Press, 2014

  11. [18]

    Brownian motion on the figure eight.J

    Ilie Grigorescu and Min Kang. Brownian motion on the figure eight.J. Theor. Prob., 15(3):817– 844, 2002

  12. [19]

    Dynamics of encoding in neuron populations: some general mathematical features.Neural Computation, 12(3):473–518, 2000

    Bruce W Knight. Dynamics of encoding in neuron populations: some general mathematical features.Neural Computation, 12(3):473–518, 2000

  13. [20]

    Dynamical models of interacting neuron populations in visual cortex

    Bruce W Knight, Dimitri Manin, and Lawrence Sirovich. Dynamical models of interacting neuron populations in visual cortex. In E.C. Gerf, editor,Symposium on Robotics and Cyber- netics: Computational Engineering in Systems Applications, pages 1–5, Cite Scientifique, Lille, Fran...

  14. [21]

    One dimensional Brownian motion with holding and jumping boundary.Math

    Yuk-J Leung. One dimensional Brownian motion with holding and jumping boundary.Math. Meth. Appl. Sci., 47(17):13051–13062, 2024

  15. [22]

    Effects of noise in excitable systems.Phys

    Benjamin Lindner, Jordi García-Ojalvo, Alexander B Neiman, and Lutz Schimansky-Geier. Effects of noise in excitable systems.Phys. Rep., 392(6):321–424, 2004. 26

  16. [23]

    Transmission of noise coded versus additive signals through a neuronal ensemble.Phys

    Benjamin Lindner and Lutz Schimansky-Geier. Transmission of noise coded versus additive signals through a neuronal ensemble.Phys. Rev. Lett., 86(14):2934–2937, 2001

  17. [24]

    Rigorous justification of the fokker–planck equations of neural networks based on an iteration perspective.SIAM Journal on Mathematical Analysis, 54(1):1270–1312, 2022

    Jian-Guo Liu, Ziheng Wang, Yuan Zhang, and Zhennan Zhou. Rigorous justification of the fokker–planck equations of neural networks based on an iteration perspective.SIAM Journal on Mathematical Analysis, 54(1):1270–1312, 2022

  18. [25]

    Population dynamics of interacting spiking neurons

    Maurizio Mattia and Paolo Del Giudice. Population dynamics of interacting spiking neurons. Phys. Rev. E, 66(5):051917, 2002

  19. [26]

    Finite-size dynamics of inhibitory and excitatory interacting spiking neurons.Phys

    Maurizio Mattia and Paolo Del Giudice. Finite-size dynamics of inhibitory and excitatory interacting spiking neurons.Phys. Rev. E, 70(5 Pt 1):052903, 2004

  20. [27]

    Low dimensional dynamics of spiking neuron networks

    Maurizio Mattia and Gianni V Vinci. Low dimensional dynamics of spiking neuron networks. Zenodo, page 5518215, 2021

  21. [29]

    Dynamics of neuronal populations: The equilibrium solution.SIAM Journal on Applied Mathematics, 60(6):2009–2028, 2000

    Ahmet Omurtag, Bruce W Knight, and Lawrence Sirovich. Dynamics of neuronal populations: The equilibrium solution.SIAM Journal on Applied Mathematics, 60(6):2009–2028, 2000

  22. [30]

    Jun Peng and Wenbo V. Li. Diffusions with holding and jumping boundary.Science China Mathematics, 56:161 – 176, 2012

  23. [31]

    Low-dimensional firing-rate dynamics for populations of renewal-type spiking neurons.Physical Review E, 102(2):022407, 2020

    Bastian Pietras, Noé Gallice, and Tilo Schwalger. Low-dimensional firing-rate dynamics for populations of renewal-type spiking neurons.Physical Review E, 102(2):022407, 2020

  24. [32]

    First-passage-time density and moments of the Ornstein- Uhlenbeck process.J

    Luigi M Ricciardi and Shunsuke Sato. First-passage-time density and moments of the Ornstein- Uhlenbeck process.J. Appl. Prob., 25(1):43–57, 1988

  25. [33]

    Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons.Phys

    Magnus J E Richardson. Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons.Phys. Rev. E, 69(5 Pt 1):051918, 2004

  26. [34]

    Firing-rate response of linear and nonlinear integrate-and-fire neurons to modulated current-based and conductance-based synaptic drive.Phys

    Magnus J E Richardson. Firing-rate response of linear and nonlinear integrate-and-fire neurons to modulated current-based and conductance-based synaptic drive.Phys. Rev. E, 76(2 Pt 1):021919, 2007

  27. [35]

    Springer Berlin Heidelberg, Berlin, Heidelberg, 1989

    Hannes Risken.The Fokker-Planck Equation, volume 18 ofSpringer Series in Synergetics. Springer Berlin Heidelberg, Berlin, Heidelberg, 1989

  28. [36]

    Modulated escape from a metastable state driven by colored noise.Phys

    Jannis Schuecker, Markus Diesmann, and Moritz Helias. Modulated escape from a metastable state driven by colored noise.Phys. Rev. E, 92(5):052119, 2015

  29. [37]

    Mind the last spike — firing rate models for mesoscopic populations of spiking neurons.Curr

    Tilo Schwalger and Anton V Chizhov. Mind the last spike — firing rate models for mesoscopic populations of spiking neurons.Curr. Opin. Neurobiol., 58:155–166, 2019

  30. [38]

    Mean-field analysis of neuronal spike dynamics.Network, 4(3):259–84, 1993

    Alessandro Treves. Mean-field analysis of neuronal spike dynamics.Network, 4(3):259–84, 1993

  31. [39]

    Cambridge University Press, 1988

    Henry Clavering Tuckwell.Introduction to theoretical neurobiology: linear cable theory and dendritic structure, volume 1. Cambridge University Press, 1988

  32. [40]

    Rosetta stone for the population dynamics of spiking neuron networks.Physical Review E, 110(3):034303, 2024

    Gianni V Vinci and Maurizio Mattia. Rosetta stone for the population dynamics of spiking neuron networks.Physical Review E, 110(3):034303, 2024. 27

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.