Reducing Neuronal Networks to Discrete Dynamics

We consider a general class of purely inhibitory and excitatory-inhibitory neuronal networks, with a general class of network architectures, and characterize the complex firing patterns that emerge. Our strategy for studying these networks is to first reduce them to a discrete model. In the discrete model, each neuron is represented as a finite number of states and there are rules for how a neuron transitions from one state to another. In this paper, we rigorously demonstrate that the continuous neuronal model can be reduced to the discrete model if the intrinsic and synaptic properties of the cells are chosen appropriately. In a companion paper [W. Just, S. Ahn, D. Terman. Minimal attractors in digraph system models of neuronal networks (preprint)], we analyse the discrete model.     

Plasticity-induced characteristic changes of pattern dynamics and the related phase transitions in small-world neuronal networks

Phase transitions widely exist in nature and occur when some control parameters are changed. In neural systems, their macroscopic states are represented by the activity states of neuron populations, and phase transitions between different activity states are closely related to corresponding functions in the brain. In particular, phase transitions to some rhythmic synchronous firing states play significant roles on diverse brain functions and disfunctions, such as encoding rhythmical external stimuli, epileptic seizure, etc. However, in previous studies, phase transitions in neuronal networks are almost driven by network parameters （e.g., external stimuli）, and there has been no investigation about the transitions between typical activity states of neuronal networks in a self-organized way by applying plastic connection weights. In this paper, we discuss phase transitions in electrically coupled and lattice-based small-world neuronal networks （LBSW networks） under spike-timing-dependent plasticity （STDP）. By applying STDP on all electrical synapses, various known and novel phase transitions could emerge in LBSW networks, particularly, the phenomenon of self-organized phase transitions （SOPTs）： repeated transitions between synchronous and asynchronous firing states. We further explore the mechanics generating SOPTs on the basis of synaptic weight dynamics.     

Voltage-controllable magnetic skyrmion dynamics for spiking neuron device applications

Voltage-controlled magnetic skyrmions have attracted special attention because they satisfy the requirements for well-controlled high-efficiency and energy saving for future skyrmion-based neuron device applications.In this work,we propose a compact leaky-integrate-fire(LIF)spiking neuron device by using the voltage-driven skyrmion dynamics in a multiferroic nanodisk structure.The skyrmion dynamics is controlled by well tailoring voltage-induced piezostrains,where the skyrmion radius can be effectively modulated by applying the piezostrain pulses.Like the biological neuron,the proposed skyrmionic neuron will accumulate a membrane potential as skyrmion radius is varied by inputting the continuous piezostrain spikes,and the skyrmion radius will return to the initial state in the absence of piezostrain.Therefore,this skyrmion radius-based membrane potential will reach a definite threshold value by the strain stimuli and then reset by removing the stimuli.Such the LIF neuronal functionality and the behaviors of the proposed skyrmionic neuron device are elucidated through the micromagnetic simulation studies.Our results may benefit the utilization of skyrmionic neuron for constructing the future energy-efficient and voltage-tunable spiking neural networks.     

Explosive phenomena in complex networks

The emergence of large-scale connectivity and synchronization are crucial to the structure, function and failure of many complex socio-technical networks. Thus, there is great interest in analyzing phase transitions to large-scale connectivity and to global synchronization, including how to enhance or delay the onset. These phenomena are traditionally studied as second-order phase transitions where, at the critical threshold, the order parameter increases rapidly but continuously. In 2009, an extremely abrupt transition was found for a network growth process where links compete for addition in an attempt to delay percolation. This observation of ‘explosive percolation’ was ultimately revealed to be a continuous transition in the thermodynamic limit, yet with very atypical finite-size scaling, and it started a surge of work on explosive phenomena and their consequences. Many related models are now shown to yield discontinuous percolation transitions and even hybrid transitions. Explosive percolation enables many other features such as multiple giant components, modular structures, discrete scale invariance and non-self-averaging, relating to properties found in many real phenomena such as explosive epidemics, electric breakdowns and the emergence of molecular life. Models of explosive synchronization provide an analytic framework for the dynamics of abrupt transitions and reveal the interplay between the distribution in natural frequencies and the network structure, with applications ranging from epileptic seizures to waking from anesthesia. Here we review the vast literature on explosive phenomena in networked systems and synthesize the fundamental connections between models and survey the application areas. We attempt to classify explosive phenomena based on underlying mechanisms and to provide a coherent overview and perspective for future research to address the many vital questions that remained unanswered.     

