Ising-like model replicating time-averaged spiking behaviour of in vitro neuronal networks.

We analyze time-averaged experimental data from in vitro activities of neuronal networks. Through a Pairwise Maximum-Entropy method, we identify through an inverse binary Ising-like model the local fields and interaction couplings which best reproduce the average activities of each neuron as well as the statistical correlations between the activities of each pair of neurons in the system. The specific information about the type of neurons is mainly stored in the local fields, while a symmetric distribution of interaction constants seems generic. Our findings demonstrate that, despite not being directly incorporated into the inference approach, the experimentally observed correlations among groups of three neurons are accurately captured by the derived Ising-like model. Within the context of the thermodynamic analogy inherent to the Ising-like models developed in this study, our findings additionally indicate that these models demonstrate characteristics of second-order phase transitions between ferromagnetic and paramagnetic states at temperatures above, but close to, unity. Considering that the operating temperature utilized in the Maximum-Entropy method is T o = 1 , this observation further expands the thermodynamic conceptual parallelism postulated in this work for the manifestation of criticality in neuronal network behavior.

Morphological Transition between Patterns Formed by Threads of Magnetic Beads.

Magnetic beads attract each other, forming chains. We push such chains into an inclined Hele-Shaw cell and discover that they spontaneously form self-similar patterns. Depending on the angle of inclination of the cell, two completely different situations emerge; namely, above the static friction angle the patterns resemble the stacking of a rope and below they look similar to a fortress from above. Moreover, locally the first pattern forms a square lattice, while the second pattern exhibits triangular symmetry. For both patterns, the size distributions of enclosed areas follow power laws. We characterize the morphological transition between the two patterns experimentally and numerically and explain the change in polarization as a competition between friction-induced buckling and gravity.

Phase transitions in a conservative game of life.

We investigate the dynamics of a conservative version of Conway's Game of Life, in which a pair consisting of a dead and a living cell can switch their states following Conway's rules but only by swapping their positions, irrespective of their mutual distance. Our study is based on square-lattice simulations as well as a mean-field calculation. As the density of dead cells is increased, we identify a discontinuous phase transition between an inactive phase, in which the dynamics freezes after a finite time, and an active phase, in which the dynamics persists indefinitely in the thermodynamic limit. Further increasing the density of dead cells leads the system back to an inactive phase via a second transition, which is continuous on the square lattice but discontinuous in the mean-field limit.

A universal approach for drainage basins.

Drainage basins are essential to Geohydrology and Biodiversity. Defining those regions in a simple, robust and efficient way is a constant challenge in Earth Science. Here, we introduce a model to delineate multiple drainage basins through an extension of the Invasion Percolation-Based Algorithm (IPBA). In order to prove the potential of our approach, we apply it to real and artificial datasets. We observe that the perimeter and area distributions of basins and anti-basins display long tails extending over several orders of magnitude and following approximately power-law behaviors. Moreover, the exponents of these power laws depend on spatial correlations and are invariant under the landscape orientation, not only for terrestrial, but lunar and martian landscapes. The terrestrial and martian results are statistically identical, which suggests that a hypothetical martian river would present similarity to the terrestrial rivers. Finally, we propose a theoretical value for the Hack's exponent based on the fractal dimension of watersheds, γ = D/2. We measure γ = 0.54 ± 0.01 for Earth, which is close to our estimation of γ ≈ 0.55. Our study suggests that Hack's law can have its origin purely in the maximum and minimum lines of the landscapes.

Elastic Backbone Defines a New Transition in the Percolation Model.

The elastic backbone is the set of all shortest paths. We found a new phase transition at p_{eb} above the classical percolation threshold at which the elastic backbone becomes dense. At this transition in 2D, its fractal dimension is 1.750±0.003, and one obtains a novel set of critical exponents ß_{eb}=0.50±0.02, γ_{eb}=1.97±0.05, and ν_{eb}=2.00±0.02, fulfilling consistent critical scaling laws. Interestingly, however, the hyperscaling relation is violated. Using Binder's cumulant, we determine, with high precision, the critical probabilities p_{eb} for the triangular and tilted square lattice for site and bond percolation. This transition describes a sudden rigidification as a function of density when stretching a damaged tissue.

