Directed percolation

Last updated

In statistical physics, directed percolation (DP) refers to a class of models that mimic filtering of fluids through porous materials along a given direction, due to the effect of gravity. Varying the microscopic connectivity of the pores, these models display a phase transition from a macroscopically permeable (percolating) to an impermeable (non-percolating) state. Directed percolation is also used as a simple model for epidemic spreading with a transition between survival and extinction of the disease depending on the infection rate.

Contents

More generally, the term directed percolation stands for a universality class of continuous phase transitions which are characterized by the same type of collective behavior on large scales. Directed percolation is probably the simplest universality class of transitions out of thermal equilibrium.

Lattice models

Realization of permeable (solid line) or closed bonds (no line) on a tilted square lattice. The cluster of sites connected by a directed path (following the arrows) to the source (circled point) is indicated in red. Bond Directed Percolation.svg
Realization of permeable (solid line) or closed bonds (no line) on a tilted square lattice. The cluster of sites connected by a directed path (following the arrows) to the source (circled point) is indicated in red.

One of the simplest realizations of DP is bond directed percolation. This model is a directed variant of ordinary (isotropic) percolation and can be introduced as follows. The figure shows a tilted square lattice with bonds connecting neighboring sites. The bonds are permeable (open) with probability and impermeable (closed) otherwise. The sites and bonds may be interpreted as holes and randomly distributed channels of a porous medium.

The difference between ordinary and directed percolation is illustrated to the right. In isotropic percolation a spreading agent (e.g. water) introduced at a particular site percolates along open bonds, generating a cluster of wet sites. Contrarily, in directed percolation the spreading agent can pass open bonds only along a preferred direction in space, as indicated by the arrow. The resulting red cluster is directed in space.

As a dynamical process

Interpreting the preferred direction as a temporal degree of freedom, directed percolation can be regarded as a stochastic process that evolves in time. In a minimal, two-parameter model [1] that includes bond and site DP as special cases, a one-dimensional chain of sites evolves in discrete time , which can be viewed as a second dimension, and all sites are updated in parallel. Activating a certain site (called initial seed) at time the resulting cluster can be constructed row by row. The corresponding number of active sites varies as time evolves.

Universal scaling behavior

The DP universality class is characterized by a certain set of critical exponents. These exponents depend on the spatial dimension . Above the so-called upper critical dimension they are given by their mean-field values while in dimensions they have been estimated numerically. Current estimates are summarized in the following table:

Critical exponents of directed percolation in dimensions (September 2006) [2]
exponent

Other examples

In two dimensions, the percolation of water through a thin tissue (such as toilet paper) has the same mathematical underpinnings as the flow of electricity through two-dimensional random networks of resistors. In chemistry, chromatography can be understood with similar models.

The propagation of a tear or rip in a sheet of paper, in a sheet of metal, or even the formation of a crack in ceramic bears broad mathematical resemblance to the flow of electricity through a random network of electrical fuses. Above a certain critical point, the electrical flow will cause a fuse to pop, possibly leading to a cascade of failures, resembling the propagation of a crack or tear. The study of percolation helps indicate how the flow of electricity will redistribute itself in the fuse network, thus modeling which fuses are most likely to pop next, and how fast they will pop, and what direction the crack may curve in.

Examples can be found not only in physical phenomena, but also in biology, neuroscience, ecology (e.g. evolution), and economics (e.g. diffusion of innovation).

Percolation can be considered to be a branch of the study of dynamical systems or statistical mechanics. In particular, percolation networks exhibit a phase change around a critical threshold.

Experimental realizations

In spite of vast success in the theoretical and numerical studies of DP, obtaining convincing experimental evidence has proved challenging. In 1999 an experiment on flowing sand on an inclined plane was identified as a physical realization of DP. [3] In 2007, critical behavior of DP was finally found in the electrohydrodynamic convection of liquid crystal, where a complete set of static and dynamic critical exponents and universal scaling functions of DP were measured in the transition to spatiotemporal intermittency between two turbulent states. [4] [5]

See also

Sources

Literature

Related Research Articles

<span class="mw-page-title-main">Phase transition</span> Physical process of transition between basic states of matter

In chemistry, thermodynamics, and other related fields, a phase transition is the physical process of transition between one state of a medium and another. Commonly the term is used to refer to changes among the basic states of matter: solid, liquid, and gas, and in rare cases, plasma. A phase of a thermodynamic system and the states of matter have uniform physical properties. During a phase transition of a given medium, certain properties of the medium change as a result of the change of external conditions, such as temperature or pressure. This can be a discontinuous change; for example, a liquid may become gas upon heating to its boiling point, resulting in an abrupt change in volume. The identification of the external conditions at which a transformation occurs defines the phase transition point.

<span class="mw-page-title-main">Percolation theory</span> Mathematical theory on behavior of connected clusters in a random graph

