Pith. sign in

REVIEW 4 major objections 4 minor 62 references

Hyperchaos and complex dynamical regimes in $N$-dimensional neuron lattices

T0 review · 4 major / 4 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read As electrical coupling grows, nonchaotic Rulkov neurons on N-dimensional lattices follow a single regime sequence from independent spiking to synchronized hyperchaos, with dimensionality shifting where chaos first appears.

desk verdict Honest N-dimensional Rulkov lattice extension; the dimensional peak-shift trend is plausible but rests on finite-time Lyapunov exponents without convergence evidence or error bars. read the letter →

arxiv 2505.03051 v2 pith:FKHRJR4G submitted 2025-05-05 nlin.CD cond-mat.dis-nnmath.DSq-bio.NC

Brandon B. Le , Dima Watkins This is my paper · ORCID
classification nlin.CDcond-mat.dis-nnmath.DSq-bio.NC MSC 37D4537M2537N25 PACS 05.45.-a05.45.Xt
keywords RulkovneuronsN-dimensionallatticeshyperchaosLyapunovexponentssynchronizedburstinglagsynchronizationelectricalcouplingnext-nearest-neighbor
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

The paper sets out to show that a minimal discrete-time neuron model, the nonchaotic Rulkov map, placed on cubic lattices in one to four spatial dimensions, develops a reproducible sequence of collective dynamical regimes as electrical coupling strength $g$ is swept from 0 to 1. In homogeneous lattices the sequence runs from uncoupled nonchaotic spiking through unsynchronized chaotic spiking and synchronized chaotic bursting to synchronized hyperchaos; heterogeneous lattices add a low-coupling 'local quasi-bursting' phase. Dimensionality and connectivity do not alter the sequence itself, but the paper argues they shift its thresholds: the first chaotic peak in the maximal Lyapunov exponent moves to larger $g$ and higher values as $N$ grows, and next-nearest-neighbor coupling delays the onset of synchronized bursting through a diluting 'destructive interference' of many neighbors. At maximal coupling, large 2D and 3D lattices display a checkerboard of adjacent neurons alternating high and low voltage each timestep, an extreme form of one-timestep lag synchronization. A reader should care because this suggests that spatial wiring, not the internal complexity of the neuron model, can generate dimension-dependent chaos and synchronization, and it gives a computationally cheap test bed for synchronization phenomena.

What carries the argument

The carrying object is the lattice state tensor $X_{ia}=(x_i,y_i)$ of Rulkov maps with electrical coupling $C_i(k)=\frac{g}{|N_i|}\sum_{j\in N_i}(x_j-x_i)$, where $|N_i|=2N$ for nearest neighbors and $4N$ for next-nearest neighbors, with periodic boundary conditions. The iteration function feeds this current into both the fast variable (through an effective shift of the slow variable) and the slow variable (through an effective shift of the excitation parameter), which is the mechanism behind the high-frequency spiking and subsequent quiescence that defines synchronized chaotic bursting. The argument is carried quantitatively by an explicitly derived tensorial Jacobian, converted to a matrix via a base-$\zeta$ index map, whose QR-based Lyapunov spectra provide the regime boundaries. Two mechanisms explain the key patterns: 'destructive interference' dilutes each pairwise interaction when a neuron has more neighbors, delaying synchronization, and the discrete-time lag in current flow reverses the voltage difference every timestep at $g=1$, producing checkerboard lag synchronization.

What would settle it

Recompute the maximal Lyapunov exponent sweep for the $N=2$ homogeneous nearest-neighbor lattice with $\zeta=8$ using orbit length $k=100000$ and fifty independent initial conditions; if the sharp drop near $g\approx 0.2$, the rightward peak shift with $N$, and the small bumps in the next-nearest-neighbor curves do not survive these longer runs, the central regime-transition claims are called into question.

Watch

Extended reading notes

Core claim

The central claim is that electrically coupled nonchaotic Rulkov neurons on an $N$-dimensional cubic lattice, with nearest-neighbor or next-nearest-neighbor coupling, pass through the same basic sequence of regimes as the coupling conductance $g$ is increased. For homogeneous parameters the regimes are uncoupled nonchaotic spiking, unsynchronized chaotic spiking, synchronized chaotic bursting, and synchronized hyperchaos; for heterogeneous parameters, a local quasi-bursting phase appears at low coupling. The paper's quantitative support is the maximal Lyapunov exponent $\lambda_1$ computed from the full Jacobian: as $N$ increases from 1 to 4, the first chaotic peak shifts rightward and upward, the descent into synchronized bursting becomes more uniform, and the final rise to synchronized hyperchaos near $g=1$ remains sharp. In next-nearest-neighbor lattices the synchronized-hyperchaos rise disappears within $g\in[0,1]$, which the paper attributes to destructive interference from the larger neighborhood. At $g=1$ in large 2D and 3D lattices, the paper identifies checkerboard antiphase spiking of adjacent neurons as extreme one-timestep lag synchronization.

Load-bearing premise

The regime classification and all transition thresholds rest on Lyapunov exponents computed from short finite orbits ($k=2000$, or $10000$ for the two-dimensional nearest-neighbor case) with no convergence curves or error bars, so numerical error could shift the reported regime boundaries.

Editorial extensions

