Catalog footprint

What is connected

26works
21topics
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

26 published item(s)

preprint2025arXiv

Complexity and dynamics of partially symmetric random neural networks

Neural circuits exhibit structured connectivity, including an overrepresentation of reciprocal connections between neuron pairs. Despite important advances, a full understanding of how such partial symmetry in connectivity shapes neural dynamics remains elusive. Here we ask how correlations between reciprocal connections in a random, recurrent neural network affect phase-space complexity, defined as the exponential proliferation rate (with network size) of the number of fixed points that accompanies the transition to chaotic dynamics. We find a striking pattern: partial anti-symmetry strongly amplifies complexity, while partial symmetry suppresses it. These opposing trends closely track changes in other measures of dynamical behavior, such as dimensionality, Lyapunov exponents, and transient path length, supporting the view that fixed-point structure is a key determinant of network dynamics. Thus, positive reciprocal correlations favor low-dimensional, slowly varying activity, whereas negative correlations promote high-dimensional, rapidly fluctuating chaotic activity. These results yield testable predictions about the link between connection reciprocity, neural dynamics and function.

preprint2020arXiv

Resolution limit revisited: community detection using generalized modularity density

Various attempts have been made in recent years to solve the Resolution Limit (RL) problem in community detection by considering variants of the modularity metric in the detection algorithms. These metrics purportedly largely mitigate the RL problem and are preferable to modularity in many realistic scenarios. However, they are not generally suitable for analyzing weighted networks or for detecting hierarchical community structure. Resolution limit problems can be complicated, though, and in particular it can be unclear when it should be considered as problem. In this paper, we introduce a metric that we call generalized modularity density $Q_g$ that eliminates the RL problem at any desired resolution and is easily extendable to study weighted and hierarchical networks. We also propose a benchmark test to quantify the resolution limit problem, examine various modularity-like metrics to show that the new metric $Q_g$ performs best, and show that $Q_g$ can identify modular structure in real-world and artificial networks that is otherwise hidden.

preprint2016arXiv

Effects of single- and multi-substituted Zn ions in doped-122 type iron-based superconductors

Recent experiments on Zn-substituted 122-type iron-based superconductors (FeSCs) at electron- and hole- doped region provide us with a testing ground for understanding the effect of Zn impurities in these systems. Our first-principle calculations of the electronic structure reveal that the Zn 3d orbitals are far below the Fermi level and chemically inactive, while the Zn 4s-orbital is partially occupied and its wave function overlapping with those 3d-orbitals of neighboring Fe-ions. This suggests that the impurity effect is originating in the Zn 4s-orbital, not its 3d-orbitals. Employing a phenomenological two-orbital lattice model for 122-FeSCs and the self-consistent Bogoliubov-de Gennes equations, we study how the Zn-impurities suppress the superconductivity in electron- and hole- doped compounds. Our obtained results qualitatively agree with the experimental measurements.

preprint2016arXiv

The Many-agent limit of the Extreme Introvert-Extrovert model

We consider a toy model of interacting extrovert and introvert agents introduced earlier by Liu et al [Europhys. Lett. {\bf 100} (2012) 66007]. The number of extroverts, and introverts is $N$ each. At each time step, we select an agent at random, and allow her to modify her state. If an extrovert is selected, she adds a link at random to an unconnected introvert. If an introvert is selected, she removes one of her links. The set of $N^2$ links evolves in time, and may be considered as a set of Ising spins on an $N \times N$ square-grid with single-spin-flip dynamics. This dynamics satisfies detailed balance condition, and the probability of different spin configurations in the steady state can be determined exactly. The effective hamiltonian has long-range multi-spin couplings that depend on the row and column sums of spins. If the relative bias of choosing an extrovert over an introvert is varied, this system undergoes a phase transition from a state with very few links to one in which most links are occupied. We show that the behavior of the system can be determined exactly in the limit of large $N$. The behavior of large fluctuations in the total numer of links near the phase transition is determined. We also discuss two variations, called egalitarian and elitist agents, when the agents preferentially add or delete links to their least/ most-connected neighbor. These shows interesting cooperative behavior.

