Die Perkolationstheorie (lat. percolare – durchsickern) ist ein Teilgebiet der Stochastik und der statistischen Physik und beschreibt das Ausbilden zusammenhängender Gebiete (Cluster) bei zufallsbedingtem Besetzen von Strukturen (Gittern).

Die Perkolationstheorie eignet sich zur Modellierung physikalischer Phänomene, die bei Phasenübergängen entstehen. Des Weiteren eignet sie sich zur Beschreibung der elektrischen Leitfähigkeit von Legierungen, der Ausbreitung von Epidemien und Waldbränden, sowie zur Beschreibung von Wachstumsmodellen. In der Geologie und Hydrologie beschreibt die Perkolation einfache Modelle zur Ausbreitung von Flüssigkeiten in porösem Gestein (siehe Perkolation (Technik)), die als anschauliche Beispiele der unten beschriebenen Clusterbildung dienen. Unterarten sind die Punktperkolation, bei der Gitterpunkte mit einer bestimmten Wahrscheinlichkeit besetzt werden, und die Kantenperkolation, bei der besetzte Punkte untereinander verbunden werden. Man kann sich beliebige zufällig erzeugte Objekte, z. B. Tröpfchen, vorstellen, die untersucht werden.

Definition

Bearbeiten

Eine Perkolation ist ein zufälliger Teilgraph eines Graphen  , genannt Referenzgraph, wobei   eine abzählbare Menge von Knoten bezeichnet und   die abzählbare Menge der Kanten. Man unterscheidet zwischen einer Verbindungs-Perkolation (englisch bond percolation) und einer Orts-Perkolation (englisch site percolation).

Die Verbindungs-Perkolation entsteht durch eine Wahrscheinlichkeitsfunktion  , die jeder Kante   entweder den Zustand

  • offen (auch besetzt genannt)   oder
  • geschlossen (auch leer genannt)  

zuweist. Die Verbindungs-Perkolation wird daher auch als Kanten-Perkolation bezeichnet.

Die Orts-Perkolation entsteht durch eine Wahrscheinlichkeitsfunktion  , die jedem Knoten   entweder den Zustand

  • offen / besetzt   oder
  • geschlossen / leer  

zuweist. Die Orts-Perkolation wird daher auch als Knoten-Perkolation bezeichnet.[1]

Die Perkolations-Konfiguration ist

 

Der resultierende Zufallsgraph bildet eine Sammlung von Clustern.

Bernoulli-Perkolation

Bearbeiten

Das einfachste Modell ist die Bernoulli-Perkolation von Broadbent und Hammersley (1957[2]). Die Wahrscheinlichkeitsfunktion folgt dabei einer Bernoulli-Verteilung, das heißt, für   (respektive  ) gilt,   tritt mit Wahrscheinlichkeit   und   mit Wahrscheinlichkeit   ein.

Als Referenzgraphen wählt man dabei ein Gitter wie zum Beispiel  , wobei  .[1]

Wir notieren mit   einen zusammenhängenden Pfad der Länge  , der im Ursprung startet. Dann notiert   die Wahrscheinlichkeit

 .

Die unendliche Verbindung wird oft als   notiert.[1]

Kritische Wahrscheinlichkeit

Bearbeiten

Für Bernoulli-Perkolationen auf transitiven unendlichen Graphen existiert eine kritische Wahrscheinlichkeit  , so dass die Wahrscheinlichkeit, dass ein unendlicher zusammenhängender Teilgraph existiert,   ist, wenn  , und  , wenn  .[1]

Die Region   bezeichnet man als subkritische Region und   als superkritische Region.

Erläuterungen

Bearbeiten

Es werden auch Kombinationen von beiden Arten der Perkolation untersucht. Seit 1957 spielen Computersimulationen eine entscheidende Rolle in der Perkolationstheorie. Obwohl viele interessante Größen heute genau bekannt sind, lieferten Simulationen in vielen Fällen die notwendige Intuition, um sie zu bestimmen, und die meisten interessanten Größen sind immer noch nur durch Simulationen insgesamt zugänglich. Perkolationsmodelle können in einer Vielzahl verschiedener Systeme verwendet werden.[3]

Geschichte

Bearbeiten

Historisch geht die Perkolationstheorie (englisch percolation theory) auf Paul Flory und Walter H. Stockmayer zurück, welche die Flory-Stockmayer-Theorie in den 1940er Jahren entwickelten, um Polymerisationsprozesse bei der Gelierung zu beschreiben. Der Polymerisationsprozess kommt durch das Aneinanderreihen von Molekülen zustande, die dadurch Makromoleküle bilden. Der Verbund solcher Makromoleküle führt zu einem Netzwerk von Verbindungen, die sich durch das ganze System ziehen können. Broadbent und Hammersley führten das moderne Konzept der Perkolation ein[4]. Ein anderes Modell zur Beschreibung von Zufallsgraphen ist das Erdős–Rényi-Modell.

