Diverse electrical responses in a network of fractional-order conductance-based excitable Morris-Lecar systems.

The diverse excitabilities of cells often produce various spiking-bursting oscillations that are found in the neural system. We establish the ability of a fractional-order excitable neuron model with Caputo's fractional derivative to analyze the effects of its dynamics on the spike train features observed in our results. The significance of this generalization relies on a theoretical framework of the model in which memory and hereditary properties are considered. Employing the fractional exponent, we first provide information about the variations in electrical activities. We deal with the 2D class I and class II excitable Morris-Lecar (M-L) neuron models that show the alternation of spiking and bursting features including MMOs & MMBOs of an uncoupled fractional-order neuron. We then extend the study with the 3D slow-fast M-L model in the fractional domain. The considered approach establishes a way to describe various characteristics similarities between fractional-order and classical integer-order dynamics. Using the stability and bifurcation analysis, we discuss different parameter spaces where the quiescent state emerges in uncoupled neurons. We show the characteristics consistent with the analytical results. Next, the Erdös-Rényi network of desynchronized mixed neurons (oscillatory and excitable) is constructed that is coupled through membrane voltage. It can generate complex firing activities where quiescent neurons start to fire. Furthermore, we have shown that increasing coupling can create cluster synchronization, and eventually it can enable the network to fire in unison. Based on cluster synchronization, we develop a reduced-order model which can capture the activities of the entire network. Our results reveal that the effect of fractional-order depends on the synaptic connectivity and the memory trace of the system. Additionally, the dynamics captures spike frequency adaptation and spike latency that occur over multiple timescales as the effects of fractional derivative, which has been observed in neural computation.

The generation of diverse traveling pulses and its solution scheme in an excitable slow-fast dynamics.

In this article, we report on the generation and propagation of traveling pulses in a homogeneous network of diffusively coupled, excitable, slow-fast dynamical neurons. The spatially extended system is modeled using the nearest neighbor coupling theory, in which the diffusion part measures the spatial distribution of coupling topology. We derive analytically the conditions for traveling wave profiles that allow the construction of the shape of traveling nerve impulses. The analytical and numerical results are used to explore the nature of propagating pulses. The symmetric or asymmetric nature of traveling pulses is characterized, and the wave velocity is derived as a function of system parameters. Moreover, we present our results for an extended excitable medium by considering a slow-fast biophysical model with a homogeneous, diffusive coupling that can exhibit various traveling pulses. The appearance of series of pulses is an interesting phenomenon from biophysical and dynamical perspective. Varying the perturbation and coupling parameters, we observe the propagation of activities with various amplitude modulations and transition phases of different wave profiles that affect the speed of pulses in certain parameter regimes. We observe different types of traveling pulses, such as envelope solitons and multi-bump solutions, and show how system parameters and coupling play a major role in the formation of different traveling pulses. Finally, we obtain the conditions for stable and unstable plane waves.

A new model for freedom of movement using connectomic analysis.

The problem of whether we can execute free acts or not is central in philosophical thought, and it has been studied by numerous scholars throughout the centuries. Recently, neurosciences have entered this topic contributing new data and insights into the neuroanatomical basis of cognitive processes. With the advent of connectomics, a more refined landscape of brain connectivity can be analysed at an unprecedented level of detail. Here, we identify the connectivity network involved in the movement process from a connectomics point of view, from its motivation through its execution until the sense of agency develops. We constructed a "volitional network" using data derived from the Brainnetome Atlas database considering areas involved in volitional processes as known in the literature. We divided this process into eight processes and used Graph Theory to measure several structural properties of the network. Our results show that the volitional network is small-world and that it contains four communities. Nodes of the right hemisphere are contained in three of these communities whereas nodes of the left hemisphere only in two. Centrality measures indicate the nucleus accumbens is one of the most connected nodes in the network. Extensive connectivity is observed in all processes except in Decision (to move) and modulation of Agency, which might correlate with a mismatch mechanism for perception of Agency.

A generic model for pandemics in networks of communities and the role of vaccination.