If this is right

  • In NN-coupled homogeneous lattices, the same four regimes appear for $N=1$ through $N=4$, so the regime sequence is a robust feature of the lattice rather than an accident of one dimension.
  • Because the first chaotic peak shifts rightward and upward with $N$, higher-dimensional lattices remain unsynchronized-chaotic over a wider coupling window and reach stronger chaos there; extrapolating the trend predicts sharper transitions in higher dimensions.
  • Next-nearest-neighbor coupling eliminates the synchronized-hyperchaos regime for $g\in[0,1]$; the lattice stays in synchronized chaotic bursting even at $g=1$.
  • Large 2D and 3D lattices show quasi-synchronization and localized synchronization, and at $g=1$ they display extreme one-timestep lag synchronization as checkerboard anti-phase spiking.
  • Heterogeneity can create local quasi-bursting at low coupling, but the phase weakens as dimension increases and appears extinguished for fully heterogeneous lattices in $N=1$ and $N=2$.

Reading between the lines

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

  • If the rightward peak shift is a genuine scaling law, then the transition coupling $g_c(N)$ should either approach a finite limit or follow a dependence like $g_c\sim 1/(2N)$ for large $N$; computing $N=5$ and $N=6$ with longer orbits would test this prediction.
  • The checkerboard one-timestep lag synchronization is a discrete-time analogue of anti-phase synchronization; in continuous-time models it may appear as half-period lag, giving a concrete signature to look for in coupled oscillator experiments or neuromorphic circuits.
  • The local quasi-bursting phase in heterogeneous lattices may be a finite-size or finite-time artifact: its valley moves toward $g=0$ as $N$ grows, so larger systems could show a monotonic chaos onset instead.
  • The paper's 'universal saturation' in strongly coupled heterogeneous lattices suggests the hyperchaotic dynamics may become extensive with dimension; computing the Kaplan-Yorke dimension per neuron across $N$ could reveal a dimension-independent density of positive Lyapunov directions.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. This manuscript studies N-dimensional cubic lattices of electrically coupled nonchaotic Rulkov neurons (N = 1,...,4, plus large 2D and 3D cases) under nearest-neighbor and next-nearest-neighbor coupling, with homogeneous and heterogeneous neuron parameters. For each setting, the authors compute Lyapunov spectra via QR factorization using a piecewise Jacobian derived in Appendix A and classify dynamical regimes as a function of coupling strength g ∈ [0,1]. The main claims are: a universal sequence of regimes—uncoupled nonchaotic spiking, unsynchronized chaotic spiking, synchronized chaotic bursting, synchronized hyperchaos—in homogeneous lattices; an additional "local quasi-bursting" regime in low-conductance heterogeneous lattices; a delay and increase of the first chaos peak under NNN coupling due to "destructive interference"; miniature "phase transitions" in NNN and higher-dimensional lattices; a rightward and upward shift of the first chaotic peak with increasing spatial dimension (Fig. 12); and emergent local, quasi-, and lag synchronization in large lattices, including a one-timestep checkerboard lag-synchronized state at g = 1. The paper includes a tensorial formulation, a Jacobian derivation, and an implementation appendix.

Significance. If confirmed, the results would extend the study of Rulkov map networks from rings to higher-dimensional lattices and would provide a systematic regime catalog for a tractable discrete-time neuronal model. The manuscript has clear strengths: the Jacobian derivation in Appendix A is explicit and correctly reduces to the ring case; the code is publicly available; and the paper makes falsifiable predictions (e.g., peak-shift trends, disappearance of local quasi-bursting at higher dimension). However, the quantitative backbone—finite-time Lyapunov exponents with k = 2000, moving-window smoothing, visual regime assignment, and single realizations for heterogeneous cases—is not yet at the standard needed to support the headline dimensional trends. The contribution is potentially significant but requires substantial numerical validation.

major comments (4)
  1. [§4.1, Fig. 12] The headline dimensional trend—the rightward and upward shift of the first chaotic peak in λ1 with increasing N—rests on Lyapunov exponents computed from orbits of only k = 2000 timesteps (Sec. 3.2, Fig. 12 caption) and on curves smoothed with moving windows of 40–60 points over 1000 values of g. The paper asserts that k = 2000 is long enough "by computer experiment," but no convergence curves or error bars are shown. Because finite-time Lyapunov exponents can misclassify transient contraction as chaos and the moving average can shift peak locations by up to about half the window width (≈0.02–0.03 in g), the peak-shift trend and the associated predictions about higher dimensions are not yet quantitatively supported. Please provide convergence curves for representative g values and smoothing-robustness checks.
  2. [§3.2, Figs. 8 and 10] The miniature "phase transitions" are inferred from small bumps in the λ1 curves that the paper itself states are "on the same order of magnitude as the variance in the Lyapunov exponent" (Sec. 3.2). The supporting visualization (Fig. 10) shows six neurons at six g values, but it does not establish that the bumps are statistically distinct from fluctuations or that they arise from coordinated threshold crossings. Without estimates of the Lyapunov-exponent variance or repeated realizations, this claim is not supported.
  3. [§2.2, Figs. 4, 8, 12] For the heterogeneous cases, a single random realization is used for each parameter distribution (Sec. 2.2). Consequently, statements such as "the local quasi-bursting phase in the N = 2 lattice has been extinguished" in the fully heterogeneous case (Sec. 4.1) or the comparison between partially and fully heterogeneous curves in Figs. 4b and 4c may be realization-specific. Ensemble averages over multiple draws of σi and αi, with error bars, are needed before these contrasts can be taken as properties of the distribution.
  4. [§4.1, Fig. 12a] The regime classification itself is performed by visual inspection of six adjacent neurons and representative snapshots (e.g., Figs. 5–7, 10, 13). Terms such as "local quasi-bursting" and boundaries like 0.2 ≲ g ≲ 0.8 are not tied to quantitative order parameters (e.g., synchronization error, burst length, fraction of silent neurons). Because the paper's central claims concern the existence and ordering of these regimes, I would ask for quantitative definitions and validation, especially for the newly introduced "local quasi-bursting" regime.