Modellbildung

Bearbeiten

Perkolationen werden auf Gittern modelliert, wobei Kristallgitter Interpretationen mathematischer Gitter sind.

Knotenperkolation (site percolation)

Bearbeiten
 
Knoten- und Kanten-Perkolation

Allgemein lässt sich ein einfaches Modell für die „Knoten-“ oder „Platzperkolation“ konstruieren:

Die Felder eines zweidimensionalen Quadratgitters werden mit einer bestimmten Wahrscheinlichkeit besetzt. Ob ein Feld besetzt wird oder leer bleibt, ist unabhängig von der Besetzung aller anderen Felder. Des Weiteren wird das Gitter als so groß angenommen, dass Randeffekte vernachlässigt werden können, und ist im Idealfall unendlich groß. Abhängig von der gegebenen Verteilung werden sich Gruppen auf dem Gitter bilden, d. h. besetzte Felder in unmittelbarer Nachbarschaft. Diese Gruppen – als Cluster bezeichnet – werden umso größer sein, je größer die Wahrscheinlichkeit zur Besetzung eines Feldes ist. Die Perkolationstheorie beschäftigt sich nun mit Eigenschaften wie Größe oder Anzahl dieser Cluster.

Wenn   die Wahrscheinlichkeit ist, dass ein Feld besetzt ist, so bilden sich mit dem Ansteigen von   größere Cluster aus. Die Perkolationsschwelle   ist definiert als der Wert von  , bei dem mindestens ein Cluster eine Größe erreicht, dass er sich durch das gesamte System erstreckt, also eine Ausdehnung auf dem Gitter von der rechten zur linken und von der oberen zur unteren Seite hat. Man sagt: „Der Cluster perkoliert durch das System.“

Kantenperkolation (bond percolation)

Bearbeiten

Das Gegenstück dazu wird „Kantenperkolation“ (englisch bond percolation) genannt.

 
Kantenperkolation in zwei Dimensionen mit Kanten-Besetzungs­wahrscheinlichkeit p = 0,51 auf einem 50×50-Ausschnitt. Es gibt einen Pfad geschlossener Kanten, der die untere mit der oberen Bildkante verbindet.
 
Kantenperkolation in einem Quadratgitter mit p = 0,3 bis p = 0,52

Ein Gitter, z. B. oben genanntes Quadratgitter, ist vollständig besetzt und von jedem Feld des Gitters bestehen vier Verbindungen zu den jeweils vier Nachbarfeldern. Nun ist mit einer Wahrscheinlichkeit   eine Verbindung zu einem Nachbarfeld geöffnet und mit einer Wahrscheinlichkeit   die Verbindung geschlossen. Diese Art der Perkolation lässt sich gut mit dem oben genannten Modell in der Geologie vergleichen: Die Hohlräume in einem porösen Gestein sind mit Wasser gefüllt und durch ein Netzwerk von Kanälen verbunden; mit einer Wahrscheinlichkeit   besteht ein Kanal zwischen zwei nächsten Nachbarn und mit einer Wahrscheinlichkeit von   besteht keiner.

Ein Cluster ist dann definiert als Gruppe von Gitterplätzen, die durch offene Kanäle verbunden sind. Auch hier ist   wieder die Perkolationsschwelle, und für   gibt es einen Cluster, der durch das gesamte System perkoliert, während ein solcher Cluster bei   nicht existiert. Die Perkolationsschwelle ist bei der Kantenperkolation niedriger als bei Systemen, welche sich entsprechend der Knotenperkolation verhalten. Das gilt für alle Gittertypen.

Gittertyp[5][6] Knotenperkolationsschwelle Kantenperkolationsschwelle
Sechseckgitter 0,6962 0,6527… = 1 − 2 sin (π/18)[7]
Quadratgitter 0,592746 0,5
Dreiecksgitter 0,5 0,34729… = 2 sin (π/18)[7]
Diamantgitter 0,43 0,388
einfach kubisches Gitter 0,3116 0,2488
BCC 1. 0,246 0,1803
FCC 0,198 0,119
Hyperkubisches Gitter (4d) 0,197 0,1601
Hyperkubisches Gitter (5d) 0,141 0,1182
Hyperkubisches Gitter (6d) 0,107 0,0942
Hyperkubisches Gitter (7d) 0,089 0,0787

