Finding network communities using modularity density

Many real-world complex networks exhibit a community structure, in which the modules correspond to actual functional units. Identifying these communities is a key challenge for scientists. A common approach is to search for the network partition that maximizes a quality function. Here, we present a detailed analysis of a recently proposed function, namely modularity density. We show that it does not incur in the drawbacks suffered by traditional modularity, and that it can identify networks without ground-truth community structure, deriving its analytical dependence on link density in generic random graphs. In addition, we show that modularity density allows an easy comparison between networks of different sizes, and we also present some limitations that methods based on modularity density may suffer from. Finally, we introduce an efficient, quadratic community detection algorithm based on modularity density maximization, validating its accuracy against theoretical predictions and on a set of benchmark networks

Data Mining a Medieval Medical Text Reveals Patterns in Ingredient Choice That Reflect Biological Activity against Infectious Agents

We used established methodologies from network science to identify patterns in medicinal ingredient combinations in a key medieval text, the 15th-century Lylye of Medicynes, focusing on recipes for topical treatments for symptoms of microbial infection. We conducted experiments screening the antimicrobial activity of selected ingredients. These experiments revealed interesting examples of ingredients that potentiated or interfered with each other’s activity and that would be useful bases for future, more detailed experiments. Our results highlight (i) the potential to use methodologies from network science to analyze medieval data sets and detect patterns of ingredient combination, (ii) the potential of interdisciplinary collaboration to reveal different aspects of the ethnopharmacology of historical medical texts, and (iii) the potential development of novel therapeutics inspired by premodern remedies in a time of increased need for new antibiotics.The pharmacopeia used by physicians and laypeople in medieval Europe has largely been dismissed as placebo or superstition. While we now recognize that some of the materia medica used by medieval physicians could have had useful biological properties, research in this area is limited by the labor-intensive process of searching and interpreting historical medical texts. Here, we demonstrate the potential power of turning medieval medical texts into contextualized electronic databases amenable to exploration by the use of an algorithm. We used established methodologies from network science to reveal patterns in ingredient selection and usage in a key text, the 15th-century Lylye of Medicynes, focusing on remedies to treat symptoms of microbial infection. In providing a worked example of data-driven textual analysis, we demonstrate the potential of this approach to encourage interdisciplinary collaboration and to shine a new light on the ethnopharmacology of historical medical texts

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.Comment: 5 pages, 3 figure

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.Comment: 5 pages, 4 figure

Synchronization in networks with multiple interaction layers

The structure of many real-world systems is best captured by networks consisting of several interaction layers. Understanding how a multilayered structure of connections affects the synchronization properties of dynamical systems evolving on top of it is a highly relevant endeavor in mathematics and physics and has potential applications in several socially relevant topics, such as power grid engineering and neural dynamics. We propose a general framework to assess the stability of the synchronized state in networks with multiple interaction layers, deriving a necessary condition that generalizes the master stability function approach. We validate our method by applying it to a network of Rössler oscillators with a double layer of interactions and show that highly rich phenomenology emerges from this. This includes cases where the stability of synchronization can be induced even if both layers would have individually induced unstable synchrony, an effect genuinely arising from the true multilayer structure of the interactions among the units in the network

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.Comment: 4 pages, 2 figure

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 and 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 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 ${\mathcal{O}}({NM})$ where N is the number of nodes and M is the number of edges

Synchronization in dynamical networks with unconstrained structure switching

We provide a rigorous solution to the problem of constructing a structural evolution for a network of coupled identical dynamical units that switches between specified topologies without constraints on their structure. The evolution of the structure is determined indirectly, from a carefully built transformation of the eigenvector matrices of the coupling Laplacians, which are guaranteed to change smoothly in time. In turn, this allows to extend the Master Stability Function formalism, which can be used to assess the stability of a synchronized state. This approach is independent from the particular topologies that the network visits, and is not restricted to commuting structures. Also, it does not depend on the time scale of the evolution, which can be faster than, comparable to, or even secular with respect to the the dynamics of the units.Comment: 8 pages, 3 figures. arXiv admin note: text overlap with arXiv:1407.074

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.Comment: 5 pages, 4 figure

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.Comment: 12 pages, 8 figure