In statistical physics and mathematics, percolation theory describes the behavior of a network when nodes or links are added. This is a geometric type of phase transition, since at a critical fraction of addition the network of small, disconnected clusters merge into significantly larger connected, so-called spanning clusters. The applications of percolation theory to materials science and in many other disciplines are discussed here and in the articles Network theory and Percolation.

<span class="mw-page-title-main">Percolation</span> Filtration of fluids through porous materials

In physics, chemistry, and materials science, percolation refers to the movement and filtering of fluids through porous materials. It is described by Darcy's law. Broader applications have since been developed that cover connectivity of many systems modeled as lattices or graphs, analogous to connectivity of lattice components in the filtration problem that modulates capacity for percolation.

In statistical mechanics, a universality class is a collection of mathematical models which share a single scale-invariant limit under the process of renormalization group flow. While the models within a class may differ dramatically at finite scales, their behavior will become increasingly similar as the limit scale is approached. In particular, asymptotic phenomena such as critical exponents will be the same for all models in the class.

<span class="mw-page-title-main">Self-organized criticality</span> Concept in physics

Self-organized criticality (SOC) is a property of dynamical systems that have a critical point as an attractor. Their macroscopic behavior thus displays the spatial or temporal scale-invariance characteristic of the critical point of a phase transition, but without the need to tune control parameters to a precise value, because the system, effectively, tunes itself as it evolves towards criticality.

<span class="mw-page-title-main">Lambda point</span>

The lambda point is the temperature at which normal fluid helium makes the transition to superfluid helium II. The lowest pressure at which He-I and He-II can coexist is the vapor−He-I−He-II triple point at 2.1768 K (−270.9732 °C) and 5.0418 kPa (0.049759 atm), which is the "saturated vapor pressure" at that temperature. The highest pressure at which He-I and He-II can coexist is the bcc−He-I−He-II triple point with a helium solid at 1.762 K (−271.388 °C), 29.725 atm (3,011.9 kPa).

Critical exponents describe the behavior of physical quantities near continuous phase transitions. It is believed, though not proven, that they are universal, i.e. they do not depend on the details of the physical system, but only on some of its general features. For instance, for ferromagnetic systems, the critical exponents depend only on:

<span class="mw-page-title-main">Percolation threshold</span> Threshold of percolation theory models

The percolation threshold is a mathematical concept in percolation theory that describes the formation of long-range connectivity in random systems. Below the threshold a giant connected component does not exist; while above it, there exists a giant component of the order of system size. In engineering and coffee making, percolation represents the flow of fluids through porous media, but in the mathematics and physics worlds it generally refers to simplified lattice models of random systems or networks (graphs), and the nature of the connectivity in them. The percolation threshold is the critical value of the occupation probability p, or more generally a critical surface for a group of parameters p1, p2, ..., such that infinite connectivity (percolation) first occurs.

The Binder parameter or Binder cumulant in statistical physics, also known as the fourth-order cumulant is defined as the kurtosis of the order parameter, s, introduced by Austrian theoretical physicist Kurt Binder. It is frequently used to determine accurately phase transition points in numerical simulations of various models.

In the context of the physical and mathematical theory of percolation, a percolation transition is characterized by a set of universal critical exponents, which describe the fractal properties of the percolating medium at large scales and sufficiently close to the transition. The exponents are universal in the sense that they only depend on the type of percolation model and on the space dimension. They are expected to not depend on microscopic details such as the lattice structure, or whether site or bond percolation is considered. This article deals with the critical exponents of random percolation.

<span class="mw-page-title-main">Active matter</span> Matter behavior at system scale

Active matter is matter composed of large numbers of active "agents", each of which consumes energy in order to move or to exert mechanical forces. Such systems are intrinsically out of thermal equilibrium. Unlike thermal systems relaxing towards equilibrium and systems with boundary conditions imposing steady currents, active matter systems break time reversal symmetry because energy is being continually dissipated by the individual constituents. Most examples of active matter are biological in origin and span all the scales of the living, from bacteria and self-organising bio-polymers such as microtubules and actin, to schools of fish and flocks of birds. However, a great deal of current experimental work is devoted to synthetic systems such as artificial self-propelled particles. Active matter is a relatively new material classification in soft matter: the most extensively studied model, the Vicsek model, dates from 1995.

<span class="mw-page-title-main">Epidemic models on lattices</span>

Classic epidemic models of disease transmission are described in Compartmental models in epidemiology. Here we discuss the behavior when such models are simulated on a lattice. Lattice models, which were first explored in the context of cellular automata, act as good first approximations of more complex spatial configurations, although they do not reflect the heterogeneity of space. Lattice-based epidemic models can also be implemented as fixed agent-based models.

The Ziff–Gulari–Barshad (ZGB) model is a simple Monte Carlo method for catalytic reactions of oxidation of carbon monoxide to carbon dioxide on a surface using Monte-Carlo methods which captures correctly the essential dynamics: the phase transition between two poisoned states (either CO2- or O-poisoned) and a steady-state in between. It is named after Robert M. Ziff, Erdogan Gulari, and Yoav Barshad, who published it in 1986.