Hindmarsh-Rose神经元阵列自发产生螺旋波的研究

采用Hindmarsh-Rose(HR)神经元模型,研究了二维神经元阵列系统从一个具有随机相位分布的初态演化最终是否能自发产生螺旋波的问题.数值结果表明:系统是否出现螺旋波与单个HR神经元的状态、系统的初态和耦合强度有关,其中单个HR神经元的振荡状态起主要作用.当单个HR神经元处于一周期振荡态时,在一定的耦合强度范围内系统都会自发出现多个螺旋波和螺旋波对,出现螺旋波与系统初态无关,只要适当选择耦合强度,在系统中可以出现单个螺旋波.当耦合强度超过某一阈值后,继续增加耦合强度,系统会呈现三种不同的动力学行为,分别与三类初态有关.系统从第一类初态演化将偶尔出现单个螺旋波,系统从第二类和第三类初态演化将分别出现间歇性全局同步振荡和振荡死亡.当单个神经元处于二周期态时,只有当系统神经元的初相位比较均匀分布时,系统才能自发出现螺旋波,而且出现螺旋波的耦合强度范围大为减少.当神经元处于更高的周期态时,系统一般不容易自发出现螺旋波.这些结果有助于人们了解大脑皮层自发产生螺旋波的机制.     

Mutual and intermittent enhancements of synchronization transitions by autaptic and synaptic delay in scale-free neuron networks

In neural networks, there exist both synaptic delays among different neurons and autaptic self-feedback delays in a neuron itself. In this paper, we study synchronization transitions induced by synaptic and autaptic delays in scale-free neuron networks, mainly exploring how these two time delays affect synchronization transitions induced by each other. It is found that the synchronization transitions induced by synaptic (autaptic) delay are intermittently enhanced when autaptic (synaptic) delay is varied. There are optimal autaptic strength and synaptic coupling strength by which the synchronization transitions induced by autaptic and synaptic delays become strongest. The underlying mechanisms are briefly discussed in terms of the relationships of autaptic delay, synaptic delay, and inter-burst interval. These results show that synaptic and autaptic delays could contribute to each other and enhance synchronization transitions in the neuronal networks. This implies that autaptic and synaptic delays could play a vital role for the information transmission in neural systems.     

Chemical memory erasure; a photochemical model

In this paper we propose an exactly solvable model of a topological insulator defined on a spin- \(\tfrac{1}{2}\) square decorated lattice. Itinerant fermions defined in the framework of the Haldane model interact via the Kitaev interaction with spin- \(\tfrac{1}{2}\) Kitaev sublattice. The presented model, whose ground state is a non-trivial topological phase, is solved exactly. We have found out that various phase transitions without gap closing at the topological phase transition point outline the separate states with different topological numbers. We provide a detailed analysis of the model’s ground-state phase diagram and demonstrate how quantum phase transitions between topological states arise. We have found that the states with both the same and different topological numbers are all separated by the quantum phase transition without gap closing. The transition between topological phases is accompanied by a rearrangement of the spin subsystem’s spectrum from band to flat-band states.     

Theory for photoinduced structural phase transitions

