Chialvo map

Last updated
Activation variable as a function of time for the chaotic regime Chaotic act variable.png
Activation variable as a function of time for the chaotic regime
Solution for the Chialvo map equations for the chaotic regime Chaotic regime.png
Solution for the Chialvo map equations for the chaotic regime
Activation variable as a function of time for the excitable regime Excitable act variable.png
Activation variable as a function of time for the excitable regime
Solution for the Chialvo map equations for the excitable regime Excitable regime.png
Solution for the Chialvo map equations for the excitable regime

The Chialvo map is a two-dimensional map proposed by Dante R. Chialvo in 1995 [1] to describe the generic dynamics of excitable systems. The model is inspired by Kunihiko Kaneko's Coupled map lattice [2] (CML) numerical approach which considers time and space as discrete variables but state as a continuous one. Later on Rulkov popularized a similar approach. [3] By using only three parameters the model is able to efficiently mimic generic neuronal dynamics in computational simulations, as single elements or as parts of inter-connected networks.

Contents

The model

The model is an iterative map where at each time step, the behavior of one neuron is updated as the following equations:

in which, is called activation or action potential variable, and is the recovery variable. The model has four parameters, is a time-dependent additive perturbation or a constant bias, is the time constant of recovery , is the activation-dependence of the recovery process and is an offset constant. The model has a rich dynamics, presenting from oscillatory  to chaotic behavior, [4] [5] as well as non trivial responses to small stochastic fluctuations. [6] [7]

Analysis

Bursting and chaos

The map is able to capture the aperiodic solutions and the bursting behavior which are remarkable in the context of neural systems. For example, for the values , and and changing b from to the system passes from oscillations to aperiodic bursting solutions.

Fixed points

Considering the case where and the model mimics the lack of ‘voltage-dependence inactivation’ for real neurons and the evolution of the recovery variable is fixed at . Therefore, the dynamics of the activation variable is basically described by the iteration of the following equations

in which as a function of has a period-doubling bifurcation structure.

Examples

Example 1

A practical implementation is the combination of neurons over a lattice, for that, it can be defined as a coupling constant for combining the neurons. For neurons in a single row, we can define the evolution of action potential on time by the diffusion of the local temperature in:

where is the time step and is the index of each neuron. For the values , , and , in absence of perturbations they are at the resting state. If we introduce a stimulus over cell 1, it induces two propagated waves circulating in opposite directions that eventually collapse and die in the middle of the ring.

Example 2

Analogous to the previous example, it's possible create a set of coupling neurons over a 2-D lattice, in this case the evolution of action potentials is given by:

where , , represent the index of each neuron in a square lattice of size , . With this example spiral waves can be represented for specific values of parameters. In order to visualize the spirals, we set the initial condition in a specific configuration and the recovery as .

Example of spiral waves for the Two-dimensional Chialvo map in 100 x 100 lattice and parameters a=0.89, b=0.6, c= 0.26, and k=0.02. SpiralChialvomap.gif
Example of spiral waves for the Two-dimensional Chialvo map in 100 x 100 lattice and parameters a=0.89, b=0.6, c= 0.26, and k=0.02.

The map can also present chaotic dynamics for certain parameter values. In the following figure we show the chaotic behavior of the variable on a square network of for the parameters , , and .

Evolution of Potential X as a function of time in a 500x500 lattice for a chaotic regime with parameters a=0.89, b=0.18, c= 0.28, and k=0.026. Evolution of Potential X as a function of time in a 500x500 lattice for a chaotic regime.gif
Evolution of Potential X as a function of time in a 500x500 lattice for a chaotic regime with parameters a=0.89, b=0.18, c= 0.28, and k=0.026.

The map can be used to simulated a nonquenched disordered lattice (as in Ref [8] ), where each map connects with four nearest neighbors on a square lattice, and in addition each map has a probability of connecting to another one randomly chosen, multiple coexisting circular excitation waves will emerge at the beginning of the simulation until spirals takes over.

