313 research outputs found
Mean field approximation of two coupled populations of excitable units
The analysis on stability and bifurcations in the macroscopic dynamics
exhibited by the system of two coupled large populations comprised of
stochastic excitable units each is performed by studying an approximate system,
obtained by replacing each population with the corresponding mean-field model.
In the exact system, one has the units within an ensemble communicating via the
time-delayed linear couplings, whereas the inter-ensemble terms involve the
nonlinear time-delayed interaction mediated by the appropriate global
variables. The aim is to demonstrate that the bifurcations affecting the
stability of the stationary state of the original system, governed by a set of
4N stochastic delay-differential equations for the microscopic dynamics, can
accurately be reproduced by a flow containing just four deterministic
delay-differential equations which describe the evolution of the mean-field
based variables. In particular, the considered issues include determining the
parameter domains where the stationary state is stable, the scenarios for the
onset and the time-delay induced suppression of the collective mode, as well as
the parameter domains admitting bistability between the equilibrium and the
oscillatory state. We show how analytically tractable bifurcations occurring in
the approximate model can be used to identify the characteristic mechanisms by
which the stationary state is destabilized under different system
configurations, like those with symmetrical or asymmetrical inter-population
couplings.Comment: 5 figure
Stability and Hopf bifurcation of controlled complex networks model with two delays
none3siopenJinde Cao, Luca Guerrini, Zunshui ChengCao, Jinde; Guerrini, Luca; Cheng, Zunshu
The role of fixed delays in neuronal rate models
Fixed delays in neuronal interactions arise through synaptic and dendritic processing. Previous work has shown that such delays, which play an important role in shaping the dynamics of networks of large numbers of spiking neurons with continuous synaptic kinetics, can be taken into account with a rate model through the addition of an explicit, fixed delay. Here we extend this work to account for arbitrary symmetric patterns of synaptic connectivity and generic nonlinear transfer functions. Specifically, we conduct a weakly nonlinear analysis of the dynamical states arising via primary instabilities of the stationary uniform state. In this way we determine analytically how the nature and stability of these states depend on the choice of transfer function and connectivity. While this dependence is, in general, nontrivial, we make use of the smallness of the ratio in the delay in neuronal interactions to the effective time constant of integration to arrive at two general observations of physiological relevance. These are: 1 - fast oscillations are always supercritical for realistic transfer functions. 2 - Traveling waves are preferred over standing waves given plausible patterns of local connectivity
Chimera states: Coexistence of coherence and incoherence in networks of coupled oscillators
A chimera state is a spatio-temporal pattern in a network of identical
coupled oscillators in which synchronous and asynchronous oscillation coexist.
This state of broken symmetry, which usually coexists with a stable spatially
symmetric state, has intrigued the nonlinear dynamics community since its
discovery in the early 2000s. Recent experiments have led to increasing
interest in the origin and dynamics of these states. Here we review the history
of research on chimera states and highlight major advances in understanding
their behaviour.Comment: 26 pages, 3 figure
Metabifurcation analysis of a mean field model of the cortex
Mean field models (MFMs) of cortical tissue incorporate salient features of
neural masses to model activity at the population level. One of the common
aspects of MFM descriptions is the presence of a high dimensional parameter
space capturing neurobiological attributes relevant to brain dynamics. We study
the physiological parameter space of a MFM of electrocortical activity and
discover robust correlations between physiological attributes of the model
cortex and its dynamical features. These correlations are revealed by the study
of bifurcation plots, which show that the model responses to changes in
inhibition belong to two families. After investigating and characterizing
these, we discuss their essential differences in terms of four important
aspects: power responses with respect to the modeled action of anesthetics,
reaction to exogenous stimuli, distribution of model parameters and oscillatory
repertoires when inhibition is enhanced. Furthermore, while the complexity of
sustained periodic orbits differs significantly between families, we are able
to show how metamorphoses between the families can be brought about by
exogenous stimuli. We unveil links between measurable physiological attributes
of the brain and dynamical patterns that are not accessible by linear methods.
They emerge when the parameter space is partitioned according to bifurcation
responses. This partitioning cannot be achieved by the investigation of only a
small number of parameter sets, but is the result of an automated bifurcation
analysis of a representative sample of 73,454 physiologically admissible sets.
Our approach generalizes straightforwardly and is well suited to probing the
dynamics of other models with large and complex parameter spaces
Dynamics of a Fractional-Order Delayed Model of COVID-19 with Vaccination Efficacy
In this study, we provide a fractional-order mathematical model that considers the effect of vaccination on COVID-19 spread dynamics. The model accounts for the latent period of intervention strategies by incorporating a time delay τ. A basic reproduction number, R0, is determined for the model, and prerequisites for endemic equilibrium are discussed. The model’s endemic equilibrium point also exhibits local asymptotic stability (under certain conditions), and a Hopf bifurcation condition is established. Different scenarios of vaccination efficacy are simulated. As a result of the vaccination efforts, the number of deaths and those affected have decreased. COVID-19 may not be effectively controlled by vaccination alone. To control infections, several non-pharmacological interventions are necessary. Based on numerical simulations and fitting to real observations, the theoretical results are proven to be effective
Dynamics and precursor signs for phase transitions in neural systems
This thesis investigates neural state transitions associated with sleep, seizure and anaesthesia. The aim is to address the question: How does a brain traverse the critical threshold between distinct cortical states, both healthy and pathological? Specifically we are interested in sub-threshold neural behaviour immediately prior to state transition. We use theoretical neural modelling (single spiking neurons, a network of these, and a mean-field continuum limit) and in vitro experiments to address this question.
Dynamically realistic equations of motion for thalamic relay neuron, reticular nuclei, cortical pyramidal and cortical interneuron in different vigilance states are developed, based on the Izhikevich spiking neuron model. A network of cortical neurons is assembled to examine the behaviour of the gamma-producing cortical network and its transition to lower frequencies due to effect of anaesthesia. Then a three-neuron model for the thalamocortical loop for sleep spindles is presented. Numerical simulations of these networks confirms spiking consistent with reported in vivo measurement results, and provides supporting evidence for precursor indicators of imminent phase transition due to occurrence of individual spindles.
To complement the spiking neuron networks, we study the Wilson–Cowan neural mass equations describing homogeneous cortical columns and a 1D spatial cluster of such columns. The abstract representation of cortical tissue by a pair of coupled integro-differential equations permits thorough linear stability, phase plane and bifurcation analyses. This model shows a rich set of spatial and temporal bifurcations marking the boundary to state transitions: saddle-node, Hopf, Turing, and mixed Hopf–Turing. Close to state transition, white-noise-induced subthreshold fluctuations show clear signs of critical slowing down with prolongation and strengthening of autocorrelations, both in time and space, irrespective of bifurcation type.
Attempts at in vitro capture of these predicted leading indicators form the last part of the thesis. We recorded local field potentials (LFPs) from cortical and hippocampal slices of mouse brain. State transition is marked by the emergence and cessation of spontaneous seizure-like events (SLEs) induced by bathing the slices in an artificial cerebral spinal fluid containing no magnesium ions. Phase-plane analysis of the LFP time-series suggests that distinct bifurcation classes can be responsible for state change to seizure. Increased variance and growth of spectral power at low frequencies (f < 15 Hz) was observed in LFP recordings prior to initiation of some SLEs. In addition we demonstrated prolongation of electrically evoked potentials in cortical tissue, while forwarding the slice to a seizing regime. The results offer the possibility of capturing leading temporal indicators prior to seizure generation, with potential consequences for understanding epileptogenesis.
Guided by dynamical systems theory this thesis captures evidence for precursor signs of phase transitions in neural systems using mathematical and computer-based modelling as well as in vitro experiments
Recommended from our members
Dynamics of neural systems with time delays
Complex networks are ubiquitous in nature. Numerous neurological diseases, such as
Alzheimer's, Parkinson's, epilepsy are caused by the abnormal collective behaviour of
neurons in the brain. In particular, there is a strong evidence that Parkinson's disease is
caused by the synchronisation of neurons, and understanding how and why such synchronisation
occurs will bring scientists closer to the design and implementation of appropriate
control to support desynchronisation required for the normal functioning of the brain. In
order to study the emergence of (de)synchronisation, it is necessary first to understand
how the dynamical behaviour of the system under consideration depends on the changes
in systems parameters. This can be done using a powerful mathematical method, called
bifurcation analysis, which allows one to identify and classify different dynamical regimes,
such as, for example, stable/unstable steady states, Hopf and fold bifurcations, and find
periodic solutions by varying parameters of the nonlinear system.
In real-world systems, interactions between elements do not happen instantaneously
due to a finite time of signal propagation, reaction times of individual elements, etc.
Moreover, time delays are normally non-constant and may vary with time. This means
that it is vital to introduce time delays in any realistic model of neural networks. In
this thesis, I consider four different models. First, in order to analyse the fundamental
properties of neural networks with time-delayed connections, I consider a system of four
coupled nonlinear delay differential equations. This model represents a neural network,
where one subsystem receives a delayed input from another subsystem. The exciting
feature of this model is the combination of both discrete and distributed time delays, where
distributed time delays represent the neural feedback between the two sub-systems, and the
discrete delays describe neural interactions within each of the two subsystems. Stability
properties are investigated for different commonly used distribution kernels, and the results
are compared to the corresponding stability results for networks with no distributed delays.
It is shown how approximations to the boundary of stability region of an equilibrium point
can be obtained analytically for the cases of delta, uniform, and gamma delay distributions.
Numerical techniques are used to investigate stability properties of the fully nonlinear
system and confirm our analytical findings.
In the second part of this thesis, I consider a globally coupled network composed of
active (oscillatory) and inactive (non-oscillatory) oscillators with distributed time delayed
coupling. Analytical conditions for the amplitude death, where the oscillations are quenched,
are obtained in terms of the coupling strength, the ratio of inactive oscillators, the width
of the uniformly distributed delay and the mean time delay for gamma distribution. The
results show that for uniform distribution, by increasing both the width of the delay distribution
and the ratio of inactive oscillators, the amplitude death region increases in the
mean time delay and the coupling strength parameter space. In the case of the gamma
distribution kernel, we find the amplitude death region in the space of the ratio of inactive
oscillators, the mean time delay for gamma distribution, and the coupling strength for
both weak and strong gamma distribution kernels.
Furthermore, I analyse a model of the subthalamic nucleus (STN)-globus palidus (GP)
network with three different transmission delays. A time-shift transformation reduces the
model to a system with two time delays, for which the existence of a unique steady
state is established. Conditions for stability of the steady state are derived in terms of
system parameters and the time delays. Numerical stability analysis is performed using
traceDDE and DDE-BIFTOOL in Matlab to investigate different dynamical regimes in
the STN-GP model, and to obtain critical stability boundaries separating stable (healthy)
and oscillatory (Parkinsonian-like) neural ring. Direct numerical simulations of the fully
nonlinear system are performed to confirm analytical findings, and to illustrate different
dynamical behaviours of the system.
Finally, I consider a ring of n neurons coupled through the discrete and distributed
time delays. I show that the amplitude death occurs in the symmetric (asymmetric) region
depending on the even (odd) number of neurons in the ring neural system. Analytical
conditions for linear stability of the trivial steady state are represented in a parameter space
of the synaptic weight of the self-feedback and the coupling strength between the connected
neurons, as well as in the space of the delayed self-feedback and the coupling strength
between the neurons. It is shown that both Hopf and steady-state bifurcations may occur
when the steady state loses its stability. Stability properties are also investigated for
different commonly used distribution kernels, such as delta function and weak gamma
distributions. Moreover, the obtained analytical results are confirmed by the numerical
simulations of the fully nonlinear system
- …