Critical phenomena and noise-induced phase transitions in neuronal networks

Controlling permeation in electrically deforming liquid crystal network films: A dynamical Landau theory

Liquid crystal networks exploit the coupling between the responsivity of liquid crystalline mesogens, e.g., to electric fields, and the (visco)elastic properties of a polymer network. Because of this, these materials have been put forward for a wide array of applications, including responsive surfaces such as artificial skins and membranes. For such applications, the desired functional response must generally be realized under strict geometrical constraints, such as provided by supported thin films. To model such settings, we present a dynamical, spatially heterogeneous Landau-type theory for electrically actuated liquid crystal network films. We find that the response of the liquid crystal network permeates the film from top to bottom, and illustrate how this affects the timescale associated with macroscopic deformation. Finally, by linking our model parameters to experimental quantities, we suggest that the permeation rate can be controlled by varying the aspect ratio of the mesogens and their degree of orientational order when crosslinked into the polymer network, for which we predict a single optimum. Our results contribute specifically to the rational design of future applications involving transport or on-demand release of molecular cargo in liquid crystal network films.

Dynamical Landau-de Gennes theory for electrically-responsive liquid crystal networks

Liquid crystal networks combine the orientational order of liquid crystals with the elastic properties of polymer networks, leading to a vast application potential in the field of responsive coatings, e.g., for haptic feedback, self-cleaning surfaces, and static and dynamic pattern formation. Recent experimental work has further paved the way toward such applications by realizing the fast and reversible surface modulation of a liquid crystal network coating upon in-plane actuation with an AC electric field [Liu, Tito, and Broer, Nat. Commun. 8, 1526 (2017)10.1038/s41467-017-01448-w]. Here, we construct a Landau-type theory for electrically-responsive liquid crystal networks and perform molecular dynamics simulations to explain the findings of these experiments and inform on rational design strategies. Qualitatively, the theory agrees with our simulations and reproduces the salient experimental features. We also provide a set of testable predictions: the aspect ratio of the nematogens, their initial orientational order when cross-linked into the polymer network, and the cross-linking fraction of the network all increase the plasticization time required for the film to macroscopically deform. We demonstrate that the dynamic response to oscillating electric fields is characterized by two resonances, which can likewise be influenced by varying these parameters, providing an experimental handle to fine-tune device design.

Aperiodic stochastic resonance in neural information processing with Gaussian colored noise

The aim of this paper is to explore the phenomenon of aperiodic stochastic resonance in neural systems with colored noise. For nonlinear dynamical systems driven by Gaussian colored noise, we prove that the stochastic sample trajectory can converge to the corresponding deterministic trajectory as noise intensity tends to zero in mean square, under global and local Lipschitz conditions, respectively. Then, following forbidden interval theorem we predict the phenomenon of aperiodic stochastic resonance in bistable and excitable neural systems. Two neuron models are further used to verify the theoretical prediction. Moreover, we disclose the phenomenon of aperiodic stochastic resonance induced by correlation time and this finding suggests that adjusting noise correlation might be a biologically more plausible mechanism in neural signal processing.

Discrete Langevin machine: Bridging the gap between thermodynamic and neuromorphic systems

A formulation of Langevin dynamics for discrete systems is derived as a class of generic stochastic processes. The dynamics simplify for a two-state system and suggest a network architecture which is implemented by the Langevin machine. The Langevin machine represents a promising approach to compute successfully quantitative exact results of Boltzmann distributed systems by LIF neurons. Besides a detailed introduction of the dynamics, different simplified models of a neuromorphic hardware system are studied with respect to a control of emerging sources of errors.

Statistical Approaches for the Analysis of Dependency Among Neurons Under Noise

Neuronal noise is a major factor affecting the communication between coupled neurons. In this work, we propose a statistical toolset to infer the coupling between two neurons under noise. We estimate these statistical dependencies from data which are generated by a coupled Hodgkin-Huxley (HH) model with additive noise. To infer the coupling using observation data, we employ copulas and information-theoretic quantities, such as the mutual information (MI) and the transfer entropy (TE). Copulas and MI between two variables are symmetric quantities, whereas TE is asymmetric. We demonstrate the performances of copulas and MI as functions of different noise levels and show that they are effective in the identification of the interactions due to coupling and noise. Moreover, we analyze the inference of TE values between neurons as a function of noise and conclude that TE is an effective tool for finding out the direction of coupling between neurons under the effects of noise.

Proper spatial heterogeneities expand the regime of scale-free behavior in a lattice of excitable elements

