Search

Mandelbrot's empirical observation that the coast of Britain is fractal has been confirmed by many authors, but it can be described by the Schramm--Loewner Evolution? Since the self-affine surface of our planet has a positive Hurst exponent, one would not expect a priori any critical behavior. Here, we investigate numerically the roughness and fractal dimension of the isoheight lines of real and artificial landscapes. Using a novel algorithm to take into account overhangs, we find that the roughness exponent of isoheight lines is consistent with unity regardless of the Hurst exponent of the rough surface. Moreover, the effective fractal dimension of the iso-height lines decays linearly with the Hurst exponent of the surface. We perform several tests to verify if the complete and accessible perimeters would follow the Schramm--Loewner Evolution and find that the left passage probability test is clearly violated, implying that coastlines violate SLE., Comment: 17 pages, 11 figures

We study the average long-time behavior of the binary opinions of a social group with peer-to-peer interactions under the influence of an external bias and a persuadable leader, a strongly-biased agent with a dynamic opinion with the intention of spreading it across the system. We use a generalized, fully-connected Ising model, with each spin representing the binary opinion of an agent at a given time and a single, super spin representing the opinion of the leader. External fields and interaction constants model the opinion bias and peer-to-peer interactions, respectively, while the temperature $T$ models an idealized social climate, representing an authoritarian regime if $T$ is low or a liberal one if $T$ is high. We derive a mean-field solution for the average magnetization $m$, the "social mood", and investigate how $m$ and the super spin magnetization vary as a function of $T$. We find that, depending on the initial conditions, due to the presence of metastable states, the sign of the average magnetization depends on the temperature. Finally, we verify that this effect is also present even if we consider only nearest-neighbor interactions within the social group.

We scrutinize the phenomenology arising from a minimal vaccination-epidemic (MVE) dynamics using three methods: mean-field approach, Monte Carlo simulations, and finite-size scaling analysis. The mean-field formulation reveals that the MVE model exhibits either a continuous or a discontinuous active-to-absorbing phase transition, accompanied by bistability and a tricritical point. However, on square lattices, we detect no signs of bistability, and we disclose that the active-to-absorbing state transition has a scaling invariance and critical exponents compatible with the continuous transition of the directed percolation universality class. Additionally, our findings indicate that the tricritical and crossover behaviors of the MVE dynamics belong to the universality class of mean-field tricritical directed percolation., Comment: 12 pages, 6 figures and 2 tables

We use the configuration model to generate networks having a degree distribution that follows a $q$-exponential, $P_q(k)=(2-q)\lambda[1-(1-q)\lambda k]^{1/(q-1)}$, for arbitrary values of the parameters $q$ and $\lambda$. We study the assortativity and the shortest path of these networks finding that the more the distribution resembles a pure power law, the less well connected are the corresponding nodes. In fact, the average degree of a nearest neighbor grows monotonically with $\lambda^{-1}$. Moreover, our results show that $q$-exponential networks are more robust against random failures and against malicious attacks than standard scale-free networks. Indeed, the critical fraction of removed nodes grows logarithmically with $\lambda^{-1}$ for malicious attacks. An analysis of the $k_s$-core decomposition shows that $q$-exponential networks have a highest $k_s$-core, that is bigger and has a larger $k_s$ than pure scale-free networks. Being at the same time well connected and robust, networks with $q$-exponential degree distribution exhibit scale-free and small-world properties, making them a particularly suitable model for application in several systems., Comment: 6 pages, 8 figures

