Fracturing ranked surfaces

Breaking universality in random sequential adsorption on a square lattice with long-range correlated defects

Jamming and percolation transitions in the standard random sequential adsorption of particles on regular lattices are characterized by a universal set of critical exponents. The universality class is preserved even in the presence of randomly distributed defective sites that are forbidden for particle deposition. However, using large-scale Monte Carlo simulations by depositing dimers on the square lattice and employing finite-size scaling, we provide evidence that the system does not exhibit such well-known universal features when the defects have spatial long-range (power-law) correlations. The critical exponents ν_{j} and ν associated with the jamming and percolation transitions, respectively, are found to be nonuniversal for strong spatial correlations and approach systematically their own universal values as the correlation strength is decreased. More crucially, we have found a difference in the values of the percolation correlation length exponent ν for a small but finite density of defects with strong spatial correlations. Furthermore, for a fixed defect density, it is found that the percolation threshold of the system, at which the largest cluster of absorbed dimers first establishes the global connectivity, gets reduced with increasing the strength of the spatial correlation.

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.

Occluding junctions as novel regulators of tissue mechanics during wound repair

Simple epithelial repair is mediated by the contraction of an actomyosin cable and cellular rearrangements at the wound edge. Carvalho et al. show that occluding junctions are required for epithelial repair by regulating these cellular rearrangements and tissue mechanical properties.

In epithelial tissues, cells tightly connect to each other through cell–cell junctions, but they also present the remarkable capacity of reorganizing themselves without compromising tissue integrity. Upon injury, simple epithelia efficiently resolve small lesions through the action of actin cytoskeleton contractile structures at the wound edge and cellular rearrangements. However, the underlying mechanisms and how they cooperate are still poorly understood. In this study, we combine live imaging and theoretical modeling to reveal a novel and indispensable role for occluding junctions (OJs) in this process. We demonstrate that OJ loss of function leads to defects in wound-closure dynamics: instead of contracting, wounds dramatically increase their area. OJ mutants exhibit phenotypes in cell shape, cellular rearrangements, and mechanical properties as well as in actin cytoskeleton dynamics at the wound edge. We propose that OJs are essential for wound closure by impacting on epithelial mechanics at the tissue level, which in turn is crucial for correct regulation of the cellular events occurring at the wound edge.

A comparison of hydrological and topological watersheds

We introduce the hydrological watershed, a watershed where water can penetrate the soil, and compare it with the topological watershed for a two-dimensional landscape. For this purpose, we measure the fractal dimension of the hydrological watershed for different penetration depths and different grid sizes. Through finite size scaling, we find that the fractal dimension is 1.31 ± 0.02 which is significantly higher than the fractal dimension of the topological watershed. This indicates that the hydrological watershed belongs to a new universality class. We also find that, as opposed to the topological watershed, the hydrodynamic watershed can exhibit disconnected islands.

Schramm-Loewner evolution and perimeter of percolation clusters of correlated random landscapes

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 ∈ [-1, 0] 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.

Abnormal grain growth mediated by fractal boundary migration at the nanoscale

Modern engineered materials are composed of space-filling grains or domains separated by a network of interfaces or boundaries. Such polycrystalline microstructures have the capacity to coarsen through boundary migration. Grain growth theories account for the topology of grains and the connectivity of the boundary network in terms of the familiar Euclidian dimension and Euler’s polyhedral formula, both of which are based on integer numbers. However, we recently discovered an unusual growth mode in a nanocrystalline Pd-Au alloy, in which grains develop complex, highly convoluted surface morphologies that are best described by a fractional dimension of ∼1.2 (extracted from the perimeters of grain cross sections). This fractal value is characteristic of a variety of domain growth scenarios—including explosive percolation, watersheds of random landscapes, and the migration of domain walls in a random field of pinning centers—which suggests that fractal grain boundary migration could be a manifestation of the same universal behavior.

The influence of statistical properties of Fourier coefficients on random Gaussian surfaces

Many examples of natural systems can be described by random Gaussian surfaces. Much can be learned by analyzing the Fourier expansion of the surfaces, from which it is possible to determine the corresponding Hurst exponent and consequently establish the presence of scale invariance. We show that this symmetry is not affected by the distribution of the modulus of the Fourier coefficients. Furthermore, we investigate the role of the Fourier phases of random surfaces. In particular, we show how the surface is affected by a non-uniform distribution of phases.

Critical Fragmentation Properties of Random Drilling: How Many Holes Need to Be Drilled to Collapse a Wooden Cube?

A solid wooden cube fragments into pieces as we sequentially drill holes through it randomly. This seemingly straightforward observation encompasses deep and nontrivial geometrical and probabilistic behavior that is discussed here. Combining numerical simulations and rigorous results, we find off-critical scale-free behavior and a continuous transition at a critical density of holes that significantly differs from classical percolation.

Two Types of Discontinuous Percolation Transitions in Cluster Merging Processes

