REVISAR
Theoretical foundations of studying
criticality in the brain
Yang Tian1,2
, Zeren Tan3, Hedong Hou4, Guoqi Li5,6, Aohua Cheng7,
Yike Qiu7, Kangyu Weng7, Chun Chen1, y Pei Sun1
1Department of Psychology & Tsinghua Laboratory of Brain and Intelligence, Universidad de Tsinghua, Beijing, Porcelana
2Laboratorio de Computación Avanzada y Almacenamiento, Instituto Central de Investigaciones, 2012 Laboratorios,
Huawei Technologies Co.. Limitado., Beijing, Porcelana
3Institute for Interdisciplinary Information Science, Universidad de Tsinghua, Beijing, Porcelana
4UFR de Mathématiques, Université de Paris, París, Francia
5Institute of Automation, Chinese Academy of Science, Beijing, Porcelana
6University of Chinese Academy of Science, Beijing, Porcelana
7Tsien Excellence in Engineering Program, School of Aerospace Engineering, Universidad de Tsinghua, Beijing, Porcelana
Palabras clave: Nonequilibrium criticality, Neural avalanches, Neural dynamics, Directed percolation
ABSTRACTO
Criticality is hypothesized as a physical mechanism underlying efficient transitions between
cortical states and remarkable information-processing capacities in the brain. Mientras
considerable evidence generally supports this hypothesis, nonnegligible controversies persist
regarding the ubiquity of criticality in neural dynamics and its role in information processing.
Validity issues frequently arise during identifying potential brain criticality from empirical data.
Además, the functional benefits implied by brain criticality are frequently misconceived or
unduly generalized. These problems stem from the nontriviality and immaturity of the physical
theories that analytically derive brain criticality and the statistic techniques that estimate brain
criticality from empirical data. To help solve these problems, we present a systematic review
and reformulate the foundations of studying brain criticality, eso es, ordinary criticality (OC),
quasi-criticality (qC), self-organized criticality (SOC), and self-organized quasi-criticality
(SOqC), using the terminology of neuroscience. We offer accessible explanations of the
physical theories and statistical techniques of brain criticality, providing step-by-step
derivations to characterize neural dynamics as a physical system with avalanches. Nosotros
summarize error-prone details and existing limitations in brain criticality analysis and suggest
possible solutions. Además, we present a forward-looking perspective on how optimizing
the foundations of studying brain criticality can deepen our understanding of various
neuroscience questions.
RESUMEN DEL AUTOR
The brain criticality hypothesis is one of the most focused and controversial topics in
neuroscience and biophysics. This research develops a unified framework to reformulate the
physics theories of four basic types of brain criticality, ordinary criticality (OC), quasi-criticality
(qC), self-organized criticality (SOC), and self-organized quasi-criticality (SOqC), into more
accessible and neuroscience-related forms. For the statistic techniques used to validate the
brain criticality hypothesis, we also present comprehensive explanations of them, summarize
their error-prone details, and suggest possible solutions. This framework may help resolve
un acceso abierto
diario
Citación: tian, y., Broncearse, Z., Hou, h., li, GRAMO.,
cheng, A., Qiu, y., Weng, K., Chen, C.,
& Sol, PAG. (2022). Teórico
foundations of studying criticality in
the brain. Neurociencia en red, 6(4),
1148–1185. https://doi.org/10.1162/netn
_a_00269
DOI:
https://doi.org/10.1162/netn_a_00269
Recibió: 2 Abril 2022
Aceptado: 12 Julio 2022
Conflicto de intereses: Los autores tienen
declaró que no hay intereses en competencia
existir.
Autor correspondiente:
Pei Sol
peisun@tsinghua.edu.cn
Editor de manejo:
Olaf Sporns
Derechos de autor: © 2022
Instituto de Tecnología de Massachusetts
Publicado bajo Creative Commons
Atribución 4.0 Internacional
(CC POR 4.0) licencia
La prensa del MIT
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
/
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
potential controversies in studying the brain criticality hypothesis, especially those arising from
the misconceptions about the theoretical foundations of brain criticality.
Neurophysics:
A branch of biophysics that develops
and uses physics theories to study the
neural system.
Conectividad cerebral:
Refers to the anatomical connectivity
formed by synaptic connections and
the functional connectivity formed
by dynamic interactions among
neuronas.
Statistical physics:
A branch of theoretical physics that
develops mathematical theories
of the characterization and
approximations of large populations
with inherently stochastic natures.
Teoría de la percolación:
Una teoría física que caracteriza
fenómenos críticos y fase
transiciones desde un modelo probabilístico y
perspectiva geométrica.
Hopfield network:
An Ising model of a neural network,
which serves as a content-
addressable memory system with
binary nodes or continuous
variables.
Information thermodynamics:
A branch of statistical physics that
develops mathematical theories of
the exchange between information
quantities and thermodynamic
quantities.
INTRODUCCIÓN
Neuroscience is dawning upon revealing physics foundations of the brain (Abbott, 2008). Ever
since the 1970s, the term neurophysics has been suggested as a term to indicate the essential
role of physics in understanding the brain (Scott, 1977). More recently, substantial progress has
been accomplished in studying brain connectivity and brain functions with statistical physics
teorías (Lynn & bassett, 2019).
For brain connectivity, physics provides insights for its emergence, organización, y
evolution. Random graphs (Betzel et al., 2016; Betzel & bassett, 2017a), percolation (Breskin,
Soriano, Moses, & Tlusty, 2006; Guo et al., 2021), and other physics theories of correlated
sistemas (Haimovici, Tagliazucchi, Balenzuela, & Chialvo, 2013; Lobo, 2005) are applied to
reveal the underlying mechanisms accounting for the origins of brain network properties.
Complex network theories act as the foundation of characterizing brain connectivity organi-
zational features (p.ej., comunidad; Betzel & bassett, 2017b; Betzel, Medalla, & bassett, 2018;
Khambhati, Sizemore, Betzel, & bassett, 2018), hub (decoración, Tononi, boly, & Kringelbach,
2015; Gong et al., 2009), and small-world (bullmore & despreciar, 2012; Deco et al., 2015; struc-
turas) and embedding attributes into physical space (Bassett y cols., 2010; Kaiser & Hilgetag,
2006). Network evolution driven by neural plasticity helps to explain the dynamics of brain
connectivity structures during information processing (Del Pozo et al., 2021; Galván, 2010;
montesco, Dayán, & Sejnowski, 1996; Roberto & Vignoud, 2021; Song, Molinero, & Abbott,
2000). For brain functions, physics presents possible explanations for the origin of information
processing capacities from collective neural activities. From single neuron dynamics models
(Gerstner, Kistler, Beneficio, & paninski, 2014), stochastic network models of neural populations
and circuits (tian, li, & Sol, 2021; tian & Sol, 2021), mean-field neural mass models of brain
regiones (David & Friston, 2003; Touboul, Wendling, Chauvel, & Faugeras, 2011), eventually to
models of entire brain networks (Hopfield, 1982; Schneidman, Berry, Segev, & Bialek, 2006),
important efforts have been devoted to characterize information-processing-related neural
dynamics across different scales. Networks with memory capacities (p.ej., Hopfield networks;
Tyulmankov, Fang, Vadaparty, & Cual, 2021), which are equivalent to Ising models under spe-
cific conditions (Lynn & bassett, 2019), have been applied to study neural information storage
and recall (Haldeman & ruega, 2005; Krotov & Hopfield, 2020), adaptation to environment
cambios (Shew et al., 2015), information transmission optimization (ruega & Plenz, 2003),
dynamic range maximization (Kinouchi & Copelli, 2006; mostrar, Cual, Petermann, roy, &
Plenz, 2009), and neural computation power (Bertschinger & Natschläger, 2004). Estos
models are further related to maximum entropy models (p.ej., specific fine-tuned Ising models)
that predict long-range correlations observed among neurons (Ganmor, Segev, & Schneidman,
2011; Schneidman et al., 2006). Además, general theories of free-energy principle (Friston,
2009, 2010; Guevara, 2021) and information thermodynamics (Capolupo, Hombre libre, &
Vitiello, 2013; Cuello & grifo, 2015; Sartori, Granjero, Sotavento, & Horowitz, 2014; tian &
Sol, 2022) are suggested as the unified foundations of perception, acción, and learning in
the brain.
If one needs to specify one of the most focused and controversial topics among all the
works mentioned above, brain criticality may be a potential candidate (ruega & Hora,
Neurociencia en red
1149
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Dynamic range:
The intensity or amplitude range of
stimulus inputs that is encoded by the
given neural dynamics.
2012). The hypothesis of the critical brain has received increasing attention in recent decades,
serving as a possible mechanism underlying various intriguing but elusive phenomena in the
cerebro. In light of our limited understanding of the complex nature of collective neural dynam-
circuitos integrados, these phenomena include, to name a few, efficient transitions between cortical states
(Fontenele et al., 2019), maximal dynamic ranges of neural responses (Antonopoulos, 2016;
Gautama, Hoang, McClanahan, grady, & mostrar, 2015; Kinouchi & Copelli, 2006; Shew et al.,
2009), optimized information transmission and representation (Antonopoulos, 2016; X. li &
Pequeño, 2012; mostrar, Cual, Yu, roy, & Plenz, 2011), and numerous other issues concerning
brain functions that we have mentioned above. One can see Beggs (2007), Chialvo (2010),
cocineros, Perdido, Brilla, and Breakspear (2017), Hesse and Gross (2014), and Shew and Plenz
(2013) for systematic reviews of the diverse function advantages implied by brain criticality
and their experimental demonstrations. From a Darwinian perspective, one potential reason
for the brain to feature criticality lay in that the most informative parts of external world prin-
cipally occur at a borderline between purely ordered and purely disordered states (información
would be trivial in a purely ordered world while it would be incomprehensible in a purely
disordered world). Becoming critical may be a potential way for the brain to adapt to the com-
plex world, where nontrivial information has a finite opportunity to occur (Bak, 2013; Chialvo,
2010). Hasta la fecha, generic features of a critical brain with the characteristics discussed above,
such as divergent correlation length, neuronal avalanches with power law behaviors, y
long-range correlations on the microscopic scale (p.ej., neural populations), have been exten-
sively observed in mathematical models in conjunction with experimental data (p.ej., ruega &
Plenz, 2003; Dalla Porta & Copelli, 2019; Oscuro, Williams-García, ruega, & Ortíz, 2021;
Gireesh & Plenz, 2008; Piedra dura, Mansvelder, & Linkenkaer-Hansen, 2014; Petermann
et al., 2009; Cabello, Piedra dura, Mansvelder, & Linkenkaer-Hansen, 2012; Cabello, van Ooyen, &
Linkenkaer-Hansen, 2008; Ponce-Alvarez, Jouary, Privat, decoración, & Sumbre, 2018; GRAMO. Scott
et al., 2014; Shew et al., 2009; Shriki et al., 2013; Tagliazucchi, Balenzuela, Fraiman, &
Chialvo, 2012; Tkačik et al., 2015).
Our work does not aim at repeatedly reviewing experimental advances concerning brain
criticality and its biological significance, given that they have been comprehensively summa-
rized by existing reviews (ruega, 2007; Chialvo, 2010; Cocchi et al., 2017; Hesse & Bruto,
2014; Muñoz, 2018; mostrar & Plenz, 2013). On the contrary, our motivation is to present a
systematic and accessible review of the theoretical methods applied to achieve these
advances, which have not received necessary attention yet.
These theoretical foundations are initially thought to be incomprehensible and irrelevant to
neurociencia. Sin embargo, practice suggests that omitting these physical and mathematical back-
grounds does not significantly improve the accessibility of studies on brain criticality. En cambio,
the lack of detailed explanations of theoretical foundations has frequently misled neuroscien-
tistas, leading to diverse confusions about the precise meaning, identification criteria, and bio-
logical corollaries of brain criticality. Como resultado, criticidad, an analytic statistical physics theory
with solid foundations, unnecessarily becomes an elusive black box for neuroscientists. A
address this issue, we use the terminology of neuroscience to present a self-contained frame-
work of brain criticality, reviewing and reformulating (1) physical theories that analytically
derive brain criticality and (2) statistic techniques that computationally estimate brain criticality
from empirical data. Given the frequent misunderstanding of neural avalanches, our discus-
sions primarily focus on brain criticality analysis on the microscopic scale of the brain. El
objectives guiding our review are tripartite: (1) explaining why brain criticality matters in the
cerebro, (2) understanding what is brain criticality and what it conveys about the brain, y (3)
confirming how to identify potential brain criticality and ensure the validity of analyses.
Neurociencia en red
1150
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
.
/
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
BRAIN CRITICALITY: GENERAL CONCEPTS
Overview of Brain Criticality
Brain criticality frequently confuses neuroscientists since too many distinct phenomena are
studied under this name without being properly classified. In this review, brain criticality refers
to a family of critical processes in neural dynamics where erratic fluctuations appear to reduce
dynamic stability. To present a systematic classification framework, we discuss three funda-
mental perspectives concerning brain criticality. Mesa 1 provides all the necessary glossaries
in comprehensible forms.
Being nonequilibrium. Primero, the brain, similar to other biological systems, generally exhibits
temporal evolution from initial states that are far away from equilibrium (Gnesotto, Mura,
Gladrow, & Broedersz, 2018; Lynn, Cornblath, Papadopoulos, Bertolero, & bassett, 2021).
These departures from equilibrium arise due to diverse endogenous causes (Gnesotto et al.,
2018; Perl et al., 2021) to break the detailed balance to support consciousness, sintiendo,
and adaptation (Lynn et al., 2021; Perl et al., 2021). Por lo tanto, potential critical phenomena
underlying neural dynamics, at least in most neural dynamics models and empirical datasets,
are basically nonequilibrium and cannot be characterized by equilibrium statistic mechanics.
In Figure 1A, we illustrate the difference between equilibrium and nonequilibrium dynamics.
Fine tuning versus self-organization. Segundo, there exist two types of general mechanisms
underlying the existence of brain criticality. One type of mechanisms either arise from the
external manipulations outside the brain (p.ej., researchers manipulate the tonic dopamine
D1-receptor stimulation; Stewart & Plenz, 2006, 2008) or adjust network topology (Kaiser &
Fine tuning:
A concept in theoretical physics,
which refers to the case where
system parameters must be precisely
manipulated in order to fit with
certain observations.
Concept
Equilibrium
Mesa 1.
Key concepts in describing brain criticality
Significado
A case where the system maximizes entropy and conserves energy simultaneously. The stationary
probability distribution Peq(·) of system states of a system at equilibrium is the Boltzmann distribution. En
equilibrium, the transition dynamics between system states c and c 0 satisfies the detailed balance
condition Peq(C) W.(c → c 0) = Peq(C 0) W.(C 0 → c), where W(· → ·) denotes the transition probability.
Nonequilibrium
A case where the system is out of equilibrium because the transition dynamics between system states breaks
the detailed balance condition. En otras palabras, the transition dynamics between states becomes
directional rather than symmetric.
Self-organization
A process where the internal complexity of a system increases without being tuned by any external
mechanism. All potentially emergent properties are created by endogenous feedback processes or other
internal factors inside the system.
Criticality
A kind of phenomena where the systems is generally close to specific critical points separating between
multiple system states. Small disturbances are sufficient to make the system experience dramatic and sharp
transitions between system states.
Quasi-criticality
A kind of phenomena where all statistical physics relations required by criticality are principally adhered by
the system but slight and inconstant deviations from perfect criticality can be seen on the actual values of
characteristic variables. These deviations robustly exist and are generally independent of data noises.
Sub-criticality
A kind of system states below criticality. They occur when the order parameter (es decir., the macroscopic
observable used to describe system states) remains at zero even with the addition of derives,
corresponding to disordered system dynamics.
Super-criticality
A kind of system states above criticality. They occur when the order parameter is positive, correspondiente a
ordered system dynamics.
Neurociencia en red
1151
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
.
/
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
.
/
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Cifra 1. Conceptual illustrations of brain criticality. (A) Difference between equilibrium and nonequilibrium dynamics in a three-state brain
(paralelo superior). Brain states are characterized by three system components. We illustrate an instance of nonequilibrium dynamics between
these states (paralelo inferior). (B) Fine tuning with exogenous mechanisms (represented by animated hands) makes the brain evolve from a
noncritical state (upper left) to the critical state (upper right). Endogenous mechanisms enable the brain to self-organize from a noncritical state
(abajo a la izquierda) to the critical state (bottom right). (C) Increasing stimulus intensity enlarges the quasi-critical region around the perfect critical
point in a quasi-critical system. (D) The approaching process to a critical point in a self-organized quasi-critical system consists of two stages. En
the first stage, the brain self-organizes from a noncritical state to a quasi-critical region based on certain endogenous mechanisms. En el
second stage, additional exogenous mechanisms are necessary to fine tune the brain to the critical point. De lo contrario, the brain just hovers
within the quasi-critical region. (mi) The difference between four types of brain criticality from the perspective of susceptibility. For standard
brain criticality (p.ej., ordinary criticality and self-organized criticality), susceptibility becomes divergent (es decir., infinite) at the critical point. Para
nonstandard brain criticality (p.ej., quasi-criticality and self-organized quasi-criticality), susceptibility is always nondivergent (es decir., finite). El
quasi-critical region is defined as a set of all control parameters where susceptibility values are no less than a specific threshold (p.ej., half-
maximum value). (F) The commonness and difference between four types of brain criticality.
Hilgetag, 2010; Rubinov, despreciar, Thivierge, & romper la lanza, 2011; S. Wang & zhou, 2012) o
belong to the top-down biological processes that globally function on neural dynamics inside
the brain (p.ej., anesthesia effects; Fontenele et al., 2019; Hahn et al., 2017; Ribeiro et al.,
2010) as well as sleep restoration effects (Meisel, Olbrich, Shriki, & Achermann, 2013). Neural
dynamics is passively fine-tuned toward or away from ordinary criticality (OC) by these
Neurociencia en red
1152
Theoretical foundations of studying criticality in the brain
exogenous mechanisms, similar to ordinary critical phenomena that require the fine tuning of
order parameters.
Another type of mechanisms includes all endogenous factors of neural dynamics (p.ej.,
neural plasticity mechanisms such as spike-timing dependent synaptic plasticity (Effenberger,
Jost, & Levina, 2015; Meisel & Bruto, 2009; espinilla & kim, 2006), short-term synaptic plasticity
(Levina, Herrmann, & Geisel, 2007, 2009), retro-synaptic signals (Hernandez-Urbina &
Herrmann, 2017), and Hebbian rules (De Arcangelis & Herrmann, 2010; De Arcangelis,
Perrone-Capano, & Herrmann, 2006), which locally function on neural dynamics as drive
and dissipation components. The interactions between these components naturally form feed-
back control loops to support the self-organization of neural dynamics toward the critical point
(ruega, 2007; Chialvo, 2010). This spontaneously emerged brain criticality, distinct from ordi-
nary critical phenomena, is conjectured as a kind of self-organized criticality (SOC) (Chialvo,
2010). In Figure 1B, we present conceptual illustrations of ordinary criticality and self-
organized criticality in the brain.
Standard versus nonstandard. Tercero, brain criticality frequently occurs in nonstandard forms
due to stimulus drives or endogenous factors. Por un lado, slight and inconstant devia-
tions from perfect brain criticality can be seen on the actual values of characteristic variables,
differentiating the characterized phenomena from the standard criticality (Oscuro et al., 2021;
Williams-García, moore, ruega, & Ortíz, 2014). Por otro lado, all statistical physics rela-
tions required by perfect brain criticality are still adhered by these actual characteristic vari-
ables, distinguishing the brain from being noncritical (Oscuro et al., 2021; Williams-García
et al., 2014).
For ordinary criticality, its nonstandard form is referred to as quasi-criticality (qC) (Oscuro
et al., 2021; Williams-García et al., 2014). Diverse mechanisms can force the brain to depart
from perfect ordinary criticality and exhibit quasi-critical neural dynamics, among which,
stimulus derive may be the most common one (Oscuro et al., 2021; Williams-García et al.,
2014). En general, sufficiently strong stimulus drives can capture or even govern neural dynam-
circuitos integrados. Similar to the situation where external inputs suppress irregular neural dynamics
(Molgedey, Schuchhardt, & Schuster, 1992), the stimuli that are too strong may evoke intense
but less changeable neural dynamics to make the brain depart from the perfect critical point
(Oscuro et al., 2021; Williams-García et al., 2014). Let us take the qC phenomenon introduced
by Fosque et al. (2021) and Williams-García et al. (2014) as an instance. Under specific con-
ditions, the actual brain state may be close to a Widom line in the three-dimensional space
defined by the stimulus intensity υ, refractory period length τ, and branching ratio κ (es decir., el
time-dependent average number of subsequent neural activities caused by a single neuron
activation event (Haldeman & ruega, 2005). The Widom line is a line of all the combinations
de (υ, t, κ) where the susceptibility of neural dynamics is maximized (Oscuro et al., 2021;
∂y
Williams-García et al., 2014). The susceptibility is defined by limx→0
∂x, where y is the neural
dynamics state and x denotes a factor that affects y. En general, one can understand suscepti-
bility as the degree to which fluctuations in the state of each neuron can propagate to
neighbored neurons (Williams-García et al., 2014). Being close to the Widom line suggests
the existence of quasi-criticality in the brain. Moving along the Widom line as the stimulus
intensity increases, the susceptibility of neural dynamics decreases, and the branching ratio
at maximal susceptibility will decrease as well (Oscuro et al., 2021; Williams-García et al.,
2014). Significant deviations from the Widom line suggest noncriticality (es decir., the subcriticality
where neural dynamics is disordered and the super-criticality where neural dynamics is
ordenado; Williams-García et al., 2014). In Figure 1C, we conceptually illustrate how stimuli
Neurociencia en red
1153
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
t
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
imply qC in the brain. In Figure 2D, the qC phenomenon in Fosque et al. (2021) and Williams-
García et al. (2014) is shown in detail.
As for self-organized criticality (SOC), its nonstandard form is defined according to statisti-
cal physics criteria. Perfect self-organized criticality only exists in conserved neural dynamics
(p.ej., see integrate-and-fire neurons analyzed by Levina et al., 2007), where system energy
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Cifra 2. Conceptual illustrations of the relations between neural avalanches and brain criticality. (A) Instances of neural avalanche, self-
sustained neural dynamics, and vanished neural dynamics. (B) The recovery processes of brain states after the same perturbation in the space
of absorbing and active phases (paralelo superior) and the space of synchronous and asynchronous phases (paralelo inferior). The recovery pro-
cesses after perturbations are relatively fast when the brain is far from the critical point or the quasi-critical region, These recovery processes
slow down when the brain is close to the critical point or the quasi-critical region. (C) The conceptual illustrations of neural dynamics when the
brain state is asynchronous, synchronous, or at ordinary criticality. (D) Without stimuli, there initially exist disordered (gray), ordered but asyn-
chronous (light blue), synchronous (verde) phases in the phase space of the brain. Stimulus inputs imply quasi-criticality in the brain. Un
increasing stimulus intensity enlarges the quasi-critical region (purple) around the Widom line (purple dashed line). (mi) The conceptual illus-
trations of how endogenous mechanisms in conserved neural dynamics can function as drive or dissipation terms to create self-organized
criticality between absorbing and active phases in the brain. (F) In the self-organized quasi-critical brain, endogenous mechanisms in non-
conserved neural dynamics only support the self-organization towards a quasi-critical region between asynchronous and synchronous phases.
Extra exogenous mechanisms are required to fine tune the brain towards the critical point.
Neurociencia en red
1154
Theoretical foundations of studying criticality in the brain
(es decir., neural activities) either conserves within the system and only dissipates at the system
boundary, or dissipates inside the system (es decir., bulk dissipation) with a dissipation rate vanish-
ing in the system size limit (Malcai, Shilo, & Biham, 2006). Under more general conditions
where neural dynamics are not conserved (p.ej., see leaky integrate-and-fire neurons analyzed
by Levina et al., 2007; millman, Mihalas, kirkwood, & Niebur, 2010; Rubinov et al., 2011;
Stepp, Plenz, & Srinivasa, 2015; where neural dynamics dissipates within the system due to
voltage leak), perfect self-organized criticality can be broken by any rate of bulk dissipation
(Bonachela, De Franciscis, Torres, & Muñoz, 2010; Bonachela & Muñoz, 2009; Buendía, di
Santo, Villegas, Burioni, & Muñoz, 2020b; de Andrade Costa, Copelli, & Kinouchi, 2015).
Stronger bulk dissipation implies larger deviations from perfect self-organized criticality (De
Arcangelis et al., 2006). Como consecuencia, the self-organization process of nonconserved neural
dynamics only make the brain hover around the critical point. Any further closeness toward
the critical point requires the fine tuning of order parameter by additional exogenous mecha-
nisms, which is different from pure self-organized criticality (Bonachela et al., 2010;
Bonachela & Muñoz, 2009; Buendía et al., 2020b; de Andrade Costa et al., 2015). This non-
conserved self-organization process is termed as self-organized quasi-criticality (SOqC)
(Bonachela & Muñoz, 2009). Similar to SOC in conserved dynamics, neural plasticity mech-
anisms, such as spike-timing-dependent synaptic plasticity (Rubinov et al., 2011), Hebbian
normas (De Arcangelis et al., 2006), short-term synaptic depression in conjunction with spike-
dependent threshold increase (Girardi-Schappo et al., 2021), and inhibitory plasticity in con-
junction with network topology (Mamá, Turrigiano, wessel, & Hengen, 2019), can serve as
underlying self-organization mechanisms of SOqC. Because purely conserved neural dynam-
ics is relatively rare in empirical data (p.ej., neural dynamics is conserved for integrate-and-fire
neuronas; Levina et al., 2007; and leaky integrate-and-fire neurons whose presynaptic inputs
are exactly equal to the sum of voltage leak and potential costs during neural spiking
(Bonachela et al., 2010), we suggest that SOqC may be more common in the brain than
SOC. In Figure 1D, we present conceptual instances of the two-stage approaching process
towards the critical point in the brain with SOqC.
Classification of brain criticality. The above discussion has presented a classification framework
of brain criticality, eso es, ordinary criticality (OC), quasi-criticality (qC), self-organized criti-
cality (SOC), and self-organized quasi-criticality (SOqC). In Figure 1E, we compare between
these four types of brain criticality in term of susceptibility. En general, susceptibility diverges at
the critical point in a brain with standard criticality (p.ej., OC and SOC) while it does not
diverge in the quasi-critical region of a brain with nonstandard criticality (p.ej., qC and SOqC).
In Figure 1F, we summarize the commonness and difference between these four types of brain
criticality discussed in our review. From a neuroscience perspective, a brain with critical neu-
ral dynamics is expected to be near the critical point and prepared for tremendous changes in
cortical states during a short duration. This intriguing property coincides with the experimen-
tally observed efficient transitions between cortical states (p.ej., Cardin, 2019; Holcman &
Tsodyks, 2006; Jercog et al., 2017; h. Sotavento, Wang, & Hudetz, 2020; Reimer et al., 2014)
y, por lo tanto, interests researchers for the potential existence of brain criticality. The impor-
tance of identifying brain criticality in neural dynamics is beyond brain criticality itself because
it implies an opportunity to explain and predict brain function characteristics by various
statistical physics theories built on nonequilibrium criticality.
Neural Avalanches and Their Phases
To identify potential nonequilibrium criticality in the brain, researchers actually characterize
neural dynamics as a physical system with absorbing states and avalanche behaviors
Neurociencia en red
1155
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
t
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
(Hinrichsen, 2000; Larremore, Carpintero, Ott, & Restrepo, 2012; Lübeck, 2004). En general,
one needs to consider the propagation of neural dynamics where neurons are either activated
(“on” state) or silent (”off” state) (Dalla Porta & Copelli, 2019). A silent neuron may be acti-
vated with a probability defined by the number of activated presynaptic neurons and the cou-
pling strength θ among neurons (p.ej., neural correlation; Franke et al., 2016). An activated
neuron spontaneously becomes silent at a constant rate (p.ej., after the refractory period;
Kinouchi & Copelli, 2006; Squire et al., 2012). These definitions naturally support to distin-
guish between different phases of neural dynamics. Here we review two kinds of phase
partition that are active in neuroscience.
Absorbing versus active. The first group of phases are absorbing and active phases (Larremore
et al., 2012). The absorbing phase refers to cases where couplings between neurons are weak
and all neurons eventually become silent (neural dynamics vanishes). Once a neural dynamics
process vanishes, it cannot reappear by itself. The brain requires new drives (p.ej., neuronas
activated spontaneously or by stimuli) to trigger new neural dynamics. The active phase, en
the other hand, correspond to cases where the “on” state propagates among neurons with
strong couplings, leading to stable self-sustained neural dynamics (p.ej., nonzero time- y
ensemble-averaged density of active neurons in the brain). In Figure 2A, we show conceptual
instances of neural avalanches, self-sustained neural dynamics, and vanished neural dynam-
circuitos integrados. Denoting ρ(t) as the density of active neurons at moment t, we can simply represent the
absorbing (Ecuación 1) and active (Ecuación 2) phases of a neural dynamics process triggered
by an active neuron at moment 0 como
ρ tð Þ ¼ 0; ∃t > 0;
ρ tð Þ > 0; ∀t > 0:
(1)
(2)
Synchronous versus asynchronous. The second group of phases are synchronous and asynchro-
nous phases (di Santo, Villegas, Burioni, & Muñoz, 2018; Fontenele et al., 2019; Girardi-
Schappo et al., 2021). As their names suggest, these two phases correspond to the situations
where synchronization emerges or disappears in neural activities, respectivamente. Synchroniza-
tion refers to the cases where “on” states appear in an oscillatory, although not strictly peri-
odic, manner. To quantify its potential existence, we can measure the variability of neural
dynamics using the coefficient of variation (CV ) (di Santo et al., 2018; Fontenele et al.,
2019; Girardi-Schappo et al., 2021) or the Kuramoto order parameter (Acebrón, Bonilla, Pérez
Vicente, Ritort, & Spigler, 2005; Arenas, Díaz-Guilera, Kurths, Moreno, & zhou, 2008). CV
can be defined from diverse perspectives, yet the most common definition is the ratio between
the standard deviation and the mean of the interspike interval length (di Santo et al., 2018;
Fontenele et al., 2019; Girardi-Schappo et al., 2021). A higher value of CV implies the
reduction of synchronization. For most neural dynamics data, an empirical choice of the
CV threshold that separates between synchronous and asynchronous phases may be ≃1
(Fontenele et al., 2019) or ≃ 3
2 (Fontenele et al., 2019). The Kuramoto order parameter
Vaya 2 [0, 1] measures the coherent degree of neural dynamics based on the Kuramoto model
of oscillators (for detailed definitions see Acebrón et al., 2005; Arenas et al., 2008). Perfect
synchronization emerges when ω = 1 and vanishes when ω = 0 (Acebrón et al., 2005;
Arenas et al., 2008).
Critical point or quasi-critical region. The boundary between these two phases is the critical
punto, at which the brain is on the edge of exhibiting self-sustained (for absorbing and active
phases) or synchronous (for synchronous and asynchronous phases) neural dynamics.
Neurociencia en red
1156
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Perturbations (p.ej., the propagation of “on” state among neurons) to the absorbing or asynchro-
nous phase do not have characteristic lifetime and size. These perturbations, referred to as
neural avalanches, are expected to exhibit power law properties in their lifetime (time differ-
ence between the first and last activation of neurons in between complete quiescent epochs)
and size (number of active neurons along with the excursion) distributions (Hesse & Bruto,
2014; Hinrichsen, 2000; Larremore et al., 2012; Lübeck, 2004). En general, the emergence
of neural avalanches implies the slowing down of neural dynamics, eso es, the brain state
recovery process toward the baseline state after fluctuations changes from fast (exponential)
to slow (ley de potencia) (Cocchi et al., 2017; Hesse & Bruto, 2014). The dynamic stability of neu-
ral dynamics is limited by the slow recovery and, por lo tanto, cannot robustly counteract per-
turbations. Como consecuencia, small perturbations initiated on the microscopic scale may still make
the brain change sharply on the macroscopic scale (Cocchi et al., 2017; Hesse & Bruto, 2014).
In Figure 2B, we conceptually illustrate how the recovery process slows down when the brain
is close to the critical point or the quasi-critical region.
General Relations Between Neural Avalanches and Brain Criticality
The relation between neural avalanches and brain criticality is frequently neglected or misun-
entendido. Neural avalanche data alone is not sufficient to determine the concrete type of brain
criticidad (es decir., OC, qC, SOC, and SOqC) unless additional information about the mechanisms
underlying neural avalanche emergence is provided (p.ej., if neural dynamics is conserved or
self-organizing). To explore a concrete type of brain criticality, researchers need to explicitly
present its definition depending on different control parameters (p.ej., the balance between
excitatory and inhibitory neurons in CROS models; Hardstone et al., 2014; Poil et al.,
2012) and order parameters (p.ej., active neuron density and synchronous degree; Dalla Porta
& Copelli, 2019). A brain criticality hypothesis without strict definitions of control and order
parameters is not informative (Cocchi et al., 2017; Girardi-Schappo, 2021). To present concep-
tual instances, we illustrate four possible critical phenomena in Figure 2, each of which cor-
responds to a concrete brain criticality type.
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
.
/
Instance of ordinary criticality. To produce ordinary criticality (OC), we can control neural
dynamics and manipulate hθi, the expectation of coupling strength θ among all neurons
(p.ej., averaged neural correlation), by some top-down and global biological effects. Estos
efectos, por ejemplo, may be anesthesia effects (p.ej., by ketamine-xylazine; Ribeiro et al.,
2010; and isoflurane; Hahn et al., 2017) or sleep restoration effects; Meisel et al., 2013).
We use the Kuramoto order parameter ω (Acebrón et al., 2005; Arenas et al., 2008) as the order
parameter to define synchronous and asynchronous phases (di Santo et al., 2018; Fontenele
et al., 2019). As hθi increases, we may see transitions from asynchronous to synchronous
phase in some situations (see a similar instance in Villegas, Moretti, & Muñoz, 2014). Uno
can see Figure 2C for conceptual illustrations.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Instance of quasi-criticality. To produce quasi-criticality (qC), we can manipulate refractory
period length τ, branching ratio κ, and stimulus intensity υ as control parameters (p.ej., control
τ and κ by pharmacological perfusion or ionic concentration adjustment; Chiappalone et al.,
2003; Shew et al., 2011). There exist a disordered phase (subcritical), an ordered but asynchro-
nous phase (supercritical), and a synchronous (quasi-periodic) phase in the space of (υ, t, κ)
(Oscuro et al., 2021; Williams-García et al., 2014). These phases can be characterized by spe-
cific order parameters related to synchronization. As υ increases, a qC phenomenon emerges
in the space, where the quasi-critical region is defined by all combinations of (υ, t, κ) cuyo
susceptibility values are at least half-maximum. Cross-over behaviors (es decir., a generalization of
Neurociencia en red
1157
Theoretical foundations of studying criticality in the brain
phase transition with finite susceptibility) emerge when the quasi-critical region has overlaps
with at least two phases (Oscuro et al., 2021; Williams-García et al., 2014). In Figure 2D, nosotros
show this qC phenomenon in details.
Instance of self-organized criticality. To study self-organized criticality (SOC), consideramos el
conserved neural dynamics generated by integrate-and-fire neurons (Levina et al., 2007). El
order parameter is active neuron density ρ, whose dynamics is controlled by parameter hθiA,
the averaged coupling strength θ between activated neurons and their postsynaptic neurons
(here A denotes the set of activated neurons). In specific cases, the considered neural dynamics
may self-organize to the critical point under the joint effects of excitatory and inhibitory neu-
ron, neural spiking processes (activation and silence), as well as neural plasticity. In Figure 2E,
we conceptually illustrate a case where these endogenous mechanisms enable the brain to
self-organize to the criticality between absorbing and active phases.
Instance of self-organized quasi-criticality. To analyze self-organized quasi-criticality (SOqC), nosotros
consider the nonconserved neural dynamics affected by two homeostatic adaptation processes,
eso es, the short-term depression of inhibition and the spike-dependent threshold increase.
These processes are controlled by y^, the maximum inhibitory coupling strength, as well as τx
and τy, the decay timescales of neural activation threshold increase and synaptic depression.
These control parameters affect neural activation threshold x and inhibition strength y to shape
neural dynamics states (p.ej., the active neuron density ρ). With appropriate x, y, and ρ, neural
avalanches with power law behaviors will occur to indicate the criticality between an asyn-
chronous phase (stochastic oscillations) and a synchronous phase (periodic oscillations).
According to Girardi-Schappo et al. (2021), x and ρ self-organize to their appropriate values
through quasi-critical fluctuations under biologically reasonable conditions (es decir., τx ≫ 1) mientras
y hovers around the expected value. Additional fine tuning of y based on exogenous mecha-
nisms are necessary to place neural dynamics at the perfect criticality. Mientras tanto, sináptico
homeostasis is discovered as constantly canceled by the variation of the activation threshold,
impeding neural dynamics from self-organizing to perfect criticality. In Figure 2F, we concep-
tually illustrate the defined SOqC phenomenon in a similar manner of Figure 2D and Figure 2E.
As for the precise description of quasi-critical fluctuations, one can see Girardi-Schappo et al.
(2021) for details.
To this point, we have conceptually introduced the phenomenological properties of brain
criticidad. To verify the hypothetical brain criticality, one needs to learn about analytic brain
criticality theories and the properties of neural avalanche predicted by them. Abajo, we pres-
ent accessible expositions of these theoretical foundations.
BRAIN CRITICALITY: PHYSICAL THEORIES
Mean-Field and Stochastic Field Theories of Brain Criticality
One of the main challenges faced by neuroscientists in studying ordinary criticality (OC),
quasi-criticality (qC), self-organized criticality (SOC), and self-organized quasi-criticality
(SOqC) is how to understand their theoretical relations (Girardi-Schappo, 2021). Overcoming
this challenge is crucial for understanding why we can verify the existence of different types of
brain criticality with certain theoretical tools. To present a concise and thorough review, nosotros
first focus on brain criticality between absorbing and active phases, where we generalize the
ideas in Bonachela and Muñoz (2009) and Buendía, di Santo, Bonachela, and Muñoz (2020a)
to present a possible framework for unification.
Neurociencia en red
1158
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
.
t
/
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
En general, brain criticality in the space of absorbing
Langevin formulation of ordinary criticality.
and active phases are related to directed percolation (Dalla Porta & Copelli, 2019), a univer-
sality class of continuous phase transitions into absorbing states (Hinrichsen, 2000; Lübeck,
2004). Aquí, a universality class can be understood as the set of all systems with the same
scaling properties (Hinrichsen, 2000; Lübeck, 2004; Sethna, Dahmen, & miers, 2001).
Directed percolation theory initially covers OC phenomena (Hinrichsen, 2000; Lübeck,
2004). Let us begin with a variant of the classic Reggeon field theory, the simplest description
of absorbing phase transitions (Henkel, Hinrichsen, Lübeck, & Pleimling, 2008). The Lange-
→
vin equation of the activity neuron field ρ(X
, t ) is defined as
(cid:3)
(cid:1)
ρ x
→; t
(cid:1)
(cid:1)
→; t
¼ a þ bν x
(cid:3)
(cid:3)
(cid:3)
(cid:1)
→; t
ρ x
(cid:1)
→; t
− cρ2 x
(cid:3)
(cid:1)
þ d∇2ρ x
→; t
(cid:3)
þ e
r
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
(cid:1)
(cid:3)
(cid:1)
→; t
→; t
σ x
ρ x
(cid:3)
;
(cid:3)
(cid:1)
ν x
→; t
(cid:1)
→; t
¼ ∇2ν x
(cid:3)
(cid:1)
→; t
þ f x
(cid:3)
(cid:1)
→; t
− g x
(cid:1)
(cid:3)
→; t
ρ x
(cid:3)
;
(3)
(4)
∂
∂t
∂
∂t
, t ) ρ(X
→ − x
, t ) = 0 under the effects of factor
→
represents spatial coordinates, a 2 ℝ, b 2 (0, ∞), C 2 (0, ∞), d 2 ℝ is the diffusion
where x
factor, and e 2 ℝ is the noise factor. Function σ(·, ·) defines a zero-mean Gaussian noise with
→ 0) δ(t − t 0), where δ(·) is the delta
→ 0, t 0)i = δ(X
→
a spatiotemporal correlation h ρ(X
función. En general, pag(·, ·) reflects the collective fluctuations in neural activities that vanish
ffiffiffiffiffiffiffiffiffiffiffiffiffiffi
→
→
ρ x→; t
in the absorbing phase ρ(X
. The term ∇2ρ(X
, t )
Þ
d
→
reflects the propagation of neural dynamics. The function ν(X
, t) defines the energy (es decir.,
→
membrane potential) that propagates according to ∇2ν(X
, t ), increases with external drives
→
, t ) ≥ 0
F (X
always hold. The initial active neuron density and energy are assumed as nonzero. It is clear
→
that a + bν(X
, t ) < 0 makes the neural dynamics eventually vanish (i.e., absorbing phase)
→
while a + bν(x
, t ) > 0 does not (es decir., active phase). Por lo tanto, we can fine tune the
→
control parameter ν(X
, t ) = 0, a
→
critical point defined by νc. The fine tuning relies on manipulating f (X
, t ) por
exogenous mechanisms.
→
, t ) to make the brain exhibit OC dynamics at a + bνc(X
→
, t ), and decreases with bulk dissipation g (X
→
, t ). Please note that ρ(X
→
, t ) ≥ 0 and ν(X
→
, t ) and g (X
pag
Langevin formulation of quasi-criticality. Then we turn to analyzing qC, whose mean-field
approximation is initially derived based on the cortical branching model (Oscuro et al.,
2021; Williams-García et al., 2014). A cortical branching model with no stimulus input
belongs to the directed percolation universality class according to the Janssen-Grassberger
conjecture (Williams-García et al., 2014). Nonzero stimulus inputs make the cortical
branching model depart from directed percolation universality class to create qC (williams-
García et al., 2014). Sin embargo, the above mean-field theory is defined in the space of
synchronous and asynchronous phases. To derive a qC phenomenon between absorbing
and active phases, we can provisionally analyze a mean-field approximation of Equations 3
y 4:
(cid:3)
(cid:1)
→; t
ρ x
(cid:1)
(cid:1)
→; t
¼ a þ bν x
(cid:3)
(cid:1)
(cid:3)
→; t
ρ x
(cid:3)
(cid:1)
→; t
− cρ2 x
(cid:3)
;
(cid:3)
(cid:1)
→; t
ν x
(cid:1)
¼ f x
→; t
(cid:3)
(cid:1)
→; t
− g x
(cid:1)
(cid:3)
→; t
ρ x
(cid:3)
;
∂
∂t
∂
∂t
(5)
(6)
→
where ∇2ρ(X
→
field assumption. We consider the cases where stimulus inputs vanish, eso es, F (X
, t ) in Equations 3 y 4 are neglected under the mean-
, t ) ≡ 0. El
→
, t ), and σ(X
→
, t ), ∇2ν(X
Neurociencia en red
1159
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
/
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
critical point between active and absorbing phase becomes νc = − a
of Equations 5 y 6 son
b. The steady state solutions
(cid:3)
→; t
(cid:1)
ρ x
(cid:1)
→; t
ν x
¼ 0;
(cid:3)
¼ r 2 0; ∞ð
Þ;
(7)
(8)
respectivamente. Por lo tanto, OC is one of the steady states of neural dynamics when there is no stim-
ulus. In the cases where stimulus inputs become increasingly strong, there exists no steady-state
solution of Equations 5 y 6 unless f x→;t
Þ
→ r 2 (0, ∞) sostiene, we can derive
→;t
Þ
→; t
f x
d
g x→; t
d
(cid:5)
→ r 2 (0, ∞). If f x→;t
→;t
(cid:1)
ρ x
→ r;
→; t
d
g x
d
d
g x
d
(9)
Þ
Þ
¼
(cid:6)
(cid:3)
Þ
Þ
(cid:3)
(cid:1)
ν x
→; t
¼
1
b
→; t
f x
d
g x→; t
d
Þ
Þ
C
1
b
− a
→
cr − a
d
Þ:
(10)
Because the critical point νc = −a
b is not necessarily a steady state, it can be disturbed by
diverse factors (p.ej., by stimuli). Unless there exist certain ideal exogenous mechanisms that
→
persistently enlarge g (X
, t ) aumenta, the fine tuning of neural dynamics
→
cannot cancel the effects of f (X
, t ). Como consecuencia, the fine tuning process may only enable
the brain to reach a quasi-critical region where the susceptibility of neural dynamics is rela-
tively large. The initial OC vanishes and is replaced by qC.
→
, t ) whenever f (X
Langevin formulation of self-organized criticality. Although SOC is treated as a rather isolated con-
cept after its first discovery in statistical physics (Bak, Espiga, & Wiesenfeld, 1987), subsequent anal-
yses demonstrate SOC as relevant with ordinary continuous phase transitions into infinitely many
absorbing states (Dickman, Muñoz, Vespignani, & Zapperi, 2000; Dickman, Vespignani, &
Zapperi, 1998; Narayan & Middleton, 1994; Sornette, Johansen, & Dornic, 1995). Específicamente,
SOC models can be subdivided into two families, which we refer to as external dynamics family
(p.ej., Bak-Sneppen model; Bak & Sneppen, 1993) and conserved field family (p.ej., sandpile models
such as Manna model; Manna, 1991; and Bak-Tang-Wiesenfeld model; Bak et al., 1987). El
second family, being the main theoretical source of studying SOC in neural dynamics, corresponds
to absorbing-state transitions since it can represent any system with conserved local dynamics and
continuous transitions to absorbing states (Dickman et al., 2000; Lübeck, 2004). Aunque el
universality class of the second family should be precisely referred to as conserved directed per-
colation, the explicit behaviors (p.ej., avalanche exponents and scaling relations) of conserved
directed percolation are similar to those of directed percolation in high-dimensional systems
(p.ej., neural dynamics) (Bonachela & Muñoz, 2008, 2009; Buendía et al., 2020a). Por lo tanto,
SOC and OC share some identification criteria in practice. To understand the connections between
→ 0
→
SOC and OC more precisely, we can consider the cases where f(X
, t) → 0 such that f x→;t
→;t
d
g x
d
Þ
Þ
(es decir., infinite separation of timescales). The steady-state solutions of Equations 5 y 6 convertirse
→; t
f x
d
g x→; t
d
(cid:5)
(cid:1)
ρ x
→; t
→ 0;
Þ
Þ
¼
(cid:6)
(cid:3)
(11)
(cid:3)
(cid:1)
→; t
ν x
¼
1
b
→; t
f x
d
g x→; t
d
Þ
Þ
C
− a
→ νc;
(12)
respectivamente. Self-organization properties are reflected by the following processes: if the brain is
→
in the absorbing phase because neural dynamics vanishes, eso es, ρ(X
, t ) → 0, Ecuación 6
→
becomes ∂
norte (X
, t ) to shift the brain toward the active phase; if the brain is in the active
∂t
phase, Ecuación 6 becomes ∂
, t ) ≪
∂t
→
, t ) to reduce neural dynamics since f (X
→
, t ) ≃ −g (X
→
, t ) = f (X
→
, t ) ρ(X
→
norte (X
Neurociencia en red
1160
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
/
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
→
, t ). These feedback control loops drive the brain to the critical point. One may be curious
gramo (X
→
, t ) → 0, is necessary for SOC since the above der-
about why energy conservation, eso es, gramo (X
→
, t ) → 0 en
ivations seem to be independent of g (X
Ecuación 14 makes the active phase no longer exist. En otras palabras, the nonconserved energy
implies a kind of continuous phase transition that does not belong to conserved directed
percolation or directed percolation when the infinite separation of timescales is satisfied.
Por lo tanto, energy conservation is necessary for SOC.
→
, t ) → 0. Later we show that the absence of g(X
Langevin formulation of self-organized quasi-criticality. As for SOqC, nonzero bulk dissipation
breaks the conservation law to generate non-Markovian components in neural dynamics
(Bonachela & Muñoz, 2009; Buendía et al., 2020a). In the ideal cases where the drive terms
(p.ej., stimulus inputs) of a sufficiently large neural dynamics system occur at an arbitrarily slow
timescale (es decir., only occur in the interval between neural avalanches), the brain exhibits pure
dynamical percolation behaviors (Buendía et al., 2020a). To understand this property, let us
→
consider a variant of Equations 5 y 6 where the dissipation term g (X
(cid:3)
(cid:1)
;
, t ) is nonnegligible
(cid:1)
→; t
− cρ2 x
(13)
(cid:3)
(cid:3)
(cid:3)
(cid:1)
→; t
ρ x
(cid:1)
→; t
ν x
(cid:3)
(cid:1)
→; t
¼ a þ bν x
(cid:1)
(cid:1)
→; t
ρ x
¼ −g x
→; t
(cid:3)
(cid:1)
(cid:3)
→; t
ρ x
(cid:3)
:
∂
∂t
∂
∂t
(cid:3)
(cid:1)
(cid:3)
→; t
ρ x
R
→
t
0 gramo(X
, t)
By integrating Equation 14 and plugging the integral into Equation 13, we can derive
(cid:3)
(cid:1)
→; t
ρ x
(cid:1)
(cid:1)
→; 0
¼ a þ bν x
∂
∂t
(cid:3)
(cid:1)
→; t
− cρ2 x
(cid:3)
(cid:1)
→; t
− bρ x
(cid:3) z
t
(cid:1)
→; t
g x
(cid:1)
(cid:3)
→; t
ρ x
(cid:3)
0
dτ;
(15)
→
, t) ρ(X
→
The non-Markovian term −bρ(X
, t)dτ in Equation 15 makes the regions already
visited by neural dynamics become more unlikely to be activated (Bonachela & Muñoz, 2009;
Buendía et al., 2020a). Por lo tanto, the pure self-sustained active phase vanishes and is replaced
by a spreading phase, where local perturbations can transiently propagate across the whole system
without reaching a self-sustained state, and a nonspreading phase, where local perturbations can
never span the entire system (Bonachela & Muñoz, 2009; Buendía et al., 2020a). The phase
transition and corresponding critical point νd > νc between spreading and nonspreading phases
belong to the universality class of dynamical percolation rather than conserved directed per-
colation (Bonachela & Muñoz, 2009; Buendía et al., 2020a). The initial neural dynamics can be
created by random shifts at moment 0 (Bonachela & Muñoz, 2009; Buendía et al., 2020a)
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
t
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
(14)
(16)
(17)
(cid:7)
→
ρ x
(cid:8)
* ; 0
→ (cid:2);
(cid:7)
→
ν x
(cid:8)
*; 0
(cid:7)
→
→ ν x
(cid:8)
*; 0
(cid:7)
→
þ h x
*; 0
(cid:8)
;
→
* is a randomly selected coordinate, and function h(·, ·) is a driving function of energy at
where x
moment 0. Every time a neural avalanche occurs after random shifts, the strong dissipation term
→
gramo(X
, t ) pushes the brain towards the subcritical phase. Como consecuencia, the brain cannot exactly
self-organize to the perfect criticality. En cambio, the brain just hovers around the critical point νd to
form a quasi-critical region, exhibiting finite fluctuations to the both sides of νd . In the more
realistic cases where the drive terms do not necessarily occur at an arbitrarily slow timescale
(es decir., can occur at an arbitrary moment), sin embargo, neural dynamics may be phenomenology
controlled by conserved directed percolation transitions and hover around the critical point.
Let us add a drive term in Equation 14
(cid:3)
(cid:3)
(cid:1)
→; t
ν x
(cid:1)
→; t
¼ f x
(cid:1)
→; t
− g x
(cid:1)
(cid:3)
→; t
ρ x
(cid:3)
:
∂
∂t
(18)
1161
Neurociencia en red
Theoretical foundations of studying criticality in the brain
Then Equation 15 becomes
(cid:1)
(cid:1)
→; t
→; t
ρ x
ρ x
(cid:3)
(cid:1)
→; 0
¼ a þ bν x
(cid:3)
(cid:1)
(cid:3)
∂
∂t
(cid:3)
(cid:1)
→; t
− cρ2 x
(cid:3)
(cid:1)
→; t
− bρ x
(cid:3)z
0
(cid:3)
t
(cid:1)
(cid:1)
→; t
f x
(cid:1)
→; t
− g x
(cid:1)
(cid:3)
ρ x
→; t
(cid:3)
(cid:3)
dτ:
(19)
→
If we can ideally fine tune the drive term f (X
, t ) to ensure that
→;t
d
f x
→;t
g x
d
Þ
Þ
→ r 2 (0, ∞), the steady-
state solutions of Equations 18 y 19 son
(cid:3)
(cid:1)
→; t
ρ x
¼
(cid:3)
(cid:1)
→; 0
ν x
¼
→ r;
Þ
Þ
→; t
f x
d
→; t
g x
d
(cid:5)
→; t
f x
d
→; t
g x
d
(cid:6)
− a
→
1
b
Þ
Þ
cr − a
Þ;
d
C
1
b
(20)
(21)
(cid:3)
(cid:1)
ν x
→; t
(cid:3)
(cid:1)
→; 0
¼ ν x
þ
z
0
t
(cid:3)
(cid:1)
(cid:1)
→; t
f x
(cid:1)
→; t
− g x
(cid:1)
(cid:3)
→; t
ρ x
(cid:3)
(cid:3)
dτ →
1
b
cr − a
Þ:
d
(22)
→
Equations 20–22 correspond to a steady state of the brain with ρ(X
, t ) → r and conserved
energía, which is similar to SOC. Por lo tanto, the brain may self-organize to a quasi-critical
region around νc, the critical point of SOC. Reaching the critical point requires ideal fine
tuning. These emerged conserved directed percolation behaviors enable scientists to recog-
nize SOqC in a similar manner of SOC in practice (es decir., when stimulus inputs can occur at
any moment) (Bonachela & Muñoz, 2009; Buendía et al., 2020a).
Summary of theoretical relations. Tomados juntos, neuroscientists can approximately verify the
existence of brain criticality in the space of absorbing and active phases with specific tools
coming from directed percolation theory. This is because OC, qC, SOC, and SOqC exhibit
or approximately exhibit directed percolation behaviors under certain conditions. The verifi-
cation may be inaccurate since the approximation holds conditionally. As for the brain
criticality between asynchronous and synchronous phases, sin embargo, the universality class
properties become rather elusive because an analytic and complete theory of synchronous
phase transitions in the brain remains absent yet (see Buendía, Villegas, Burioni, & Muñoz,
2021; di Santo et al., 2018; for early attempts). Although some behaviors of absorbing phase
transitions can be observed in synchronous phase transitions (p.ej., see Buendía et al., 2021; di
Santo et al., 2018; Fontenele et al., 2019; Girardi-Schappo et al., 2021), there also exist
numerous differences between them (p.ej., see Buendía et al., 2021; Fontenele et al., 2019;
Girardi-Schappo et al., 2021). As suggested by Dalla Porta and Copelli (2019), it remains
elusive if directed percolation properties are applicable, at least conditionally applicable, a
analyzing synchronous phase transitions. More explorations are necessary in the future.
There are numerous properties of brain criticality predicted by directed percolation theory,
among which, neural avalanche exponents (the power law exponents of lifetime and size dis-
tributions), scaling relation, universal collapse shape, and slow decay of autocorrelation are
applicable in both analytic derivations and statistical estimations from empirical data. Estos
properties are our main focus. Por conveniencia, we summarize important glossaries and sym-
bol conventions before we discuss theoretical details (Mesa 2).
Neural Avalanche Exponents
As we have mentioned above, neural avalanches are expected to exhibit power law properties
in their lifetime and size distributions when the brain is at the critical point (Hinrichsen, 2000;
Larremore et al., 2012; Lübeck, 2004). Por lo tanto, it is pivotal to confirm the detailed values of
Neurociencia en red
1162
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
.
/
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Variable
t
S
A
hS(t )i
hS(t |t )i
PT (t )
PS (s)
a
b
γ
h(·)
Aquellos(·, ·)
χ
ξ
Mesa 2. Glossaries and symbol conventions
Significado
The lifetime of the neural avalanche
The size of the neural avalanche
The area of the neural avalanche
The averaged size of neural avalanches with lifetime T
The averaged time-dependent avalanche size at moment t during neural avalanches with the lifetime T
The probability distribution of neural avalanche lifetime
The probability distribution of neural avalanche size
Power law exponent of the neural avalanche lifetime distribution PT (t ) ∝ t−α
Power law exponent of the neural avalanche size distribution PS (s) ∝ s−β
Power exponent of the neural avalanche area A ∝ T γ
Universal scaling function
Autocorrelation function
Power law decay rate of autocorrelation
Exponential decay rate of autocorrelation
Nota. Please note that Table 2 mainly contains important glossaries with fixed symbol definitions. There are many symbols uncovered by Table 2 since they are
only used for mathematical derivations.
neural avalanche exponents. To analytically derive these exponents, one can consider critical
branching process (di Santo, Villegas, Burioni, & Muñoz, 2017; García-Pelayo, Salazar, &
Schieve, 1993; Gros, 2010; harris, 1963; Otter, 1949), neural field theory (robinson,
2021), mean-field Abelian sandpile models (Janowsky & Laberge, 1993; D. S. Sotavento, Goh,
Kahng, & kim, 2004), and avalanches in networks (Larremore et al., 2012). The key idea to
derive neural avalanche exponents shared by these existing theories is to confirm the explicit
forms of PT (t ) and PS (s), the probability distributions of the lifetime and size of neural ava-
lanches, under ideal conditions (p.ej., when the maximum lifetime and size are unlimited
and can be infinitely large). In real cases where lifetime and size are restricted because the
brain is a finite system, slight deviations from idea values may be observed but theoretical
derivations of neural avalanche exponents principally hold.
To present accessible expositions, we consider a critical branching process in Equations
23–34 to describe related backgrounds. More importantly, we present a novel and simple idea
to calculate target exponents in the context of neuroscience in Box 1. Abstractly, one can
define P(norte, t ) as the probability for an active neuron at moment t to activate n postsynaptic
neurons subsequently and define Z(norte, t ) as the probability of finding n active neurons at
moment t. Mientras tanto, one denotes
F x; t
d
Þ ¼
G x; t
d
Þ ¼
X∞
n¼0
X∞
n¼0
P n; t
d
Þx n;
Z n; t
d
Þx n
(23)
(24)
1163
Neurociencia en red
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
/
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
as the corresponding generating functions (Fristedt & Gray, 2013; Rao & Swift, 2006). Entonces,
one can readily see the recursion relation
G x; t
d
Þ ¼
X∞
Z n; t − δt
d
ÞF x; t − δt
d
Þn;
n¼0
¼ G F x; t − δt
d
d
Þ; t − δt
Þ;
(25)
(26)
where δt denotes the minimum time step. Ecuación 26 implies that branching processes are
Markovian. Similarmente, one can measure the expectations
μ tð Þ ¼
ϕ tð Þ ¼
∂
∂x
∂
∂x
F x; t
d
Þjx¼1;
G x; t
d
Þjx¼1
to derive another recursion relation
ϕ tð Þ ¼
∂
∂x
F x; t − 1
d
Þjx¼1
∂
∂x
G x; t − 1
d
Þjx¼1;
¼ ϕ t − δ
d
Þμ t − δ
d
Þ;
¼
Yt−δ
τ¼0
μ τð Þ:
(27)
(28)
(29)
(30)
(31)
Note that Equation 31 is derived from the fact that ϕ(0) = 1 (one neuron is activated at moment
0 to trigger neural avalanches). Please see Marković and Gros (2014) for more explanations of
Equations 25–31. Assuming that ϕ(t ) scales as exp( λt ) for large t, we know that ϕ(t ) converges
a 0 given a negative Lyapunov exponent λ (the branching process is subcritical; di Santo et al.,
2017; García-Pelayo et al., 1993; Gros, 2010; harris, 1963; Otter, 1949) and diverges with a
positive Lyapunov exponent λ (the branching process is supercritical; di Santo et al., 2017;
García-Pelayo et al., 1993; Gros, 2010; harris, 1963; Otter, 1949). Here λ can be defined
according to Equation 32
(cid:5)
λ ¼ lim
t→∞
ln
(cid:6)
ϕ tð Þ
¼ lim
t→∞
1
t
1
t
Xt−δ
τ¼0
ln μ τð Þ
d
Þ:
(32)
If the branching process is homogeneous, namely P(norte, t) =P(norte), z(norte, t ) = Z(norte), μ(t) = μ, and ϕ(t) =
ϕ for every moment τ, then μ = 1 is the condition for the branching process to be critical. A
relate these results with neural avalanches, one only need to consider the avalanche size S =
(cid:2)t z(t ), where z(t ) ∼ Z denotes the number of active neurons at moment t, and the avalanche
life time T = min{t | z(t ) > 0 and z(t + δt ) = 0}. It has been analytically proved that in terms of
fixed environments and a Poisson generating function F one can derive (Otter, 1949)
PS sð Þ ∼ s−3=2μs−1 exp s 1 − μ
d
½
(cid:2);
Þ
PT tð Þ ∼ t −2μt−1 exp t 1 − μ
d
½
(cid:2):
Þ
(33)
(34)
In the case with μ = 1, one can obtain PS (s ) ∼ s−3/2 and PT (t ) ∼ t−2, the power law
distributions of neural avalanche size and neural avalanche life time (di Santo et al., 2017;
García-Pelayo et al., 1993; Gros, 2010; harris, 1963; Janowsky & Laberge, 1993; Jung, Le,
Sotavento, & Sotavento, 2020; Larremore et al., 2012; D. S. Lee et al., 2004; Lombardi, Herrmann, &
De Arcangelis, 2017; Otter, 1949; robinson, 2021), from Equations 33 y 34.
Neurociencia en red
1164
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
.
/
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
The derivations of avalanche exponents α = 2 and β = 3
2 are nontrivial. Sin embargo, pocos
neuroscience studies elaborate on these details, impeding researchers from understanding
the theoretical foundations of brain criticality in the brain. The importance of these deriva-
tions is beyond the detailed values of avalanche exponents since they reveal the fundamental
properties of neural dynamics (Cocchi et al., 2017; di Santo et al., 2017; Girardi-Schappo,
2021). In Box 1, we sketch an original idea to derive these avalanche exponents in the
terminology of neuroscience. In Figure 3A, we present graphical illustrations of our idea in
Box 1.
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
/
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Box 1. Derivations of neural avalanche exponents.
Consider a time-continuous neural dynamics process, where an active neuron implies three
possibilities: becoming absorbed with probability ς, activating another neuron with probabil-
ity η, or remaining effect-free with probability 1 − (ς + η). In critical states, we have ς = η
(García-Pelayo et al., 1993). We define An(t ) as the probability for n active neurons to exist
at t * + t given that 1 active neuron exists at t *. Assuming the independence of neuron acti-
vación, tenemos
An tð Þ ¼
An1 tð Þ … Ank tð Þ:
X
n1þ…þnk ¼n
If An(t ), norte 2 + admits a Maclaurin expansion An(t ) = ant + oh(t 2) (when n ≠ 1) or An(t ) = ant +
1 + oh(t 2) (when n = 1) where an = dAn(0)/dt, we can readily derive a0 = a2 = ς and a1 = −2ς
(García-Pelayo et al., 1993). Mientras tanto, we can know
An t þ dt
d
Þ − An tð Þ ¼
X∞
ak An−k tð Þdt:
k¼0
Ecuaciones 15 y 16 readily lead to
X∞
X∞
X
Yk
ak
k¼0
n¼0
n1þ…þnk ¼n−k
i¼1
!
Ani tð Þ
x n ¼
X∞
ak W x; t
d
Þk ;
k¼0
∂
∂t
W x; t
d
Þ ¼
P∞
n¼0 An(t )xn, X 2 [0, 1] denotes the generating function. Applying a trick
where W(X, t ) =
introduced in García-Pelayo et al. (1993), we define H(X) = ∂
∂tW(X, 0), which naturally leads to
∂
∂tW(X, t ) = H(W.(X, t )). Mientras tanto, h(X) = ς(1 − x )2 can be derived based on a0, a1, and a2
(García-Pelayo et al., 1993). Tomados juntos, tenemos
(35)
(36)
(37)
∂
∂t
W x; t
d
Þ ¼ ς 1 − W x; t
d
d
2:
Þ
Þ
(38)
Note that the initial condition is W(X, 0) = x since one neuron is activated at t *. Solving
Ecuación 38, we derive that
W x; t
d
Þ ¼
ς 1 − x
d
ς 1 − x
d
Þt
Þt þ 1
:
(39)
Por lo tanto, we have A0(t ) = W(0, t ) = ςt
ςtþ1, supporting a calculation of lifetime distribution PT(t )
lim
t →∞
PT tð Þ ¼ lim
t →∞
d
dt
W. 0; t
d
Þ ∼ t −2:
(40)
Following García-Pelayo et al. (1993), harris (1963), and Otter (1949), one can similarly
calculate
PS sð Þ ∼ s−3
2:
lim
s→∞
(41)
1165
Neurociencia en red
Theoretical foundations of studying criticality in the brain
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
.
/
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Cifra 3. Conceptual illustrations of the neural avalanche properties predicted by analytic theories of brain criticality. (A) Illustrations of the
framework to derive neural avalanche exponents in Box 3. (B) Illustrations of the scaling relation satisfied by neural avalanches under mean-
field assumptions. (C) Illustrations of the universal collapse shape of neural avalanches. The unscaled plot of t versus hS(t | t )i (upper left) y
the scaled plot hS(t | t )iT1−γ versus t
T and T 1−γ respectively serve as scale factors on x-axis
and y-axis to create a universal collapse shape. Mientras tanto, the symmetric collapse shape in Equation 51 (abajo a la izquierda) and the asymmetric
collapse shape controlled by skewness c in Equation 52 (bottom right) are also presented. (D) Autocorrelations and their decays in critical (izquierda)
and noncritical (bien) cases are shown. Autocorrelations are calculated after a ti 2 [0, t ) is randomly selected as a reference (upper left and
bien). Mientras tanto, the autocorrelation decays measured on tj 2 [de, t] is fitted to derive the coefficients χ (x-axis corresponds to ln(tj −ti
t )) and ξ
(x-axis corresponds to tj −ti
t ) en la ecuación 53 and Equation 55 (bottom left and right). One can see that autocorrelations in the critical case have
slower decays (smaller χ and ξ) than those in the noncritical case.
t (upper right) are shown for comparison. Here terms 1
There are three important things to remember. Primero, the lifetime exponent α = 2 and size
exponent β = 3
2 can only be treated as ideal exponents under mean-field assumptions of directed
percolation. There are numerous factors, such as granularity, network topology, and neural
dynamics variability, that can be considered in derivations to affect the detailed values of ava-
lanche exponents (Bonachela & Muñoz, 2009; Girardi-Schappo, 2021). En mesa 3, we summa-
rize the possible intervals of α and β in empirical neural data. Segundo, α and β alone are not
sufficient to verify the existence of brain criticality. Even when the actual values of α and β in
empirical data are exactly equal to theoretical predictions, they may still not satisfy the scaling
relation and universal collapse. Mientras tanto, as we shall discuss later, estimating α and β in prac-
tice is statistically error prone. Tercero, one cannot confirm or disprove a detail type of brain crit-
icality based on α and β unless additional information is provided. Although four types of brain
criticality exhibit (p.ej., OC) or approximately exhibit (p.ej., qC, SOC, and SOqC) directed perco-
lation behaviors under certain conditions, these preconditions are difficult to verify in practice.
Scaling Relation
In the previous section, we discuss how the neural avalanche lifetime and size distributions
exhibit power law properties when the brain is at the critical point (Hinrichsen, 2000;
Neurociencia en red
1166
Theoretical foundations of studying criticality in the brain
Mesa 3. Neural avalanche exponents with scaling relation in empirical data
Data type
LFP recordings filtered by spike sorting
Observed interval of α
a 2 [1.35, 2.67]
Observed interval of β
b 2 [1.3, 2.5]
Observed interval of γ
γ 2 [1.16, 1.48]
LFP recordings with thresholding
Ca and voltage imaging
a 2 [1.82, 2.84]
a 2 [2.15, 3.5]
b 2 [1.57, 2.59]
b 2 [1.5, 2.3]
γ 2 [1.12, 1.39]
γ 2 [1.75, 2.5]
Nota. The data is acquired from Girardi-Schappo (2021), dónde 45 experimental observations of neuronal avalanches reported by 30 studies are summarized.
These observations can be classified according to the recording techniques of neural avalanches. Detailed data classification criteria (p.ej., details of spike
sorting and thresholding) can be seen in Girardi-Schappo (2021). We only include the data where α, b, and γ are all recorded and satisfy the scaling relation
en la ecuación 45. For LFP recordings filtered by spike sorting, included observations are reported by Carvalho et al. (2021), Fontenele et al. (2019), Fosque et al.
(2021), Ma et al. (2019), Mariani et al. (2021), and Senzai et al. (2019). For LFP recordings with thresholding, included observations come from Mariani et al.
(2021) and Shew et al. (2015). For Ca and voltage imaging, observations are provided by Ponce-Alvarez et al. (2018) and Yaghoubi et al. (2018). A pesar de
numerous studies report neural avalanches in whole-brain imaging (p.ej., MEG, M/EEG, and invasive ECoG), these studies either do not report three exponents
together (Palva et al., 2013; Shriki et al., 2013; Zhigalov et al., 2015) or have not observed the scaling relation in Equation 45 (Varley et al., 2020). One can see
Girardi-Schappo (2021) for a summary of these results.
Larremore et al., 2012; Lübeck, 2004). Apart from lifetime T and size S, there are several other
quantities that characterize neural avalanches, such as area A (number of distinct active neu-
ron, measured as A ≃ hS(t )i where the expectation h·i is averaged across all neural ava-
lanches with the same lifetime T ) and radius exponent R (radius of gyration) (Lübeck, 2004;
Lübeck & Heger, 2003). En general, the corresponding probability distributions of these four
quantities decay algebraically
PX xð Þ∝x −λX ;
(42)
where random variable X 2 {S, t, A, R} can be an arbitrary quantity to characterize neural
avalanches. The avalanche exponent λX is defined according to the selected meaning of X
(p.ej., λT = 2 and λT = 3
2 under mean filed assumptions). Assuming that variables {S, t, A, R}
scale as a power of each other
X 0 ∝ X ψ
X0 X ; ∀X; X 0 2 S; t ; A; R
F
gramo;
we can derive the scaling relation from Equations 42 y 43
ψ
X 0X ¼
λX − 1
λX 0 − 1
:
If we let X 0 = A and X = T, we can specify Equation 44 como
γ ¼
α − 1
β − 1
;
(43)
(44)
(45)
where PT (t ) ∝ t−α, PS (s) ∝ s−β, and A ∝ T γ. Ecuación 45 leads to γ = 2 in the mean-field theory
of directed percolation. En mesa 3, one can see the possible interval of γ in empirical neural
datos. Ecuación 45 is widely used as a criterion to verify if the brain is at the critical point in
neuroscience studies (p.ej., Dalla Porta & Copelli, 2019; Fontenele et al., 2019; norte. Friedman
et al., 2012; Ponce-Alvarez et al., 2018). Once the scaling relation is confirmed among
observed neural avalanche exponents, it indicates key features of the universality class (please
note that α, b, and γ should be derived independently). For neuroscientists, the importance of
Ecuación 45 lies in that it provides extra verification of the validity of estimated neural ava-
lanche exponents. This verification is necessary given that neural avalanche exponent estima-
tion is frequently inaccurate (Fontenele et al., 2019). In Figure 3B, we illustrate the scaling
relation in Equation 45 under mean-field assumptions.
Neurociencia en red
1167
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
In Lübeck (2004), one can further learn about how brain criticality is mapped to a directed
percolation transition characterized by ordinary critical exponents. Mientras tanto, one can see
how to connect these neural avalanche exponents with second-order phase transition expo-
nents (Lübeck & Heger, 2003).
Universal Collapse Shape
Universal collapse with an implicit scaling function. Apart from the scaling relation discussed
arriba, the average temporal shape of bursts, a fundamental signature of avalanches
(Baldassarri, Colaiori, & Castellano, 2003; Laurson et al., 2013; Papanikolaou et al., 2011),
can also be used to verify the existence of brain criticality in a more precise manner. Este
approach has been previously applied on diverse physical systems, such as plastically deform-
ing crystals (Laurson & Alava, 2006) and Barkhausen noise (Mehta, Mills, Dahmen, & Sethna,
2002; Papanikolaou et al., 2011), and was recently introduced into neuroscience (Dalla Porta
& Copelli, 2019; Fontenele et al., 2019; norte. Friedman et al., 2012; Pausch, Garcia-Millan, &
Pruessner, 2020); Ponce-Alvarez et al., 2018). To understand this approach, let us step back to
the power relation in Equation 43 and specify that X 0 = A and X = T. These settings naturally
lead to
A ≃ S Tð Þ
h
i ≡
z
t
0
S t jTð
h
idt ∝ T γ;
Þ
(46)
where hS(t |t )i measures the averaged time-dependent avalanche size during an avalanche
and the expectation h·i is averaged across all neural avalanches with the same lifetime T. Equa-
ción 46 can be readily reformulated as
S t jTð
h
i ∝ T γ−1:
Þ
(47)
The general form of Equation 47 is usually given by Baldassarri et al. (2003), Laurson et al.
(2013), and Papanikolaou et al. (2011)
S tjTð
h
i ¼ T γ−1
Þ
h
;
(48)
(cid:5) (cid:6)
t
t
where H(·) denotes a universal scaling function. When the brain is at the critical point, all data
of hS (t | t )iT 1−γ is expected to collapse onto H(·) with reasonable errors (Baldassarri et al.,
2003; Laurson et al., 2013; Papanikolaou et al., 2011). Here the terminology “collapse onto”
means that all data generally exhibits a similar pattern in a plot of hS (t | t )iT 1−γ vs. t
t (p.ej., todo
data follows function H(·)). Mientras tanto, scaling function H(·) is expected to be a parabolic
función (Baldassarri et al., 2003; Laurson et al., 2013; Papanikolaou et al., 2011). By testing
these properties, neuroscientists can verify whether the brain is at criticality (p.ej., Dalla Porta &
Copelli, 2019; Fontenele et al., 2019; norte. Friedman et al., 2012; Ponce-Alvarez et al., 2018). En
Figura 3C, we graphically illustrate these properties.
Universal collapse with an explicit scaling function. Under specific conditions, researchers can
further consider an explicit form of scaling function H(·) (Laurson et al., 2013). Assuming that
the early-time growth of neural avalanches averagely follows a power law of time, one can
derive that hS (t | t )i ∝ t κ for certain t
≤ ε ≪ 1. Mientras tanto, one knows that hS ( εT | t )i ∝ T γ−1
t
should hold according to Equation 28. To ensure these two properties, one needs to have
hS (εT | t )i ∝ (εT )κ ∝ T γ−1, which readily leads to κ = γ − 1. Based on these derivations, uno
can know
h
S t jTð
i ∝ t γ−1;
Þ
t ≪ T :
(49)
1168
Neurociencia en red
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
.
/
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
To find an explicit form of H(·) that satisfies Equations 48 y 49, one can consider a
possible answer (Laurson et al., 2013)
(cid:5) (cid:6)
t
t
(cid:5)
1 − t
t
(cid:10)γ−1
(50)
t
t
h
¼
(cid:6)
(cid:9)
;
t )γ−1
which can be analytically derived by multiplying Equation 49 por (1 − t
is a term to characterize the deceleration at the ends of neural avalanches (Laurson et al.,
2013). Because γ = 2 is expected for critical neural avalanches under mean-field assump-
ciones, Ecuaciones 48 y 50 imply that
t )γ−1. Aquí (1 − t
S t jTð
h
Þ
(cid:6)
(cid:5)
i ∝ t 1 − t
t
:
(51)
This result is consistent with the prediction by the ABBM model in the limit of vanishing
drive rate and demagnetizing factor (norte. Friedman et al., 2012; Laurson et al., 2013).
A potential limitation of Equation 51 in applications lies in its internal symmetry property
(Laurson et al., 2013). Although avalanches under mean-field frameworks have a symmetric
average shape (norte. Friedman et al., 2012), it does not mean that symmetry generally holds in
real complex systems (Laurson et al., 2013). Applying Equation 48 on neural data, investigadores
may observe a nonstandard parabolic function H(·) with specific skewness. This does not
necessarily mean that neural dynamics is not at criticality. When neural avalanches are time
irreversible (this is generally true in the brain since the detailed balance of neural dynamics is
frequently broken; Lynn et al., 2021), one can consider small temporal asymmetry in the
collapse shape (Laurson et al., 2013). To characterize potential asymmetry, one can add a
correction term controlled by skewness degree c into Equation 51
S t jTð
h
Þ
(cid:5)
i ∝ t 1 − t
t
(cid:6)
(cid:9)
(cid:5)
(cid:6)
(cid:10)
1 − c
−
1
2
t
t
:
(52)
If c = 0, then Equation 52 reduces to Equation 51. De lo contrario, neural avalanches can have a
temporally asymmetric collapse shape with a positive (c > 0) or negative (C < 0) skewness
(Laurson et al., 2013). We suggest that Equation 52 may be more applicable to real data of
neural dynamics. In Figure 3C, we show examples of Equations 51 and 52.
Slow Decay of Autocorrelation
In applications, researchers can also consider a more practical verification of the potential
brain criticality. When the brain is at the critical point, a slow decay of autocorrelation is
expected to occur in neural avalanches, corresponding to long-range correlations (Dalla Porta
& Copelli, 2019; Erdos, Kruger, & Renfrew, 2018; Schaworonkow, Blythe, Kegeles, Curio, &
Nikulin, 2015; Smit et al., 2011). This slow decay property is initially derived from the power
law decay of autocorrelation, which can be analytically derived as a part of the scaling rela-
tion if ordinary critical exponents of directed percolation transition are considered (for details
see Girardi-Schappo, 2021). The power law decay is expressed as
(cid:6)
(cid:9)
(cid:10)
(cid:8)
(cid:8)
(cid:7)
Cov S ti jTð
Cov S ti jTð
ð
(cid:7)
Þ; S tj jT
Þ; S ti jTð
Þ
Þ
ln
(cid:5)
¼ −χln tj − ti
T
þ r;
(53)
where ti 2 [0, T ) is used as a reference and tj 2 [ti, T] traverses the entire interval
(Schaworonkow et al., 2015; Smit et al., 2011). According to the Wiener–Khinchin theorem,
coefficient χ is related to S(f ), the power spectrum of neural avalanches (notion f denotes
frequency) (Bak et al., 1987; Girardi-Schappo, 2021; Linkenkaer-Hansen, Nikouline, Palva,
Network Neuroscience
1169
l
D
o
w
n
o
a
d
e
d
f
r
o
m
h
t
t
p
:
/
/
d
i
r
e
c
t
.
m
i
t
.
/
t
/
e
d
u
n
e
n
a
r
t
i
c
e
-
p
d
l
f
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
n
e
n
_
a
_
0
0
2
6
9
p
d
.
/
t
f
b
y
g
u
e
s
t
t
o
n
0
7
S
e
p
e
m
b
e
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
& Ilmoniemi, 2001). One may expect S(f ) ∼ f −υ at the critical point, where χ = 1 − υ (Bak
et al., 1987; Girardi-Schappo, 2021; Linkenkaer-Hansen et al., 2001). The power law decay of
autocorrelation in Equation 53 breaks down when υ = 1, leading to infinitely long temporal
correlations. Therefore, χ 2 [0, ∞) in Equation 53 is expected to be sufficiently small. Certainly,
the actual value of χ may not be perfectly zero in empirical data. For instance, χ 2 [0.58 ±
0.23, 0.73 ± 0.31] and χ 2 [0.52 ± 0.35, 0.81 ± 0.32] are observed in spontaneous alpha
oscillations in MEG and EEG data, respectively (Linkenkaer-Hansen et al., 2001).
Apart from verifying power law decay directly, one can also consider the exponential
decay, which is active in neuroscience as well (Miller & Wang, 2006; Pausch et al., 2020;
Wilting & Priesemann, 2019). The exponential decay can described by
(cid:7)
Cov S ti jTð
(cid:7)
Þ; S tj jT
(cid:8)
(cid:8)
(cid:7)
¼ −ξ Cov S ti jTð
(cid:7)
Þ; S tj jT
(cid:8)
(cid:8)
:
∂
∂t
Equation 54 directly leads to
(cid:9)
ln
(cid:7)
Cov S ti jTð
Cov S ti jTð
ð
(cid:7)
Þ; S tj jT
Þ; S ti jTð
(cid:8)
(cid:8)
(cid:10)
Þ
Þ
(cid:6)
(cid:5)
¼ −ξ tj − ti
T
þ r:
(54)
(55)
l
D
o
w
n
o
a
d
e
d
f
r
o
m
h
t
t
p
:
/
/
d
i
r
e
c
t
.
m
i
t
.
/
/
t
e
d
u
n
e
n
a
r
t
i
c
e
-
p
d
l
f
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
n
e
n
_
a
_
0
0
2
6
9
p
d
.
t
/
f
b
y
g
u
e
s
t
t
o
n
0
7
S
e
p
e
m
b
e
r
2
0
2
3
The exponential decay can be seen in the dynamics with short-term correlations (i.e., corre-
lations have a characteristic time scale). Mathematically, the exponential decay can be related
R ∞
to power law decay in a form of x−y = Γ(y )
0 zy−1 exp(−xz)dz, where Γ(·) denotes the Gamma
function. When ξ 2 [0, ∞) is sufficiently small, Equation 55 can be treated as a looser criterion
that approximately verifies the slow decay of autocorrelation and may be more applicable to
nonstandard brain criticality (e.g., qC and SOqC) (Wilting & Priesemann, 2019). Despite of its
practicality, this looser criterion should be used with caution since it is not analytically derived
from criticality theories.
In Figure 3D, we illustrate examples of autocorrelation slow decay in critical cases and
compare them with noncritical cases. Compared with other properties previously discussed,
a slow autocorrelation decay can be readily verified by conventional data fitting. However,
we need to note that one should not confirm or reject the possibility of brain criticality only
based on the decay characteristic of autocorrelation in Equations 53–55. This is because
Equations 53–55 only serve as the approximate descriptions of long-range correlations at
criticality. The strict criterion χ, ξ → 0 is rarely seen in empirical data while the determination
of whether χ and ξ are sufficiently small in the looser criterion is relatively subjective.
In summary, we have reviewed the physical foundations of identifying and characterizing
criticality in the brain. Based on these analytic derivations, we attempt to present systematic
explanations of what is brain criticality and how to identify potential criticality in neural
dynamics. Nevertheless, physical theories alone are not sufficient to support neuroscience
studies because the implementation of these theories on empirical data is even more challeng-
ing than the theories themselves. To overcome these challenges, one needs to learn about
statistic techniques to computationally estimate brain criticality from empirical data.
BRAIN CRITICALITY: STATISTIC TECHNIQUES
While most properties of neural avalanches analytically predicted by the physical theories of
brain criticality can be estimated by conventional statistic techniques, there exist several
properties that frequently imply serious validity issues and deserve special attention. Below,
we discuss them in detail.
Network Neuroscience
1170
Theoretical foundations of studying criticality in the brain
Estimating Neural Avalanche Exponents
Perhaps the estimation of neural avalanche exponents from empirical data is the most error-
prone step in brain criticality analysis. The least-square approach is abused in fitting power
law data and frequently derives highly inaccurate results (Clauset, Shalizi, & Newman, 2009;
Virkar & Clauset, 2014). To derive neural avalanche exponents α and β in Equation 25 with
reasonable errors, one needs to consider the maximum likelihood estimation (MLE)
approach and corresponding statistic tests (see MLE on unbinned data; Clauset et al.,
2009; and binned data; Virkar & Clauset, 2014). Taking the avalanche size distribution as
an instance, the estimator ^β of distribution exponent β is expected to maximize the log-
likelihood function
L βð Þ ¼ −n ln ζ β; s0
ð
½
(cid:2) − β
Þ
Xn
i¼1
ln sið Þ;
L βð Þ ¼ n β − 1
ð
Þ ln b0 þ
Xk
(cid:7)
hi ln b
i¼1
1−β
i
− b
1−β
iþ1
(cid:8)
:
(56)
(57)
Here Equation 56 and Equation 57 denote the log-likelihood functions on unbinned and
binned data, respectively. Function ζ(·, ·) is the generalized zeta function (Bauke, 2007;
Clauset et al., 2009). Notion s denotes avalanche size samples in Equation 33 and Equation 41
(Clauset et al., 2009). Notion b denotes bin boundaries defined on these samples and h
counts the number of samples within each bin (Virkar & Clauset, 2014). Notions s0 and b0
are the lower cutoffs of unbinned and binned power law distributions (Clauset et al., 2009;
Virkar & Clauset, 2014). They are necessary because few empirical data exhibits power law
properties on the entire distribution (Clauset et al., 2009). Notions n and k measure the num-
bers of samples and bins above cutoffs, respectively (Clauset et al., 2009; Virkar & Clauset,
2014). To estimate ^β precisely, researchers are suggested to follow several indispensable steps
(Clauset et al., 2009; Virkar & Clauset, 2014): (1) for each potential choice of s0 or b0, estimate
the power law model on the distribution tail above the cutoff. Compute the Kolmogorov–
Smirnov (KS) goodness-of-fit statistic between the cumulative probability distributions of
power law model and empirical data. Find the ideal choice of s0 or b0 that minimizes KS sta-
tistic; (2) derive the corresponding estimator ^β and KS statistic based on the chosen cutoff; (3)
use the semiparametric bootstrap to generate numerous synthetic data distributions that follow
the estimated power law model above the cutoff but follow the empirical data distribution
below the cutoff. Estimate new power law models on these synthetic data distributions and
measure the goodness-of-fit by KS statistic. Define a p value, the fraction of these KS statistics
in step 3 that are no less than the KS statistic in step 2. Rule out the estimated power law model
in steps 1–2 if p < 0.1 (conservative criterion). Apart from these necessary steps, one can fur-
ther consider Vuong’s likelihood ratio test for alternative distribution checking (Clauset et al.,
2009; Virkar & Clauset, 2014; Vuong, 1989) and information loss measurement of binning
approach (Virkar & Clauset, 2014). During the above process, we measure the goodness-of-
fit by KS statistic instead of the well-known χ2 statistic because the latter has less statistic power
(Bauke, 2007; Clauset et al., 2009; Virkar & Clauset, 2014). Meanwhile, KS statistic is mea-
sured on cumulative probability distributions rather than probability distributions to control the
effects of extreme values in empirical data (Clauset et al., 2009; Virkar & Clauset, 2014).
Except for the above approach (Clauset et al., 2009; Virkar & Clauset, 2014), one can also
consider the BIC method (for unbinned data) (Schwarz, 1978) and the RT method (for binned
data) (Reiss & Thomas, 2007) for comparisons. In practice, the approaches proposed by
Clauset et al. are more robust (Clauset et al., 2009; Virkar & Clauset, 2014) and have attracted
Network Neuroscience
1171
l
D
o
w
n
o
a
d
e
d
f
r
o
m
h
t
t
p
:
/
/
d
i
r
e
c
t
.
m
i
t
.
/
t
/
e
d
u
n
e
n
a
r
t
i
c
e
-
p
d
l
f
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
n
e
n
_
a
_
0
0
2
6
9
p
d
.
t
/
f
b
y
g
u
e
s
t
t
o
n
0
7
S
e
p
e
m
b
e
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
numerous follow-up studies for improvements (e.g., Deluca & Corral, 2013; Marshall et al.,
2016; Yu, Klaus, Yang, & Plenz, 2014).
Estimating Universal Collapse Shape
Another error-prone step is the calculation and evaluation of the universal collapse shape,
which is closely related to scaling relation analysis. Deriving the collapse shape from empir-
ical data may be problematic because the goodness evaluation of collapse shape is rather sub-
jective (e.g., depends on personal opinions about whether all data follows function H(·) in
Equation 48) in most cases (Marshall et al., 2016). Although important efforts have been
devoted to quantify if a given dataset exhibits shape collapse (Bhattacharjee & Seno, 2001;
Shaukat & Thivierge, 2016), common approaches in practice still depend on specific shape
collapse algorithms that search potential scaling parameters (e.g., γ in Equation 48) in a data-
driven manner (Marshall et al., 2016). In these algorithms, thresholding on neural avalanches
before analyzing the shape collapse is a standard preprocessing scheme to control noises (e.g.,
set an avalanche size threshold and remove all data below the threshold) (Marshall et al.,
2016; Papanikolaou et al., 2011). While experimental noises are partly limited, unexpected
excursions of scaling parameters away from theoretical predictions may occur after threshold-
ing as well (Laurson, Illa, & Alava, 2009). To our best knowledge, the effects of thresholding on
brain criticality analysis are nonnegligible. Although being highly practical, thresholding may
lead to significant transient effects to cloud the true scaling property (Villegas, di Santo,
Burioni, & Muñoz, 2019). Therefore, any qualitative evaluation of collapse shape after thresh-
olding is questionable regardless of its practicability. Although an ideal approach requires fur-
ther explorations, we suggest researchers to consider the following methods: (1) estimate γ by
area fitting (e.g., follow Equation 47 in scaling relation analysis) and collapse shape fitting (e.g.,
follow Equation 48 in collapse shape analysis), respectively; (2) compare between γ derived by
these two kinds of fitting and measure the difference. Search for a threshold that minimizes the
difference (e.g., makes variation amplitude < 1%) and maintains a reasonable sample size
(e.g., maintains > 80% muestras); (3) given the chosen threshold and corresponding γ, measure
the difference (p.ej., the dynamic time warping; Keogh & Pazzani, 2001) between hS(t |t )iT1−γ
derived on neural avalanches with different lifetime T in the plot of hS(t |t )iT1−γ vs. t
t. Denote
the shape collapse error as the averaged difference. Combining these three steps, investigadores
may partly avoid the errors implied by subjective judgment. Similar ideas can be seen in
Marshall et al. (2016).
Estimating the Slow Decay of Autocorrelation
Finalmente, the analysis of slow decay of autocorrelation is also error-prone in practice. A pesar de
this approach is practical and has been extensively applied (p.ej., Pausch et al., 2020; Wilting
& Priesemann, 2019), the criterion to determine if the decay is truly slow (p.ej., χ > 0 in Equa-
ción 54 and ξ > 0 en la ecuación 54 are sufficiently small) remains ambiguous. A fixed criterion
(p.ej., χ, ξ < 0.5) may serve as an explicit condition of a slow decay. However, this presupposed
criterion may deviate from real situations. For instance, the baseline of decay rate in a non-
critical brain may be essentially high (e.g., χ, ξ > 10). Even though the decay rate drops sig-
nificantly when the brain becomes critical (p.ej., χ, ξ ≃ 1), the presupposed criterion is still
unsatisfied and leads to unnecessary controversies on criticality hypothesis. Given that ξ is
principally independent from spatial subsampling on neurons or brain regions at criticality
(Pausch et al., 2020), we suggest researchers to consider the following approaches: (1) do spa-
tial subsampling in both critical and non-ritical brains to derive two groups of χ or ξ (one group
for criticality and another group for noncriticality); (2) use appropriate statistic tests (p.ej.,
Neurociencia en red
1172
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
choose t test; Kanji, 2006; Kolmogorov–Smirnov test; Berger & zhou, 2014; or Wilcoxon-
Mann-Whitney test; Fay & Proschan, 2010; according to sample distribution properties) to ver-
ify if two groups of χ or ξ belong to different distributions. Test if the expectation and variance
of χ or ξ drops significantly from the noncritical group to the critical group according to certain
effect sizes.
En resumen, statistic techniques bridge between brain criticality theory and empirical data.
Sin embargo, misconception and misuse of statistical analyses of neural avalanche properties still
occasionally appear in practice. Although existing techniques remain imperfect in brain crit-
icality analysis, we wish that our discussion may inspire future studies.
BRAIN CRITICALITY AND OTHER NEUROSCIENCE THEORIES
Ever since brain criticality was introduced into neuroscience, it is frequently speculated as
contradictory with other traditional neuroscience hypotheses, such as the conjectured hier-
archical processing characteristic of neural information (Felleman & VanEssen, 1991) y el
asynchronous-irregular characteristic of neural dynamics (p.ej., neurons spike independently
in Poisson manners; Burns & Webb, 1976; Softky & Koch, 1993; piedra, Gossen, & jones,
2005). Mientras tanto, the differences between brain criticality and scale-free neural dynamics
(Chialvo, 2010; Él, 2014; Martinello et al., 2017) are frequently neglected. Before we put an
end to our review, we discuss the relations between brain criticality and these neuroscience
teorías.
Brain Criticality and Hierarchical Processing
The hierarchical processing of neural information (Felleman & VanEssen, 1991) is initially
speculated to contradict critical neural dynamics since hierarchical topology has not been
used as an explicit condition to analytically derive criticality (p.ej., see derivations in di Santo
et al., 2017; García-Pelayo et al., 1993; Gros, 2010; harris, 1963; Janowsky & Laberge, 1993;
Larremore et al., 2012; D. S. Lee et al., 2004; Otter, 1949; robinson, 2021). On the contrary,
random graphs without strict hierarchical structures seem to be more widespread in criticality
derivations. Recientemente, this speculation has been challenged by numerous discoveries of the
facilitation effects of hierarchical modular structures on critical dynamics (mi. j. Friedman &
Landsberg, 2013; Kaiser & Hilgetag, 2010; Rubinov et al., 2011; S. Wang & zhou, 2012).
Mientras tanto, computational analysis suggests that information transmission in standard feed-
forward networks is maximized by critical neural dynamics (ruega & Plenz, 2003). Parallel
to neuroscience, a recent machine-learning study empirically observes and analytically dem-
onstrates that artificial neural networks, a kind of hierarchical structure, self-organize to criti-
cality during learning (Katsnelson, Vanchurin, & Westerhout, 2021). Por lo tanto, brain criticality
is not necessarily contradictory with hierarchical information processing, yet more analyses
are required to understand how brain criticality affects hierarchical processing schemes.
Brain Criticality and Asynchronous-Irregular Characteristic
Brain criticality and the asynchronous-irregular (AI) characteristic may correspond to distinct
encoding schemes in the brain (Girardi-Schappo et al., 2021; Wilting & Priesemann, 2019).
While AI characteristic can minimize redundancy (Atick, 1992; Barlow et al., 1961; Campana &
Sejnowski, 1997; Van Hateren & van der Schaaf, 1998) to improve neural encoding (Van
Vreeswijk & Sompolinsky, 1996), brain criticality may optimize encoding performance
utilizing a series of reverberations of neural activities (Bertschinger & Natschläger, 2004;
Boedecker, Obst, Lizier, Mayer, & Asada, 2012; Del Papa, Priesemann, & Triesch, 2017;
Neurociencia en red
1173
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Haldeman & ruega, 2005; Kinouchi & Copelli, 2006; mostrar & Plenz, 2013; X. R. Wang, Lizier,
& Prokopenko, 2011). The coexistence of empirical evidence of AI and brain criticality
characteristics initially confuses researchers since these characteristics are hypothesized as
contradictory to each other (Girardi-Schappo et al., 2021; Wilting & Priesemann, 2019). En
experimentos, AI characteristic is supported by small correlations between the spike rates of
different neurons in cortical microcircuits (cohen & Kohn, 2011; Ecker et al., 2010) and expo-
nential distributions of interspike intervals (Carandini & stevens, 2004; Kara, Reinagel, & Reid,
2000) while brain criticality characteristic is observed in neural dynamics recorded from
multiple species (p.ej., awake monkeys; Petermann et al., 2009, anesthetized rats; Gireesh &
Plenz, 2008, slices of rat cortices; ruega & Plenz, 2003; Shew et al., 2009, and humans; Cabello
et al., 2008). A recent study demonstrates that cortical spikes may propagate at somewhere
between perfect criticality (p.ej., OC or SOC depending on whether underlying mechanisms
are exogenous or endogenous) and full irregularity (Wilting & Priesemann, 2019), similar to
the cases of qC and SOqC. Mientras tanto, it is known that stimulus drives suppress irregularity in
neural activities (Molgedey et al., 1992). These results imply that brain criticality may not nec-
essarily contradict AI characteristic. On the contrary, they may coexist when stimulus drives
are too weak to disturb brain criticality (p.ej., OC or SOC) and suppress AI characteristics. En
our previous discussions, we have analytically proven that neural avalanche exponents, el
fundamental properties of brain criticality, can still be derived under the condition of indepen-
dent neuron activation, a key feature of AI characteristics (Wilting & Priesemann, 2019). Este
result suggests that brain criticality and AI characteristics do not contradict each other. Como para
the case where stimulus drives are nonnegligible, a recent study presents an elegant theory to
prove that two homeostatic adaptation mechanisms (es decir., the short-term depression of inhibi-
tion and the spike-dependent threshold increase) enable synaptic excitation/inhibition bal-
ance, AI characteristics, and SOqC to appear simultaneously in the same neural dynamics
(Girardi-Schappo et al., 2021). Similarmente, it is suggested that neural dynamics with criticality
or with AI characteristics can be generated by the same neural populations if the synaptic
excitation/inhibition balance is fine tuned appropriately (j. li & mostrar, 2020).
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Brain Criticality and Power Law Behaviors in Neural Dynamics
Neural dynamics with power law behaviors is a necessary but insufficient condition of brain
criticidad. This property is frequently neglected in practice. Power law behaviors are wide-
spread in nature because they can be generated by diverse mechanisms, such as exponential
curve summation and preferential attachment (Mitzenmacher, 2004; Reed & abrazos, 2002). Él
has been reported that the aggregate behaviors of noncritical stochastic systems may also cre-
ate scale-free dynamics within a limited range (Touboul & Destexhe, 2010, 2017). In the brain,
the generic scale-free properties can be implied by neutral dynamics, a kind of dynamics
where the population size of neutral individuals (or dynamically homogeneous individuals)
does not tend to increase or decrease after adding a new individual that is neutral to existing
unos (see neutral theories for further explanations; Blythe & McKane, 2007; Liggett, 2006).
This generic property can generate power law neural avalanches without criticality (Martinello
et al., 2017). Mientras tanto, bistability phenomena, a kind of fine-tuned or self-organized
discontinuous phase transitions with limit cycles rather than critical points, can also create
neural dynamics with power law properties (Buendía et al., 2020a; Cocchi et al., 2017; di
Santo, Burioni, Vezzani, & Muñoz, 2016). Como consecuencia, we emphasize that neural avalanche
exponents alone are insufficient to prove or disprove any brain criticality hypothesis. Estos
power law exponents are meaningless for brain criticality hypothesis unless they satisfy the
scaling relation.
Discontinuous phase transition:
A kind of phase transition where the
first-order derivative of the order
parameter diverges when transition
happens.
Neurociencia en red
1174
Theoretical foundations of studying criticality in the brain
BRAIN CRITICALITY: CONCLUSIONS ON CURRENT PROGRESSES
AND LIMITATIONS
Given what have been reviewed above, we arrive at a point to conclude on the current pro-
gresses and limitations in establishing theoretical foundations of different types of brain criti-
cality, eso es, ordinary criticality (OC), quasi-criticality (qC), self-organized criticality (SOC),
and self-organized quasi-criticality (SOqC). As we have suggested, an inescapable cause of
various controversies is the nontriviality of physical theories that analytically derive brain crit-
icality and statistical techniques that estimate brain criticality from empirical data. Immoderate
omitting of these theoretical foundations, especially their imperfection, in practice may lead to
confusions on the precise meaning, identification criteria, and biological corollaries of brain
criticidad. To address these problems, we have introduced the mainstream theoretical founda-
tions of brain criticality, reformulated them in the terminology of neuroscience, and discussed
their error-prone details.
Thanks to the increasing efforts devoted to improving theoretical frameworks of criticality in
the brain, researchers have seen substantial progresses in explaining various important neuro-
science problems, including but not limited to efficient cortical state transitions (Fontenele
et al., 2019), dynamic range maximization in neural responses (Kinouchi & Copelli, 2006;
Shew et al., 2009), and optimization of information transmission and representation (mostrar
et al., 2011). These advances have been comprehensively reviewed by existing works (ruega,
2007; Chialvo, 2010; Cocchi et al., 2017; Hesse & Bruto, 2014; Muñoz, 2018; mostrar & Plenz,
2013) and are not discussed in detail in our review. The benefits of studying brain criticality, como
we have suggested, lay in the possibility to analyze brain function characteristics with numer-
ous statistical physics theories relevant to brain criticality, such as directed percolation
(Hinrichsen, 2000; Lübeck, 2004), conserved directed percolation (Bonachela et al., 2010;
Bonachela & Muñoz, 2008), and dynamical percolation theories (Bonachela et al., 2010; Steif,
2009). These theories characterize the brain as a physical system with avalanche behaviors,
enabling researchers to analyze various propagation, synchronization, and correlation prop-
erties of neural dynamics (p.ej., continuous phase transitions). These properties intrinsically
shape neural information processing (p.ej., encoding; Bertschinger & Natschläger, 2004;
Boedecker et al., 2012; Del Papa et al., 2017; Haldeman & ruega, 2005; Kinouchi & Copelli,
2006; mostrar & Plenz, 2013; X. R. Wang y cols., 2011; transmission; Shew et al., 2011, y
memory; Haldeman & ruega, 2005; Krotov & Hopfield, 2020) and can be readily recorded
in neuroscience experiments. Por lo tanto, the nonequilibrium dynamic processes and potential
criticality defined by statistical physics theories are highly applicable to characterizing brain
funciones. As we have discussed in Figure 2, researchers can consider diverse brain criticality
phenomena in neural dynamics by defining different control (p.ej., the balance between excit-
atory and inhibitory neurons; Hardstone et al., 2014; Poil et al., 2012) and order (p.ej., active
neuron density; Dalla Porta & Copelli, 2019) parámetros, corresponding to multifarious bio-
logical mechanisms underlying neural dynamics (p.ej., synaptic depression; Levina et al.,
2007). Mientras tanto, the definition of neural avalanches can flexibly change from neural spikes
and local field potentials to global cortical oscillations. The flexibility of brain criticality and
neural avalanche definitions enables researchers to analyze different functional properties on
distinct organizational levels in the brain.
The limited theoretical foundations of brain criticality in the brain, sin embargo, have become
irreconcilable with their increasingly widespread applications. Although the analytic theories
of brain criticality have solid physics backgrounds, they needlessly become black boxes for
neuroscientists in practice. Por un lado, the details of brain criticality theory frequently
Neurociencia en red
1175
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
/
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
experience immoderate neglecting in neuroscience studies. Por otro lado, to our best
conocimiento, there is no accessible and systematic introduction of the statistical physics foun-
dations of brain criticality in the terminology of neuroscience yet. These obstacles severely
impede neuroscientists from comprehensively understanding brain criticality, eventually
motivating us to present this review. When we turn to bridging between brain criticality
theories and experiments, one can find nonnegligible gaps separating between theories
and experiments. Although numerous biological factors (p.ej., neural plasticity; De Arcangelis
et al., 2006; Levina et al., 2007, 2009; membrane potential leakage; Levina et al., 2007;
Millman et al., 2010; Rubinov et al., 2011; Stepp et al., 2015; retro-synaptic signals;
Hernandez-Urbina & Herrmann, 2017; spatial heterogeneity; Girardi-Schappo, Bortolotto,
Gonsalves, Pinto, & Tragtenberg, 2016; Moretti & Muñoz, 2013; and refractory period;
Oscuro et al., 2021; Williams-García et al., 2014) have been considered in brain criticality
characterization, existing theories more or less suffer from deviations from actual neural sys-
tem properties. Por ejemplo, the requirements of conserved neural dynamics and an infinite
timescale separation between the dissipation and drive processes required by SOC may not
be biologically realistic (Muñoz, 2018). The implicit requirement of a sufficiently large system
size by the mean-field theories of brain criticality may not always be satisfied during neural
avalanche recording, implying nonnegligible finite size effects (Girardi-Schappo, 2021).
Mientras tanto, to precisely verify the existence of a detailed type of brain criticality (p.ej., estafa-
firm the actual universality class) in empirical neural data is principally infeasible. As we have
explained, the common criteria used for brain criticality hypothesis verification, such as neu-
ral avalanche exponents (Bauke, 2007; Clauset et al., 2009; Deluca & Corral, 2013; marshall
et al., 2016; Yu et al., 2014; scaling relation; Lübeck, 2004; Lübeck & Heger, 2003; universal
collapse shape; Bhattacharjee & Seno, 2001; Laurson et al., 2009; Marshall et al., 2016;
Papanikolaou et al., 2011; and slow decay of auto-correlation; Pausch et al., 2020; Wilting
& Priesemann, 2019), are derived according to directed percolation theory under mean-field
suposiciones. Among four types of brain criticality in absorbing phase transitions, only OC
originally belongs to directed the percolation universality class, while qC, SOC and SOqC
conditionally exhibit directed percolation behaviors. In most cases, one can only verify if
the brain is plausibly at criticality (p.ej., whether neural avalanches obey universal collapse
and have the power law exponents that satisfy the scaling relation). When observed neural
avalanche exponents depart from their mean-field approximation results but still satisfy the
scaling relation, there may exist an OC phenomenon affected by nonmean-field factors (p.ej.,
network topology; Girardi-Schappo, 2021) or exist a certain qC, SOC, or SOqC phenomenon
caused by diverse mechanisms. Additional information of neural dynamics properties is inev-
itably required to determine the category belonging to the hypothesized brain criticality,
which poses daunting challenges to neuroscience experiment designs. Además, the poten-
tial validity issues of applying the theoretical tools derived from directed percolation theory to
verify brain criticality in synchronous phase transitions deserve special attention (for similar
opinions see Dalla Porta & Copelli, 2019). It remains controversial if absorbing and synchro-
nous phase transitions robustly share specific features (see reported similarities; Buendía
et al., 2021; di Santo et al., 2018; Fontenele et al., 2019; Girardi-Schappo et al., 2021;
and differences; Buendía et al., 2021; Fontenele et al., 2019; Girardi-Schappo et al.,
2021). Any speculated relations between these two kinds of critical phenomena should be
tested with caution. Además, statistic techniques to estimate and verify brain criticality
from empirical data are yet imperfect. The estimation of some properties of neural
avalanches is error prone in practice and may lead to serious validity issues. Although we
suggest compromised solutions to these issues, more optimal approaches are required in
future studies.
Neurociencia en red
1176
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
.
/
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
BRAIN CRITICALITY: SUGGESTIONS OF FUTURE DIRECTION
We submit that this review not only summarizes the latest developments in the field of study-
ing criticality in the brain, but also serves as a blueprint for further explorations. Abajo, nosotros
offer concrete recommendations of future directions.
Primero, we suggest researchers carefully rethink the theoretical foundations of criticality in the
cerebro. Immoderately omitting these foundations in neuroscience needlessly muddies an
already complex scientific field and leads to potential validity issues. While we have presented
a self-contained framework of brain criticality to characterize neural dynamics as a physical
system with avalanches, plentiful details are uncovered in this article (p.ej., the Landau–
Ginzburg theory; di Santo et al., 2018) because the statistical physics theories of brain criticality
are essentially grand. We recommend researchers to further improve our work and explore a
more accessible and systematic reformulation of related physics theories, such as directed per-
colation, conserved directed percolation, dynamic percolation, and nonequilibrium dynamics,
in the context of neuroscience. Además, we note that these theories are not initially proposed
for brain analysis. It is normal to see gaps between these theories and real situations of the
cerebro. We urge researchers to develop new variants of criticality formalism that is more appli-
cable to the brain or even explore new universality classes of continuous phase transitions.
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
t
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
.
/
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Segundo, neuroscience is in urgent need of new physical theories and statistical techniques
to bridge between brain criticality hypotheses and experiments. Although existing theories and
techniques have become increasingly widespread and cover most of the pivotal details of
brain criticality, there remain various limitations, as we have suggested. Específicamente, we sug-
gest five potential directions to resolve these problems: (1) combine brain criticality theories
with large-scale neural dynamics recording or computation to include more realistic biological
details into brain criticality theories and establish a closer connection with experimental obser-
vaciones; (2) try to summarize, standardize, and subdivide these theories according to the con-
crete biological meanings of brain criticality phenomena, prerequisites of model definitions,
and scopes of application—try to avoid abusing or misusing of different brain criticality the-
ories; (3) develop open-source toolboxes of theoretical models and statistical techniques to
routinize brain criticality analysis in neuroscience studies (one can see existing efforts to
achieve this objective; Marshall et al., 2016); (4) establish open-source, multispecies, y
large-scale datasets of neural dynamics recorded from both critical and noncritical
brains—validate different statistic techniques of brain criticality estimation and testing on these
datasets and, more importantly, confirm appropriate baselines to define the criteria of brain
criticality identification (see notable contributions in Girardi-Schappo, 2021); (5) explorar
new nonequilibrium statistical physics theories for synchronous phase transitions or analyti-
cally verify the theoretical validity of directed percolation formulation of synchronous phase
transitions.
Tercero, parallel to neuroscience, the discoveries of critical phenomena in other learning and
computation systems also merit attention. Learning or computing at the edge of chaos has
been proven as a mechanism to optimize the performance of learners (p.ej., recurrent neural
redes; Bertschinger & Natschläger, 2004). The well-known residual connections can con-
trol the performance degradation of artificial neural networks because they enable networks to
self-organize to criticality between stability and chaos to preserve gradient information flows
(Cual & Schoenholz, 2017). It is recently demonstrated that any artificial neural network gen-
erally self-organizes to criticality during the learning process (Katsnelson et al., 2021). En el
future, it would be interesting to explore whether information processing processes in brains
and artificial neural networks can be universally characterized by a unified criticality theory.
Neurociencia en red
1177
Theoretical foundations of studying criticality in the brain
En general, we anticipate the potential of well-validated studies of criticality in the brain to
greatly deepen our understanding of neural dynamics characteristics and their roles in neural
procesamiento de información. Laying solid theoretical foundations of studies is the most effective and
indispensable path to contributing to this booming research area.
EXPRESIONES DE GRATITUD
Authors are grateful for discussions and assistance of Drs. Yaoyuan Wang and Ziyang Zhang
from the Laboratory of Advanced Computing and Storage, Instituto Central de Investigaciones, 2012
Laboratorios, Huawei Technologies Co.. Limitado., Beijing, 100084, Porcelana.
CONTRIBUCIONES DE AUTOR
Yang Tian: Conceptualización; Análisis formal; Investigación; Metodología; Visualización;
Escritura – borrador original; Escritura – revisión & edición. Zeren Tan: Análisis formal; Metodología.
Hedong Hou: Análisis formal; Metodología. Guoqi Li: Validación; Escritura – revisión & edición.
Aohua Cheng: Análisis formal; Escritura – revisión & edición. Yike Qiu: Validación; Writing –
revisar & edición. Kangyu Weng: Validación; Escritura – revisión & edición. Chun Chen:
Validación; Escritura – revisión & edición. Pei Sol: Conceptualización; Administración de proyecto;
Supervisión; Validación; Escritura – borrador original; Escritura – revisión & edición.
INFORMACIÓN DE FINANCIACIÓN
Pei Sol, The Artificial and General Intelligence Research Program of Guo Qiang Research
Institute at Tsinghua University, Award ID: 2020GQG1017.
REFERENCIAS
Abbott, l. F. (2008). Theoretical neuroscience rising. Neurona, 60(3),
489–495. https://doi.org/10.1016/j.neuron.2008.10.019,
PubMed: 18995824
Acebrón, j. A., Bonilla, l. l., Pérez Vicente, C. J., Ritort, F., &
Spigler, R. (2005). The Kuramoto model: A simple paradigm for
synchronization phenomena. Reseñas de Física Moderna, 77,
137–185. https://doi.org/10.1103/RevModPhys.77.137
Antonopoulos, C. GRAMO. (2016). Dynamic range in the C. elegans brain net-
trabajar. Chaos: An Interdisciplinary Journal of Nonlinear Science, 26(1),
013102. https://doi.org/10.1063/1.4939837, PubMed: 26826854
Arenas, A., Díaz-Guilera, A., Kurths, J., Moreno, y., & zhou, C.
(2008). Synchronization in complex networks. Physics Reports,
469(3), 93–153. https://doi.org/10.1016/j.physrep.2008.09.002
Atick, j. j. (1992). Could information theory provide an ecological
theory of sensory processing? Red: Computation in Neural
Sistemas, 3(2), 213–251. https://doi.org/10.1088/0954-898X_3_2_009
Bak, PAG. (2013). How nature works: The science of self-organized
criticidad. Nueva York, Nueva York: Ciencia Springer & Medios comerciales.
Bak, PAG., & Sneppen, k. (1993). Punctuated equilibrium and criticality in
a simple model of evolution. Physical Review Letters, 71(24),
4083–4086. https://doi.org/10.1103/ PhysRevLett.71.4083,
PubMed: 10055149
Bak, PAG., Espiga, C., & Wiesenfeld, k. (1987). Self-organized criticality:
An explanation of the 1/f noise. Physical Review Letters, 59(4),
381–384. https://doi.org/10.1103/PhysRevLett.59.381, PubMed:
10035754
Baldassarri, A., Colaiori, F., & Castellano, C. (2003). Average shape
of a fluctuation: Universality in excursions of stochastic pro-
cesses. Physical Review Letters, 90(6), 060601. https://doi.org
/10.1103/PhysRevLett.90.060601, PubMed: 12633282
Barlow, h. B., et al. (1961). Possible principles underlying the trans-
formation of sensory messages. Sensory Communication, 1(1).
bassett, D. S., Greenfield, D. l., Meyer-Lindenberg, A., Weinberger,
D. r., moore, S. w., & bullmore, mi. t. (2010). Efficient physical
embedding of topologically complex information processing
networks in brains and computer circuits. PLoS computacional
Biología, 6(4), e1000748. https://doi.org/10.1371/journal.pcbi
.1000748, PubMed: 20421990
Bauke, h. (2007). Parameter estimation for power-law distributions
by maximum likelihood methods. The European Physical Journal
B, 58(2), 167–173. https://doi.org/10.1140/epjb/e2007-00219-y
ruega, j. METRO. (2007). How to build a critical mind. Nature Physics,
3(12), 835. https://doi.org/10.1038/nphys799
ruega, j. METRO., & Plenz, D. (2003). Neuronal avalanches in neocor-
tical circuits. Revista de neurociencia, 23(35), 11167–11177.
https://doi.org/10.1523/ JNEUROSCI.23-35-11167.2003,
PubMed: 14657176
ruega, j. METRO., & Hora, norte. (2012). Being critical of criticality in the
cerebro. Frontiers in Physiology, 3, 163. https://doi.org/10.3389
/fphys.2012.00163, PubMed: 22701101
Campana, A. J., & Sejnowski, t. j. (1997). The “independent components”
of natural scenes are edge filters. Vision Research, 37(23),
Neurociencia en red
1178
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
t
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
3327–3338. https://doi.org/10.1016/S0042-6989(97)00121-1,
PubMed: 9425547
Berger, V. w., & zhou, Y. (2014). Kolmogorov–Smirnov test: Encima-
vista. Wiley Statsref: Statistics reference online. https://doi.org/10
.1002/9781118445112.stat06558
Bertschinger, NORTE., & Natschläger, t. (2004). Real-time computa-
tion at the edge of chaos in recurrent neural networks. Neural
Cálculo, 16(7), 1413–1436. https://doi.org/10.1162
/089976604323057443, PubMed: 15165396
Betzel, R. F., Avena-Koenigsberger, A., Goñi, J., Él, y., De Reus,
METRO. A., Griffa, A., … Sporns, oh. (2016). Generative models of the
human connectome. NeuroImagen, 124, 1054–1064. https://doi
.org/10.1016/j.neuroimage.2015.09.041, PubMed: 26427642
Betzel, R. F., & bassett, D. S. (2017a). Generative models for net-
work neuroscience: Prospects and promise. Journal of the Royal
Society Interface, 14(136), 20170623. https://doi.org/10.1098/rsif
.2017.0623, PubMed: 29187640
Betzel, R. F., & bassett, D. S. (2017b). Multi-scale brain networks.
NeuroImagen, 160, 73–83. https://doi.org/10.1016/j.neuroimage
.2016.11.006, PubMed: 27845257
Betzel, R. F., Medalla, j. D., & bassett, D. S. (2018). Diversity of
meso-scale architecture in human and non-human connectomes.
Comunicaciones de la naturaleza, 9(1), 1–14. https://doi.org/10.1038
/s41467-017-02681-z, PubMed: 29367627
Bhattacharjee, S. METRO., & Seno, F. (2001). A measure of data collapse
for scaling. Journal of Physics A: Mathematical and General,
34(33), 6375. https://doi.org/10.1088/0305-4470/34/33/302
Blythe, R. A., & McKane, A. j. (2007). Stochastic models of evolu-
tion in genetics, ecology and linguistics. Journal of Statistical
Mechanics: Theory and Experiment, 2007(7), P07018. https://
doi.org/10.1088/1742-5468/2007/07/P07018
Boedecker, J., Obst, o., Lizier, j. T., Mayer, norte. METRO., & Asada, METRO.
(2012). Information processing in echo state networks at the edge
of chaos. Theory in Biosciences, 131(3), 205–213. https://doi.org
/10.1007/s12064-011-0146-8, PubMed: 22147532
Bonachela, j. A., De Franciscis, S., Torres, j. J., & Muñoz, METRO. A.
(2010). Self-organization without conservation: Are neuronal
avalanches generically critical? Revista de mecánica estadística:
Theory and Experiment, 2010(2), P02015. https://doi.org/10
.1088/1742-5468/2010/02/P02015
Bonachela, j. A., & Muñoz, METRO. A. (2008). Confirming and extend-
ing the hypothesis of universality in sandpiles. Physical Review E,
78(4), 041102. https://doi.org/10.1103/ PhysRevE.78.041102,
PubMed: 18999374
Bonachela, j. A., & Muñoz, METRO. A. (2009). Self-organization without
conservation: True or just apparent scale-invariance? Diario de
Statistical Mechanics: Theory and Experiment, 2009(9), P09009.
https://doi.org/10.1088/1742-5468/2009/09/P09009
Breskin, I., Soriano, J., Moses, MI., & Tlusty, t. (2006). Percolation in
living neural networks. Physical Review Letters, 97(18), 188102.
https://doi.org/10.1103/ PhysRevLett.97.188102, PubMed:
17155581
Buendía, v., di Santo, S., Bonachela, j. A., & Muñoz, METRO. A. (2020a).
Feedback mechanisms for self-organization to the edge of a
phase transition. Frontiers in Physics, 8, 333. https://doi.org/10
.3389/fphy.2020.00333
Buendía, v., di Santo, S., Villegas, PAG., Burioni, r., & Muñoz, METRO. A.
(2020b). Self-organized bistability and its possible relevance for
brain dynamics. Physical Review Research, 2(1), 013318. https://
doi.org/10.1103/PhysRevResearch.2.013318
Buendía, v., Villegas, PAG., Burioni, r., & Muñoz, METRO. A. (2021).
Hybrid-type synchronization transitions: Where incipient oscilla-
ciones, scale-free avalanches, and bistability live together. Físico
Review Research, 3(2), 023224. https://doi.org/10.1103
/PhysRevResearch.3.023224
bullmore, MI., & despreciar, oh. (2012). The economy of brain network
organización. Naturaleza Reseñas Neurociencia, 13(5), 336–349.
https://doi.org/10.1038/nrn3214, PubMed: 22498897
Burns, B. D., & Webb, A. (1976). The spontaneous activity of
neurones in the cat’s cerebral cortex. Proceedings of the Royal
Society of London. Serie B. Ciencias Biologicas, 194(1115),
211–223. https://doi.org/10.1098/rspb.1976.0074, PubMed:
11486
Capolupo, A., Hombre libre, W.. J., & Vitiello, GRAMO. (2013). Dissipation of
‘dark energy’ by cortex in knowledge retrieval. Physics of Life
Reseñas, 10(1), 85–94. https://doi.org/10.1016/j.plrev.2013.01
.001, PubMed: 23333569
Carandini, METRO., & stevens, C. (2004). Amplification of trial-to-trial
response variability by neurons in visual cortex. Biología PLoS,
2(9), e264. https://doi.org/10.1371/journal.pbio.0020264,
PubMed: 15328535
Cardin, j. A. (2019). Functional flexibility in cortical circuits. Actual
Opinion in Neurobiology, 58, 175–180. https://doi.org/10.1016/j
.conb.2019.09.008, PubMed: 31585330
Carvalho, t. T., Fontenele, A. J., Girardi-Schappo, METRO., Feliciano, T.,
Aguiar, l. A., silva, t. PAG., … Copelli, METRO. (2021). Subsampled
directed-percolation models explain scaling relations experimen-
tally observed in the brain. Frontiers in Neural Circuits, 14,
576727. https://doi.org/10.3389/fncir.2020.576727, PubMed:
33519388
Chialvo, D. R. (2010). Emergent complex neural dynamics. Naturaleza
Physics, 6(10), 744–750. https://doi.org/10.1038/nphys1803
Chiappalone, METRO., Vato, A., Tedesco, METRO. B. Marcoli, METRO., Davide, F.,
& Martinoia, S. (2003). Networks of neurons coupled to micro-
electrode arrays: A neuronal sensory system for pharmacological
applications. Biosensors and Bioelectronics, 18(5–6), 627–634.
https://doi.org/10.1016/S0956-5663(03)00041-1, PubMed:
12706572
cláusula, A., shalizi, C. r., & Hombre nuevo, METRO. mi. (2009). Power-law dis-
tributions in empirical data. SIAM Review, 51(4), 661–703.
https://doi.org/10.1137/070710111
cocineros, l., Perdido, l. l., Brilla, A., & romper la lanza, METRO. (2017). Crit-
icality in the brain: A synthesis of neurobiology, models and cog-
nition. Progress in Neurobiology, 158, 132–152. https://doi.org
/10.1016/j.pneurobio.2017.07.002, PubMed: 28734836
cohen, METRO. r., & Kohn, A. (2011). Measuring and interpreting neu-
ronal correlations. Neurociencia de la naturaleza, 14(7), 811–819. https://
doi.org/10.1038/nn.2842, PubMed: 21709677
Cuello, GRAMO., & grifo, j. (2015). Brain activity and cognition: A
connection from thermodynamics and information theory. Fron-
tiers in Psychology, 6, 818. https://doi.org/10.3389/fpsyg.2015
.00818, PubMed: 26136709
Dalla Porta, l., & Copelli, METRO. (2019). Modeling neuronal ava-
lanches and long-range temporal correlations at the emergence
of collective oscillations: Continuously varying exponents mimic
m/eeg results. Biología Computacional PLoS, 15(4), e1006924.
Neurociencia en red
1179
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
https://doi.org/10.1371/journal.pcbi.1006924, PubMed:
30951525
David, o., & Friston, k. j. (2003). A neural mass model for MEG/
EEG: Coupling and neuronal dynamics. NeuroImagen, 20(3),
1743–1755. https://doi.org/10.1016/j.neuroimage.2003.07.015,
PubMed: 14642484
de Andrade Costa, A., Copelli, METRO., & Kinouchi, oh. (2015). Can
dynamical synapses produce true self-organized criticality?
Revista de mecánica estadística: Theory and Experiment,
2015(6), P06004. https://doi.org/10.1088/1742-5468/2015/06
/P06004
De Arcangelis, l., & Herrmann, h. j. (2010). Learning as a phenom-
enon occurring in a critical state. Actas del Nacional
Academia de Ciencias, 107(9), 3977–3981. https://doi.org/10
.1073/pnas.0912289107, PubMed: 20160107
De Arcangelis, l., Perrone-Capano, C., & Herrmann, h. j. (2006).
Self-organized criticality model for brain plasticity. Físico
Review Letters, 96(2), 028107. https://doi.org/10.1103
/PhysRevLett.96.028107, PubMed: 16486652
decoración, GRAMO., Tononi, GRAMO., boly, METRO., & Kringelbach, METRO. l. (2015).
Rethinking segregation and integration: Contributions of
whole-brain modelling. Naturaleza Reseñas Neurociencia, 16(7),
430–439. https://doi.org/10.1038/nrn3963, PubMed: 26081790
Del Papa, B., Priesemann, v., & Triesch, j. (2017). Criticality meets
aprendiendo: Criticality signatures in a self-organizing recurrent neu-
ral network. Más uno, 12(5), e0178683. https://doi.org/10.1371
/diario.pone.0178683, PubMed: 28552964
Del Pozo, S. METRO., Laufs, h., Bonhomme, v., Laureys, S., Balenzuela,
PAG., & Tagliazucchi, mi. (2021). Unconsciousness reconfigures
modular brain network dynamics. Chaos: An Interdisciplinary
Journal of Nonlinear Science, 31(9), 093117. https://doi.org/10
.1063/5.0046047, PubMed: 34598477
Deluca, A., & Corral, Á. (2013). Fitting and goodness-of-fit test of
non-truncated and truncated power-law distributions. Acta Geo-
physica, 61(6), 1351–1394. https://doi.org/10.2478/s11600-013
-0154-9
Dickman, r., Muñoz, METRO. A., Vespignani, A., & Zapperi, S. (2000). Paths
to self-organized criticality. Brazilian Journal of Physics, 30(1),
27–41. https://doi.org/10.1590/S0103-97332000000100004
Dickman, r., Vespignani, A., & Zapperi, S. (1998). Self-organized
criticality as an absorbing-state phase transition. Physical Review
mi, 57(5), 5095. https://doi.org/10.1103/PhysRevE.57.5095
di Santo, S., Burioni, r., Vezzani, A., & Muñoz, METRO. A. (2016). Self-
organized bistability associated with first-order phase transitions.
Physical Review Letters, 116(24), 240601. https://doi.org/10
.1103/PhysRevLett.116.240601, PubMed: 27367373
di Santo, S., Villegas, PAG., Burioni, r., & Muñoz, METRO. A. (2017). Simple
unified view of branching process statistics: Random walks in
balanced logarithmic potentials. Physical Review E, 95(3),
032115. https://doi.org/10.1103/PhysRevE.95.032115, PubMed:
28415350
di Santo, S., Villegas, PAG., Burioni, r., & Muñoz, METRO. A. (2018). Lan-
dau–Ginzburg theory of cortex dynamics: Scale-free avalanches
emerge at the edge of synchronization. Actas de la
Academia Nacional de Ciencias, 115(7), E1356–E1365. https://
doi.org/10.1073/pnas.1712989115, PubMed: 29378970
Ecker, A. S., Berens, PAG., Keliris, GRAMO. A., Bethge, METRO., Logothetis, norte. K.,
& Tolias, A. S. (2010). Decorrelated neuronal firing in cortical
microcircuits. Ciencia, 327(5965), 584–587. https://doi.org/10
.1126/science.1179867, PubMed: 20110506
Effenberger, F., Jost, J., & Levina, A. (2015). Self-organization in bal-
anced state networks by STDP and homeostatic plasticity. PLoS
Biología Computacional, 11(9), e1004420. https://doi.org/10.1371
/journal.pcbi.1004420, PubMed: 26335425
Erdos, l., Kruger, T., & Renfrew, D. (2018). Power law decay for
systems of randomly coupled differential equations. SIAM Journal
on Mathematical Analysis, 50(3), 3271–3290. https://doi.org/10
.1137/17M1143125
Fay, METRO. PAG., & Proschan, METRO. A. (2010). Wilcoxon-Mann-Whitney or
prueba t? On assumptions for hypothesis tests and multiple interpre-
tations of decision rules. Statistics Surveys, 4, 1–39. https://doi
.org/10.1214/09-SS051, PubMed: 20414472
Felleman, D. J., & VanEssen, D. C. (1991). Distributed hierarchical
processing in the primate cerebral cortex. Corteza cerebral, 1(1),
1–47. https://doi.org/10.1093/cercor/1.1.1-a, PubMed: 1822724
Fontenele, A. J., de Vasconcelos, norte. A., Feliciano, T., Aguiar, l. A.,
Soares-Cunha, C., Coimbra, B., … Copelli, METRO. (2019). Criticality
between cortical states. Physical Review Letters, 122(20),
208101. https://doi.org/10.1103/ PhysRevLett.122.208101,
PubMed: 31172737
Oscuro, l. J., Williams-García, R. v., ruega, j. METRO., & Ortíz, GRAMO.
(2021). Evidence for quasicritical brain dynamics. Físico
Review Letters, 126(9), 098101. https://doi.org/10.1103
/PhysRevLett.126.098101, PubMed: 33750159
Franke, F., Fiscella, METRO., Sevelev, METRO., Roska, B., Hierlemann, A., &
da Silveira, R. A. (2016). Structures of neural correlation and how
they favor coding. Neurona, 89(2), 409–422. https://doi.org/10
.1016/j.neuron.2015.12.037, PubMed: 26796692
Friedman, mi. J., & Landsberg, A. S. (2013). Hierarchical networks,
power laws, and neuronal avalanches. Chaos: An Interdisciplin-
ary Journal of Nonlinear Science, 23(1), 013135. https://doi.org
/10.1063/1.4793782, PubMed: 23556972
Friedman, NORTE., Ito, S., Brinkman, B. A., Shimono, METRO., DeVille, R. l.,
Dahmen, k. A., … Butler, t. C. (2012). Universal critical dynam-
ics in high resolution neuronal avalanche data. Physical Review
Letras, 108(20), 208102. https://doi.org/10.1103/ PhysRevLett
.108.208102, PubMed: 23003192
Fristedt, B. MI., & Gray, l. F. (2013). A modern approach to proba-
bility theory. Nueva York, Nueva York: Ciencia Springer & Medios comerciales.
Friston, k. (2009). The free-energy principle: A rough guide to the
cerebro? Tendencias en Ciencias Cognitivas, 13(7), 293–301. https://doi
.org/10.1016/j.tics.2009.04.005, PubMed: 19559644
Friston, k. (2010). The free-energy principle: A unified brain theory?
Naturaleza Reseñas Neurociencia, 11(2), 127–138. https://doi.org/10
.1038/nrn2787, PubMed: 20068583
Galván, A. (2010). Neural plasticity of development and learning.
Mapeo del cerebro humano, 31(6), 879–890. https://doi.org/10.1002
/hbm.21029, PubMed: 20496379
Ganmor, MI., Segev, r., & Schneidman, mi. (2011). Sparse low-order
interaction network underlies a highly correlated and learnable
neural population code. Proceedings of the National Academy
of Sciences, 108(23), 9679–9684. https://doi.org/10.1073/pnas
.1019641108, PubMed: 21602497
García-Pelayo, r., Salazar, I., & Schieve, W.. C. (1993). A branching
process model for sand avalanches. Journal of Statistical Physics,
72(1), 167–187. https://doi.org/10.1007/BF01048045
Neurociencia en red
1180
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Gautama, S. h., Hoang, t. T., McClanahan, K., grady, S. K., & mostrar,
W.. l. (2015). Maximizing sensory dynamic range by tuning the
cortical state to criticality. Biología Computacional PLoS, 11(12),
e1004576. https://doi.org/10.1371/journal.pcbi.1004576,
PubMed: 26623645
Gerstner, w., Kistler, W.. METRO., Beneficio, r., & paninski, l. (2014). nuevo-
ronal dynamics: From single neurons to networks and models of
cognition. Cambridge, Reino Unido: Prensa de la Universidad de Cambridge. https://
doi.org/10.1017/CBO9781107447615
Girardi-Schappo, METRO. (2021). Brain criticality beyond avalanches:
Open problems and how to approach them. Journal of Physics:
Complexity, 2(3), 031003. https://doi.org/10.1088/2632-072X
/ac2071
Girardi-Schappo, METRO., Bortolotto, GRAMO. S., Gonsalves, j. J., Pinto, l. T.,
& Tragtenberg, METRO. h. (2016). Griffiths phase and long-range cor-
relations in a biologically motivated visual cortex model. ciencia-
tific Reports, 6, 29561. https://doi.org/10.1038/srep29561,
PubMed: 27435679
Girardi-Schappo, METRO., Galera, mi. F., Carvalho, t. T., Brochini, l.,
Kamiji, norte. l., Roque, A. C., & Kinouchi, oh. (2021). A unified
theory of E/I synaptic balance, quasicritical neuronal avalanches
and asynchronous irregular spiking. Journal of Physics: Com-
plejidad, 2(4), 045001. https://doi.org/10.1088/2632-072X/ac2792
Gireesh, mi. D., & Plenz, D. (2008). Neuronal avalanches organize
as nested theta-and beta/gamma-oscillations during development
of cortical layer 2/3. Actas de la Academia Nacional de
Ciencias, 105(21), 7576–7581. https://doi.org/10.1073/pnas
.0800537105, PubMed: 18499802
Gnesotto, F. S., Mura, F., Gladrow, J., & Broedersz, C. PAG. (2018).
Broken detailed balance and non-equilibrium dynamics in
living systems: A review. Reports on Progress in Physics, 81(6),
066601. https://doi.org/10.1088/1361-6633/aab3ed, PubMed:
29504517
Gong, GRAMO., Él, y., Concha, l., Lebel, C., Bruto, D. w., evans, A. C.,
& Beaulieu, C. (2009). Mapping anatomical connectivity patterns
of human cerebral cortex using in vivo diffusion tensor imaging
tractography. Corteza cerebral, 19(3), 524–536. https://doi.org/10
.1093/cercor/bhn102, PubMed: 18567609
Gros, C. (2010). Complex and adaptive dynamical systems. Berlina,
Alemania: Saltador. https://doi.org/10.1007/978-3-642-04706-0
Guevara, R. (2021). Synchronization, free energy and the embryo-
genesis of the cortex. Physics of Life Reviews, 36, 5–6. https://doi
.org/10.1016/j.plrev.2020.11.006, PubMed: 33348117
guo, S., Chen, X., Liu, y., Kang, r., Liu, T., & li, D. (2021). Perco-
lation analysis of brain structural network. Frontiers in Physics, 9,
345. https://doi.org/10.3389/fphy.2021.698077
Hahn, GRAMO., Ponce-Alvarez, A., Monier, C., Benvenuti, GRAMO., Kumar, A.,
Chavane, F., … Frégnac, Y. (2017). Spontaneous cortical activity
is transiently poised close to criticality. PLoS computacional
Biología, 13(5), e1005543. https://doi.org/10.1371/journal.pcbi
.1005543, PubMed: 28542191
Haimovici, A., Tagliazucchi, MI., Balenzuela, PAG., & Chialvo, D. R.
(2013). Brain organization into resting state networks emerges
at criticality on a model of the human connectome. Físico
Review Letters, 110(17), 178101. https://doi.org/10.1103
/PhysRevLett.110.178101, PubMed: 23679783
Haldeman, C., & ruega, j. METRO. (2005). Critical branching captures
activity in living neural networks and maximizes the number
of metastable states. Physical Review Letters, 94(5), 058101.
https://doi.org/10.1103/ PhysRevLett.94.058101, PubMed:
15783702
Piedra dura, r., Mansvelder, h. D., & Linkenkaer-Hansen, k. (2014).
The neuronal network oscillation as a critical phenomenon. En
Criticality in neural systems (páginas. 293–316). Weinheim, Alemania:
wiley. https://doi.org/10.1002/9783527651009.ch13
harris, t. mi. (1963). The theory of branching processes (volumen. 6).
Berlina, Alemania: Saltador. https://doi.org/10.1007/978-3-642
-51866-9
Él, B. j. (2014). Scale-free brain activity: Past, present, and future.
Tendencias en Ciencias Cognitivas, 18(9), 480–487. https://doi.org/10
.1016/j.tics.2014.04.003, PubMed: 24788139
Henkel, METRO., Hinrichsen, h., Lübeck, S., & Pleimling, METRO. (2008).
Non-equilibrium phase transitions (volumen. 1). Berlina, Alemania:
Saltador.
Hernandez-Urbina, v., & Herrmann, j. METRO. (2017). Self-organized
criticality via retro-synaptic signals. Frontiers in Physics, 4, 54.
https://doi.org/10.3389/fphy.2016.00054
Hesse, J., & Bruto, t. (2014). Self-organized criticality as a funda-
mental property of neural systems. Frontiers in Systems Neurosci-
ence, 8, 166. https://doi.org/10.3389/fnsys.2014.00166,
PubMed: 25294989
Hinrichsen, h. (2000). Non-equilibrium critical phenomena and
phase transitions into absorbing states. Advances in Physics,
49(7), 815–958. https://doi.org/10.1080/00018730050198152
Holcman, D., & Tsodyks, METRO. (2006). The emergence of up and
down states in cortical networks. Biología Computacional PLoS,
2(3), e23. https://doi.org/10.1371/journal.pcbi.0020023,
PubMed: 16557293
Hopfield, j. j. (1982). Neural networks and physical systems with
emergent collective computational abilities. Actas de la
Academia Nacional de Ciencias, 79(8), 2554–2558. https://doi
.org/10.1073/pnas.79.8.2554, PubMed: 6953413
Janowsky, S. A., & Laberge, C. A. (1993). Exact solutions for a
mean-field abelian sandpile. Journal of Physics A: Matemático
and General, 26(19), L973. https://doi.org/10.1088/0305-4470
/26/19/001
Jercog, D., Roxin, A., Bartho, PAG., Luczak, A., Compte, A., & de la
Rocha, j. (2017). Up-down cortical dynamics reflect state transi-
tions in a bistable network. Elife, 6, e22425. https://doi.org/10
.7554/eLife.22425, PubMed: 28826485
Jung, NORTE., Le, q. A., Sotavento, K.-E., & Sotavento, j. W.. (2020). Avalanche size
distribution of an integrate-and-fire neural model on complex
redes. Chaos: An Interdisciplinary Journal of Nonlinear
Ciencia, 30(6), 063118. https://doi.org/10.1063/5.0008767,
PubMed: 32611110
Kaiser, METRO., & Hilgetag, C. C. (2006). Nonoptimal component place-
mento, but short processing paths, due to long-distance projec-
tions in neural systems. Biología Computacional PLoS, 2(7), e95.
https://doi.org/10.1371/journal.pcbi.0020095, PubMed:
16848638
Kaiser, METRO., & Hilgetag, C. C. (2010). Optimal hierarchical modular
topologies for producing limited sustained activation of neural
redes. Frontiers in Neuroinformatics, 4, 8. https://doi.org/10
.3389/fninf.2010.00008, PubMed: 20514144
Kanji, GRAMO. k. (2006). 100 statistical tests. Londres, Reino Unido: Sage. https://
doi.org/10.4135/9781849208499
Neurociencia en red
1181
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
.
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Kara, PAG., Reinagel, PAG., & Reid, R. C. (2000). Low response variability
in simultaneously recorded retinal, thalamic, and cortical neu-
ron. Neurona, 27(3), 635–646. https://doi.org/10.1016/S0896
-6273(00)00072-6, PubMed: 11055444
Katsnelson, METRO. I., Vanchurin, v., & Westerhout, t. (2021). Self-
organized criticality in neural networks. arXiv preprint
arXiv:2107.03402. https://doi.org/10.48550/arXiv.2107.03402
Keogh, mi. J., & Pazzani, METRO. j. (2001). Derivative dynamic time
warping. En Actas de la 2001 SIAM international confer-
ence on data mining (páginas. 1–11). https://doi.org/10.1137/1
.9781611972719.1
Khambhati, A. NORTE., Sizemore, A. MI., Betzel, R. F., & bassett, D. S.
(2018). Modeling and interpreting mesoscale network dynamics.
NeuroImagen, 180, 337–349. https://doi.org/10.1016/j
.neuroimage.2017.06.029, PubMed: 28645844
Kinouchi, o., & Copelli, METRO. (2006). Optimal dynamical range of
excitable networks at criticality. Nature Physics, 2(5), 348–351.
https://doi.org/10.1038/nphys289
Krotov, D., & Hopfield, j. (2020). Large associative memory prob-
lem in neurobiology and machine learning. arXiv preprint
arXiv:2008.06996. https://doi.org/10.48550/arXiv.2008.06996
Larremore, D. B., Carpintero, METRO. y., Ott, MI., & Restrepo, j. GRAMO. (2012).
Statistical properties of avalanches in networks. Physical Review
mi, 85(6), 066131. https://doi.org/10.1103/PhysRevE.85.066131,
PubMed: 23005186
Laurson, l., & Alava, METRO. j. (2006). 1/f noise and avalanche scaling
in plastic deformation. Physical Review E, 74(6), 066106. https://
doi.org/10.1103/PhysRevE.74.066106, PubMed: 17280120
Laurson, l., Illa, X., & Alava, METRO. j. (2009). The effect of thresholding
on temporal avalanche statistics. Revista de mecánica estadística:
Theory and Experiment, 2009(1), P01019. https://doi.org/10
.1088/1742-5468/2009/01/P01019
Laurson, l., Illa, X., Santucci, S., Tore Tallakstad, K., Måløy, k. J., &
Alava, METRO. j. (2013). Evolution of the average avalanche shape
with the universality class. Comunicaciones de la naturaleza, 4(1), 1–6.
https://doi.org/10.1038/ncomms3927, PubMed: 24352571
Sotavento, D. S., Goh, k. I., Kahng, B., & kim, D. (2004). Branching pro-
cess approach to avalanche dynamics on complex networks.
Journal of the Korean Physical Society, 44(1), 633–637. https://
doi.org/10.3938/jkps.44.633
Sotavento, h., Wang, S., & Hudetz, A. GRAMO. (2020). State-dependent
cortical unit activity reflects dynamic brain state transitions in
Revista de neurociencia, 40(49), 9440–9454.
anesthesia.
https://doi.org/10.1523/ JNEUROSCI.0601-20.2020, PubMed:
33122389
Levina, A., Herrmann, j. METRO., & Geisel, t. (2007). Dynamical synap-
ses causing self-organized criticality in neural networks. Naturaleza
Physics, 3(12), 857–860. https://doi.org/10.1038/nphys758
Levina, A., Herrmann, j. METRO., & Geisel, t. (2009). Phase transitions
towards criticality in a neural system with adaptive interactions.
Physical Review Letters, 102(11), 118110. https://doi.org/10.1103
/PhysRevLett.102.118110, PubMed: 19392248
li, J., & mostrar, W.. l. (2020). Tuning network dynamics from critical-
ity to an asynchronous state. Biología Computacional PLoS, 16(9),
e1008268. https://doi.org/10.1371/journal.pcbi.1008268,
PubMed: 32986705
li, X., & Pequeño, METRO. (2012). Neuronal avalanches of a self-organized
neural network with active-neuron-dominant structure. Chaos:
An Interdisciplinary Journal of Nonlinear Science, 22(2), 023104.
https://doi.org/10.1063/1.3701946, PubMed: 22757511
Liggett, t. METRO. (2006). Interacting particle systems. Nueva York, Nueva York:
Ciencia Springer & Medios comerciales.
Linkenkaer-Hansen, K., Nikouline, V. v., Palva, j. METRO., & Ilmoniemi,
R. j. (2001). Long-range temporal correlations and scaling behav-
ior in human brain oscillations. Revista de neurociencia, 21(4),
1370–1377. https://doi.org/10.1523/ JNEUROSCI.21-04-01370
.2001, PubMed: 11160408
Lombardi, F., Herrmann, h. J., & De Arcangelis, l. (2017). Balance
of excitation and inhibition determines 1/f power spectrum in
neuronal networks. Chaos: An Interdisciplinary Journal of Non-
linear Science, 27(4), 047402. https://doi.org/10.1063/1
.4979043, PubMed: 28456161
Lübeck, S. (2004). Universal scaling behavior of non-equilibrium
phase transitions. International Journal of Modern Physics B,
1 8 ( 3 1 norte 3 2 ) , 3 9 7 7 – 4 1 1 8 . h t t p s : / / d o i . o r g / 1 0 . 1 1 4 2
/S0217979204027748
Lübeck, S., & Heger, PAG. (2003). Universal finite-size scaling behav-
ior and universal dynamical scaling behavior of absorbing phase
transitions with a conserved field. Physical Review E, 68(5),
056102. https://doi.org/10.1103/PhysRevE.68.056102, PubMed:
14682841
Lynn, C. w., & bassett, D. S. (2019). The physics of brain network
estructura, function and control. Nature Reviews Physics, 1(5),
318–332. https://doi.org/10.1038/s42254-019-0040-8
Lynn, C. w., Cornblath, mi. J., Papadopoulos, l., Bertolero, METRO. A., &
bassett, D. S. (2021). Broken detailed balance and entropy pro-
duction in the human brain. Proceedings of the National Acad-
emy of Sciences, 118(47), e2109889118. https://doi.org/10.1073
/pnas.2109889118, PubMed: 34789565
Mamá, Z., Turrigiano, GRAMO. GRAMO., wessel, r., & Hengen, k. B. (2019).
Cortical circuit dynamics are homeostatically tuned to criticality
in vivo. Neurona, 104(4), 655–664. https://doi.org/10.1016/j
.neuron.2019.08.031, PubMed: 31601510
Malcai, o., Shilo, y., & Biham, oh. (2006). Dissipative sandpile
models with universal exponents. Physical Review E, 73(5),
056125. https://doi.org/10.1103/PhysRevE.73.056125, PubMed:
16803016
Manna, S. S. (1991). Two-state model of self-organized criticality.
Journal of Physics A: Mathematical and General, 24(7), L363.
https://doi.org/10.1088/0305-4470/24/7/009
Mariani, B., Nicoletti, GRAMO., Bisio, METRO., Maschietto, METRO., Oboe, r.,
Suweis, S., & Vassanelli, S. (2021). Beyond resting state neuronal
avalanches in the somatosensory barrel cortex. bioRxiv. https://
doi.org/10.1101/2021.05.13.444047
Marković, D., & Gros, C. (2014). Power laws and self-organized
criticality in theory and nature. Physics Reports, 536(2), 41–74.
https://doi.org/10.1016/j.physrep.2013.11.002
marshall, NORTE., Hora, norte. METRO., bennett, NORTE., Ripp, METRO., Lautzenhiser, MI.,
& ruega, j. METRO. (2016). Analysis of power laws, shape collapses,
and neural complexity: New techniques and MATLAB support
via the NCC toolbox. Frontiers in Physiology, 7, 250. https://doi
.org/10.3389/fphys.2016.00250, PubMed: 27445842
Martinello, METRO., Hidalgo, J., Maritan, A., di Santo, S., Plenz, D., &
Muñoz, METRO. A. (2017). Neutral theory and scale-free neural
dinámica. Physical Review X, 7(4), 041071. https://doi.org/10
.1103/PhysRevX.7.041071
Neurociencia en red
1182
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
/
t
.
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Mehta, A. PAG., Mills, A. C., Dahmen, k. A., & Sethna, j. PAG. (2002).
Universal pulse shape scaling function and exponents: Critical
test for avalanche models applied to Barkhausen noise. Físico
Review E, 65(4), 046139. https://doi.org/10.1103/PhysRevE.65
.046139, PubMed: 12005958
Meisel, C., & Bruto, t. (2009). Adaptive self-organization in a real-
istic neural network model. Physical Review E, 80(6), 061917.
https://doi.org/10.1103/ PhysRevE.80.061917, PubMed:
20365200
Meisel, C., Olbrich, MI., Shriki, o., & Achermann, PAG. (2013). Fading
signatures of critical brain dynamics during sustained wakefulness
in humans. Revista de neurociencia, 33(44), 17363–17372.
https://doi.org/10.1523/ JNEUROSCI.1516-13.2013, PubMed:
24174669
Molinero, PAG., & Wang, X.-J. (2006). Power-law neuronal fluctuations in
a recurrent network model of parametric working memory. Jour-
nal of Neurophysiology, 95(2), 1099–1114. https://doi.org/10
.1152/jn.00491.2005, PubMed: 16236788
millman, D., Mihalas, S., kirkwood, A., & Niebur, mi. (2010). Self-
organized criticality occurs in non-conservative neuronal net-
works during ‘up’ states. Nature Physics, 6(10), 801–805.
https://doi.org/10.1038/nphys1757, PubMed: 21804861
Mitzenmacher, METRO. (2004). A brief history of generative models for
power law and lognormal distributions. Internet Mathematics, 1(2),
226–251. https://doi.org/10.1080/15427951.2004.10129088
Molgedey, l., Schuchhardt, J., & Schuster, h. GRAMO. (1992). Suppress-
ing chaos in neural networks by noise. Physical Review Letters,
69(26), 3717–3719. https://doi.org/10.1103/ PhysRevLett.69
.3717, PubMed: 10046895
montesco, PAG. r., Dayán, PAG., & Sejnowski, t. j. (1996). A framework
for mesencephalic dopamine systems based on predictive
Hebbian learning. Revista de neurociencia, 16(5), 1936–1947.
https://doi.org/10.1523/ JNEUROSCI.16-05-01936.1996,
PubMed: 8774460
Moretti, PAG., & Muñoz, METRO. A. (2013). Griffiths phases and the
stretching of criticality in brain networks. Nature Communica-
ciones, 4(1), 1–10. https://doi.org/10.1038/ncomms3521,
PubMed: 24088740
Muñoz, METRO. A. (2018). Coloquio: Criticality and dynamical scal-
ing in living systems. Reseñas de Física Moderna, 90(3), 031001.
https://doi.org/10.1103/RevModPhys.90.031001
Narayan, o., & Middleton, A. A. (1994). Avalanches and the renor-
malization group for pinned charge-density waves. Físico
Review B, 49(1), 244–256. https://doi.org/10.1103/PhysRevB.49
.244, PubMed: 10009280
Otter, R. (1949). The multiplicative process. The Annals of Mathe-
matical Statistics, 20(2), 206–224. https://doi.org/10.1214/aoms
/1177730031
Palva, j. METRO., Zhigalov, A., Hirvonen, J., Korhonen, o., Linkenkaer-
Hansen, K., & Palva, S. (2013). Neuronal long-range temporal
correlations and avalanche dynamics are correlated with
behavioral scaling laws. Proceedings of the National Academy
of Sciences, 110(9), 3585–3590. https://doi.org/10.1073/pnas
.1216855110, PubMed: 23401536
Papanikolaou, S., Bohn, F., verano, R. l., Durin, GRAMO., Zapperi, S., &
Sethna, j. PAG. (2011). Universality beyond power laws and the
average avalanche shape. Nature Physics, 7(4), 316–320.
https://doi.org/10.1038/nphys1884
Pausch, J., Garcia-Millan, r., & Pruessner, GRAMO. (2020). Time-
dependent branching processes: A model of oscillating neuronal
avalanches. Informes Científicos, 10(1), 13678. https://doi.org/10
.1038/s41598-020-69705-5, PubMed: 32792658
Perl, Y. S., Bocaccio, h., Pallavicini, C., Pérez-Ipiña, I., Laureys, S.,
Laufs, h., … Tagliazucchi, mi. (2021). Nonequilibrium brain
dynamics as a signature of consciousness. Physical Review E,
104(1), 014411. https://doi.org/10.1103/PhysRevE.104.014411,
PubMed: 34412335.
Petermann, T., Thiagarajan, t. C., Lebedev, METRO. A., Nicolelis, METRO. A.,
Chialvo, D. r., & Plenz, D. (2009). Spontaneous cortical activity in
awake monkeys composed of neuronal avalanches. Actas
of the National Academy of Sciences, 106(37), 15921–15926.
https://doi.org/10.1073/pnas.0904089106, PubMed: 19717463
Cabello, S.-S., Piedra dura, r., Mansvelder, h. D., & Linkenkaer-
Hansen, k. (2012). Critical-state dynamics of avalanches and
oscillations jointly emerge from balanced excitation/inhibition in
neuronal networks. Revista de neurociencia, 32(29), 9817–9823.
https://doi.org/10.1523/ JNEUROSCI.5990-11.2012, PubMed:
22815496
Cabello, S.-S., van Ooyen, A., & Linkenkaer-Hansen, k. (2008). Ava-
lanche dynamics of human brain oscillations: Relation to critical
branching processes and temporal correlations. Human Brain
Cartografía, 29(7), 770–777. https://doi.org/10.1002/hbm.20590,
PubMed: 18454457
Ponce-Alvarez, A., Jouary, A., Privat, METRO., decoración, GRAMO., & Sumbre, GRAMO.
(2018). Whole-brain neuronal activity displays crackling noise
dinámica. Neurona, 100 (6), 1446–1459. https://doi.org/10.1016
/j.neuron.2018.10.045, PubMed: 30449656
Rao, METRO. METRO., & Swift, R. j. (2006). Probability theory with applica-
ciones (volumen. 582). Nueva York, Nueva York: Ciencia Springer & Negocio
Media.
Reed, W.. J., & abrazos, B. D. (2002). From gene families and genera
to incomes and internet file sizes: Why power laws are so com-
mon in nature. Physical Review E, 66(6), 067103. https://doi.org
/10.1103/PhysRevE.66.067103, PubMed: 12513446
Reimer, J., Froudarakis, MI., Cadwell, C. r., Yatsenko, D., Denfield,
GRAMO. h., & Tolias, A. S. (2014). Pupil fluctuations track fast switch-
ing of cortical states during quiet wakefulness. Neurona, 84(2),
355–362. https://doi.org/10.1016/j.neuron.2014.09.033,
PubMed: 25374359
Reiss, R.-D., & tomás, METRO. (2007). Statistical analysis of extreme
valores: With applications to insurance, finance, hydrology and
other fields. Nueva York, Nueva York: Ciencia Springer & Medios comerciales.
Ribeiro, t. l., Copelli, METRO., Caixeta, F., Belchior, h., Chialvo, D. r.,
Nicolelis, METRO. A., & Ribeiro, S. (2010). Spike avalanches exhibit
universal dynamics across the sleep-wake cycle. Más uno,
5(11), e14129. https://doi.org/10.1371/journal.pone.0014129,
PubMed: 21152422
Roberto, PAG., & Vignoud, GRAMO. (2021). Stochastic models of neural
synaptic plasticity: A scaling approach. SIAM Journal on Applied
Matemáticas, 81(6), 2362–2386. https://doi.org/10.1137
/20M1382891
robinson, PAG. (2021). Neural field theory of neural avalanche expo-
nents. Cibernética biológica, 115(3), 237–243. https://doi.org/10
.1007/s00422-021-00875-9, PubMed: 33939016
Rubinov, METRO., despreciar, o., Thivierge, J.-P., & romper la lanza, METRO. (2011).
Neurobiologically realistic determinants of self-organized
Neurociencia en red
1183
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
t
/
/
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
t
.
/
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
criticality in networks of spiking neurons. PLoS computacional
Biología, 7(6), e1002038. https://doi.org/10.1371/journal.pcbi
.1002038, PubMed: 21673863
Sartori, PAG., Granjero, l., Sotavento, C. F., & Horowitz, j. METRO. (2014).
Thermodynamic costs of information processing in sensory adap-
tation. Biología Computacional PLoS, 10(12), e1003974. https://
doi.org/10.1371/journal.pcbi.1003974, PubMed: 25503948
Schaworonków, NORTE., Blythe, D. A., Kegeles, J., Curio, GRAMO., & Nikulin,
V. V. (2015). Power-law dynamics in neuronal and behavioral
data introduce spurious correlations. Mapeo del cerebro humano,
36(8), 2901–2914. https://doi.org/10.1002/ hbm.22816,
PubMed: 25930148
Schneidman, MI., Berry, METRO. J., Segev, r., & Bialek, W.. (2006). Weak
pairwise correlations imply strongly correlated network states in
a neural population. Naturaleza, 440(7087), 1007–1012. https://doi
.org/10.1038/nature04701, PubMed: 16625187
Schwarz, GRAMO. (1978). Estimating the dimension of a model. El
Annals of Statistics, 6(2), 461–464. https://doi.org/10.1214/aos
/1176344136
Scott, A. (1977). Neurophysics. Nueva York, Nueva York: John Wiley & Sons.
Scott, GRAMO., Fagerholm, mi. D., Mutoh, h., Sanguijuela, r., Sharp, D. J.,
mostrar, W.. l., & Knöpfel, t. (2014). Voltage imaging of waking
mouse cortex reveals emergence of critical neuronal dynamics.
Revista de neurociencia, 34(50), 16611–16620. https://doi.org
/10.1523/JNEUROSCI.3474-14.2014, PubMed: 25505314
Senzai, y., Fernandez-Ruiz, A., & Buzsáki, GRAMO. (2019). Layer-specific
physiological features and interlaminar interactions in the
primary visual cortex of the mouse. Neurona, 101(3), 500–513.
https://doi.org/10.1016/j.neuron.2018.12.009, PubMed:
30635232
Sethna, j. PAG., Dahmen, k. A., & miers, C. R. (2001). Crackling
ruido. Naturaleza, 410(6825), 242–250. https://doi.org/10.1038
/35065675, PubMed: 11258379
Shaukat, A., & Thivierge, J.-P. (2016). Statistical evaluation of
waveform collapse reveals scale-free properties of neuronal
avalanches. Frontiers in Computational Neuroscience, 10, 29.
https://doi.org/10.3389/fncom.2016.00029, PubMed: 27092071
mostrar, W.. l., Clawson, W.. PAG., Pobst, J., Karimipanah, y., Wright,
norte. C., & wessel, R. (2015). Adaptation to sensory input tunes
visual cortex to criticality. Nature Physics, 11(8), 659–663.
https://doi.org/10.1038/nphys3370
mostrar, W.. l., & Plenz, D. (2013). The functional benefits of critical-
ity in the cortex. The Neuroscientist, 19(1), 88–100. https://doi
.org/10.1177/1073858412445487, PubMed: 22627091
mostrar, W.. l., Cual, h., Petermann, T., roy, r., & Plenz, D.
(2009). Neuronal avalanches imply maximum dynamic range in
cortical networks at criticality. Revista de neurociencia, 29(49),
15595–15600. https://doi.org/10.1523/ JNEUROSCI.3864-09
.2009, PubMed: 20007483
mostrar, W.. l., Cual, h., Yu, S., roy, r., & Plenz, D. (2011). Informa-
tion capacity and transmission are maximized in balanced corti-
cal networks with neuronal avalanches. Revista de neurociencia,
31(1), 55–63. https://doi.org/10.1523/ JNEUROSCI.4637-10
.2011, PubMed: 21209189
espinilla, C.-W., & kim, S. (2006). Self-organized criticality and
scale-free properties in emergent functional neural networks.
Physical Review E, 74(4), 045101. https://doi.org/10.1103
/PhysRevE.74.045101, PubMed: 17155118
Shriki, o., Alstott, J., Carver, F., Holroyd, T., Henson, R. NORTE., Herrero,
METRO. l., … Plenz, D. (2013). Neuronal avalanches in the resting
MEG of the human brain. Revista de neurociencia, 33(16),
7079–7090. https://doi.org/10.1523/JNEUROSCI.4286-12.2013,
PubMed: 23595765
Smit, D. J., de Geus, mi. J., van de Nieuwenhuijzen, METRO. MI., camioneta
Beijsterveldt, C. MI., van Baal, GRAMO. C. METRO., Mansvelder, h. D., …
Linkenkaer-Hansen, k. (2011). Scale-free modulation of
resting-state neuronal oscillations reflects prolonged brain
maturation in humans. Revista de neurociencia, 31(37),
13128–13136. https://doi.org/10.1523/ JNEUROSCI.1678-11
.2011, PubMed: 21917796
Softky, W.. r., & Koch, C. (1993). The highly irregular firing of cor-
tical cells is inconsistent with temporal integration of random
EPSPs. Revista de neurociencia, 13(1), 334–350. https://doi.org
/10.1523/JNEUROSCI.13-01-00334.1993, PubMed: 8423479
Song, S., Molinero, k. D., & Abbott, l. F. (2000). Competitive Hebbian
learning through spike-timing-dependent synaptic plasticity.
Neurociencia de la naturaleza, 3(9), 919–926. https://doi.org/10.1038
/78829, PubMed: 10966623
Sornette, D., Johansen, A., & Dornic, I. (1995). Cartografía
self-organized criticality onto criticality. Journal de Physique I,
5(3), 325–335. https://doi.org/10.1051/jp1:1995129
Squire, l., Iceberg, D., Bloom, F. MI., Du Lac, S., Ghosh, A., & Spitzer,
norte. C. (2012). Fundamental neuroscience. Cambridge, MAMÁ:
Prensa académica.
Steif, j. mi. (2009). A survey of dynamical percolation. In Fractal
geometry and stochastics IV (páginas. 145–174). Bostón, MAMÁ: Saltador.
https://doi.org/10.1007/978-3-0346-0030-9_5
piedra, R. B., Gossen, mi. r., & jones, k. mi. (2005). Neuronal vari-
capacidad: Noise or part of the signal? Naturaleza Reseñas Neurociencia,
6(5), 389–397. https://doi.org/10.1038/nrn1668, PubMed:
15861181
Stepp, NORTE., Plenz, D., & Srinivasa, norte. (2015). Synaptic plasticity
enables adaptive self-tuning critical networks. PLoS Computa-
tional Biology, 11(1), e1004043. https://doi.org/10.1371/journal
.pcbi.1004043, PubMed: 25590427
Stewart, C. v., & Plenz, D. (2006). Inverted-U profile of dopamine–
NMDA-mediated spontaneous avalanche recurrence in superfi-
cial layers of rat prefrontal cortex. Revista de neurociencia,
26(31), 8148–8159. https://doi.org/10.1523/JNEUROSCI.0723
-06.2006, PubMed: 16885228
Stewart, C. v., & Plenz, D. (2008). Homeostasis of neuronal ava-
lanches during postnatal cortex development in vitro. Diario
of Neuroscience Methods, 169(2), 405–416. https://doi.org/10
.1016/j.jneumeth.2007.10.021, PubMed: 18082894
Tagliazucchi, MI., Balenzuela, PAG., Fraiman, D., & Chialvo, D. R.
(2012). Criticality in large-scale brain fMRI dynamics unveiled
by a novel point process analysis. Frontiers in Physiology, 3, 15.
https://doi.org/10.3389/fphys.2012.00015, PubMed: 22347863
tian, y., li, GRAMO., & Sol, PAG. (2021). Bridging the information and
dynamics attributes of neural activities. Physical Review Research,
3(4), 043085. https://doi.org/10.1103/PhysRevResearch.3.043085
tian, y., & Sol, PAG. (2021). Characteristics of the neural coding of
causality. Physical Review E, 103(1), 012406. https://doi.org/10
.1103/PhysRevE.103.012406, PubMed: 33601638
tian, y., & Sol, PAG. (2022). Information thermodynamics of encoding
and encoders. Chaos: An Interdisciplinary Journal of Nonlinear
Neurociencia en red
1184
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
.
/
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
Theoretical foundations of studying criticality in the brain
Ciencia, 32(6), 063109. https://doi.org/10.1063/5.0068115,
PubMed: 35778156
Tkačik, GRAMO., Mora, T., Marre, o., Amodei, D., Palmer, S. MI., Berry,
METRO. J., & Bialek, W.. (2015). Thermodynamics and signatures of
criticality in a network of neurons. Actas del Nacional
Academia de Ciencias, 112(37), 11508–11513. https://doi.org/10
.1073/pnas.1514188112, PubMed: 26330611
Touboul, J., & Destexhe, A. (2010). Can power-law scaling and
neuronal avalanches arise from stochastic dynamics? PLoS
Uno, 5(2), e8982. https://doi.org/10.1371/journal.pone
.0008982, PubMed: 20161798
Touboul, J., & Destexhe, A. (2017). Power-law statistics and univer-
sal scaling in the absence of criticality. Physical Review E, 95(1),
012413. https://doi.org/10.1103/PhysRevE.95.012413, PubMed:
28208383
Touboul, J., Wendling, F., Chauvel, PAG., & Faugeras, oh. (2011).
Neural mass activity, bifurcations, and epilepsy. Neural Compu-
tation, 23(12), 3232–3286. https://doi.org/10.1162/ NECO_a
_00206, PubMed: 21919787.
Tyulmankov, D., Fang, C., Vadaparty, A., & Cual, GRAMO. R. (2021).
Biological key-value memory networks. Avances en Neurología
Sistemas de procesamiento de información, 34.
Van Hateren, j. h., & van der Schaaf, A. (1998). Independent com-
ponent filters of natural images compared with simple cells in
primary visual cortex. Proceedings of the Royal Society of
Londres. Serie B: Ciencias Biologicas, 265(1394), 359–366.
https://doi.org/10.1098/rspb.1998.0303, PubMed: 9523437
Van Vreeswijk, C., & Sompolinsky, h. (1996). Chaos in neuronal
networks with balanced excitatory and inhibitory activity.
Ciencia, 274(5293), 1724–1726. https://doi.org/10.1126
/science.274.5293.1724, PubMed: 8939866
Varley, t. F., despreciar, o., Chip, A., & ruega, j. (2020). Differential
effects of propofol and ketamine on critical brain dynamics. PLoS
Biología Computacional, 16(12), e1008418. https://doi.org/10
.1371/journal.pcbi.1008418, PubMed: 33347455.
Villegas, PAG., di Santo, S., Burioni, r., & Muñoz, METRO. A. (2019). Time-
series thresholding and the definition of avalanche size. Físico
Review E, 100(1), 012133. https://doi.org/10.1103/PhysRevE.100
.012133, PubMed: 31499802
Villegas, PAG., Moretti, PAG., & Muñoz, METRO. A. (2014). Frustrated hierar-
chical synchronization and emergent complexity in the human
connectome network. Informes Científicos, 4, 5990. https://doi.org
/10.1038/srep05990, PubMed: 25103684
Obras, y., & cláusula, A. (2014). Power-law distributions in binned
empirical data. The Annals of Applied Statistics, 8(1), 89–119.
https://doi.org/10.1214/13-AOAS710
Vuong, q. h. (1989). Likelihood ratio tests for model selection and
non-nested hypotheses. Econometrica: Journal of the Economet-
ric Society, 57(2), 307–333. https://doi.org/10.2307/1912557
Wang, S., & zhou, C. (2012). Hierarchical modular structure
enhances the robustness of self-organized criticality in neural
redes. New Journal of Physics, 14(2), 023005. https://doi
.org/10.1088/1367-2630/14/2/023005
Wang, X. r., Lizier, j. T., & Prokopenko, METRO. (2011). Fisher informa-
tion at the edge of chaos in random boolean networks. Artificial
Life, 17(4), 315–329. https://doi.org/10.1162/artl_a_00041,
PubMed: 21762019
Williams-García, R. v., moore, METRO., ruega, j. METRO., & Ortíz, GRAMO. (2014).
Quasicritical brain dynamics on a nonequilibrium Widom line.
Physical Review E, 90(6), 062714. https://doi.org/10.1103
/PhysRevE.90.062714, PubMed: 25615136
Wilting, J., & Priesemann, V. (2019). Between perfectly critical and
fully irregular: A reverberating model captures and predicts cor-
tical spike propagation. Corteza cerebral, 29(6), 2759–2770.
https://doi.org/10.1093/cercor/bhz049, PubMed: 31008508
Lobo, F. (2005). Symmetry, multistability, and long-range interac-
tions in brain development. Physical Review Letters, 95(20),
208701. https://doi.org/10.1103/ PhysRevLett.95.208701,
PubMed: 16384113
Yaghoubi, METRO., de Graaf, T., Orlandi, j. GRAMO., Girotto, F., Colicos,
METRO. A., & Davidsen, j. (2018). Neuronal avalanche dynamics indi-
cates different universality classes in neuronal cultures. Científico
Informes, 8(1), 3417. https://doi.org/10.1038/s41598-018-21730-1,
PubMed: 29467426
Cual, GRAMO., & Schoenholz, S. (2017). Mean field residual networks:
On the edge of chaos. Avances en el procesamiento de información neuronal
Sistemas, 30.
Yu, S., Klaus, A., Cual, h., & Plenz, D. (2014). Scale-invariant neu-
ronal avalanche dynamics and the cut-off in size distributions.
Más uno, 9(6), e99761. https://doi.org/10.1371/journal.pone
.0099761, PubMed: 24927158
Zhigalov, A., Arnulfo, GRAMO., Nobili, l., Palva, S., & Palva, j. METRO. (2015).
Relationship of fast- and slow-timescale neuronal dynamics in
human MEG and SEEG. Revista de neurociencia, 35(13),
5385–5396. https://doi.org/10.1523/JNEUROSCI.4880-14.2015,
PubMed: 25834062
Neurociencia en red
1185
yo
D
oh
w
norte
oh
a
d
mi
d
F
r
oh
metro
h
t
t
pag
:
/
/
d
i
r
mi
C
t
.
metro
i
t
.
/
/
t
mi
d
tu
norte
mi
norte
a
r
t
i
C
mi
–
pag
d
yo
F
/
/
/
/
6
4
1
1
4
8
2
0
5
9
7
7
7
norte
mi
norte
_
a
_
0
0
2
6
9
pag
d
.
/
t
F
b
y
gramo
tu
mi
s
t
t
oh
norte
0
7
S
mi
pag
mi
metro
b
mi
r
2
0
2
3