Signatures of criticality, such as power law scaling of observables, have been empirically found in a plethora of real-life settings, including biological systems. The presence of critical states is believed to have many functional advantages and is associated with optimal operational abilities. Typically, critical dynamics arises in the proximity of phase transition points between absorbing disordered states (subcriticality) and ordered active regimes (supercriticality) and requires a high degree of fine tuning to emerge, which is unlikely to occur in real biological systems. In the present study we propose a rather simple, and biologically relevant mechanism that profoundly expands the critical-like region. In particular, by means of numerical simulation we show that incorporating spatial heterogeneities into the square lattice of map-based excitable oscillators broadens the parameter space in which the distribution of excitation wave sizes follows closely a power law. Most importantly, this behavior is only observed if the spatial profile exhibits intermediate-sized patches with similar excitability levels, whereas for large and small spatial clusters only marginal widening of the critical state is detected. Furthermore, it turned out that the presence of spatial disorder in general amplifies the size of excitation waves, whereby the relatively highest contributions are observed in the proximity of the critical point. We argue that the reported mechanism is of particular importance for excitable systems with local interactions between individual elements.

Distinct dynamical behavior in Erdős-Rényi networks, regular random networks, ring lattices, and all-to-all neuronal networks

Neuronal network dynamics depends on network structure. In this paper we study how network topology underpins the emergence of different dynamical behaviors in neuronal networks. In particular, we consider neuronal network dynamics on Erdős-Rényi (ER) networks, regular random (RR) networks, ring lattices, and all-to-all networks. We solve analytically a neuronal network model with stochastic binary-state neurons in all the network topologies, except ring lattices. Given that apart from network structure, all four models are equivalent, this allows us to understand the role of network structure in neuronal network dynamics. While ER and RR networks are characterized by similar phase diagrams, we find strikingly different phase diagrams in the all-to-all network. Neuronal network dynamics is not only different within certain parameter ranges, but it also undergoes different bifurcations (with a richer repertoire of bifurcations in ER and RR compared to all-to-all networks). This suggests that local heterogeneity in the ratio between excitation and inhibition plays a crucial role on emergent dynamics. Furthermore, we also observe one subtle discrepancy between ER and RR networks, namely, ER networks undergo a neuronal activity jump at lower noise levels compared to RR networks, presumably due to the degree heterogeneity in ER networks that is absent in RR networks. Finally, a comparison between network oscillations in RR networks and ring lattices shows the importance of small-world properties in sustaining stable network oscillations.

A Basic Phase Diagram of Neuronal Dynamics

The extreme complexity of the brain has attracted the attention of neuroscientists and other researchers for a long time. More recently, the neuromorphic hardware has matured to provide a new powerful tool to study neuronal dynamics. Here, we study neuronal dynamics using different settings on a neuromorphic chip built with flexible parameters of neuron models. Our unique setting in the network of leaky integrate-and-fire (LIF) neurons is to introduce a weak noise environment. We observed three different types of collective neuronal activities, or phases, separated by sharp boundaries, or phase transitions. From this, we construct a rudimentary phase diagram of neuronal dynamics and demonstrate that a noise-induced chaotic phase (N-phase), which is dominated by neuronal avalanche activity (intermittent aperiodic neuron firing), emerges in the presence of noise and its width grows with the noise intensity. The dynamics can be manipulated in this N-phase. Our results and comparison with clinical data is consistent with the literature and our previous work showing that healthy brain must reside in the N-phase. We argue that the brain phase diagram with further refinement may be used for the diagnosis and treatment of mental disease and also suggest that the dynamics may be manipulated to serve as a means of new information processing (e.g., for optimization). Neuromorphic chips, similar to the one we used but with a variety of neuron models, may be used to further enhance the understanding of human brain function and accelerate the development of neuroscience research.

Postoperative seizure outcome-guided machine learning for interictal electrocorticography in neocortical epilepsy

The objective of this study was to introduce a new machine learning guided by outcome of resective epilepsy surgery defined as the presence/absence of seizures to improve data mining for interictal pathological activities in neocortical epilepsy. Electrocorticographies for 39 patients with medically intractable neocortical epilepsy were analyzed. We separately analyzed 38 frequencies from 0.9 to 800 Hz including both high-frequency activities and low-frequency activities to select bands related to seizure outcome. An automatic detector using amplitude-duration-number thresholds was used. Interictal electrocorticography data sets of 8 min for each patient were selected. In the first training data set of 20 patients, the automatic detector was optimized to best differentiate the seizure-free group from not-seizure-free-group based on ranks of resection percentages of activities detected using a genetic algorithm. The optimization was validated in a different data set of 19 patients. There were 16 (41%) seizure-free patients. The mean follow-up duration was 21 ± 11 mo (range, 13-44 mo). After validation, frequencies significantly related to seizure outcome were 5.8, 8.4-25, 30, 36, 52, and 75 among low-frequency activities and 108 and 800 Hz among high-frequency activities. Resection for 5.8, 8.4-25, 108, and 800 Hz activities consistently improved seizure outcome. Resection effects of 17-36, 52, and 75 Hz activities on seizure outcome were variable according to thresholds. We developed and validated an automated detector for monitoring interictal pathological and inhibitory/physiological activities in neocortical epilepsy using a data-driven approach through outcome-guided machine learning. NEW & NOTEWORTHY Outcome-guided machine learning based on seizure outcome was used to improve detections for interictal electrocorticographic low- and high-frequency activities. This method resulted in better separation of seizure outcome groups than others reported in the literature. The automatic detector can be trained without human intervention and no prior information. It is based only on objective seizure outcome data without relying on an expert's manual annotations. Using the method, we could find and characterize pathological and inhibitory activities.