Many systems in Nature exhibit avalanche dynamics with scale-free features. A general scaling theory has been proposed for critical avalanche profiles in crackling noise, predicting the collapse onto a universal avalanche shape, as well as the scaling behaviour of the activity power spectrum as Brown noise. Recently, much attention has been given to the profile of neuronal avalanches, measured in neuronal systems in vitro and in vivo. Although a universal profile was evidenced, confirming the validity of the general scaling theory, the parallel study of the power spectrum scaling under the same conditions was not performed. The puzzling observation is that in the majority of healthy neuronal systems the power spectrum exhibits a behaviour close to $1/f$, rather than Brown, noise. Here we perform a numerical study of the scaling behaviour of avalanche shape and power spectrum for a model of integrate and fire neurons with a short-term plasticity parameter able to tune the system to criticality. We confirm that, at criticality, the average avalanche size and the avalanche profile fulfill the general avalanche scaling theory. However, the power spectrum consistently exhibits Brown noise behaviour, for both fully excitatory networks and systems with 30\% inhibitory networks. Conversely, a behaviour closer to $1/f$ noise is observed in systems slightly off-criticality. Results suggest that the power spectrum is a good indicator to determine how close neuronal activity is to criticality., Comment: 9 pages, 7 figures

The two-dimensional Loewner exploration process is generalized to the case where the random force is self-similar with positively correlated increments. We model this random force by a fractional Brownian motion with Hurst exponent $H\geq \frac{1}{2}\equiv H_{\text{BM}}$, where $H_{\text{BM}}$ stands for the one-dimensional Brownian motion. By manipulating the deterministic force, we design a scale-invariant equation describing self-similar traces which lack conformal invariance. The model is investigated in terms of the "input diffusivity parameter" $\kappa$, which coincides with the one of the ordinary Schramm-Loewner evolution (SLE) at $H=H_{\text{BM}}$. In our numerical investigation, we focus on the scaling properties of the traces generated for $\kappa=2,3$, $\kappa=4$ and $\kappa=6,8$ as the representatives, respectively, of the dilute phase, the transition point and the dense phase of the ordinary SLE. The resulting traces are shown to be scale-invariant. Using two equivalent schemes, we extract the fractal dimension, $D_f(H)$, of the traces which decrease monotonically with increasing $H$, reaching $D_f=1$ at $H=1$ for all $\kappa$ values. The left passage probability (LPP) test demonstrates that, for $H$ values not far from the uncorrelated case (small $\epsilon_H\equiv \frac{H-H_{\text{BM}}}{H_{\text{BM}}}$) the prediction of the ordinary SLE is applicable with an effective diffusivity parameter $\kappa_{\text{eff}}$. Not surprisingly, the $\kappa_{\text{eff}}$'s do not fulfill the prediction of SLE for the relation between $D_f(H)$ and the diffusivity parameter.

Magnetic beads attract each other forming chains. We pushed such chains into an inclined Hele-Shaw cell and discovered 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., Comment: 5 pages, 6 figures

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., Comment: To appear in Physical Review E

Molecular dynamics simulations of a three dimensional relativistic gas with a soft potential are conducted with different interactions and particle masses. For all cases the velocity distribution agrees numerically with the J\"uttner distribution. We show how the relativistic gas can be coupled to a thermostat to simulate the canonical ensemble at a given temperature T. The behaviour of the thermostat is investigated as a function of the thermal inertia and its appropriate range is determined by evaluating the kinetic energy fluctuations.

A Monte Carlo algorithm to simulate the isothermal recrystallization process of snow is presented. The snow metamorphism is approximated by two mass redistribution processes, surface diffusion and sublimation-deposition. The algorithm is justified and its parametrization is determined. The simulation results are compared to experimental data, in particular, the temporal evolution of the specific surface area and the ice thickness. We find that the two effects of surface diffusion and sublimation-deposition can accurately model many aspects of the isothermal metamorphism of snow. Furthermore, it is shown that sublimation-deposition is the dominant contribution for temperatures close to the melting point, whereas surface diffusion dominates at temperatures far below the melting point. A simple approximation of gravitational compaction is implemented to simulate density change., Comment: 6 pages, 5 figures, 2 tables

