Metastable brain waves

Shaping dynamical neural computations using spatiotemporal constraints

Dynamics play a critical role in computation. The principled evolution of states over time enables both biological and artificial networks to represent and integrate information to make decisions. In the past few decades, significant multidisciplinary progress has been made in bridging the gap between how we understand biological versus artificial computation, including how insights gained from one can translate to the other. Research has revealed that neurobiology is a key determinant of brain network architecture, which gives rise to spatiotemporally constrained patterns of activity that underlie computation. Here, we discuss how neural systems use dynamics for computation, and claim that the biological constraints that shape brain networks may be leveraged to improve the implementation of artificial neural networks. To formalize this discussion, we consider a natural artificial analog of the brain that has been used extensively to model neural computation: the recurrent neural network (RNN). In both the brain and the RNN, we emphasize the common computational substrate atop which dynamics occur-the connectivity between neurons-and we explore the unique computational advantages offered by biophysical constraints such as resource efficiency, spatial embedding, and neurodevelopment.

NHSMM-MAR-sdNC: A novel data-driven computational framework for state-dependent effective connectivity analysis

The brain exhibits intrinsic dynamics characterized by spontaneous spatiotemporal reorganization of neural activity or metastability, which is associated closely with functional integration and segregation. Compared to dynamic functional connectivity, state-dependent effective connectivity (i.e., dynamic effective connectivity) is more suitable for exploring the metastability as its ability to infer causalities between brain regions. However, methods for state-dependent effective connectivity are scarce and urgently needed. In this study, a novel data-driven computational framework, named NHSMM-MAR-sdNC integrating nonparametric hidden semi-Markov model combined with multivariate autoregressive model and state-dependent new causality, is proposed to investigate the state-dependent effective connectivity. The framework is not constrained by any biological assumptions. Furthermore, state number can be inferred from the observed data directly and the state duration distributions will be estimated explicitly rather than restricted by geometric form, which overcomes limitations of hidden Markov model. Experimental results of synthetic data show that the framework can identify the state number adaptively and the state-dependent causality networks accurately. The dynamics of state-related causality networks are also revealed by the new method on real-world resting-state fMRI data. Our method provides a new data-driven computational framework for identifying state-dependent effective connectivity, which will facilitate the identification and assessment of metastability and itinerant dynamics of the brain.

A network control theory pipeline for studying the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

Discovering the effective connectome of the brain with dynamic Bayesian DAG learning

Understanding the complex mechanisms of the brain can be unraveled by extracting the Dynamic Effective Connectome (DEC). Recently, score-based Directed Acyclic Graph (DAG) discovery methods have shown significant improvements in extracting the causal structure and inferring effective connectivity. However, learning DEC through these methods still faces two main challenges: one with the fundamental impotence of high-dimensional dynamic DAG discovery methods and the other with the low quality of fMRI data. In this paper, we introduce Bayesian Dynamic DAG learning with M-matrices Acyclicity characterization (BDyMA) method to address the challenges in discovering DEC. The presented dynamic DAG enables us to discover direct feedback loop edges as well. Leveraging an unconstrained framework in the BDyMA method leads to more accurate results in detecting high-dimensional networks, achieving sparser outcomes, making it particularly suitable for extracting DEC. Additionally, the score function of the BDyMA method allows the incorporation of prior knowledge into the process of dynamic causal discovery which further enhances the accuracy of results. Comprehensive simulations on synthetic data and experiments on Human Connectome Project (HCP) data demonstrate that our method can handle both of the two main challenges, yielding more accurate and reliable DEC compared to state-of-the-art and traditional methods. Additionally, we investigate the trustworthiness of DTI data as prior knowledge for DEC discovery and show the improvements in DEC discovery when the DTI data is incorporated into the process.

The direction of theta and alpha travelling waves modulates human memory processing

To support a range of behaviours, the brain must flexibly coordinate neural activity across widespread brain regions. One potential mechanism for this coordination is a travelling wave, in which a neural oscillation propagates across the brain while organizing the order and timing of activity across regions. Although travelling waves are present across the brain in various species, their potential functional relevance has remained unknown. Here, using rare direct human brain recordings, we demonstrate a distinct functional role for travelling waves of theta- and alpha-band (2-13 Hz) oscillations in the cortex. Travelling waves propagate in different directions during separate cognitive processes. In episodic memory, travelling waves tended to propagate in a posterior-to-anterior direction during successful memory encoding and in an anterior-to-posterior direction during recall. Because travelling waves of oscillations correspond to local neuronal spiking, these patterns indicate that rhythmic pulses of activity move across the brain in different directions for separate behaviours. More broadly, our results suggest a fundamental role for travelling waves and oscillations in dynamically coordinating neural connectivity, by flexibly organizing the timing and directionality of network interactions across the cortex to support cognition and behaviour.

Human connectome topology directs cortical traveling waves and shapes frequency gradients

Traveling waves and neural oscillation frequency gradients are pervasive in the human cortex. While the direction of traveling waves has been linked to brain function and dysfunction, the factors that determine this direction remain elusive. We hypothesized that structural connectivity instrength gradients - defined as the gradually varying sum of incoming connection strengths across the cortex - could shape both traveling wave direction and frequency gradients. We confirm the presence of instrength gradients in the human connectome across diverse cohorts and parcellations. Using a cortical network model, we demonstrate how these instrength gradients direct traveling waves and shape frequency gradients. Our model fits resting-state MEG functional connectivity best in a regime where instrength-directed traveling waves and frequency gradients emerge. We further show how structural subnetworks of the human connectome generate opposing wave directions and frequency gradients observed in the alpha and beta bands. Our findings suggest that structural connectivity instrength gradients affect both traveling wave direction and frequency gradients.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

An Agent-Based Model to Reproduce the Boolean Logic Behaviour of Neuronal Self-Organised Communities through Pulse Delay Modulation and Generation of Logic Gates