The slogan "nobody is safe until everybody is safe" is a dictum to raise awareness that in an interconnected world, pandemics, such as COVID-19, require a global approach. Motivated by the ongoing COVID-19 pandemic, we model here the spread of a virus in interconnected communities and explore different vaccination scenarios, assuming that the efficacy of the vaccination wanes over time. We start with susceptible populations and consider a susceptible-vaccinated-infected-recovered model with unvaccinated ("Bronze"), moderately vaccinated ("Silver"), and very-well-vaccinated ("Gold") communities, connected through different types of networks via a diffusive linear coupling for local spreading. We show that when considering interactions in "Bronze"-"Gold" and "Bronze"-"Silver" communities, the "Bronze" community is driving an increase in infections in the "Silver" and "Gold" communities. This shows a detrimental, unidirectional effect of non-vaccinated to vaccinated communities. Regarding the interactions between "Gold," "Silver," and "Bronze" communities in a network, we find that two factors play a central role: the coupling strength in the dynamics and network density. When considering the spread of a virus in Barabási-Albert networks, infections in "Silver" and "Gold" communities are lower than in "Bronze" communities. We find that the "Gold" communities are the best in keeping their infection levels low. However, a small number of "Bronze" communities are enough to give rise to an increase in infections in moderately and well-vaccinated communities. When studying the spread of a virus in dense Erdos-Rényi and sparse Watts-Strogatz and Barabási-Albert networks, the communities reach the disease-free state in the dense Erdos-Rényi networks, but not in the sparse Watts-Strogatz and Barabási-Albert networks. However, we also find that if all these networks are dense enough, all types of communities reach the disease-free state. We conclude that the presence of a few unvaccinated or partially vaccinated communities in a network can increase significantly the rate of infected population in other communities. This reveals the necessity of a global effort to facilitate access to vaccines for all communities.

Dynamical analysis of the infection status in diverse communities due to COVID-19 using a modified SIR model.

In this article, we model and study the spread of COVID-19 in Germany, Japan, India and highly impacted states in India, i.e., in Delhi, Maharashtra, West Bengal, Kerala and Karnataka. We consider recorded data published in Worldometers and COVID-19 India websites from April 2020 to July 2021, including periods of interest where these countries and states were hit severely by the pandemic. Our methodology is based on the classic susceptible-infected-removed (SIR) model and can track the evolution of infections in communities, i.e., in countries, states or groups of individuals, where we (a) allow for the susceptible and infected populations to be reset at times where surges, outbreaks or secondary waves appear in the recorded data sets, (b) consider the parameters in the SIR model that represent the effective transmission and recovery rates to be functions of time and (c) estimate the number of deaths by combining the model solutions with the recorded data sets to approximate them between consecutive surges, outbreaks or secondary waves, providing a more accurate estimate. We report on the status of the current infections in these countries and states, and the infections and deaths in India and Japan. Our model can adapt to the recorded data and can be used to explain them and importantly, to forecast the number of infected, recovered, removed and dead individuals, as well as it can estimate the effective infection and recovery rates as functions of time, assuming an outbreak occurs at a given time. The latter information can be used to forecast the future basic reproduction number and together with the forecast on the number of infected and dead individuals, our approach can further be used to suggest the implementation of intervention strategies and mitigation policies to keep at bay the number of infected and dead individuals. This, in conjunction with the implementation of vaccination programs worldwide, can help reduce significantly the impact of the spread around the world and improve the wellbeing of people.

Spatiotemporal characteristics in systems of diffusively coupled excitable slow-fast FitzHugh-Rinzel dynamical neurons.

In this paper, we study an excitable, biophysical system that supports wave propagation of nerve impulses. We consider a slow-fast, FitzHugh-Rinzel neuron model where only the membrane voltage interacts diffusively, giving rise to the formation of spatiotemporal patterns. We focus on local, nonlinear excitations and diverse neural responses in an excitable one- and two-dimensional configuration of diffusively coupled FitzHugh-Rinzel neurons. The study of the emerging spatiotemporal patterns is essential in understanding the working mechanism in different brain areas. We derive analytically the coefficients of the amplitude equations in the vicinity of Hopf bifurcations and characterize various patterns, including spirals exhibiting complex geometric substructures. Furthermore, we derive analytically the condition for the development of antispirals in the neighborhood of the bifurcation point. The emergence of broken target waves can be observed to form spiral-like profiles. The spatial dynamics of the excitable system exhibits two- and multi-arm spirals for small diffusive couplings. Our results reveal a multitude of neural excitabilities and possible conditions for the emergence of spiral-wave formation. Finally, we show that the coupled excitable systems with different firing characteristics participate in a collective behavior that may contribute significantly to irregular neural dynamics.