After their re-emergence in the last decades, dengue fever and other vector-borne diseases are a potential threat to the lives of millions of people. Based on a data set of dengue cases in the Brazilian city of Fortaleza, collected from 2011 to 2016, we study the spatio-temporal characteristics of dengue outbreaks to characterize epidemic and non-epidemic years. First, we identify regions that show a high prevalence of dengue cases and mosquito larvae in different years and also analyze their corresponding correlations. Our results show that the characteristic correlation length of the epidemic is of the order of the system size, suggesting that factors such as citizen mobility may play a major role as a drive for spatial spreading of vector-borne diseases. Inspired by this observation, we perform a mean-field estimation of the basic reproduction number and find that our estimated values agree well with the values reported for other regions, pointing towards similar underlying spreading mechanisms. These findings provide insights into the spreading characteristics of dengue in densely populated areas and should be of relevance for the design of improved disease containment strategies., Comment: 25 pages, 8 figures

We propose a method for the simulation of particle fragmentation based on the calculation of the energy landscape inside the particle. The landscape of strain energy is calculated in terms of internal stress using the principles of damage and fracture mechanics. Numerical calculation of the landscape s ridges is used to determine the breakage criterion as well as the shape of the postbreakage fragments. This method provides a physical-based understanding of the breakage effect in granular material.

Stroke is one of the main causes of human disabilities. Experimental observations indicate that several mechanisms are activated during the recovery of functional activity after a stroke. Here we unveil how the brain recovers by explaining the role played by three mechanisms: Plastic adaptation, hyperexcitability and synaptogenesis. We consider two different damages in a neural network: A diffuse damage that simply causes the reduction of the effective system size and a localized damage, a stroke, that strongly alters the spontaneous activity of the system. Recovery mechanisms observed experimentally are implemented both separately and in a combined way. Interestingly, each mechanism contributes to the recovery to a limited extent. Only the combined application of all three together is able to recover the spontaneous activity of the undamaged system. This explains why the brain triggers independent mechanisms, whose cooperation is the fundamental ingredient for the system recovery.

Recent political campaigns have demonstrated how technologies are used to boost election outcomes by microtargeting voters. We propose and analyze a framework which analyzes how political activists use technologies to target voters. Voters are represented as nodes of a network. Political activists reach out locally to voters and try to convince them. Depending on their technological advantage and budget, political activists target certain regions in the network where their activities are able to generate the largest vote-share gains. Analytically and numerically, we quantify vote-share gains and savings in terms of budget and number of activists from employing superior targeting technologies compared to traditional campaigns. Moreover, we demonstrate that the technological precision must surpass a certain threshold in order to lead to a vote-share gain or budget advantage. Finally, by calibrating the technology parameters to the recent U.S. presidential election, we show that a pure targeting technology advantage is consistent with Trump winning against Clinton., Comment: 25 pages, 7 figures

We investigate the effect of sequentiallydisrupting the shortest path of percolation clusters at criticality by comparing it with the shortest alternative path. We measure the difference in length and the enclosed area between the two paths. The sequential approach allows to study spatial correlations. We find the lengths of the segments of successively constant differences in length to be uncorrelated. Simultaneously, we study the distance between red bonds. We find the probability distributions for the enclosed areas A, the differences in length $\Delta l$, and the lengths between the redbonds $l_r$ to follow power law distributions. Using maximum likelihood estimation and extrapolation we find the exponents $\beta$ = 1.38 $\pm$ 0.03 for $\Delta l$, $\alpha$ = 1.186 $\pm$ 0.008 for A and $\delta$ = 1.64 $\pm$ 0.025 for thedistribution of $l_r$.

Based on the numerical solution of the quantum lattice Boltzmann method in curved space, we predict the onset of a quantized alternating current on curved graphene sheets. Such numerical prediction is verified analytically via a set of semi-classical equations relating the Berry curvature to real space curvature. The proposed quantised oscillating current on curved graphene could form the basis for the implementation of quantum information processing algorithms., Comment: 7 pages 2 figures. arXiv admin note: substantial text overlap with arXiv:1902.03892

We report on the buckling and subsequent collapse of orthotropic elastic spherical shells under volume and pressure control. Going far beyond what is known for isotropic shells, a rich morphological phase space with three distinct regimes emerges upon variation of shell slenderness and degree of orthotropy. Our extensive numerical simulations are in agreement with experiments using fabricated polymer shells. The shell buckling pathways and corresponding strain energy evolution are shown to depend strongly on material orthotropy. We find surprisingly robust orthotropic structures with strong similarities to stomatocytes and tricolpate pollen grains, suggesting that the shape of several of Nature's collapsed shells could be understood from the viewpoint of material orthotropy., Comment: 7 pages, 5 figures

