Mean-field solution of the neural dynamics in a Greenberg-Hastings model with excitatory and inhibitory units.

We present a mean-field solution of the dynamics of a Greenberg-Hastings neural network with both excitatory and inhibitory units. We analyze the dynamical phase transitions that appear in the stationary state as the model parameters are varied. Analytical solutions are compared with numerical simulations of the microscopic model defined on a fully connected network. We found that the stationary state of this system exhibits a first-order dynamical phase transition (with the associated hysteresis) when the fraction of inhibitory units f is smaller than some critical value f_{t}≲1/2, even for a finite system. Moreover, any solution for f<1/2 can be mapped to a solution for purely excitatory systems (f=0). In finite systems, when the system is dominated by inhibition (f>f_{t}), the first-order transition is replaced by a pseudocritical one, namely a continuous crossover between regions of low and high activity that resembles the finite size behavior of a continuous phase transition order parameter. However, in the thermodynamic limit (i.e., infinite-system-size limit), we found that f_{t}→1/2 and the activity for the inhibition dominated case (f≥f_{t}) becomes negligible for any value of the parameters, while the first-order transition between low- and high-activity phases for f

Scale-free correlations in the dynamics of a small (N∼10000) cortical network.

The advent of novel optogenetics technology allows the recording of brain activity with a resolution never seen before. The characterization of these very large data sets offers new challenges as well as unique theory-testing opportunities. Here we discuss whether the spatial and temporal correlations of the collective activity of thousands of neurons are tangled as predicted by the theory of critical phenomena. The analysis shows that both the correlation length ξ and the correlation time τ scale as predicted as a function of the system size. With some peculiarities that we discuss, the analysis uncovers evidence consistent with the view that the large-scale brain cortical dynamics corresponds to critical phenomena.

Finite-size correlation behavior near a critical point: A simple metric for monitoring the state of a neural network.

In this article, a correlation metric κ_{c} is proposed for the inference of the dynamical state of neuronal networks. κ_{C} is computed from the scaling of the correlation length with the size of the observation region, which shows qualitatively different behavior near and away from the critical point of a continuous phase transition. The implementation is first studied on a neuronal network model, where the results of this new metric coincide with those obtained from neuronal avalanche analysis, thus well characterizing the critical state of the network. The approach is further tested with brain optogenetic recordings in behaving mice from a publicly available database. Potential applications and limitations for its use with currently available optical imaging techniques are discussed.

Tricritical behavior in a neural model with excitatory and inhibitory units.

While the support for the relevance of critical dynamics to brain function is increasing, there is much less agreement on the exact nature of the advocated critical point. Thus, a considerable number of theoretical efforts are currently concentrated on which mechanisms and what type(s) of transition can be exhibited by neuronal network models. In that direction, the present work describes the effect of incorporating a fraction of inhibitory neurons on the collective dynamics. As we show, this results in the appearance of a tricritical point for highly connected networks and a nonzero fraction of inhibitory neurons.

Emergent dynamics in a robotic model based on the Caenorhabditis elegans connectome.

We analyze the neural dynamics and their relation with the emergent actions of a robotic vehicle that is controlled by a neural network numerical simulation based on the nervous system of the nematode Caenorhabditis elegans. The robot interacts with the environment through a sensor that transmits the information to sensory neurons, while motor neurons outputs are connected to wheels. This is enough to allow emergent robot actions in complex environments, such as avoiding collisions with obstacles. Working with robotic models makes it possible to simultaneously keep track of the dynamics of all the neurons and also register the actions of the robot in the environment in real time, while avoiding the complex technicalities of simulating a real environment. This allowed us to identify several relevant features of the neural dynamics associated with the emergent actions of the robot, some of which have already been observed in biological worms. These results suggest that some basic aspects of behaviors observed in living beings are determined by the underlying structure of the associated neural network.

Box scaling as a proxy of finite size correlations.