The human brain is arguably the most complex "machine" to ever exist. Its detailed functioning is yet to be fully understood, let alone modelled. Neurological processes have logical signal-processing and biophysical aspects, and both affect the brain's structure, functioning and adaptation. Mathematical approaches based on both information and graph theory have been extensively used in an attempt to approximate its biological functioning, along with Artificial Intelligence frameworks inspired by its logical functioning. In this article, an approach to model some aspects of the brain learning and signal processing is presented, mimicking the metastability and backpropagation found in the real brain while also accounting for neuroplasticity. Several simulations are carried out with this model to demonstrate how dynamic neuroplasticity, neural inhibition and neuron migration can reshape the brain's logical connectivity to synchronise signal processing and obtain certain target latencies. This work showcases the importance of dynamic logical and biophysical remodelling in brain plasticity. Combining mathematical (agents, graph theory, topology and backpropagation) and biomedical ingredients (metastability, neuroplasticity and migration), these preliminary results prove complex brain phenomena can be reproduced-under pertinent simplifications-via affordable computations, which can be construed as a starting point for more ambitiously accurate simulations.

A modular and adaptable analysis pipeline to compare slow cerebral rhythms across heterogeneous datasets

Neuroscience is moving toward a more integrative discipline where understanding brain function requires consolidating the accumulated evidence seen across experiments, species, and measurement techniques. A remaining challenge on that path is integrating such heterogeneous data into analysis workflows such that consistent and comparable conclusions can be distilled as an experimental basis for models and theories. Here, we propose a solution in the context of slow-wave activity (<1 Hz), which occurs during unconscious brain states like sleep and general anesthesia and is observed across diverse experimental approaches. We address the issue of integrating and comparing heterogeneous data by conceptualizing a general pipeline design that is adaptable to a variety of inputs and applications. Furthermore, we present the Collaborative Brain Wave Analysis Pipeline (Cobrawap) as a concrete, reusable software implementation to perform broad, detailed, and rigorous comparisons of slow-wave characteristics across multiple, openly available electrocorticography (ECoG) and calcium imaging datasets.

Multifractal long-range dependence pattern of functional magnetic resonance imaging in the human brain at rest

Long-range dependence is a prevalent phenomenon in various biological systems that characterizes the long-memory effect of temporal fluctuations. While recent research suggests that functional magnetic resonance imaging signal has fractal property, it remains unknown about the multifractal long-range dependence pattern of resting-state functional magnetic resonance imaging signals. The current study adopted the multifractal detrended fluctuation analysis on highly sampled resting-state functional magnetic resonance imaging scans to investigate long-range dependence profile associated with the whole-brain voxels as specific functional networks. Our findings revealed the long-range dependence's multifractal properties. Moreover, long-term persistent fluctuations are found for all stations with stronger persistency in whole-brain regions. Subsets with large fluctuations contribute more to the multifractal spectrum in the whole brain. Additionally, we found that the preprocessing with band-pass filtering provided significantly higher reliability for estimating long-range dependence. Our validation analysis confirmed that the optimal pipeline of long-range dependence analysis should include band-pass filtering and removal of daily temporal dependence. Furthermore, multifractal long-range dependence characteristics in healthy control and schizophrenia are different significantly. This work has provided an analytical pipeline for the multifractal long-range dependence in the resting-state functional magnetic resonance imaging signal. The findings suggest differential long-memory effects in the intrinsic functional networks, which may offer a neural marker finding for understanding brain function and pathology.

Shaping dynamical neural computations using spatiotemporal constraints

Dynamics play a critical role in computation. The principled evolution of states over time enables both biological and artificial networks to represent and integrate information to make decisions. In the past few decades, significant multidisciplinary progress has been made in bridging the gap between how we understand biological versus artificial computation, including how insights gained from one can translate to the other. Research has revealed that neurobiology is a key determinant of brain network architecture, which gives rise to spatiotemporally constrained patterns of activity that underlie computation. Here, we discuss how neural systems use dynamics for computation, and claim that the biological constraints that shape brain networks may be leveraged to improve the implementation of artificial neural networks. To formalize this discussion, we consider a natural artificial analog of the brain that has been used extensively to model neural computation: the recurrent neural network (RNN). In both the brain and the RNN, we emphasize the common computational substrate atop which dynamics occur-the connectivity between neurons-and we explore the unique computational advantages offered by biophysical constraints such as resource efficiency, spatial embedding, and neurodevelopment.

The coupling between the spatial and temporal scales of neural processes revealed by a joint time-vertex connectome spectral analysis

Brain oscillations are produced by the coordinated activity of large groups of neurons and different rhythms are thought to reflect different modes of information processing. These modes, in turn, are known to occur at different spatial scales. Nevertheless, how these rhythms support different spatial modes of information processing at the brain scale is not yet fully understood. Here we use "Joint Time-Vertex Spectral Analysis" to characterize the joint spectral content of brain activity both in time (temporal frequencies) and in space over the connectivity graph (spatial connectome harmonics). This method allows us to characterize the relationship between spatially localized or distributed neural processes on one side and their respective temporal frequency bands in source-reconstructed M/EEG signals. We explore this approach on two different datasets, an auditory steady-state response (ASSR) and a visual grating task. Our results suggest that different information processing mechanisms are carried out at different frequency bands: while spatially distributed activity (which may also be interpreted as integration) specifically occurs at low temporal frequencies (alpha and theta) and low graph spatial frequencies, localized electrical activity (i.e., segregation) is observed at high temporal frequencies (high and low gamma) over restricted high spatial graph frequencies. Crucially, the estimated contribution of the distributed and localized neural activity predicts performance in a behavioral task, demonstrating the neurophysiological relevance of the joint time-vertex spectral representation.

Using network control theory to study the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Low-frequency repetitive transcranial magnetic stimulation alters the individual functional dynamical landscape

Repetitive transcranial magnetic stimulation (rTMS) is a noninvasive approach to modulate brain activity and behavior in humans. Still, how individual resting-state brain dynamics after rTMS evolves across different functional configurations is rarely studied. Here, using resting state fMRI data from healthy subjects, we aimed to examine the effects of rTMS to individual large-scale brain dynamics. Using Topological Data Analysis based Mapper approach, we construct the precise dynamic mapping (PDM) for each participant. To reveal the relationship between PDM and canonical functional representation of the resting brain, we annotated the graph using relative activation proportion of a set of large-scale resting-state networks (RSNs) and assigned the single brain volume to corresponding RSN-dominant or a hub state (not any RSN was dominant). Our results show that (i) low-frequency rTMS could induce changed temporal evolution of brain states; (ii) rTMS didn't alter the hub-periphery configurations underlined resting-state brain dynamics; and (iii) the rTMS effects on brain dynamics differ across the left frontal and occipital lobe. In conclusion, low-frequency rTMS significantly alters the individual temporo-spatial dynamics, and our finding further suggested a potential target-dependent alteration of brain dynamics. This work provides a new perspective to comprehend the heterogeneous effect of rTMS.