Percolation is a paradigmatic model in disordered systems and has been applied to various natural phenomena. The percolation transition is known as one of the most robust continuous transitions. However, recent extensive studies have revealed that a few models exhibit a discontinuous percolation transition (DPT) in cluster merging processes. Unlike the case of continuous transitions, understanding the nature of discontinuous phase transitions requires a detailed study of the system at hand, which has not been undertaken yet for DPTs. Here we examine the cluster size distribution immediately before an abrupt increase in the order parameter of DPT models and find that DPTs induced by cluster merging kinetics can be classified into two types. Moreover, the type of DPT can be determined by the key characteristic of whether the cluster kinetic rule is homogeneous with respect to the cluster sizes. We also establish the necessary conditions for each type of DPT, which can be used effectively when the discontinuity of the order parameter is ambiguous, as in the explosive percolation model.

Loop-erased random walk on a percolation cluster is compatible with Schramm-Loewner evolution

We study the scaling limit of a planar loop-erased random walk (LERW) on the percolation cluster, with occupation probability p≥p(c). We numerically demonstrate that the scaling limit of planar LERW(p) curves, for all p>p(c), can be described by Schramm-Loewner evolution (SLE) with a single parameter κ that is close to the normal LERW in a Euclidean lattice. However, our results reveal that the LERW on critical incipient percolation clusters is compatible with SLE, but with another diffusivity coefficient κ. Several geometrical tests are applied to ascertain this. All calculations are consistent with SLE(κ), where κ=1.732±0.016. This value of the diffusivity coefficient is outside the well-known duality range 2≤κ≤8. We also investigate how the winding angle of the LERW(p) crosses over from Euclidean to fractal geometry by gradually decreasing the value of the parameter p from 1 to p(c). For finite systems, two crossover exponents and a scaling relation can be derived. This finding should, to some degree, help us understand and predict the existence of conformal invariance in disordered and fractal landscapes.

Loop-erased random walk on a percolation cluster: crossover from Euclidean to fractal geometry

We study loop-erased random walk (LERW) on the percolation cluster, with occupation probability p ≥ p_{c}, in two and three dimensions. We find that the fractal dimensions of LERW_{p} are close to normal LERW in a Euclidean lattice, for all p>p_{c}. However, our results reveal that LERW on critical incipient percolation clusters is fractal with d_{f}=1.217 ± 0.002 for d=2 and 1.43 ± 0.02 for d=3, independent of the coordination number of the lattice. These values are consistent with the known values for optimal path exponents in strongly disordered media. We investigate how the behavior of the LERW_{p} crosses over from Euclidean to fractal geometry by gradually decreasing the value of the parameter p from 1 to p_{c}. For finite systems, two crossover exponents and a scaling relation can be derived. This work opens up a theoretical window regarding the diffusion process on fractal and random landscapes.

Retention capacity of correlated surfaces

We extend the water retention model [C. L. Knecht et al., Phys. Rev. Lett. 108, 045703 (2012)] to correlated random surfaces. We find that the retention capacity of discrete random landscapes is strongly affected by spatial correlations among the heights. This phenomenon is related to the emergence of power-law scaling in the lake volume distribution. We also solve the uncorrelated case exactly for a small lattice and present bounds on the retention of uncorrelated landscapes.

Crossover behavior of conductivity in a discontinuous percolation model

When conducting bonds are occupied randomly in a two-dimensional square lattice, the conductivity of the system increases continuously as the density of those conducting bonds exceeds the percolation threshold. Such a behavior is well known in percolation theory; however, the conductivity behavior has not been studied yet when the percolation transition is discontinuous. Here we investigate the conductivity behavior through a discontinuous percolation model evolving under a suppressive external bias. Using effective medium theory, we analytically calculate the conductivity behavior as a function of the density of conducting bonds. The conductivity function exhibits a crossover behavior from a drastically to a smoothly increasing function beyond the percolation threshold in the thermodynamic limit. The analytic expression fits well our simulation data.

Percolation with long-range correlated disorder

Long-range power-law correlated percolation is investigated using Monte Carlo simulations. We obtain several static and dynamic critical exponents as functions of the Hurst exponent H, which characterizes the degree of spatial correlation among the occupation of sites. In particular, we study the fractal dimension of the largest cluster and the scaling behavior of the second moment of the cluster size distribution, as well as the complete and accessible perimeters of the largest cluster. Concerning the inner structure and transport properties of the largest cluster, we analyze its shortest path, backbone, red sites, and conductivity. Finally, bridge site growth is also considered. We propose expressions for the functional dependence of the critical exponents on H.

Percolation model with continuously varying exponents

