Search

The stability of complex networks, from power grids to biological systems, is crucial for their proper functioning. It is thus important to control such systems to maintain or restore their stability. Traditional approaches rely on real-time state measurements for feedback control, but this can be challenging in many real-world systems, such as the brain, due to their complex and dynamic nature. This paper utilizes vibrational control -- an open-loop strategy -- to regulate network stability. Unlike conventional methods targeting network nodes, our approach focuses on manipulating network edges through vibrational inputs. We establish sufficient graph-theoretic conditions for vibration-induced functional modifications of network edges and stabilization of network systems as a whole. Additionally, we provide methods for designing effective vibrational control inputs and validate our theoretical findings through numerical simulations., Comment: 12 pages, 9 figures

Epilepsy is one of the most common neurological disorders globally, affecting millions of individuals. Despite significant advancements, the precise mechanisms underlying this condition remain largely unknown, making accurately predicting and preventing epileptic seizures challenging. In this paper, we employ a bistable model, where a stable equilibrium and a stable limit cycle coexist, to describe epileptic dynamics. The equilibrium captures normal steady-state neural activity, while the stable limit cycle signifies seizure-like oscillations. The noise-driven switch from the equilibrium to the limit cycle characterizes the onset of seizures. The differences in the regions of attraction of these two stable states distinguish epileptic brain dynamics from healthy ones. We analytically construct the regions of attraction for both states. Further, using the notion of input-to-state stability, we theoretically show how the regions of attraction influence the stability of the system subject to external perturbations. Generalizing the bistable system into coupled networks, we also find the role of network parameters in shaping the regions of attraction. Our findings shed light on the intricate interplay between brain networks and epileptic activity, offering mechanistic insights into potential avenues for more predictable treatments., Comment: 6 pages, 4 figures, submitted to IEEE CDC 2024

Cluster synchronization is of paramount importance for the normal functioning of numerous technological and natural systems. Deviations from normal cluster synchronization patterns are closely associated with various malfunctions, such as neurological disorders in the brain. Therefore, it is crucial to restore normal system functions by stabilizing the appropriate cluster synchronization patterns. Most existing studies focus on designing controllers based on state measurements to achieve system stabilization. However, in many real-world scenarios, measuring system states, such as neuronal activity in the brain, poses significant challenges, rendering the stabilization of such systems difficult. To overcome this challenge, in this paper, we employ an open-loop control strategy, vibrational control, which does not requires any state measurements. We establish some sufficient conditions under which vibrational inputs stabilize cluster synchronization. Further, we provide a tractable approach to design vibrational control. Finally, numerical experiments are conducted to demonstrate our theoretical findings., Comment: Submitted to Open Journal of Control Systems

Many natural and man-made network systems need to maintain certain patterns, such as working at equilibria or limit cycles, to function properly. Thus, the ability to stabilize such patterns is crucial. Most of the existing studies on stabilization assume that network systems states can be measured online so that feedback control strategies can be used. However, in many real-world scenarios, systems states, e.g., neuronal activity in the brain, are often difficult to measure. In this paper, we take this situation into account and study the stabilization problem of linear network systems with an open-loop control strategy (vibrational control). We derive a graph-theoretic sufficient condition for structural vibrational stabilizability, under which network systems can always be stabilized. We further provide an approach to select the locations in the network for control placement and design corresponding vibrational inputs to stabilize systems that satisfy this condition. Finally, we provide some numerical results that demonstrate the validity of our theoretical findings.

Cluster synchronization underlies various functions in the brain. Abnormal patterns of cluster synchronization are often associated with neurological disorders. Deep brain stimulation (DBS) is a neurosurgical technique used to treat several brain diseases, which has been observed to regulate neuronal synchrony patterns. Despite its widespread use, the mechanisms of DBS remain largely unknown. In this paper, we hypothesize that DBS plays a role similar to vibrational control since they both highly rely on high-frequency excitation to function. Under the framework of Kuramoto-oscillator networks, we study how vibrations introduced to network connections can stabilize cluster synchronization. We derive some sufficient conditions and also provide an effective approach to design vibrational control. Also, a numerical example is presented to demonstrate our theoretical findings., Comment: 7 pages, 4 figures

To design control strategies that predictably manipulate a system's behavior, it is first necessary to understand how the system's structure relates to its response. Many complex systems can be represented as multilayer networks whose response and control can be studied in the framework of network control theory. It remains unknown how a system's layered architecture dictates its control properties, particularly when control signals can only access the system through a single input layer. Here, we use the framework of linear control theory to probe the control properties of a duplex network with directed interlayer links. We determine the manner in which the structural properties of layers and relative interlayer arrangement together dictate the system's response. For this purpose, we calculate the exact expression of optimal control energy in terms of layer spectra and the relative alignment between the eigenmodes of the input layer and the deeper target layer. For a range of numerically constructed duplex networks, we then calculate the control properties of the two layers as a function of target-layer densities and layer topology. The alignment of layer eigenmodes emerges as an important parameter that sets the cost and routing of the optimal energy for control. We understand these results in a simplified limit of a single-mode approximation, and we build metrics to characterize the routing of optimal control energy through the eigenmodes of each layer. Our analytical and numerical results together provide insights into the relationships between structure and control in duplex networks. More generally, they serve as a platform for future work designing optimal \emph{interlayer} control strategies., Comment: 14 pages, 5 figures, 19 pages supplementary material, 10 supplementary figures