Gender differences in scientific collaborations: Women are more egalitarian than men.

By analyzing a unique dataset of more than 270,000 scientists, we discovered substantial gender differences in scientific collaborations. While men are more likely to collaborate with other men, women are more egalitarian. This is consistently observed over all fields and regardless of the number of collaborators a scientist has. The only exception is observed in the field of engineering, where this gender bias disappears with increasing number of collaborators. We also found that the distribution of the number of collaborators follows a truncated power law with a cut-off that is gender dependent and related to the gender differences in the number of published papers. Considering interdisciplinary research, our analysis shows that men and women behave similarly across fields, except in the case of natural sciences, where women with many collaborators are more likely to have collaborators from other fields.

Optimal array of sand fences.

Sand fences are widely applied to prevent soil erosion by wind in areas affected by desertification. Sand fences also provide a way to reduce the emission rate of dust particles, which is triggered mainly by the impacts of wind-blown sand grains onto the soil and affects the Earth's climate. Many different types of fence have been designed and their effects on the sediment transport dynamics studied since many years. However, the search for the optimal array of fences has remained largely an empirical task. In order to achieve maximal soil protection using the minimal amount of fence material, a quantitative understanding of the flow profile over the relief encompassing the area to be protected including all employed fences is required. Here we use Computational Fluid Dynamics to calculate the average turbulent airflow through an array of fences as a function of the porosity, spacing and height of the fences. Specifically, we investigate the factors controlling the fraction of soil area over which the basal average wind shear velocity drops below the threshold for sand transport when the fences are applied. We introduce a cost function, given by the amount of material necessary to construct the fences. We find that, for typical sand-moving wind velocities, the optimal fence height (which minimizes this cost function) is around 50 cm, while using fences of height around 1.25 m leads to maximal cost.

Mesoscopic approach to subcritical fatigue crack growth.

We investigate a model for fatigue crack growth in which damage accumulation is assumed to follow a power law of the local stress amplitude, a form that can be generically justified on the grounds of the approximately self-similar aspect of microcrack distributions. Our aim is to determine the relation between model ingredients and the Paris exponent governing subcritical crack-growth dynamics at the macroscopic scale, starting from a single small notch propagating along a fixed line. By a series of analytical and numerical calculations, we show that, in the absence of disorder, there is a critical damage-accumulation exponent γ, namely γ_{c}=2, separating two distinct regimes of behavior for the Paris exponent m. For γ>γ_{c}, the Paris exponent is shown to assume the value m=γ, a result that proves robust against the separate introduction of various modifying ingredients. Explicitly, we deal here with (i) the requirement of a minimum stress for damage to occur, (ii) the presence of disorder in local damage thresholds, and (iii) the possibility of crack healing. On the other hand, in the regime γ<γ_{c}, the Paris exponent is seen to be sensitive to the different ingredients added to the model, with rapid healing or a high minimum stress for damage leading to m=2 for all γ<γ_{c}, in contrast with the linear dependence m=6-2γ observed for very long characteristic healing times in the absence of a minimum stress for damage. Upon the introduction of disorder on the local fatigue thresholds, which leads to the possible appearance of multiple cracks along the propagation line, the Paris exponent tends to m≈4 for γ≲2 while retaining the behavior m=γ for γ≳4.

Stochastic Loewner evolution relates anomalous diffusion and anisotropic percolation.

We disclose the origin of anisotropic percolation perimeters in terms of the stochastic Loewner evolution (SLE) process. Precisely, our results from extensive numerical simulations indicate that the perimeters of multilayered and directed percolation clusters at criticality are the scaling limits of the Loewner evolution of an anomalous Brownian motion, being superdiffusive and subdiffusive, respectively. The connection between anomalous diffusion and fractal anisotropy is further tested by using long-range power-law correlated time series (fractional Brownian motion) as the driving functions in the evolution process. The fact that the resulting traces are distinctively anisotropic corroborates our hypothesis. Under the conceptual framework of SLE, our study therefore reveals different perspectives for mathematical and physical interpretations of non-Markovian processes in terms of anisotropic paths at criticality and vice versa.

