Диффузионная неустойчивость · модель Шнакенберга

Узор Тьюринга

Узор Тьюринга — это рисунок из пятен или полос, возникающий в почти однородной смеси двух реагирующих веществ: крошечные возмущения усиливаются, потому что одно вещество распространяется намного быстрее другого.

Узоры Тьюринга в этой визуализации

На этой странице узоры Тьюринга возникают в безразмерной модели Шнакенберга с активатором uu и топливом vv:

∂u∂t=a−u+u2v+∇2u,∂v∂t=b−u2v+d ∇2v.\frac{\partial u}{\partial t} = a - u + u^2 v + \nabla^2 u, \qquad \frac{\partial v}{\partial t} = b - u^2 v + d\,\nabla^2 v .

Активатор и топливо поступают со скоростями aa и bb. Член u2vu^2 v отвечает реакции 2u+v→3u2u + v \to 3u, в которой активатор превращает топливо в дополнительный активатор; сам активатор также распадается. Топливо диффундирует в dd раз быстрее активатора. Состав фиксирован: a=0.05a = 0.05 и b=1.6b = 1.6, меняется только dd. Здесь можно увидеть возникновение узора при переходе порога, сравнить его с длинами волн, рост которых предсказывает теория, и добавить активатор кистью.

Насыщенность цвета показывает количество активатора или топлива, выбранного в разделе Показать. Диапазон подстраивается автоматически и указан в углу. Обе части малого графика используют одну ось длин волн. Вверху — теория: как быстро растёт небольшое возмущение каждой длины волны при данном dd. Внизу — измерение: сколько каждой длины волны сейчас содержится в изображении. Растущая полоса выделена фоном в обеих частях, поэтому можно наблюдать, как измеренные столбцы поднимаются внутри предсказанной теорией полосы.

Почему распространение может создавать узор? Реакция уравновешивается в однородной смеси u0=a+b=1.65u_0 = a + b = 1.65 и v0=b/(a+b)2≈0.588v_0 = b/(a+b)^2 \approx 0.588 и сама по себе там устойчива: в хорошо перемешанном сосуде малое возмущение затухает. Диффузия сама по себе тоже лишь выравнивает распределение.

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

Для малых возмущений около однородной смеси линейная теория описывает это точно. Волна cos⁡(k⋅x)\cos(\mathbf k \cdot \mathbf x) с длиной 2π/k2\pi/k сохраняет форму при диффузии, поэтому каждая длина волны растёт или затухает независимо. Топливо должно распространяться быстрее активатора, но одного этого мало: при данном составе скорость роста впервые достигает нуля при

dc≈17.66,2πkc≈10.03.d_c \approx 17.66, \qquad \frac{2\pi}{k_c} \approx 10.03 .

Ниже dcd_c все возмущения затухают. Выше растёт полоса длин волн около 10, расширяющаяся с увеличением dd. Сосуд допускает лишь определённые волны, поэтому чуть выше dcd_c в этой полосе может не оказаться ни одной допустимой волны. На сетке страницы критическая длина волны равна 9,86; малый график и индикатор показывают сеточные значения. Линейная теория говорит, какие длины волн растут первыми, а не каким будет конечный узор.

Линейная теория выбирает длину волны, а не форму. Компоненты нужной длины волны могут складываться в полосы, шестиугольную решётку пятен или ячеистую структуру отверстий. Нелинейные члены решают, что вырастет и сохранится. В этой модели три волны под углами 120° способны усиливать друг друга, что благоприятствует шестиугольникам вблизи порога.

В большом сосуде расстояния задаёт химия, а не сосуд: больший сосуд вмещает больше пятен, а не более крупные пятна. Форма важнее всего там, где ширина области примерно равна длине волны. Пример Туловище и хвост показывает пятна в широком туловище и поперечные полосы на узком хвосте.

Алан Тьюринг предложил механизм в 1952 году в работе «The Chemical Basis of Morphogenesis», а Юрген Шнакенберг ввёл эту модель в 1979 году.

На что обратить внимание

  • Преобладающая длина волны — не расстояние между пятнами. Шестиугольная решётка пятен образуется из трёх волн длины LL под углами 120°, а ближайшие пятна находятся на расстоянии 2L/3≈1.15 L2L/\sqrt3 \approx 1.15\,L.
  • Не для всякого узора животных доказано, что это узор Тьюринга. Экспериментальное подтверждение есть в химии и для нескольких случаев в биологии.

Связанные понятия