Oscillatory activity is ubiquitous in natural and engineered network systems. The interaction scheme underlying interdependent oscillatory components governs the emergence of network-wide patterns of synchrony that regulate and enable complex functions. Yet, understanding, and ultimately harnessing, the structure-function relationship in oscillator networks remains an outstanding challenge of modern science. Here, we address this challenge by presenting a principled method to prescribe exact and robust functional configurations from local network interactions through optimal tuning of the oscillators' parameters. To quantify the behavioral synchrony between coupled oscillators, we introduce the notion of functional pattern, which encodes the pairwise relationships between the oscillators' phases. Our procedure is computationally efficient and provably correct, accounts for constrained interaction types, and allows to concurrently assign multiple desired functional patterns. Further, we derive algebraic and graph-theoretic conditions to guarantee the feasibility and stability of target functional patterns. These conditions provide an interpretable mapping between the structural constraints and their functional implications in oscillator networks. As a proof of concept, we apply the proposed method to replicate empirically recorded functional relationships from cortical oscillations in a human brain, and to redistribute the active power flow in different models of electrical grids.

Across all scales of the physical world, dynamical systems can often be usefully represented as abstract networks that encode the system's units and inter-unit interactions. Understanding how physical rules shape the topological structure of those networks can clarify a system's function and enhance our ability to design, guide, or control its behavior. In the emerging area of quantum network science, a key challenge lies in distinguishing between the topological properties that reflect a system's underlying physics and those that reflect the assumptions of the employed conceptual model. To elucidate and address this challenge, we study networks that represent non-equilibrium quantum-electronic transport through quantum antidot devices -- an example of an open, mesoscopic quantum system. The network representations correspond to two different models of internal antidot states: a single-particle, non-interacting model and an effective model for collective excitations including Coulomb interactions. In these networks, nodes represent accessible energy states and edges represent allowed transitions. We find that both models reflect spin conservation rules in the network topology through bipartiteness and the presence of only even-length cycles. The models diverge, however, in the minimum length of cycle basis elements, in a manner that depends on whether electrons are considered to be distinguishable. Furthermore, the two models reflect spin-conserving relaxation effects differently, as evident in both the degree distribution and the cycle-basis length distribution. Collectively, these observations serve to elucidate the relationship between network structure and physical constraints in quantum-mechanical models. More generally, our approach underscores the utility of network science in understanding the dynamics and control of quantum systems., Comment: 37 pages, including supplementary material

From spiking activity in neuronal networks to force chains in granular materials, the behavior of many real-world systems depends on a network of both strong and weak interactions. These interactions give rise to complex and higher-order system behaviors, and are encoded using data as the network's edges. However, distinguishing between true weak edges and low-weight edges caused by noise remains a challenge. We address this problem by examining the higher-order structure of noisy, weak edges added to model networks. We find that the structure of low-weight, noisy edges varies according to the topology of the model network to which it is added. By investigating this variation more closely, we see that at least three qualitative classes of noise structure emerge. Furthermore, we observe that the structure of noisy edges contains enough model-specific information to classify the model networks with moderate accuracy. Finally, we offer network generation rules that can drive different types of structure in added noisy edges. Our results demonstrate that noise does not present as a monolithic nuisance, but rather as a nuanced, topology-dependent, and even useful entity in characterizing higher-order network interactions. Hence, we provide an alternate approach to noise management by embracing its role in such interactions., Comment: 15 pages, 6 figures (main text)

Probing the developing neural circuitry in Caenorhabditis elegans has enhanced our understanding of nervous systems. The C. elegans connectome, like those of other species, is characterized by a rich club of densely connected neurons embedded within a small-world architecture. This organization of neuronal connections, captured by quantitative network statistics, provides insight into the system's capacity to perform integrative computations. Yet these network measures are limited in their ability to detect weakly connected motifs, such as topological cavities, that may support the systems capacity to perform segregated computations. We address this limitation by using persistent homology to track the evolution of topological cavities in the growing C. elegans connectome throughout neural development, and assess the degree to which the growing connectomes topology is resistant to biological noise. We show that the developing connectome topology is both relatively robust to changes in neuron birth times and not captured by similar growth models. Additionally, we quantify the consequence of a neurons specific birth time and ask if this metric tracks other biological properties of neurons. Our results suggest that the connectomes growing topology is a robust feature of the developing connectome that is distinct from other network properties, and that the growing topology is particularly sensitive to the exact birth times of a small set of predominantly motor neurons. By utilizing novel measurements that track biological features, we anticipate that our study will be helpful in the construction of more accurate models of neuronal development in C. elegans