Structural stability of interaction networks against negative external fields

We explore structural stability of weighted and unweighted networks of positively interacting agents against a negative external field. We study how the agents support the activity of each other to confront the negative field, which suppresses the activity of agents and can lead to collapse of the whole network. The competition between the interactions and the field shape the structure of stable states of the system. In unweighted networks (uniform interactions) the stable states have the structure of k-cores of the interaction network. The interplay between the topology and the distribution of weights (heterogeneous interactions) impacts strongly the structural stability against a negative field, especially in the case of fat-tailed distributions of weights. We show that apart from critical slowing down there is also a critical change in the system structure that precedes the network collapse. The change can serve as an early warning of the critical transition. To characterize changes of network structure we develop a method based on statistical analysis of the k-core organization and so-called "corona" clusters belonging to the k-cores.

Neuronal network model of interictal and recurrent ictal activity

We propose a neuronal network model which undergoes a saddle node on an invariant circle bifurcation as the mechanism of the transition from the interictal to the ictal (seizure) state. In the vicinity of this transition, the model captures important dynamical features of both interictal and ictal states. We study the nature of interictal spikes and early warnings of the transition predicted by this model. We further demonstrate that recurrent seizures emerge due to the interaction between two networks.

Approximate-master-equation approach for the Kinouchi-Copelli neural model on networks

In this work, we use the approximate-master-equation approach to study the dynamics of the Kinouchi-Copelli neural model on various networks. By categorizing each neuron in terms of its state and also the states of its neighbors, we are able to uncover how the coupled system evolves with respective to time by directly solving a set of ordinary differential equations. In particular, we can easily calculate the statistical properties of the time evolution of the network instantaneous response, the network response curve, the dynamic range, and the critical point in the framework of the approximate-master-equation approach. The possible usage of the proposed theoretical approach to other spreading phenomena is briefly discussed.

Optimal Information Representation and Criticality in an Adaptive Sensory Recurrent Neuronal Network

Recurrent connections play an important role in cortical function, yet their exact contribution to the network computation remains unknown. The principles guiding the long-term evolution of these connections are poorly understood as well. Therefore, gaining insight into their computational role and into the mechanism shaping their pattern would be of great importance. To that end, we studied the learning dynamics and emergent recurrent connectivity in a sensory network model based on a first-principle information theoretic approach. As a test case, we applied this framework to a model of a hypercolumn in the visual cortex and found that the evolved connections between orientation columns have a "Mexican hat" profile, consistent with empirical data and previous modeling work. Furthermore, we found that optimal information representation is achieved when the network operates near a critical point in its dynamics. Neuronal networks working near such a phase transition are most sensitive to their inputs and are thus optimal in terms of information representation. Nevertheless, a mild change in the pattern of interactions may cause such networks to undergo a transition into a different regime of behavior in which the network activity is dominated by its internal recurrent dynamics and does not reflect the objective input. We discuss several mechanisms by which the pattern of interactions can be driven into this supercritical regime and relate them to various neurological and neuropsychiatric phenomena.

Brain Performance versus Phase Transitions

We here illustrate how a well-founded study of the brain may originate in assuming analogies with phase-transition phenomena. Analyzing to what extent a weak signal endures in noisy environments, we identify the underlying mechanisms, and it results a description of how the excitability associated to (non-equilibrium) phase changes and criticality optimizes the processing of the signal. Our setting is a network of integrate-and-fire nodes in which connections are heterogeneous with rapid time-varying intensities mimicking fatigue and potentiation. Emergence then becomes quite robust against wiring topology modification—in fact, we considered from a fully connected network to the Homo sapiens connectome—showing the essential role of synaptic flickering on computations. We also suggest how to experimentally disclose significant changes during actual brain operation.

Noise-enhanced nonlinear response and the role of modular structure for signal detection in neuronal networks

We show that sensory noise can enhance the nonlinear response of neuronal networks, and when delivered together with a weak signal, it improves the signal detection by the network. We reveal this phenomenon in neuronal networks that are in a dynamical state preceding a saddle-node bifurcation corresponding to the appearance of sustained network oscillations. In this state, even a weak subthreshold pulse can evoke a large-amplitude oscillation of neuronal activity. The signal-to-noise ratio reaches a maximum at an optimum level of sensory noise, manifesting stochastic resonance (SR) at the population level. We demonstrate SR by use of simulations and numerical integration of rate equations in a cortical model. Using this model, we mimic the experiments of Gluckman et al. [Phys. Rev. Lett. 77, 4098 (1996)PRLTAO0031-900710.1103/PhysRevLett.77.4098] that have given evidence of SR in mammalian brain. We also study neuronal networks in which neurons are grouped in modules and every module works in the regime of SR. We find that even a few modules can strongly enhance the reliability of signal detection in comparison with the case when a modular organization is absent.