Example of spiral waves for the Two-dimensional Chialvo map in an annealed random network starting from a 128 x 128 lattice with probability of rewiring
p
=
0.25
{\displaystyle p=0.25}
and parameters a=0.89, b=0.6, c= 0.26, and k=0.02. ESPIRALES2.gif
Example of spiral waves for the Two-dimensional Chialvo map in an annealed random network starting from a 128 x 128 lattice with probability of rewiring and parameters a=0.89, b=0.6, c= 0.26, and k=0.02.

Chaotic and periodic behavior for a neuron

For a neuron, in the limit of , the map becomes 1D, since converges to a constant. If the parameter is scanned in a range, different orbits will be seen, some periodic, others chaotic, that appear between two fixed points, one at  ; and the other close to the value of (which would be the regime excitable).

Evolution of
x
{\displaystyle x}
as a function of parameter
b
{\displaystyle b}
for a Chialvo map neuron. Parameters:
a
=
0.89
{\displaystyle a=0.89}
,
c
=
0.28
{\displaystyle c=0.28}
,
k
=
0.026
{\displaystyle k=0.026}
, and
b
{\displaystyle b}
from
0.16
{\displaystyle 0.16}
to
0.4
{\displaystyle 0.4}
. Chialvo map for an Neuron.png
Evolution of as a function of parameter for a Chialvo map neuron. Parameters: , , , and from to .
GIF: Evolution of
x
{\displaystyle x}
as a function of parameter
b
{\displaystyle b}
for a Chialvo map neuron. Parameters:
a
=
0.89
{\displaystyle a=0.89}
,
c
=
0.28
{\displaystyle c=0.28}
,
k
=
0.026
{\displaystyle k=0.026}
, and
b
{\displaystyle b}
from
0.16
{\displaystyle 0.16}
to
0.4
{\displaystyle 0.4}
. Chialvomal aneuron.gif
GIF: Evolution of as a function of parameter for a Chialvo map neuron. Parameters: , , , and from to .

Related Research Articles

The Ising model, named after the physicists Ernst Ising and Wilhelm Lenz, is a mathematical model of ferromagnetism in statistical mechanics. The model consists of discrete variables that represent magnetic dipole moments of atomic "spins" that can be in one of two states. The spins are arranged in a graph, usually a lattice, allowing each spin to interact with its neighbors. Neighboring spins that agree have a lower energy than those that disagree; the system tends to the lowest energy but heat disturbs this tendency, thus creating the possibility of different structural phases. The model allows the identification of phase transitions as a simplified model of reality. The two-dimensional square-lattice Ising model is one of the simplest statistical models to show a phase transition.

In probability theory and mathematical physics, a random matrix is a matrix-valued random variable—that is, a matrix in which some or all elements are random variables. Many important properties of physical systems can be represented mathematically as matrix problems. For example, the thermal conductivity of a lattice can be computed from the dynamical matrix of the particle-particle interactions within the lattice.

<span class="mw-page-title-main">Hodgkin–Huxley model</span> Describes how neurons transmit electric signals

The Hodgkin–Huxley model, or conductance-based model, is a mathematical model that describes how action potentials in neurons are initiated and propagated. It is a set of nonlinear differential equations that approximates the electrical engineering characteristics of excitable cells such as neurons and muscle cells. It is a continuous-time dynamical system.

<span class="mw-page-title-main">Lorenz system</span> System of ordinary differential equations with chaotic solutions

The Lorenz system is a system of ordinary differential equations first studied by mathematician and meteorologist Edward Lorenz. It is notable for having chaotic solutions for certain parameter values and initial conditions. In particular, the Lorenz attractor is a set of chaotic solutions of the Lorenz system. In popular media the "butterfly effect" stems from the real-world implications of the Lorenz attractor, namely that several different initial chaotic conditions evolve in phase space in a way that never repeats, so all chaos is unpredictable. This underscores that chaotic systems can be completely deterministic and yet still be inherently unpredictable over long periods of time. Because chaos continually increases in systems, we cannot predict the future of systems well. E.g., even the small flap of a butterfly’s wings could set the world on a vastly different trajectory, such as by causing a hurricane. The shape of the Lorenz attractor itself, when plotted in phase space, may also be seen to resemble a butterfly.