Pattern selection by material aging: Modeling chemical gardens in two and three dimensions

Chemical gardens are complex, often macroscopic, structures formed by precipitation reactions. Their thin walls compartmentalize the system and adjust in size and shape if the volume of the interior reactant solution is increased by osmosis or active injection. Spatial confinement to a thin layer is known to result in various patterns including self-extending filaments and flower-like patterns organized around a continuous, expanding front. Here, we describe a cellular automaton model for this type of self-organization, in which each lattice site is occupied by one of the two reactants or the precipitate. Reactant injection causes the random replacement of precipitate and generates an expanding near-circular precipitate front. If this process includes an age bias favoring the replacement of fresh precipitate, thin-walled filaments arise and grow-like in the experiments-at the leading tip. In addition, the inclusion of a buoyancy effect allows the model to capture various branched and unbranched chemical garden shapes in two and three dimensions. Our results provide a model of chemical garden structures and highlight the importance of temporal changes in the self-healing membrane material.

Interacting spiral wave patterns underlie complex brain dynamics and are related to cognitive processing

The large-scale activity of the human brain exhibits rich and complex patterns, but the spatiotemporal dynamics of these patterns and their functional roles in cognition remain unclear. Here by characterizing moment-by-moment fluctuations of human cortical functional magnetic resonance imaging signals, we show that spiral-like, rotational wave patterns (brain spirals) are widespread during both resting and cognitive task states. These brain spirals propagate across the cortex while rotating around their phase singularity centres, giving rise to spatiotemporal activity dynamics with non-stationary features. The properties of these brain spirals, such as their rotational directions and locations, are task relevant and can be used to classify different cognitive tasks. We also demonstrate that multiple, interacting brain spirals are involved in coordinating the correlated activations and de-activations of distributed functional regions; this mechanism enables flexible reconfiguration of task-driven activity flow between bottom-up and top-down directions during cognitive processing. Our findings suggest that brain spirals organize complex spatiotemporal dynamics of the human brain and have functional correlates to cognitive processing.

Geometric constraints on human brain function

The anatomy of the brain necessarily constrains its function, but precisely how remains unclear. The classical and dominant paradigm in neuroscience is that neuronal dynamics are driven by interactions between discrete, functionally specialized cell populations connected by a complex array of axonal fibres1-3. However, predictions from neural field theory, an established mathematical framework for modelling large-scale brain activity4-6, suggest that the geometry of the brain may represent a more fundamental constraint on dynamics than complex interregional connectivity7,8. Here, we confirm these theoretical predictions by analysing human magnetic resonance imaging data acquired under spontaneous and diverse task-evoked conditions. Specifically, we show that cortical and subcortical activity can be parsimoniously understood as resulting from excitations of fundamental, resonant modes of the brain's geometry (that is, its shape) rather than from modes of complex interregional connectivity, as classically assumed. We then use these geometric modes to show that task-evoked activations across over 10,000 brain maps are not confined to focal areas, as widely believed, but instead excite brain-wide modes with wavelengths spanning over 60 mm. Finally, we confirm predictions that the close link between geometry and function is explained by a dominant role for wave-like activity, showing that wave dynamics can reproduce numerous canonical spatiotemporal properties of spontaneous and evoked recordings. Our findings challenge prevailing views and identify a previously underappreciated role of geometry in shaping function, as predicted by a unifying and physically principled model of brain-wide dynamics.

Scale-free dynamics in the core-periphery topography and task alignment decline from conscious to unconscious states

Scale-free physiological processes are ubiquitous in the human organism. Resting-state functional MRI studies observed the loss of scale-free dynamics under anesthesia. In contrast, the modulation of scale-free dynamics during task-related activity remains an open question. We investigate scale-free dynamics in the cerebral cortex's unimodal periphery and transmodal core topography in rest and task states during three conscious levels (awake, sedation, and anesthesia) complemented by computational modelling (Stuart-Landau model). The empirical findings demonstrate that the loss of the brain's intrinsic scale-free dynamics in the core-periphery topography during anesthesia, where pink noise transforms into white noise, disrupts the brain's neuronal alignment with the task's temporal structure. The computational model shows that the stimuli's scale-free dynamics, namely pink noise distinguishes from brown and white noise, also modulate task-related activity. Together, we provide evidence for two mechanisms of consciousness, temporo-spatial nestedness and alignment, suggested by the Temporo-Spatial Theory of Consciousness (TTC).

Noise improves the association between effects of local stimulation and structural degree of brain networks

Stimulation to local areas remarkably affects brain activity patterns, which can be exploited to investigate neural bases of cognitive function and modify pathological brain statuses. There has been growing interest in exploring the fundamental action mechanisms of local stimulation. Nevertheless, how noise amplitude, an essential element in neural dynamics, influences stimulation-induced brain states remains unknown. Here, we systematically examine the effects of local stimulation by using a large-scale biophysical model under different combinations of noise amplitudes and stimulation sites. We demonstrate that noise amplitude nonlinearly and heterogeneously tunes the stimulation effects from both regional and network perspectives. Furthermore, by incorporating the role of the anatomical network, we show that the peak frequencies of unstimulated areas at different stimulation sites averaged across noise amplitudes are highly positively related to structural connectivity. Crucially, the association between the overall changes in functional connectivity as well as the alterations in the constraints imposed by structural connectivity with the structural degree of stimulation sites is nonmonotonically influenced by the noise amplitude, with the association increasing in specific noise amplitude ranges. Moreover, the impacts of local stimulation of cognitive systems depend on the complex interplay between the noise amplitude and average structural degree. Overall, this work provides theoretical insights into how noise amplitude and network structure jointly modulate brain dynamics during stimulation and introduces possibilities for better predicting and controlling stimulation outcomes.

Functional re-organization of hippocampal-cortical gradients during naturalistic memory processes