Itinerant Conductance in Fuse-Antifuse Networks.

We report on a novel dynamic phase in electrical networks, in which current channels perpetually change in time. This occurs when the elementary units of the network are fuse-antifuse devices, namely, become insulators within a certain finite interval of local applied voltages. As a consequence, the macroscopic current exhibits temporal fluctuations which increase with system size. We determine the conditions under which this exotic situation appears by establishing a phase diagram as a function of the applied field and the size of the insulating window. Besides its obvious application as a versatile electronic device, due to its rich variety of behaviors, this network model provides a possible description for particle-laden flow through porous media leading to dynamical clogging and reopening of the local channels in the pore space.

Mandala networks: ultra-small-world and highly sparse graphs.

The increasing demands in security and reliability of infrastructures call for the optimal design of their embedded complex networks topologies. The following question then arises: what is the optimal layout to fulfill best all the demands? Here we present a general solution for this problem with scale-free networks, like the Internet and airline networks. Precisely, we disclose a way to systematically construct networks which are robust against random failures. Furthermore, as the size of the network increases, its shortest path becomes asymptotically invariant and the density of links goes to zero, making it ultra-small world and highly sparse, respectively. The first property is ideal for communication and navigation purposes, while the second is interesting economically. Finally, we show that some simple changes on the original network formulation can lead to an improved topology against malicious attacks.

Incipient speciation with biased gene flow between two lineages of the Western Diamondback Rattlesnake (Crotalus atrox).

We used mitochondrial DNA sequence data from 151 individuals to estimate population genetic structure across the range of the Western Diamondback Rattlesnake (Crotalus atrox), a widely distributed North American pit viper. We also tested hypotheses of population structure using double-digest restriction site associated DNA (ddRADseq) data, incorporating thousands of nuclear genome-wide SNPs from 42 individuals. We found strong mitochondrial support for a deep divergence between eastern and western C. atrox populations, and subsequent intermixing of these populations in the Inter-Pecos region of the United States and Mexico. Our nuclear RADseq data also identify these two distinct lineages of C. atrox, and provide evidence for nuclear admixture of eastern and western alleles across a broad geographic region. We identified contrasting patterns of mitochondrial and nuclear genetic variation across this genetic fusion zone that indicate partially restricted patterns of gene flow, which may be due to either pre- or post-zygotic isolating mechanisms. The failure of these two lineages to maintain complete genetic isolation, and evidence for partially-restricted gene flow, imply that these lineages were in the early stages of speciation prior to secondary contact.

Explosive electric breakdown due to conducting-particle deposition on an insulating substrate.

We introduce a theoretical model to investigate the electric breakdown of a substrate on which highly conducting particles are adsorbed and desorbed with a probability that depends on the local electric field. We find that, by tuning the relative strength q of this dependence, the breakdown can change from continuous to explosive. Precisely, in the limit in which the adsorption probability is the same for any finite voltage drop, we can map our model exactly onto the q-state Potts model and thus the transition to a jump occurs at q = 4. In another limit, where the adsorption probability becomes independent of the local field strength, the traditional bond percolation model is recovered. Our model is thus an example of a possible experimental realization exhibiting a truly discontinuous percolation transition.

Numerical modeling of the wind flow over a transverse dune.

Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee - the separation bubble - displays a surprisingly strong dependence on the wind shear velocity, u: it is nearly independent of u for shear velocities within the range between 0.2 m/s and 0.8 m/s but increases linearly with u for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if u is larger than approximately 0.39 m/s, whereas a larger value of u (about 0.49 m/s) is required to initiate this reverse transport.

Overspill avalanching in a dense reservoir network.

Sustainability of communities, agriculture, and industry is strongly dependent on an effective storage and supply of water resources. In some regions the economic growth has led to a level of water demand that can only be accomplished through efficient reservoir networks. Such infrastructures are not always planned at larger scale but rather made by farmers according to their local needs of irrigation during droughts. Based on extensive data from the upper Jaguaribe basin, one of the world's largest system of reservoirs, located in the Brazilian semiarid northeast, we reveal that surprisingly it self-organizes into a scale-free network exhibiting also a power-law in the distribution of the lakes and avalanches of discharges. With a new self-organized-criticality-type model we manage to explain the novel critical exponents. Implementing a flow model we are able to reproduce the measured overspill evolution providing a tool for catastrophe mitigation and future planning.