In the field of mathematical modeling, a radial basis function network is an artificial neural network that uses radial basis functions as activation functions. The output of the network is a linear combination of radial basis functions of the inputs and neuron parameters. Radial basis function networks have many uses, including function approximation, time series prediction, classification, and system control. They were first formulated in a 1988 paper by Broomhead and Lowe, both researchers at the Royal Signals and Radar Establishment.

Synchronization of chaos is a phenomenon that may occur when two or more dissipative chaotic systems are coupled.

In computational neuroscience, the Wilson–Cowan model describes the dynamics of interactions between populations of very simple excitatory and inhibitory model neurons. It was developed by Hugh R. Wilson and Jack D. Cowan and extensions of the model have been widely used in modeling neuronal populations. The model is important historically because it uses phase plane methods and numerical solutions to describe the responses of neuronal populations to stimuli. Because the model neurons are simple, only elementary limit cycle behavior, i.e. neural oscillations, and stimulus-dependent evoked responses are predicted. The key findings include the existence of multiple stable states, and hysteresis, in the population response.

<span class="mw-page-title-main">Biological neuron model</span> Mathematical descriptions of the properties of certain cells in the nervous system

Biological neuron models, also known as a spiking neuron models, are mathematical descriptions of neurons. In particular, these models describe how the voltage potential across the cell membrane changes over time. In an experimental setting, stimulating neurons with an electrical current generates an action potential, that propagates down the neuron's axon. This axon can branch out and connect to a large number of downstream neurons at sites called synapses. At these synapses, the spike can cause release of a biochemical substance (neurotransmitter), which in turn can change the voltage potential of downstream neurons, potentially leading to spikes in those downstream neurons, thus propagating the signal. As many as 85% of neurons in the neocortex, the outermost layer of the mammalian brain, consists of excitatory pyramidal neurons, and each pyramidal neuron receives tens of thousands of inputs from other neurons. Thus, spiking neurons are a major information processing unit of the nervous system.

<span class="mw-page-title-main">Light front quantization</span> Technique in computational quantum field theory

The light-front quantization of quantum field theories provides a useful alternative to ordinary equal-time quantization. In particular, it can lead to a relativistic description of bound systems in terms of quantum-mechanical wave functions. The quantization is based on the choice of light-front coordinates, where plays the role of time and the corresponding spatial coordinate is . Here, is the ordinary time, is one Cartesian coordinate, and is the speed of light. The other two Cartesian coordinates, and , are untouched and often called transverse or perpendicular, denoted by symbols of the type . The choice of the frame of reference where the time and -axis are defined can be left unspecified in an exactly soluble relativistic theory, but in practical calculations some choices may be more suitable than others.

<span class="mw-page-title-main">Hindmarsh–Rose model</span> Of the spiking-bursting behavior of a neuron

The Hindmarsh–Rose model of neuronal activity is aimed to study the spiking-bursting behavior of the membrane potential observed in experiments made with a single neuron. The relevant variable is the membrane potential, x(t), which is written in dimensionless units. There are two more variables, y(t) and z(t), which take into account the transport of ions across the membrane through the ion channels. The transport of sodium and potassium ions is made through fast ion channels and its rate is measured by y(t), which is called the spiking variable. z(t) corresponds to an adaptation current, which is incremented at every spike, leading to a decrease in the firing rate. Then, the Hindmarsh–Rose model has the mathematical form of a system of three nonlinear ordinary differential equations on the dimensionless dynamical variables x(t), y(t), and z(t). They read:

A coupled map lattice (CML) is a dynamical system that models the behavior of nonlinear systems. They are predominantly used to qualitatively study the chaotic dynamics of spatially extended systems. This includes the dynamics of spatiotemporal chaos where the number of effective degrees of freedom diverges as the size of the system increases.

The dynamical systems approach to neuroscience is a branch of mathematical biology that utilizes nonlinear dynamics to understand and model the nervous system and its functions. In a dynamical system, all possible states are expressed by a phase space. Such systems can experience bifurcation as a function of its bifurcation parameters and often exhibit chaos. Dynamical neuroscience describes the non-linear dynamics at many levels of the brain from single neural cells to cognitive processes, sleep states and the behavior of neurons in large-scale neuronal simulation.