Ein Cluster der Größe   entsteht, wenn benachbarte   Plätze nebeneinander besetzt sind, die von zwei leeren Plätzen begrenzt werden. Wenn die Größe des Clusters gegen unendlich geht, können die Effekte der Ränder ignoriert werden. Die Wahrscheinlichkeit, dass ein Platz auf der linken Seite des Clusters besetzt ist, ist dann  .

Für   ist die Wahrscheinlichkeit, dass ein beliebiger Platz zu einem endlichen Cluster gehört, gleich der Wahrscheinlichkeit  , dass der Platz besetzt ist. Weil die Wahrscheinlichkeit, dass eine beliebiger Platz zu einem  -Cluster gehört, gegeben ist durch  , erhält man für   mithilfe der geometrische Reihe

 ,

also  . Für die mittlere Größe   des Clusters ergibt sich  .

Die Stärke   des unendlichen Clusters ist die Wahrscheinlichkeit, dass ein beliebiger Platz zum Cluster gehört. Die Stärke misst, wie groß der Anteil der Plätze des unendlichen Gitters ist, die zum unendlichen Cluster gehören, und wird Ordnungsparameter (englisch order parameter) genannt. Das Phänomen, dass der Ordnungsparameter für   größer als 0 wird, ist als Phasenübergang bekannt und   wird als kritische Besetzungswahrscheinlichkeit (englisch critical occupation probability) bezeichnet.

Ein besetzter Platz gehört entweder zum unendlichen Cluster oder zu einem endlichen Cluster. Daher gilt

 

für alle  . Daraus folgt

 

Der Term   ist höchstens von der Größenordnung   (siehe Landau-Symbole). Die Summe läuft über alle endlichen Clustergrößen   und schließt das unendliche Cluster aus. Für   ist  .[8]

Gerichtete Perkolation

Bearbeiten

Die gerichtete Perkolation (englisch directed percolation) lässt sich anschaulich mit einer Kaffeemaschine (englisch coffee percolator) oder porösem Gestein erklären.

Anhand der Kantenperkolation wird der Unterschied zwischen „normaler“ bzw. isotroper Perkolation und der gerichteten Perkolation klar.

 
Gerichtete Perkolation

Wenn Wasser auf ein poröses Medium gegossen wird, stellt sich die Frage, ob das Medium durchdrungen werden kann, d. h., ob es einen Kanal von der Oberseite zur Unterseite des Mediums gibt, oder ob das Wasser vom Medium absorbiert wird. Die Wahrscheinlichkeit, dass das Wasser auf einen offenen Kanal trifft, ist wie bei einer isotropen Perkolation gegeben durch  . Im Gegensatz zur isotropen Perkolation existiert jedoch eine gegebene Vorzugsrichtung: Wasser in porösem Gestein und ebenso in der Kaffeemaschine bewegt sich in die Richtung, die durch die Gravitation bestimmt wird. Die Perkolationsschwelle   ist bei der gerichteten Perkolation größer als bei der isotropen Perkolation.

Mikrokanonisches und kanonisches Ensemble

Bearbeiten

Um Finite-Size-Effekte und den Phasenübergang zu studieren, müssen die betreffenden Observables für einen großen Bereich der Wahrscheinlichkeit   bestimmt werden, wenn nicht für eine kontinuierliche Nachbarschaft der Mindestwahrscheinlichkeit   des gesamten Einheitsintervalls. Diese Aufgabe kann erheblich vereinfacht werden, indem man vom kanonischen zum mikrokanonischen Ensemble wechselt.

Wenn jede der   Kanten des Gitters mit Wahrscheinlichkeit   besetzt ist und mit Wahrscheinlichkeit   leer ist, spricht man auch von einem kanonischen Ensemble. Wenn genau   Kanten des Gitters besetzt sind, spricht man auch von einem mikrokanonischen Ensemble. Wenn   ein Observable im kanonischen Ensemble ist, dann wird das entsprechende Observable im mikrokanonischen Ensemble mit   bezeichnet. Das gewünschte Observable kann durch eine Faltung mit der Binomialverteilung erhalten werden. Diese kann wie folgt geschrieben werden:

 

Unter Verwendung des mikrokanonischen Ensembles ist es somit möglich,   für alle   auf einmal zu bestimmen und dies, indem nur   Werte gemessen werden.[3]

Anwendung in der Epidemiologie

Bearbeiten