Functional connectivity (FC) can be represented as a network, and is frequently used to better understand the neural underpinnings of complex tasks such as motor imagery (MI) detection in brain-computer interfaces (BCIs). However, errors in the estimation of connectivity can affect the detection performances. In this work, we address the problem of denoising common connectivity estimates to improve the detectability of different connectivity states. Specifically, we propose a denoising algorithm that acts on the network graph Laplacian, which leverages recent graph signal processing results. Further, we derive a novel formulation of the Jensen divergence for the denoised Laplacian under different states. Numerical simulations on synthetic data show that the denoising method improves the Jensen divergence of connectivity patterns corresponding to different task conditions. Furthermore, we apply the Laplacian denoising technique to brain networks estimated from real EEG data recorded during MI-BCI experiments. Using our novel formulation of the J-divergence, we are able to quantify the distance between the FC networks in the motor imagery and resting states, as well as to understand the contribution of each Laplacian variable to the total J-divergence between two states. Experimental results on real MI-BCI EEG data demonstrate that the Laplacian denoising improves the separation of motor imagery and resting mental states, and shortens the time interval required for connectivity estimation. We conclude that the approach shows promise for the robust detection of connectivity states while being appealing for implementation in real-time BCI applications., Comment: This work has been submitted to the IEEE for possible publication

Phase-amplitude coupling (PAC) describes the phenomenon where the power of a high-frequency oscillation evolves with the phase of a low-frequency one. We propose a model that explains the emergence of PAC in two commonly-accepted architectures in the brain, namely, a high-frequency neural oscillation driven by an external low-frequency input and two interacting local oscillations with distinct, locally-generated frequencies. We further propose an interconnection structure for brain regions and demonstrate that low-frequency phase synchrony can integrate high-frequency activities regulated by local PAC and control the direction of information flow across distant regions., Comment: 6 pages, 5 figures

Many complex networks depend upon biological entities for their preservation. Such entities, from human cognition to evolution, must first encode and then replicate those networks under marked resource constraints. Networks that survive are those that are amenable to constrained encoding, or, in other words, are compressible. But how compressible is a network? And what features make one network more compressible than another? Here we answer these questions by modeling networks as information sources before compressing them using rate-distortion theory. Each network yields a unique rate-distortion curve, which specifies the minimal amount of information that remains at a given scale of description. A natural definition then emerges for the compressibility of a network: the amount of information that can be removed via compression, averaged across all scales. Analyzing an array of real and model networks, we demonstrate that compressibility increases with two common network properties: transitivity (or clustering) and degree heterogeneity. These results indicate that hierarchical organization -- which is characterized by modular structure and heterogeneous degrees -- facilitates compression in complex networks. Generally, our framework sheds light on the interplay between a network's structure and its capacity to be compressed, enabling investigations into the role of compression in shaping real-world networks., Comment: 18 pages, 10 figures

Brain-computer interfaces (BCIs) constitute a promising tool for communication and control. However, mastering non-invasive closed-loop systems remains a learned skill that is difficult to develop for a non-negligible proportion of users. The involved learning process induces neural changes associated with a brain network reorganization that remains poorly understood. To address this inter-subject variability, we adopted a multilayer approach to integrate brain network properties from electroencephalographic (EEG) and magnetoencephalographic (MEG) data resulting from a four-session BCI training program followed by a group of healthy subjects. Our method gives access to the contribution of each layer to multilayer network that tends to be equal with time. We show that regardless the chosen modality, a progressive increase in the integration of somatosensory areas in the alpha band was paralleled by a decrease of the integration of visual processing and working memory areas in the beta band. Notably, only brain network properties in multilayer network correlated with future BCI scores in the alpha2 band: positively in somatosensory and decision-making related areas and negatively in associative areas. Our findings cast new light on neural processes underlying BCI training. Integrating multimodal brain network properties provides new information that correlates with behavioral performance and could be considered as a potential marker of BCI learning., Comment: This is the version of the article before editing, as submitted by an author to the Journal of Neural Engineering. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online athttp://iopscience.iop.org/article/10.1088/1741-2552/abef39

Philosophers of science have long postulated how collective scientific knowledge grows. Empirical validation has been challenging due to limitations in collecting and systematizing large historical records. Here, we capitalize on the largest online encyclopedia to formulate knowledge as growing networks of articles and their hyperlinked inter-relations. We demonstrate that concept networks grow not by expanding from their core but rather by creating and filling knowledge gaps, a process which produces discoveries that are more frequently awarded Nobel prizes than others. Moreover, we operationalize paradigms as network modules to reveal a temporal signature in structural stability across scientific subjects. In a network formulation of scientific discovery, data-driven conditions underlying breakthroughs depend just as much on identifying uncharted gaps as on advancing solutions within scientific communities.

White matter microstructure underpins cognition and function in the human brain through the facilitation of neuronal communication, and the non-invasive characterization of this structure remains a research frontier in the neuroscience community. Efforts to assess white matter microstructure, however, are hampered by the sheer amount of information needed for characterization. Current techniques within neuroimaging deal with this problem by representing white matter features with single scalars that are often not easy to interpret. Here, we address these issues by introducing tools from materials science for the characterization of white matter microstructure. We investigate structure on a mesoscopic scale by analyzing its homogeneity and determining which regions of the brain are structurally homogeneous, or "crystalline" in the context of materials science. We find that crystallinity provides novel information and varies across the brain along interpretable lines of anatomical difference, with highest homogeneity in regions adjacent to the corpus callosum. Furthermore, crystallinity is highly reliable, and yet also varies between individuals, making it a potentially useful biomarker to examine individual differences in white matter along various dimensions. We also parcellate white matter into "crystal grains," or contiguous sets of voxels of high structural similarity, and find overlap with other white matter parcellations. Finally, we characterize the shapes of local white matter signatures through another tool from materials science - bond-orientational order parameters - and locate fiber crossings and fascicles with this new technique. Our results provide new means of assessing white matter microstructure on multiple length scales, and open new avenues of future inquiry. We hope that this work fosters future productive dialogue between the fields of soft matter and neuroscience.

