Spotlight #57 – Хаос, теорія мереж та морфогенез

Більшість механічних систем є нелінійними: пружини жорсткішають при великій амплітуді, маятники сповільнюються поблизу вертикальної рівноваги, а металева балка під стисненням має два стійких положення спокою. Осцилятор Дуффінга охоплює всі ці явища в єдиному рівнянні з кубічним пружинним доданком:

ẍ + δẋ + αx + βx³ = γ cos(ωt)

Для β > 0 і α < 0 потенціальна енергія V(x) = ½αx² + ¼βx⁴ утворює подвійну яму: дві стійкі рівноваги в точках x = ±√(|α|/β), розділені нестійким піком при x = 0. Без збудження (γ=0) осцилятор просто осідає в тій ямі, з якої почав. За помірного збудження осцилятор може непередбачувано перестрибувати між ямами — це визначальна ознака хаосу.

Подвоєння періоду і шлях до хаосу

У міру зростання амплітуди збудження γ від нуля відгук системи проходить каскад подвоєння періоду: період-1 → період-2 → період-4 → період-8 → … → хаос, при послідовних значеннях параметра, що підпорядковуються сталій Фейгенбаума δ ≈ 4,669. Той самий каскад з'являється в будь-якому гладкому унімодальному відображенні і став одним із перших кількісних містків між різними хаотичними системами.

Як читати перерізи Пуанкаре

Переріз Пуанкаре фіксує стан фазового простору (x, ẋ) при заданій фазі збудження, зводячи неперервний 3D-потік до 2D-відображення.

• Одна точка — орбіта періоду 1: система точно повторюється кожен цикл збудження T = 2π/ω.
• n точок — орбіта періоду n: система повторюється через n циклів збудження (субгармонічний резонанс n:1).
• Фрактальна хмара — хаос: сам переріз Пуанкаре має фрактальну структуру із самоподібним шаруванням на всіх масштабах.
• Замкнена крива — квазіперіодичність: дві несумірні частоти, що часто трапляється поблизу меж стійкості.

У симуляції можна спостерігати еволюцію перерізу Пуанкаре в реальному часі. Спробуйте збільшити γ від 0,28 (дві точки) до 0,38 (чотири точки) і до 0,50 (хаотична хмара), щоб безпосередньо простежити каскад подвоєння періоду.

Досліджуйте в бібліотеці:

🌀 Осцилятор Дуффінга 🌿 Діаграма біфуркації 🎲 Гра в хаос ⚖️ Подвійний маятник

До 1990-х років стало зрозуміло, що багато реальних мереж — граф маршрутизаторів Інтернету, мережі білок-білкових взаємодій, графи спільного авторства, енергомережі — мають дві разючі спільні властивості: дуже коротку середню довжину шляху (кожен перебуває в межах «шести рукостискань») і дуже високу локальну кластеризацію (ваші друзі, як правило, знайомі між собою). Класичні випадкові графи мали короткі шляхи, але низьку кластеризацію; регулярні ґратки мали високу кластеризацію, але довгі шляхи. Альберт-Ласло Барабаші та Река Альберт (1999), а також Дункан Уоттс і Стівен Строгац (1998) кожен знайшов свою частину головоломки.

Ердеш–Реньї G(N,p) — випадковий базовий рівень

Найпростіша модель: кожне з N(N−1)/2 можливих ребер включається незалежно з імовірністю p. Розподіл ступенів є біноміальним (пуассонівським для великих N). Фазовий перехід перколяції відбувається при критичній імовірності p_c = 1/N: нижче неї всі компоненти малі (O(log N)); вище неї раптово з'являється єдина «гігантська компонента» розміру O(N). Модель Ердеша–Реньї передбачає пуассонівський розподіл ступенів — але реальні мережі мають значно більше вузлів-хабів високого ступеня, ніж це передбачає модель.

Барабаші–Альберт — преференційне приєднання і степеневі закони

Модель BA вирощує мережу, додаючи вузли по одному. Кожен новий вузол приєднує m ребер до наявних вузлів, але не випадково: імовірність з'єднання з вузлом i пропорційна його поточному ступеню k_i («багаті стають багатшими»). Цей механізм преференційного приєднання породжує безмасштабний розподіл ступенів P(k) ~ k^−3, з довгим степеневим хвостом: невелика кількість вузлів-хабів дуже високого ступеня (як-от Google чи великі аеропорти) співіснує з більшістю вузлів низького ступеня. Безмасштабні мережі стійкі до випадкових збоїв (випадкові відмови рідко зачіпають хаби), але вразливі до цілеспрямованих атак на хаби.