This work analyzes a percolation model on the diamond hierarchical lattice (DHL), where the percolation transition is retarded by the inclusion of a probability of erasing specific connected structures. It has been inspired by the recent interest on the existence of other universality classes of percolation models. The exact scale invariance and renormalization properties of DHL leads to recurrence maps, from which analytical expressions for the critical exponents and precise numerical results in the limit of very large lattices can be derived. The critical exponents ν and β of the investigated model vary continuously as the erasing probability changes. An adequate choice of the erasing probability leads to the result ν=∞, like in some phase transitions involving vortex formation. The percolation transition is continuous, with β>0, but β can be as small as desired. The modified percolation model turns out to be equivalent to the Q→1 limit of a Potts model with specific long range interactions on the same lattice.

Exact evaluation of the cutting path length in a percolation model on a hierarchical network

This work presents an approach to evaluate the exact value of the fractal dimension of the cutting path d(f)(CP) on hierarchical structures with finite order of ramification. Our approach is based on a renormalization group treatment of the universality class of watersheds. By making use of the self-similar property, we show that d(f)(CP) depends only on the average cutting path (CP) of the first generation of the structure. For the simplest Wheastone hierarchical lattice (WHL), we present a mathematical proof. For a larger WHL structure, the exact value of d(f)(CP) is derived based on a computer algorithm that identifies the length of all possible CP's of the first generation.

Avoiding a spanning cluster in percolation models

When dynamics in a system proceeds under suppressive external bias, the system can undergo an abrupt phase transition, as can happen when an epidemic spreads. Recently, an explosive percolation (EP) model was introduced to understand such phenomena. The order of the EP transition has not been clarified in a unified framework covering low-dimensional systems and the mean-field limit. We introduce a stochastic model in which a rule for dynamics is designed to avoid the formation of a spanning cluster through competitive selection in Euclidean space. We use heuristic arguments to show that in the thermodynamic limit and depending on a control parameter, the EP transition can be either continuous or discontinuous if d < d(c) and is always continuous if d ≥ d(c), where d(c) is the spatial dimension and d is the upper critical dimension.

Fracturing highly disordered materials

We investigate the role of disorder on the fracturing process of heterogeneous materials by means of a two-dimensional fuse network model. Our results in the extreme disorder limit reveal that the backbone of the fracture at collapse, namely, the subset of the largest fracture that effectively halts the global current, has a fractal dimension of 1.22 ± 0.01. This exponent value is compatible with the universality class of several other physical models, including optimal paths under strong disorder, disordered polymers, watersheds and optimal path cracks on uncorrelated substrates, hulls of explosive percolation clusters, and strands of invasion percolation fronts. Moreover, we find that the fractal dimension of the largest fracture under extreme disorder, d(f) = 1.86 ± 0.01, is outside the statistical error bar of standard percolation. This discrepancy is due to the appearance of trapped regions or cavities of all sizes that remain intact till the entire collapse of the fuse network, but are always accessible in the case of standard percolation. Finally, we quantify the role of disorder on the structure of the largest cluster, as well as on the backbone of the fracture, in terms of a distinctive transition from weak to strong disorder characterized by a new crossover exponent.

Watersheds are Schramm-Loewner evolution curves

We show that in the continuum limit watersheds dividing drainage basins are Schramm-Loewner evolution (SLE) curves, being described by one single parameter κ. Several numerical evaluations are applied to ascertain this. All calculations are consistent with SLE(κ), with κ = 1.734 ± 0.005, being the only known physical example of an SLE with κ<2. This lies outside the well-known duality conjecture, bringing up new questions regarding the existence and reversibility of dual models. Furthermore, it constitutes a strong indication for conformal invariance in random landscapes and suggests that watersheds likely correspond to a logarithmic conformal field theory with a central charge c ≈ -7/2.

How to share underground reservoirs

Many resources, such as oil, gas, or water, are extracted from porous soils and their exploration is often shared among different companies or nations. We show that the effective shares can be obtained by invading the porous medium simultaneously with various fluids. Partitioning a volume in two parts requires one division surface while the simultaneous boundary between three parts consists of lines. We identify and characterize these lines, showing that they form a fractal set consisting of a single thread spanning the medium and a surrounding cloud of loops. While the spanning thread has fractal dimension 1.55 ± 0.03, the set of all lines has dimension 1.69 ± 0.02. The size distribution of the loops follows a power law and the evolution of the set of lines exhibits a tricritical point described by a crossover with a negative dimension at criticality.

Corrections to scaling for watersheds, optimal path cracks, and bridge lines

We study the corrections to scaling for the mass of the watershed, the bridge line, and the optimal path crack in two and three dimensions (2D and 3D). We disclose that these models have numerically equivalent fractal dimensions and leading correction-to-scaling exponents. We conjecture all three models to possess the same fractal dimension, namely, d(f) =1.2168 ± 0.0005 in 2D and d(f) = 2.487 ± 0.003 in 3D, and the same exponent of the leading correction, Ω = 0.9 ± 0.1 and Ω=1.0 ± 0.1, respectively. The close relations between watersheds, optimal path cracks in the strong disorder limit, and bridge lines are further supported by either heuristic or exact arguments.