Catalog footprint

What is connected

95works
53topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

95 published item(s)

preprint2026arXiv

A disease-spread model on hypergraphs with distinct droplet and aerosol transmission modes

We examine the spread of an infectious disease, such as one that is caused by a respiratory virus, with two distinct modes of transmission. To do this, we consider a susceptible--infected--susceptible (SIS) disease on a hypergraph, which allows us to incorporate the effects of both dyadic (i.e., pairwise) and polyadic (i.e., group) interactions on disease propagation. This disease can spread either via large droplets through direct social contacts, which we associate with edges (i.e., hyperedges of size 2), or via infected aerosols in the environment through hyperedges of size at least 3 (i.e., polyadic interactions). We derive mean-field approximations of our model for two types of hypergraphs, and we obtain threshold conditions that characterize whether the disease dies out or becomes endemic. Additionally, we numerically simulate our model and a mean-field approximation of it to examine the impact of various factors, such as hyperedge size (when the size is uniform), hyperedge-size distribution (when the sizes are nonuniform), and hyperedge-recovery rates (when the sizes are nonuniform) on the disease dynamics.

preprint2026arXiv

Graph energy as a measure of community detectability in networks

A key challenge in network science is the detection of communities, which are sets of nodes in a network that are densely connected internally but sparsely connected to the rest of the network. A fundamental result in community detection is the existence of a nontrivial threshold for community detectability on sparse graphs that are generated by the planted partition model (PPM). Below this so-called ``detectability limit'', no community-detection method can perform better than random chance. Spectral methods for community detection fail before this detectability limit because the eigenvalues corresponding to the eigenvectors that are relevant for community detection can be absorbed by the bulk of the spectrum. One can bypass the detectability problem by using special matrices, like the non-backtracking matrix, but this requires one to consider higher-dimensional matrices. In this paper, we show that the difference in graph energy between a PPM and an Erdős--Rényi (ER) network has a distinct transition at the detectability threshold even for the adjacency matrices of the underlying networks. The graph energy is based on the full spectrum of an adjacency matrix, so our result suggests that standard graph matrices still allow one to separate the parameter regions with detectable and undetectable communities.

preprint2022arXiv

A Majority-Vote Model On Multiplex Networks with Community Structure

We investigate a majority-vote model on two-layer multiplex networks with community structure. In our majority-vote model, the edges on each layer encode one type of social relationship and an individual changes their opinion based on the majority opinions of their neighbors in each layer. To capture the fact that different relationships often have different levels of importance, we introduce a layer-preference parameter, which determines the probability of a node to adopt an opinion when the node's neighborhoods on the two layers have different majority opinions. We construct our networks so that each node is a member of one community on each layer, and we consider situations in which nodes tend to have more connections with nodes from the same community than with nodes from different communities. We study the influence of the layer-preference parameter, the intralayer communities, and interlayer membership correlation on the steady-state behavior of our model using both direct numerical simulations and a mean-field approximation. We find three different types of steady-state behavior: a fully-mixed state, consensus states, and polarized states. We demonstrate that a stronger interlayer community correlation makes polarized steady states reachable for wider ranges of the other model parameters. We also show that different values of the layer-preference parameter result in qualitatively different phase diagrams for the mean opinions at steady states.

preprint2022arXiv

Analysis of Spatial and Spatiotemporal Anomalies Using Persistent Homology: Case Studies with COVID-19 Data

We develop a method for analyzing spatial and spatiotemporal anomalies in geospatial data using topological data analysis (TDA). To do this, we use persistent homology (PH), which allows one to algorithmically detect geometric voids in a data set and quantify the persistence of such voids. We construct an efficient filtered simplicial complex (FSC) such that the voids in our FSC are in one-to-one correspondence with the anomalies. Our approach goes beyond simply identifying anomalies; it also encodes information about the relationships between anomalies. We use vineyards, which one can interpret as time-varying persistence diagrams (which are an approach for visualizing PH), to track how the locations of the anomalies change with time. We conduct two case studies using spatially heterogeneous COVID-19 data. First, we examine vaccination rates in New York City by zip code at a single point in time. Second, we study a year-long data set of COVID-19 case rates in neighborhoods of the city of Los Angeles.

preprint2022arXiv

Lonely individuals process the world in idiosyncratic ways

Loneliness is detrimental to well-being and is often accompanied by self-reported feelings of not being understood by others. What contributes to such feelings in lonely people? We used functional magnetic resonance imaging (fMRI) of 66 participants to unobtrusively measure the relative alignment of people's mental processing of naturalistic stimuli and tested whether or not lonely people actually process the world in idiosyncratic ways. We found evidence for such idiosyncrasy: lonely individuals' neural responses were dissimilar to their peers, particularly in regions of the default-mode network in which similar responses have been associated with shared perspectives and subjective understanding. These relationships persisted when controlling for demographic similarities, objective social isolation, and participants' friendships with each other. Our findings suggest the possibility that being surrounded by people who see the world differently from oneself, even if one is friends with them, may be a risk factor for loneliness.

preprint2022arXiv

Mixed Logit Models and Network Formation

The study of network formation is pervasive in economics, sociology, and many other fields. In this paper, we model network formation as a `choice' that is made by nodes in a network to connect to other nodes. We study these `choices' using discrete-choice models, in which an agent chooses between two or more discrete alternatives. We employ the `repeated-choice' (RC) model to study network formation. We argue that the RC model overcomes important limitations of the multinomial logit (MNL) model, which gives one framework for studying network formation, and that it is well-suited to study network formation. We also illustrate how to use the RC model to accurately study network formation using both synthetic and real-world networks. Using edge-independent synthetic networks, we also compare the performance of the MNL model and the RC model. We find that the RC model estimates the data-generation process of our synthetic networks more accurately than the MNL model. In a patent citation network, which forms sequentially, we present a case study of a qualitatively interesting scenario -- the fact that new patents are more likely to cite older, more cited, and similar patents -- for which employing the RC model yields interesting insights.

preprint2021arXiv

Counterparty Credit Limits: The Impact of a Risk-Mitigation Measure on Everyday Trading

A counterparty credit limit (CCL) is a limit that is imposed by a financial institution to cap its maximum possible exposure to a specified counterparty. CCLs help institutions to mitigate counterparty credit risk via selective diversification of their exposures. In this paper, we analyze how CCLs impact the prices that institutions pay for their trades during everyday trading. We study a high-quality data set from a large electronic trading platform in the foreign exchange spot market, which enables institutions to apply CCLs. We find empirically that CCLs had little impact on the vast majority of trades in this data. We also study the impact of CCLs using a new model of trading. By simulating our model with different underlying CCL networks, we highlight that CCLs can have a major impact in some situations.

preprint2021arXiv

Tie-decay networks in continuous time and eigenvector-based centralities

Network theory is a useful framework for studying interconnected systems of interacting entities. Many networked systems evolve continuously in time, but most existing methods for the analysis of time-dependent networks rely on discrete or discretized time. In this paper, we propose an approach for studying networks that evolve in continuous time by distinguishing between \emph{interactions}, which we model as discrete contacts, and \emph{ties}, which encode the strengths of relationships as functions of time. To illustrate our tie-decay network formalism, we adapt the well-known PageRank centrality score to our tie-decay framework in a mathematically tractable and computationally efficient way. We apply this framework to a synthetic example and then use it to study a network of retweets during the 2012 National Health Service controversy in the United Kingdom. Our work also provides guidance for similar generalizations of other tools from network theory to continuous-time networks with tie decay, including for applications to streaming data.

preprint2020arXiv

A Model for the Influence of Media on the Ideology of Content in Online Social Networks

Many people rely on online social networks as sources of news and information, and the spread of media content with ideologies across the political spectrum influences online discussions and impacts actions offline. To examine the impact of media in online social networks, we generalize bounded-confidence models of opinion dynamics by incorporating media accounts as influencers in a network. We quantify partisanship of content with a continuous parameter on an interval, and we formulate higher-dimensional generalizations to incorporate content quality and increasingly nuanced political positions. We simulate our model with one and two ideological dimensions, and we use the results of our simulations to quantify the "entrainment" of content from non-media accounts to the ideologies of media accounts in a network. We maximize media impact in a social network by tuning the number of media accounts that promote the content and the number of followers of the accounts. Using numerical computations, we find that the entrainment of the ideology of content spread by non-media accounts to media ideology depends on a network's structural features, including its size, the mean number of followers of its nodes, and the receptiveness of its nodes to different opinions. We then introduce content quality --- a key novel contribution of our work --- into our model. We incorporate multiple media sources with ideological biases and quality-level estimates that we draw from real media sources and demonstrate that our model can produce distinct communities ("echo chambers") that are polarized in both ideology and quality. Our model provides a step toward understanding content quality and ideology in spreading dynamics, with ramifications for how to mitigate the spread of undesired content and promote the spread of desired content.

preprint2020arXiv

A unified framework for equivalences in social networks

A key concern in network analysis is the study of social positions and roles of actors in a network. The notion of "position" refers to an equivalence class of nodes that have similar ties to other nodes, whereas a "role" is an equivalence class of compound relations that connect the same pairs of nodes. An open question in network science is whether it is possible to simultaneously perform role and positional analysis. Motivated by the principle of functoriality in category theory we propose a new method that allows to tie role and positional analysis together. We illustrate our methods on two well-studied data sets in network science.

preprint2020arXiv

Connecting the Dots: Discovering the "Shape" of Data

Scientists use a mathematical subject called 'topology' to study the shapes of objects. An important part of topology is counting the numbers of pieces and holes in objects, and people use this information to group objects into different types. For example, a doughnut has the same number of holes and the same number of pieces as a teacup with one handle, but it is different from a ball. In studies that resemble activities like "connect the dots", scientists use ideas from topology to study the shape of data. Data can take many possible forms: a picture made of dots, a large collection of numbers from a scientific experiment, or something else. The approach in these studies is called 'topological data analysis', and it has been used to study the branching structures of veins in leaves, how people vote in elections, flight patterns in models of bird flocking, and more. Scientists can take data on the way veins branch on leaves and use topological data analysis to divide the leaves into different groups and discover patterns that may otherwise be hard to find.

preprint2020arXiv

Disease Detectives: Using Mathematics to Forecast the Spread of Infectious Diseases

The COVID-19 pandemic has led to significant changes in how people are currently living their lives. To determine how to best reduce the effects of the pandemic and start reopening societies, governments have drawn insights from mathematical models of the spread of infectious diseases. In this article, we give an introduction to a family of mathematical models (called "compartmental models") and discuss how the results of analyzing these models influence government policies and human behavior, such as encouraging mask wearing and physical distancing to help slow the spread of the disease.

