Quasicycles in the stochastic hybrid Morris-Lecar neural model

Conceptual foundations of physiological regulation incorporating the free energy principle and self-organized criticality

Bettinger, J. S., K. J. Friston. Conceptual Foundations of Physiological Regulation incorporating the Free Energy Principle & Self-Organized Criticality. NEUROSCI BIOBEHAV REV 23(x) 144-XXX, 2022. Since the late nineteen-nineties, the concept of homeostasis has been contextualized within a broader class of "allostatic" dynamics characterized by a wider-berth of causal factors including social, psychological and environmental entailments; the fundamental nature of integrated brain-body dynamics; plus the role of anticipatory, top-down constraints supplied by intrinsic regulatory models. Many of these evidentiary factors are integral in original descriptions of homeostasis; subsequently integrated; and/or cite more-general operating principles of self-organization. As a result, the concept of allostasis may be generalized to a larger category of variational systems in biology, engineering and physics in terms of advances in complex systems, statistical mechanics and dynamics involving heterogenous (hierarchical/heterarchical, modular) systems like brain-networks and the internal milieu. This paper offers a three-part treatment. 1) interpret "allostasis" to emphasize a variational and relational foundation of physiological stability; 2) adapt the role of allostasis as "stability through change" to include a "return to stability" and 3) reframe the model of homeostasis with a conceptual model of criticality that licenses the upgrade to variational dynamics.

A universal description of stochastic oscillators

Many systems in physics, chemistry, and biology exhibit oscillations with a pronounced random component. Such stochastic oscillations can emerge via different mechanisms, for example, linear dynamics of a stable focus with fluctuations, limit-cycle systems perturbed by noise, or excitable systems in which random inputs lead to a train of pulses. Despite their diverse origins, the phenomenology of random oscillations can be strikingly similar. Here, we introduce a nonlinear transformation of stochastic oscillators to a complex-valued function [Formula: see text](x) that greatly simplifies and unifies the mathematical description of the oscillator's spontaneous activity, its response to an external time-dependent perturbation, and the correlation statistics of different oscillators that are weakly coupled. The function [Formula: see text] (x) is the eigenfunction of the Kolmogorov backward operator with the least negative (but nonvanishing) eigenvalue λ1 = μ1 + iω1. The resulting power spectrum of the complex-valued function is exactly given by a Lorentz spectrum with peak frequency ω1 and half-width μ1; its susceptibility with respect to a weak external forcing is given by a simple one-pole filter, centered around ω1; and the cross-spectrum between two coupled oscillators can be easily expressed by a combination of the spontaneous power spectra of the uncoupled systems and their susceptibilities. Our approach makes qualitatively different stochastic oscillators comparable, provides simple characteristics for the coherence of the random oscillation, and gives a framework for the description of weakly coupled oscillators.

Low Cost Digital Implementation of Hybrid FitzHugh Nagumo-Morris Lecar Neuron Model Considering Electromagnetic Flux Coupling

Digital realization of neuron models, especially implementation on a field programmable gate array (FPGA), is one of the key objectives of neuromorphic research, because the effective hardware realization of the biological neural networks plays a crucial role in implementing the behaviors of the brain for future applications. In this paper, a hybrid FitzHugh Nagumo-Morris Lecar (FNML) neuron model with electromagnetic flux coupling is considered, and two multiplierless piecewise linear (PWL) models, which have similar behaviors to the biological neuron, are presented. A comparison between digital implementation results of the original FNML and PWL models illustrates that, the PWL1 model provides a 65% speed-up with an overall saving (in FPGA resources) of 66.2%, and the PWL2 model yields a 71% speed-up with an overall saving of 78.2%.

Isostables for Stochastic Oscillators