The quantum Hall effect has been predicted and discovered in various condensed-matter systems. A promising quantum material for such topological effects is graphene. We report the numerical observation of a curvature-induced spin-Hall effect in a monolayer graphene M\"obius strip. The solution of the Dirac equation on the nontrivial and non-Euclidean manifold reveals that despite the absence of a Hall current, a spin-Hall current is a natural consequence for such a topology, as predicted from symmetry arguments., Comment: arXiv admin note: text overlap with arXiv:1810.02102

We investigate the mechanical behavior of a confined granular packing of irregular polyhedral particles under repeated heating and cooling cycles by means of numerical simulations with the Non-Smooth Contact Dynamics method. Assuming a homogeneous temperature distribution as well as constant temperature rate, we study the effect of the container shape, and coefficients of thermal expansions on the pressure buildup at the confining walls and the density evolution. We observe that small changes in the opening angle of the confinement can lead to a drastic peak pressure reduction. Furthermore, the displacement fields over several thermal cycles are obtained and we discover the formation of convection cells inside the granular material having the shape of a torus. The root mean square of the vorticity is then calculated from the displacement fields and a quadratic dependency on the ratio of thermal expansion coefficients is established.

We study the fragment size distributions after crushing of single and many particles under uniaxial compression inside a cylindrical container by means of numerical simulations. Under the assumption that breaking goes through the bulk of the particle we obtain the size distributions of fragments for both cases after large displacements. For the single particle crushing, this fragmentation mechanism produces a log-normal size distribution, which deviates from the power-law distribution of fragment sizes for the packed bed. We show that as the breaking process evolves, a power-law dependency on the displacement is present for the single grain, while for the many grains system, the distribution converges to a steady state. We further investigate the force networks and the average coordination number as a function of the particle size, which gives inside about the origin of the power-law distributions for the granular assembly under uniaxial compression., Comment: 12 figures

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, $\gamma=D/2$. We measure $\gamma=0.54 \pm 0.01$ for Earth, which is close to our estimation of $\gamma \approx 0.55$. Our study suggests that Hack's law can have its origin purely in the maximum and minimum lines of the landscapes., Comment: 20 pages, 6 Figures, and 1 Table

We study the discrete energy spectrum of curved graphene sheets in the presence of a magnetic field. The shifting of the Landau levels is determined for complex and realistic geometries of curved graphene sheets. The energy levels follow a similar square root dependence on the energy quantum number as for rippled and flat graphene sheets. The Landau levels are shifted towards lower energies proportionally to the average deformation and the effect is larger compared to a simple uni-axially rippled geometry. Furthermore, the resistivity of wrinkled graphene sheets is calculated for different average space curvatures and shown to obey a linear relation. The study is carried out with a quantum lattice Boltzmann method, solving the Dirac equation on curved manifolds., Comment: 6 pages, 4 figures, 27th International Conference on Discrete Simulation of Fluid Dynamics

The existence of power-law distributions is only a first requirement in the validation of the critical behavior of a system. Long-range spatio-temporal correlations are fundamental for the spontaneous neuronal activity to be the expression of a system acting close to a critical point. This chapter focuses on temporal correlations and avalanche dynamics in the spontaneous activity of cortex slice cultures and in the resting fMRI BOLD signal. Long-range correlations are investigated by means of the scaling of power spectra and of Detrended Fluctuations Analysis. The existence of 1/f decay in the power spectrum, as well as of power-law scaling in the root mean square fluctuations function for the appropriate balance of excitation and inhibition suggests that long-range temporal correlations are distinctive of "healthy brains". The corresponding temporal organization of neuronal avalanches can be dissected by analyzing the distribution of inter-event times between successive events. In rat cortex slice cultures this distribution exhibits a non-monotonic behavior, not usually found in other natural processes. Numerical simulations provide evidences that this behavior is a consequence of the alternation between states of high and low activity, leading to a dynamic balance between excitation and inhibition that tunes the system at criticality. In this scenario, inter-times show a peculiar relation with avalanche sizes, resulting in a hierarchical structure of avalanche sequences. Large avalanches correspond to low-frequency oscillations, and trigger cascades of smaller avalanches that are part of higher frequency rhythms. The self-regulated balance of excitation and inhibition observed in cultures is confirmed at larger scales, i.e. on fMRI data from resting brain activity, and appears to be closely related to critical features of avalanche activity., Comment: Conference paper