preprint2020arXiv

Fitting In and Breaking Up: A Nonlinear Version of Coevolving Voter Models

We investigate a nonlinear version of coevolving voter models, in which node states and network structure update as a coupled stochastic dynamical process. Most prior work on coevolving voter models has focused on linear update rules with fixed and homogeneous rewiring and adopting probabilities. By contrast, in our nonlinear version, the probability that a node rewires or adopts is a function of how well it "fits in" within its neighborhood. To explore this idea, we incorporate a parameter $σ$ that represents the fraction of neighbors of an updating node that share its opinion state. In an update, with probability $σ^q$ (for some nonlinearity parameter $q$), the updating node rewires; with complementary probability $1-σ^q$, the updating node adopts a new opinion state. We study this mechanism using three rewiring schemes: after an updating node deletes a discordant edge, it then either (1) "rewires-to-random" by choosing a new neighbor in a random process; (2) "rewires-to-same" by choosing a new neighbor in a random process from nodes that share its state; or (3) "rewires-to-none" by not rewiring at all (akin to "unfriending" on social media). We compare our nonlinear coevolving model to several existing linear models, and we find in our model that initial network topology plays a larger role in the dynamics and the choice of rewiring mechanism plays a smaller role. A particularly interesting feature of our model is that, under certain conditions, the opinion state that is held initially by a minority of the nodes can effectively spread to almost every node in a network if the minority nodes view themselves as the majority. In light of this observation, we relate our results to recent work on the majority illusion in social networks.

preprint2020arXiv

Forecasting elections using compartmental models of infection

Forecasting elections -- a challenging, high-stakes problem -- is the subject of much uncertainty, subjectivity, and media scrutiny. To shed light on this process, we develop a method for forecasting elections from the perspective of dynamical systems. Our model borrows ideas from epidemiology, and we use polling data from United States elections to determine its parameters. Surprisingly, our general model performs as well as popular forecasters for the 2012 and 2016 U.S. races for president, senators, and governors. Although contagion and voting dynamics differ, our work suggests a valuable approach to elucidate how elections are related across states. It also illustrates the effect of accounting for uncertainty in different ways, provides an example of data-driven forecasting using dynamical systems, and suggests avenues for future research on political elections. We conclude with our forecasts for the senatorial and gubernatorial races on 6~November 2018, which we posted on 5 November 2018.

preprint2020arXiv

Inference of Edge Correlations in Multilayer Networks

Many recent developments in network analysis have focused on multilayer networks, which one can use to encode time-dependent interactions, multiple types of interactions, and other complications that arise in complex systems. Like their monolayer counterparts, multilayer networks in applications often have mesoscale features, such as community structure. A prominent type of method for inferring such structures is the employment of multilayer stochastic block models (SBMs). A common (but {potentially} inadequate) assumption of these models is the sampling of edges in different layers independently, conditioned on the community labels of the nodes. In this paper, we relax this assumption of independence by incorporating edge correlations into an SBM-like model. We derive maximum-likelihood estimates of the key parameters of our model, and we propose a measure of layer correlation that reflects the similarity between connectivity patterns in different layers. Finally, we explain how to use correlated models for edge "prediction" (i.e., inference) in multilayer networks. By taking into account edge correlations, prediction accuracy improves both in synthetic networks and in a temporal network of shoppers who are connected to previously-purchased grocery products.

preprint2020arXiv

Migration Networks: Applications of Network Analysis to Macroscale Migration Patterns

An emerging area of research is the study of macroscale migration patterns as a network of nodes that represent places (e.g., countries, cities, and rural areas) and edges that encode migration ties that connect those places. In this chapter, we first review advances in the study of migration networks and recent work that has employed network analysis to examine such networks at different geographical scales. In our discussion, we focus in particular on global scale migration networks. We then propose ways to leverage network analysis in concert with digital technologies and online geolocated data to examine the structure and dynamics of migration networks. The implementation of such approaches for studying migration networks faces many challenges, including ethical ones, methodological ones, socio-technological ones (e.g., data availability and reuse), and research reproducibility. We detail these challenges, and we then consider possible ways of linking digital geolocated data to administrative and survey data as a way of harnessing new technologies to construct increasingly realistic migration networks (e.g., using multiplex networks). We also briefly discuss new methods (e.g., multilayer network analysis) in network analysis and adjacent fields (e.g., machine learning) that can help advance understanding of macroscale patterns of migration.

preprint2020arXiv

Random walks and diffusion on networks

Random walks are ubiquitous in the sciences, and they are interesting from both theoretical and practical perspectives. They are one of the most fundamental types of stochastic processes; can be used to model numerous phenomena, including diffusion, interactions, and opinions among humans and animals; and can be used to extract information about important entities or dense groups of entities in a network. Random walks have been studied for many decades on both regular lattices and (especially in the last couple of decades) on networks with a variety of structures. In the present article, we survey the theory and applications of random walks on networks, restricting ourselves to simple cases of single and non-adaptive random walkers. We distinguish three main types of random walks: discrete-time random walks, node-centric continuous-time random walks, and edge-centric continuous-time random walks. We first briefly survey random walks on a line, and then we consider random walks on various types of networks. We extensively discuss applications of random walks, including ranking of nodes (e.g., PageRank), community detection, respondent-driven sampling, and opinion models such as voter models.

preprint2020arXiv

Social Network Analysis for Social Neuroscientists

Although social neuroscience is concerned with understanding how the brain interacts with its social environment, prevailing research in the field has primarily considered the human brain in isolation, deprived of its rich social context. Emerging work in social neuroscience that leverages tools from network analysis has begun to pursue this issue, advancing knowledge of how the human brain influences and is influenced by the structures of its social environment. In this paper, we provide an overview of key theory and methods in network analysis (especially for social systems) as an introduction for social neuroscientists who are interested in relating individual cognition to the structures of an individual's social environments. We also highlight some exciting new work as examples of how to productively use these tools to investigate questions of relevance to social neuroscientists. We include tutorials to help with practical implementation of the concepts that we discuss. We conclude by highlighting a broad range of exciting research opportunities for social neuroscientists who are interested in using network analysis to study social systems.

preprint2020arXiv

Spatial Applications of Topological Data Analysis: Cities, Snowflakes, Random Structures, and Spiders Spinning Under the Influence

Spatial networks are ubiquitous in social, geographical, physical, and biological applications. To understand the large-scale structure of networks, it is important to develop methods that allow one to directly probe the effects of space on structure and dynamics. Historically, algebraic topology has provided one framework for rigorously and quantitatively describing the global structure of a space, and recent advances in topological data analysis (TDA) have given scholars a new lens for analyzing network data. In this paper, we study a variety of spatial networks -- including both synthetic and natural ones -- using novel topological methods that we recently developed for analyzing spatial networks. We demonstrate that our methods are able to capture meaningful quantities, with specifics that depend on context, in spatial networks and thereby provide useful insights into the structure of those networks, including a novel approach for characterizing them based on their topological structures. We illustrate these ideas with examples of synthetic networks and dynamics on them, street networks in cities, snowflakes, and webs spun by spiders under the influence of various psychotropic substances.

preprint2020arXiv

Spatial Strength Centrality and the Effect of Spatial Embeddings on Network Architecture

For many networks, it is useful to think of their nodes as being embedded in a latent space, and such embeddings can affect the probabilities for nodes to be adjacent to each other. In this paper, we extend existing models of synthetic networks to spatial network models by first embedding nodes in Euclidean space and then modifying the models so that progressively longer edges occur with progressively smaller probabilities. We start by extending a geographical fitness model by employing Gaussian-distributed fitnesses, and we then develop spatial versions of preferential attachment and configuration models. We define a notion of "spatial strength centrality" to help characterize how strongly a spatial embedding affects network structure, and we examine spatial strength centrality on a variety of real and synthetic networks.

preprint2020arXiv

Topological Data Analysis of Task-Based fMRI Data from Experiments on Schizophrenia

We use methods from computational algebraic topology to study functional brain networks, in which nodes represent brain regions and weighted edges encode the similarity of fMRI time series from each region. With these tools, which allow one to characterize topological invariants such as loops in high-dimensional data, we are able to gain understanding into low-dimensional structures in networks in a way that complements traditional approaches that are based on pairwise interactions. In the present paper, we use persistent homology to analyze networks that we construct from task-based fMRI data from schizophrenia patients, healthy controls, and healthy siblings of schizophrenia patients. We thereby explore the persistence of topological structures such as loops at different scales in these networks. We use persistence landscapes and persistence images to create output summaries from our persistent-homology calculations, and we study the persistence landscapes and images using $k$-means clustering and community detection. Based on our analysis of persistence landscapes, we find that the members of the sibling cohort have topological features (specifically, their 1-dimensional loops) that are distinct from the other two cohorts. From the persistence images, we are able to distinguish all three subject groups and to determine the brain regions in the loops (with four or more edges) that allow us to make these distinctions.

preprint2020arXiv

Tunable Eigenvector-Based Centralities for Multiplex and Temporal Networks