The functional organization of the hippocampus mirrors that of the cortex, changing smoothly along connectivity gradients and abruptly at inter-areal boundaries. Hippocampal-dependent cognitive processes require flexible integration of these hippocampal gradients into functionally related cortical networks. To understand the cognitive relevance of this functional embedding, we acquired fMRI data while participants viewed brief news clips, either containing or lacking recently familiarized cues. Participants were 188 healthy mid-life adults and 31 adults with mild cognitive impairment (MCI) or Alzheimer's disease (AD). We employed a recently developed technique - connectivity gradientography - to study gradually changing patterns of voxel to whole brain functional connectivity and their sudden transitions. We observed that functional connectivity gradients of the anterior hippocampus map onto connectivity gradients across the default mode network during these naturalistic stimuli. The presence of familiar cues in the news clips accentuates a stepwise transition across the boundary from the anterior to the posterior hippocampus. This functional transition is shifted in the posterior direction in the left hippocampus of individuals with MCI or AD. These findings shed new light on the functional integration of hippocampal connectivity gradients into large-scale cortical networks, how these adapt with memory context and how these change in the presence of neurodegenerative disease.

Uncovering the Interplay of Oscillatory Processes During Creative Problem Solving: A Dynamic Modeling Approach

To solve a new problem, people spontaneously engage multiple cognitive processes. Previous work has identified a diverse set of oscillatory components critical at different stages of creative problem solving. In this project, we use hidden state modeling to untangle the roles of oscillation processes over time as people solve puzzles. Building on earlier work, we further developed analytical methods, such as incorporating source separating techniques and identifying the optimal number of states using cross-validation. We extracted brain states characterized by spatio-spectral topographies from time-resolved EEG spectral powers. The data driven approach allowed us to infer the dynamic, trial-by-trial, state sequences, and provided a comprehensive depiction of how various oscillation components interact recurrently throughout the trial. The properties of the states suggest their dissociable cognitive functions. For example, we identified three states with dominant activation in alpha bands but having distinct spatial distributions. People were differentially engaged in these states depending on the stages (e.g., onset or response) and outcomes of the trials (solved with insight or analysis). The current approach, applicable to many tasks requiring extended trial duration, can potentially reconcile findings from previous EEG studies and drive new hypotheses to further our understanding of the complex creative process.

Ion gradient-driven bifurcations of a multi-scale neuronal model

Metabolic limitations within the brain frequently arise in the context of aging and disease. As the largest consumers of energy within the brain, ion pumps that maintain the neuronal membrane potential are the most affected when energy supply becomes limited. To characterize the effects of such limitations, we analyze the ion gradients present in a conductance-based (Morris-Lecar) neural mass model. We show the existence and locations of Neimark-Sacker and period-doubling bifurcations in the sodium, calcium, and potassium reversal potentials and demonstrate that these bifurcations form physiologically relevant bounds of ion gradient variability. Within these bounds, we show how depolarization of the gradients causes decreased neural activity. We also show that the depolarization of ion gradients decreases inter-regional coherence, causing a shift in the critical point at which the coupling occurs and thereby inducing loss of synchrony between regions. In this way, we show that the Larter-Breakspear model captures ion gradient variability present at the microscale level and propagates these changes to the macroscale effects such as those observed in human neuroimaging studies.

Alteration of resting-state network dynamics in autism spectrum disorder based on leading eigenvector dynamics analysis

Background

Neurobiological models to explain the vulnerability of autism spectrum disorders (ASDs) are scarce, and previous resting-state functional magnetic resonance imaging (rs-fMRI) studies mostly examined static functional connectivity (FC). Given that FC constantly evolves, it is critical to probe FC dynamic differences in ASD patients.

Methods

We characterized recurring phase-locking (PL) states during rest in 45 ASD patients and 47 age- and sex-matched healthy controls (HCs) using Leading Eigenvector Dynamics Analysis (LEiDA) and probed the organization of PL states across different fine grain sizes.

Results

Our results identified five different groups of discrete resting-state functional networks, which can be defined as recurrent PL state overtimes. Specifically, ASD patients showed an increased probability of three PL states, consisting of the visual network (VIS), frontoparietal control network (FPN), default mode network (DMN), and ventral attention network (VAN). Correspondingly, ASD patients also showed a decreased probability of two PL states, consisting of the subcortical network (SUB), somatomotor network (SMN), FPN, and VAN.

Conclusion

Our findings suggested that the temporal reorganization of brain discrete networks was closely linked to sensory to cognitive systems of the brain. Our study provides new insights into the dynamics of brain networks and contributes to a deeper understanding of the neurological mechanisms of ASD.

Dynamical interactions reconfigure the gradient of cortical timescales

The functional organization of the brain is usually presented with a back-to-front gradient of timescales, reflecting regional specialization with sensory areas (back) processing information faster than associative areas (front), which perform information integration. However, cognitive processes require not only local information processing but also coordinated activity across regions. Using magnetoencephalography recordings, we find that the functional connectivity at the edge level (between two regions) is also characterized by a back-to-front gradient of timescales following that of the regional gradient. Unexpectedly, we demonstrate a reverse front-to-back gradient when nonlocal interactions are prominent. Thus, the timescales are dynamic and can switch between back-to-front and front-to-back patterns.

Modulations of local synchrony over time lead to resting-state functional connectivity in a parsimonious large-scale brain model

Biophysical models of large-scale brain activity are a fundamental tool for understanding the mechanisms underlying the patterns observed with neuroimaging. These models combine a macroscopic description of the within- and between-ensemble dynamics of neurons within a single architecture. A challenge for these models is accounting for modulations of within-ensemble synchrony over time. Such modulations in local synchrony are fundamental for modeling behavioral tasks and resting-state activity. Another challenge comes from the difficulty in parametrizing large scale brain models which hinders researching principles related with between-ensembles differences. Here we derive a parsimonious large scale brain model that can describe fluctuations of local synchrony. Crucially, we do not reduce within-ensemble dynamics to macroscopic variables first, instead we consider within and between-ensemble interactions similarly while preserving their physiological differences. The dynamics of within-ensemble synchrony can be tuned with a parameter which manipulates local connectivity strength. We simulated resting-state static and time-resolved functional connectivity of alpha band envelopes in models with identical and dissimilar local connectivities. We show that functional connectivity emerges when there are high fluctuations of local and global synchrony simultaneously (i.e. metastable dynamics). We also show that for most ensembles, leaning towards local asynchrony or synchrony correlates with the functional connectivity with other ensembles, with the exception of some regions belonging to the default-mode network.

Spontaneous neuronal avalanches as a correlate of access consciousness