Iterated filtering algorithms are a tool for maximum likelihood inference on partially observed dynamical systems. Stochastic perturbations to the unknown parameters are used to explore the parameter space. Applying sequential Monte Carlo to this extended model results in the selection of the parameter values that are more consistent with the data. Appropriately constructed procedures, iterating with successively diminished perturbations, converge to the maximum likelihood estimate. Iterated filtering methods have so far been used most extensively to study infectious disease transmission dynamics. Case studies include cholera, Ebola virus, influenza, malaria, HIV, pertussis, poliovirus and measles. Other areas which have been proposed to be suitable for these methods include ecological dynamics and finance.

<span class="mw-page-title-main">Theta model</span>

The theta model, or Ermentrout–Kopell canonical model, is a biological neuron model originally developed to mathematically describe neurons in the animal Aplysia. The model is particularly well-suited to describe neural bursting, which is characterized by periodic transitions between rapid oscillations in the membrane potential followed by quiescence. This bursting behavior is often found in neurons responsible for controlling and maintaining steady rhythms such as breathing, swimming, and digesting. Of the three main classes of bursting neurons, the theta model describes parabolic bursting, which is characterized by a parabolic frequency curve during each burst.

<span class="mw-page-title-main">Rulkov map</span> Map used to model a biological neuron

The Rulkov map is a two-dimensional iterated map used to model a biological neuron. It was proposed by Nikolai F. Rulkov in 2001. The use of this map to study neural networks has computational advantages because the map is easier to iterate than a continuous dynamical system. This saves memory and simplifies the computation of large neural networks.

In the fields of dynamical systems and control theory, a fractional-order system is a dynamical system that can be modeled by a fractional differential equation containing derivatives of non-integer order. Such systems are said to have fractional dynamics. Derivatives and integrals of fractional orders are used to describe objects that can be characterized by power-law nonlocality, power-law long-range dependence or fractal properties. Fractional-order systems are useful in studying the anomalous behavior of dynamical systems in physics, electrochemistry, biology, viscoelasticity and chaotic systems.

The spike response model (SRM) is a spiking neuron model in which spikes are generated by either a deterministic or a stochastic threshold process. In the SRM, the membrane voltage V is described as a linear sum of the postsynaptic potentials (PSPs) caused by spike arrivals to which the effects of refractoriness and adaptation are added. The threshold is either fixed or dynamic. In the latter case it increases after each spike. The SRM is flexible enough to account for a variety of neuronal firing pattern in response to step current input. The SRM has also been used in the theory of computation to quantify the capacity of spiking neural networks; and in the neurosciences to predict the subthreshold voltage and the firing times of cortical neurons during stimulation with a time-dependent current stimulation. The name Spike Response Model points to the property that the two important filters and of the model can be interpreted as the response of the membrane potential to an incoming spike (response kernel , the PSP) and to an outgoing spike (response kernel , also called refractory kernel). The SRM has been formulated in continuous time and in discrete time. The SRM can be viewed as a generalized linear model (GLM) or as an (integrated version of) a generalized integrate-and-fire model with adaptation.

Tau functions are an important ingredient in the modern mathematical theory of integrable systems, and have numerous applications in a variety of other domains. They were originally introduced by Ryogo Hirota in his direct method approach to soliton equations, based on expressing them in an equivalent bilinear form.

In computational and mathematical biology, a biological lattice-gas cellular automaton (BIO-LGCA) is a discrete model for moving and interacting biological agents, a type of cellular automaton. The BIO-LGCA is based on the lattice-gas cellular automaton (LGCA) model used in fluid dynamics. A BIO-LGCA model describes cells and other motile biological agents as point particles moving on a discrete lattice, thereby interacting with nearby particles. Contrary to classic cellular automaton models, particles in BIO-LGCA are defined by their position and velocity. This allows to model and analyze active fluids and collective migration mediated primarily through changes in momentum, rather than density. BIO-LGCA applications include cancer invasion and cancer progression.

