1. Введение: Что, если молоко, налитое в кофе, начнёт пульсировать само по себе?

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

Однако в физической химии существует реакция, которая эффектно опровергает эту бытовую интуицию: реакция Белоусова — Жаботинского (реакция БЖ). При смешивании нескольких прозрачных водных реагентов в стакане раствор спонтанно и с точностью метронома десятки раз меняет цвет туда и обратно; а если налить его тонким неперемешиваемым слоем в чашку Петри, он сам по себе порождает расходящиеся концентрические кольца (пейсмейкеры) и вращающиеся архимедовы спиральные волны.

1951 год: «Еретическое» открытие, отвергнутое советской наукой

История началась в 1951 году в Москве. Советский химик и биофизик Борис Павлович Белоусов (1893–1970), работавший в лаборатории биофизики Минздрава СССР, искал неорганический аналог биологического цикла трикарбоновых кислот (цикла Кребса). Он растворил в разбавленной серной кислоте бромат калия ($\text{KBrO}_3$) как окислитель, лимонную кислоту как органический субстрат и сульфат церия ($\text{Ce}^{4+}/\text{Ce}^{3+}$) как окислительно-восстановительный катализатор. Вместо монотонного перехода от жёлтого ($\text{Ce}^{4+}$) к бесцветному ($\text{Ce}^{3+}$) раствор Белоусова более часа периодически пульсировал между жёлтым и бесцветным состояниями!

Когда Белоусов в 1951 и 1957 годах направил рукописи в ведущие советские химические журналы, они были категорически отвергнуты. Вердикт рецензентов гласил:

«Утверждение о периодических колебаниях гомогенной химической реакции туда и обратно через точку равновесия прямо противоречит Второму началу термодинамики. Наблюдения автора — это артефакт плохого перемешивания, пыли или пузырьков газа».

Глубоко оскорблённый 64-летний учёный опубликовал в 1958 году лишь двухстраничный реферат в малоизвестном «Сборнике рефератов по радиационной медицине» и навсегда оставил эту тему.

1961 год: Реабилитация Анатолием Жаботинским и красно-синее биение

В 1961 году профессор МГУ Симон Эльевич Шноль передал рецепт молодому аспиранту Анатолию Марковичу Жаботинскому (1938–2008). Жаботинский заменил лимонную кислоту на малоновую кислоту ($\text{CH}_2(\text{COOH})_2$) и ввёл редокс-индикатор ферроин ($[\text{Fe}(\text{phen})_3]^{2+}/[\text{Fe}(\text{phen})_3]^{3+}$), благодаря чему раствор стал ярко пульсировать между кирпично-красным ($[\text{Fe}(\text{phen})_3]^{2+}$) и небесно-голубым ($[\text{Fe}(\text{phen})_3]^{3+}$) цветами. В 1980 году Белоусову (посмертно), Жаботинскому, Заикину, Кринскому и Иваницкому была присуждена Ленинская премия.

2. Разрешение термодинамического парадокса: почему реакция БЖ не нарушает Второе начало

Монотонное убывание энергии Гиббса и соотношения Онсагера вблизи равновесия

В закрытой системе при постоянных температуре $T$ и давлении $P$ Второе начало требует неотрицательного производства внутренней энтропии $\sigma_S \equiv \frac{d_i S}{dt} \ge 0$, поэтому энергия Гиббса $G = H - TS$ убывает монотонно:

где $v_{\rho}$ — скорость реакции $\rho$, а $A_{\rho} \equiv -\sum_i \nu_{i\rho} \mu_i = -\left(\frac{\partial G}{\partial \xi_{\rho}}\right)_{T,P}$ — химическое сродство Де Дондера.

В линейной области вблизи равновесия ($|A_{\rho}| \ll RT$) микроскопическая обратимость во времени требует принципа детального равновесия и соотношений взаимности Онсагера (Ларс Онсагер, Нобелевская премия 1968 г.), в силу которых матрица кинетических коэффициентов строго симметрична ($L_{\alpha\beta} = L_{\beta\alpha}$). Поэтому матрица Якоби $\mathbf{J}_{\text{eq}}$ подобна вещественной симметричной отрицательно полуопределённой матрице, и все её собственные значения чисто вещественны и неположительны ($\lambda_k \in \mathbb{R}, \lambda_k \le 0, \text{Im}(\lambda_k) = 0$). Вблизи термодинамического равновесия невозможны даже затухающие колебания!