minor comments (4)
  1. [Fig. 10 caption] The caption says "three select values of g" but the figure contains six panels; please correct the caption.
  2. [Sec. 2.1] The notation t = k for discrete time is later reused for the Lyapunov orbit length k; this is not an error but could be clarified to avoid confusion in Sec. 3.2.
  3. [Ref. [45]] Reference [45] appears to misstate the title of the McCulloch–Pitts paper; the canonical title is "A logical calculus of the ideas immanent in nervous activity."
  4. [Sec. 4.1] The claim that λ1 "increases monotonically with dimension" around g ≈ 0.2 and g ≳ 0.9 is not obvious from the overlaid curves in Fig. 12a because the N = 2, 3, 4 curves nearly overlap; a zoomed inset or separate panels would help substantiate this statement.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: Lyapunov spectra are computed from the stated lattice equations with no fitted parameters; the only self-citations are non-load-bearing comparisons.

full rationale

The paper's central quantities are produced by direct numerical integration of the model defined in Eqs. (2), (3), (8), and (9). The Jacobian tensor (Eq. (14)) is derived in Appendix A by differentiating those same equations, and the Lyapunov exponents are computed from that Jacobian via QR factorization. There are no free parameters fitted to the regime structure, no regime-defining observable is inserted back into the dynamics, and the 'predictions' for large lattices and higher dimensions are extrapolations tested against new simulations rather than quantities fitted to the data. The regime labels (unsynchronized chaotic spiking, synchronized chaotic bursting, etc.) are qualitative classifications of the observed λ1 curves and voltage traces; classifying data does not make the dynamics circular. The only self-citations are to the author's prior ring-lattice study (Ref. [36]) and to a companion analysis (Ref. [39]), used as baselines or comparisons; none is used to force the main result. The paper's own caveats about finite-time Lyapunov exponents (k = 2000), moving-average smoothing, and bumps 'on the same order of magnitude as the variance in the Lyapunov exponent' are numerical robustness concerns, not circularity. Accordingly, the derivation is self-contained and the circularity score is minimal.

Assumptions & free parameters 8 free parameters · 7 assumptions · 3 invented entities

The central claims depend on standard map-model assumptions plus several computational choices made by the authors. No parameters are fitted to target outcomes, but the scan resolution, orbit length, smoothing window, realization seed, and lattice side length shape the regime boundaries reported, so they are listed as free parameters rather than as fitted constants.

free parameters (8)
  • slow-variable timescale mu = 0.001
    Chosen standard value to make y a slow variable; Eq. (2).
  • current response coefficients beta_c and sigma_c = beta_c=sigma_c=1
    Set for all neurons to incorporate both current effects; Sec. 2.1.
  • coupling strength scan g = 5000 values in [0,1] for small 2D lattices; 1000 values for N-dimensional lattices
    Main control parameter; all regime classifications are read from lambda1(g).
  • neuron parameter ranges sigma_i and alpha_i = sigma=-0.5, alpha=4.5; sigma in U(-1.5,-0.5), alpha in U(4.25,4.75)
    Chosen standard ranges that produce regular spiking, chaotic spiking, or bursting; Sec. 2.2.
  • initial conditions x_i(0) and y_i(0) = x_i(0) in U(-1,1); y_i(0)=-3.25
    Single random draw per lattice, with no specified seed; heterogeneity results depend on this draw.
  • orbit length k for Lyapunov exponents = 10000 for 2D NN; 2000 for NNN and N-dimensional simulations
    Convergence is asserted by computer experiment; shorter orbits introduce more error.
  • moving-average smoothing window = 60 points for N=1; 40 points for N=2,3,4
    Used to produce Fig. 12; the authors state the N=1 divergence is an artifact of this smoothing.
  • lattice side length zeta = 8 and 300 in 2D; 4 in N-dimensional; 50 in large 3D
    Lattice sizes chosen for computational feasibility; small-lattice extrapolation to large lattices relies on locality.