Dirac particles have been notoriously difficult to confine. Implementing a curved space Dirac equation solver based on the quantum Lattice Boltzmann method, we show that curvature in a 2-D space can confine a portion of a charged, mass-less Dirac fermion wave-packet. This is equivalent to a finite probability of confining the Dirac fermion within a curved space region. We propose a general power law expression for the probability of confinement with respect to average spatial curvature for the studied geometry., Comment: 10 pages 8 figures

Sand fences are often erected to reduce wind speed, prevent aeolian soil erosion, and induce sand deposition and dune formation in areas affected by sand encroachment and desertification. However, the search for the most efficient array of fences by means of field experiments alone poses a challenging task given that field experiments are affected by weather conditions. Here we apply Computational Fluid Dynamic simulations to investigate the three-dimensional structure of the turbulent wind field over an array of fences of different sizes, porosity and spacing. Our goal is to perform a quantitative analysis of this structure in the absence of saltation or suspension transport, as first step toward the development of a continuum simulation of the aeolian soil in presence of the fences. We find that the area of soil protected against direct aerodynamic entrainment has two regimes, depending on the spacing $L_x$ between the fences. When $L_x$ is smaller than a critical value $L_{xc}$, the wake zones associated with each fence are inter-connected (regime A), while these wake zones appear separated from each other (regime B) when $L_x$ exceeds this critical value of spacing. The system undergoes a second order phase transition at $L_x = L_{xc}$, with the cross-wind width of the protected zone scaling with $\left[{1-{L}_x}/{L}_{xc}\right]^\beta$ in regime A, with $\beta \approx 0.32$. Our findings have implication for a better understanding of aeolian transport in the presence of sand fences, as well as to develop optimization strategies for measures to protect soils from wind erosion., Comment: 35 pages, 13 figures

We investigate the behaviour of the shortest path on a directed two-dimensional square lattice for bond percolation at the critical probability $p_c$ . We observe that flipping an edge lying on the shortest path has a non-local effect in the form of power-law distributions for both the differences in shortest path lengths and for the minimal enclosed areas. Using maximum likelihood estimation and extrapolation we find the exponents $\alpha = 1.36 \pm 0.01$ for the path length differences and $\beta = 1.186 \pm 0.001$ for the enclosed areas.

We investigate percolation on a randomly directed lattice, an intermediate between standard percolation and directed percolation, focusing on the isotropic case in which bonds on opposite directions occur with the same probability. We derive exact results for the percolation threshold on planar lattices, and present a conjecture for the value the percolation-threshold for in any lattice. We also identify presumably universal critical exponents, including a fractal dimension, associated with the strongly-connected components both for planar and cubic lattices. These critical exponents are different from those associated either with standard percolation or with directed percolation., Comment: 10 pages, 11 figures

Introducing a model to study deposition and erosion of single particles at microscopic scale, we investigate the clogging and erosive processes in a pore. The particle diameter, concentration, and adhesive forces rule the way particles are deposited, and therefore, characterize the clogging process. We study the hydraulic pressure that induces erosive bursts and conclude that this pressure depends linearly on the deposited volume and inversely on the pores' diameter. While cohesion does not play an important role for erosive bursts, the adhesion is the main force initiating clogging and when overcome by the hydraulic pressure, erosive bursts are triggered. Finally, we show how the magnitude of erosive bursts depends on the pore length, particle diameter and pore size., Comment: more pictures than pages