The psychosis spectrum is associated with structural dysconnectivity concentrated in transmodal association cortex. However, understanding of this pathophysiology has been limited by an exclusive focus on the direct connections to a region. Using Network Control Theory, we measured variation in both direct and indirect structural connections to a region to gain new insights into the pathophysiology of the psychosis spectrum. We used psychosis symptom data and structural connectivity in 1,068 youths aged 8 to 22 years from the Philadelphia Neurodevelopmental Cohort. Applying a Network Control Theory metric called average controllability, we estimated each brain region's capacity to leverage its direct and indirect structural connections to control linear brain dynamics. Next, using non-linear regression, we determined the accuracy with which average controllability could predict negative and positive psychosis spectrum symptoms in out-of-sample testing. We also compared prediction performance for average controllability versus strength, which indexes only direct connections to a region. Finally, we assessed how the prediction performance for psychosis spectrum symptoms varied over the functional hierarchy spanning unimodal to transmodal cortex. Average controllability outperformed strength at predicting positive psychosis spectrum symptoms, demonstrating that indexing indirect structural connections to a region improved prediction performance. Critically, improved prediction was concentrated in association cortex for average controllability, whereas prediction performance for strength was uniform across the cortex, suggesting that indexing indirect connections is crucial in association cortex. Examining inter-individual variation in direct and indirect structural connections to association cortex is crucial for accurate prediction of positive psychosis spectrum symptoms.

In networks of coupled oscillators, it is of interest to understand how interaction topology affects synchronization. Many studies have gained key insights into this question by studying the classic Kuramoto oscillator model on static networks. However, new questions arise when network structure is time-varying or when the oscillator system is multistable, the latter of which can occur when an inertial term is added to the Kuramoto model. While the consequences of evolving topology and multistability on collective behavior have been examined separately, real-world systems such as gene regulatory networks and the brain can exhibit these properties simultaneously. How does the rewiring of network connectivity affect synchronization in systems with multistability, where different paths of network evolution may differentially impact system dynamics? To address this question, we study the effects of time-evolving network topology on coupled Kuramoto oscillators with inertia. We show that hysteretic synchronization behavior occurs when the network density of coupled inertial oscillators is slowly varied as the dynamics evolve. Moreover, we find that certain fixed-density rewiring schemes induce significant changes to the level of global synchrony, and that these changes remain after the network returns to its initial configuration and are robust to a wide range of network perturbations. Our findings suggest that the specific progression of network topology, in addition to its initial or final static structure, can play a considerable role in modulating the collective behavior of systems evolving on complex networks.

When presented with information of any type, from music to language to mathematics, the human mind subconsciously arranges it into a network. A network puts pieces of information like musical notes, syllables or mathematical concepts into context by linking them together. These networks help our minds organize information and anticipate what is coming. Here we present two questions about network building. 1) Can humans more easily learn some types of networks than others? 2) Do humans find some links between ideas more surprising than others? The answer to both questions is "Yes," and we explain why. The findings provide much-needed insight into the ways that humans learn about the networked world around them. Moreover, the study paves the way for future efforts seeking to optimize how information is presented to accelerate human learning., Comment: 11 pages, 3 figures

Primates display a marked ability to learn habits in uncertain and dynamic environments. The associated perceptions and actions of such habits engage distributed neural circuits. Yet, precisely how such circuits support the computations necessary for habit learning remain far from understood. Here we construct a formal theory of network energetics to account for how changes in brain state produce changes in sequential behavior. We exercise the theory in the context of multi-unit recordings spanning the caudate nucleus, prefrontal cortex, and frontal eyefields of female macaque monkeys engaged in 60-180 sessions of a free scan task that induces motor habits. The theory relies on the determination of effective connectivity between recording channels, and on the stipulation that a brain state is taken to be the trial-specific firing rate across those channels. The theory then predicts how much energy will be required to transition from one state into another, given the constraint that activity can spread solely through effective connections. Consistent with the theory's predictions, we observed smaller energy requirements for transitions between more similar and more complex trial saccade patterns, and for sessions characterized by less entropic selection of saccade patterns. Using a virtual lesioning approach, we demonstrate the resilience of the observed relationships between minimum control energy and behavior to significant disruptions in the inferred effective connectivity. Our theoretically principled approach to the study of habit learning paves the way for future efforts examining how behavior arises from changing patterns of activity in distributed neural circuitry., Comment: Main Text: 17 pages and 6 figures; Supplement Text: 9 pages and 8 figures