A general concept for photoinduced structural phase transitions is developed in terms of the hidden multistability of the ground state and the proliferations of optically excited states. Taking the ionic→neutral (I - N) phase transition in the organic charge transfer crystal, TTF—CA, as a typical example for this type of transition, we, at first, theoretically show an adiabatic path of this transition, which starts from a single charge transfer exciton in the ionic phase, but finally reaches a neutral domain with a macroscopic size. In connection with this I—N transition, the concept of the initial condition sensitivity is also developed so as to clarify experimentally observed nonlinear characteristics of this material. In the next, using a more simplified model for the many-exciton system, we theoretically study the early time quantum dynamics of the exciton proliferation, which finally results in the domain formation of a large number of ex-citons. For this purpose, we derive a stepwise iterative equation to describe the exciton proliferation, and clarify the origin of the initial condition sensitivity.     

Studies of phase return map and symbolic dynamics in a periodically driven Hodgkin–Huxley neuron

How neuronal spike trains encode external information is a hot topic in neurodynamics studies.In this paper,we investigate the dynamical states of the Hodgkin–Huxley neuron under periodic forcing.Depending on the parameters of the stimulus,the neuron exhibits periodic,quasiperiodic and chaotic spike trains.In order to analyze these spike trains quantitatively,we use the phase return map to describe the dynamical behavior on a one-dimensional(1D)map.According to the monotonicity or discontinuous point of the 1D map,the spike trains are transformed into symbolic sequences by implementing a coarse-grained algorithm—symbolic dynamics.Based on the ordering rules of symbolic dynamics,the parameters of the external stimulus can be measured in high resolution with finite length symbolic sequences.A reasonable explanation for why the nervous system can discriminate or cognize the small change of the external signals in a short time is also presented.     

Bifurcation in neuronal networks with hub structure

We investigate bifurcations in neuronal networks with a hub structure. It is known that hubs play a leading role in characterizing the network dynamical behavior. However, the dynamics of hubs or star-coupled systems is not well understood. Here, we study rather subnetworks with a star-like configuration. This coupled system is an important motif in complex networks. Thus, our study is a basic step for understanding structure formation in large networks. We use the Morris-Lecar neuron with class I and class II excitabilities as a node. Homogeneous (coupling the same class neurons) and heterogeneous (coupling different class neurons) cases are considered for both excitatory and inhibitory coupling. For the homogeneous system class II neurons are suitable for achieving both complete and cluster synchronization in excitatory and inhibitory coupling, respectively. For the heterogeneous system with inhibitory coupling, the class I hub neuron has a wider parameter region of synchronous firings than the class II hub. Moreover, the class I hub neuron with the excitatory synapse gives rise to bifurcations of synchronized states and multi-stability (coexistence of a few different states) is observed.     

Complex transitions between spike,burst or chaos synchronization states in coupled neurons with coexisting bursting patterns

We investigated the synchronization dynamics of a coupled neuronal system composed of two identical Chay model neurons. The Chay model showed coexisting period-1 and period-2 bursting patterns as a parameter and initial values are varied. We simulated multiple periodic and chaotic bursting patterns with non-(NS), burst phase(BS), spike phase(SS),complete(CS), and lag synchronization states. When the coexisting behavior is near period-2 bursting, the transitions of synchronization states of the coupled system follows very complex transitions that begins with transitions between BS and SS, moves to transitions between CS and SS, and to CS. Most initial values lead to the CS state of period-2 bursting while only a few lead to the CS state of period-1 bursting. When the coexisting behavior is near period-1 bursting, the transitions begin with NS, move to transitions between SS and BS, to transitions between SS and CS, and then to CS. Most initial values lead to the CS state of period-1 bursting but a few lead to the CS state of period-2 bursting. The BS was identified as chaos synchronization. The patterns for NS and transitions between BS and SS are insensitive to initial values. The patterns for transitions between CS and SS and the CS state are sensitive to them. The number of spikes per burst of non-CS bursting increases with increasing coupling strength. These results not only reveal the initial value- and parameterdependent synchronization transitions of coupled systems with coexisting behaviors, but also facilitate interpretation of various bursting patterns and synchronization transitions generated in the nervous system with weak coupling strength.     