We present a boundary-spheropolygon element method (BSEM), that combines the boundary integral method (BIM) and the spheropolygon-based discrete element method (SEM). The interaction between particles is simulated via the SEM, and the sub-particle stress (stress inside the grains) is calculated by BIM. The framework of BSEM is presented. Then the accuracy and efficiency of the method are analysed by comparison with both analytical solutions and a well-established finite element method (ABAQUS). The results demonstrate that BSEM could efficiently provide instant sub-particle stress for irregular particles with an optimized compromise between computational time and accuracy. The effect of particles aspect ratio, coordination number and heterogeneity on the sub-particle stress are discussed through parametric studies. Key conclusions on particle breakage are derived based on the analysis of the distribution of the sub-particle tensile stress. The simulation results suggest that BSEM could overcome most of the disadvantages of existing numerical methods and must be used for advanced simulations of particle breakage

We study two new models of two particle species invading a surface from opposite sides. Collisions of particles of different species lead to the formation of congestion fronts. One of the models implements a reversible process whereas in the other model the congestion front forms irreversibly. For both models we find that the congestion fronts are self-affine but with different roughness exponents. For low densities the system does not congest and we find a phase transition between a phase of freely moving particles and a congestion phase.

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\pm 0.003$, and one obtains a novel set of critical exponents $\beta_{eb} = 0.50\pm 0.02$, $\gamma_{eb} = 1.97\pm 0.05$, and $\nu_{eb} = 2.00\pm 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., Comment: 5 pages, 7 figures

We present a novel method for fluid structure interaction (FSI) simulations where an original 2nd-order curved space lattice Boltzmann fluid solver (LBM) is coupled to a finite element method (FEM) for thin shells. The LBM can work independently on a standard lattice in curved coordinates without the need for interpolation, re-meshing or an immersed boundary. The LBM distribution functions are transformed dynamically under coordinate change. In addition, force and momentum can be calculated on the nodes exactly in any geometry. Furthermore, the FEM shell is a complete numerical tool with implementations such as growth, self-contact and strong external forces. We show resolution convergent error for standard tests under metric deformation. Mass and volume conservation, momentum transfer, boundary-slip and pressure maintenance are verified through specific examples. Additionally, a brief deformation stability analysis is carried out. Next, we study the interaction of a square fluid flow channel to a deformable shell. Finally, we simulate a flag at moderate Reynolds number, air flow channel. The scheme is limited to small deformations of O(10%) relative to domain size, by improving its stability the method can be naturally extended to multiple applications without further implementations., Comment: 18 pages, 19 figures

How big is the risk that a few initial failures of nodes in a network amplify to large cascades that span a substantial share of all nodes? Predicting the final cascade size is critical to ensure the functioning of a system as a whole. Yet, this task is hampered by uncertain or changing parameters and missing information. In infinitely large networks, the average cascade size can often be well estimated by established approaches building on local tree approximations and mean field approximations. Yet, as we demonstrate, in finite networks, this average does not even need to be a likely outcome. Instead, we find broad and even bimodal cascade size distributions. This phenomenon persists for system sizes up to $10^{7}$ and different cascade models, i.e. it is relevant for most real systems. To show this, we derive explicit closed-form solutions for the full probability distribution of the final cascade size. We focus on two topological limit cases, the complete network representing a dense network with a very narrow degree distribution, and the star network representing a sparse network with a inhomogeneous degree distribution. Those topologies are of great interest, as they either minimize or maximize the average cascade size and are common motifs in many real world networks., Comment: systemic risk, finite size effects, cascades, networks

Political campaigns involve, in the simplest case, two competing campaign groups which try to obtain a majority of votes. We propose a novel mathematical framework to study political campaign dynamics on social networks whose constituents are either political activists or persuadable individuals. Activists are convinced and do not change their opinion and they are able to move around in the social network to motivate persuadable individuals to vote according to their opinion. We describe the influence of the complex interplay between the number of activists, political clout, budgets, and campaign costs on the campaign result. We also identify situations where the choice of one campaign group to send a certain number of activists already pre-determines their victory. Moreover, we show that a candidate's advantage in terms of political clout can overcome a substantial budget disadvantage or a lower number of activists, as illustrated by the US presidential election 2016., Comment: 23 pages, 10 figures