The Vicsek model is a mathematical model used to describe active matter. One motivation of the study of active matter by physicists is the rich phenomenology associated to this field. Collective motion and swarming are among the most studied phenomena. Within the huge number of models that have been developed to catch such behavior from a microscopic description, the most famous is the model introduced by Tamás Vicsek et al. in 1995.

The Kibble–Zurek mechanism (KZM) describes the non-equilibrium dynamics and the formation of topological defects in a system which is driven through a continuous phase transition at finite rate. It is named after Tom W. B. Kibble, who pioneered the study of domain structure formation through cosmological phase transitions in the early universe, and Wojciech H. Zurek, who related the number of defects it creates to the critical exponents of the transition and to its rate—to how quickly the critical point is traversed.

Robustness, the ability to withstand failures and perturbations, is a critical attribute of many complex systems including complex networks.

<span class="mw-page-title-main">Bianconi–Barabási model</span>

The Bianconi–Barabási model is a model in network science that explains the growth of complex evolving networks. This model can explain that nodes with different characteristics acquire links at different rates. It predicts that a node's growth depends on its fitness and can calculate the degree distribution. The Bianconi–Barabási model is named after its inventors Ginestra Bianconi and Albert-László Barabási. This model is a variant of the Barabási–Albert model. The model can be mapped to a Bose gas and this mapping can predict a topological phase transition between a "rich-get-richer" phase and a "winner-takes-all" phase.

In statistical mechanics, Lee–Yang theory, sometimes also known as Yang–Lee theory, is a scientific theory which seeks to describe phase transitions in large physical systems in the thermodynamic limit based on the properties of small, finite-size systems. The theory revolves around the complex zeros of partition functions of finite-size systems and how these may reveal the existence of phase transitions in the thermodynamic limit.

<span class="mw-page-title-main">Amnon Aharony</span> Physicist at Ben Gurion University in Israel

Amnon Aharony is an Israeli Professor (Emeritus) of Physics in the School of Physics and Astronomy at Tel Aviv University, Israel and in the Physics Department of Ben Gurion University of the Negev, Israel. After years of research on statistical physics, his current research focuses on condensed matter theory, especially in mesoscopic physics and spintronics. He is a member of the Israel Academy of Sciences and Humanities, a Foreign Honorary Member of the American Academy of Arts and Sciences and of several other academies. He also received several prizes, including the Rothschild Prize in Physical Sciences, and the Gunnar Randers Research Prize, awarded every other year by the King of Norway.

Replica cluster move in condensed matter physics refers to a family of non-local cluster algorithms used to simulate spin glasses. It is an extension of the Swendsen-Wang algorithm in that it generates non-trivial spin clusters informed by the interaction states on two replicas instead of just one. It is different from the replica exchange method, as it performs a non-local update on a fraction of the sites between the two replicas at the same temperature, while parallel tempering directly exchanges all the spins between two replicas at different temperature. However, the two are often used alongside to achieve state-of-the-art efficiency in simulating spin-glass models.

References

  1. Domany, Eytan; Kinzel, Wolfgang (23 July 1984). "Equivalence of Cellular Automata to Ising Models and Directed Percolation". Physical Review Letters. 53 (4): 311.
  2. Hinrichsen, Haye (2000). "Non-equilibrium critical phenomena and phase transitions into absorbing states". Advances in Physics. 49 (7): 815–958. arXiv: cond-mat/0001070 . Bibcode:2000AdPhy..49..815H. doi:10.1080/00018730050198152. ISSN   0001-8732. S2CID   119106856.
  3. Hinrichsen, Haye; Jimenez-Dalmaroni, Andrea; Rozov, Yadin; Domany, Eytan (13 December 1999). "Flowing Sand: A Physical Realization of Directed Percolation". Physical Review Letters. 83 (24): 4999.
  4. Takeuchi, Kazumasa A.; Kuroda, Masafumi; Chaté, Hugues; Sano, Masaki (2007-12-05). "Directed Percolation Criticality in Turbulent Liquid Crystals". Physical Review Letters. 99 (23). American Physical Society (APS): 234503. arXiv: 0706.4151 . Bibcode:2007PhRvL..99w4503T. doi:10.1103/physrevlett.99.234503. ISSN   0031-9007. PMID   18233372. S2CID   6723341.
  5. Takeuchi, Kazumasa A.; Kuroda, Masafumi; Chaté, Hugues; Sano, Masaki (2009-11-16). "Experimental realization of directed percolation criticality in turbulent liquid crystals". Physical Review E. 80 (5). American Physical Society (APS): 051116. arXiv: 0907.4297 . Bibcode:2009PhRvE..80e1116T. doi:10.1103/physreve.80.051116. ISSN   1539-3755. PMID   20364956. S2CID   16205917.

Sources