振動反応（BZ反応など）を反応速度論で解く

問いは木下が出したもので、答えは Claude Code と調べて書いたものです。

Q. 信号反応・ヨウ素時計反応（2026-10-18）の続きで、振動反応（色や濃度が周期的に行ったり来たりする反応）を反応速度論から解くと？

A. まず前提の整理から。熱力学的には、閉じた系の反応は自由エネルギーが単調に減るので「平衡のまわりを行き来する」ことはできない（第二法則。平衡付近の濃度振動は禁止）。振動反応が起きるのは、反応物の大きな在庫が単調に減り続ける途中で、中間体の濃度だけが振動するから。見た目は往復でも、全体の自由エネルギーは一方通行で坂を下っている。振動するのは系全体でなく一部の変数（中間体・触媒の酸化状態）だけ、というのが肝。

振動に必要な3条件：①平衡から遠い（反応物が豊富／流通系で供給される）、②自己触媒（正のフィードバック）＝ある中間体が自分自身の生成を速める、③遅れて効く負のフィードバック＝増えた中間体が別の経路で自分を抑える。②で急増、③で急減、これが交互に起きて振動する。数学的には不動点が不安定化してリミットサイクル（閉じた軌道）に落ちること。

BZ反応（ベロウソフ–ジャボチンスキー反応）の中身：酸性でマロン酸を臭素酸 $\ce{BrO3^-}$ が酸化する反応を、金属触媒（セリウム $\ce{Ce^{3+}/Ce^{4+}}$ やフェロイン）が仲立ちする。触媒の酸化状態で色が変わり（$\ce{Ce^{3+}}$ 無色↔$\ce{Ce^{4+}}$ 黄、フェロインなら赤↔青）、周期的に色が切り替わる。機構（FKN 機構）は臭化物イオン $\ce{Br^-}$ の多寡で切り替わる3つの過程からなる：

• 過程A（$\ce{Br^-}$ が多いとき）：$\ce{BrO3^-}$ が $\ce{Br^-}$ を消費していく。$\ce{Br^-}$ が減って閾値を割ると —
• 過程B（$\ce{Br^-}$ が少ないとき）：亜臭素酸 $\ce{HBrO2}$ の自己触媒的な増殖 $\ce{BrO3^- + HBrO2 + H^+ -> 2 HBrO2 + \text{(酸化)}}$ が暴走し、触媒が一気に $\ce{Ce^{4+}}$ へ酸化される（色が変わる）。
• 過程C：$\ce{Ce^{4+}}$ がマロン酸を酸化しながら $\ce{Br^-}$ を再生する。$\ce{Br^-}$ が増えると過程Aへ戻る。

$\ce{Br^-}$ が「スイッチ」、$\ce{HBrO2}$ が自己触媒の「アクセル」、$\ce{Ce^{4+}}$ 経由の $\ce{Br^-}$ 再生が「遅れて効くブレーキ」。この3役が条件②③そのもの。

最小の数理モデルで振動を出す。実際の FKN を3変数に縮めたのがオレゴネーターだが、振動が「なぜ起きるか」は2変数モデルで十分見える。

(1) 最初の化学振動子モデル＝ロトカ機構（Lotka 1920）。$\ce{A + X -> 2X}$（速度 $k_1[\ce{A}]x$）、$\ce{X + Y -> 2Y}$（$k_2xy$）、$\ce{Y -> P}$（$k_3y$）で $[\ce{A}]$ 一定とすると、これは捕食者–被食者（Lotka–Volterra）と同じ式になり振動する。ただし不動点が中心（中立安定）なので、少し揺らすと振幅が勝手に変わる脆い振動。本物の BZ は外乱に強い（振幅が決まった）リミットサイクルなので、これでは不足。

(2) リミットサイクルをきれいに出すブラッセレーター。$\ce{A -> X}$、$\ce{2X + Y -> 3X}$（3次の自己触媒）、$\ce{B + X -> Y + D}$、$\ce{X -> E}$。無次元化すると

$$\frac{dx}{dt}=a-(b+1)x+x^2y,\qquad \frac{dy}{dt}=bx-x^2y$$

不動点は $x^*=a,\ y^*=b/a$。ヤコビ行列は

$$J=\begin{pmatrix} b-1 & a^2 \\ -b & -a^2 \end{pmatrix},\quad \operatorname{tr}J=b-1-a^2,\quad \det J=a^2>0$$

$\det J>0$ は常に成り立つので、不動点の安定性はトレースの符号だけで決まる。$\operatorname{tr}J$ が $0$ を正へ横切る $b=1+a^2$ がホップ分岐で、

$$b>1+a^2 \ \Longrightarrow\ \text{不動点は不安定} \ \Longrightarrow\ \text{リミットサイクル（持続振動）}$$

となる。有界な正領域に不安定な不動点が1つだけなので、ポアンカレ–ベンディクソンの定理が閉軌道の存在を保証する。$\operatorname{tr}J$ の $+2xy=2b$ の項（＝自己触媒 $\ce{2X+Y->3X}$ の寄与）がトレースを正へ押し上げており、振動の正体は自己触媒による不動点の不安定化だと式の上で見える。

時計反応・信号反応との対比：ヨウ素時計は「在庫が尽きた1回だけ」急変する（1発の閾値）。信号反応は2つの速度の綱引きで、$\ce{O2}$ を使い切れば止まる（振れば戻せるが自発的には往復しない）。振動反応はフィードバックで自発的に何度も往復する点が本質的に違う。空間に広げると濃度の波（BZ の渦巻きパターン）になり、これは反応拡散系（チューリング・パターンと同じ枠組み）。別系統の振動反応にブリッグス–ラウシャー反応（BR反応、$\ce{IO3^-}$＋過酸化水素＋マロン酸＋でんぷん、無色↔黄↔青黒を往復）がある。関連：[[auto-add-futoshita-gimon]]