Searching for biomarkers has been a chief pursuit of the field of psychiatry. Toward this end, studies have catalogued candidate resting-state biomarkers in nearly all forms of mental disorder. However, it is becoming increasingly clear that these biomarkers lack specificity, limiting their capacity to yield clinical impact. We discuss three avenues of research that are overcoming this limitation: (i) the adoption of transdiagnostic research designs, which involve studying and explicitly comparing multiple disorders from distinct diagnostic axes of psychiatry; (ii) dimensional models of psychopathology that map the full spectrum of symptomatology and that cut across traditional disorder boundaries; and (iii) modeling individuals' unique functional connectomes throughout development. We provide a framework for tying these subfields together that draws on tools from machine learning and network science., Comment: Keywords: Psychiatry; biomarker; transdiagnostic; dimensional psychopathology; normative modeling; fingerprint; network science

Throughout life, we might seek a calling, companions, skills, entertainment, truth, self-knowledge, beauty, and edification. The practice of curiosity can be viewed as an extended and open-ended search for valuable information with hidden identity and location in a complex space of interconnected information. Despite its importance, curiosity has been challenging to computationally model because the practice of curiosity often flourishes without specific goals, external reward, or immediate feedback. Here, we show how network science, statistical physics, and philosophy can be integrated into an approach that coheres with and expands the psychological taxonomies of specific-diversive and perceptual-epistemic curiosity. Using this interdisciplinary approach, we distill functional modes of curious information seeking as searching movements in information space. The kinesthetic model of curiosity offers a vibrant counterpart to the deliberative predictions of model-based reinforcement learning. In doing so, this model unearths new computational opportunities for identifying what makes curiosity curious.

Complex systems thinking is applied to a wide variety of domains, from neuroscience to computer science and economics. The wide variety of implementations has resulted in two key challenges: the progenation of many domain-specific strategies that are seldom revisited or questioned, and the siloing of ideas within a domain due to inconsistency of complex systems language. In this work we offer basic, domain-agnostic language in order to advance towards a more cohesive vocabulary. We use this language to evaluate each step of the complex systems analysis pipeline, beginning with the system and data collected, then moving through different mathematical formalisms for encoding the observed data (i.e. graphs, simplicial complexes, and hypergraphs), and relevant computational methods for each formalism. At each step we consider different types of \emph{dependencies}; these are properties of the system that describe how the existence of one relation among the parts of a system may influence the existence of another relation. We discuss how dependencies may arise and how they may alter interpretation of results or the entirety of the analysis pipeline. We close with two real-world examples using coauthorship data and email communications data that illustrate how the system under study, the dependencies therein, the research question, and choice of mathematical representation influence the results. We hope this work can serve as an opportunity of reflection for experienced complexity scientists, as well as an introductory resource for new researchers., Comment: 50 pages, 18 figures

Living systems break detailed balance at small scales, consuming energy and producing entropy in the environment in order to perform molecular and cellular functions. However, it remains unclear how broken detailed balance manifests at macroscopic scales, and how such dynamics support higher-order biological functions. Here we present a framework to quantify broken detailed balance by measuring entropy production in macroscopic systems. We apply our method to the human brain, an organ whose immense metabolic consumption drives a diverse range of cognitive functions. Using whole-brain imaging data, we demonstrate that the brain nearly obeys detailed balance when at rest, but strongly breaks detailed balance when performing physically and cognitively demanding tasks. Using a dynamic Ising model, we show that these large-scale violations of detailed balance can emerge from fine-scale asymmetries in the interactions between elements, a known feature of neural systems. Together, these results suggest that violations of detailed balance are vital for cognition, and provide a general tool for quantifying entropy production in macroscopic systems., Comment: 19 pages, 15 figures

The ability to store and manipulate information is a hallmark of computational systems. Whereas computers are carefully engineered to represent and perform mathematical operations on structured data, neurobiological systems perform analogous functions despite flexible organization and unstructured sensory input. Recent efforts have made progress in modeling the representation and recall of information in neural systems. However, precisely how neural systems learn to modify these representations remains far from understood. Here we demonstrate that a recurrent neural network (RNN) can learn to modify its representation of complex information using only examples, and we explain the associated learning mechanism with new theory. Specifically, we drive an RNN with examples of translated, linearly transformed, or pre-bifurcated time series from a chaotic Lorenz system, alongside an additional control signal that changes value for each example. By training the network to replicate the Lorenz inputs, it learns to autonomously evolve about a Lorenz-shaped manifold. Additionally, it learns to continuously interpolate and extrapolate the translation, transformation, and bifurcation of this representation far beyond the training data by changing the control signal. Finally, we provide a mechanism for how these computations are learned, and demonstrate that a single network can simultaneously learn multiple computations. Together, our results provide a simple but powerful mechanism by which an RNN can learn to manipulate internal representations of complex information, allowing for the principled study and precise design of RNNs., Comment: 7 main text figures, 3 supplementary figures

The nature of yield in amorphous materials under stress has yet to be fully elucidated. In particular, understanding how microscopic rearrangement gives rise to macroscopic structural and rheological signatures in disordered systems is vital for the prediction and characterization of yield and for the study of how memory is stored in disordered materials. Here, we investigate the evolution of local structural homogeneity on an individual particle level in amorphous jammed two-dimensional systems under oscillatory shear, and relate this evolution to rearrangement, memory, and macroscale rheological measurements. We identify a new structural metric, crystalline shielding, that is predictive of rearrangement propensity and the structural volatility of individual particles under shear. We use this metric to identify localized regions of the system in which the material's memory of its preparation is preserved. Our results contribute to a growing understanding of how local structure relates to dynamic response and memory in disordered systems.