Detecting synchronization in coupled stochastic ecosystem networks

Instantaneous phase difference, synchronization index and mutual information are considered in order to detect phase transitions, collective behaviours and synchronization phenomena that emerge for different levels of diffusive and reactive activity in stochastic networks. The network under investigation is a spatial 2D lattice which serves as a substrate for Lotka-Volterra dynamics with 3rd order nonlinearities. Kinetic Monte Carlo simulations demonstrate that the system spontaneously organizes into a number of asynchronous local oscillators, when only nearest neighbour interactions are considered. In contrast, the oscillators can be correlated, phase synchronized and completely synchronized when introducing different interactivity rules (diffusive or reactive) for nearby and distant species. The quantitative measures of synchronization show that long distance diffusion coupling induces phase synchronization after a well defined transition point, while long distance reaction coupling induces smeared phase synchronization.     

Dislocation Coupling-Induced Transition of Synchronization in Two-Layer Neuronal Networks

The mutual coupling between neurons in a realistic neuronal system is much complex, and a two-layer neuronal network is designed to investigate the transition of electric activities of neurons. The Hindmarsh-Rose neuron model is used to describe the local dynamics of each neuron, and neurons in the two-layer networks are coupled in dislocated type. The coupling intensity between two-layer networks, and the coupling ratio (Pro), which defines the percentage involved in the coupling in each layer, are changed to observe the synchronization transition of collective behaviors in the two-layer networks. It is found that the two-layer networks of neurons becomes synchronized with increasing the coupling intensity and coupling ratio (Pro) beyond certain thresholds. An ordered wave in the first layer is useful to wake up the rest state in the second layer, or suppress the spatiotemporal state in the second layer under coupling by generating target wave or spiral waves. And the scheme of dislocation coupling can be used to suppress spatiotemporal chaos and excite quiescent neurons.     

Instability and Death of Spiral Wave in a Two-Dimensional Array of Hindmarsh-Rose Neurons

Spiral wave could be observed in the excitable media, the neurons are often excitable within appropriate parameters. The appearance and formation of spiral wave in the cardiac tissue is linked to monomorphic ventricular tachycardia that can denervate into polymorphic tachycardia and ventricular fibrillation. The neuronal system often consists of a large number of neurons with complex connections. In this paper, we theoretically study the transition from spiral wave to spiral turbulence and homogeneous state （death of spiral wave） in two-dimensional array of the Hindmarsh-Rose neuron with completely nearest-neighbor connections. In our numerical studies, a stable rotating spiral wave is developed and selected as the initial state, then the bifurcation parameters are changed to different values to observe the transition from spiral wave to homogeneous state, breakup of spiral wave and weak change of spiral wave, respectively. A statistical factor of synchronization is defined with the mean field theory to analyze the transition from spiral wave to other spatial states, and the snapshots of the membrane potentials of all neurons and time series of mean membrane potentials of all neurons are also plotted to discuss the change of spiral wave. It is found that the sharp changing points in the curve for factor of synchronization vs. bifurcation parameter indicate sudden transition from spiral wave to other states. And the results are independent of the number of neurons we used.     

Hybrid phase transitions of spreading dynamics in multiplex networks

In this paper, we study the spreading dynamics of social behaviors and focus on heterogenous responses of individuals depending on whether they realize the spreading or not. We model the system with a two-layer multiplex network, in which one layer describes the spreading of social behaviors and the other layer describes the diffusion of the awareness about the spreading. We use the susceptible-infected-susceptible (SIS) model to describe the dynamics of an individual if it is unaware of the spreading of the behavior. While when an individual is aware of the spreading of the social behavior its dynamics will follow the threshold model, in which an individual will adopt a behavior only when the fraction of its neighbors who have adopted the behavior is above a certain threshold. We find that such heterogenous reactions can induce intriguing dynamical properties. The dynamics of the whole network may exhibit hybrid phase transitions with the coexistence of continuous phase transition and bi-stable states. Detailed study of how the diffusion of the awareness influences the spreading dynamics of social behavior is provided. The results are supported by theoretical analysis.     