assumptions (7)
  • domain assumption The Rulkov map is a valid phenomenological model of neuronal spiking, bursting, and chaotic behavior at the parameter values used.
    Sec. 2.1 and Fig. 1; the entire study interprets lattice behaviors as neuron dynamics.
  • domain assumption Coupling current into neuron i is well described by the mean voltage difference Ci = g/|N_i| sum (xj - xi), with sigma_c = beta_c = 1.
    Eqs. (4), (8), and (9); this linear diffusive coupling is asserted, not derived from biophysics.
  • domain assumption Periodic boundary conditions and uniform random draws of parameters and initial conditions produce representative dynamics.
    Sec. 2.2; no seed is reported, so the reported curves are single realizations.
  • standard math QR-factorization of the piecewise Jacobian over finite orbits yields Lyapunov exponents that describe the asymptotic dynamics.
    Sec. 3.1 and Appendix A; assumes differentiability except on measure-zero boundaries.
  • ad hoc to paper k=2000 timesteps is sufficient for Lyapunov exponent convergence for NNN and N>2 lattices.
    Sec. 3.2 says by computer experiment without showing convergence data, and also says the reduced orbit length adds error.
  • ad hoc to paper Moving-window smoothing does not change the qualitative trends in lambda1(g).
    Sec. 4.1 invokes a moving average and later states the N=1 divergence is an artifact of the smoothing, so smoothing can affect observed features.
  • ad hoc to paper Plotting six adjacent neurons and hand-picked snapshots reveals the dynamical regime of the whole lattice.
    Secs. 3.1 and 3.3; no synchronization order parameter or clustering metric is used to verify representativeness.
invented entities (3)
  • local quasi-bursting dynamical regime
    purpose: Classify a low-coupling phase in heterogeneous lattices where neighborhoods oscillate between near-spiking and silence without full synchronization.
    Defined by visual inspection of selected time series (Figs. 6b, 7b, 13); no quantitative order parameter or external signature distinguishes it from ordinary uncoupled chaos.
  • destructive interference mechanism
    purpose: Explain the delayed onset of synchronized bursting as more neighbors dilute each neuron's individual influence.
    Presented as a physical explanation in Sec. 3.2, but no direct measurement of interference is given; it is an analogy consistent with the averaging in Eq. (9).
  • miniature phase transitions
    purpose: Describe coordinated threshold crossings that add spikes to silent neurons and produce bumps in lambda1(g).
    A metaphor for stepwise changes in spike count per burst; no thermodynamic order parameter or finite-size scaling is computed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Hyperchaos and complex dynamical regimes in $N$-dimensional neuron lattices." pith.science (2026). https://pith.science/paper/FKHRJR4G

@misc{pith2026250503051,
  author       = {Pith},
  title        = {Pith review of: Hyperchaos and complex dynamical regimes in $N$-dimensional neuron lattices},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FKHRJR4G}},
  note         = {Machine review of arXiv:2505.03051}
}
abstract

We study the dynamics of $N$-dimensional lattices of nonchaotic Rulkov neurons coupled with a flow of electrical current. We consider both nearest-neighbor and next-nearest-neighbor couplings, homogeneous and heterogeneous neurons, and small and large lattices over a wide range of electrical coupling strengths. As the coupling strength is varied, the neurons exhibit a number of complex dynamical regimes, including unsynchronized chaotic spiking, local quasi-bursting, synchronized chaotic bursting, and synchronized hyperchaos. For lattices in higher spatial dimensions, we discover dynamical effects arising from the "destructive interference" of many connected neurons and miniature "phase transitions" from coordinated spiking threshold crossings. In large two- and three-dimensional neuron lattices, we observe emergent dynamics such as local synchronization, quasi-synchronization, and lag synchronization. These results illustrate the rich dynamics that emerge from coupled neurons in multiple spatial dimensions, highlighting how dimensionality, connectivity, and heterogeneity critically shape the collective behavior of neuronal systems.

Figures

Figures reproduced from arXiv: 2505.03051 by the authors.

