Про моделі випадкового росту

Ця симуляція порівнює дві класичні моделі випадкового росту кластерів: модель Едена та дифузійно-обмежену агрегацію (ДЛА). Обидві починаються з однієї зернової частинки в центрі сітки і ростуть назовні, але механізм росту кардинально відрізняється — утворюючи разюче різні геометричні структури.

Модель Едена додає частинки рівномірно випадковим чином у будь-яку порожню клітину, суміжну з існуючим кластером. Оскільки кожна комірка периметра має однакову ймовірність бути вибраною, виступи не мають переваги і кластер залишається компактним, наближаючись до кругового диска. Флуктуації поверхні належать до класу універсальності Кардара-Парізі-Чжана (KPZ).

ДЛА вводить дифузію: кожна нова частинка стартує далеко від кластера і здійснює випадкове блукання, поки не торкнеться кластера, де вона прилипає. Це створює екранування — кінчики та гілки перехоплюють блукачів до того, як ті можуть досягти фіордів (внутрішніх зазорів), змушуючи кластер розгалужуватися у самоподібний фрактал із виміром D≈1.71 у 2D.

Фрактальний вимір за методом підрахунку боксів та радіус інерції Rg обчислюються в реальному часі та відображаються на бічній панелі. Для ДЛА D має наближатися до ~1.7; для Едена — до ~2.0 (компактний диск).

Часті запитання

• Що таке дифузійно-обмежена агрегація (ДЛА)?
  Дифузійно-обмежена агрегація (ДЛА) — модель росту кластерів, запропонована Уіттеном і Сандером у 1981 році. У початковій точці фіксується зернова частинка. Потім частинки випускаються по одній із випадкового місця далеко від кластера і здійснюють випадкове блукання (броунівський рух), поки не торкнуться кластера, де вони прилипають назавжди. Утворений кластер є фрактальним деревом із характерною розгалуженою структурою та фрактальним виміром D≈1.71 у 2D. Кластери ДЛА схожі на сніжинки, мінеральні дендрити або шляхи блискавок.
• Що таке модель Едена?
  Модель Едена (запропонована Мюрреєм Еденом у 1961 році) вирощує кластер, випадково вибираючи одну клітину з периметра кластера (набір порожніх клітин, суміжних із кластером) і додаючи її. На відміну від ДЛА, тут немає дальнодіючої дифузії — ймовірність росту рівномірна по периметру незалежно від форми. Кластери Едена є компактними та приблизно круглими з фрактальною шорсткістю поверхні. Поверхня належить до класу універсальності Кардара-Парізі-Чжана (KPZ), що робить її актуальною для росту поверхні у тонких плівках, поширення бактеріальних колоній і росту пухлин.
• Чому ДЛА утворює фрактальні гілки, а Еден — компактний кластер?
  Ключова відмінність — екранування. У ДЛА випадкові блукачі здалеку набагато частіше вперше потрапляють на виступаючий кінчик кластера, ніж проникають у фіорд (внутрішній зазор). Це створює позитивний зворотний зв'язок: кінчики ростуть швидше, оскільки більше піддаються впливу блукачів, тоді як фіорди екрановані від них. Ця нестабільність екранування підсилює будь-який невеликий виступ у гілку, потім підгілки, утворюючи фрактальне дерево. Ріст Едена не має екранування: кожна комірка периметра має однакову ймовірність, тому виступи не ростуть переважно і кластер залишається компактним.
• Як фрактальний вимір вимірюється методом підрахунку боксів?

  Метод підрахунку боксів вимірює фрактальний вимір D, накладаючи сітки різного розміру клітин ε і підраховуючи N(ε) = кількість клітин, що містять хоча б одну частинку кластера. Для фракталу N(ε) ∝ ε^(−D), тому D = −lim(ε→0) log N(ε) / log ε. На практиці: вибирають розміри ε = 1,2,4,8,...,32; рахують заповнені бокси для кожного; апроксимують лінією лог-лог графік N vs 1/ε; нахил є D. Для ДЛА у 2D D ≈ 1.71. Для кластерів Едена (суцільний диск), D → 2.0. Для фрактальної кривої 1 < D < 2.