The scaling of correlations as a function of size provides important hints to understand critical phenomena on a variety of systems. Its study in biological structures offers two challenges: usually they are not of infinite size, and, in the majority of cases, dimensions can not be varied at will. Here we discuss how finite-size scaling can be approximated in an experimental system of fixed and relatively small extent, by computing correlations inside of a reduced field of view of various widths (we will refer to this procedure as "box-scaling"). A relation among the size of the field of view, and measured correlation length, is derived at, and away from, the critical regime. Numerical simulations of a neuronal network, as well as the ferromagnetic 2D Ising model, are used to verify such approximations. Numerical results support the validity of the heuristic approach, which should be useful to characterize relevant aspects of critical phenomena in biological systems.

Similar local neuronal dynamics may lead to different collective behavior.

This report is concerned with the relevance of the microscopic rules that implement individual neuronal activation, in determining the collective dynamics, under variations of the network topology. To fix ideas we study the dynamics of two cellular automaton models, commonly used, rather in-distinctively, as the building blocks of large-scale neuronal networks. One model, due to Greenberg and Hastings (GH), can be described by evolution equations mimicking an integrate-and-fire process, while the other model, due to Kinouchi and Copelli (KC), represents an abstract branching process, where a single active neuron activates a given number of postsynaptic neurons according to a prescribed "activity" branching ratio. Despite the apparent similarity between the local neuronal dynamics of the two models, it is shown that they exhibit very different collective dynamics as a function of the network topology. The GH model shows qualitatively different dynamical regimes as the network topology is varied, including transients to a ground (inactive) state, continuous and discontinuous dynamical phase transitions. In contrast, the KC model only exhibits a continuous phase transition, independently of the network topology. These results highlight the importance of paying attention to the microscopic rules chosen to model the interneuronal interactions in large-scale numerical simulations, in particular when the network topology is far from a mean-field description. One such case is the extensive work being done in the context of the Human Connectome, where a wide variety of types of models are being used to understand the brain collective dynamics.

Controlling a complex system near its critical point via temporal correlations.

Many complex systems exhibit large fluctuations both across space and over time. These fluctuations have often been linked to the presence of some kind of critical phenomena, where it is well known that the emerging correlation functions in space and time are closely related to each other. Here we test whether the time correlation properties allow systems exhibiting a phase transition to self-tune to their critical point. We describe results in three models: the 2D Ising ferromagnetic model, the 3D Vicsek flocking model and a small-world neuronal network model. We demonstrate that feedback from the autocorrelation function of the order parameter fluctuations shifts the system towards its critical point. Our results rely on universal properties of critical systems and are expected to be relevant to a variety of other settings.

Universal and nonuniversal neural dynamics on small world connectomes: A finite-size scaling analysis.

Evidence of critical dynamics has been found recently in both experiments and models of large-scale brain dynamics. The understanding of the nature and features of such a critical regime is hampered by the relatively small size of the available connectome, which prevents, among other things, the determination of its associated universality class. To circumvent that, here we study a neural model defined on a class of small-world networks that share some topological features with the human connectome. We find that varying the topological parameters can give rise to a scale-invariant behavior either belonging to the mean-field percolation universality class or having nonuniversal critical exponents. In addition, we find certain regions of the topological parameter space where the system presents a discontinuous, i.e., noncritical, dynamical phase transition into a percolated state. Overall, these results shed light on the interplay of dynamical and topological roots of the complex brain dynamics.

Three-state model with competing antiferromagnetic and pairing interactions.

Motivated by the rich phase diagram of the high-temperature superconductors, we introduce a pseudospin model with three state variables which can be interpreted as two states (spin ±1/2) particles and holes. The Hamiltonian has a term which favors antiferromagnetism and an additional competing interaction which favors bonding between pairs of antiparallel spins mediated by holes. For low concentration of holes the dominant interaction between particles has antiferromagnetic character, leading to an antiferromagnetic phase in the temperature-hole concentration phase diagram, qualitatively similar to the antiferromagnetic phase of doped Mott insulators. For growing concentration of holes antiferromagnetic order is weakened and a phase with a different kind of order mediated by holes appears. This last phase has the form of a dome in the T-hole concentration plane. The whole phase diagram resembles those of some families of high-T_{c} superconductors. We compute the phase diagram in the mean-field approximation and characterize the different phase transitions through Monte Carlo simulations.

Nrf2 stabilization prevents critical oxidative damage in Down syndrome cells.