Influence of Autapses on Synchronization in Neural Networks With Chemical Synapses.

A great deal of research has been devoted on the investigation of neural dynamics in various network topologies. However, only a few studies have focused on the influence of autapses, synapses from a neuron onto itself via closed loops, on neural synchronization. Here, we build a random network with adaptive exponential integrate-and-fire neurons coupled with chemical synapses, equipped with autapses, to study the effect of the latter on synchronous behavior. We consider time delay in the conductance of the pre-synaptic neuron for excitatory and inhibitory connections. Interestingly, in neural networks consisting of both excitatory and inhibitory neurons, we uncover that synchronous behavior depends on their synapse type. Our results provide evidence on the synchronous and desynchronous activities that emerge in random neural networks with chemical, inhibitory and excitatory synapses where neurons are equipped with autapses.

Labyrinth chaos: Revisiting the elegant, chaotic, and hyperchaotic walks.

Labyrinth chaos was discovered by Otto Rössler and René Thomas in their endeavor to identify the necessary mathematical conditions for the appearance of chaotic and hyperchaotic motion in continuous flows. Here, we celebrate their discovery by considering a single labyrinth walk system and an array of coupled labyrinth chaos systems that exhibit complex, chaotic behavior, reminiscent of chimera-like states, a peculiar synchronization phenomenon. We discuss the properties of the single labyrinth walk system and review the ability of coupled labyrinth chaos systems to exhibit chimera-like states due to the unique properties of their space-filling, chaotic trajectories, which amounts to elegant, hyperchaotic walks. Finally, we discuss further implications in relation to the labyrinth walk system by showing that even though it is volume-preserving, it is not force-conservative.

Dynamic tracking with model-based forecasting for the spread of the COVID-19 pandemic.

In this paper, a susceptible-infected-removed (SIR) model has been used to track the evolution of the spread of COVID-19 in four countries of interest. In particular, the epidemic model, that depends on some basic characteristics, has been applied to model the evolution of the disease in Italy, India, South Korea and Iran. The economic, social and health consequences of the spread of the virus have been cataclysmic. Hence, it is imperative that mathematical models can be developed and used to compare published datasets with model predictions. The predictions estimated from the presented methodology can be used in both the qualitative and quantitative analysis of the spread. They give an insight into the spread of the virus that the published data alone cannot, by updating them and the model on a daily basis. We show that by doing so, it is possible to detect the early onset of secondary spikes in infections or the development of secondary waves. We considered data from March to August, 2020, when different communities were affected severely and demonstrate predictions depending on the model's parameters related to the spread of COVID-19 until the end of December, 2020. By comparing the published data with model results, we conclude that in this way, it may be possible to reflect better the success or failure of the adequate measures implemented by governments and authorities to mitigate and control the current pandemic.

A SIR model assumption for the spread of COVID-19 in different communities.

In this paper, we study the effectiveness of the modelling approach on the pandemic due to the spreading of the novel COVID-19 disease and develop a susceptible-infected-removed (SIR) model that provides a theoretical framework to investigate its spread within a community. Here, the model is based upon the well-known susceptible-infected-removed (SIR) model with the difference that a total population is not defined or kept constant per se and the number of susceptible individuals does not decline monotonically. To the contrary, as we show herein, it can be increased in surge periods! In particular, we investigate the time evolution of different populations and monitor diverse significant parameters for the spread of the disease in various communities, represented by China, South Korea, India, Australia, USA, Italy and the state of Texas in the USA. The SIR model can provide us with insights and predictions of the spread of the virus in communities that the recorded data alone cannot. Our work shows the importance of modelling the spread of COVID-19 by the SIR model that we propose here, as it can help to assess the impact of the disease by offering valuable predictions. Our analysis takes into account data from January to June, 2020, the period that contains the data before and during the implementation of strict and control measures. We propose predictions on various parameters related to the spread of COVID-19 and on the number of susceptible, infected and removed populations until September 2020. By comparing the recorded data with the data from our modelling approaches, we deduce that the spread of COVID-19 can be under control in all communities considered, if proper restrictions and strong policies are implemented to control the infection rates early from the spread of the disease.