• Що таке радіус інерції кластера?

  Радіус інерції Rg = √(⟨r²⟩ − ⟨r⟩²), де r — відстань кожної частинки від центру мас кластера. Він вимірює просторовий розмір кластера. Для компактного кругового кластера з N частинок Rg ∝ √N (2D диск: Rg = R/√2). Для фракталу ДЛА з виміром D, Rg ∝ N^(1/D). Вимірювання Rg як функції від N дає ще один спосіб оцінити D: нахил log Rg vs log N у лог-лог графіку дорівнює 1/D.
• Які реальні явища моделює ДЛА?

  ДЛА та подібні моделі широко зустрічаються в природі. Мінеральні дендрити (відкладення оксиду марганцю у вигляді папоротей у каменях) формуються шляхом дифузії іонів, подібної до ДЛА. Електроосадження цинку виробляє фрактальні дерева цинку в розчині. Блискавки створюють розгалужені канали розряду, що слідують шляху найменшого опору (подібно до ДЛА). Пальцеподібні явища в'язкості (один флюїд витісняє інший у пористій породі) виробляють пальці Саффмана-Тейлора, подібні до ДЛА у 2D. Розгалуження рук сніжинок, ріст аксонів нейронів і деякі схеми інвазії пухлин демонструють фрактальну геометрію класу ДЛА.
• Чому фрактальний вимір ДЛА приблизно дорівнює 1.71?

  Точний фрактальний вимір ДЛА в 2D чисельно обчислений як D ≈ 1.710 ± 0.004, підтверджений симуляціями мільярдів частинок. Точне теоретичне передбачення залишається одним із відкритих питань статистичної фізики. Евристично D визначається балансом між дифузією (яка згладжує, прагнучи до D=2) та екрануванням (яке створює гілки, зменшуючи D нижче 2). У 3D вимір ДЛА D ≈ 2.50. У дуже великих вимірах d, D → d−1, коли кластер стає гіпердеревом.
• Як ефективно симулюється випадкове блукання в ДЛА?

  Наївне випадкове блукання, де частинка рухається по одному кроку за ітерацію, є дуже повільним: частинка, запущена з радіуса R=100, в середньому робить ~R² ≈ 10,000 кроків лише для досягнення кластера. Тут використовуються дві оптимізації. По-перше, 'коло смерті': якщо блукач виходить за максимальний радіус (радіус кластера + запас безпеки), його знищують і перезапускають. По-друге, 'позарешітковий стрибок': коли блукач далеко від будь-якої частинки кластера, він може стрибнути на більший випадковий крок без ризику перескочити через кластер, значно зменшуючи кількість кроків.
• Що таке клас Кардара-Парізі-Чжана (KPZ) і чому Еден належить до нього?

  Клас універсальності Кардара-Парізі-Чжана (KPZ) описує стохастичний ріст поверхні: шорстка межа розділу кластера Едена належить до цього класу. Ключові ознаки: ширина межі w росте як t^β із KPZ показником β=1/3 у 2D, а кореляції висота-висота спадають із показником χ=1/2. Це означає, що поверхні Едена, випадкове осадження з релаксацією, краї бактеріальних колоній, плями від кави та навіть 1D транспорт TASEP — усі вони належать до одного класу універсальності — глибокий зв'язок у нерівноважній статистичній механіці.
• Як кількість частинок впливає на фрактальний візерунок?

  У ДЛА фрактальна структура стає чіткішою зі збільшенням кількості частинок. Малі кластери (N<100) виглядають як неправильні павутиноподібні форми. При N≈1000 ієрархія розгалуження стає видимою. При N≈10,000 самоподібна структура чітко простежується на 2–3 декадах масштабу. Фрактальний вимір D можна надійно оцінити лише для N>500 частинок. У кластерах Едена форма наближається до кола зі збільшенням N (за центральною граничною теоремою), але шорсткість поверхні зберігається — вона ніколи не стає ідеально гладкою, оскільки нові флуктуації завжди додаються на периметрі.