Activating function

Last updated

The activating function is a mathematical formalism that is used to approximate the influence of an extracellular field on an axon or neurons. [1] [2] [3] [4] [5] [6] It was developed by Frank Rattay and is a useful tool to approximate the influence of functional electrical stimulation (FES) or neuromodulation techniques on target neurons. [7] It points out locations of high hyperpolarization and depolarization caused by the electrical field acting upon the nerve fiber. As a rule of thumb, the activating function is proportional to the second-order spatial derivative of the extracellular potential along the axon.

Contents

Equations

In a compartment model of an axon, the activating function of compartment n, , is derived from the driving term of the external potential, or the equivalent injected current

,

where is the membrane capacity, the extracellular voltage outside compartment relative to the ground and the axonal resistance of compartment .

The activating function represents the rate of membrane potential change if the neuron is in resting state before the stimulation. Its physical dimensions are V/s or mV/ms. In other words, it represents the slope of the membrane voltage at the beginning of the stimulation. [8]

Following McNeal's [9] simplifications for long fibers of an ideal internode membrane, with both membrane capacity and conductance assumed to be 0 the differential equation determining the membrane potential for each node is:

,

where is the constant fiber diameter, the node-to-node distance, the node length the axomplasmatic resistivity, the capacity and the ionic currents. From this the activating function follows as:

.

In this case the activating function is proportional to the second order spatial difference of the extracellular potential along the fibers. If and then:

.

Thus is proportional to the second order spatial differential along the fiber.

Interpretation

Positive values of suggest a depolarization of the membrane potential and negative values a hyperpolarization of the membrane potential.

Related Research Articles

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols , (where is the nabla operator), or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

<span class="mw-page-title-main">Chemiosmosis</span> Electrochemical principle that enables cellular respiration

Chemiosmosis is the movement of ions across a semipermeable membrane bound structure, down their electrochemical gradient. An important example is the formation of adenosine triphosphate (ATP) by the movement of hydrogen ions (H+) across a membrane during cellular respiration or photosynthesis.

<span class="mw-page-title-main">Backpropagation</span> Optimization algorithm for artificial neural networks

In machine learning, backpropagation is a widely used algorithm for training feedforward artificial neural networks or other parameterized networks with differentiable nodes. It is an efficient application of the Leibniz chain rule (1673) to such networks. It is also known as the reverse mode of automatic differentiation or reverse accumulation, due to Seppo Linnainmaa (1970). The term "back-propagating errors" was introduced in 1962 by Frank Rosenblatt, but he did not know how to implement this, although Henry J. Kelley had a continuous precursor of backpropagation already in 1960 in the context of control theory.

<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 characteristics of excitable cells such as neurons and muscle cells. It is a continuous-time dynamical system.

<span class="mw-page-title-main">Wavelet transform</span> Mathematical technique used in data compression and analysis

In mathematics, a wavelet series is a representation of a square-integrable function by a certain orthonormal series generated by a wavelet. This article provides a formal, mathematical definition of an orthonormal wavelet and of the integral wavelet transform.

In neuroscience, single-unit recordings provide a method of measuring the electro-physiological responses of a single neuron using a microelectrode system. When a neuron generates an action potential, the signal propagates down the neuron as a current which flows in and out of the cell through excitable membrane regions in the soma and axon. A microelectrode is inserted into the brain, where it can record the rate of change in voltage with respect to time. These microelectrodes must be fine-tipped, high-impedance conductors; they are primarily glass micro-pipettes, metal microelectrodes made of platinum, tungsten, iridium or even iridium oxide. Microelectrodes can be carefully placed close to the cell membrane, allowing the ability to record extracellularly.

In neurobiology, the length constant (λ) is a mathematical constant used to quantify the distance that a graded electric potential will travel along a neurite via passive electrical conduction. The greater the value of the length constant, the farther the potential will travel. A large length constant can contribute to spatial summation—the electrical addition of one potential with potentials from adjacent areas of the cell.