Over one third of the estimated 3 million people with epilepsy in the US are medication resistant. Responsive neurostimulation from chronically implanted electrodes provides a promising treatment option and alternative to resective surgery. However, determining personalized optimal stimulation parameters, including when and where to intervene to guarantee a positive patient outcome, is a major open challenge. Network neuroscience and control theory offer useful tools that may guide improvements in parameter selection for control of anomalous neural activity. Here we use a novel method to characterize dynamic controllability across consecutive effective connectivity (EC) networks based on regularized partial correlations between implanted electrodes during the onset, propagation, and termination phases of thirty-four seizures. We estimate regularized partial correlation adjacency matrices from one-second time windows of intracranial electrocorticography recordings using the Graphical Least Absolute Shrinkage and Selection Operator (GLASSO). Average and modal controllability metrics calculated from each resulting EC network track the time-varying controllability of the brain on an evolving landscape of conditionally dependent network interactions. We show that average controllability increases throughout a seizure and is negatively correlated with modal controllability throughout. Furthermore, our results support the hypothesis that the energy required to drive the brain to a seizure-free state from an ictal state is smallest during seizure onset; yet, we find that applying control energy at electrodes in the seizure onset zone may not always be energetically favorable. Our work suggests that a low-complexity model of time-evolving controllability may offer new insights for developing and improving control strategies targeting seizure suppression., Comment: Main text is 7 pages and contains 4 figures, supplement is 17 pages and contains 9 figures and one table

Our ability to manipulate the behavior of complex networks depends on the design of efficient control algorithms and, critically, on the availability of an accurate and tractable model of the network dynamics. While the design of control algorithms for network systems has seen notable advances in the past few years, knowledge of the network dynamics is a ubiquitous assumption that is difficult to satisfy in practice, especially when the network topology is large and, possibly, time-varying. In this paper we overcome this limitation, and develop a data-driven framework to control a complex dynamical network optimally and without requiring any knowledge of the network dynamics. Our optimal controls are constructed using a finite set of experimental data, where the unknown complex network is stimulated with arbitrary and possibly random inputs. In addition to optimality, we show that our data-driven formulas enjoy favorable computational and numerical properties even compared to their model-based counterpart. Although our controls are provably correct for networks with linear dynamics, we also characterize their performance against noisy experimental data and in the presence of nonlinear dynamics, as they arise when mitigating cascading failures in power-grid networks and when manipulating neural activity in brain networks.

Cortical regions without direct neuronal connections have been observed to exhibit synchronized dynamics. A recent empirical study has further revealed that such regions that share more common neighbors are more likely to behave coherently. To analytically investigate the underlying mechanisms, we consider that a set of n oscillators, which have no direct connections, are linked through m intermediate oscillators (called mediators), forming a complete bipartite network structure. Modeling the oscillators by the Kuramoto-Sakaguchi model, we rigorously prove that mediated remote synchronization, i.e., synchronization between those n oscillators that are not directly connected, becomes more robust as the number of mediators increases. Simulations are also carried out to show that our theoretical findings can be applied to other general and complex networks.

How an individual's unique brain connectivity determines that individual's cognition, behavior, and risk for pathology is a fundamental question in basic and clinical neuroscience. In seeking answers, many have turned to machine learning, with some noting the particular promise of deep neural networks in modelling complex non-linear functions. However, it is not clear that complex functions actually exist between brain connectivity and behavior, and thus if deep neural networks necessarily outperform simpler linear models, or if their results would be interpretable. Here we show that, across 52 subject measures of cognition and behavior, deep neural networks fit to each brain region's connectivity outperform linear regression, particularly for the brain's connector hubs--regions with diverse brain connectivity--whereas the two approaches perform similarly when fit to brain systems. Critically, averaging deep neural network predictions across brain regions results in the most accurate predictions, demonstrating the ability of deep neural networks to easily model the various functions that exists between regional brain connectivity and behavior, carving the brain at its joints. Finally, we shine light into the black box of deep neural networks using multislice network models. We determined that the relationship between connector hubs and behavior is best captured by modular deep neural networks. Our results demonstrate that both simple and complex relationships exist between brain connectivity and behavior, and that deep neural networks can fit both. Moreover, deep neural networks are particularly powerful when they are first fit to the various functions of a system independently and then combined. Finally, deep neural networks are interpretable when their architectures are structurally characterized using multislice network models.

Recent advances in computational models of signal propagation and routing in the human brain have underscored the critical role of white matter structure. A complementary approach has utilized the framework of network control theory to better understand how white matter constrains the manner in which a region or set of regions can direct or control the activity of other regions. Despite the potential for both of these approaches to enhance our understanding of the role of network structure in brain function, little work has sought to understand the relations between them. Here, we seek to explicitly bridge computational models of communication and principles of network control in a conceptual review of the current literature. By drawing comparisons between communication and control models in terms of the level of abstraction, the dynamical complexity, the dependence on network attributes, and the interplay of multiple spatiotemporal scales, we highlight the convergence of and distinctions between the two frameworks. Based on the understanding of the intertwined nature of communication and control in human brain networks, this work provides an integrative perspective for the field and outlines exciting directions for future work., Comment: 23 pages, 4 figures