Thomas and Lindner [P. J. Thomas and B. Lindner, Phys. Rev. Lett. 113, 254101 (2014).PRLTAO0031-900710.1103/PhysRevLett.113.254101], defined an asymptotic phase for stochastic oscillators as the angle in the complex plane made by the eigenfunction, having a complex eigenvalue with a least negative real part, of the backward Kolmogorov (or stochastic Koopman) operator. We complete the phase-amplitude description of noisy oscillators by defining the stochastic isostable coordinate as the eigenfunction with the least negative nontrivial real eigenvalue. Our results suggest a framework for stochastic limit cycle dynamics that encompasses noise-induced oscillations.

Noise induced suppression of spiral waves in a hybrid FitzHugh-Nagumo neuron with discontinuous resetting

A modified FitzHugh-Nagumo neuron model with sigmoid function-based recovery variable is considered with electromagnetic flux coupling. The dynamical properties of the proposed neuron model are investigated, and as the excitation current becomes larger, the number of fixed points decreases to one. The bifurcation plots are investigated to show the chaotic and periodic regimes for various values of excitation current and parameters. A N×N network of the neuron model is constructed to study the wave propagation and wave re-entry phenomena. Investigations are conducted to show that for larger flux coupling values, the spiral waves are suppressed, but for such values of the flux coupling, the individual nodes are driven into periodic regimes. By introducing Gaussian noise as an additional current term, we showed that when noise is introduced for the entire simulation time, the dynamics of the nodes are largely altered while the noise exposure for 200-time units will not alter the dynamics of the nodes completely.

Firing statistics in the bistable regime of neurons with homoclinic spike generation

Neuronal voltage dynamics of regularly firing neurons typically has one stable attractor: either a fixed point (like in the subthreshold regime) or a limit cycle that defines the tonic firing of action potentials (in the suprathreshold regime). In two of the three spike onset bifurcation sequences that are known to give rise to all-or-none type action potentials, however, the resting-state fixed point and limit cycle spiking can coexist in an intermediate regime, resulting in bistable dynamics. Here, noise can induce switches between the attractors, i.e., between rest and spiking, and thus increase the variability of the spike train compared to neurons with only one stable attractor. Qualitative features of the resulting spike statistics depend on the spike onset bifurcations. This paper focuses on the creation of the spiking limit cycle via the saddle-homoclinic orbit (HOM) bifurcation and derives interspike interval (ISI) densities for a conductance-based neuron model in the bistable regime. The ISI densities of bistable homoclinic neurons are found to be unimodal yet distinct from the inverse Gaussian distribution associated with the saddle-node-on-invariant-cycle bifurcation. It is demonstrated that for the HOM bifurcation the transition between rest and spiking is mainly determined along the downstroke of the action potential-a dynamical feature that is not captured by the commonly used reset neuron models. The deduced spike statistics can help to identify HOM dynamics in experimental data.

Multiple Firing Patterns in Coupled Hindmarsh-Rose Neurons with a Nonsmooth Memristor

A model is introduced by coupling two three-dimensional Hindmarsh-Rose models with the help of a nonsmooth memristor. The firing patterns dependent on the external forcing current are explored, which undergo a process from adding-period to chaos. The stability of equilibrium points of the considered model is investigated via qualitative analysis, from which it can be gained that the model has diversity in the number and stability of equilibrium points for different coupling coefficients. The coexistence of multiple firing patterns relative to initial values is revealed, which means that the referred model can appear various firing patterns with the change of the initial value. Multiple firing patterns of the addressed neuron model induced by different scales are uncovered, which suggests that the discussed model has a multiscale effect for the nonzero initial value.

Phase descriptions of a multidimensional Ornstein-Uhlenbeck process