In information theory, information dimension is an information measure for random vectors in Euclidean space, based on the normalized entropy of finely quantized versions of the random vectors. This concept was first introduced by Alfréd Rényi in 1959.

<span class="mw-page-title-main">Chronaxie</span> Electrophysiology metric

Chronaxie is the minimum time required for an electric current double the strength of the rheobase to stimulate a muscle or a neuron. Rheobase is the lowest intensity with indefinite pulse duration which just stimulated muscles or nerves. Chronaxie is dependent on the density of voltage-gated sodium channels in the cell, which affect that cell’s excitability. Chronaxie varies across different types of tissue: fast-twitch muscles have a lower chronaxie, slow-twitch muscles have a higher one. Chronaxie is the tissue-excitability parameter that permits choice of the optimum stimulus pulse duration for stimulation of any excitable tissue. Chronaxie (c) is the Lapicque descriptor of the stimulus pulse duration for a current of twice rheobasic (b) strength, which is the threshold current for an infinitely long-duration stimulus pulse. Lapicque showed that these two quantities (c,b) define the strength-duration curve for current: I = b(1+c/d), where d is the pulse duration. However, there are two other electrical parameters used to describe a stimulus: energy and charge. The minimum energy occurs with a pulse duration equal to chronaxie. Minimum charge (bc) occurs with an infinitely short-duration pulse. Choice of a pulse duration equal to 10c requires a current of only 10% above rheobase (b). Choice of a pulse duration of 0.1c requires a charge of 10% above the minimum charge (bc).

<span class="mw-page-title-main">Cable theory</span> Mathematical model of a dendrite

Classical cable theory uses mathematical models to calculate the electric current along passive neurites, particularly the dendrites that receive synaptic inputs at different sites and times. Estimates are made by modeling dendrites and axons as cylinders composed of segments with capacitances and resistances combined in parallel. The capacitance of a neuronal fiber comes about because electrostatic forces are acting through the very thin lipid bilayer. The resistance in series along the fiber is due to the axoplasm's significant resistance to movement of electric charge.

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.

<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 the properties of certain cells in the nervous system that generate sharp electrical potentials across their cell membrane, roughly one millisecond in duration, called action potentials or spikes. Since spikes are transmitted along the axon and synapses from the sending neuron to many other neurons, spiking neurons are considered to be a major information processing unit of the nervous system. Spiking neuron models can be divided into different categories: the most detailed mathematical models are biophysical neuron models that describe the membrane voltage as a function of the input current and the activation of ion channels. Mathematically simpler are integrate-and-fire models that describe the membrane voltage as a function of the input current and predict the spike times without a description of the biophysical processes that shape the time course of an action potential. Even more abstract models only predict output spikes as a function of the stimulation where the stimulation can occur through sensory input or pharmacologically. This article provides a short overview of different spiking neuron models and links, whenever possible to experimental phenomena. It includes deterministic and probabilistic models.

In neurophysiology, several mathematical models of the action potential have been developed, which fall into two basic types. The first type seeks to model the experimental data quantitatively, i.e., to reproduce the measurements of current and voltage exactly. The renowned Hodgkin–Huxley model of the axon from the Loligo squid exemplifies such models. Although qualitatively correct, the H-H model does not describe every type of excitable membrane accurately, since it considers only two ions, each with only one type of voltage-sensitive channel. However, other ions such as calcium may be important and there is a great diversity of channels for all ions. As an example, the cardiac action potential illustrates how differently shaped action potentials can be generated on membranes with voltage-sensitive calcium channels and different types of sodium/potassium channels. The second type of mathematical model is a simplification of the first type; the goal is not to reproduce the experimental data, but to understand qualitatively the role of action potentials in neural circuits. For such a purpose, detailed physiological models may be unnecessarily complicated and may obscure the "forest for the trees". The FitzHugh–Nagumo model is typical of this class, which is often studied for its entrainment behavior. Entrainment is commonly observed in nature, for example in the synchronized lighting of fireflies, which is coordinated by a burst of action potentials; entrainment can also be observed in individual neurons. Both types of models may be used to understand the behavior of small biological neural networks, such as the central pattern generators responsible for some automatic reflex actions. Such networks can generate a complex temporal pattern of action potentials that is used to coordinate muscular contractions, such as those involved in breathing or fast swimming to escape a predator.