Вдали от равновесия: нарушение симметрии и диссипативные структуры Пригожина

Разгадка парадокса состоит в том, что реакция БЖ колеблется не вокруг точки термодинамического равновесия, а протекает вдали от равновесия ($A_{\rho} \gg RT$) в сильно нелинейном режиме!

Вдали от равновесия автокаталитическая обратная связь нарушает симметрию Онсагера, матрица Якоби $\mathbf{J}(\mathbf{c}^*)$ становится несимметричной и может иметь комплексно-сопряжённые собственные значения $\lambda_{\pm} = \mu \pm i\omega$ ($\omega \neq 0$). Бельгийский физикохимик Илья Пригожин (Нобелевская премия 1977 г.) показал, что вдали от равновесия диссипация энергии ($\sigma_S > 0$) становится источником пространственно-временной самоорганизации — диссипативных структур:

1. В открытом проточном реакторе (CSTR): Приток свежих реагентов компенсирует производство энтропии ($dS = d_e S + d_i S = 0$), поддерживая незатухающие колебания бесконечно долго.
2. В закрытом стакане (ступенчатый спуск $G(t)$): Высококонцентрированные основные реагенты ($\text{BrO}_3^-$ и малоновая кислота, $\sim 10^{-1}\text{ M}$) расходуются строго монотонно ($\frac{dG}{dt} \le 0$). Замкнутые петли в пространстве концентраций описывают только следовые интермедиаты ($\text{HBrO}_2, \text{Br}^-, \text{Ce}^{4+}$, $\sim 10^{-7}\text{–}10^{-4}\text{ M}$) — в точности как водяное колесо в водопаде, которое вращается сотни раз, пока вода течёт строго сверху вниз.

3. Механизм FKN (1972): три сопряжённых процесса и автокаталитический переключатель

В 1972 году Ричард Филд, Эндре Кёрёш и Ричард Нойес расшифровали 18-стадийный механизм FKN окисления и бромирования малоновой кислоты броматом:

Цикл переключается между процессами A, B и C благодаря конкуренции между активатором — бромистой кислотой ($\text{HBrO}_2$) и ингибитором — бромид-ионом ($\text{Br}^-$):

  flowchart TD
    subgraph ProcessA["Процесс A: Нерадикальное расходование Br⁻ (Восстановленное состояние: Красный / Бесцветный)"]
        A1["Концентрация бромида высока: [Br⁻] > [Br⁻]_crit"]
        A2["Ионная реакция: BrO₃⁻ + 2Br⁻ + 3H⁺ → 3HOBr"]
        A3["HOBr бромирует малоновую кислоту в BrMA; [Br⁻] падает"]
        A1 --> A2 --> A3
    end

    subgraph Switch1["Автокаталитический переключатель зажигания"]
        S1["[Br⁻] падает ниже порога [Br⁻]_crit! Тормоз снят"]
    end

    subgraph ProcessB["Процесс B: Автокаталитический взрыв и окисление катализатора (Окисленное состояние: Синий / Жёлтый)"]
        B1["Радикальная стадия: BrO₃⁻ + HBrO₂ + H⁺ → 2BrO₂• + H₂O"]
        B2["Окисление катализатора: 2BrO₂• + 2Ce³⁺ + 2H⁺ → 2HBrO₂ + 2Ce⁴⁺"]
        B3["Положительная обратная связь: 1 молекула HBrO₂ рождает 2 HBrO₂!"]
        B4["Скачок цвета от Красного (Fe²⁺/Ce³⁺) к Синему (Fe³⁺/Ce⁴⁺)"]
        B1 --> B2 --> B3 --> B4
    end

    subgraph ProcessC["Процесс C: Запаздывающая отрицательная обратная связь и сброс"]
        C1["Накопленный Ce⁴⁺ окисляет малоновую и броммалоновую кислоту (BrMA)"]
        C2["Ce⁴⁺ → Ce³⁺ с выделением бромид-ионов (f Br⁻)"]
        C3["[Br⁻] превышает верхний порог и гасит HBrO₂"]
        C1 --> C2 --> C3
    end

    A3 --> S1 --> B1
    B4 --> C1
    C3 --> A1

