🌸 Випадковий Ріст — Модель Едена та Фрактал ДЛА

Дифузійно-обмежена агрегація (ДЛА) — це модель росту кластера, введена Віттеном і Сандером у 1981 році. Зародкова частинка фіксується в початку координат. Потім частинки випускаються по одній з випадкової позиції далеко й здійснюють випадкове блукання (броунівський рух), доки не торкнуться кластера, де назавжди приклеюються. Отриманий кластер — це фрактальне дерево з характерною розгалуженою структурою й фрактальною розмірністю D≈1.71 у 2D. Кластери ДЛА виглядають як сніжинки, мінеральні дендрити або шляхи блискавкового розряду.

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

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

Підрахунок коробок вимірює фрактальну розмірність D, накладаючи сітки з різними розмірами комірок ε і рахуючи N(ε) = скільки комірок містять принаймні одну частинку кластера. Для фракталу N(ε) ∝ ε^(−D), тож D = −lim(ε→0) log N(ε) / log ε. На практиці: оберіть розміри ε = 1,2,4,8,...,64 пікселів; порахуйте заповнені комірки для кожного; підберіть пряму до логарифмічного графіка N проти 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 проти log N на логарифмічному графіку дорівнює 1/D.

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

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

Наївне випадкове блукання, де частинка рухається на один крок за ітерацію, дуже повільне: частинка, запущена з радіуса R=100, робить у середньому ~R² ≈ 10 000 кроків, щоб просто досягти кластера. Тут використовуються дві оптимізації. По-перше, «коло смерті»: якщо блукач виходить за межі максимального радіуса (радіус кластера + запас безпеки), його вбивають і перезапускають. По-друге, «стрибок поза ґраткою»: коли блукач далеко від будь-якої частинки кластера (найближча частинка кластера на відстані d), він може стрибнути на випадкову відстань до d−1 без ризику проскочити кластер, значно зменшуючи кількість кроків.

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

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