🌸 Wzrost losowy — model Edena i fraktal DLA

O modelach losowego wzrostu

Ta symulacja porównuje dwa klasyczne modele losowego wzrostu klastrów: model Edena oraz dyfuzyjnie ograniczoną agregację (DLA). Oba zaczynają się od pojedynczej cząstki zalążkowej w centrum siatki i rosną na zewnątrz, ale mechanizm wzrostu różni się dramatycznie — dając uderzająco odmienne struktury geometryczne.

Model Edena dodaje cząstki równomiernie losowo do dowolnej pustej komórki sąsiadującej z istniejącym klastrem. Ponieważ każde miejsce na obwodzie jest równie prawdopodobne do wyboru, wypustki nie mają przewagi, a klaster pozostaje zwarty, zbliżając się do okrągłego dysku w miarę wzrostu. Fluktuacje powierzchni należą do klasy uniwersalności Kardara-Parisiego-Zhanga (KPZ).

DLA wprowadza dyfuzję: każda nowa cząstka zaczyna daleko od klastra i wykonuje błądzenie losowe, aż dotknie klastra, po czym się przykleja. Tworzy to ekranowanie — końcówki i gałęzie przechwytują nadchodzące cząstki błądzące, zanim zdążą one dotrzeć do fiordów (wewnętrznych luk), powodując rozgałęzianie się klastra w samopodobny fraktal o wymiarze D≈1,71 w 2D.

Wymiar fraktalny liczony metodą zliczania pudełek (box-counting) oraz promień żyracji Rg są obliczane w czasie rzeczywistym i wyświetlane na panelu bocznym. Dla DLA, D powinno dążyć do ~1,7; dla Edena — do ~2,0 (zwarty dysk).

Najczęściej zadawane pytania

• Czym jest dyfuzyjnie ograniczona agregacja (DLA)?
  Dyfuzyjnie ograniczona agregacja (DLA) to model wzrostu klastra wprowadzony przez Wittena i Sandera w 1981 roku. Cząstka zalążkowa jest ustalona w początku układu. Następnie cząstki są wypuszczane jedna po drugiej z losowej pozycji daleko od klastra i wykonują błądzenie losowe (ruch Browna), dopóki nie dotkną klastra, gdzie przyklejają się na stałe. Powstały klaster to fraktalne drzewo o charakterystycznej strukturze rozgałęzień i wymiarze fraktalnym D≈1,71 w 2D. Klastry DLA wyglądają jak płatki śniegu, dendryty mineralne lub ścieżki wyładowań piorunowych.
• Czym jest model Edena?
  Model Edena (zaproponowany przez Murraya Edena w 1961 roku) rozrasta klaster, losowo wybierając jedną komórkę z obwodu klastra (zbioru pustych komórek sąsiadujących z klastrem) i dodając ją. W przeciwieństwie do DLA, nie ma tu dyfuzji dalekiego zasięgu — prawdopodobieństwo wzrostu jest jednolite na całym obwodzie, niezależnie od kształtu. Klastry Edena są zwarte i mniej więcej okrągłe, z fraktalną chropowatością powierzchni. Powierzchnia należy do klasy uniwersalności Kardara-Parisiego-Zhanga (KPZ), co czyni ją istotną dla wzrostu powierzchni w cienkich warstwach, rozprzestrzeniania się kolonii bakteryjnych i wzrostu guzów.
• Dlaczego DLA daje fraktalne rozgałęzienia, a Eden daje zwarty klaster?
  Kluczową różnicą jest ekranowanie. W DLA cząstki błądzące losowo z daleka znacznie częściej trafiają najpierw na wystający czubek klastra, niż przenikają do fiordu (wewnętrznej luki). Tworzy to dodatnie sprzężenie zwrotne: czubki rosną szybciej, ponieważ są bardziej wyeksponowane na nadchodzące cząstki błądzące, podczas gdy fiordy są ekranowane przed cząstkami. Ta niestabilność ekranowania wzmacnia każdą małą wypustkę w gałąź, a następnie w podgałęzie, tworząc fraktalne drzewo. Wzrost w modelu Edena nie ma ekranowania: każde miejsce na obwodzie ma równe prawdopodobieństwo, więc wypustki nie rosną preferencyjnie, a klaster pozostaje zwarty.
• Jak mierzy się wymiar fraktalny metodą zliczania pudełek?

  Zliczanie pudełek (box counting) mierzy wymiar fraktalny D, nakładając siatki o różnych rozmiarach komórek ε i licząc N(ε) = ile komórek zawiera co najmniej jedną cząstkę klastra. Dla fraktala, N(ε) ∝ ε^(−D), więc D = −lim(ε→0) log N(ε) / log ε. W praktyce: wybiera się rozmiary ε = 1,2,4,8,...,64 piksele; liczy się zajęte pudełka dla każdego z nich; dopasowuje się linię do wykresu log-log N w funkcji 1/ε; nachylenie to D. Dla DLA w 2D, D ≈ 1,71. Dla klastrów Edena (lity dysk), D → 2,0. Dla krzywej fraktalnej, 1 < D < 2.