preprint2015arXiv

A unified description of superconducting pairing symmetry in electron-doped Fe-based-122 compounds

The pairing symmetry is examined in highly electron-doped Ba(Fe$_{1-x}$Co$_x$As)$_2$ and A$_y$Fe$_2$Se$_2$ (with A=K, Cs) compounds, with similar crystallographic and electronic band structures. Starting from a phenomenological two-orbital model, we consider nearest-neighbor and next-nearest-neighbor intraorbital pairing interactions on the Fe square lattice. In this model, we find a unified description of the evolution from $s_\pm$-wave pairing ($2.0 < n \lesssim 2.4$) to $d$-wave pairing ($2.4 \lesssim n \lesssim 2.5$) as a function of electron filling. In the crossover region a novel time-reversal symmetry breaking state with $s_\pm+id$ pairing symmetry emerges. This minimal model offers an overall picture of the evolution of superconductivity with electron doping for both $s_\pm$-wave [Ba(Fe$_{1-x}$Co$_x$As)$_2$] and $d$-wave [A$_y$Fe$_2$Se$_2$] pairing, as long as the dopants only play the role of a charge reservoir. However, the situation is more complicated for Ba(Fe$_{1-x}$Co$_x$As)$_2$. A real-space study further shows that when the impurity scattering effects of Co dopants are taken into account, the superconductivity is completely suppressed for $n > 2.4$. This preempts any observation of $d$-wave pairing in this compound, in contrast to A$_y$Fe$_2$Se$_2$.

preprint2015arXiv

Exact sampling of graphs with prescribed degree correlations

Many real-world networks exhibit correlations between the node degrees. For instance, in social networks nodes tend to connect to nodes of similar degree. Conversely, in biological and technological networks, high-degree nodes tend to be linked with low-degree nodes. Degree correlations also affect the dynamics of processes supported by a network structure, such as the spread of opinions or epidemics. The proper modelling of these systems, i.e., without uncontrolled biases, requires the sampling of networks with a specified set of constraints. We present a solution to the sampling problem when the constraints imposed are the degree correlations. In particular, we develop an efficient and exact method to construct and sample graphs with a specified joint-degree matrix, which is a matrix providing the number of edges between all the sets of nodes of a given degree, for all degrees, thus completely specifying all pairwise degree correlations, and additionally, the degree sequence itself. Our algorithm always produces independent samples without backtracking. The complexity of the graph construction algorithm is O(NM) where N is the number of nodes and M is the number of edges.

preprint2015arXiv

Networks with preferred degree: A mini-review and some new results

Since their inception about a decade ago, dynamic networks which adapt to the state of the nodes have attracted much attention. One simple case of such an adaptive dynamics is a model of social networks in which individuals are typically comfortable with a certain number of contacts, i.e., preferred degrees. This paper is partly a review of earlier work of single homogeneous systems and ones with two interacting networks, and partly a presentation of some new results. In general, the dynamics does not obey detailed balance and the stationary distributions are not known analytically. A particular limit of the latter is a system of extreme introverts and extroverts - the XIE model. Remarkably, in this case, the detailed balance condition is satisfied, the exact distribution and an effective Hamiltonian can be found explicitly. Further, the model exhibits a phase transition in which the total number of links in the system - a macroscopically interesting quantity, displays an extreme Thouless effect. We show that in the limit of large populations and away from the transition, the model reduces to one with non-interacting agents of the majority subgroup. We determine the nature of fluctuations near the transition. We also introduce variants of the model where the agents show preferential attachment or detachment. There are significant changes to the degree distributions in the steady state, some of which can be understood by theoretical arguments and some remain to be explored. Many intriguing questions are posed, providing some food for thought and avenues for future research.