Figure 1
Figure 1. Examples of the dynamics of individual Rulkov neurons. [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Graphs of the response of the fast variable [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Visualization of a two-dimensional electrically coupled neuron lattice with [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (13 more)
Figure 4
Figure 4. Figure 4: Graphs of the maximal Lyapunov exponent λ1 for 5000 electrical coupling strength values g between 0 and 1 in the homogeneous, partially heterogeneous, and fully heterogeneous cases of an N = 2-dimensional lattice of NN electrically coupled Rulkov neurons with ζ = 8 neu…
Figure 5
Figure 5. Figure 5: Graphs of the voltages x(k) of six adjacent homogeneous NN coupled neurons for four select values of g. The system is an N = 2-dimensional lattice of NN electrically coupled homogeneous Rulkov neurons with ζ = 8 neurons per side. over a range of 5000 electrical couplin…
Figure 6
Figure 6. Figure 6: Graphs of the voltages x(k) of six adjacent partially heterogeneous NN coupled neurons for five select values of g. The system is an N = 2-dimensional lattice of NN electrically coupled partially heterogeneous Rulkov neurons with ζ = 8 neurons per side. electrically se…
Figure 7
Figure 7. Figure 7: Graphs of the voltages x(k) of six adjacent fully heterogeneous NN coupled neurons for five select values of g (the same neuron locations and g values as [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]
Figure 8
Figure 8. Figure 8: Graphs of the maximal Lyapunov exponent λ1 for 5000 electrical coupling strength values g between 0 and 1 in the homogeneous, partially heterogeneous, and fully heterogeneous cases of an N = 2-dimensional lattice of NNN electrically coupled Rulkov neurons with ζ = 8 ne…
Figure 9
Figure 9. Figure 9: Graphs of the voltages x(k) of six adjacent homogeneous NNN coupled neurons for three select values of g. The system is an N = 2-dimensional lattice of NNN electrically coupled homogeneous Rulkov neurons with ζ = 8 neurons per side. (a) g = 0.25, λ1 ≈ 0.104 (b) g = 0.2…
Figure 10
Figure 10. Figure 10: Graphs of the voltages x(k) of six adjacent fully heterogeneous NNN coupled neurons for three select values of g. The system is an N = 2-dimensional lattice of NNN electrically coupled fully heterogeneous Rulkov neurons with ζ = 8 neurons per side. 15 [PITH_FULL_IMAG…
Figure 11
Figure 11. Figure 11: Dynamics of an N = 2-dimensional large lattice of NN electrically coupled homogeneous Rulkov neurons with ζ = 300 neurons per side. The top graphs show the voltages x(k) of six adjacent neurons in the lattice. The bottom graphs show snapshots of the voltages of all th…
Figure 12
Figure 12. Figure 12: Overlaid graphs of the maximal Lyapunov exponent [PITH_FULL_IMAGE:figures/full_fig_p019_12.png]
Figure 13
Figure 13. Figure 13: Graphs of the voltages x(k) of six adjacent partially heterogeneous NN coupled neurons with g = 0.045. These neurons are in the local quasi-bursting regime, which exists in this lattice with ζ = 4 in N = 2, 3, and 4 spatial dimensions. As we have now come to expect, t…
Figure 14
Figure 14. Figure 14: Dynamics of an N = 3-dimensional large lattice of NN electrically coupled homogeneous Rulkov neurons with ζ = 50 neurons per side. The neurons are in the synchronized chaotic bursting regime with electrical coupling strength g = 0.5. Each row displays a snapshot of th…
Figure 15
Figure 15. Figure 15: Dynamics of an N = 3-dimensional large lattice of NN electrically coupled homogeneous Rulkov neurons with ζ = 50 neurons per side. The neurons are in the (quasi-)synchronized hyperchaotic regime with electrical coupling strength g = 1. Each row displays a snapshot of …
Figure 16
Figure 16. Figure 16: Dynamics of an N = 3-dimensional large lattice of NN electrically coupled homogeneous Rulkov neurons with ζ = 50 neurons per side and electrical coupling strength g = 1. These graphs show the existence of extreme lag synchronization in the (quasi-)synchronized hyperch…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

62 extracted references · 60 canonical work pages

  1. [36]

    Chaotic dynamics and fractal geometry in ring lattice systems of nonchaotic Rulkov neurons,

    B. B. Le, “Chaotic dynamics and fractal geometry in ring lattice systems of nonchaotic Rulkov neurons,” arXiv:2412.12134 [nlin.CD] , 2025

    arXiv 2025

  2. [1]

    A review for dynamics in neuron and neuronal network,

    J. Ma and J. Tang, “A review for dynamics in neuron and neuronal network,” Nonlinear Dyanamics , vol. 89, pp. 1569–1578, 2017. 28

    work page 2017

  3. [2]

    A quantitative description of membrane current and its application to conduction and excitation in nerve,

    A. L. Hodgkin and A. F. Huxley, “A quantitative description of membrane current and its application to conduction and excitation in nerve,” The Journal of Physiology , vol. 117, no. 4, pp. 500–544, 1952

    work page 1952

  4. [3]

    A model of neuronal bursting using three coupled first order differential equations,

    J. L. Hindmarsh and R. M. Rose, “A model of neuronal bursting using three coupled first order differential equations,” Proceedings of the Royal Society London Series B Biological Sciences , vol. 221, p. 87–102, 1984

    work page 1984

  5. [4]

    Suppressing spiral wave turbulence in a simple fractional-order discrete neuron map using impulse triggering,

    K. Rajagopal, S. Panahi, M. Chen, S. Jafari, and B. Bao, “Suppressing spiral wave turbulence in a simple fractional-order discrete neuron map using impulse triggering,” Fractals, vol. 29, no. 8, p. 2140030, 2021

    work page 2021

  6. [5]

    Simple model of spiking neurons,

    E. Izhikevich, “Simple model of spiking neurons,” IEEE Transactions on Neural Networks , vol. 14, no. 6, pp. 1569–1572, 2003

    work page 2003

  7. [6]

    Chaos in a three-variable model of an excitable cell,

    T. R. Chay, “Chaos in a three-variable model of an excitable cell,” Physica D: Nonlinear Phenomena , vol. 16, no. 2, pp. 233–242, 1985

    work page 1985

  8. [7]

    Mathematical model of an identified stomatogastric ganglion neuron,

    F. Buchholtz, J. Golowasch, I. R. Epstein, and E. Marder, “Mathematical model of an identified stomatogastric ganglion neuron,” Journal of Neurophysiology , vol. 67, no. 2, pp. 332–340, 1992

    work page 1992

Show all 62 references
  1. [8]

    Can the thermodynamic Hodgkin–Huxley model of voltage-dependent conductance extrapolate for temperature?,

    M. D. Forrest, “Can the thermodynamic Hodgkin–Huxley model of voltage-dependent conductance extrapolate for temperature?,” Computation, vol. 1, p. 47–60, 2014

    2014

  2. [9]

    Poisson-Boltzmann-Nernst-Planck model,

    Q. Zheng and G. W. Wei, “Poisson-Boltzmann-Nernst-Planck model,” Journal of Chemical Physics , vol. 134, no. 19, 2011

    2011

  3. [10]

    Classification of bursting mappings,

    E. M. Izhikevich and F. Hoppensteadt, “Classification of bursting mappings,” International Journal of Bifurca- tion and Chaos , vol. 14, no. 11, pp. 3847–3854, 2004

    2004

  4. [11]

    Generic excitable dynamics on a two-dimensional map,

    D. R. Chialvo, “Generic excitable dynamics on a two-dimensional map,” Chaos, Solitons and Fractals , vol. 5, no. 3-4, pp. 461–479, 1995

    1995

  5. [12]

    On a response characteristic of a mathematical neuron model,

    J. Nagumo and S. Sato, “On a response characteristic of a mathematical neuron model,” Kybernetik, vol. 10, pp. 155–164, 1972

    1972

  6. [13]

    Discrete memristive neuron model and its interspike interval-encoded application in image encryption,

    H. Bao, Z. Hua, W. Liu, and B. Bao, “Discrete memristive neuron model and its interspike interval-encoded application in image encryption,” Science China Technological Sciences , vol. 64, p. 2281–2291, 2021

    2021

  7. [14]

    Map-based models in neuronal dynamics,

    B. Ibarz, J. M. Casado, and M. A. F. Sanju´ an, “Map-based models in neuronal dynamics,” Physics Reports, vol. 501, no. 1-2, pp. 1–74, 2011

    2011

  8. [15]

    Modeling of spiking-bursting neural behavior using two-dimensional map,

    N. F. Rulkov, “Modeling of spiking-bursting neural behavior using two-dimensional map,” Physical Review E , vol. 65, 2002

    2002

  9. [16]

    Regularization of synchronized chaotic bursts,

    N. F. Rulkov, “Regularization of synchronized chaotic bursts,” Physical Review Letters , vol. 86, no. 1, pp. 183– 186, 2001

    2001

  10. [17]

    Stability and chaos of Rulkov map-based neuron network with electrical synapse,

    C. Wang and H. Cao, “Stability and chaos of Rulkov map-based neuron network with electrical synapse,” Communications in Nonlinear Science and Numerical Simulation , vol. 20, no. 2, 2015

    2015

  11. [18]

    Controlling the bursting size in the two-dimensional Rulkov model,

    J. L´ opez, M. Coccolo, R. Cape´ ans, and M. A. Sanju´ an, “Controlling the bursting size in the two-dimensional Rulkov model,” Communications in Nonlinear Science and Numerical Simulation , vol. 120, 2023

    2023

  12. [19]

    Symbolic analysis of bursting dynamical regimes of Rulkov neural networks,

    R. Budzinski, S. Lopes, and C. Masoller, “Symbolic analysis of bursting dynamical regimes of Rulkov neural networks,” Neurocomputing, vol. 441, 2021

    2021

  13. [20]

    Asymmetric coupling of nonchaotic Rulkov neurons: Fractal attractors, quasimultistability, and final state sensitivity,

    B. B. Le, “Asymmetric coupling of nonchaotic Rulkov neurons: Fractal attractors, quasimultistability, and final state sensitivity,” Physical Review E , vol. 111, p. 034201, 2025

    2025

  14. [21]

    Intermittent evolution routes to the periodic or the chaotic orbits in Rulkov map featured,

    P. Ge and H. Cao, “Intermittent evolution routes to the periodic or the chaotic orbits in Rulkov map featured,” Chaos, vol. 31, no. 9, 2021. 29

    2021

  15. [22]

    Hidden and self-excited firing activities of an improved Rulkov neuron, and its application in information patterns,

    Z. T. Njitacke, C. N. Takembo, G. Sani, N. Marwan, R. Yamapi, and J. Awrejcewicz, “Hidden and self-excited firing activities of an improved Rulkov neuron, and its application in information patterns,”Nonlinear Dynamics, vol. 112, 2024

    2024

  16. [23]

    Extreme multi-stability and microchaos of fractional-order memristive Rulkov neuron model considering magnetic induction and its digital watermarking application,

    D. Ding, Y. Niu, Z. Yang, J. Wang, W. Wang, M. Wang, and F. Jin, “Extreme multi-stability and microchaos of fractional-order memristive Rulkov neuron model considering magnetic induction and its digital watermarking application,” Nonlinear Dynamics, vol. 112, 2024

    2024

  17. [24]

    Stability and synchronization of coupled Rulkov map-based neurons with chemical synapses,

    D. Hu and H. Cao, “Stability and synchronization of coupled Rulkov map-based neurons with chemical synapses,” Communications in Nonlinear Science and Numerical Simulation , vol. 35, pp. 105–122, 2016

    2016

  18. [25]

    Complete synchronization of coupled Rulkov neuron networks,

    H. Sun and H. Cao, “Complete synchronization of coupled Rulkov neuron networks,” Nonlinear Dynamics , vol. 84, p. 2423–2434, 2016

    2016

  19. [26]

    Chaotic phase synchronization in scale-free networks of bursting neurons,

    C. A. S. Batista, A. M. Batista, J. A. C. de Pontes, R. L. Viana, and S. R. Lopes, “Chaotic phase synchronization in scale-free networks of bursting neurons,” Phys. Rev. E , vol. 76, p. 016218, 2007

    2007

  20. [27]

    Dynamics of discrete memristor-based Rulkov neuron,

    L. J. Liu, Y. H. Qin, and D. Q. Wei, “Dynamics of discrete memristor-based Rulkov neuron,” IEEE Access, vol. 10, pp. 72051–72056, 2022

    2022

  21. [28]

    Chaos in the Rulkov neuron model based on Marotto’s theorem,

    P. Ge and H. Cao, “Chaos in the Rulkov neuron model based on Marotto’s theorem,” International Journal of Bifurcation and Chaos , vol. 31, no. 15, p. 2150233, 2021

    2021

  22. [29]

    Noise induced complexity: patterns and collective phenomena in a small-world neuronal network,

    Y. Zheng, Q. Wang, and M. F. Danca, “Noise induced complexity: patterns and collective phenomena in a small-world neuronal network,” Cognitive Neurodynamics, vol. 8, p. 143–149, 2014

    2014

  23. [30]

    Neural excitability, spiking and bursting,

    E. M. Izhikevich, “Neural excitability, spiking and bursting,” International Journal of Bifurcation and Chaos , vol. 10, no. 6, pp. 1171–1266, 1999

    1999

  24. [31]

    Which model to use for cortical spiking neurons?,

    E. Izhikevich, “Which model to use for cortical spiking neurons?,” IEEE Transactions on Neural Networks , vol. 15, no. 5, pp. 1063–1070, 2004

    2004

  25. [32]

    Enhancing synchronization in chaotic oscillators by induced heterogeneity,

    R. Banerjee, B. K. Bera, D. Ghosh, and S. K. Dana, “Enhancing synchronization in chaotic oscillators by induced heterogeneity,” The European Physical Journal Special Topics , vol. 226, p. 1893–1902, 2017

    1902

  26. [33]

    Synchronization in a network of model neurons,

    M. P. K. Jampa, A. R. Sonawane, P. M. Gade, and S. Sinha, “Synchronization in a network of model neurons,” Physical Review E , vol. 75, 2007

    2007

  27. [34]

    Enhancement of neuronal coherence by diversity in coupled Rulkov-map models,

    H. Chen, J. Zhang, and J. Liu, “Enhancement of neuronal coherence by diversity in coupled Rulkov-map models,” Physica A: Statistical Mechanics and its Applications , vol. 387, no. 5, pp. 1071–1076, 2008

    2008

  28. [35]

    Synchronization of slow-fast systems,

    I. Omelchenko, M. Rosenblum, and A. Pikovsky, “Synchronization of slow-fast systems,” The European Physical Journal Special Topics , vol. 191, pp. 3–14, 2011

    2011

  29. [37]

    Coexisting attractors and multistate noise-induced intermittency in a cycle ring of Rulkov neurons,

    I. A. Bashkirtseva, A. N. Pisarchik, and L. B. Ryashko, “Coexisting attractors and multistate noise-induced intermittency in a cycle ring of Rulkov neurons,” Mathematics, vol. 11, no. 3, p. 597, 2023

    2023

  30. [38]

    Lag synchronization in an unidirectional ring of memristive neurons,

    V. Vijayan, H. Natiq, S. Momani, V.-T. Pham, and M. Perc, “Lag synchronization in an unidirectional ring of memristive neurons,” The European Physical Journal Special Topics , 2025

    2025

  31. [39]

    Exploring geometrical properties of chaotic systems through an analysis of the Rulkov neuron maps,

    B. B. Le and N. A. Gandhi, “Exploring geometrical properties of chaotic systems through an analysis of the Rulkov neuron maps,” arXiv:2406.08385 [nlin.CD] , 2024

    2024 arXiv

  32. [40]

    On next-nearest-neighbor interaction in linear chain. I,

    C. K. Majumdar and D. K. Ghosh, “On next-nearest-neighbor interaction in linear chain. I,” Journal of Math- ematical Physics, vol. 10, no. 8, pp. 1388–1398, 1969

    1969

  33. [41]

    Phase diagrams and critical behavior in Ising square lattices with nearest- and next-nearest-neighbor interactions,

    K. Binder and D. P. Landau, “Phase diagrams and critical behavior in Ising square lattices with nearest- and next-nearest-neighbor interactions,” Physical Review B , vol. 21, pp. 1941–1962, 1980. 30

    1941

  34. [42]

    Synchronization of pulse-coupled biological oscillators,

    R. E. Mirollo and S. H. Strogatz, “Synchronization of pulse-coupled biological oscillators,” SIAM Journal on Applied Mathematics, vol. 50, no. 6, pp. 1645–1662, 1990

    1990

  35. [43]

    Bursting synchronization in non-locally coupled maps,

    J. de Pontes, R. Viana, S. Lopes, C. Batista, and A. Batista, “Bursting synchronization in non-locally coupled maps,” Physica A: Statistical Mechanics and its Applications , vol. 387, no. 16, 2008

    2008

  36. [44]

    Ergodic theory of chaos and strange attractors,

    J.-P. Eckmann and D. Ruelle, “Ergodic theory of chaos and strange attractors,” Reviews of Modern Physics , vol. 57, no. 3, pp. 617–656, 1985

    1985

  37. [45]

    Neural networks and physical systems with emergent collective computational abilities,

    W. S. McCulloch and W. Pitts, “Neural networks and physical systems with emergent collective computational abilities,” Bulletin of Mathematical Biophysics , vol. 5, p. 115–133, 1943

    1943

  38. [46]

    Neural networks and physical systems with emergent collective computational abilities,

    J. J. Hopfield, “Neural networks and physical systems with emergent collective computational abilities,” Pro- ceedings of the National Academy of Sciences , vol. 92, pp. 6655–6662, 1982

    1982

  39. [47]

    Cluster synchronization in networks of coupled nonidentical dynamical systems,

    W. Lu, B. Liu, and T. Chen, “Cluster synchronization in networks of coupled nonidentical dynamical systems,” Chaos, vol. 20, no. 1, p. 013120, 2010

    2010

  40. [48]

    Cluster synchronization in oscillatory networks,

    V. N. Belykh, G. V. Osipov, V. S. Petrov, J. A. K. Suykens, and J. Vandewalle, “Cluster synchronization in oscillatory networks,” Chaos, vol. 18, no. 3, p. 037106, 2008

    2008

  41. [49]

    Cluster synchronization in an ensemble of neurons interacting through chemical synapses,

    M. Yoshioka, “Cluster synchronization in an ensemble of neurons interacting through chemical synapses,” Phys- ical Review E , vol. 71, p. 061914, 2005

    2005

  42. [50]

    Quasi-synchronization of heterogeneous dynamic networks via distributed impulsive control: Error estimation, optimization and design,

    W. He, F. Qian, J. Lam, G. Chen, Q.-L. Han, and J. Kurths, “Quasi-synchronization of heterogeneous dynamic networks via distributed impulsive control: Error estimation, optimization and design,” Automatica, vol. 62, pp. 249–262, 2015

    2015

  43. [51]

    Lag quasi-synchronization of coupled delayed systems with parameter mismatch,

    W. He, F. Qian, Q.-L. Han, and J. Cao, “Lag quasi-synchronization of coupled delayed systems with parameter mismatch,” IEEE Transactions on Circuits and Systems I: Regular Papers , vol. 58, no. 6, pp. 1345–1357, 2011

    2011

  44. [52]

    Emergence of local synchronization in neuronal networks with adaptive couplings,

    S. Chakravartula, P. Indic, B. Sundaram, and T. Killingback, “Emergence of local synchronization in neuronal networks with adaptive couplings,” PLOS One , vol. 12, no. 6, 2017

    2017

  45. [53]

    Local synchronization in complex networks of coupled oscillators,

    J. Stout, M. Whiteway, E. Ott, M. Girvan, and T. M. Antonsen, “Local synchronization in complex networks of coupled oscillators,” Chaos, vol. 21, no. 2, p. 025109, 2011

    2011

  46. [54]

    Rapid local synchronization of action potentials: toward computation with coupled integrate-and-fire neurons,

    J. J. Hopfield and A. V. Herz, “Rapid local synchronization of action potentials: toward computation with coupled integrate-and-fire neurons,” Proceedings of the National Academy of Sciences , vol. 92, no. 15, pp. 6655– 6662, 1995

    1995

  47. [55]

    Lag synchronization in time-delayed systems,

    E. M. Shahverdiev, S. Sivaprakasam, and K. A. Shore, “Lag synchronization in time-delayed systems,” Physics Letters A, vol. 292, no. 6, pp. 320–324, 2002

    2002

  48. [56]

    From phase to lag synchronization in coupled chaotic oscillators,

    M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, “From phase to lag synchronization in coupled chaotic oscillators,” Physical Review Letters , vol. 78, pp. 4193–4196, 1997

    1997

  49. [57]

    Lag synchronization of hyperchaos with application to secure communica- tions,

    C. Li, X. Liao, and K.-W. Wong, “Lag synchronization of hyperchaos with application to secure communica- tions,” Chaos, Solitons and Fractals , vol. 23, no. 1, pp. 183–193, 2005

    2005

  50. [58]

    Lag synchronization of hyperchaotic complex nonlinear systems,

    G. M. Mahmoud and E. E. Mahmoud, “Lag synchronization of hyperchaotic complex nonlinear systems,” Nonlinear Dynamics, vol. 67, p. 1613–1622, 2012

    2012

  51. [59]

    Neural synchrony in brain disorders: Relevance for cognitive dysfunctions and pathophysiology,

    P. J. Uhlhaas and W. Singer, “Neural synchrony in brain disorders: Relevance for cognitive dysfunctions and pathophysiology,” Neuron, vol. 52, no. 1, pp. 155–168, 2006

    2006

  52. [60]

    On the nature of seizure dynamics,

    V. K. Jirsa, W. C. Stacey, P. P. Quilichini, A. I. Ivanov, and C. Bernard, “On the nature of seizure dynamics,” Brain, vol. 137, no. 8, pp. 2210–2230, 2014

    2014

  53. [61]

    Memory and information processing in neuromorphic systems,

    G. Indiveri and S.-C. Liu, “Memory and information processing in neuromorphic systems,” Proceedings of the IEEE, vol. 103, no. 8, pp. 1379–1397, 2015

    2015

  54. [62]

    https://github.com/brandon-bd-le/ND-Lattices . 31

Pith tools

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