The human brain displays rich communication dynamics that are thought to be particularly well-reflected in its marked community structure. Yet, the precise relationship between community structure in structural brain networks and the communication dynamics that can emerge therefrom is not well-understood. In addition to offering insight into the structure-function relationship of networked systems, such an understanding is a critical step towards the ability to manipulate the brain's large-scale dynamical activity in a targeted manner. We investigate the role of community structure in the controllability of structural brain networks. At the region level, we find that certain network measures of community structure are sometimes statistically correlated with measures of linear controllability. However, we then demonstrate that this relationship depends on the distribution of network edge weights. We highlight the complexity of the relationship between community structure and controllability by performing numerical simulations using canonical graph models with varying mesoscale architectures and edge weight distributions. Finally, we demonstrate that weighted subgraph centrality, a measure rooted in the graph spectrum, and which captures higher-order graph architecture, is a stronger and more consistent predictor of controllability. Our study contributes to an understanding of how the brain's diverse mesoscale structure supports transient communication dynamics., Comment: 32 pages, 7 figures in main text, and 22 pages, 11 figures in supplement

A fundamental understanding of synchronized behavior in multi-agent systems can be acquired by studying analytically tractable Kuramoto models. However, such models typically diverge from many real systems whose dynamics evolve under non-negligible resource constraints. Here we construct a system of coupled Kuramoto oscillators that consume or produce resources as a function of their oscillation frequency. At high coupling, we observe strongly synchronized dynamics, whereas at low coupling we observe independent oscillator dynamics, as expected from standard Kuramoto models. For intermediate coupling, which typically induces a partially synchronized state, we empirically observe that (and theoretically explain why) the system can exist in either (i) a state in which the order parameter oscillates in time, or (ii) a state in which multiple synchronization states are simultaneously stable. Whether (i) or (ii) occurs depends upon whether the oscillators consume or produce resources, respectively. Relevant for systems as varied as coupled neurons and social groups, our study lays important groundwork for future efforts to develop quantitative predictions of synchronized dynamics for systems embedded in environments marked by sparse resources., Comment: 8 pages, 5 figure, submitted to Physical Review E

At the macroscale, the brain operates as a network of interconnected neuronal populations, which display rhythmic dynamics that support interareal communication. Understanding how stimulation of a particular brain area impacts such concerted activity is important for gaining basic insights into brain function and for developing neuromodulation as a therapeutic tool. However, it remains difficult to predict the downstream effects of focal stimulation. Specifically, little is known about how the collective oscillatory regime of network activity may affect the outcomes of regional perturbations on cooperative dynamics. Here, we combine connectome data and biophysical modeling to begin filling these gaps. By tuning parameters that control the collective dynamics of the network, we identify distinct states of simulated brain activity, and investigate how the distributed effects of stimulation manifest in different states. When baseline oscillations are weak, the stimulated area exhibits enhanced power and frequency, and due to network interactions, nearby regions develop phase locked activity in the excited frequency band. Importantly, we find that focal stimulation also causes more distributed modifications to network coherence at regions' baseline oscillation frequencies, and that these effects are better predicted by functional rather than structural connectivity. In contrast, when the network operates in a regime of stronger endogenous oscillations, stimulation causes only slight shifts in power and frequency, and network averaged changes in coherence are more homogenous across the choice of the stimulated area. In sum, this work builds upon and extends previous computational studies investigating the impacts of stimulation, and highlights that both the stimulation site, and, crucially, the regime of brain network dynamics, can influence the network wide responses to local perturbations., Comment: Main Text: 32 pages, 10 figures; Supplement: 14 pages, 10 figures

In systems neuroscience, most models posit that brain regions communicate information under constraints of efficiency. Yet, evidence for efficient communication in structural brain networks characterized by hierarchical organization and highly connected hubs remains sparse. The principle of efficient coding proposes that the brain transmits maximal information in a metabolically economical or compressed form to improve future behavior. To determine how structural connectivity supports efficient coding, we develop a theory specifying minimum rates of message transmission between brain regions to achieve an expected fidelity, and we test five predictions from the theory based on random walk communication dynamics. In doing so, we introduce the metric of compression efficiency, which quantifies the trade-off between lossy compression and transmission fidelity in structural networks. In a large sample of youth (n = 1,042; age 8-23 years), we analyze structural networks derived from diffusion weighted imaging and metabolic expenditure operationalized using cerebral blood flow. We show that structural networks strike compression efficiency trade-offs consistent with theoretical predictions. We find that compression efficiency prioritizes fidelity with development, heightens when metabolic resources and myelination guide communication, explains advantages of hierarchical organization, links higher input fidelity to disproportionate areal expansion, and shows that hubs integrate information by lossy compression. Lastly, compression efficiency is predictive of behavior--beyond a conventional metric--for cognitive domains including executive function, memory, complex reasoning, and social cognition. Our findings elucidate how macroscale connectivity supports efficient coding, and serve to foreground communication processes that utilize random walk dynamics constrained by network connectivity.