Transverse instability of dunes.

The simplest type of dune is the transverse one, which propagates with invariant profile orthogonally to a fixed wind direction. Here we show, by means of numerical simulations, that transverse dunes are unstable with respect to along-axis perturbations in their profile and decay on the bedrock into barchan dunes. Any forcing modulation amplifies exponentially with growth rate determined by the dune turnover time. We estimate the distance covered by a transverse dune before fully decaying into barchans and identify the patterns produced by different types of perturbation.

Oil displacement through a porous medium with a temperature gradient.

We investigate the effect of a temperature gradient on oil recovery in a two-dimensional pore-network model. The oil viscosity depends on temperature as µ(o) [Please see text] e(B/T), where B is a physicochemical parameter, depending on the type of oil, and T is the temperature. A temperature gradient is applied across the medium in the flow direction. Initially, the porous medium is saturated with oil, and then another fluid is injected. We have considered two cases representing different injection strategies. In the first case, the invading fluid viscosity is constant (finite viscosity ratio), while in the second one, the invading fluid is inviscid (infinite viscosity ratio). Our results show that for the case of finite viscosity ratio, recovery increases with ΔT independent of strength or sign of the gradient. For an infinite viscosity ratio, a positive temperature gradient is necessary to enhance recovery. Moreover, we show that for ΔT>0, the percentage of oil recovery generally decreases (increases) with B for a finite (infinite) viscosity ratio. Finally, we also extend our results for infinite viscosity ratio to a three-dimensional porous media geometry.

Transport on exploding percolation clusters.

We propose a simple generalization of the explosive percolation process [Achlioptas et al., Science 323, 1453 (2009)], and investigate its structural and transport properties. In this model, at each step, a set of q unoccupied bonds is randomly chosen. Each of these bonds is then associated with a weight given by the product of the cluster sizes that they would potentially connect, and only that bond among the q set which has the smallest weight becomes occupied. Our results indicate that, at criticality, all finite-size scaling exponents for the spanning cluster, the conducting backbone, the cutting bonds, and the global conductance of the system, change continuously and significantly with q. Surprisingly, we also observe that systems with intermediate values of q display the worst conductive performance. This is explained by the strong inhibition of loops in the spanning cluster, resulting in a substantially smaller associated conducting backbone.

Tactical voting in plurality elections.

How often will elections end in landslides? What is the probability for a head-to-head race? Analyzing ballot results from several large countries rather anomalous and yet unexplained distributions have been observed. We identify tactical voting as the driving ingredient for the anomalies and introduce a model to study its effect on plurality elections, characterized by the relative strength of the feedback from polls and the pairwise interaction between individuals in the society. With this model it becomes possible to explain the polarization of votes between two candidates, understand the small margin of victories frequently observed for different elections, and analyze the polls' impact in American, Canadian, and Brazilian ballots. Moreover, the model reproduces, quantitatively, the distribution of votes obtained in the Brazilian mayor elections with two, three, and four candidates.

Dune formation under bimodal winds.

The study of dune morphology represents a valuable tool in the investigation of planetary wind systems--the primary factor controlling the dune shape is the wind directionality. However, our understanding of dune formation is still limited to the simplest situation of unidirectional winds: There is no model that solves the equations of sand transport under the most common situation of seasonally varying wind directions. Here we present the calculation of sand transport under bimodal winds using a dune model that is extended to account for more than one wind direction. Our calculations show that dunes align longitudinally to the resultant wind trend if the angle(w) between the wind directions is larger than 90 degrees. Under high sand availability, linear seif dunes are obtained, the intriguing meandering shape of which is found to be controlled by the dune height and by the time the wind lasts at each one of the two wind directions. Unusual dune shapes including the "wedge dunes" observed on Mars appear within a wide spectrum of bimodal dune morphologies under low sand availability.