Decades of research have advanced our understanding of the biophysical mechanisms underlying consciousness. However, an overarching framework bridging between models of consciousness and the large-scale organization of spontaneous brain activity is still missing. Based on the observation that spontaneous brain activity dynamically switches between epochs of segregation and large-scale integration of information, we hypothesize a brain-state dependence of conscious access, whereby the presence of either segregated or integrated states marks distinct modes of information processing. We first review influential works on the neuronal correlates of consciousness, spontaneous resting-state brain activity and dynamical system theory. Then, we propose a test experiment to validate our hypothesis that conscious access occurs in aperiodic cycles, alternating windows where new incoming information is collected but not experienced, to punctuated short-lived integration events, where conscious access to previously collected content occurs. In particular, we suggest that the integration events correspond to neuronal avalanches, which are collective bursts of neuronal activity ubiquitously observed in electrophysiological recordings. If confirmed, the proposed framework would link the physics of spontaneous cortical dynamics, to the concept of ignition within the global neuronal workspace theory, whereby conscious access manifest itself as a burst of neuronal activity.

How the Brain Becomes the Mind: Can Thermodynamics Explain the Emergence and Nature of Emotions?

The neural systems' electric activities are fundamental for the phenomenology of consciousness. Sensory perception triggers an information/energy exchange with the environment, but the brain's recurrent activations maintain a resting state with constant parameters. Therefore, perception forms a closed thermodynamic cycle. In physics, the Carnot engine is an ideal thermodynamic cycle that converts heat from a hot reservoir into work, or inversely, requires work to transfer heat from a low- to a high-temperature reservoir (the reversed Carnot cycle). We analyze the high entropy brain by the endothermic reversed Carnot cycle. Its irreversible activations provide temporal directionality for future orientation. A flexible transfer between neural states inspires openness and creativity. In contrast, the low entropy resting state parallels reversible activations, which impose past focus via repetitive thinking, remorse, and regret. The exothermic Carnot cycle degrades mental energy. Therefore, the brain's energy/information balance formulates motivation, sensed as position or negative emotions. Our work provides an analytical perspective of positive and negative emotions and spontaneous behavior from the free energy principle. Furthermore, electrical activities, thoughts, and beliefs lend themselves to a temporal organization, an orthogonal condition to physical systems. Here, we suggest that an experimental validation of the thermodynamic origin of emotions might inspire better treatment options for mental diseases.

The expanding horizons of network neuroscience: From description to prediction and control

The field of network neuroscience has emerged as a natural framework for the study of the brain and has been increasingly applied across divergent problems in neuroscience. From a disciplinary perspective, network neuroscience originally emerged as a formal integration of graph theory (from mathematics) and neuroscience (from biology). This early integration afforded marked utility in describing the interconnected nature of neural units, both structurally and functionally, and underscored the relevance of that interconnection for cognition and behavior. But since its inception, the field has not remained static in its methodological composition. Instead, it has grown to use increasingly advanced graph-theoretic tools and to bring in several other disciplinary perspectives-including machine learning and systems engineering-that have proven complementary. In doing so, the problem space amenable to the discipline has expanded markedly. In this review, we discuss three distinct flavors of investigation in state-of-the-art network neuroscience: (i) descriptive network neuroscience, (ii) predictive network neuroscience, and (iii) a perturbative network neuroscience that draws on recent advances in network control theory. In considering each area, we provide a brief summary of the approaches, discuss the nature of the insights obtained, and highlight future directions.

Precision dynamical mapping using topological data analysis reveals a hub-like transition state at rest

In the absence of external stimuli, neural activity continuously evolves from one configuration to another. Whether these transitions or explorations follow some underlying arrangement or lack a predictable ordered plan remains to be determined. Here, using fMRI data from highly sampled individuals (~5 hours of resting-state data per individual), we aimed to reveal the rules that govern transitions in brain activity at rest. Our Topological Data Analysis based Mapper approach characterized a highly visited transition state of the brain that acts as a switch between different neural configurations to organize the spontaneous brain activity. Further, while the transition state was characterized by a uniform representation of canonical resting-state networks (RSNs), the periphery of the landscape was dominated by a subject-specific combination of RSNs. Altogether, we revealed rules or principles that organize spontaneous brain activity using a precision dynamics approach.

Brain Complexities and Anesthesia: Their Meaning and Measurement

A complex system is often associated with emergence of new phenomena from the interactions between the system's components. General anesthesia reduces brain complexity and so inhibits the emergence of consciousness. An understanding of complexity is necessary for the interpretation of brain monitoring algorithms. Complexity indices capture the "difficulty" of understanding brain activity over time and/or space. Complexity-entropy plots reveal the types of complexity indices and their balance of randomness and structure. Lempel-Ziv complexity is a common index of temporal complexity for single-channel electroencephalogram containing both power spectral and nonlinear effects, revealed by phase-randomized surrogate data. Computing spatial complexities involves forming a connectivity matrix and calculating the complexity of connectivity patterns. Spatiotemporal complexity can be estimated in multiple ways including temporal or spatial concatenation, estimation of state switching, or integrated information. This article illustrates the concept and application of various complexities by providing working examples; a website with interactive demonstrations has also been created.

Metastable oscillatory modes emerge from synchronization in the brain spacetime connectome

A rich repertoire of oscillatory signals is detected from human brains with electro- and magnetoencephalography (EEG/MEG). However, the principles underwriting coherent oscillations and their link with neural activity remain under debate. Here, we revisit the mechanistic hypothesis that transient brain rhythms are a signature of metastable synchronization, occurring at reduced collective frequencies due to delays between brain areas. We consider a system of damped oscillators in the presence of background noise - approximating the short-lived gamma-frequency oscillations generated within neuronal circuits - coupled according to the diffusion weighted tractography between brain areas. Varying the global coupling strength and conduction speed, we identify a critical regime where spatially and spectrally resolved metastable oscillatory modes (MOMs) emerge at sub-gamma frequencies, approximating the MEG power spectra from 89 healthy individuals at rest. Further, we demonstrate that the frequency, duration, and scale of MOMs - as well as the frequency-specific envelope functional connectivity - can be controlled by global parameters, while the connectome structure remains unchanged. Grounded in the physics of delay-coupled oscillators, these numerical analyses demonstrate how interactions between locally generated fast oscillations in the connectome spacetime structure can lead to the emergence of collective brain rhythms organized in space and time.

Generative Models of Brain Dynamics