1. Процесс A (Бромидный тормоз): Пока $[\text{Br}^-] > [\text{Br}^-]_{\text{crit}}$, ионы $\text{Br}^-$ мгновенно перехватывают $\text{HBrO}_2$ по реакции $\text{HBrO}_2 + \text{Br}^- + \text{H}^+ \xrightarrow{k_2} 2\text{HOBr}$ с образованием броммалоновой кислоты ($\text{BrMA}$). Катализатор остаётся восстановленным (красным).
2. Процесс B (Автокаталитический взрыв): Как только $[\text{Br}^-] < [\text{Br}^-]_{\text{crit}} = \frac{k_5}{k_2}[\text{BrO}_3^-] \approx 5 \times 10^{-7}\text{ M}$, запускается радикальная цепь через $\text{BrO}_2^\bullet$: $$\text{BrO}_3^- + \mathbf{HBrO_2} + 3\text{H}^+ + 2\text{Ce}^{3+} \longrightarrow \mathbf{2\,\text{HBrO}_2} + 2\text{Ce}^{4+} + \text{H}_2\text{O}$$ Каждая молекула $\text{HBrO}_2$ порождает две молекулы $\text{HBrO}_2$! Концентрация $[\text{HBrO}_2]$ взлетает на четыре порядка, окисляя катализатор ($\text{Ce}^{3+} \to \text{Ce}^{4+}$, переход в синий цвет), и ограничивается диспропорционированием $2\text{HBrO}_2 \xrightarrow{k_4} \text{BrO}_3^- + \text{HOBr} + \text{H}^+$.
3. Процесс C (Запаздывающая отрицательная обратная связь): Накопленный $\text{Ce}^{4+}$ окисляет броммалоновую кислоту, восстанавливаясь до $\text{Ce}^{3+}$ и высвобождая $f$ ионов $\text{Br}^-$, что гасит автокатализ и возвращает систему к Процессу A.

4. Математическая теория орегонатора: релаксационные колебания, бифуркация Хопфа и предельный цикл

В 1974 году Филд и Нойес свели механизм FKN к системе из трёх безразмерных ОДУ орегонатора (Oregonator) для $x \propto [\text{HBrO}_2]$, $y \propto [\text{Br}^-]$, $z \propto [\text{Ce}^{4+}]$:

• Разделение времён и релаксационные колебания: Поскольку $\varepsilon' \sim 10^{-4} \ll \varepsilon \sim 10^{-2} \ll 1$, квазистационарное приближение $y \approx \frac{f z}{x + q}$ сводит систему к плоскости $(x, z)$, где нуль-изоклина $z = \frac{x(1-x)(x+q)}{f(x-q)}$ имеет кубическую S-образную (или N-образную) складчатую форму. Траектория медленно ползёт вдоль устойчивых ветвей и скачком перепрыгивает между ними в точках срыва, порождая пилообразные релаксационные колебания.
• Бифуркация Хопфа и теорема Пуанкаре — Бендиксона: В стационарной точке $(x^*, z^*)$ всегда $\det(\mathbf{J}) > 0$. Когда стехиометрический фактор $f$ лежит в окне колебаний ($0{,}5 \lesssim f \lesssim 1 + \sqrt{2} \approx 2{,}414$), точка $(x^*, z^*)$ попадает на неустойчивую среднюю ветвь с $\text{Tr}(\mathbf{J}) > 0$. Переход собственных значений $\lambda_{\pm} = \mu(f) \pm i\omega(f)$ через мнимую ось порождает бифуркацию Андронова — Хопфа, а теорема Пуанкаре — Бендиксона гарантирует существование глобально притягивающего предельного цикла.

