Реакционно-диффузионные системы

Цели урока

• Понять механизм Тьюринговой неустойчивости и условие D_v >> D_u
• Вывести дисперсионное соотношение и найти критическое волновое число
• Численно моделировать систему Gray-Scott с явной разностной схемой
• Связать математическую теорию с реальными биологическими паттернами

Предварительные знания

• PDE второго порядка
• Линейная алгебра (собственные значения)
• Устойчивость ОДУ

• Краевые задачи PDE

Почему диффузия - обычно стабилизирующий процесс - может дестабилизировать равновесие и породить структуру из хаоса?

• Паттерны на шерсти гепарда и полосы зебры - реакционно-диффузионная динамика
• Формирование пальцев в эмбриогенезе управляется системой BMP-Noggin
• Ориентационные колонки зрительной коры - нейронный аналог модели Тьюринга
• Neural Cellular Automata (Google Brain) обучают реакционно-диффузионные правила для self-repair

Алан Тьюринг и морфогенез

В 1952 году, преследуемый британскими властями, Тьюринг публикует 'The Chemical Basis of Morphogenesis' - статью, которую мало кто понял тогда и которая изменила биологию через 30 лет. Тьюринг предложил: химические вещества (морфогены) диффундируют и реагируют, создавая паттерны. Парадокс Тьюринга - диффузия дестабилизирует - оставался теоретическим до 1990 года, когда химики создали первую систему с контролируемой диффузией CIMA. В 2012 году генетики напрямую подтвердили тьюринговский механизм в нёбе мыши.

Неустойчивость Тьюринга и условие паттернообразования

В 1952 году Алан Тьюринг предложил математическую модель морфогенеза. Суть парадокса: однородное равновесие устойчиво без диффузии, но диффузия - обычно стабилизирующий процесс - его дестабилизирует. Именно это противоречие и порождает паттерны на коже рыб, шерсти леопардов, раковинах моллюсков.

В 2012 году эксперименты на мышах подтвердили: узор на нёбе определяется реакционно-диффузионной системой с FGF и Shh в роли активатора и ингибитора. Масштаб полос регулируется концентрацией FGF.

Длина волны паттерна L ~ 2*pi/k_c. Чем больше соотношение d = D_v/D_u, тем мельче паттерн. Это позволяет моделировать и крупные пятна гепарда, и мелкую сетку рыбы-зебры.

Модель Gray-Scott и численное моделирование

Модель Gray-Scott - классический пример химической реакции-диффузии с богатым фазовым портретом: при разных параметрах F (feed rate) и k (kill rate) возникают пятна, полосы, лабиринты, движущиеся пятна и хаотические режимы. Именно эта модель стала стандартом для визуализации реакционно-диффузионной динамики.

Явная схема Эйлера устойчива при dt <= dx^2 / (2*max(Du, Dv)). При нарушении этого условия решение взрывается за несколько шагов. Для больших сеток и долгих симуляций используют неявные или полунеявные схемы.

Паттерны в биологии и приложения в нейросетях

Реакционно-диффузионные системы описывают не только химию. Паттерны зрительной коры, ориентационные колонки V1 моделируются теми же уравнениями. Нейронные поля Вильсона-Коуэна - прямой аналог: возбуждение распространяется быстро (малая диффузия), торможение - медленнее (большая диффузия).

Нейронные клеточные автоматы (Neural Cellular Automata) от Google Brain (2020) обучают реакционно-диффузионные правила, чтобы клетки восстанавливали целевой паттерн после повреждения. Каждая клетка - крошечная нейросеть, читающая локальных соседей - прямая аналогия с химическим сигналингом.

Анализ устойчивости и фазовые портреты

Полный анализ реакционно-диффузионной системы включает: нахождение равновесий, линеаризацию, анализ спектра оператора диффузии-реакции, построение дисперсионной кривой. Критическое волновое число k_c задаёт масштаб доминирующего паттерна - первый нарастающий пространственный режим.

Реальные паттерны не строго синусоидальны - нелинейность насыщает нарастающие моды и стабилизирует амплитуду. Слабонелинейный анализ (уравнение Гинзбурга-Ландау) описывает этот переход от экспоненциального роста к устойчивому паттерну.

Связи с другими областями

Реакционно-диффузионные системы - пересечение нелинейной динамики, биологии развития и вычислительной математики.

• Метод конечных элементов — Связанная тема
• Бифуркационный анализ — Связанная тема
• Neural Cellular Automata — Связанная тема
• Уравнение Гинзбурга-Ландау — Связанная тема

Итоги

• Тьюрингова неустойчивость: равновесие устойчиво без диффузии, но неустойчиво к пространственным возмущениям при D_v >> D_u
• Критическое волновое число k_c задаёт масштаб паттерна; длина волны L_c = 2*pi/k_c
• Модель Gray-Scott воспроизводит пятна, полосы, лабиринты при разных параметрах F и k
• Биологические паттерны от нёба мыши до зрительной коры описываются теми же уравнениями

Вопросы для размышления

• Как изменится паттерн, если увеличить соотношение D_v/D_u в 10 раз?
• Почему нейронные поля Вильсона-Коуэна математически аналогичны системам Тьюринга?
• Какой биологический смысл у параметра F (feed rate) в модели Gray-Scott?

Связанные уроки

• de-25-fem — МКЭ решает реакционно-диффузионные системы на произвольных областях
• de-29-einstein-equations — Нелинейные PDE с пространственными паттернами - общий класс
• de-23-pde-bvp — Краевые задачи PDE - основа для реакционно-диффузионных систем