Like many scientific disciplines, neuroscience has increasingly attempted to confront pervasive gender imbalances within the field. While much of the conversation has centered around publishing and conference participation, recent research in other fields has called attention to the prevalence of gender bias in citation practices. Because of the downstream effects that citations can have on visibility and career advancement, understanding and eliminating gender bias in citation practices is vital for addressing inequity in a scientific community. In this study, we sought to determine whether there is evidence of gender bias in the citation practices of neuroscientists. Using data from five top neuroscience journals, we find that reference lists tend to include more papers with men as first and last author than would be expected if gender were not a factor in referencing. Importantly, we show that this overcitation of men and undercitation of women is driven largely by the citation practices of men, and is increasing over time as the field becomes more diverse. We develop a co-authorship network to assess homophily in researchers' social networks, and we find that men tend to overcite men even when their social networks are representative. We discuss possible mechanisms and consider how individual researchers might address these findings in their own practices.

Whether listening to overlapping conversations in a crowded room or recording the simultaneous electrical activity of millions of neurons, the natural world abounds with sparse measurements of complex overlapping signals that arise from dynamical processes. While tools that separate mixed signals into linear sources have proven necessary and useful, the underlying equational forms of most natural signals are unknown and nonlinear. Hence, there is a need for a framework that is general enough to extract sources without knowledge of their generating equations, and flexible enough to accommodate nonlinear, even chaotic, sources. Here we provide such a framework, where the sources are chaotic trajectories from independently evolving dynamical systems. We consider the mixture signal as the sum of two chaotic trajectories, and propose a supervised learning scheme that extracts the chaotic trajectories from their mixture. Specifically, we recruit a complex dynamical system as an intermediate processor that is constantly driven by the mixture. We then obtain the separated chaotic trajectories based on this intermediate system by training the proper output functions. To demonstrate the generalizability of this framework \emph{in silico}, we employ a tank of water as the intermediate system, and show its success in separating two-part mixtures of various chaotic trajectories. Finally, we relate the underlying mechanism of this method to the state-observer problem. This relation provides a quantitative theory that explains the performance of our method, such as why separation is difficult when two source signals are trajectories from the same chaotic system.

The extraction of brain functioning features is a crucial step in the definition of brain-computer interfaces (BCIs). In the last decade, functional connectivity (FC) estimators have been increasingly explored based on their ability to capture synchronization between multivariate brain signals. However, the underlying neurophysiological mechanisms and the extent to which they can improve performance in BCI-related tasks, is still poorly understood. To address this gap in knowledge, we considered a group of 20 healthy subjects during an EEG-based hand motor imagery (MI) task. We studied two well-established FC estimators, i.e. spectral- and imaginary-coherence, and investigated how they were modulated by the MI task. We characterized the resulting FC networks by extracting the strength of connectivity of each EEG sensor and compared the discriminant power with respect to standard power spectrum features. At the group level, results showed that while spectral-coherence based network features were increasing the controlateral motor area, those based on imaginary-coherence were decreasing. We demonstrated that this opposite, but complementary, behavior was respectively determined by the increase in amplitude and phase synchronization between the brain signals. At the individual level, we proved that including these network connectivity features in the classification of MI mental states led to an overall improvement in accuracy. Taken together, our results provide fresh insights into the oscillatory mechanisms subserving brain network changes during MI and offer new perspectives to improve BCI performance.

Network neuroscience represents the brain as a collection of regions and inter-regional connections. Given its ability to formalize systems-level models, network neuroscience has generated unique explanations of neural function and behavior. The mechanistic status of these explanations and how they can contribute to and fit within the field of neuroscience as a whole has received careful treatment from philosophers. However, these philosophical contributions have not yet reached many neuroscientists. Here we complement formal philosophical efforts by providing an applied perspective from and for neuroscientists. We discuss the mechanistic status of the explanations offered by network neuroscience and how they contribute to, enhance, and interdigitate with other types of explanations in neuroscience. In doing so, we rely on philosophical work concerning the role of causality, scale, and mechanisms in scientific explanations. In particular, we make the distinction between an explanation and the evidence supporting that explanation, and we argue for a scale-free nature of mechanistic explanations. In the course of these discussions, we hope to provide a useful applied framework in which network neuroscience explanations can be exercised across scales and combined with other fields of neuroscience to gain deeper insights into the brain and behavior., Comment: This article is part a forthcoming Topics in Cognitive Science Special Issue: "Levels of Explanation in Cognitive Science: From Molecules to Culture," Matteo Colombo and Markus Knauff (Topic Editors)

Humans communicate, receive, and store information using sequences of items -- from words in a sentence or notes in music to abstract concepts in lectures and books. The networks formed by these items (nodes) and the sequential transitions between them (edges) encode important structural features of human communication and knowledge. But how do humans learn the networks of probabilistic transitions that underlie sequences of items? Moreover, what do people's internal maps of these networks look like? Here, we introduce graph learning, a growing and interdisciplinary field focused on studying how humans learn and represent networks in the world around them. We begin by describing established results from statistical learning showing that humans are adept at detecting differences in the transition probabilities between items in a sequence. We next present recent experiments that directly control for differences in transition probabilities, demonstrating that human behavior also depends critically on the abstract network structure of transitions. Finally, we present computational models that researchers have proposed to explain the effects of network structure on human behavior and cognition. Throughout, we highlight a number of exciting open questions in the study of graph learning that will require creative insights from cognitive scientists and network scientists alike., Comment: 9 pages, 6 figures