Здесь а – некоторый эффективный коэффициент, зависящий от значения управляющего параметра α. Система становится неустойчивой при а > 0.
186
3.5.1. Простейшая система с двумя стационарными состояниями, бифуркацией и неустойчивостью стационарного состояния
Рассмотрим кинетически необратимый брутто-процесс
R → P,
осуществляющийся по нелинейной схеме
ε1
R + Хε2
2 Х, (3.22) Х → P,
где R и P – исходный реагент и конечный продукт соответственно, а Х – интермедиат-автокатализатор. Введение кинетической необратимости во второй стадии означает, что эта стадия далека от термодинамического равновесия и поэтому можно ожидать неустойчивости стационарных состояний. Проверим это.
Для эволюции концентрации интермедиата
|
|
|
d[X] |
|
~ ~ ~2 |
~ |
|
|
|
|
|
|
dt |
= ε1(R X − X |
) − ε2 X . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отсюда видно, что стационарному состоянию отвечают два |
||||||||
решения: |
|
|
|
|
|
|
|
||
|
~ |
= 0; |
|
|
|
|
|
|
|
|
1) X |
|
|
|
|
|
|
||
|
~ |
~ |
ε2 |
|
|
|
|
|
|
|
2) X |
= R − ε1 . |
|
|
~ |
|
|
|
|
|
Физический смысл имеет только решение X ≥ 0. |
|
|
||||||
|
|
|
|
~ |
~ |
ε1 существует только од- |
|||
|
Таким образом, при R < Rкр ≡ ε2 |
||||||||
|
|
|
|
|
~ |
~ |
~ |
= ε2 |
ε1 |
но действительное решение X = 0 . Однако при R < Rкр |
|||||||||
возможны |
два действительных |
решения: |
~ |
= 0 |
и |
||||
X1 |
|||||||||
~ |
~ |
~ |
|
|
|
~ |
~ |
|
|
X2 |
= R −Rкр |
> 0 (рис. 3.2). Таким образом, точка R = Rкр явля- |
|||||||
ется точкой бифуркации (раздвоения) стационарных решений.
186
|
|
|
|
|
|
|
Рис. 3.2. Зависимость термоди- |
||
|
|
|
|
|
|
|
намического напора интерме- |
||
|
|
|
|
|
|
|
диата Х для стационарных со- |
||
|
|
|
|
|
|
|
стояний схемы (3.22) от управ- |
||
|
|
|
|
|
|
|
ляющего параметра – термоди- |
||
|
|
|
|
|
|
|
намического напора исходного |
||
|
|
|
|
|
|
|
|
~ |
~ |
|
|
|
|
|
|
|
реагента R . Точка Rкр = ε2 ε1 |
||
|
|
|
|
|
|
|
– точка бифуркации стационар- |
||
|
|
|
|
|
|
|
ных состояний |
|
|
Покажем, |
что в условиях существования двух решений |
||||||||
|
~ |
~ |
решение |
~ |
= 0 является неустойчивым. |
||||
(т. е. при R > Rкр ) |
X1 |
||||||||
|
|
|
~ ~ ~ |
2 |
|
~ |
~ |
~ ~ |
~ |
P = ∑viAi = ε1(R X − X |
)(lnRR −ln X) + ε2 X(ln X −lnP) . |
||||||||
~ |
i |
|
|
|
|
|
|
|
|
= 0 в стационарном состоянии Р1 = 0. |
|
||||||||
При X |
|
||||||||
~ |
~ |
−(ε2 |
ε1) > 0 в стационарном состоянии Р2 > 0. |
||||||
При X |
= R |
||||||||
Рассмотрим флуктуацию скорости диссипации энергии при
~
флуктуации X из стационарного значения.
δP = ∑δvi δAi = δv1 δA1 + δv2 δA2 =
|
i |
|
|
|
|
|
|
|
|
|
|
|
~ ~ |
~ |
~ |
~ |
~ |
~ |
|
|
|
|
|
= ε1(RδX − 2nXδX) |
(−δln X) + ε2δX δln X = |
|
|
|
|||||||
= {− ε1(R − 2X) + ε2 |
}(δX δln X)= {ε1(R − |
2X)+ ε2 } |
|
~ 2 |
. |
|
|||||
|
~ |
|
|||||||||
~ |
~ |
|
~ |
~ |
~ |
~ |
|
(δX) |
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
1. Для стационарного решения X = 0 |
|
|
|
|
|
|
|||||
|
|
|
~ |
~ |
~ |
|
|
|
|
|
|
|
δP = (−ε1R + ε2 )(δX |
δln X) . |
|
|
|
|
|
||||
Таким образом, при |
~ |
< ε2 |
ε1 δP |
> 0 |
и рассматриваемое |
||||||
R |
|||||||||||
стационарное |
решение |
устойчиво. |
Однако |
при |
~ |
|
ε1 |
||||
R > ε2 |
|||||||||||
δP < 0, т. е. в этой ситуации стационарное состояние неустойчиво.
187
2. Для стационарного решения |
~ ~ |
ε2 |
||||
X = R − |
ε1 |
|||||
|
|
~ |
ε2 |
~ |
~ |
|
|
δP = { − ε1( −R + |
ε |
) + ε2 |
} (δX δln X) = |
||
|
|
~ |
|
1 |
~ |
|
|
|
|
~ |
|
||
|
|
= {ε1R − ε2 |
} (δX δln X) > 0, |
|
||
~ |
ε2 |
это состояние устойчиво. |
|
|||
т. е. при R > |
ε |
|
||||
|
1 |
|
|
|
|
|
Существенно, что неустойчивость стационарного состояния была получена вследствие предположения о необратимости реакции получения конечного продукта Р.
Вывод о наличии бифуркации в автокаталитических системах приводит, в частности, к выводу о возможности существования аналога биологического естественного отбора даже в абиогенных автокаталитических системах.
Действительно, пусть существует сразу несколько популяций автокатализаторов Хi, получаемых~ из одной и той же
«пищи» R. При уменьшении значения R происходит полное «вымирание» всех автокатализаторов с наибольшими
~ |
~ |
Rкрi = ε2i |
ε1i < R . При этом при последующем увеличении ко- |
~
личества «пищи» (т. е. при росте R ) вымершие популяции не восстанавливаются, поскольку нет затравок для этих автокатализаторов. Однако если возможна мутация~автокатализато-
ров с уменьшением критического параметра Rкрi, то при уве-
~
личении R популяция этого автокатализатора начнет расти.
3.5.2. Множественность стационарных состояний. Пример простейшей «триггерной» системы с тремя стационарными состояниями
Приведенный выше пример демонстрирует возможность наличия у открытой системы двух, т. е. множественности стационарных состояний. Множественность стационарных состояний может быть и более сложной.
188
Рассмотрим кинетически необратимый брутто-процесс
R → P,
который протекает по схеме |
|
ε1 |
|
|||
R + 2 Y |
3 Y, |
|||||
|
|
|
||||
|
|
|
||||
|
|
ε2 |
(2.23) |
|||
|
Y → |
P, |
||||
где Y – реакционный интермедиат.
Кинетика изменения концентрации этого интермедиата
описывается уравнением |
|
|
|
|
|
|||
d[Y] |
~ ~ |
2 |
~ |
3 |
~ |
~ ~ ~ |
2 |
~ |
dt |
= ε1(R Y |
|
− Y |
|
) − ε2Y ≡ (ε1R Y − ε1Y |
|
− ε2 )Y . |
|
|
|
|
|
|
|
|
|
|
Очевидно, что в данной ситуации имеются три стационарных решения:
|
|
|
|
|
|
|
|
|
|
~ |
|
= 0, |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
Y1 |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
~ |
|
|
1 |
~ |
|
1 |
|
~ |
2 |
|
ε2 |
|
|
||||||
|
|
|
|
|
Y2,3 = |
|
|
R ± |
|
|
|
|
R |
|
− |
|
|
. |
|
|||||
|
|
|
|
|
2 |
4 |
|
ε1 |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Это |
значит, что |
при |
~ |
|
|
~ |
≡ 2 |
|
|
ε2 |
|
у |
|
системы |
существует |
|||||||||
R |
< Rкр |
|
|
ε1 |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|||
только |
одно |
стационарное состояние |
|
однако при |
||||||||||||||||||||
с Y = 0 , |
||||||||||||||||||||||||
~ |
~ |
|
ε2 |
у |
системы |
существуют |
|
три действительных |
||||||||||||||||
R > Rкр ≡ 2 |
ε1 |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
решения: |
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
= 0 , |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
Y1 |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
~ |
|
|
1 |
|
~ |
|
1 |
|
~ |
2 |
|
ε2 |
|
|
|
||||
|
|
|
|
|
Y2 |
= |
|
|
R − |
|
|
|
|
R |
|
|
− |
|
|
, |
|
|||
|
|
|
|
|
|
2 |
|
4 |
|
|
ε1 |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
~ |
|
|
1 ~ |
|
1 ~ |
2 |
|
ε2 |
|
|
|
||||||||
|
|
|
|
|
Y3 |
= |
|
|
R + |
|
|
|
|
R |
|
|
− |
|
|
, |
|
|||
|
|
|
|
|
|
2 |
|
4 |
|
|
ε1 |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
189