<span class="mw-page-title-main">Synthetic nervous system</span> Computational neuroscience model

Synthetic Nervous System (SNS) is a computational neuroscience model that may be developed with the Functional Subnetwork Approach (FSA) to create biologically plausible models of circuits in a nervous system. The FSA enables the direct analytical tuning of dynamical networks that perform specific operations within the nervous system without the need for global optimization methods like genetic algorithms and reinforcement learning. The primary use case for a SNS is system control, where the system is most often a simulated biomechanical model or a physical robotic platform. An SNS is a form of a neural network much like artificial neural networks (ANNs), convolutional neural networks (CNN), and recurrent neural networks (RNN). The building blocks for each of these neural networks is a series of nodes and connections denoted as neurons and synapses. More conventional artificial neural networks rely on training phases where they use large data sets to form correlations and thus “learn” to identify a given object or pattern. When done properly this training results in systems that can produce a desired result, sometimes with impressive accuracy. However, the systems themselves are typically “black boxes” meaning there is no readily distinguishable mapping between structure and function of the network. This makes it difficult to alter the function, without simply starting over, or extract biological meaning except in specialized cases. The SNS method differentiates itself by using details of both structure and function of biological nervous systems. The neurons and synapse connections are intentionally designed rather than iteratively changed as part of a learning algorithm.

References

  1. Chialvo, Dante R. (1995-03-01). "Generic excitable dynamics on a two-dimensional map". Chaos, Solitons & Fractals. Nonlinear Phenomena in Excitable Physiological Systems. 5 (3): 461–479. Bibcode:1995CSF.....5..461C. doi:10.1016/0960-0779(93)E0056-H. ISSN   0960-0779.
  2. "Coupled map lattice", Wikipedia, 2022-07-23, retrieved 2022-09-11
  3. Rulkov, Nikolai F. (2002-04-10). "Modeling of spiking-bursting neural behavior using two-dimensional map". Physical Review E. 65 (4): 041922. arXiv: nlin/0201006 . Bibcode:2002PhRvE..65d1922R. doi:10.1103/PhysRevE.65.041922. PMID   12005888. S2CID   1998912.
  4. Pilarczyk, Paweł; Signerska-Rynkowska, Justyna; Graff, Grzegorz (2022-09-07). "Topological-numerical analysis of a two-dimensional discrete neuron model". arXiv: 2209.03443 [math.DS].
  5. Wang, Fengjuan; Cao, Hongjun (2018-03-01). "Mode locking and quasiperiodicity in a discrete-time Chialvo neuron model". Communications in Nonlinear Science and Numerical Simulation. 56: 481–489. Bibcode:2018CNSNS..56..481W. doi:10.1016/j.cnsns.2017.08.027. ISSN   1007-5704.
  6. Chialvo, Dante R.; Apkarian, A. Vania (1993-01-01). "Modulated noisy biological dynamics: Three examples". Journal of Statistical Physics. 70 (1): 375–391. Bibcode:1993JSP....70..375C. doi:10.1007/BF01053974. ISSN   1572-9613. S2CID   121830779.
  7. Bashkirtseva, Irina; Ryashko, Lev; Used, Javier; Seoane, Jesús M.; Sanjuán, Miguel A. F. (2023-01-01). "Noise-induced complex dynamics and synchronization in the map-based Chialvo neuron model". Communications in Nonlinear Science and Numerical Simulation. 116: 106867. Bibcode:2023CNSNS.11606867B. doi: 10.1016/j.cnsns.2022.106867 . ISSN   1007-5704. S2CID   252140483.
  8. Sinha, Sitabhra; Saramäki, Jari; Kaski, Kimmo (2007-07-09). "Emergence of self-sustained patterns in small-world excitable media". Physical Review E. 76 (1): 015101. arXiv: cond-mat/0701121 . Bibcode:2007PhRvE..76a5101S. doi:10.1103/PhysRevE.76.015101. ISSN   1539-3755. PMID   17677522. S2CID   11714109.