Transition disorder–order–disorder in reaction-diffusion biophysical system

Full Text

Введение. В последние годы наблюдается все возрастающий интерес к изучению влияния внешних шумов на открытые нелинейные пространственно распределенные системы, что подтверждается появлением в ведущих научных изданиях и журналах сравнительно большого числа работ, посвященных этим исследованиям, например [1–7]. Актуальность такого рода исследований объясняется широчайшим спектром приложений таких систем к различным областям знаний, а также тем очевидным фактом, что шумы являются неотъемлемой частью окружающего нас мира и часто играют решающую роль в изменении режимов поведения динамических систем. В работе [8] авторами развита теория, позволяющая с единой точки зрения концепции параметров порядка провести последовательное и детальное изучение процессов образования пространственных диссипативных структур (паттернов), спонтанно возникающих в открытых нелинейных пространственно распределенных системах с внешними шумами как в окрестности, так и вдали от точки перехода. Настоящая работа является продолжением исследований, начатых авторами в [8–11], и посвящена изучению шумоиндуцированных упорядочивающего и разупорядочивающего фазовых переходов, возникающих в конкретной биофизической системе в области параметров детерминированной системы, соответствующих бифуркации Тьюринга. 81 М а к с и м о в В. В. 1. Математическая модель. Математическая модель изучаемой биофизической системы описывается уравнениями [12–14] ∂x1 ax1 x2 = rx1 (1 − x1 ) − + D1 2 x1 , ∂t 1 + bx1 ax1 x2 g 2 x2 ∂x2 2 = − mx2 − f + D2 ∂t 1 + bx1 1 + h2 x2 2 2 (1) x2 , где x1 , x2 — функции состояния; параметры r, a, b, m, g, h, f , D1 и D2 описаны подробно в [12, 14]. Исследование локальной динамики и бифуркационный анализ системы (1) представлены в [13, 14]. Введём безразмерные время τ = rt, координаты x = x r/D1 и представим параметры m/r и a/r в виде m/r = (m0 /r0 )(1 + f1 (x , τ )), a/r = (a0 /r0 )(1 + f2 (x , τ )). Здесь m0 , r0 , a0 — пространственно-временные средние соответствующих параметров, случайные однородные изотропные гауссовы поля fi (x , τ ) определяют пространственно-временные флуктуации этих параметров и имеют нулевые средние и корреляционные функции вида K fj (r, t), fj (r , t ) = Φj ( r − r ) exp(−ktj t − t )δjj , где Φj (|r − r |) = θj exp(−kf j |r − r |), θj — интенсивности флуктуаций, kf j и ktj определяют их характерные пространственный и временной масштабы. Принимая во внимание внешний шум, получим ∂x1 a0 x1 x2 = x1 (1 − x1 ) − (1 + f2 (x , τ )) + D1 2 x1 , ∂τ r0 1 + bx1 ∂x2 a0 x1 x2 m0 = (1 + f2 (x , τ )) − (1 + f1 (x , τ ))x2 − ∂τ r0 1 + bx1 r0 g 2 x2 2 f + D2 2 x2 . − 1 + h2 x2 2 (2) 2. Результаты аналитического исследования. В работе [8] было получено уравнение Фоккера—Планка для критического параметра порядка стохастических систем реакционно-диффузионного типа. В этом разделе представлены результаты аналитического исследования системы (2), полученные на основании этого уравнения при следующих значениях параметров: r0 = 1, a0 = 8, g = 1,434, f = 0,093, h = 0,857, b = 11,905, m0 = 0,490, rf 1 = = rf 2 = 1, θ1 = 2,4 · 10−5 , критическое значение контрольного параметра D = 135. Известно, что в докритической области состояние любой системы является однородным статистически стационарным (беспорядок) и этому состоянию соответствует одномодальная плотность распределения вероятности. Возникновение неоднородного статистически стационарного состояния (порядок) проявляется в расщеплении максимума плотности распределения вероятности на два симметричных. На рис. 1 и 2 представлено изменение плотности стационарного распределения вероятности амплитуды критической моды ξkc с увеличением интенсивности шума θ2 при переходе через точку бифуркации детерминированной системы (1). 82 Переход «беспорядок – порядок – беспорядок» . . . Рис. 1. Стационарная плотность распределения вероятности значений амплитуды критической моды системы (2) в докритической области для двух значений интенсивности шума θ2 : a) θ2 = 6,4 · 10−5 , б) θ2 = 3,2 · 10−5 Рис. 1 иллюстрирует плотность распределения вероятности амплитуды критической моды в докритической области. При малых шумах стационарная плотность вероятности близка к δ-функции, и среднее и наиболее вероятное значения модуля амплитуды критической моды (параметра порядка) совпадают и равны нулю, т.е. однородное статистически стационарное состояние системы является наиболее вероятным (см. рис. 1, a). При увеличении интенсивности шума (см. рис. 1, б) происходит деформация кривой стационарной плотности вероятности: максимум, не смещаясь, значительно уменьшается, при этом основание кривой расширяется. Так как полученное распределение вероятности не является Гауссовым, среднее значение модуля амплитуды становится отличным от наиболее вероятного. Таким образом, несмотря на отсутствие расщепления одномодальной плотности вероятности на бимодальную, среднее значение параметра порядка становится отличным от нуля и следует ожидать возникновения неоднородного статистически стационарного состояния (порядка), вероятность которого мала. Очевидно, что в рассматриваемом случае возникновение нового состояния носит случайный характер. Объяснить это можно так. Формирование структур в докритической области может происходить на сильных неоднородностях среды. Такую неоднородность может создать сильная (крупномасштабная) флуктуация, вероятность которой невелика и зависит от параметров шума. Таким образом, при длительном наблюдении за эволюцией системы и подходящих параметрах внешнего шума можно ожидать формирования таких случайных неоднородностей, которые вызовут возникновение структур в области их нахождения. Рис. 2 демонстрирует стационарную плотность вероятности амплитуды критической моды системы (2) в закритической области при различных интенсивностях шума. Бимодальные распределения соответствуют существованию диссипативных структур. При этом наиболее вероятное значение и математическое ожидание параметра порядка становится отличным от нуля. Рис. 2 ясно показывает, что при увеличении интенсивности шума происходит 83 М а к с и м о в В. В. Рис. 2. Стационарная плотность распределения вероятности значений амплитуды критической моды системы (2) в закритической области для пяти значений интенсивности шума θ2 : a) θ2 = 0, б) θ2 = 6,4 · 10−5 , в) θ2 = 6,4 · 10−3 , г) θ2 = 1,6 · 10−1 , д) θ2 = 6,4 · 10−1 постепенное слияние максимумов, и при некотором критическом значении интенсивности шума вновь возникает одномодальная плотность. При этом |ξkc mp | = 0, а |ξkc | = 0. Система (2) переходит в состояние сильно нерегулярного поведения (беспорядок). Таким образом, полученное изменение плотности стационарного распределения вероятности амплитуды критической моды свидетельствует о существовании в системе (2) фазового перехода «беспорядок – порядок – беспорядок». Обсуждавшееся выше изменение статистически стационарных среднего и наиболее вероятного значений модуля амплитуды 84 Переход «беспорядок – порядок – беспорядок» . . . критической моды, соответствующих изменению плотностей, изображенных на рис. 2, демонстрирует рис. 3. На рис. 4 представлена граница шумоиндуцированного фазового перехода «порядок – беспорядок» для системы (2), предсказанная с применением подхода, развитого в [8]. Следует отметить, что при приближении к детерминированной точке перехода даже очень малые флуктуации будут способствовать потере устойчивости неоднородного состояния и вызывать неупорядоченное состояние. Определим относительные флуктуации параметра порядка (восприимчи2 вость) в соответствии с работой [15] как χ = [ ξkc − ξkc 2 ]/θ2 и его кумулянт 2 второго порядка κ2 = ξkc / ξkc 2 . На рис. 5 изображены графики χ и κ2 как функции интенсивности шума θ2 . Наличие максимумов восприимчивости ясно показывает увеличение флуктуаций в окрестности двух критических точек. Качественно вид кривых κ2 и χ, полученных для системы (2), совпадает в соответствующей области с аналогичными кривыми, полученными численно в работе [15] для системы, в которой наблюдается чисто индуцированный шумом переход. Это свидетельствует о том, что в системе (2) шумоиндуцированный переход «порядок – беспорядок» имеет ту же природу. Рис. 3. Статистически стационарныe среднее |ξkc | (сплошная линия) и наиболее вероятное |ξkc mp | (штриховая линия) значения модуля амплитуды критической моды в зависимости от интенсивности шума θ2 (закритическая область) Рис. 4. Граница перехода «порядок – беспорядок» на плоскости параметров D и θ2 85 М а к с и м о в В. В. Рис. 5. Восприимчивость χ (сплошная линия) и кумулянт второго порядка κ2 (штриховая линия) как функции интенсивности шума θ2 Заключение. В работе исследована эволюция двухкомпонентной биофизической системы реакционно-диффузионного типа в поле внешних флуктуаций, имеющая важное теоретическое и практическое значение. Изучено поведение статистически стационарных характеристик первого и второго порядков амплитуды критической моды этой системы. Определена граница перехода «порядок – беспорядок». Исследование восприимчивости и кумулянта второго порядка абсолютного значения амплитуды критической моды позволяет сделать вывод о том, что шумоиндуцированные упорядочивающий и разупорядочивающий переходы, связанные с образованием диссипативных структур в этой системе, аналогичны чисто индуцированным шумом переходам и являются неравновесными фазовыми переходами второго рода. Результаты численного моделирования эволюции системы (2) [11] имеют удовлетворительное количественное соответствие результатам аналитического исследования в окрестности точки перехода «беспорядок – порядок».

About the authors

Valerii V Maksimov

Samara State Transport University

Email: vvmaksimov52@mail.ru
(Ph. D. (Techn.)), Dept. of Higher Mathematics 18, First Bezimyanniy per., Samara, 443066, Russia

References

Supplementary files