• Czym jest promień żyracji klastra?

  Promień żyracji Rg = √(⟨r²⟩ − ⟨r⟩²), gdzie r to odległość każdej cząstki od środka masy klastra. Mierzy on rozciągłość przestrzenną klastra. Dla zwartego, okrągłego klastra złożonego z N cząstek, Rg ∝ √N (dysk 2D: Rg = R/√2). Dla fraktala DLA o wymiarze D, Rg ∝ N^(1/D). Pomiar Rg jako funkcji N daje inny sposób oszacowania D: nachylenie log Rg w funkcji log N na wykresie log-log jest równe 1/D.
• Jakie rzeczywiste zjawiska modeluje DLA?

  DLA i wzrost typu DLA występuje szeroko w naturze. Dendryty mineralne (osady tlenku manganu przypominające paprocie w skale) powstają wskutek dyfuzji jonów podobnej do DLA. Elektroosadzanie cynku daje fraktalne drzewa cynkowe w roztworze. Pioruny tworzą rozgałęzione kanały wyładowań podążające ścieżką najmniejszego oporu (podobnie jak DLA). Palcowanie lepkościowe (viscous fingering) (jeden płyn wypierający inny w porowatej skale) daje palce Saffmana-Taylora, które są podobne do DLA w geometrii 2D. Rozgałęzienia ramion płatków śniegu, wzrost aksonów neuronów oraz niektóre wzorce inwazji nowotworowej wykazują geometrię fraktalną klasy DLA.
• Dlaczego wymiar fraktalny DLA wynosi w przybliżeniu 1,71?

  Dokładny wymiar fraktalny DLA w 2D obliczono numerycznie na D ≈ 1,710 ± 0,004, co potwierdzają symulacje miliardów cząstek. Dokładna teoretyczna predykcja pozostaje jednym z otwartych problemów fizyki statystycznej. Heurystycznie, D jest determinowane przez równowagę między dyfuzją (która wygładza, dążąc do D=2) a ekranowaniem (które tworzy rozgałęzienia, obniżając D poniżej 2). W 3D wymiar DLA wynosi D ≈ 2,50. W bardzo wysokich wymiarach d, D → d−1, gdy klaster staje się hiperdrzewem.
• Jak efektywnie symulowane jest błądzenie losowe w DLA?

  Naiwne błądzenie losowe, w którym cząstka wykonuje jeden krok na iterację, jest bardzo powolne: cząstka wystrzelona z promienia R=100 potrzebuje średnio ~R² ≈ 10 000 kroków tylko po to, by dotrzeć do klastra. Zastosowano tu dwie optymalizacje. Pierwsza to "krąg śmierci": jeśli cząstka błądząca przekroczy maksymalny promień (promień klastra + margines bezpieczeństwa), zostaje usunięta i uruchomiona od nowa. Druga to "skoki poza siatką": gdy cząstka błądząca jest daleko od jakiejkolwiek cząstki klastra (najbliższa cząstka klastra w odległości d), może wykonać losowy skok o odległość do d−1 bez ryzyka przeskoczenia klastra, co drastycznie zmniejsza liczbę kroków.
• Czym jest klasa Kardara-Parisiego-Zhanga (KPZ) i dlaczego Eden do niej należy?

  Klasa uniwersalności Kardara-Parisiego-Zhanga (KPZ) opisuje stochastyczny wzrost powierzchni: chropowaty interfejs klastra Edena należy do tej klasy. Kluczowe sygnatury to: szerokość interfejsu w rośnie jak t^β z wykładnikiem KPZ β=1/3 w 2D, a korelacje wysokość-wysokość zanikają z wykładnikiem χ=1/2. Oznacza to, że powierzchnie Edena, losowe osadzanie z relaksacją, krawędzie kolonii bakteryjnych, plamy po kawie, a nawet jednowymiarowy transport cząstek TASEP należą do tej samej klasy uniwersalności — głębokiego związku w nierównowagowej mechanice statystycznej.
• Jak liczba cząstek wpływa na wzór fraktalny?

  W DLA struktura fraktalna staje się wyraźniejsza wraz ze wzrostem liczby cząstek. Małe klastry (N<100) wyglądają jak nieregularne, pajęcze kształty. Przy N≈1000 widoczna staje się hierarchia rozgałęzień. Przy N≈10 000 samopodobna struktura jest wyraźna w zakresie 2–3 rzędów wielkości skali. Wymiar fraktalny D można wiarygodnie oszacować tylko dla N>500 cząstek. W klastrach Edena kształt klastra zbliża się do koła wraz ze wzrostem N (zgodnie z centralnym twierdzeniem granicznym), ale chropowatość powierzchni się utrzymuje — nigdy nie staje się idealnie gładka, ponieważ na obwodzie zawsze pojawiają się nowe fluktuacje.