preprint2015arXiv

Quantifying randomness in real networks

Represented as graphs, real networks are intricate combinations of order and disorder. Fixing some of the structural properties of network models to their values observed in real networks, many other properties appear as statistical consequences of these fixed observables, plus randomness in other respects. Here we employ the $dk$-series, a complete set of basic characteristics of the network structure, to study the statistical dependencies between different network properties. We consider six real networks---the Internet, US airport network, human protein interactions, technosocial web of trust, English word network, and an fMRI map of the human brain---and find that many important local and global structural properties of these networks are closely reproduced by $dk$-random graphs whose degree distributions, degree correlations, and clustering are as in the corresponding real network. We discuss important conceptual, methodological, and practical implications of this evaluation of network randomness, and release software to generate $dk$-random graphs.

preprint2014arXiv

Fast and accurate determination of modularity and its effect size

We present a fast spectral algorithm for community detection in complex networks. Our method searches for the partition with the maximum value of the modularity via the interplay of several refinement steps that include both agglomeration and division. We validate the accuracy of the algorithm by applying it to several real-world benchmark networks. On all these, our algorithm performs as well or better than any other known polynomial scheme. This allows us to extensively study the modularity distribution in ensembles of Erdős-Rényi networks, producing theoretical predictions for means and variances inclusive of finite-size corrections. Our work provides a way to accurately estimate the effect size of modularity, providing a $z$-score measure of it and enabling a more informative comparison of networks with different numbers of nodes and links.

preprint2014arXiv

Mesoscopic structures and the Laplacian spectra of random geometric graphs

We investigate the Laplacian spectra of random geometric graphs (RGGs). The spectra are found to consist of both a discrete and a continuous part. The discrete part is a collection of Dirac delta peaks at integer values roughly centered around the mean degree. The peaks are mainly due to the existence of mesoscopic structures that occur far more abundantly in RGGs than in non-spatial networks. The probability of certain mesoscopic structures is analytically calculated for one-dimensional RGGs and they are shown to produce integer-valued eigenvalues that comprise a significant fraction of the spectrum, even in the large network limit. A phenomenon reminiscent of Bose-Einstein condensation in the appearance of zero eigenvalues is also found.

preprint2014arXiv

Symmetry in Critical Random Boolean Network Dynamics

Using Boolean networks as prototypical examples, the role of symmetry in the dynamics of heterogeneous complex systems is explored. We show that symmetry of the dynamics, especially in critical states, is a controlling feature that can be used both to greatly simplify analysis and to characterize different types of dynamics. Symmetry in Boolean networks is found by determining the frequency at which the various Boolean output functions occur. There are classes of functions that consist of Boolean functions that behave similarly. These classes are orbits of the controlling symmetry group. We find that the symmetry that controls the critical random Boolean networks is expressed through the frequency by which output functions are utilized by nodes that remain active on dynamical attractors. This symmetry preserves canalization, a form of network robustness. We compare it to a different symmetry known to control the dynamics of an evolutionary process that allows Boolean networks to organize into a critical state. Our results demonstrate the usefulness and power of using the symmetry of the behavior of the nodes to characterize complex network dynamics, and introduce a novel approach to the analysis of heterogeneous complex systems.

preprint2012arXiv

Magnetic and superconducting structures near twin boundaries in low doped Fe-pnictides

The effects of twin boundaries (TBs) on the complex interaction between magnetism and superconductivity in slightly electron-doped Ba(Ca)(FeAs)$_2$ superconductors are investigated. The spatial distributions of the magnetic, superconducting and charge density orders near two different types of TBs are calculated. We find that TBs corresponding to a 90$^\circ$ lattice rotation in the a-b plane enable magnetic domain walls to form with only a small effective Coulomb interaction between valance electrons, and that superconductivity is enhanced at such TBs. Contrastingly, we find that superconductivity is suppressed at TBs corresponding to an asymmetrical placement of As atoms with respect to the Fe atoms in the a-b plane.