Motivated by the fact that many physical landscapes are characterized by long-range height-height correlations that are quantified by the Hurst exponent H, we investigate the statistical properties of the iso-height lines of correlated surfaces in the framework of Schramm-Loewner evolution (SLE). We show numerically that in the continuum limit the external perimeter of a percolating cluster of correlated surfaces with H between -1 and zero is statistically equivalent to SLE curves. Our results suggest that the external perimeter also retains the Markovian properties, confirmed by the absence of time correlations in the driving function and the fact that the latter is Gaussian distributed for any specific time. We also confirm that for all H the variance of the winding angle grows logarithmically with size.

We propose a targeted intervention protocol where recovery is restricted to individuals that have the least number of infected neighbours. Our recovery strategy is highly efficient on any kind of network, since epidemic outbreaks are minimal when compared to the baseline scenario of spontaneous recovery. In the case of spatially embedded networks, we find that an epidemic stays strongly spatially confined with a characteristic length scale undergoing a random walk. We demonstrate numerically and analytically that this dynamics leads to an epidemic spot with a flat surface structure and a radius that grows linearly with the spreading rate., Comment: 6 pages, 5 figures

In this paper, we develop and characterize the fully dissipative Lattice Boltzmann method for ultra-relativistic fluids in two dimensions using three equilibrium distribution functions: Maxwell-J\"uttner, Fermi-Dirac and Bose-Einstein. Our results stem from the expansion of these distribution functions up to fifth order in relativistic polynomials. We also obtain new Gaussian quadratures for square lattices that preserve the spatial resolution. Our models are validated with the Riemann problem and the limitations of lower order expansions to calculate higher order moments are shown. The kinematic viscosity and the thermal conductivity are numerically obtained using the Taylor-Green vortex and the Fourier flow respectively and these transport coefficients are compared with the theoretical prediction from Grad's theory. In order to compare different expansion orders, we analyze the temperature and heat flux fields on the time evolution of a hot spot.

Erosion and deposition during flow through porous media can lead to large erosive bursts that manifest as jumps in permeability and pressure loss. Here we reveal that the cause of these bursts is the re-opening of clogged pores when the pressure difference between two opposite sites of the pore surpasses a certain threshold. We perform numerical simulations of flow through porous media and compare our predictions to experimental results, recovering with excellent agreement shape and power-law distribution of pressure loss jumps, and the behavior of the permeability jumps as function of particle concentration. Furthermore, we find that erosive bursts only occur for pressure gradient thresholds within the range of two critical values, independent on how the flow is driven. Our findings provide a better understanding of sudden sand production in oil wells and breakthrough in filtration., Comment: 7 pages, 8 figures

Recent studies have proposed that the diffusion of messenger molecules, such as monoamines, can mediate the plastic adaptation of synapses in supervised learning of neural networks. Based on these findings we developed a model for neural learning, where the signal for plastic adaptation is assumed to propagate through the extracellular space. We investigate the conditions allowing learning of Boolean rules in a neural network. Even fully excitatory networks show very good learning performances. Moreover, the investigation of the plastic adaptation features optimizing the performance suggests that learning is very sensitive to the extent of the plastic adaptation and the spatial range of synaptic connections.

In particle-laden flows through porous media, porosity and permeability are significantly affected by the deposition and erosion of particles. Experiments show that the permeability evolution of a porous medium with respect to a particle suspension is not smooth, but rather exhibits significant jumps followed by longer periods of continuous permeability decrease. Their origin seems to be related to internal flow path reorganization by avalanches of deposited material due to erosion inside the porous medium. We apply neutron tomography to resolve the spatio-temporal evolution of the pore space during clogging and unclogging to prove the hypothesis of flow path reorganization behind the permeability jumps. This mechanistic understanding of clogging phenomena is relevant for a number of applications from oil production to filters or suffosion as the mechanisms behind sinkhole formation., Comment: 18 pages, 9 figures