Für die Anwendung in der Epidemiologie kann man sich für den Erreger empfängliche Individuen als Knoten und Kontakt zwischen den Individuen als Kanten vorstellen. Ab einer bestimmten Populationsdichte, wenn es genügend Kontakte zwischen den Individuen gibt, würde die Perkolationsschwelle überschritten, sich also große, zusammenhängende Cluster bilden, die zu einer Ausbreitung des Erregers auf größere Bereiche der Population führen. Empirisch wurde die Existenz einer solchen Perkolationsschwelle anhand der Großen Rennmaus gezeigt, deren Kolonien unterschiedliche Populationsdichten aufweisen.[9]

Anwendungen im Alltag

Bearbeiten

Im täglichen Leben kommen viele perkolationsartige Phasenübergänge vor, z. B. das „Puddingproblem“ (Gel-Bildung[10]), das „Sahnesteif-Problem“ und das Problem der „Verklumpung“. In allen Fällen geht die Wirkung erst bei Überschreiten eines kritischen Wertes des ursächlichen Parameters gegen das erwünschte oder unerwünschte Maximum, und zwar meist nach einem Potenzgesetz mit einem kritischen Exponenten, wobei die maximale Wirkung bei Überschreiten des kritischen Wertes zunächst sehr rasch ansteigt. Durch chemische Zusätze, etwa Pudding- oder „Sahnesteif“-Pulver, kann man den kritischen Wert herabsetzen, ohne allerdings das Prinzip zu ändern.

Siehe auch

Bearbeiten

Literatur

Bearbeiten
  • P. J. Flory: Thermodynamics of High Polymer Solutions. In: Journal of Chemical Physics. 9, Nr. 8, August 1941, S. 660.
  • P. J. Flory: Thermodynamics of high polymer solutions. In: J. Chem. Phys. 10, 1942, S. 51–61.
  • W. H. Stockmayer: Theory of molecular size distribution and gel formation in branched polymers. In: J. Chem. Phys. 11, 1943, S. 45–55.
  • D. Stauffer, A. Aharony: Introduction to Percolation Theory. Taylor and Fransis, London 1994.
  • D. Achlioptas u. a.: Explosive Percolation in Random Networks. In: Science. 2009.
  • A. Bunde, H. E. Roman: Gesetzmäßigkeiten der Unordnung. In: Physik in unserer Zeit. 27, 1996, S. 246–256.
  • Vincent Beffara, Vladas Sidoravicius: Percolation. In: Encyclopedia of Mathematical Physics. Elsevier, 2006. Arxiv

Einzelnachweise

Bearbeiten
  1. a b c d Hugo Duminil-Copin: Sixty years of percolation. Hrsg.: arXiv. 2017, doi:10.48550/ARXIV.1712.04651, arxiv:1712.04651 [abs].
  2. Simon Broadbent und John Hammersley: Percolation processes I. Crystals and mazes. In: Mathematical Proceedings of the Cambridge Philosophical Society. Band 53, Nr. 3, 1957, ISSN 0305-0041, S. 629–641, doi:10.1017/S0305004100032680.
  3. a b Universität Bonn, Jan Hasenbusch, Matthias Wilhelm, : Project in Computational Physics, Percolation
  4. S. R. Broadbent, J. M. Hammersley: Percolation processes. In: Mathematical Proceedings of the Cambridge Philosophical Society. 53. Jahrgang, Nr. 3, 2008, ISSN 0305-0041, S. 629, doi:10.1017/S0305004100032680, bibcode:1957PCPS...53..629B.
  5. Werte entnommen aus uni-stuttgart.de, Skript (27. April 2005): Simulationsmethoden (Memento vom 21. September 2004 im Internet Archive) (PDF; 1,3 MB), S. 41–52.
  6. Humboldt-Universität zu Berlin: Notes on percolation theory
  7. a b M. F. Sykes, J. W. Essam: Exact critical percolation probabilities for site and bond problems in two dimensions. In: Journal of Mathematical Physics. Band 5, Nr. 8, 1964, S. 1117–1127, doi:10.1063/1.1704215, bibcode:1964JMP.....5.1117S (englisch, bme.hu [PDF]).
  8. Dr. Kim Christensen, Imperial College London: Percolation Theory
  9. S. Davis, P. Trapman, H. Leirs, M. Begon, J. a. P. Heesterbeek: The abundance threshold for plague as a critical percolation phenomenon. In: Nature. Band 454, Nr. 7204, Juli 2008, ISSN 1476-4687, S. 634–637, doi:10.1038/nature07053 (nature.com [abgerufen am 10. September 2020]).
  10. Siehe etwa die Dissertation von Markus Lechtenfeld im Fachbereich Chemie der Universität Duisburg, zum Thema Auswertung der rheologischen und optischen Untersuchungen während der Gelierung des Systems Gelatine/Wasser mit Hilfe der Perkolationstheorie. Duisburg 2001, (online)