This review article gives a high-level overview of the approaches across different scales of organization and levels of abstraction. The studies covered in this paper include fundamental models in computational neuroscience, nonlinear dynamics, data-driven methods, as well as emergent practices. While not all of these models span the intersection of neuroscience, AI, and system dynamics, all of them do or can work in tandem as generative models, which, as we argue, provide superior properties for the analysis of neuroscientific data. We discuss the limitations and unique dynamical traits of brain data and the complementary need for hypothesis- and data-driven modeling. By way of conclusion, we present several hybrid generative models from recent literature in scientific machine learning, which can be efficiently deployed to yield interpretable models of neural dynamics.

Uncovering hidden resting state dynamics: A new perspective on auditory verbal hallucinations

In the absence of sensory stimulation, the brain transits between distinct functional networks. Network dynamics such as transition patterns and the time the brain stays in each network link to cognition and behavior and are subject to much investigation. Auditory verbal hallucinations (AVH), the temporally fluctuating unprovoked experience of hearing voices, are associated with aberrant resting state network activity. However, we lack a clear understanding of how different networks contribute to aberrant activity over time. An accurate characterization of latent network dynamics and their relation to neurocognitive changes necessitates methods that capture the sub-second temporal fluctuations of the networks' functional connectivity signatures. Here, we critically evaluate the assumptions and sensitivity of several approaches commonly used to assess temporal dynamics of brain connectivity states in M/EEG and fMRI research, highlighting methodological constraints and their clinical relevance to AVH. Identifying altered brain connectivity states linked to AVH can facilitate the detection of predictive disease markers and ultimately be valuable for generating individual risk profiles, differential diagnosis, targeted intervention, and treatment strategies.

On the intersection between data quality and dynamical modelling of large-scale fMRI signals

Large-scale dynamics of the brain are routinely modelled using systems of nonlinear dynamical equations that describe the evolution of population-level activity, with distinct neural populations often coupled according to an empirically measured structural connectivity matrix. This modelling approach has been used to generate insights into the neural underpinnings of spontaneous brain dynamics, as recorded with techniques such as resting state functional MRI (fMRI). In fMRI, researchers have many degrees of freedom in the way that they can process the data and recent evidence indicates that the choice of pre-processing steps can have a major effect on empirical estimates of functional connectivity. However, the potential influence of such variations on modelling results are seldom considered. Here we show, using three popular whole-brain dynamical models, that different choices during fMRI preprocessing can dramatically affect model fits and interpretations of findings. Critically, we show that the ability of these models to accurately capture patterns in fMRI dynamics is mostly driven by the degree to which they fit global signals rather than interesting sources of coordinated neural dynamics. We show that widespread deflections can arise from simple global synchronisation. We introduce a simple two-parameter model that captures these fluctuations and performs just as well as more complex, multi-parameter biophysical models. From our combined analyses of data and simulations, we describe benchmarks to evaluate model fit and validity. Although most models are not resilient to denoising, we show that relaxing the approximation of homogeneous neural populations by more explicitly modelling inter-regional effective connectivity can improve model accuracy at the expense of increased model complexity. Our results suggest that many complex biophysical models may be fitting relatively trivial properties of the data, and underscore a need for tighter integration between data quality assurance and model development.

Spatiotemporal Patterns of Adaptation-Induced Slow Oscillations in a Whole-Brain Model of Slow-Wave Sleep

During slow-wave sleep, the brain is in a self-organized regime in which slow oscillations (SOs) between up- and down-states travel across the cortex. While an isolated piece of cortex can produce SOs, the brain-wide propagation of these oscillations are thought to be mediated by the long-range axonal connections. We address the mechanism of how SOs emerge and recruit large parts of the brain using a whole-brain model constructed from empirical connectivity data in which SOs are induced independently in each brain area by a local adaptation mechanism. Using an evolutionary optimization approach, good fits to human resting-state fMRI data and sleep EEG data are found at values of the adaptation strength close to a bifurcation where the model produces a balance between local and global SOs with realistic spatiotemporal statistics. Local oscillations are more frequent, last shorter, and have a lower amplitude. Global oscillations spread as waves of silence across the undirected brain graph, traveling from anterior to posterior regions. These traveling waves are caused by heterogeneities in the brain network in which the connection strengths between brain areas determine which areas transition to a down-state first, and thus initiate traveling waves across the cortex. Our results demonstrate the utility of whole-brain models for explaining the origin of large-scale cortical oscillations and how they are shaped by the connectome.

Cortical traveling waves reflect state-dependent hierarchical sequencing of local regions in the human connectome network

Recent human studies using electrocorticography have demonstrated that alpha and theta band oscillations form traveling waves on the cortical surface. According to neural synchronization theories, the cortical traveling waves may group local cortical regions and sequence them by phase synchronization; however these contributions have not yet been assessed. This study aimed to evaluate the functional contributions of traveling waves using connectome-based network modeling. In the simulation, we observed stable traveling waves on the entire cortical surface wherein the topographical pattern of these phases was substantially correlated with the empirically obtained resting-state networks, and local radial waves also appeared within the size of the empirical networks (< 50 mm). Importantly, individual regions in the entire network were instantaneously sequenced by their internal frequencies, and regions with higher intrinsic frequency were seen in the earlier phases of the traveling waves. Based on the communication-through-coherence theory, this phase configuration produced a hierarchical organization of each region by unidirectional communication between the arbitrarily paired regions. In conclusion, cortical traveling waves reflect the intrinsic frequency-dependent hierarchical sequencing of local regions, global traveling waves sequence the set of large-scale cortical networks, and local traveling waves sequence local regions within individual cortical networks.

Novel Cognitive Functions Arise at the Convergence of Macroscale Gradients

Functions in higher-order brain regions are the source of extensive debate. Although past trends have been to describe the brain-especially posterior cortical areas-in terms of a set of functional modules, a new emerging paradigm focuses on the integration of proximal functions. In this review, we synthesize emerging evidence that a variety of novel functions in the higher-order brain regions are due to convergence: convergence of macroscale gradients brings feature-rich representations into close proximity, presenting an opportunity for novel functions to arise. Using the TPJ as an example, we demonstrate that convergence is enabled via three properties of the brain: (1) hierarchical organization, (2) abstraction, and (3) equidistance. As gradients travel from primary sensory cortices to higher-order brain regions, information becomes abstracted and hierarchical, and eventually, gradients meet at a point maximally and equally distant from their sensory origins. This convergence, which produces multifaceted combinations, such as mentalizing another person's thought or projecting into a future space, parallels evolutionary and developmental characteristics in such regions, resulting in new cognitive and affective faculties.