When searching for a target within an image our brain can adopt different strategies, but which one does it choose? This question can be answered by tracking the motion of the eye while it executes the task. Following many individuals performing various search tasks we distinguish between two competing strategies. Motivated by these findings, we introduce a model that captures the interplay of the search strategies and allows us to create artificial eye-tracking trajectories, which could be compared to the experimental ones. Identifying the model parameters allows us to quantify the strategy employed in terms of ensemble averages, characterizing each experimental cohort. In this way we can discern with high sensitivity the relation between the visual landscape and the average strategy, disclosing how small variations in the image induce changes in the strategy.

The $1/f$-like decay observed in the power spectrum of electro-physiological signals, along with scale-free statistics of the so-called neuronal avalanches, constitute evidences of criticality in neuronal systems. Recent in vitro studies have shown that avalanche dynamics at criticality corresponds to some specific balance of excitation and inhibition, thus suggesting that this is a basic feature of the critical state of neuronal networks. In particular, a lack of inhibition significantly alters the temporal structure of the spontaneous avalanche activity and leads to an anomalous abundance of large avalanches. Here we study the relationship between network inhibition and the scaling exponent $\beta$ of the power spectral density (PSD) of avalanche activity in a neuronal network model inspired in Self-Organized Criticality (SOC). We find that this scaling exponent depends on the percentage of inhibitory synapses and tends to the value $\beta = 1$ for a percentage of about 30%. More specifically, $\beta$ is close to $2$, namely brownian noise, for purely excitatory networks and decreases towards values in the interval $[1,1.4]$ as the percentage of inhibitory synapses ranges between 20 and 30%, in agreement with experimental findings. These results indicate that the level of inhibition affects the frequency spectrum of resting brain activity and suggest the analysis of the PSD scaling behavior as a possible tool to study pathological conditions.

Particle detachment bursts during the flow of suspensions through porous media are a phenomenon that can severely affect the efficiency of deep bed filters. Despite the relevance in several industrial fields, little is known about the statistical properties and the temporal organization of these events. We present experiments of suspensions of deionized water carrying quartz particles pushed with a peristaltic pump through a filter of glass beads measuring simultaneously pressure drop, flux and suspension solid fraction. We find that the burst size distribution scales consistently with a power-law, suggesting that we are in the presence of a novel experimental realization of a self-organized critical system. Temporal correlations are present in the time series, alike in other phenomena as earthquakes or neuronal activity bursts, and also an analog to Omori's law can be shown. The understanding of bursts statistics could provide novel insights in different fields, e.g. in filter and petroleum industries., Comment: 7 pages, 9 figures

We provide numerical evidence that a Kelvin-Helmholtz instability occurs in the Dirac fluid of electrons in graphene and can be detected in current experiments. This instability appears for electrons in the viscous regime passing though a micrometer-scale obstacle and affects measurements on the time scale of nanoseconds. A possible realization with a needle-shaped obstacle is proposed to produce and detect this instability by measuring the electric potential difference between contact points located before and after the obstacle. We also show that, for our setup, the Kelvin-Helmholtz instability leads to the formation of whirlpools similar to the ones reported in Bandurin et al. [Science 351, 1055 (2016)]. To perform the simulations, we develop a lattice Boltzmann method able to recover the full dissipation in a fluid of massless particles., Comment: 12 pages, 12 figures

Magnetic nanoparticles in a colloidal solution self-assemble in various aligned structures, which has a profound influence on the flow behavior. However, the precise role of the microstructure in the development of the rheological response has not been reliably quantified. We investigate the self-assembly of dipolar colloids in simple shear using hybrid molecular dynamics and multi-particle collision dynamics simulations with explicit coarse-grained hydrodynamics; conduct simulated rheometric studies and apply micromechanical models to produce master curves, showing evidence of the universality of the structural behavior governed by the competition of the bonding (dipolar) and erosive (thermal and/or hydrodynamic) stresses. The simulations display viscosity changes across several orders of magnitude in fair quantitative agreement with various literature sources, substantiating the universality of the approach, which seems to apply generally across vastly different length scales and a broad range of physical systems.

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., Comment: 10 pages