Emergence of Mixed Mode Oscillations in Random Networks of Diverse Excitable Neurons: The Role of Neighbors and Electrical Coupling.

In this paper, we focus on the emergence of diverse neuronal oscillations arising in a mixed population of neurons with different excitability properties. These properties produce mixed mode oscillations (MMOs) characterized by the combination of large amplitudes and alternate subthreshold or small amplitude oscillations. Considering the biophysically plausible, Izhikevich neuron model, we demonstrate that various MMOs, including MMBOs (mixed mode bursting oscillations) and synchronized tonic spiking appear in a randomly connected network of neurons, where a fraction of them is in a quiescent (silent) state and the rest in self-oscillatory (firing) states. We show that MMOs and other patterns of neural activity depend on the number of oscillatory neighbors of quiescent nodes and on electrical coupling strengths. Our results are verified by constructing a reduced-order network model and supported by systematic bifurcation diagrams as well as for a small-world network. Our results suggest that, for weak couplings, MMOs appear due to the de-synchronization of a large number of quiescent neurons in the networks. The quiescent neurons together with the firing neurons produce high frequency oscillations and bursting activity. The overarching goal is to uncover a favorable network architecture and suitable parameter spaces where Izhikevich model neurons generate diverse responses ranging from MMOs to tonic spiking.

Bistable Firing Pattern in a Neural Network Model.

Excessively high, neural synchronization has been associated with epileptic seizures, one of the most common brain diseases worldwide. A better understanding of neural synchronization mechanisms can thus help control or even treat epilepsy. In this paper, we study neural synchronization in a random network where nodes are neurons with excitatory and inhibitory synapses, and neural activity for each node is provided by the adaptive exponential integrate-and-fire model. In this framework, we verify that the decrease in the influence of inhibition can generate synchronization originating from a pattern of desynchronized spikes. The transition from desynchronous spikes to synchronous bursts of activity, induced by varying the synaptic coupling, emerges in a hysteresis loop due to bistability where abnormal (excessively high synchronous) regimes exist. We verify that, for parameters in the bistability regime, a square current pulse can trigger excessively high (abnormal) synchronization, a process that can reproduce features of epileptic seizures. Then, we show that it is possible to suppress such abnormal synchronization by applying a small-amplitude external current on > 10% of the neurons in the network. Our results demonstrate that external electrical stimulation not only can trigger synchronous behavior, but more importantly, it can be used as a means to reduce abnormal synchronization and thus, control or treat effectively epileptic seizures.

Hyperchaos & labyrinth chaos: Revisiting Thomas-Rössler systems.

We consider a multidimensional extension of Thomas-Rössler systems, that was inspired by Thomas' earlier work on biological feedback circuits, and we report on our first results that shows its ability to sustain spatio-temporal behaviour reminiscent of chimera states. The novelty here is that its underlying mechanism is based on "chaotic walks" discovered by Thomas during the course of his investigations on what he called Labyrinth Chaos. We briefly review the main properties of these systems and their chaotic and hyperchaotic dynamics and discuss the simplest way of coupling, necessary for this spatio-temporal behaviour that allows the emergence of complex dynamical behaviours. We also recall Thomas' memorable influence and interaction with the authors as we dedicate this work to his memory.

Evaluating performance of neural codes in model neural communication networks.