Mounting evidence implicates chronic oxidative stress as a critical driver of the aging process. Down syndrome (DS) is characterized by a complex phenotype, including early senescence. DS cells display increased levels of reactive oxygen species (ROS) and mitochondrial structural and metabolic dysfunction, which are counterbalanced by sustained Nrf2-mediated transcription of cellular antioxidant response elements (ARE). Here, we show that caspase 3/PKCδdependent activation of the Nrf2 pathway in DS and Dp16 (a mouse model of DS) cells is necessary to protect against chronic oxidative damage and to preserve cellular functionality. Mitochondria-targeted catalase (mCAT) significantly reduced oxidative stress, restored mitochondrial structure and function, normalized replicative and wound healing capacity, and rendered the Nrf2-mediated antioxidant response dispensable. These results highlight the critical role of Nrf2/ARE in the maintenance of DS cell homeostasis and validate mitochondrial-specific interventions as a key aspect of antioxidant and antiaging therapies.

Mitochondrial network complexity emerges from fission/fusion dynamics.

Mitochondrial networks exhibit a variety of complex behaviors, including coordinated cell-wide oscillations of energy states as well as a phase transition (depolarization) in response to oxidative stress. Since functional and structural properties are often interwinded, here we characterized the structure of mitochondrial networks in mouse embryonic fibroblasts using network tools and percolation theory. Subsequently we perturbed the system either by promoting the fusion of mitochondrial segments or by inducing mitochondrial fission. Quantitative analysis of mitochondrial clusters revealed that structural parameters of healthy mitochondria laid in between the extremes of highly fragmented and completely fusioned networks. We confirmed our results by contrasting our empirical findings with the predictions of a recently described computational model of mitochondrial network emergence based on fission-fusion kinetics. Altogether these results offer not only an objective methodology to parametrize the complexity of this organelle but also support the idea that mitochondrial networks behave as critical systems and undergo structural phase transitions.

Magnetization reversal in mixed ferrite-chromite perovskites with non magnetic cation on the A-site.

In this work, we have performed Monte Carlo simulations in a classical model for RFe1-x Cr x O3 with R = Y and Lu, comparing the numerical simulations with experiments and mean field calculations. In the analyzed compounds, the antisymmetric exchange or Dzyaloshinskii-Moriya (DM) interaction induced a weak ferromagnetism due to a canting of the antiferromagnetically ordered spins. This model is able to reproduce the magnetization reversal (MR) observed experimentally in a field cooling process for intermediate x values and the dependence with x of the critical temperatures. We also analyzed the conditions for the existence of MR in terms of the strength of DM interactions between Fe(3+) and Cr(3+) ions with the x values variations.

Finite-temperature phase diagram of ultrathin magnetic films without external fields.

We analyze the finite-temperature phase diagram of ultrathin magnetic films by introducing a mean-field theory, valid in the low-anisotropy regime, i.e., close to the spin reorientation transition. The theoretical results are compared with Monte Carlo simulations carried out on a microscopic Heisenberg model. Connections between the finite-temperature behavior and the ground-state properties of the system are established. Several properties of the stripe pattern, such as the presence of canted states, the stripe width variation phenomenon, and the associated magnetization profiles, are also analyzed.

Nonequilibrium structures and slow dynamics in a two-dimensional spin system with competing long-range and short-range interactions.

We present a lattice spin model that mimics a system of interacting particles through a short-range repulsive potential and a long-range attractive power-law decaying potential. We perform a detailed analysis of the general equilibrium phase diagram of the model at finite temperature, showing that the only possible equilibrium phases are the ferromagnetic and the antiferromagnetic ones. We then study the nonequilibrium behavior of the model after a quench to subcritical temperatures, in the antiferromagnetic region of the phase diagram region, where the pair interaction potential behaves in the same qualitative way as in a Lennard-Jones gas. We find that even in the absence of quenched disorder or geometric frustration, the competition between interactions gives rise to nonequilibrium disordered structures at low enough temperatures that strongly slow down the relaxation of the system. This nonequilibrium state presents several features characteristic of glassy systems such as subaging, nontrivial fuctuation dissipation relations, and possible logarithmic growth of free-energy barriers to coarsening.

Emergent self-organized complex network topology out of stability constraints.