Characterizing the importances (i.e., centralities) of nodes in social, biological, and technological networks is a core topic in both network science and data science. We present a linear-algebraic framework that generalizes eigenvector-based centralities, including PageRank and hub/authority scores, to provide a common framework for two popular classes of multilayer networks: multiplex networks (which have layers that encode different types of relationships) and temporal networks (in which the relationships change over time). Our approach involves the study of joint, marginal, and conditional "supracentralities" that one can calculate from the dominant eigenvector of a supracentrality matrix [Taylor et al., 2017], which couples centrality matrices that are associated with individual network layers. We extend this prior work (which was restricted to temporal networks with layers that are coupled by adjacent-in-time coupling) by allowing the layers to be coupled through a (possibly asymmetric) interlayer-adjacency matrix $\tilde{\bf A}$, where the entry $\tilde{A}_{tt'} \geq 0$ encodes the coupling between layers $t$ and $t'$. Our framework provides a unifying foundation for centrality analysis of multiplex and temporal networks; it also illustrates a complicated dependency of the supracentralities on the topology and weights of interlayer coupling. By scaling $\tilde{\bf A}$ by an interlayer-coupling strength $ω\ge0$ and developing a singular perturbation theory for the limits of weak ($ω\to0^+$) and strong coupling ($ω\to\infty$), we also reveal an interesting dependence of supracentralities on the dominant left and right eigenvectors of $\tilde{\bf A}$.

preprint2019arXiv

A Framework for the Construction of Generative Models for Mesoscale Structure in Multilayer Networks

Multilayer networks allow one to represent diverse and coupled connectivity patterns --- e.g., time-dependence, multiple subsystems, or both --- that arise in many applications and which are difficult or awkward to incorporate into standard network representations. In the study of multilayer networks, it is important to investigate mesoscale (i.e., intermediate-scale) structures, such as dense sets of nodes known as communities, to discover network features that are not apparent at the microscale or the macroscale. The ill-defined nature of mesoscale structure and its ubiquity in empirical networks make it crucial to develop generative models that can produce the features that one encounters in empirical networks. Key purposes of such generative models include generating synthetic networks with empirical properties of interest, benchmarking mesoscale-detection methods and algorithms, and inferring structure in empirical multilayer networks. In this paper, we introduce a framework for the construction of generative models for mesoscale structures in multilayer networks. Our framework provides a standardized set of generative models, together with an associated set of principles from which they are derived, for studies of mesoscale structures in multilayer networks. It unifies and generalizes many existing models for mesoscale structures in fully-ordered (e.g., temporal) and unordered (e.g., multiplex) multilayer networks. One can also use it to construct generative models for mesoscale structures in partially-ordered multilayer networks (e.g., networks that are both temporal and multiplex). Our framework has the ability to produce many features of empirical multilayer networks, and it explicitly incorporates a user-specified dependency structure between layers.

preprint2016arXiv

A Local Perspective on Community Structure in Multilayer Networks

The analysis of multilayer networks is among the most active areas of network science, and there are now several methods to detect dense "communities" of nodes in multilayer networks. One way to define a community is as a set of nodes that trap a diffusion-like dynamical process (usually a random walk) for a long time. In this view, communities are sets of nodes that create bottlenecks to the spreading of a dynamical process on a network. We analyze the local behavior of different random walks on multiplex networks (which are multilayer networks in which different layers correspond to different types of edges) and show that they have very different bottlenecks that hence correspond to rather different notions of what it means for a set of nodes to be a good community. This has direct implications for the behavior of community-detection methods that are based on these random walks.

preprint2016arXiv

Detection of Core-Periphery Structure in Networks Using Spectral Methods and Geodesic Paths

We introduce several novel and computationally efficient methods for detecting "core--periphery structure" in networks. Core--periphery structure is a type of mesoscale structure that includes densely-connected core vertices and sparsely-connected peripheral vertices. Core vertices tend to be well-connected both among themselves and to peripheral vertices, which tend not to be well-connected to other vertices. Our first method, which is based on transportation in networks, aggregates information from many geodesic paths in a network and yields a score for each vertex that reflects the likelihood that a vertex is a core vertex. Our second method is based on a low-rank approximation of a network's adjacency matrix, which can often be expressed as a tensor-product matrix. Our third approach uses the bottom eigenvector of the random-walk Laplacian to infer a coreness score and a classification into core and peripheral vertices. We also design an objective function to (1) help classify vertices into core or peripheral vertices and (2) provide a goodness-of-fit criterion for classifications into core versus peripheral vertices. To examine the performance of our methods, we apply our algorithms to both synthetically-generated networks and a variety of networks constructed from real-world data sets.

preprint2016arXiv

Eigenvector-Based Centrality Measures for Temporal Networks

Numerous centrality measures have been developed to quantify the importances of nodes in time-independent networks, and many of them can be expressed as the leading eigenvector of some matrix. With the increasing availability of network data that changes in time, it is important to extend such eigenvector-based centrality measures to time-dependent networks. In this paper, we introduce a principled generalization of network centrality measures that is valid for any eigenvector-based centrality. We consider a temporal network with N nodes as a sequence of T layers that describe the network during different time windows, and we couple centrality matrices for the layers into a supra-centrality matrix of size NTxNT whose dominant eigenvector gives the centrality of each node i at each time t. We refer to this eigenvector and its components as a joint centrality, as it reflects the importances of both the node i and the time layer t. We also introduce the concepts of marginal and conditional centralities, which facilitate the study of centrality trajectories over time. We find that the strength of coupling between layers is important for determining multiscale properties of centrality, such as localization phenomena and the time scale of centrality changes. In the strong-coupling regime, we derive expressions for time-averaged centralities, which are given by the zeroth-order terms of a singular perturbation expansion. We also study first-order terms to obtain first-order-mover scores, which concisely describe the magnitude of nodes' centrality changes over time. As examples, we apply our method to three empirical temporal networks: the United States Ph.D. exchange in mathematics, costarring relationships among top-billed actors during the Golden Age of Hollywood, and citations of decisions from the United States Supreme Court.

preprint2016arXiv

Heterogeneity of Global and Local Connectivity in Spatial Network Structures of World Migration

We examine world migration as a social-spatial network of countries connected via movements of people. We assess how multilateral migratory relationships at global, regional, and local scales coexist ("glocalization"), divide ("polarization"), or form an interconnected global system ("globalization"). To do this, we decompose the world migration network (WMN) into communities---sets of countries with denser than expected migration connections---and characterize their pattern of local (i.e., intracommunity) and global (i.e., intercommunity) connectivity. We distinguish community signatures---"cave", "biregional", and "bridging"---with distinct migration patterns, spatial network structures, temporal dynamics, and underlying antecedents. Cave communities are tightly-knit, enduring structures that tend to channel local migration between contiguous countries; biregional communities are likely to merge migration between two distinct geographic regions (e.g., North Africa and Europe); and bridging communities have hub-and-spoke structures that tend to emerge dynamically from globe-spanning movements. We find that world migration is neither globally interconnected nor reproduces the geographic boundaries as drawn on a world map but involves a heterogeneous interplay of global and local tendencies in different network regions. We discuss the implications of our results for the understating of variability in today's transnational mobility patterns and migration opportunities across the globe.

preprint2016arXiv

Quasi-Centralized Limit Order Books

A quasi-centralized limit order book (QCLOB) is a limit order book (LOB) in which financial institutions can only access the trading opportunities offered by counterparties with whom they possess sufficient bilateral credit. We perform an empirical analysis of a recent, high-quality data set from a large electronic trading platform that utilizes QCLOBs to facilitate trade. We find many significant differences between our results and those widely reported for other LOBs. We also uncover a remarkable empirical universality: although the distributions describing order flow and market state vary considerably across days, a simple, linear rescaling causes them to collapse onto a single curve. Motivated by this finding, we propose a semi-parametric model of order flow and market state in a QCLOB on a single trading day. Our model provides similar performance to that of parametric curve-fitting techniques, while being simpler to compute and faster to implement.

preprint2016arXiv

Scattering of Waves by Impurities in Precompressed Granular Chains

We study scattering of waves by impurities in strongly precompressed granular chains. We explore the linear scattering of plane waves and identify a closed-form expression for the reflection and transmission coefficients for the scattering of the waves from both a single impurity and a double impurity. For single-impurity chains, we show that, within the transmission band of the host granular chain, high-frequency waves are strongly attenuated (such that the transmission coefficient vanishes as the wavenumber $k\rightarrow \pmπ$), whereas low-frequency waves are well-transmitted through the impurity. For double-impurity chains, we identify a resonance --- enabling full transmission at a particular frequency --- in a manner that is analogous to the Ramsauer--Townsend (RT) resonance from quantum physics. We also demonstrate that one can tune the frequency of the RT resonance to any value in the pass band of the host chain. We corroborate our theoretical predictions both numerically and experimentally, and we directly observe complete transmission for frequencies close to the RT resonance frequency. Finally, we show how this RT resonance can lead to the existence of reflectionless modes even in granular chains (including disordered ones) with multiple double impurities.

preprint2016arXiv

The Topological "Shape" of Brexit

Persistent homology is a method from computational algebraic topology that can be used to study the "shape" of data. We illustrate two filtrations --- the weight rank clique filtration and the Vietoris--Rips (VR) filtration --- that are commonly used in persistent homology, and we apply these filtrations to a pair of data sets that are both related to the 2016 European Union "Brexit" referendum in the United Kingdom. These examples consider a topical situation and give useful illustrations of the strengths and weaknesses of these methods.

preprint2015arXiv

Dynamical Systems on Networks: A Tutorial

We give a tutorial for the study of dynamical systems on networks. We focus especially on "simple" situations that are tractable analytically, because they can be very insightful and provide useful springboards for the study of more complicated scenarios. We briefly motivate why examining dynamical systems on networks is interesting and important, and we then give several fascinating examples and discuss some theoretical results. We also briefly discuss dynamical systems on dynamical (i.e., time-dependent) networks, overview software implementations, and give an outlook on the field.

preprint2015arXiv

Estimating inter-event time distributions from finite observation periods in communication networks

A diverse variety of processes --- including recurrent disease episodes, neuron firing, and communication patterns among humans --- can be described using inter-event time (IET) distributions. Many such processes are ongoing, although event sequences are only available during a finite observation window. Because the observation time window is more likely to begin or end during long IETs than during short ones, the analysis of such data is susceptible to a bias induced by the finite observation period. In this paper, we illustrate how this length bias is born and how it can be corrected without assuming any particular shape for the IET distribution. To do this, we model event sequences using stationary renewal processes, and we formulate simple heuristics for determining the severity of the bias. To illustrate our results, we focus on the example of empirical communication networks, which are temporal networks that are constructed from communication events. The IET distributions of such systems guide efforts to build models of human behavior, and the variance of IETs is very important for estimating the spreading rate of information in networks of temporal interactions. We analyze several well-known data sets from the literature, and we find that the resulting bias can lead to systematic underestimates of the variance in the IET distributions and that correcting for the bias can lead to qualitatively different results for the tails of the IET distributions.

preprint2015arXiv

Extraction of Force-Chain Network Architecture in Granular Materials Using Community Detection

Force chains form heterogeneous physical structures that can constrain the mechanical stability and acoustic transmission of granular media. However, despite their relevance for predicting bulk properties of materials, there is no agreement on a quantitative description of force chains. Consequently, it is difficult to compare the force-chain structures in different materials or experimental conditions. To address this challenge, we treat granular materials as spatially-embedded networks in which the nodes (particles) are connected by weighted edges that represent contact forces. We use techniques from community detection, which is a type of clustering, to find sets of closely connected particles. By using a geographical null model that is constrained by the particles' contact network, we extract chain-like structures that are reminiscent of force chains. We propose three diagnostics to measure these chain-like structures, and we demonstrate the utility of these diagnostics for identifying and characterizing classes of force-chain network architectures in various materials. To illustrate our methods, we describe how force-chain architecture depends on pressure for two very different types of packings: (1) ones derived from laboratory experiments and (2) ones derived from idealized, numerically-generated frictionless packings. By resolving individual force chains, we quantify statistical properties of force-chain shape and strength, which are potentially crucial diagnostics of bulk properties (including material stability). These methods facilitate quantitative comparisons between different particulate systems, regardless of whether they are measured experimentally or numerically.

preprint2015arXiv

Heterogeneous, Weakly Coupled Map Lattices

Coupled map lattices (CMLs) are often used to study emergent phenomena in nature. It is typically assumed (unrealistically) that each component is described by the same map, and it is important to relax this assumption. In this paper, we characterize periodic orbits and the laminar regime of type-I intermittency in heterogeneous weakly coupled map lattices (HWCMLs). We show that the period of a cycle in an HWCML is preserved for arbitrarily small coupling strengths even when an associated uncoupled oscillator would experience a period-doubling cascade. Our results characterize periodic orbits both near and far from saddle--node bifurcations, and we thereby provide a key step for examining the bifurcation structure of heterogeneous CMLs.

preprint2015arXiv

MuxViz: A Tool for Multilayer Analysis and Visualization of Networks

Multilayer relationships among entities and information about entities must be accompanied by the means to analyze, visualize, and obtain insights from such data. We present open-source software (muxViz) that contains a collection of algorithms for the analysis of multilayer networks, which are an important way to represent a large variety of complex systems throughout science and engineering. We demonstrate the ability of muxViz to analyze and interactively visualize multilayer data using empirical genetic, neuronal, and transportation networks. Our software is available at https://github.com/manlius/muxViz.

preprint2015arXiv

Null Models for Community Detection in Spatially-Embedded, Temporal Networks

In the study of networks, it is often insightful to use algorithms to determine mesoscale features such as "community structure", in which densely connected sets of nodes constitute "communities" that have sparse connections to other communities. The most popular way of detecting communities algorithmically is to optimize the quality function known as modularity. When optimizing modularity, one compares the actual connections in a (static or time-dependent) network to the connections obtained from a random-graph ensemble that acts as a null model. The communities are then the sets of nodes that are connected to each other densely relative to what is expected from the null model. Clearly, the process of community detection depends fundamentally on the choice of null model, so it is important to develop and analyze novel null models that take into account appropriate features of the system under study. In this paper, we investigate the effects of using null models that take incorporate spatial information, and we propose a novel null model based on the radiation model of population spread. We also develop novel synthetic spatial benchmark networks in which the connections between entities are based on distance or flux between nodes, and we compare the performance of both static and time-dependent radiation null models to the standard ("Newman-Girvan") null model for modularity optimization and a recently-proposed gravity null model. In our comparisons, we use both the above synthetic benchmarks and time-dependent correlation networks that we construct using countrywide dengue fever incidence data for Peru. We also evaluate a recently-proposed correlation null model, which was developed specifically for correlation networks that are constructed from time series, on the epidemic-correlation data.

preprint2015arXiv

Numerical Methods for the Computation of the Confluent and Gauss Hypergeometric Functions

The two most commonly used hypergeometric functions are the confluent hypergeometric function and the Gauss hypergeometric function. We review the available techniques for accurate, fast, and reliable computation of these two hypergeometric functions in different parameter and variable regimes. The methods that we investigate include Taylor and asymptotic series computations, Gauss-Jacobi quadrature, numerical solution of differential equations, recurrence relations, and others. We discuss the results of numerical experiments used to determine the best methods, in practice, for each parameter and variable regime considered. We provide 'roadmaps' with our recommendation for which methods should be used in each situation.

preprint2015arXiv

Structure of Triadic Relations in Multiplex Networks

Recent advances in the study of networked systems have highlighted that our interconnected world is composed of networks that are coupled to each other through different "layers" that each represent one of many possible subsystems or types of interactions. Nevertheless, it is traditional to aggregate multilayer networks into a single weighted network in order to take advantage of existing tools. This is admittedly convenient, but it is also extremely problematic, as important information can be lost as a result. It is therefore important to develop multilayer generalizations of network concepts. In this paper, we analyze triadic relations and generalize the idea of transitivity to multiplex networks. By focusing on triadic relations, which yield the simplest type of transitivity, we generalize the concept and computation of clustering coefficients to multiplex networks. We show how the layered structure of such networks introduces a new degree of freedom that has a fundamental effect on transitivity. We compute multiplex clustering coefficients for several real multiplex networks and illustrate why one must take great care when generalizing standard network concepts to multiplex networks. We also derive analytical expressions for our clustering coefficients for ensemble averages of networks in a family of random multiplex networks. Our analysis illustrates that social networks have a strong tendency to promote redundancy by closing triads at every layer and that they thereby have a different type of multiplex transitivity from transportation networks, which do not exhibit such a tendency. These insights are invisible if one only studies aggregated networks.

preprint2015arXiv

Superdiffusive Transport and Energy Localization in Disordered Granular Crystals

We study the spreading of initially localized excitations in 1D disordered granular crystals. We thereby investigate localization phenomena in strongly nonlinear systems, which we demonstrate to be fundamentally different from localization in linear and weakly nonlinear systems. We compare wave dynamics in chains with 3 different types of disorder: an uncorrelated (Anderson-like) disorder and 2 types of correlated disorders (random dimer arrangements), and for 2 types of initial conditions: displacement excitations and velocity excitations. For strongly precompressed chains, the dynamics depend strongly on the initial condition. For displacement excitations, the long-time behavior of the second moment $\tilde{m}_2$ has oscillations that depend on the type of disorder, with a complex trend that differs markedly from a power law and which is particularly evident for an Anderson disorder. For velocity excitations, we find a scaling $\tilde{m}_2\sim t^γ$ (for a constant $γ$) for all 3 types of disorder. For weakly precompressed (strongly nonlinear) chains, $\tilde{m}_2$ and the inverse participation ratio $P^{-1}$ satisfy $\tilde{m}_2\sim t^γ$ and $P^{-1}\sim t^{-η}$, and the dynamics is superdiffusive for all examined cases. When precompression is strong, the IPR decreases slowly for all 3 types of disorder, and we observe a partial localization around the core and the leading edge of the wave. For an Anderson disorder, displacement perturbations lead to localization of energy primarily in the core, and velocity perturbations cause the energy to be divided between the core and the leading edge. This localization phenomenon does not occur in the sonic-vacuum regime, which yields the surprising result that the energy is no longer contained in strongly nonlinear waves but instead is spread across many sites. In this regime, the exponents are very similar in all cases.

preprint2015arXiv

The Long Memory of Order Flow in the Foreign Exchange Spot Market

We study the long memory of order flow for each of three liquid currency pairs on a large electronic trading platform in the foreign exchange (FX) spot market. Due to the extremely high levels of market activity on the platform, and in contrast to existing empirical studies of other markets, our data enables us to perform statistically stable estimation without needing to aggregate data from different trading days. We find strong evidence of long memory, with a Hurst exponent of approximately 0.7, for each of the three currency pairs and on each trading day in our sample. We repeat our calculations using data that spans different trading days, and we find no significant differences in our results. We test and reject the hypothesis that the apparent long memory of order flow is an artifact caused by structural breaks, in favour of the alternative hypothesis of true long memory. We therefore conclude that the long memory of order flow in the FX spot market is a robust empirical property that persists across daily boundaries.

preprint2015arXiv

Topological data analysis of contagion maps for examining spreading processes on networks

Social and biological contagions are influenced by the spatial embeddedness of networks. Historically, many epidemics spread as a wave across part of the Earth's surface; however, in modern contagions long-range edges -- for example, due to airline transportation or communication media -- allow clusters of a contagion to appear in distant locations. Here we study the spread of contagions on networks through a methodology grounded in topological data analysis and nonlinear dimension reduction. We construct "contagion maps" that use multiple contagions on a network to map the nodes as a point cloud. By analyzing the topology, geometry, and dimensionality of manifold structure in such point clouds, we reveal insights to aid in the modeling, forecast, and control of spreading processes. Our approach highlights contagion maps also as a viable tool for inferring low-dimensional structure in networks.

preprint2014arXiv

A Simple Generative Model of Collective Online Behaviour

Human activities increasingly take place in online environments, providing novel opportunities for relating individual behaviours to population-level outcomes. In this paper, we introduce a simple generative model for the collective behaviour of millions of social networking site users who are deciding between different software applications. Our model incorporates two distinct components: one is associated with recent decisions of users, and the other reflects the cumulative popularity of each application. Importantly, although various combinations of the two mechanisms yield long-time behaviour that is consistent with data, the only models that reproduce the observed temporal dynamics are those that strongly emphasize the recent popularity of applications over their cumulative popularity. This demonstrates---even when using purely observational data without experimental design---that temporal data-driven modelling can effectively distinguish between competing microscopic mechanisms, allowing us to uncover new aspects of collective online behaviour.

preprint2014arXiv

Convergence Time Towards Periodic Orbits in Discrete Dynamical Systems

We investigate the convergence towards periodic orbits in discrete dynamical systems. We examine the probability that a randomly chosen point converges to a particular neighborhood of a periodic orbit in a fixed number of iterations, and we use linearized equations to examine the evolution near that neighborhood. The underlying idea is that points of stable periodic orbit are associated with intervals. We state and prove a theorem that details what regions of phase space are mapped into these intervals (once they are known) and how many iterations are required to get there. We also construct algorithms that allow our theoretical results to be implemented successfully in practice.

preprint2014arXiv

Cross-Linked Structure of Network Evolution

We study the temporal co-variation of network co-evolution via the cross-link structure of networks, for which we take advantage of the formalism of hypergraphs to map cross-link structures back to network nodes. We investigate two sets of temporal network data in detail. In a network of coupled nonlinear oscillators, hyperedges that consist of network edges with temporally co-varying weights uncover the driving co-evolution patterns of edge weight dynamics both within and between oscillator communities. In the human brain, networks that represent temporal changes in brain activity during learning exhibit early co-evolution that then settles down with practice, and subsequent decreases in hyperedge size are consistent with emergence of an autonomous subgraph whose dynamics no longer depends on other parts of the network. Our results on real and synthetic networks give a poignant demonstration of the ability of cross-link structure to uncover unexpected co-evolution attributes in both real and synthetic dynamical systems. This, in turn, illustrates the utility of analyzing cross-links for investigating the structure of temporal networks.

preprint2014arXiv

Density-based and transport-based core-periphery structures in networks

Networks often possess mesoscale structures, and studying them can yield insights into both structure and function. It is most common to study community structure, but numerous other types of mesoscale structures also exist. In this paper, we examine core-periphery structures based on both density and transport. In such structures, core network components are well-connected both among themselves and to peripheral components, which are not well-connected to anything. We examine core-periphery structures in a wide range of examples of transportation, social, and financial networks---including road networks in large urban areas, a rabbit warren, a dolphin social network, a European interbank network, and a migration network between counties in the United States. We illustrate that a recently developed transport-based notion of node coreness is very useful for characterizing transportation networks. We also generalize this notion to examine core versus peripheral edges, and we show that the resulting diagnostic is also useful for transportation networks. To examine the properties of transportation networks further, we develop a family of generative models of roadlike networks. We illustrate the effect of the dimensionality of the embedding space on transportation networks, and we demonstrate that the correlations between different measures of coreness can be very different for different types of networks.

preprint2014arXiv

Dynamics on Modular Networks with Heterogeneous Correlations

We develop a new ensemble of modular random graphs in which degree-degree correlations can be different in each module and the inter-module connections are defined by the joint degree-degree distribution of nodes for each pair of modules. We present an analytical approach that allows one to analyze several types of binary dynamics operating on such networks, and we illustrate our approach using bond percolation, site percolation, and the Watts threshold model. The new network ensemble generalizes existing models (e.g., the well-known configuration model and LFR networks) by allowing a heterogeneous distribution of degree-degree correlations across modules, which is important for the consideration of nonidentical interacting networks.

preprint2014arXiv

Matchmaker, Matchmaker, Make Me a Match: Migration of Populations via Marriages in the Past

The study of human mobility is both of fundamental importance and of great potential value. For example, it can be leveraged to facilitate efficient city planning and improve prevention strategies when faced with epidemics. The newfound wealth of rich sources of data---including banknote flows, mobile phone records, and transportation data---has led to an explosion of attempts to characterize modern human mobility. Unfortunately, the dearth of comparable historical data makes it much more difficult to study human mobility patterns from the past. In this paper, we present an analysis of long-term human migration, which is important for processes such as urbanization and the spread of ideas. We demonstrate that the data record from Korean family books (called "jokbo") can be used to estimate migration patterns via marriages from the past 750 years. We apply two generative models of long-term human mobility to quantify the relevance of geographical information to human marriage records in the data, and we find that the wide variety in the geographical distributions of the clans poses interesting challenges for the direct application of these models. Using the different geographical distributions of clans, we quantify the "ergodicity" of clans in terms of how widely and uniformly they have spread across Korea, and we compare these results to those obtained using surname data from the Czech Republic. To examine population flow in more detail, we also construct and examine a population-flow network between regions. Based on the correlation between ergodicity and migration in Korea, we identify two different types of migration patterns: diffusive and convective. We expect the analysis of diffusive versus convective effects in population flows to be widely applicable to the study of mobility and migration patterns across different cultures.

preprint2014arXiv

Multilayer Networks

In most natural and engineered systems, a set of entities interact with each other in complicated patterns that can encompass multiple types of relationships, change in time, and include other types of complications. Such systems include multiple subsystems and layers of connectivity, and it is important to take such "multilayer" features into account to try to improve our understanding of complex systems. Consequently, it is necessary to generalize "traditional" network theory by developing (and validating) a framework and associated tools to study multilayer systems in a comprehensive fashion. The origins of such efforts date back several decades and arose in multiple disciplines, and now the study of multilayer networks has become one of the most important directions in network science. In this paper, we discuss the history of multilayer networks (and related concepts) and review the exploding body of work on such networks. To unify the disparate terminology in the large body of recent work, we discuss a general framework for multilayer networks, construct a dictionary of terminology to relate the numerous existing concepts to each other, and provide a thorough discussion that compares, contrasts, and translates between related notions such as multilayer networks, multiplex networks, interdependent networks, networks of networks, and many others. We also survey and discuss existing data sets that can be represented as multilayer networks. We review attempts to generalize single-layer-network diagnostics to multilayer networks. We also discuss the rapidly expanding research on multilayer-network models and notions like community structure, connected components, tensor decompositions, and various types of dynamical processes on multilayer networks. We conclude with a summary and an outlook.

preprint2014arXiv

Think Locally, Act Locally: The Detection of Small, Medium-Sized, and Large Communities in Large Networks

It is common in the study of networks to investigate meso-scale features to try to gain an understanding of network structure and function. For example, numerous algorithms have been developed to try to identify "communities," which are typically construed as sets of nodes with denser connections internally than with the remainder of a network. In this paper, we adopt a complementary perspective that "communities" are associated with bottlenecks of locally-biased dynamical processes that begin at seed sets of nodes, and we employ several different community-identification procedures (using diffusion-based and geodesic-based dynamics) to investigate community quality as a function of community size. Using several empirical and synthetic networks, we identify several distinct scenarios for ``size-resolved community structure'' that can arise in real (and realistic) networks. Depending on which scenario holds, one may or may not be able to successfully identify ``good'' communities in a given network, the manner in which different small communities fit together to form meso-scale network structures can be very different, and processes such as viral propagation and information diffusion can exhibit very different dynamics.In addition, our results suggest that, for many large realistic networks, the output of locally-biased methods that focus on communities that are centered around a given seed node might have better conceptual grounding and greater practical utility than the output of global community-detection methods. They also illustrate subtler structural properties that are important to consider in the development of better benchmark networks to test methods for community detection. [Note: Because of space limitations in the arXiv's abstract field, this is an abridged version of the paper's abstract.]

preprint2013arXiv

A Method Based on Total Variation for Network Modularity Optimization using the MBO Scheme

The study of network structure is pervasive in sociology, biology, computer science, and many other disciplines. One of the most important areas of network science is the algorithmic detection of cohesive groups of nodes called "communities". One popular approach to find communities is to maximize a quality function known as {\em modularity} to achieve some sort of optimal clustering of nodes. In this paper, we interpret the modularity function from a novel perspective: we reformulate modularity optimization as a minimization problem of an energy functional that consists of a total variation term and an $\ell_2$ balance term. By employing numerical techniques from image processing and $\ell_1$ compressive sensing -- such as convex splitting and the Merriman-Bence-Osher (MBO) scheme -- we develop a variational algorithm for the minimization problem. We present our computational results using both synthetic benchmark networks and real data.

preprint2013arXiv

Dark Solitary Waves in a Class of Collisionally Inhomogeneous Bose-Einstein Condensates

We study the structure, stability, and dynamics of dark solitary waves in parabolically trapped, collisionally inhomogeneous Bose-Einstein condensates (BECs) with spatially periodic variations of the scattering length. This collisional inhomogeneity yields a nonlinear lattice, which we tune from a small-amplitude, approximately sinusoidal structure to a periodic sequence of densely spaced spikes. We start by investigating time-independent inhomogeneities, and we subsequently examine the dynamical response when one starts with a collisionally homogeneous BEC and then switches on an inhomogeneity either adiabatically or nonadiabatically. Using Bogoliubov-de Gennes linearization as well as direct numerical simulations of the Gross-Pitaevskii equation, we observe dark solitary waves, which can become unstable through oscillatory or exponential instabilities. We find a critical wavelength of the nonlinear lattice that is comparable to the healing length. Near this value, the fundamental eigenmode responsible for the stability of the dark solitary wave changes its direction of movement as a function of the strength of the nonlinearity. When it increases, it collides with other eigenmodes, leading to oscillatory instabilities; when it decreases, it collides with the origin and becomes imaginary, illustrating that the instability mechanism is fundamentally different in wide-well versus narrow-well lattices. When starting from a collisionally homogeneous setup and switching on inhomogeneities, we find that dark solitary waves are preserved generically for aligned lattices. We briefly examine the time scales for the onset of solitary-wave oscillations in this scenario.

preprint2013arXiv

Decentralized Routing on Spatial Networks with Stochastic Edge Weights

We investigate algorithms to find short paths in spatial networks with stochastic edge weights. Our formulation of the problem of finding short paths differs from traditional formulations because we specifically do not make two of the usual simplifying assumptions: (1) we allow edge weights to be stochastic rather than deterministic; and (2) we do not assume that global knowledge of a network is available. We develop a decentralized routing algorithm that provides en route guidance for travelers on a spatial network with stochastic edge weights without the need to rely on global knowledge about the network. To guide a traveler, our algorithm uses an estimation function that evaluates cumulative arrival probability distributions based on distances between pairs of nodes. The estimation function carries a notion of proximity between nodes and thereby enables routing without global knowledge. In testing our decentralized algorithm, we define a criterion that allows one to discriminate among arrival probability distributions, and we test our algorithm and this criterion using both synthetic and real networks.

preprint2013arXiv

Discriminating Power of Centrality Measures

The calculation of centrality measures is common practice in the study of networks, as they attempt to quantify the importance of individual vertices, edges, or other components. Different centralities attempt to measure importance in different ways. In this paper, we examine a conjecture posed by E. Estrada regarding the ability of several measures to distinguish the vertices of networks. Estrada conjectured that if all vertices of a graph have the same subgraph centrality, then all vertices must also have the same degree, eigenvector, closeness, and betweenness centralities. We provide a counterexample for the latter two centrality measures and propose a revised conjecture.

preprint2013arXiv

Formation of Dominance Relationships via Strategy Updating in an Asymmetric Hawk-Dove Game

We develop a model to describe the development of dominance relations between social animals as they use past experiences to inform future interactions. Using the game-theoretic framework of a Hawk-Dove game with asymmetric resource-holding potentials (RHPs), we derive a simple model that describes the social interactions of animals that compete for resources. We then consider a game-playing strategy in which animals acquire information about their RHP asymmetry from the results of their previous contests and subsequently use their asymmetry assessment to inform their behavior in future contests. We examine how directly incorporating the fact that animals have incomplete information in their interactions can lead to outcomes that differ from what would be expected if one considers the situation as a static game in which the animals have perfect information about the asymmetry size. We thereby obtain results that are consistent with observations of dominance-hierarchy formation in social animals. We also discuss how increased interactions between animals can speed up the asymmetry assessment process and how this can explain why aggression can sometimes decrease as the number of interactions between animals increases over time.

preprint2013arXiv

Limit Order Books

Limit order books (LOBs) match buyers and sellers in more than half of the world's financial markets. This survey highlights the insights that have emerged from the wealth of empirical and theoretical studies of LOBs. We examine the findings reported by statistical analyses of historical LOB data and discuss how several LOB models provide insight into certain aspects of the mechanism. We also illustrate that many such models poorly resemble real LOBs and that several well-established empirical facts have yet to be reproduced satisfactorily. Finally, we identify several key unresolved questions about LOBs.

preprint2013arXiv

Mathematical Formulation of Multi-Layer Networks

A network representation is useful for describing the structure of a large variety of complex systems. However, most real and engineered systems have multiple subsystems and layers of connectivity, and the data produced by such systems is very rich. Achieving a deep understanding of such systems necessitates generalizing "traditional" network theory, and the newfound deluge of data now makes it possible to test increasingly general frameworks for the study of networks. In particular, although adjacency matrices are useful to describe traditional single-layer networks, such a representation is insufficient for the analysis and description of multiplex and time-dependent networks. One must therefore develop a more general mathematical framework to cope with the challenges posed by multi-layer complex systems. In this paper, we introduce a tensorial framework to study multi-layer networks, and we discuss the generalization of several important network descriptors and dynamical processes --including degree centrality, clustering coefficients, eigenvector centrality, modularity, Von Neumann entropy, and diffusion-- for this framework. We examine the impact of different choices in constructing these generalizations, and we illustrate how to obtain known results for the special cases of single-layer and multiplex networks. Our tensorial approach will be helpful for tackling pressing problems in multi-layer complex systems, such as inferring who is influencing whom (and by which media) in multichannel social networks and developing routing techniques for multimodal transportation systems.

preprint2013arXiv

Multi-Stage Complex Contagions

The spread of ideas across a social network can be studied using complex contagion models, in which agents are activated by contact with multiple activated neighbors. The investigation of complex contagions can provide crucial insights into social influence and behavior-adoption cascades on networks. In this paper, we introduce a model of a multi-stage complex contagion on networks. Agents at different stages --- which could, for example, represent differing levels of support for a social movement or differing levels of commitment to a certain product or idea --- exert different amounts of influence on their neighbors. We demonstrate that the presence of even one additional stage introduces novel dynamical behavior, including interplay between multiple cascades, that cannot occur in single-stage contagion models. We find that cascades --- and hence collective action --- can be driven not only by high-stage influencers but also by low-stage influencers.

preprint2013arXiv

Noise-Induced Synchronization, Desynchronization, and Clustering in Globally Coupled Nonidentical Oscillators

We study ensembles of globally coupled, nonidentical phase oscillators subject to correlated noise, and we identify several important factors that cause noise and coupling to synchronize or desychronize a system. By introducing noise in various ways, we find a novel estimate for the onset of synchrony of a system in terms of the coupling strength, noise strength, and width of the frequency distribution of its natural oscillations. We also demonstrate that noise alone is sufficient to synchronize nonidentical oscillators. However, this synchrony depends on the first Fourier mode of a phase-sensitivity function, through which we introduce common noise into the system. We show that higher Fourier modes can cause desychronization due to clustering effects, and that this can reinforce clustering caused by different forms of coupling. Finally, we discuss the effects of noise on an ensemble in which antiferromagnetic coupling causes oscillators to form two clusters in the absence of noise.

preprint2013arXiv

Random Walks on Stochastic Temporal Networks

In the study of dynamical processes on networks, there has been intense focus on network structure -- i.e., the arrangement of edges and their associated weights -- but the effects of the temporal patterns of edges remains poorly understood. In this chapter, we develop a mathematical framework for random walks on temporal networks using an approach that provides a compromise between abstract but unrealistic models and data-driven but non-mathematical approaches. To do this, we introduce a stochastic model for temporal networks in which we summarize the temporal and structural organization of a system using a matrix of waiting-time distributions. We show that random walks on stochastic temporal networks can be described exactly by an integro-differential master equation and derive an analytical expression for its asymptotic steady state. We also discuss how our work might be useful to help build centrality measures for temporal networks.

preprint2013arXiv

Robust Detection of Dynamic Community Structure in Networks

We describe techniques for the robust detection of community structure in some classes of time-dependent networks. Specifically, we consider the use of statistical null models for facilitating the principled identification of structural modules in semi-decomposable systems. Null models play an important role both in the optimization of quality functions such as modularity and in the subsequent assessment of the statistical validity of identified community structure. We examine the sensitivity of such methods to model parameters and show how comparisons to null models can help identify system scales. By considering a large number of optimizations, we quantify the variance of network diagnostics over optimizations (`optimization variance') and over randomizations of network structure (`randomization variance'). Because the modularity quality function typically has a large number of nearly-degenerate local optima for networks constructed using real data, we develop a method to construct representative partitions that uses a null model to correct for statistical noise in sets of partitions. To illustrate our results, we employ ensembles of time-dependent networks extracted from both nonlinear oscillators and empirical neuroscience data.

preprint2013arXiv

Solitary Matter Waves in Combined Linear and Nonlinear Potentials: Detection, Stability, and Dynamics

We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous interactions and outline an experimental realization of compensating linear and nonlinear potentials that can yield constant-density solutions. We illustrate how the presence of a step in the nonlinearity coefficient can only be revealed dynamically and consider, in particular, how to reveal it by exploiting the inhomogeneity of the sound speed with a defect-dragging experiment. We conduct computational experiments and observe the spontaneous emergence of dark solitary waves. We use effective-potential theory to perform a detailed analytical investigation of the existence and stability of solitary waves in this setting, and we corroborate these results computationally using a Bogoliubov-de Gennes linear stability analysis. We find that dark solitary waves are unstable for all step widths, whereas bright solitary waves can become stable through a symmetry-breaking bifurcation as one varies the step width. Using phase-plane analysis, we illustrate the scenarios that permit this bifurcation and explore the dynamical outcomes of the interaction between the solitary wave and the step.

preprint2013arXiv

Task-Based Core-Periphery Organisation of Human Brain Dynamics

As a person learns a new skill, distinct synapses, brain regions, and circuits are engaged and change over time. In this paper, we develop methods to examine patterns of correlated activity across a large set of brain regions. Our goal is to identify properties that enable robust learning of a motor skill. We measure brain activity during motor sequencing and characterize network properties based on coherent activity between brain regions. Using recently developed algorithms to detect time-evolving communities, we find that the complex reconfiguration patterns of the brain's putative functional modules that control learning can be described parsimoniously by the combined presence of a relatively stiff temporal core that is composed primarily of sensorimotor and visual regions whose connectivity changes little in time and a flexible temporal periphery that is composed primarily of multimodal association regions whose connectivity changes frequently. The separation between temporal core and periphery changes over the course of training and, importantly, is a good predictor of individual differences in learning success. The core of dynamically stiff regions exhibits dense connectivity, which is consistent with notions of core-periphery organization established previously in social networks. Our results demonstrate that core-periphery organization provides an insightful way to understand how putative functional modules are linked. This, in turn, enables the prediction of fundamental human capacities, including the production of complex goal-directed behavior.

preprint2013arXiv

Two-Particle Circular Billiards Versus Randomly Perturbed One-Particle Circular Billiards

We study a two-particle circular billiard containing two finite-size circular particles that collide elastically with the billiard boundary and with each other. Such a two-particle circular billiard provides a clean example of an "intermittent" system. This billiard system behaves chaotically, but the time scale on which chaos manifests can become arbitrarily long as the sizes of the confined particles become smaller. The finite-time dynamics of this system depends on the relative frequencies of (chaotic) particle-particle collisions versus (integrable) particle-boundary collisions, and investigating these dynamics is computationally intensive because of the long time scales involved. To help improve understanding of such two-particle dynamics, we compare the results of diagnostics used to measure chaotic dynamics for a two-particle circular billiard with those computed for two types of one-particle circular billiards in which a confined particle undergoes random perturbations. Importantly, such one-particle approximations are much less computationally demanding than the original two-particle system, and we expect them to yield reasonable estimates of the extent of chaotic behavior in the two-particle system when the sizes of confined particles are small. Our computations of recurrence-rate coefficients, finite-time Lyapunov exponents, and autocorrelation coefficients support this hypothesis and suggest that studying randomly perturbed one-particle billiards has the potential to yield insights into the aggregate properties of two-particle billiards, which are difficult to investigate directly without enormous computation times (especially when the sizes of the confined particles are small).

preprint2012arXiv

Accuracy of Mean-Field Theory for Dynamics on Real-World Networks

Mean-field analysis is an important tool for understanding dynamics on complex networks. However, surprisingly little attention has been paid to the question of whether mean-field predictions are accurate, and this is particularly true for real-world networks with clustering and modular structure. In this paper, we compare mean-field predictions to numerical simulation results for dynamical processes running on 21 real-world networks and demonstrate that the accuracy of the theory depends not only on the mean degree of the networks but also on the mean first-neighbor degree. We show that mean-field theory can give (unexpectedly) accurate results for certain dynamics on disassortative real-world networks even when the mean degree is as low as 4.

preprint2012arXiv

Dynamic Network Centrality Summarizes Learning in the Human Brain

We study functional activity in the human brain using functional Magnetic Resonance Imaging and recently developed tools from network science. The data arise from the performance of a simple behavioural motor learning task. Unsupervised clustering of subjects with respect to similarity of network activity measured over three days of practice produces significant evidence of `learning', in the sense that subjects typically move between clusters (of subjects whose dynamics are similar) as time progresses. However, the high dimensionality and time-dependent nature of the data makes it difficult to explain which brain regions are driving this distinction. Using network centrality measures that respect the arrow of time, we express the data in an extremely compact form that characterizes the aggregate activity of each brain region in each experiment using a single coefficient, while reproducing information about learning that was discovered using the full data set. This compact summary allows key brain regions contributing to centrality to be visualized and interpreted. We thereby provide a proof of principle for the use of recently proposed dynamic centrality measures on temporal network data in neuroscience.

preprint2012arXiv

Multislice Modularity Optimization in Community Detection and Image Segmentation

Because networks can be used to represent many complex systems, they have attracted considerable attention in physics, computer science, sociology, and many other disciplines. One of the most important areas of network science is the algorithmic detection of cohesive groups (i.e., "communities") of nodes. In this paper, we algorithmically detect communities in social networks and image data by optimizing multislice modularity. A key advantage of modularity optimization is that it does not require prior knowledge of the number or sizes of communities, and it is capable of finding network partitions that are composed of communities of different sizes. By optimizing multislice modularity and subsequently calculating diagnostics on the resulting network partitions, it is thereby possible to obtain information about network structure across multiple system scales. We illustrate this method on data from both social networks and images, and we find that optimization of multislice modularity performs well on these two tasks without the need for extensive problem-specific adaptation. However, improving the computational speed of this method remains a challenging open problem.

preprint2012arXiv

Taxonomies of Networks

The study of networks has grown into a substantial interdisciplinary endeavour that encompasses myriad disciplines in the natural, social, and information sciences. Here we introduce a framework for constructing taxonomies of networks based on their structural similarities. These networks can arise from any of numerous sources: they can be empirical or synthetic, they can arise from multiple realizations of a single process, empirical or synthetic, or they can represent entirely different systems in different disciplines. Since the mesoscopic properties of networks are hypothesized to be important for network function, we base our comparisons on summaries of network community structures. While we use a specific method for uncovering network communities, much of the introduced framework is independent of that choice. After introducing the framework, we apply it to construct a taxonomy for 746 individual networks and demonstrate that our approach usefully identifies similar networks. We also construct taxonomies within individual categories of networks, and in each case we expose non-trivial structure. For example we create taxonomies for similarity networks constructed from both political voting data and financial data. We also construct network taxonomies to compare the social structures of 100 Facebook networks and the growth structures produced by different types of fungi.

preprint2012arXiv

The Extraordinary SVD

The singular value decomposition (SVD) is a popular matrix factorization that has been used widely in applications ever since an efficient algorithm for its computation was developed in the 1970s. In recent years, the SVD has become even more prominent due to a surge in applications and increased computational memory and speed. To illustrate the vitality of the SVD in data analysis, we highlight three of its lesser-known yet fascinating applications: the SVD can be used to characterize political positions of Congressmen, measure the growth rate of crystals in igneous rock, and examine entanglement in quantum computation. We also discuss higher-dimensional generalizations of the SVD, which have become increasingly crucial with the newfound wealth of multidimensional data and have launched new research initiatives in both theoretical and applied mathematics. With its bountiful theory and applications, the SVD is truly extraordinary.

preprint2012arXiv

The Influence of Network Topology on Sound Propagation in Granular Materials

Granular materials, whose features range from the particle scale to the force-chain scale to the bulk scale, are usually modeled as either particulate or continuum materials. In contrast with either of these approaches, network representations are natural for the simultaneous examination of microscopic, mesoscopic, and macroscopic features. In this paper, we treat granular materials as spatially-embedded networks in which the nodes (particles) are connected by weighted edges obtained from contact forces. We test a variety of network measures for their utility in helping to describe sound propagation in granular networks and find that network diagnostics can be used to probe particle-, curve-, domain-, and system-scale structures in granular media. In particular, diagnostics of meso-scale network structure are reproducible across experiments, are correlated with sound propagation in this medium, and can be used to identify potentially interesting size scales. We also demonstrate that the sensitivity of network diagnostics depends on the phase of sound propagation. In the injection phase, the signal propagates systemically, as indicated by correlations with the network diagnostic of global efficiency. In the scattering phase, however, the signal is better predicted by meso-scale community structure, suggesting that the acoustic signal scatters over local geographic neighborhoods. Collectively, our results demonstrate how the force network of a granular system is imprinted on transmitted waves.

preprint2011arXiv

A Mathematical Model for the Dynamics and Synchronization of Cows

We formulate a mathematical model for daily activities of a cow (eating, lying down, and standing) in terms of a piecewise affine dynamical system. We analyze the properties of this bovine dynamical system representing the single animal and develop an exact integrative form as a discrete-time mapping. We then couple multiple cow "oscillators" together to study synchrony and cooperation in cattle herds. We comment on the relevant biology and discuss extensions of our model. With this abstract approach, we not only investigate equations with interesting dynamics but also develop interesting biological predictions. In particular, our model illustrates that it is possible for cows to synchronize \emph{less} when the coupling is increased.

preprint2011arXiv

Dynamic reconfiguration of human brain networks during learning

Human learning is a complex phenomenon requiring flexibility to adapt existing brain function and precision in selecting new neurophysiological activities to drive desired behavior. These two attributes -- flexibility and selection -- must operate over multiple temporal scales as performance of a skill changes from being slow and challenging to being fast and automatic. Such selective adaptability is naturally provided by modular structure, which plays a critical role in evolution, development, and optimal network function. Using functional connectivity measurements of brain activity acquired from initial training through mastery of a simple motor skill, we explore the role of modularity in human learning by identifying dynamic changes of modular organization spanning multiple temporal scales. Our results indicate that flexibility, which we measure by the allegiance of nodes to modules, in one experimental session predicts the relative amount of learning in a future session. We also develop a general statistical framework for the identification of modular architectures in evolving systems, which is broadly applicable to disciplines where network adaptability is crucial to the understanding of system performance.

preprint2011arXiv

Generalized Master Equations for Non-Poisson Dynamics on Networks

The traditional way of studying temporal networks is to aggregate the dynamics of the edges to create a static weighted network. This implicitly assumes that the edges are governed by Poisson processes, which is not typically the case in empirical temporal networks. Consequently, we examine the effects of non-Poisson inter-event statistics on the dynamics of edges, and we apply the concept of a generalized master equation to the study of continuous-time random walks on networks. We show that the equation reduces to the standard rate equations when the underlying process is Poisson and that the stationary solution is determined by an effective transition matrix whose leading eigenvector is easy to calculate. We discuss the implications of our work for dynamical processes on temporal networks and for the construction of network diagnostics that take into account their nontrivial stochastic nature.

preprint2011arXiv

Party Polarization in Congress: A Network Science Approach

We measure polarization in the United States Congress using the network science concept of modularity. Modularity provides a conceptually-clear measure of polarization that reveals both the number of relevant groups and the strength of inter-group divisions without making restrictive assumptions about the structure of the party system or the shape of legislator utilities. We show that party influence on Congressional blocs varies widely throughout history, and that existing measures underestimate polarization in periods with weak party structures. We demonstrate that modularity is a significant predictor of changes in majority party and that turnover is more prevalent at medium levels of modularity. We show that two variables related to modularity, called `divisiveness' and `solidarity,' are significant predictors of reelection success for individual House members. Our results suggest that modularity can serve as an early warning of changing group dynamics, which are reflected only later by changes in party labels.

preprint2011arXiv

Social Structure of Facebook Networks

We study the social structure of Facebook "friendship" networks at one hundred American colleges and universities at a single point in time, and we examine the roles of user attributes - gender, class year, major, high school, and residence - at these institutions. We investigate the influence of common attributes at the dyad level in terms of assortativity coefficients and regression models. We then examine larger-scale groupings by detecting communities algorithmically and comparing them to network partitions based on the user characteristics. We thereby compare the relative importances of different characteristics at different institutions, finding for example that common high school is more important to the social organization of large institutions and that the importance of common major varies significantly between institutions. Our calculations illustrate how microscopic and macroscopic perspectives give complementary insights on the social organization at universities and suggest future studies to investigate such phenomena further.

preprint2011arXiv

Temporal Evolution of Financial Market Correlations

We investigate financial market correlations using random matrix theory and principal component analysis. We use random matrix theory to demonstrate that correlation matrices of asset price changes contain structure that is incompatible with uncorrelated random price changes. We then identify the principal components of these correlation matrices and demonstrate that a small number of components accounts for a large proportion of the variability of the markets that we consider. We then characterize the time-evolving relationships between the different assets by investigating the correlations between the asset price time series and principal components. Using this approach, we uncover notable changes that occurred in financial markets and identify the assets that were significantly affected by these changes. We show in particular that there was an increase in the strength of the relationships between several different markets following the 2007--2008 credit and liquidity crisis.

preprint2010arXiv

Community Structure in the United Nations General Assembly

We study the community structure of networks representing voting on resolutions in the United Nations General Assembly. We construct networks from the voting records of the separate annual sessions between 1946 and 2008 in three different ways: (1) by considering voting similarities as weighted unipartite networks; (2) by considering voting similarities as weighted, signed unipartite networks; and (3) by examining signed bipartite networks in which countries are connected to resolutions. For each formulation, we detect communities by optimizing network modularity using an appropriate null model. We compare and contrast the results that we obtain for these three different network representations. In so doing, we illustrate the need to consider multiple resolution parameters and explore the effectiveness of each network representation for identifying voting groups amidst the large amount of agreement typical in General Assembly votes.

preprint2010arXiv

Community Structure in Time-Dependent, Multiscale, and Multiplex Networks

Network science is an interdisciplinary endeavor, with methods and applications drawn from across the natural, social, and information sciences. A prominent problem in network science is the algorithmic detection of tightly-connected groups of nodes known as communities. We developed a generalized framework of network quality functions that allowed us to study the community structure of arbitrary multislice networks, which are combinations of individual networks coupled through links that connect each node in one network slice to itself in other slices. This framework allows one to study community structure in a very general setting encompassing networks that evolve over time, have multiple types of links (multiplexity), and have multiple scales.

preprint2010arXiv

Comparing Community Structure to Characteristics in Online Collegiate Social Networks

We study the structure of social networks of students by examining the graphs of Facebook "friendships" at five American universities at a single point in time. We investigate each single-institution network's community structure and employ graphical and quantitative tools, including standardized pair-counting methods, to measure the correlations between the network communities and a set of self-identified user characteristics (residence, class year, major, and high school). We review the basic properties and statistics of the pair-counting indices employed and recall, in simplified notation, a useful analytical formula for the z-score of the Rand coefficient. Our study illustrates how to examine different instances of social networks constructed in similar environments, emphasizes the array of social forces that combine to form "communities," and leads to comparative observations about online social lives that can be used to infer comparisons about offline social structures. In our illustration of this methodology, we calculate the relative contributions of different characteristics to the community structure of individual universities and subsequently compare these relative contributions at different universities, measuring for example the importance of common high school affiliation to large state universities and the varying degrees of influence common major can have on the social structure at different universities. The heterogeneity of communities that we observe indicates that these networks typically have multiple organizing factors rather than a single dominant one.

preprint2010arXiv

Competition for Popularity in Bipartite Networks

We present a dynamical model for rewiring and attachment in bipartite networks in which edges are added between nodes that belong to catalogs that can either be fixed in size or growing in size. The model is motivated by an empirical study of data from the video rental service Netflix, which invites its users to give ratings to the videos available in its catalog. We find that the distribution of the number of ratings given by users and that of the number of ratings received by videos both follow a power law with an exponential cutoff. We also examine the activity patterns of Netflix users and find bursts of intense video-rating activity followed by long periods of inactivity. We derive ordinary differential equations to model the acquisition of edges by the nodes over time and obtain the corresponding time-dependent degree distributions. We then compare our results with the Netflix data and find good agreement. We conclude with a discussion of how catalog models can be used to study systems in which agents are forced to choose, rate, or prioritize their interactions from a very large set of options.

preprint2010arXiv

Dynamical Clustering of Exchange Rates

We use techniques from network science to study correlations in the foreign exchange (FX) market over the period 1991--2008. We consider an FX market network in which each node represents an exchange rate and each weighted edge represents a time-dependent correlation between the rates. To provide insights into the clustering of the exchange rate time series, we investigate dynamic communities in the network. We show that there is a relationship between an exchange rate's functional role within the market and its position within its community and use a node-centric community analysis to track the time dynamics of this role. This reveals which exchange rates dominate the market at particular times and also identifies exchange rates that experienced significant changes in market role. We also use the community dynamics to uncover major structural changes that occurred in the FX market. Our techniques are general and will be similarly useful for investigating correlations in other markets.

preprint2010arXiv

Intrinsic Energy Localization through Discrete Gap Breathers in One-Dimensional Diatomic Granular Crystals

We present a systematic study of the existence and stability of discrete breathers that are spatially localized in the bulk of a one-dimensional chain of compressed elastic beads that interact via Hertzian contact. The chain is diatomic, consisting of a periodic arrangement of heavy and light spherical particles. We examine two families of discrete gap breathers: (1) an unstable discrete gap breather that is centered on a heavy particle and characterized by a symmetric spatial energy profile and (2) a potentially stable discrete gap breather that is centered on a light particle and is characterized by an asymmetric spatial energy profile. We investigate their existence, structure, and stability throughout the band gap of the linear spectrum and classify them into four regimes: a regime near the lower optical band edge of the linear spectrum, a moderately discrete regime, a strongly discrete regime that lies deep within the band gap of the linearized version of the system, and a regime near the upper acoustic band edge. We contrast discrete breathers in anharmonic FPU-type diatomic chains with those in diatomic granular crystals, which have a tensionless interaction potential between adjacent particles, and highlight in that the asymmetric nature of the latter interaction potential may lead to a form of hybrid bulk-surface localized solutions.

preprint2010arXiv

Nonlinear Waves in Disordered Diatomic Granular Chains

We investigate the propagation and scattering of highly nonlinear waves in disordered granular chains composed of diatomic (two-mass) units of spheres that interact via Hertzian contact. Using ideas from statistical mechanics, we consider each diatomic unit to be a "spin", so that a granular chain can be viewed as a spin chain composed of units that are each oriented in one of two possible ways. Experiments and numerical simulations both reveal the existence of two different mechanisms of wave propagation: In low-disorder chains, we observe the propagation of a solitary pulse with exponentially decaying amplitude. Beyond a critical level of disorder, the wave amplitude instead decays as a power law, and the wave transmission becomes insensitive to the level of disorder. We characterize the spatio-temporal structure of the wave in both propagation regimes and propose a simple theoretical interpretation for such a transition. Our investigation suggests that an elastic spin chain can be used as a model system to investigate the role of heterogeneities in the propagation of highly nonlinear waves.

preprint2010arXiv

Revisiting Date and Party Hubs: Novel Approaches to Role Assignment in Protein Interaction Networks

The idea of 'date' and 'party' hubs has been influential in the study of protein-protein interaction networks. Date hubs display low co-expression with their partners, whilst party hubs have high co-expression. It was proposed that party hubs are local coordinators whereas date hubs are global connectors. Here we show that the reported importance of date hubs to network connectivity can in fact be attributed to a tiny subset of them. Crucially, these few, extremely central, hubs do not display particularly low expression correlation, undermining the idea of a link between this quantity and hub function. The date/party distinction was originally motivated by an approximately bimodal distribution of hub co-expression; we show that this feature is not always robust to methodological changes. Additionally, topological properties of hubs do not in general correlate with co-expression. Thus, we suggest that a date/party dichotomy is not meaningful and it might be more useful to conceive of roles for protein-protein interactions rather than individual proteins. We find significant correlations between interaction centrality and the functional similarity of the interacting proteins.

preprint2010arXiv

The Function of Communities in Protein Interaction Networks at Multiple Scales

Background: If biology is modular then clusters, or communities, of proteins derived using only protein interaction network structure should define protein modules with similar biological roles. We investigate the link between biological modules and network communities in yeast and its relationship to the scale at which we probe the network. Results: Our results demonstrate that the functional homogeneity of communities depends on the scale selected, and that almost all proteins lie in a functionally homogeneous community at some scale. We judge functional homogeneity using a novel test and three independent characterizations of protein function, and find a high degree of overlap between these measures. We show that a high mean clustering coefficient of a community can be used to identify those that are functionally homogeneous. By tracing the community membership of a protein through multiple scales we demonstrate how our approach could be useful to biologists focusing on a particular protein. Conclusions: We show that there is no one scale of interest in the community structure of the yeast protein interaction network, but we can identify the range of resolution parameters that yield the most functionally coherent communities, and predict which communities are most likely to be functionally homogeneous.

preprint2010arXiv

The unreasonable effectiveness of tree-based theory for networks with clustering

We demonstrate that a tree-based theory for various dynamical processes yields extremely accurate results for several networks with high levels of clustering. We find that such a theory works well as long as the mean intervertex distance $\ell$ is sufficiently small - i.e., as long as it is close to the value of $\ell$ in a random network with negligible clustering and the same degree-degree correlations. We confirm this hypothesis numerically using real-world networks from various domains and on several classes of synthetic clustered networks. We present analytical calculations that further support our claim that tree-based theories can be accurate for clustered networks provided that the networks are "sufficiently small" worlds.

preprint2009arXiv

Communities in Networks

We survey some of the concepts, methods, and applications of community detection, which has become an increasingly important area of network science. To help ease newcomers into the field, we provide a guide to available methodology and open problems, and discuss why scientists from diverse backgrounds are interested in these problems. As a running theme, we emphasize the connections of community detection to problems in statistical physics and computational optimization.

preprint2009arXiv

Dynamic communities in multichannel data: An application to the foreign exchange market during the 2007--2008 credit crisis

We study the cluster dynamics of multichannel (multivariate) time series by representing their correlations as time-dependent networks and investigating the evolution of network communities. We employ a node-centric approach that allows us to track the effects of the community evolution on the functional roles of individual nodes without having to track entire communities. As an example, we consider a foreign exchange market network in which each node represents an exchange rate and each edge represents a time-dependent correlation between the rates. We study the period 2005-2008, which includes the recent credit and liquidity crisis. Using dynamical community detection, we find that exchange rates that are strongly attached to their community are persistently grouped with the same set of rates, whereas exchange rates that are important for the transfer of information tend to be positioned on the edges of communities. Our analysis successfully uncovers major trading changes that occurred in the market during the credit crisis.

preprint2009arXiv

Experimental Results Related to Discrete Nonlinear Schrödinger Equations

In this chapter, we discuss experiments that realize the discrete nonlinear Schrödinger (DNLS) equations. The relevance of such descriptions arises from the competition of three common features: nonlinearity, dispersion, and a medium to large level of (periodic, quasiperiodic, or random) discreteness in space. DNLS equations have been especially prevalent in atomic and molecular physics in the study of Bose-Einstein condensates in optical lattices or superlattices; and in nonlinear optics in the description of pulse propagation in waveguide arrays and photorefractive crystals. New experiments in both nonlinear optics and Bose-Einstein condensation provide new challenges for DNLS models, and DNLS and related equations have also recently been used to make important predictions in novel physical settings such as the study of composite metamaterials and arrays of superconducting devices.

preprint2009arXiv

Localized Breathing Modes in Granular Crystals with Defects

We investigate nonlinear localized modes at light-mass impurities in a one-dimensional, strongly-compressed chain of beads under Hertzian contacts. Focusing on the case of one or two such "defects", we analyze the problem's linear limit to identify the system eigenfrequencies and the linear defect modes. We then examine the bifurcation of nonlinear defect modes from their linear counterparts and study their linear stability in detail. We identify intriguing differences between the case of impurities in contact and ones that are not in contact. We find that the former bears similarities to the single defect case, whereas the latter features symmetry-breaking bifurcations with interesting static and dynamic implications.

preprint2009arXiv

Mutually-Antagonistic Interactions in Baseball Networks

We formulate the head-to-head matchups between Major League Baseball pitchers and batters from 1954 to 2008 as a bipartite network of mutually-antagonistic interactions. We consider both the full network and single-season networks, which exhibit interesting structural changes over time. We find interesting structure in the network and examine their sensitivity to baseball's rule changes. We then study a biased random walk on the matchup networks as a simple and transparent way to compare the performance of players who competed under different conditions and to include information about which particular players a given player has faced. We find that a player's position in the network does not correlate with his success in the random walker ranking but instead has a substantial effect on its sensitivity to changes in his own aggregate performance.

preprint2008arXiv

Dissipative Solitary Waves in Granular Crystals

We provide a quantitative characterization of dissipative effects in one-dimensional granular crystals. We use the propagation of highly nonlinear solitary waves as a diagnostic tool and develop optimization schemes that allow one to compute the relevant exponents and prefactors of the dissipative terms in the equations of motion. We thereby propose a quantitatively-accurate extension of the Hertzian model that encompasses dissipative effects via a discrete Laplacian of the velocities. Experiments and computations with steel, brass, and polytetrafluoroethylene reveal a {\em common} dissipation exponent with a material-dependent prefactor.

preprint2006arXiv

Random Walker Ranking for NCAA Division I-A Football

Each December, college football fans and pundits across America debate which two teams should meet in the NCAA Division I-A National Championship game. The Bowl Championship Series (BCS) standings employed to select the teams invited to this game are intended to provide an unequivocal #1 v. #2 game for the championship; however, this selection process has itself been highly controversial in recent years. The computer algorithms that constitute one part of the BCS standings often act as lightning rods for the controversy, in part because they are inadequately explained to the public. We present an alternative algorithm that is simply explained yet remains effective at ranking the best teams. We define a ranking in terms of biased random walkers on the graph formed by the schedule of games played, with two teams (vertices) connected by an edge if they played each other. Each random walker moves from team to team by selecting a game and "voting" for its winner with probability p, tracing out a never-ending path motivated by the "my team beat your team" argument. We study the statistical properties of a collection of such walkers, relate the rankings to the community structure of the underlying network, and demonstrate the results for recent NCAA Division I-A seasons. We also discuss the algorithm's asymptotic behavior, illustrated with some analytically tractable cases for round-robin tournaments, and discuss possible generalizations.

preprint2005arXiv

Dynamics and Manipulation of Matter-Wave Solitons in Optical Superlattices

We analyze the existence and stability of bright, dark, and gap matter-wave solitons in optical superlattices. Then, using these properties, we show that (time-dependent) ``dynamical superlattices'' can be used to controllably place, guide, and manipulate these solitons. In particular, we use numerical experiments to displace solitons by turning on a secondary lattice structure, transfer solitons from one location to another by shifting one superlattice substructure relative to the other, and implement solitonic ``path-following'', in which a matter wave follows the time-dependent lattice substructure into oscillatory motion.

preprint2004arXiv

A Graphical User Interface to Simulate Classical Billiard Systems

Classical billiards constitute an important class of dynamical systems. They have not only been in used in mathematical disciplines such as ergodic theory, but their properties demonstrate fundamental physical phenomena that can be observed in laboratory settings. This document provides instructions for a Matlab module that simulates classical billiard systems. It is intended to be used as both a research and teaching tool. At present, the program efficiently simulates tables that are constructed entirely from line segments and elliptical arcs. It functions less reliably for tables with more complex boundary components. The program and documentation can be downloaded from \textit{http://www.math.gatech.edu/$\sim$mason/papers/}.