Information needs to be appropriately encoded to be reliably transmitted over physical media. Similarly, neurons have their own codes to convey information in the brain. Even though it is well-known that neurons exchange information using a pool of several protocols of spatio-temporal encodings, the suitability of each code and their performance as a function of network parameters and external stimuli is still one of the great mysteries in neuroscience. This paper sheds light on this by modeling small-size networks of chemically and electrically coupled Hindmarsh-Rose spiking neurons. We focus on a class of temporal and firing-rate codes that result from neurons' membrane-potentials and phases, and quantify numerically their performance estimating the Mutual Information Rate, aka the rate of information exchange. Our results suggest that the firing-rate and interspike-intervals codes are more robust to additive Gaussian white noise. In a network of four interconnected neurons and in the absence of such noise, pairs of neurons that have the largest rate of information exchange using the interspike-intervals and firing-rate codes are not adjacent in the network, whereas spike-timings and phase codes (temporal) promote large rate of information exchange for adjacent neurons. If that result would have been possible to extend to larger neural networks, it would suggest that small microcircuits would preferably exchange information using temporal codes (spike-timings and phase codes), whereas on the macroscopic scale, where there would be typically pairs of neurons not directly connected due to the brain's sparsity, firing-rate and interspike-intervals codes would be the most efficient codes.

Opinion formation in multiplex networks with general initial distributions.

We study opinion dynamics over multiplex networks where agents interact with bounded confidence. Namely, two neighbouring individuals exchange opinions and compromise if their opinions do not differ by more than a given threshold. In literature, agents are generally assumed to have a homogeneous confidence bound. Here, we study analytically and numerically opinion evolution over structured networks characterised by multiple layers with respective confidence thresholds and general initial opinion distributions. Through rigorous probability analysis, we show analytically the critical thresholds at which a phase transition takes place in the long-term consensus behaviour, over multiplex networks with some regularity conditions. Our results reveal the quantitative relation between the critical threshold and initial distribution. Further, our numerical simulations illustrate the consensus behaviour of the agents in network topologies including lattices and, small-world and scale-free networks, as well as for structure-dependent convergence parameters accommodating node heterogeneity. We find that the critical thresholds for consensus tend to agree with the predicted upper bounds in Theorems 4 and 5 in this paper. Finally, our results indicate that multiplexity hinders consensus formation when the initial opinion configuration is within a bounded range and, provide insight into information diffusion and social dynamics in multiplex systems modeled by networks.

Inference of financial networks using the normalised mutual information rate.

In this paper, we study data from financial markets, using the normalised Mutual Information Rate. We show how to use it to infer the underlying network structure of interrelations in the foreign currency exchange rates and stock indices of 15 currency areas. We first present the mathematical method and discuss its computational aspects, and apply it to artificial data from chaotic dynamics and to correlated normal-variates data. We then apply the method to infer the structure of the financial system from the time-series of currency exchange rates and stock indices. In particular, we study and reveal the interrelations among the various foreign currency exchange rates and stock indices in two separate networks, of which we also study their structural properties. Our results show that both inferred networks are small-world networks, sharing similar properties and having differences in terms of assortativity. Importantly, our work shows that global economies tend to connect with other economies world-wide, rather than creating small groups of local economies. Finally, the consistent interrelations depicted among the 15 currency areas are further supported by a discussion from the viewpoint of economics.

Maintaining extensivity in evolutionary multiplex networks.

In this paper, we explore the role of network topology on maintaining the extensive property of entropy. We study analytically and numerically how the topology contributes to maintaining extensivity of entropy in multiplex networks, i.e. networks of subnetworks (layers), by means of the sum of the positive Lyapunov exponents, HKS, a quantity related to entropy. We show that extensivity relies not only on the interplay between the coupling strengths of the dynamics associated to the intra (short-range) and inter (long-range) interactions, but also on the sum of the intra-degrees of the nodes of the layers. For the analytically treated networks of size N, among several other results, we show that if the sum of the intra-degrees (and the sum of inter-degrees) scales as Nθ+1, θ > 0, extensivity can be maintained if the intra-coupling (and the inter-coupling) strength scales as N-θ, when evolution is driven by the maximisation of HKS. We then verify our analytical results by performing numerical simulations in multiplex networks formed by electrically and chemically coupled neurons.