In physics and engineering, the time constant, usually denoted by the Greek letter τ (tau), is the parameter characterizing the response to a step input of a first-order, linear time-invariant (LTI) system. The time constant is the main characteristic unit of a first-order LTI system.

The bidomain model is a mathematical model to define the electrical activity of the heart. It consists in a continuum (volume-average) approach in which the cardiac microstructure is defined in terms of muscle fibers grouped in sheets, creating a complex three-dimensional structure with anisotropical properties. Then, to define the electrical activity, two interpenetrating domains are considered, which are the intracellular and extracellular domains, representing respectively the space inside the cells and the region between them.

A depolarizing prepulse (DPP) is an electrical stimulus that causes the potential difference measured across a neuronal membrane to become more positive or less negative, and precedes another electrical stimulus. DPPs may be of either the voltage or current stimulus variety and have been used to inhibit neural activity, selectively excite neurons, and increase the pain threshold associated with electrocutaneous stimulation.

In quantum information theory, the classical capacity of a quantum channel is the maximum rate at which classical data can be sent over it error-free in the limit of many uses of the channel. Holevo, Schumacher, and Westmoreland proved the following least upper bound on the classical capacity of any quantum channel :

Exponential integrate-and-fire models are compact and computationally efficient nonlinear spiking neuron models with one or two variables. The exponential integrate-and-fire model was first proposed as a one-dimensional model. The most prominent two-dimensional examples are the adaptive exponential integrate-and-fire model and the generalized exponential integrate-and-fire model. Exponential integrate-and-fire models are widely used in the field of computational neuroscience and spiking neural networks because of (i) a solid grounding of the neuron model in the field of experimental neuroscience, (ii) computational efficiency in simulations and hardware implementations, and (iii) mathematical transparency.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

References

  1. Rattay, F. (1986). "Analysis of Models for External Stimulation of Axons". IEEE Transactions on Biomedical Engineering (10): 974–977. doi:10.1109/TBME.1986.325670. S2CID   33053720.
  2. Rattay, F. (1988). "Modeling the excitation of fibers under surface electrodes". IEEE Transactions on Biomedical Engineering. 35 (3): 199–202. doi:10.1109/10.1362. PMID   3350548. S2CID   27312507.
  3. Rattay, F. (1989). "Analysis of models for extracellular fiber stimulation". IEEE Transactions on Biomedical Engineering. 36 (7): 676–682. doi:10.1109/10.32099. PMID   2744791. S2CID   42935757.
  4. Rattay, F. (1990). Electrical Nerve Stimulation: Theory, Experiments and Applications . Wien, New York: Springer. pp.  264. ISBN   3-211-82247-X.
  5. Rattay, F. (1998). "Analysis of the electrical excitation of CNS neurons". IEEE Transactions on Biomedical Engineering. 45 (6): 766–772. doi:10.1109/10.678611. PMID   9609941. S2CID   789370.
  6. Rattay, F. (1999). "The basic mechanism for the electrical stimulation of the nervous system". Neuroscience. 89 (2): 335–346. doi:10.1016/S0306-4522(98)00330-3. PMID   10077317. S2CID   41408689.
  7. Danner, S.M.; Wenger, C.; Rattay, F. (2011). Electrical stimulation of myelinated axons. Saarbrücken: VDM. p. 92. ISBN   978-3-639-37082-9.
  8. Rattay, F.; Greenberg, R.J.; Resatz, S. (2003). "Neuron modeling". Handbook of Neuroprosthetic Methods. CRC Press. ISBN   978-0-8493-1100-0.
  9. McNeal, D. R. (1976). "Analysis of a Model for Excitation of Myelinated Nerve". IEEE Transactions on Biomedical Engineering. BME-23 (4): 329–337. doi:10.1109/TBME.1976.324593. PMID   1278925. S2CID   22334434.