Whole-brain propagating patterns in human resting-state brain activities

Repetitive propagating activities in resting-state brain activities have been widely observed in various species and regions. Because they resemble the preceding brain activities during tasks, they are assumed to reflect past experiences embedded in neuronal circuits. "Whole-brain" propagating activities may also reflect a process that integrates information distributed over the entire brain, such as visual and motor information. Here we reveal whole-brain propagating activities from human resting-state magnetoencephalography (MEG) and electroencephalography (EEG) data. We simultaneously recorded the MEGs and EEGs and estimated the source currents from both measurements. Then using our recently proposed algorithm, we extracted repetitive spatiotemporal patterns from the source currents. The estimated patterns consisted of multiple frequency components, each of which transiently exhibited the frequency-specific resting-state networks (RSNs) of functional MRIs (fMRIs), such as the default mode and sensorimotor networks. A simulation test suggested that the spatiotemporal patterns reflected the phase alignment of the multiple frequency oscillators induced by the propagating activities along the anatomical connectivity. These results argue that whole-brain propagating activities transiently exhibited multiple RSNs in their multiple frequency components, suggesting that they reflected a process to integrate the information distributed over the frequencies and networks.

Criticality Creates a Functional Platform for Network Transitions Between Internal and External Processing Modes in the Human Brain

Continuous switching between internal and external modes in the brain appears important for generating models of the self and the world. However, how the brain transitions between these two modes remains unknown. We propose that a large synchronization fluctuation of brain networks, emerging only near criticality (i.e., a balanced state between order and disorder), spontaneously creates temporal windows with distinct preferences for integrating the network's internal information or for processing external stimuli. Using a computational model, electroencephalography (EEG) analysis, and functional magnetic resonance imaging (fMRI) analysis during alterations of consciousness in humans, we report that synchronized and incoherent networks, respectively, bias toward internal and external information with specific network configurations. In the brain network model and EEG-based network, the network preferences are the most prominent at criticality and in conscious states associated with the bandwidth 4-12 Hz, with alternating functional network configurations. However, these network configurations are selectively disrupted in different states of consciousness such as general anesthesia, psychedelic states, minimally conscious states, and unresponsive wakefulness syndrome. The network preference for internal information integration is only significant in conscious states and psychedelic states, but not in other unconscious states, suggesting the importance of internal information integration in maintaining consciousness. The fMRI co-activation pattern analysis shows that functional networks that are sensitive to external stimuli-such as default mode, dorsal attentional, and frontoparietal networks-are activated in incoherent states, while insensitive networks, such as global activation and deactivation networks, are dominated in highly synchronized states. We suggest that criticality produces a functional platform for the brain's capability for continuous switching between two modes, which is crucial for the emergence of consciousness.

neurolib: A Simulation Framework for Whole-Brain Neural Mass Modeling

neurolib is a computational framework for whole-brain modeling written in Python. It provides a set of neural mass models that represent the average activity of a brain region on a mesoscopic scale. In a whole-brain network model, brain regions are connected with each other based on biologically informed structural connectivity, i.e., the connectome of the brain. neurolib can load structural and functional datasets, set up a whole-brain model, manage its parameters, simulate it, and organize its outputs for later analysis. The activity of each brain region can be converted into a simulated BOLD signal in order to calibrate the model against empirical data from functional magnetic resonance imaging (fMRI). Extensive model analysis is made possible using a parameter exploration module, which allows one to characterize a model’s behavior as a function of changing parameters. An optimization module is provided for fitting models to multimodal empirical data using evolutionary algorithms. neurolib is designed to be extendable and allows for easy implementation of custom neural mass models, offering a versatile platform for computational neuroscientists for prototyping models, managing large numerical experiments, studying the structure–function relationship of brain networks, and for performing in-silico optimization of whole-brain models.

Intensity patterns at the peaks of brain activity in fMRI and PET are highly correlated with neural models of spatial integration

Spatial integration during the brain's cognitive activity prompts changes in energy used by different neuroglial populations. Nevertheless, the organisation of such integration in 3D ‐brain activity remains undescribed from a quantitative standpoint. In response, we applied a cross‐correlation between brain activity and integrative models, which yielded a deeper understanding of information integration in functional brain mapping. We analysed four datasets obtained via fundamentally different neuroimaging techniques (functional magnetic resonance imaging [fMRI] and positron emission tomography [PET]) and found that models of spatial integration with an increasing input to each step of integration were significantly more correlated with brain activity than models with a constant input to each step of integration. In addition, marking the voxels with the maximal correlation, we found exceptionally high intersubject consistency with the initial brain activity at the peaks. Our method demonstrated for the first time that the network of peaks of brain activity is organised strictly according to the models of spatial integration independent of neuroimaging techniques. The highest correlation with models integrating an increasing at each step input suggests that brain activity reflects a network of integrative processes where the results of integration in some neuroglial populations serve as an input to other neuroglial populations.

Dissecting whole-brain conduction delays through MRI microstructural measures

Network models based on structural connectivity have been increasingly used as the blueprint for large-scale simulations of the human brain. As the nodes of this network are distributed through the cortex and interconnected by white matter pathways with different characteristics, modeling the associated conduction delays becomes important. The goal of this study is to estimate and characterize these delays directly from the brain structure. To achieve this, we leveraged microstructural measures from a combination of advanced magnetic resonance imaging acquisitions and computed the main determinants of conduction velocity, namely axonal diameter and myelin content. Using the model proposed by Rushton, we used these measures to calculate the conduction velocity and estimated the associated delays using tractography. We observed that both the axonal diameter and conduction velocity distributions presented a rather constant trend across different connection lengths, with resulting delays that scale linearly with the connection length. Relying on insights from graph theory and Kuramoto simulations, our results support the approximation of constant conduction velocity but also show path- and region-specific differences.

Determination of Dynamic Brain Connectivity via Spectral Analysis