New mechanism for non-trivial intra-molecular vibrational dynamics

We investigate the time evolution process of one selected (initially prepared by optical pumping) vibrational molecular state S, coupled to all other intra-molecular vibrational states R of the same molecule, and also to its environment Q. Molecular states forming the first reservoir R are characterized by a discrete dense spectrum, whereas the environment reservoir Q states form a continuous spectrum. Assuming the equidistant reservoir R states we find the exact analytical solution of the quantum dynamic equations. S-Q and R-Q couplings yield to spontaneous decay of the S and R states, whereas S-R exchange leads to recurrence cycles and Loschmidt echo at frequencies of S-R transitions and double resonances at the interlevel reservoir R transitions. Due to these couplings the system S time evolution is not reduced to a simple exponential relaxation. We predict various regimes of the system S dynamics, ranging from exponential decay to irregular damped oscillations. Namely, we show that there are possible four dynamic regimes of the evolution: (i) independent of the environment Q exponential decay suppressing backward R - S transitions, (ii) Loschmidt echo regime, (iii) incoherent dynamics with multicomponent Loschmidt echo, when the system state is exchanged its energy with many states of the reservoir, (iv) cycle mixing regime, when long time system dynamics looks as a random-like. We suggest applications of our results for interpretation of femtosecond vibration spectra of large molecules and nano-systems.     

液态金属的初始状态对凝固微结构影响的模拟研究

采用分子动力学方法对液态金属Ａｌ在不同的初始状态下，以相同急冷速率凝固的过程进行模拟跟踪研究，发现：在玻璃转变温度Ｔｇ以上（即过冷液态）时，系统的微结构组态情况基本一致，相差甚微；但在Ｔｇ以下时，不同的初始液态微结构对其固态微结构有明显的影响；且在Ｔｇ处各种固态微结构之间的差别发生突发性的变化。这一结果对于深入理解液－固微结构之间的转变关系，具有一定的重要意义。     

Percolation on Interdependent Networks with a Fraction of Antagonistic Interactions

Recently, the percolation transition has been characterized on interacting networks both in presence of interdependent interactions and in presence of antagonistic interactions. Here we characterize the phase diagram of the percolation transition in two Poisson interdependent networks with a percentage q of antagonistic nodes. We show that this system can present a bistability of the steady state solutions, and both discontinuous and continuous phase transitions. In particular, we observe a bistability of the solutions in some regions of the phase space also for a small fraction of antagonistic interactions 0<q<0.4. Moreover, we show that a fraction q>q _c =2/3 of antagonistic interactions is necessary to strongly reduce the region in phase-space in which both networks are percolating. This last result suggests that interdependent networks are robust to the presence of antagonistic interactions. Our approach can be extended to multiple networks, and to complex boolean rules for regulating the percolation phase transition.     

Signal-to-noise analysis of Hopfield neural networks with a formulation of the dynamics in terms of transition probabilities

We study Hopfield neural networks with infinite connectivity using signal-to-noise analysis with a formulation of the dynamics in terms of transition probabilities. We focus our study on models where the strongest effects of the feedback correlations appear. We introduce an analysis of the path probability of one neuron in order to obtain the contribution of all feedback correlations to the dynamics of this neuron. In this way, we obtain a complete theory for dynamics with order parameter equations identical to those obtained with general functional analysis for finite temperature. In the first part of this work, we present our method in the fully connected Little-Hopfield neural network. We obtain, in a simple and direct way, the order parameter equations found by general functional analysis. In the second part, the exposed method is applied to the fully connected Ashkin-Teller neural network. It is shown that the retrieval quality is improved by introducing four-spin couplings. Simulation results support our theoretical predictions.