preprint2012arXiv

Robust Detection of Hierarchical Communities from Escherichia coli Gene Expression Data

Determining the functional structure of biological networks is a central goal of systems biology. One approach is to analyze gene expression data to infer a network of gene interactions on the basis of their correlated responses to environmental and genetic perturbations. The inferred network can then be analyzed to identify functional communities. However, commonly used algorithms can yield unreliable results due to experimental noise, algorithmic stochasticity, and the influence of arbitrarily chosen parameter values. Furthermore, the results obtained typically provide only a simplistic view of the network partitioned into disjoint communities and provide no information of the relationship between communities. Here, we present methods to robustly detect coregulated and functionally enriched gene communities and demonstrate their application and validity for Escherichia coli gene expression data. Applying a recently developed community detection algorithm to the network of interactions identified with the context likelihood of relatedness (CLR) method, we show that a hierarchy of network communities can be identified. These communities significantly enrich for gene ontology (GO) terms, consistent with them representing biologically meaningful groups. Further, analysis of the most significantly enriched communities identified several candidate new regulatory interactions. The robustness of our methods is demonstrated by showing that a core set of functional communities is reliably found when artificial noise, modeling experimental noise, is added to the data. We find that noise mainly acts conservatively, increasing the relatedness required for a network link to be reliably assigned and decreasing the size of the core communities, rather than causing association of genes into new communities.

preprint2011arXiv

All scale-free networks are sparse

We study the realizability of scale free-networks with a given degree sequence, showing that the fraction of realizable sequences undergoes two first-order transitions at the values 0 and 2 of the power-law exponent. We substantiate this finding by analytical reasoning and by a numerical method, proposed here, based on extreme value arguments, which can be applied to any given degree distribution. Our results reveal a fundamental reason why large scale-free networks without constraints on minimum and maximum degree must be sparse.

preprint2011arXiv

Canalization in the Critical States of Highly Connected Networks of Competing Boolean Nodes

Canalization is a classic concept in Developmental Biology that is thought to be an important feature of evolving systems. In a Boolean network it is a form of network robustness in which a subset of the input signals control the behavior of a node regardless of the remaining input. It has been shown that Boolean networks can become canalized if they evolve through a frustrated competition between nodes. This was demonstrated for large networks in which each node had K=3 inputs. Those networks evolve to a critical steady-state at the boarder of two phases of dynamical behavior. Moreover, the evolution of these networks was shown to be associated with the symmetry of the evolutionary dynamics. We extend these results to the more highly connected K>3 cases and show that similar canalized critical steady states emerge with the same associated dynamical symmetry, but only if the evolutionary dynamics is biased toward homogeneous Boolean functions.

preprint2011arXiv

Collectively optimal routing for congested traffic limited by link capacity

We show that the capacity of a complex network that models a city street grid to support congested traffic can be optimized by using routes that collectively minimize the maximum ratio of betweenness to capacity in any link. Networks with a heterogeneous distribution of link capacities and with a heterogeneous transport load are considered. We find that overall traffic congestion and average travel times can be significantly reduced by a judicious use of slower, smaller capacity links.

preprint2010arXiv

Anomalous ordering in inhomogeneously strained materials

We study a continuous quasi-two-dimensional order-disorder phase transition that occurs in a simple model of a material that is inhomogeneously strained due to the presence of dislocation lines. Performing Monte Carlo simulations of different system sizes and using finite size scaling, we measure critical exponents describing the transition of beta=0.18\pm0.02, gamma=1.0\pm0.1, and alpha=0.10\pm0.02. Comparable exponents have been reported in a variety of physical systems. These systems undergo a range of different types of phase transitions, including structural transitions, exciton percolation, and magnetic ordering. In particular, similar exponents have been found to describe the development of magnetic order at the onset of the pseudogap transition in high-temperature superconductors. Their common universal critical exponents suggest that the essential physics of the transition in all of these physical systems is the same as in our model. We argue that the nature of the transition in our model is related to surface transitions, although our model has no free surface.