5. Химические волны в неперемешиваемой среде: УРЧП реакции-диффузии и спиральные волны

В тонком слое раствора БЖ в чашке Петри локальная кинетика $\mathbf{R}(\mathbf{c})$ объединяется со вторым законом диффузии Фика в систему уравнений в частных производных (УРЧП) реакции-диффузии:

1. Возбудимая среда и формула Лютера: Автоволны распространяются без затухания с постоянной скоростью, определяемой формулой Лютера (1906): $$v_0 = 2 \sqrt{D_X \, k_{\text{auto}}} = 2 \sqrt{D_X \, k_5 [\text{BrO}_3^-][\text{H}^+]} \approx 5\text{–}10\text{ мм/мин}$$
2. Взаимная аннигиляция и фазовая сингулярность спиральных волн: Сразу за синим фронтом волны тянется невозбудимый рефрактерный хвост с высокой концентрацией $\text{Br}^-$. Поэтому при лобовом столкновении две химические автоволны полностью гасят (аннигилируют) друг друга! Если разорвать кольцевой фронт волны стеклянной палочкой, свободный конец закручивается вокруг собственного рефрактерного хвоста в самоподдерживающуюся архимедову спиральную волну вокруг топологической фазовой сингулярности ($\oint_C \nabla \phi \cdot d\mathbf{l} = \pm 2\pi n$), подчиняясь уравнению эйконала $v_n = v_0 - D\kappa$.

6. От живых систем к умным материалам: структуры Тьюринга, фибрилляция сердца и самоходные гели

  flowchart TD
    RD["Неравновесная система реакции-диффузии: ∂c/∂t = D∇²c + R(c)"]

    subgraph Temporal["Временная самоорганизация (Перемешивание / D = 0)"]
        T1["Бифуркация Хопфа и колебания предельного цикла"]
        T2["Биология: Колебания гликолиза (NADH), циркадные ритмы, клеточный цикл"]
        T1 --> T2
    end

    subgraph Excitable["Пространственно-временные автоволны (D_активатора ≈ D_ингибитора)"]
        E1["Концентрические волны и вращающиеся спиральные волны (фазовая сингулярность)"]
        E2["Биология: Волны цАМФ у Dictyostelium, потенциал действия миокарда и фибрилляция"]
        E1 --> E2
    end

    subgraph Turing["Стационарные пространственные структуры (D_ингибитора >> D_активатора)"]
        S1["Диффузионная неустойчивость (Бифуркация Тьюринга)"]
        S2["Биология: Окраска животных (пятна леопарда, полоски зебры) и морфогенез"]
        S1 --> S2
    end

    subgraph SmartMat["Биомиметическая мягкая робототехника и хемокомпьютинг"]
        M1["Самоколеблющиеся гидрогели (Рё Ёсида): Шагающие искусственные мышцы без батарей"]
        M2["Светочувствительные логические вентили Ru(bpy)₃²⁺ и поиск пути в лабиринте"]
    end

    RD --> Temporal
    RD --> Excitable
    RD --> Turing
    Excitable --> SmartMat

• Морфогенез Алана Тьюринга (1952): Когда ингибитор диффундирует намного быстрее активатора ($D_v \gg D_u$), возникают стационарные пятна и полосы (структуры Тьюринга), экспериментально полученные в микроэмульсиях БЖ-АОТ.
• Фибрилляция желудочков сердца: Сердечная мышца (миокард) — это электрически возбудимая среда. Разрыв фронта электрической волны на ишемическом участке порождает вращающиеся спиральные волны (3D-скролл-волны), распад которых ведёт к смертельной фибрилляции желудочков. Разряд дефибриллятора (АНД) одновременно переводит весь миокард в рефрактерное состояние, уничтожая фазовые сингулярности.
• Самоколеблющиеся гидрогели (Рё Ёсида): Ковалентное встраивание катализатора $\text{Ru(bpy)}_3^{2+}$ в термочувствительную полимерную сетку (PNIPAAm) создаёт самоходные гидрогели, которые преобразуют редокс-колебания БЖ в перистальтические сокращения и ползут подобно гусенице без внешних источников тока.