Although most networks in nature exhibit complex topologies, the origins of such complexity remain unclear. We propose a general evolutionary mechanism based on global stability. This mechanism is incorporated into a model of a growing network of interacting agents in which each new agent's membership in the network is determined by the agent's effect on the network's global stability. It is shown that out of this stability constraint complex topological properties emerge in a self-organized manner, offering an explanation for their observed ubiquity in biological networks.

Nonequilibrium characterization of spinodal points using short time dynamics.

Although intuitively appealing, the concept of spinodal is rigorously defined only in systems with infinite range interactions (mean-field systems). In short-range systems, a pseudospinodal can be defined by extrapolation of metastable measurements, but the point itself is not reachable because it lies beyond the metastability limit. In this work we show that a sensible definition of spinodal points can be obtained through the short time dynamical behavior of the system deep inside the metastable phase by looking for a point where the system shows critical behavior. We show that spinodal points obtained by this method agree both with the thermodynamical spinodal point in mean-field systems and with the pseudospinodal point obtained by extrapolation of metaequilibrium behavior in short-range systems. With this definition, a practical determination can be achieved without regard for equilibration issues.

Comparable ecological dynamics underlie early cancer invasion and species dispersal, involving self-organizing processes.

Occupancy of new habitats through dispersion is a central process in nature. In particular, long-distance dispersal is involved in the spread of species and epidemics, although it has not been previously related with cancer invasion, a process that involves cell spreading to tissues far away from the primary tumour. Using simulations and real data we show that the early spread of cancer cells is similar to the species individuals spread and we suggest that both processes are represented by a common spatio-temporal signature of long-distance dispersal and subsequent local proliferation. This signature is characterized by a particular fractal geometry of the boundaries of patches generated, and a power-law scaled, disrupted patch size distribution. In contrast, invasions involving only dispersal but not subsequent proliferation ("physiological invasions") like trophoblast cells invasion during normal human placentation did not show the patch size power-law pattern. Our results are consistent under different temporal and spatial scales, and under different resolution levels of analysis. We conclude that the scaling properties are a hallmark and a direct result of long-distance dispersal and proliferation, and that they could reflect homologous ecological processes of population self-organization during cancer and species spread. Our results are significant for the detection of processes involving long-range dispersal and proliferation like cancer local invasion and metastasis, biological invasions and epidemics, and for the formulation of new cancer therapeutical approaches.

Interplay between coarsening and nucleation in an Ising model with dipolar interactions.

We study the dynamical behavior of a square lattice Ising model with exchange and dipolar interactions by means of Monte Carlo simulations. After a sudden quench to low temperatures, we find that the system may undergo a coarsening process where stripe phases with different orientations compete, or alternatively it can relax initially to a metastable nematic phase and then decay to the equilibrium stripe phase through nucleation. We measure the distribution of equilibration times for both processes and compute their relative probability of occurrence as a function of temperature and system size. This peculiar relaxation mechanism is due to the strong metastability of the nematic phase, which goes deep into the low-temperature stripe phase. We also measure quasiequilibrium autocorrelations in a wide range of temperatures. They show a distinct decay to a plateau that we identify as due to a finite fraction of frozen spins in the nematic phase. We find indications that the plateau is a finite-size effect. Relaxation times as a function of temperature in the metastable region show super-Arrhenius behavior, suggesting a possible glassy behavior of the system at low temperatures.

Long-term ordering kinetics of the two-dimensional q-state Potts model.

We studied the nonequilibrium dynamics of the q-state Potts model in the square lattice, after a quench to subcritical temperatures. By means of a continuous time Monte Carlo algorithm (nonconserved order parameter dynamics) we analyzed the long term behavior of the energy and relaxation time for a wide range of quench temperatures and system sizes. For q>4 we found the existence of different dynamical regimes, according to quench temperature range. At low (but finite) temperatures and very long times the Lifshitz-Allen-Cahn domain growth behavior is interrupted with finite probability when the system gets stuck in highly symmetric nonequilibrium metastable states, which induce activation in the domain growth, in agreement with early predictions of Lifshitz [JETP 42, 1354 (1962)]. Moreover, if the temperature is very low, the system always gets stuck at short times in highly disordered metastable states with finite lifetime, which have been recently identified as glassy states. The finite size scaling properties of the different relaxation times involved, as well as their temperature dependency, are analyzed in detail.