preprint2010arXiv

Depth-dependent ordering, two-length-scale phenomena and crossover behavior in a crystal featuring a skin-layer with defects

Structural defects in a crystal are responsible for the "two length-scale" behavior, in which a sharp central peak is superimposed over a broad peak in critical diffuse X-ray scattering. We have previously measured the scaling behavior of the central peak by scattering from a near-surface region of a V2H crystal, which has a first-order transition in the bulk. As the temperature is lowered toward the critical temperature, a crossover in critical behavior is seen, with the temperature range nearest to the critical point being characterized by mean field exponents. Near the transition, a small two-phase coexistence region is observed. The values of transition and crossover temperatures decay with depth. An explanation of these experimental results is here proposed by means of a theory in which edge dislocations in the near-surface region occur in walls oriented in the two directions normal to the surface. The strain caused by the dislocation lines causes the ordering in the crystal to occur as growth of roughly cylindrically shaped regions. After the regions have reached a certain size, the crossover in the critical behavior occurs, and mean field behavior prevails. At a still lower temperature, the rest of the material between the cylindrical regions orders via a weak first-order transition.

preprint2010arXiv

Efficient and exact sampling of simple graphs with given arbitrary degree sequence

Uniform sampling from graphical realizations of a given degree sequence is a fundamental component in simulation-based measurements of network observables, with applications ranging from epidemics, through social networks to Internet modeling. Existing graph sampling methods are either link-swap based (Markov-Chain Monte Carlo algorithms) or stub-matching based (the Configuration Model). Both types are ill-controlled, with typically unknown mixing times for link-swap methods and uncontrolled rejections for the Configuration Model. Here we propose an efficient, polynomial time algorithm that generates statistically independent graph samples with a given, arbitrary, degree sequence. The algorithm provides a weight associated with each sample, allowing the observable to be measured either uniformly over the graph ensemble, or, alternatively, with a desired distribution. Unlike other algorithms, this method always produces a sample, without back-tracking or rejections. Using a central limit theorem-based reasoning, we argue, that for large N, and for degree sequences admitting many realizations, the sample weights are expected to have a lognormal distribution. As examples, we apply our algorithm to generate networks with degree sequences drawn from power-law distributions and from binomial distributions.

preprint2010arXiv

Phase Diagram for a 2-D Two-Temperature Diffusive XY Model

Using Monte Carlo simulations, we determine the phase diagram of a diffusive two-temperature XY model. When the two temperatures are equal the system becomes the equilibrium XY model with the continuous Kosterlitz-Thouless (KT) vortex-antivortex unbinding phase transition. When the two temperatures are unequal the system is driven by an energy flow through the system from the higher temperature heat-bath to the lower temperature one and reaches a far-from-equilibrium steady state. We show that the nonequilibrium phase diagram contains three phases: A homogenous disordered phase and two phases with long range, spin-wave order. Two critical lines, representing continuous phase transitions from a homogenous disordered phase to two phases of long range order, meet at the equilibrium the KT point. The shape of the nonequilibrium critical lines as they approach the KT point is described by a crossover exponent of phi = 2.52 \pm 0.05. Finally, we suggest that the transition between the two phases with long-range order is first-order, making the KT-point where all three phases meet a bicritical point.

preprint2009arXiv

Depth-dependent critical behavior in V2H

Using X-ray diffuse scattering, we investigate the critical behavior of an order-disorder phase transition in a defective "skin-layer" of V2H. In the skin-layer, there exist walls of dislocation lines oriented normal to the surface. The density of dislocation lines within a wall decreases continuously with depth. We find that, because of this inhomogeneous distribution of defects, the transition effectively occurs at a depth-dependent local critical temperature. A depth-dependent scaling law is proposed to describe the corresponding critical ordering behavior.