Уоттс–Строгац — режим малого світу

Починаючи з регулярної кільцевої ґратки — кожен вузол з'єднаний зі своїми k найближчими сусідами — кожне ребро перепідключається до випадкової цілі з імовірністю β.

• β = 0: регулярна ґратка — високий коефіцієнт кластеризації C ≈ 3(k−2)/4(k−1), але велика середня довжина шляху L ~ N/2k.
• β = 1: випадковий граф — низька кластеризація C ≈ k/N, коротка довжина шляху L ~ ln(N)/ln(k).
• β ≈ 0,01–0,1: оптимальна зона малого світу — майже регулярна кластеризація, але майже випадкова довжина шляху. Цей режим відповідає корі головного мозку людини, нейронному коннектому C. elegans, енергомережі Заходу США та топології доменів Інтернету.

Досліджуйте в бібліотеці:

🌐 Мережева наука 🐜 Мурашина колонія 🐦 Зграя птахів 🚌 Кластеризація автобусів

У 1952 році той самий Алан Тюрінг, який зламав код «Енігми» та заклав основи інформатики, опублікував працю, яку сам вважав найважливішою у своєму житті: «Хімічні основи морфогенезу». Центральне питання: як сферично-симетрична запліднена яйцеклітина перетворюється на асиметричний організм з окремими органами, віссю тіла та поверхневими візерунками? Відповідь Тюрінга була контрінтуїтивною — дифузія може спричиняти нестабільність.

Нестабільність Тюрінга

Розгляньмо дві хімічні речовини, що змішуються в однорідному середовищі. Зазвичай дифузія згладжує будь-який градієнт концентрації — це стабілізувальна сила. Тюрінг показав, що якщо одна речовина (активатор) каталізує власне виробництво, а інша (інгібітор) виробляється активатором і водночас дифундує швидше, однорідна рівновага стає нестійкою до малих просторових флуктуацій. Ці флуктуації зростають і насичуються до плямистих, періодичних візерунків — плям, смуг, гексагональних масивів — які повністю визначаються співвідношенням D_inhibitor / D_activator та кінетикою реакції.

Жодна генетична програма не задає «пляму в позиції (x,y)». Візерунок самоорганізується з випадкового шуму через фізичну нестабільність. Це один із найглибших прикладів спонтанного порушення симетрії в біології.

Динаміка Грея-Скотта та карта (F, k)

Модель Грея-Скотта є чистою реалізацією нестабільності Тюрінга:

∂U/∂t = Dᵤ∇²U  −  UV²  +  F(1 − U)
∂V/∂t = Dᵥ∇²V  +  UV²  −  (F + k)V

Площина параметрів (F, k) містить багатий атлас типів візерунків, експериментально картографований Пірсоном (Pearson, 1993). П'ять найяскравіших режимів доступні як готові пресети:

• Плями (F=0,055, k=0,062) — щільна гексагональна упаковка плям активатора. Аналог візерунків гепарда.
• Смуги (F=0,060, k=0,058) — паралельні смуги активатора. Аналог смуг зебри та смуг рифових риб.
• Лабіринт (F=0,029, k=0,057) — з'єднані лабіринтоподібні канали, подібні до візерунків на деяких тропічних мушлях.
• Корал (F=0,037, k=0,060) — деревоподібні структури із зубчастими краями, що нагадують коралові поліпи.
• Мітоз (F=0,028, k=0,053) — плями, що самовідтворюються, видовжуючись і розділяючись, — прямий аналог кінематики поділу клітин.

Експериментальні та біологічні докази

Механізм Тюрінга залишався суперечливим протягом 50 років, оскільки конкретні морфогени не були відомі. Докази накопичувалися починаючи з 1990-х років:

• Сигналізація Nodal–Lefty в пігментації рибок даніо (Kondo & Asai, 1995) — перше пряме експериментальне підтвердження.
• Сигналізація WNT–DKK контролює розташування волосяних фолікулів (Sick et al., 2006).
• Зворотний зв'язок FGF–SPROUTY визначає розташування пальців під час розвитку кінцівок (Miura et al., 2006).
• Візерунки складок на піднебінні мишей мають тюрінгівську природу (Economou et al., 2012).

Досліджуйте в бібліотеці:

🔬 Дифузія Тюрінга ⬛ Клітинні автомати 🌀 Реакція Білоусова-Жаботинського 🌿 Папороть Барнслі