Stochastic oscillators play a prominent role in different fields of science. Their simplified description in terms of a phase has been advocated by different authors using distinct phase definitions in the stochastic case. One notion of phase that we put forward previously, the asymptotic phase of a stochastic oscillator, is based on the eigenfunction expansion of its probability density. More specifically, it is given by the complex argument of the eigenfunction of the backward operator corresponding to the least-negative eigenvalue. Formally, besides the "backward" phase, one can also define the "forward" phase as the complex argument of the eigenfunction of the forward Kolomogorov operator corresponding to the least-negative eigenvalue. Until now, the intuition about these phase descriptions has been limited. Here we study these definitions for a process that is analytically tractable, the two-dimensional Ornstein-Uhlenbeck process with complex eigenvalues. For this process, (i) we give explicit expressions for the two phases; (ii) we demonstrate that the isochrons are always the spokes of a wheel but that (iii) the spacing of these isochrons (their angular density) is different for backward and forward phases; (iv) we show that the isochrons of the backward phase are completely determined by the deterministic part of the vector field, whereas the forward phase also depends on the noise matrix; and (v) we demonstrate that the mean progression of the backward phase in time is always uniform, whereas this is not true for the forward phase except in the rotationally symmetric case. We illustrate our analytical results for a number of qualitatively different cases.

Noise induced escape in one-population and two-population stochastic neural networks with internal states

In the present paper, the escapes from the basins of fixed points induced by intrinsic noise are investigated in both one- and two-population stochastic hybrid neural networks. In the weak noise limit, the quasipotentials are computed through the application of WKB approximation to the original hybrid system and the results of quasi-steady-state (QSS) diffusion approximation. It is seen that the two results are consistent with each other within the neighborhood of a fixed point and an obvious discrepancy arises in the other area, of which the reason is then explored and revealed. Furthermore, the relationship between the fluctuational paths and the relaxational ones is analyzed, based on which some specific results for the hybrid system is obtained. Besides, for the two-population model, the phenomenon of saddle point avoidance is investigated by using both theoretical and numerical methods. Finally, the topological structure of Lagrangian manifold is analyzed, and its particular features and something analogous to the stochastic differential equation are found according to the accuracy of QSS within the vicinity of the saddle point.

Stochastic Hybrid Systems in Cellular Neuroscience

We review recent work on the theory and applications of stochastic hybrid systems in cellular neuroscience. A stochastic hybrid system or piecewise deterministic Markov process involves the coupling between a piecewise deterministic differential equation and a time-homogeneous Markov chain on some discrete space. The latter typically represents some random switching process. We begin by summarizing the basic theory of stochastic hybrid systems, including various approximation schemes in the fast switching (weak noise) limit. In subsequent sections, we consider various applications of stochastic hybrid systems, including stochastic ion channels and membrane voltage fluctuations, stochastic gap junctions and diffusion in randomly switching environments, and intracellular transport in axons and dendrites. Finally, we describe recent work on phase reduction methods for stochastic hybrid limit cycle oscillators.

A variational method for analyzing limit cycle oscillations in stochastic hybrid systems

Many systems in biology can be modeled through ordinary differential equations, which are piece-wise continuous, and switch between different states according to a Markov jump process known as a stochastic hybrid system or piecewise deterministic Markov process (PDMP). In the fast switching limit, the dynamics converges to a deterministic ODE. In this paper, we develop a phase reduction method for stochastic hybrid systems that support a stable limit cycle in the deterministic limit. A classic example is the Morris-Lecar model of a neuron, where the switching Markov process is the number of open ion channels and the continuous process is the membrane voltage. We outline a variational principle for the phase reduction, yielding an exact analytic expression for the resulting phase dynamics. We demonstrate that this decomposition is accurate over timescales that are exponential in the switching rate ϵ^-1. That is, we show that for a constant C, the probability that the expected time to leave an O(a) neighborhood of the limit cycle is less than T scales as T exp (-Ca/ϵ).

Minimum Action Path Theory Reveals the Details of Stochastic Transitions Out of Oscillatory States

Cell state determination is the outcome of intrinsically stochastic biochemical reactions. Transitions between such states are studied as noise-driven escape problems in the chemical species space. Escape can occur via multiple possible multidimensional paths, with probabilities depending nonlocally on the noise. Here we characterize the escape from an oscillatory biochemical state by minimizing the Freidlin-Wentzell action, deriving from it the stochastic spiral exit path from the limit cycle. We also use the minimized action to infer the escape time probability density function.