Spectral analysis based on neural field theory is used to analyze dynamic connectivity via methods based on the physical eigenmodes that are the building blocks of brain dynamics. These approaches integrate over space instead of averaging over time and thereby greatly reduce or remove the temporal averaging effects, windowing artifacts, and noise at fine spatial scales that have bedeviled the analysis of dynamical functional connectivity (FC). The dependences of FC on dynamics at various timescales, and on windowing, are clarified and the results are demonstrated on simple test cases, demonstrating how modes provide directly interpretable insights that can be related to brain structure and function. It is shown that FC is dynamic even when the brain structure and effective connectivity are fixed, and that the observed patterns of FC are dominated by relatively few eigenmodes. Common artifacts introduced by statistical analyses that do not incorporate the physical nature of the brain are discussed and it is shown that these are avoided by spectral analysis using eigenmodes. Unlike most published artificially discretized “resting state networks” and other statistically-derived patterns, eigenmodes overlap, with every mode extending across the whole brain and every region participating in every mode—just like the vibrations that give rise to notes of a musical instrument. Despite this, modes are independent and do not interact in the linear limit. It is argued that for many purposes the intrinsic limitations of covariance-based FC instead favor the alternative of tracking eigenmode coefficients vs. time, which provide a compact representation that is directly related to biophysical brain dynamics.

The Radically Embodied Conscious Cybernetic Bayesian Brain: From Free Energy to Free Will and Back Again

Drawing from both enactivist and cognitivist perspectives on mind, I propose that explaining teleological phenomena may require reappraising both "Cartesian theaters" and mental homunculi in terms of embodied self-models (ESMs), understood as body maps with agentic properties, functioning as predictive-memory systems and cybernetic controllers. Quasi-homuncular ESMs are suggested to constitute a major organizing principle for neural architectures due to their initial and ongoing significance for solutions to inference problems in cognitive (and affective) development. Embodied experiences provide foundational lessons in learning curriculums in which agents explore increasingly challenging problem spaces, so answering an unresolved question in Bayesian cognitive science: what are biologically plausible mechanisms for equipping learners with sufficiently powerful inductive biases to adequately constrain inference spaces? Drawing on models from neurophysiology, psychology, and developmental robotics, I describe how embodiment provides fundamental sources of empirical priors (as reliably learnable posterior expectations). If ESMs play this kind of foundational role in cognitive development, then bidirectional linkages will be found between all sensory modalities and frontal-parietal control hierarchies, so infusing all senses with somatic-motoric properties, thereby structuring all perception by relevant affordances, so solving frame problems for embodied agents. Drawing upon the Free Energy Principle and Active Inference framework, I describe a particular mechanism for intentional action selection via consciously imagined (and explicitly represented) goal realization, where contrasts between desired and present states influence ongoing policy selection via predictive coding mechanisms and backward-chained imaginings (as self-realizing predictions). This embodied developmental legacy suggests a mechanism by which imaginings can be intentionally shaped by (internalized) partially-expressed motor acts, so providing means of agentic control for attention, working memory, imagination, and behavior. I further describe the nature(s) of mental causation and self-control, and also provide an account of readiness potentials in Libet paradigms wherein conscious intentions shape causal streams leading to enaction. Finally, I provide neurophenomenological handlings of prototypical qualia including pleasure, pain, and desire in terms of self-annihilating free energy gradients via quasi-synesthetic interoceptive active inference. In brief, this manuscript is intended to illustrate how radically embodied minds may create foundations for intelligence (as capacity for learning and inference), consciousness (as somatically-grounded self-world modeling), and will (as deployment of predictive models for enacting valued goals).

FNS allows efficient event-driven spiking neural network simulations based on a neuron model supporting spike latency

Neural modelling tools are increasingly employed to describe, explain, and predict the human brain's behavior. Among them, spiking neural networks (SNNs) make possible the simulation of neural activity at the level of single neurons, but their use is often threatened by the resources needed in terms of processing capabilities and memory. Emerging applications where a low energy burden is required (e.g. implanted neuroprostheses) motivate the exploration of new strategies able to capture the relevant principles of neuronal dynamics in reduced and efficient models. The recent Leaky Integrate-and-Fire with Latency (LIFL) spiking neuron model shows some realistic neuronal features and efficiency at the same time, a combination of characteristics that may result appealing for SNN-based brain modelling. In this paper we introduce FNS, the first LIFL-based SNN framework, which combines spiking/synaptic modelling with the event-driven approach, allowing us to define heterogeneous neuron groups and multi-scale connectivity, with delayed connections and plastic synapses. FNS allows multi-thread, precise simulations, integrating a novel parallelization strategy and a mechanism of periodic dumping. We evaluate the performance of FNS in terms of simulation time and used memory, and compare it with those obtained with neuronal models having a similar neurocomputational profile, implemented in NEST, showing that FNS performs better in both scenarios. FNS can be advantageously used to explore the interaction within and between populations of spiking neurons, even for long time-scales and with a limited hardware configuration.

Maximum Entropy Principle Underlies Wiring Length Distribution in Brain Networks

A brain network comprises a substantial amount of short-range connections with an admixture of long-range connections. The portion of long-range connections in brain networks is observed to be quantitatively dissimilar across species. It is hypothesized that the length of connections is constrained by the spatial embedding of brain networks, yet fundamental principles that underlie the wiring length distribution remain unclear. By quantifying the structural diversity of a brain network using Shannon's entropy, here we show that the wiring length distribution across multiple species-including Drosophila, mouse, macaque, human, and C. elegans-follows the maximum entropy principle (MAP) under the constraints of limited wiring material and the spatial locations of brain areas or neurons. In addition, by considering stochastic axonal growth, we propose a network formation process capable of reproducing wiring length distributions of the 5 species, thereby implementing MAP in a biologically plausible manner. We further develop a generative model incorporating MAP, and show that, for the 5 species, the generated network exhibits high similarity to the real network. Our work indicates that the brain connectivity evolves to be structurally diversified by maximizing entropy to support efficient interareal communication, providing a potential organizational principle of brain networks.

Mean-Field Models for EEG/MEG: From Oscillations to Waves

Neural mass models have been used since the 1970s to model the coarse-grained activity of large populations of neurons. They have proven especially fruitful for understanding brain rhythms. However, although motivated by neurobiological considerations they are phenomenological in nature, and cannot hope to recreate some of the rich repertoire of responses seen in real neuronal tissue. Here we consider a simple spiking neuron network model that has recently been shown to admit an exact mean-field description for both synaptic and gap-junction interactions. The mean-field model takes a similar form to a standard neural mass model, with an additional dynamical equation to describe the evolution of within-population synchrony. As well as reviewing the origins of this next generation mass model we discuss its extension to describe an idealised spatially extended planar cortex. To emphasise the usefulness of this model for EEG/MEG modelling we show how it can be used to uncover the role of local gap-junction coupling in shaping large scale synaptic waves.