preprint2009arXiv

Improved community structure detection using a modified fine tuning strategy

The community structure of a complex network can be determined by finding the partitioning of its nodes that maximizes modularity. Many of the proposed algorithms for doing this work by recursively bisecting the network. We show that this unduely constrains their results, leading to a bias in the size of the communities they find and limiting their effectivness. To solve this problem, we propose adding a step to the existing algorithms that does not increase the order of their computational complexity. We show that, if this step is combined with a commonly used method, the identified constraint and resulting bias are removed, and its ability to find the optimal partitioning is improved. The effectiveness of this combined algorithm is also demonstrated by using it on real-world example networks. For a number of these examples, it achieves the best results of any known algorithm.

preprint2007arXiv

Finite size effects and symmetry breaking in the evolution of networks of competing Boolean nodes

The effects of the finite size of the network on the evolutionary dynamics of a Boolean network are analyzed. In the model considered, Boolean networks evolve via a competition between nodes that punishes those in the majority. It is found that finite size networks evolve in a fundamentally different way than infinitely large networks do. The symmetry of the evolutionary dynamics of infinitely large networks that selects for canalizing Boolean functions is broken in the evolutionary dynamics of finite size networks. In finite size networks there is an additional selection for input inverting Boolean functions that output a value opposite to the majority of input values. These results are revealed through an empirical study of the model that calculates the frequency of occurrence of the different possible Boolean functions. Classes of functions are found to occur with the same frequency. Those classes depend on the symmetry of the evolutionary dynamics and correspond to orbits of the relevant symmetry group. The empirical results match analytic results, determined by utilizing Polya's theorem, for the number of orbits expected in both finite size and infinitely large networks. The reason for the symmetry breaking in the evolutionary dynamics is found to be due to the need for nodes in finite size networks to behave differently in order to cooperate so that the system collectively performs as well as possible. The results suggest that both finite size effects and symmetry are important for understanding the evolution of real-world complex networks, including genetic regulatory networks.

preprint2007arXiv

Hurst Exponents, Markov Processes, and Fractional Brownian motion

There is much confusion in the literature over Hurst exponents. Recently, we took a step in the direction of eliminating some of the confusion. One purpose of this paper is to illustrate the difference between fBm on the one hand and Gaussian Markov processes where H not equal to 1/2 on the other. The difference lies in the increments, which are stationary and correlated in one case and nonstationary and uncorrelated in the other. The two- and one-point densities of fBm are constructed explicitly. The two-point density doesn't scale. The one-point density is identical with that for a Markov process with H not 1/2. We conclude that both Hurst exponents and histograms for one point densities are inadequate for deducing an underlying stochastic dynamical system from empirical data.

preprint1999arXiv

Braided Rivers and Superconducting Vortex Avalanches

Magnetic vortices intermittently flow through preferred channels when they are forced in or out of a superconductor. We study this behavior using a cellular model, and find that the vortex flow can make braided rivers strikingly similar to aerial photographs of braided fluvial rivers, such as the Brahmaputra. By developing an analysis technique suitable for characterizing a self-affine (multi)fractal, the scaling properties of the braided vortex rivers in the model are compared with those of braided fluvial rivers. We suggest that avalanche dynamics leads to braiding in both cases.

preprint1998arXiv

A Simple Model of Superconducting Vortex Avalanches

We introduce a simple lattice model of superconducting vortices driven by repulsive interactions through a random pinning potential. The model describes the behavior at the scale of the London length lambda or larger. It self-organizes to a critical state, characterized by a constant flux density gradient, where the activity takes place in terms of avalanches spanning all length scales up to the system size. We determine scaling relations as well as four universal critical exponents for avalanche moments and durations: tau = 1.63 +/- 0.02, D = 2.7 +/- 0.1, z = 1.5 +/- 0.1, and tau_t = 2.13 +/- 0.14, for the system driven at the boundary.