第 8 章：BZ 振荡反应（BZ Oscillating Reactions）

8.1 Belousov 反应与 Field–Körös–Noyes（FKN）模型（Belousov Reaction and the Field–Körös–Noyes (FKN) Model）

Belousov–Zhabotinskii 反应（简称 BZ 反应）由俄罗斯生化学家 Boris Belousov 于 1951 年发现，最初因期刊编辑的轻蔑拒绝而未正式发表——当时的教条是化学反应不可能振荡。Belousov 最终在 1959 年的一个不出名的俄罗斯医学会议文集中发表了一则简短摘要。他所观察的是催化剂离子浓度比值的振荡：在原始 Belousov 反应中催化剂是铈离子，参与柠檬酸被溴酸根氧化的过程，振荡通过 Ce³⁺ 与 Ce⁴⁺ 之间的颜色变化显现；用亚铁邻菲咯啉（ferroin）作催化剂时颜色变化更显著——Fe²⁺ 态呈砖红色，Fe³⁺ 态呈亮蓝色。Zhabotinskii (1964) 继续研究该反应，后来合称 Belousov–Zhabotinskii 反应或简称 BZ 反应。当反应物同时可以扩散时，可以形成各种复杂的空间花样，正是这些空间花样维持了生物学家与物理学家的持久兴趣，并推动了关于生物图式与形态生成中时空自组织基本问题的研究。Belousov 的开创性工作直到 1980 年才以追授形式获得列宁奖。该反应是真正具有意义的振荡反应；Tyson (1994) 给出简洁而较新近的综述，Winfree (1984) 对其历史给出有趣记述。

值得特别指出的是该反应在数学上与若干真实生物振子的相似性。BZ 反应的近似方程与黏菌 Dictyostelium discoideum 自组织过程中细胞周期信号传导的现实模型相同（见 Goldbeter 1996）；这种类比已被 Winfree (1987) 创造性地用于建模心室中的三维电活动，被 Tyson (1991) 用于研究细胞周期的分子生物学。该反应绝非仅是学术好奇心驱动的玩具系统。

现在所谓 BZ 反应泛指一类化学振荡反应：有机物在酸性介质中被溴酸根氧化，由金属离子（典型为铈或亚铁邻菲咯啉）催化。原 Belousov (1951) 反应以铈为催化剂氧化丙二酸，铈在 Ce³⁺ 与 Ce⁴⁺ 两态间循环；实验观察到铈离子浓度的持续周期振荡。Field 等 (1972) 详细测量了溴离子 [Br⁻] 与铈离子浓度比 [Ce⁴⁺]/[Ce³⁺] 的时间变化。图 8.1 表明振荡可被定性理解为两段过程 I 与 II 的交替：[Br⁻] 高时（A 点附近）过程 I 主导，Br⁻ 沿 AB 段被消耗，铈主要为 Ce³⁺ 态；当 [Br⁻] 进一步下降并越过临界值 B 时快速跌至低值 C，过程 II 接管，Ce³⁺ 被氧化为 Ce⁴⁺；在过程 II 中 Ce⁴⁺ 与 Br⁻ 又逆向反应回到 Ce³⁺ 态，使 [Br⁻] 沿 C→D→A 回升到高值，过程 I 重新主导，如此往复。这种定性的"开关"描述并不足以证明系统会真的振荡——两段过程也可能仅走到某种共存的稳态；事实上确实存在某些参数范围使系统走至稳态而非振荡。沿 BC 与 DA 的快速变化是弛豫振子（relaxation oscillator）的典型特征——极限环的某些部分被快速扫过；这种结构为用奇异摄动方法解析估计周期提供了契机（第 8.4 节讨论）。

Field 和 Noyes (1974) 将整个反应机制合理压缩为 5 个关键反应并用质量作用律写出 3 种化学物的动力学方程，称为 FKN 模型（又因相关研究在俄勒冈大学完成而被称为 "Oregonator"）。关键物种为 HBrO₂（记为 X）、Br⁻（Y）、Ce⁴⁺（Z）、BrO₃⁻（A）与 HOBr（P），5 步反应的近似形式为：A + Y →(k₁) X + P；X + Y →(k₂) 2P；A + X →(k₃) 2X + 2Z；2X →(k₄) A + P；Z →(k₅) fY。其中前两步对应过程 I，后三步对应过程 II。Field 与 Noyes 实际使用的形式是 B + Z → fY + Q，B 代表有机酸，Q 是 CO₂；[B] 视为常数并入速率常数。把 [A] 视为常数（溴酸根充足供应）、忽略 [P] 的细节，由质量作用律即得无量纲前的三阶反应动力学方程组：

dx/dt = k₁ a y − k₂ x y + k₃ a x − k₄ x², dy/dt = −k₁ a y − k₂ x y + f k₅ z, dz/dt = 2 k₃ a x − k₅ z.

振荡行为强烈依赖于参数：若 k₅ = 0，y 单调衰减到零，无振荡；若 k₅ 很大、f = 0.5 时第三与第五步合并为 A + X → 2X + Y，系统退化为双分子两物种机制；Hanusse (1972) 证明双分子两物种系统不可能振荡，因此 (f, k₅) 平面上必存在一个不能振荡的参数域。

8.2 FKN 模型的线性稳定性分析与极限环解的存在（FKN Model and Limit Cycle Solutions）

分析系统 (8.3) 的合理方式是引入无量纲化。Murray (1977) 与 Tyson (1982, 1985) 给出不同的无量纲化方式与化学解释；这里采用 Tyson (1985) 的标度，使典型速率常数取到实验值。引入标度 x = x/x₀、y = y/y₀、z = z/z₀、t = t/t₀，其中 x₀ = k₃ a / k₄ ≈ 1.2 × 10⁻⁷ M、y₀ = k₃ a / k₂ ≈ 6 × 10⁻⁷ M、z₀ = 2(k₃ a)² / (k₄ k₅) ≈ 5 × 10⁻³ M、t₀ = 1/k₅ ≈ 50 s，并定义小参数 ε = k₅ / (k₃ a) ≈ 5 × 10⁻⁵、δ = k₄ k₅ / (k₂ k₃ a) ≈ 2 × 10⁻⁴、q = k₁ k₄ / (k₂ k₃) ≈ 8 × 10⁻⁴、f ≈ 0.5。代入 (8.3) 并去掉星号得无量纲方程组：

ε dx/dt = qy − xy + x(1 − x), δ dy/dt = −qy − xy + 2 f z, dz/dt = x − z.

三个小参数 ε、δ、q 是后续分析的关键所在。

分析系统 (8.3) 的标准做法是把它写成向量形式 dr/dt = F(r; ε, δ, q, f)，其三个分量依次是 ε⁻¹(qy − xy + x − x²)、δ⁻¹(−qy − xy + 2f z)、x − z。这样写的好处是把小参数 ε 与 δ 从时间导数中显式分离出来。线性稳定性分析的程序是标准的（详见第 3 章）：先求非负稳态或多个正稳态；接着确定线性化矩阵的特征值；最后寻找一个"约束集"——有限封闭曲面 S，其上 n · dr/dt < 0（n 是单位外法向），使 S 内 t₀ 时刻的解在所有 t > t₀ 永远留在 S 内。Murray (1974) 对原始（与 (8.6) 略有差异）的 FKN 方程做过此分析，本节沿用他的方法。

令 (8.5) 左端为零解代数方程组，得非负稳态。除 (0, 0, 0) 之外还存在正稳态 (xₛ, yₛ, zₛ)，其中 zₛ = xₛ、yₛ = 2f xₛ / (q + xₛ)、xₛ 由二次方程解出，表达式为 xₛ = (1/2){(1 − 2f − q) + [(1 − 2f − q)² + 4q(1 + 2f)]^(1/2)}，并带正号取根（保证 xₛ > 0）；另一非零稳态有负分量故舍去。

把 (8.5) 在 (0, 0, 0) 附近线性化得稳定性矩阵 A，行列式 |A − λI| = 0 展开成特征方程 λ³ + λ²(1 + qδ⁻¹ − ε⁻¹) − λ[ε⁻¹(1 + qδ⁻¹) − qδ⁻¹] − q(1 + 2f) / (εδ) = 0。考察此三次方程 λ ≥ 0 时的图像或利用三次方程根之积为 q(1 + 2f) / (εδ) > 0 这一事实，都可得出至少存在一个正实根的结论；故原点 (0, 0, 0) 始终线性不稳定——它不会成为 BZ 反应长时间行为的归宿。

把 (8.5) 在正稳态 (xₛ, yₛ, zₛ) 附近线性化，经一段冗长但初等的代数化简，特征方程可写成 λ³ + Aλ² + Bλ + C = 0，其中各项系数为：A = 1 + q + xₛ/δ + E/ε，B = (q + xₛ)/δ + E/ε + (q + xₛ)E/(εδ) + yₛ(q − xₛ)/(εδ)，C = (q + xₛ)E/(εδ) − 2f(q − xₛ)/(εδ) + yₛ(q − xₛ)/(εδ) = (xₛ² + q(2f + 1)) / (εδ)。这里引入了辅助量 E = 2xₛ + yₛ − 1，用 (8.7) 化简后 E = (xₛ² + q(xₛ + 2f)) / (q + xₛ) > 0。由 E > 0 推出 A > 0；由 C 的最简形式 (xₛ² + q(2f + 1)) / (εδ) > 0 推出 C > 0；B 的符号未定，可正可负。由 Descartes 符号法则（见附录 B）得知 λ³ + Aλ² + Bλ + C = 0 在 (8.8) 中至少有一个负实特征根。全部特征根实部为负的充要条件是 Routh–Hurwitz 条件（见附录 B）AB − C > 0。把 A、B、C 的具体表达式代入，AB − C 化简为 φ(δ, f, ε) = (N δ² + M δ + L) / δ² > 0，其中 N = [(xₛ² + q(xₛ + 2f))(1 + E/ε)] / [ε(q + xₛ)] > 0、L = (q + xₛ){(q + xₛ) + [xₛ(1 − q − 4f) + 2q(1 + 3f)] / ε}，M 是 f, q, ε, xₛ 的函数（M 在后续分析中不需，故 Murray 未给出）。故正稳态 (8.7) 线性不稳定的条件是参数 (δ, f, ε) 落在 φ(δ, f, ε) < 0 的区域内，分岔曲面即 φ(δ, f, ε) = 0。

当 |B| ≫ 1 的极端情形可读出特征值的渐近行为：若 B > 0 且很大，λ ∼ −C/B（实负）、−A/2 ± i√B（实部 −A/2 < 0），稳态线性稳定；若 B < 0 且绝对值很大，λ ∼ C/|B|（正实部）、±√|B|（含正实部），稳态线性不稳定。当参数使 B = C/A 时特征方程恰有根 λ = −A 与 ±i√B，这是分岔的临界位置；稍微偏离 B = (C/A) − ω（0 < ω ≪ 1）时，展开可知 λ = ±i(C/A)^(1/2) + O(ω) 的 O(ω) 项有正实部——增长振荡的失稳形式。Hopf 分岔定理（Strogatz 1994 给出）的条件在分岔面 φ = 0 的邻域内满足，故在 φ < 0 侧的小邻域内系统存在小振幅极限环，周期 T = 2π / (C/A)^(1/2)。不过 (8.4) 给出的实验参数下振幅并非小振幅，因此该周期公式主要是教学价值；真正实用的是分岔面 φ(δ, f, ε) = 0 的位置。

由 (8.10) 知 N δ² + M δ + L = 0 是关于 δ 的二次方程，N > 0。线性不稳定要求 δ 落在 (8.15) 给出的区间 0 < δ < (−M + √(M² − 4LN)) / (2N)。要求该区间非空则右端必须为正，这等价于 L < 0；用 (8.4) 中 ε 极小的性质对 L < 0 做一阶展开（即保留 ε⁻¹ 项令其为零，O(1) 项保留）得到 (1 − 4f − q){(1 − 2f − q) + [(1 − 2f − q)² + 4q(1 + 2f)]^(1/2)} + 4q(1 + 3f) < 0，这是 f 与 q 的代数方程。把 q = 8 × 10⁻⁴ 视作极小参数分两段求解：若 (1 − 2f − q) > 0 即 2f < 1 情形直接令 q → 0 得 ₁f_c ≈ 1/4；若 2f > 1 情形则需小心做 O(q) 展开（因为 q 出现在分母的 (1 − 2f − q) 中），令 (1 − 2f − q) 的绝对值倒数展开后得到 ₂f_c ≈ (1 + √2)/2。故正稳态线性不稳定的 f 范围为 ₁f_c ≈ 1/4 < f < ₂f_c ≈ (1 + √2)/2。最后，对每个固定的 ε，δ 关于 f 的分岔曲线由 (8.18) 即 (8.15) 的右端给出，其中 ₁f_c, ₂f_c 由 (8.16) 确定，L, M, N 由 (8.10) 与 (8.9) 通过 f, q, ε 表达。

Murray 在此强调这一分岔面是"实用"对象：实验测得参数 (8.4) 后直接代入 φ 公式即可判断正稳态是否失稳。这是连接抽象分析（特征方程、Hopf 定理）与具体化学体系（BZ 振子）的桥梁。整节的分析是为第 8.3 节的全局稳定性讨论和第 8.5 节的弛豫振子近似做数学准备——前者的"约束集"需要先知道正稳态的位置，后者的"零斜线"分岔形态直接由 (δ, f) 的取值决定。

8.3 FKN 模型的非局部稳定性（Nonlocal Stability of the FKN Model）

第 8.2 节证明了对每个 ε，若 δ 与 f 落在适当区域内，正稳态线性不稳定——且当 (δ, f) 接近分岔曲线时失稳形式为增长振荡。在不稳定的整个参数区域内还须考虑全局稳定性。三阶系统没有 Poincaré–Bendixson 型的定理（这与二阶相平面情形截然不同——二阶情形下 Poincaré–Bendixson 给出周期解存在的拓扑判据），但周期解的存在要求系统有一个约束集 S（有限封闭曲面，n · dr/dt < 0）。约束集加单一不稳定稳态（且失稳为增长振荡）虽不严格证明极限环存在，但给出足够理由继续推进。三阶系统是可能存在混沌的——BZ 反应中已观察到混沌行为（Scott 1991）；Lorenz (1963) 模型是三方程流体系统，展示混沌。Hastings 与 Murray (1975) 对 FKN 模型 (8.5) 给出存在极限环周期解的严格证明，并提出一个确定周期轨迹一般路径的程序；该程序对相当广泛的一类反馈控制系统也适用（Hastings 等 1977）。Murray 强调"rigourous proof" 是关键——这不是定性图像，而是带完整数学细节的存在性证明。

构造约束集的最简单形状是一个矩形盒 S，由 6 个面 x = x₁, x = x₂; y = y₁, y = y₂; z = z₁, z = z₂ 围成正稳态 (xₛ, yₛ, zₛ) 所在的区域。逐面检查 n · dr/dt < 0 的条件。设 i, j, k 是 x, y, z 方向的单位法向。取 x₁, x₂ 满足 0 < x₁ < xₛ < x₂。x = x₁ 面上 n = −i，(8.19) 要求 −i · dr/dt |{x=x₁} = −dx/dt | < 0；由 (8.5) 第一式展开并用 ε⁻¹ 系数去掉得 [qy − xy + x − x²]{x=x₁} > 0。注意到 0 < q ≪ 1，若设 x₁ = O(q)，则该不等式化为 y(q − x₁) + x₁ − x²，对实验 q ≪ 1、0 < f < 1，xₛ 严格落在 (0, 1) 内），q = x₁ < xₛ < 1 在 q < 1 时成立。实验参数下这些条件都满足。} ≈ y(q − x₁) + x₁ > 0（对所有 y₁ ≤ y ≤ y₂ 取该不等式的最坏情形）。当 y 落到下界 y₁ = O(q²) 附近时，最坏情形要求 x₁ > 0 即可；考虑到 x₁² 已是 O(q²) 阶可忽略，约束的"主导"项是 x₁ > 0。故 x 的下界可取 x₁ = q 作为一阶近似；在 x₁ = q 处 −i · dr/dt = −q(1 − q) / ε < 0（q < 1 时成立）。x = x₂ 面上 n = i，要求 [y(q − x) + x − x²]_{x=x₂} < 0；若取 x₂ = 1，得 i · dr/dt = ε⁻¹ y(q − 1) < 0（q < 1 且 y > 0 时成立）。由 (8.7) 给出的 xₛ 表达式做代数验证（xₛ 取 (1/2){(1 − 2f − q) + [(1 − 2f − q)² + 4q(1 + 2f)]^(1/2)

z 方向：在 z = z₁ 面（z₁ < zₛ）n = −k，要求 −k · dr/dt = −(x − z) |{z=z₁} < 0；由于 S 内 x ≥ x₁ = q，z 的自然下界取 z₁ = q（严格来说 z₁ 应略小于 q，因为 x₁ = q）。z = z₂ 面（z₂ > zₛ）n = k，要求 (x − z) | < 0；因 S 内 x ≤ 1，可取 z₂ = 1（严格来说 z₂ 应略大于 1）。这是 z 方向约束集的"几何"——在 S 内 z 必落在 [q, 1] 区间。这两个面给出的不等式都是单变量线性关系，所以是所有面里最容易处理的。

y 方向：在 y = y₁ 面（y₁ < yₛ）n = −j，要求 (8.19) 化为 [y(q + x) − 2f z]{y=y₁} < 0，即 y₁ < 2f z / (q + x) 对 S 内所有 q ≤ x ≤ 1、q ≤ z ≤ 1 成立。右端 2f z / (q + x) 的极小值出现在 x 取最大、z 取最小时（因为分母递增、分子与 z 成正比），等于 2f · z_min / (q + x_max) = 2f q / (q + 1)。故 y 的下界可取 y₁ = 2f q / (q + 1)。y = y₂ 面（y₂ > yₛ）n = j，要求 2f z − y(q + x) | < 0，即 y₂ > 2f z / (q + x) 对 S 内所有 q ≤ x ≤ 1、q ≤ z ≤ 1 成立。右端 2f z / (q + x) 的极大值出现在 x 取最小、z 取最大时，等于 2f · z_max / (q + x_min) = 2f · 1 / (q + q) = f / q。故 y 的上界可取 y₂ = f / q。用 (8.4) 的典型 q 与 f 值验证 y₁ < yₛ < y₂。综上矩形盒 S 的 6 个面为 x = q, x = 1; y = 2f q / (q + 1), y = f / q; z = q, z = 1；只要 f 与 q 满足一定不等式（实验参数下满足）正稳态 (8.7) 即被围于 S 内，且 S 的 6 个面上 n · dr/dt < 0 均成立，故 S 内的解永不出 S。这是极限环存在的关键定性证据。Hastings 与 Murray (1975) 的工作正是基于这种矩形盒约束集给出严格证明。图 8.2 把该模型在 f = 0.3, δ = 1/3, q = 5 × 10⁻³, ε = 0.01 下的数值极限环与 Field 等 (1972) 的实测振荡做了对比，吻合良好。Murray 强调"a rectangular box" 是最简的形状——可能不是最紧的，但足以给出存在性证明；若需要更紧的界可改用其他形状（例如分段曲面的圆台型）。

矩形盒 S 可被进一步细化以获得更紧的界——例如把面改为非平面或用分段函数拟合 f(x) 行为——但极限环最终只能由数值或下文 (8.4 节) 的渐近方法确定——这里只须展示约束集存在即可。Murray 强调存在性证明 + 数值验证两路并行，是研究复杂反应-反应扩散体系的范式：纯拓扑/解析论证给出"为什么"会有振荡，数值验证给出"具体形态"是什么。本节作为第 8.2 节的"补充"在数学上把局部不稳定与全局有界性合在一起，为第 8.5 节中"约束集存在 ⇒ 极限环存在"这一逻辑链补完最后一步。

8.4 弛豫振子：BZ 反应的近似（Relaxation Oscillators: Approximation for the Belousov–Zhabotinskii Reaction）

重新审视图 8.1 可见，极限环的某些段被迅速扫过。Br⁻ 的轨迹尤其明显：沿 DA 突然上升、沿 BC 突然下降。极限环的某些部分相对其他部分被快速扫过的振子称为弛豫振子；这在数学上要求微分方程组的关键位置存在小参数。考察最简单的弛豫振子模型 ε dx/dt = y − f(x)、dy/dt = −x（0 < ε ≪ 1），其中 f(x) 连续且 f(x) → ±∞ 当 x → ±∞。f(x) = (1/3)x³ − x 的情形是经典 Van der Pol 振子。系统 (8.21) 是典型奇异摄动问题：ε 乘以某一导数项。图 8.3(a) 描绘 (8.21) 的典型极限环相平面轨迹，相应的 x(t) 解如图 8.3(b)。

从 (8.21) 的第一式可知：除了 y ≈ f(x) 的邻域外，x 以 O(1/ε) 量级快速变化。参照图 8.3(a)，沿 DA 与 BC x(t) 变化快速；这两段在 x(t) 上对应图 8.3(b) 中急剧上升与下降的窄区间。在这些段上应取 τ = t/ε 作为独立变量，(8.21) 的第二式在 ε → 0 时变为 dy/dτ = −εx ⇒ y ≈ 常数——这正与图 8.3(a) 上 DA 与 BC 的水平段对应。在零斜线 y = f(x) 上的 AB 与 CD 段，第二式化为 f′(x) dx/dt ≈ −x，可积出 x 关于 t 的隐式表达。若 f(x) 取 Van der Pol 三次式（或者可以合理地分段线性近似），则可以精确积出。极限环周期 T 的 O(1) 估计主要来自 AB 与 CD 段所用时间，DA 与 BC 段所用时间为 o(1)。可证 T 的渐近修正项量级为 O(ε^(2/3))；因此 O(1) 近似已经足够。Murray 这里强调"asymptotic limit cycle period of (8.21) has a correction of O(ε^(2/3))"——这是 1970 年代以来奇异摄动理论的经典结论，本书不展开推导。

具体到 Van der Pol 振子 f(x) = (1/3)x³ − x，零斜线 y = f(x) 与图 8.3(a) 的极限环形状类似——只是原点被平移到 DB 中点（不是通常意义下的 (0,0)），x(t) 解关于 x = 0 轴对称（也对应图 8.3(b) 的水平轴）。从 A 至 B 与从 C 至 D 段积分 (8.22) 可得 A 在 (2, 2/3)、B 在 (1, −2/3)；C 与 D 由对称性位于 (−1, 2/3) 与 (−2, −2/3)。取 t = 0 在 A，由对称性得 O(1) 周期 T 满足 (1/2) ∫{B}^{A} dx/x = −∫ dx/x 在 f(x) = (1/3)x³ − x 情形下可以显式做出（因零点斜线 y = f(x) 在 (−2, −2/3) 至 (2, 2/3) 间单调递增且有简单代数关系），这是 T = 3 − 2 ln 2 这一漂亮解析结果能写出来的原因。}^{T/2} dt 给出 T = 3 − 2 ln 2。即使 f′(x)/x 不可直接积分，零斜线 y = f(x) 的代数关系也可直接给出极限环变量 x(t) 的极值（对应于 A 与 C）以及 y(t) 的极值（对应于 D/A 与 B/C）。Murray 强调这是 (8.22) 的关键观察：周期估计只需对零斜线做代数与积分，避开对快速段 DA、BC 的显式处理。具体的积分 ∫_{B}^{A

比较图 8.1 的 Br⁻ 时间曲线与图 8.3(b) 的时间解，可见 BZ 振子确实属于弛豫振子，可用类似技巧作分析极限分析。Murray 在这里做了关键的物理化学对应：图 8.1 中 Br⁻ 的 [Br⁻] 在 DA 段（即 t_A → t_D → t_A）突然升高至峰值，对应零斜线 y ≈ f(x) 上的"上分支"附近的一个快速水平运动；BC 段的快速下跌则对应"下分支"附近的水平运动。这种"水平段 + 慢段"的模式是所有弛豫振子的共同指纹，不仅 Van der Pol、BZ 振子有，第 7 章的神经动作电位发放、糖酵解振荡等生物振子也都呈现这一结构。下面沿用 Tyson (1976, 1977) 的方法。回到 FKN 机制的无量纲形式 (8.24)（与 (8.5) 同一系统），注意 ε ≪ δ；因此可在 ε dx/dt ≈ 0 的近似下降低方程阶数，得到代数关系 0 = qy − xy + x(1 − x) ⇒ x = x(y) = (1/2){(1 − y) + [(1 − y)² + 4qy]^(1/2)}。把 x(y) 代回 (8.24) 即得降阶后的二阶微分方程组：

δ dy/dt = 2 f z − y [x(y) + q], dz/dt = x(y) − z.

(8.26) 可以在 (y, z) 相平面上完整分析：求稳态、做线性稳定性分析、构造约束集，从而导出极限环解存在的参数条件。这正是下一节要做的事。降阶的合理性来自奇异摄动：ε ≪ 1 意味着 x 的响应时间比 y, z 快得多，于是在 y, z 的时间尺度上 x 已"瞬时"跟踪到由 (8.25) 给出的准稳态 x(y)。Murray 用 1984 年自己的教材《Asymptotic Analysis》做参考，这是他常用的摄动分析范式。重要的是 (8.26) 的"order reduction"并不是无代价：原系统的相空间是 3 维，降阶后变成 2 维；我们丢掉了 ε 方向上的快速变化信息（DA、BC 段），但保留了关键的"缓慢"动力学（AB、CD 段），后者正好对应 BZ 振子中 [Br⁻] 与 [Ce⁴⁺] 缓慢演化的"化学平台期"。Murray 这里用词是"effectively collapse"——第三与第五步反应（(8.2) 中 A + X → 2X + 2Z 与 Z → fY）在 k₅ 很大时合并为单步 A + X → 2X + Y；这与把 ε dx/dt 视为 0 是同一件事的两面。

8.5 BZ 反应弛豫模型极限环振荡分析（Analysis of a Relaxation Model for Limit Cycle Oscillations in the Belousov–Zhabotinskii Reaction）

由 (8.26) 的 z 零斜线 z = x(y)，在 0 < q ≪ 1 下有渐近展开：x(y) ≈ (1 − y) 对 q ≪ 1 − y ≤ 1，x(y) ≈ qy / (y − 1) 对 q ≪ y − 1。y 零斜线 z = y[x(y) + q] / 2 在三段范围内有近似：当 q ≪ 1 − y ≪ 1 时 z ≈ y(1 − y) / (2f)；当 q ≪ y − 1 时 z ≈ y[q y / (y − 1) + q] / (2f)；当 y ≫ 1 时 z ≈ qy / f。z 零斜线关于 y 单调递减。y 零斜线在 f 取不同值时形态不同：局部极大点为 z_max = z_D = z_A = 1/(8f)，y_max = 1/2；局部极小点为 z_min = z_B = z_C，在 y_min = (2 + √2)/2 处取到，z_min = q(3 + 2√2) / (2f)。

图 8.4(a)–(c) 是 0 < q ≪ 1 下 (8.26) 约化 BZ 模型的三种零斜线位形，由 (8.27)、(8.28)、(8.31) 联立得到 f 范围。对应关系为：f < 1/4 对应图 8.4(b)（局部极大在稳态右侧）；1/4 < f < (1 + √2)/2 对应图 8.4(a)；f > (1 + √2)/2 对应图 8.4(c)。第 7 章已指出图 8.4(a) 是可能存在极限环振荡的位形；故 (8.33) 给出的 f 范围内、稳态不稳定的条件下约化 BZ 模型 (8.26) 有极限环解。

把图 8.4(a) 与图 8.3(a) 的弛豫极限环振子作类比，可写出 (8.26) 在 0 < δ ≪ 1 下的 O(1) 周期 T ≈ ∫AB dt + ∫_CD dt = ∫ [dz/dt]⁻¹ dz + ∫}^{z_A{z_C}^{z_D} [dz/dt]⁻¹ dz。具体地 T_AB = ∫ [x(y) − z]⁻¹ dz、T_CD = ∫}^{z_A{z_C}^{z_D} [x(y) − z]⁻¹ dz；为便于计算把 z 对 y 的函数代入。沿 AB 段的大部分有 x(y) ≈ qy / (y − 1) ∼ q（y 大时），由此得 T_AB = ∫ (q − z)⁻¹ dz = ln[(z_A − q) / (z_B − q)] ∼ −ln[4(3 − 2f + 2√2) q]（q ≪ 1）。注意 x(y) 在 y → ∞ 时才渐进为 q，所以这是一个 T_AB 的上界，下界由在 x(y) 中将 q 替换为 q(1 + √2) 得到，即 ln[4(3 − 2f + 2√2) q] < T_AB < −ln[4(3 − 2f + 2√2)(1 + √2) q]。}^{z_A

CD 段则用 x(y) ≈ 1 − y、z ≈ y(1 − y) / (2f) 代入；改换积分变量为 y，并利用 dz = [(1 − 2y) / (2f)] dy，整理得 T_CD ≈ [(4f − 1) / (2f − 1)] ln[2^(1/(4f−1)) · (4f − 1) / (4f)]。周期 T = T_AB + T_CD（量纲乘以 t₀ 还原）。至于弛豫振子周期内变量的时程：取 t = t_A = 0、z(0) = z_A，AB 段由 dz/dt ≈ q − z 给出 z(t) = z_A e⁻ᵗ + q(1 − e⁻ᵗ)，z 在 O(1) 时间内从 z_A 指数衰减到 z_B = O(q)；BC 段所用时间可忽略。CD 段有 ln[(1 − y) / ((2f − y)^(4f−1))] = K + (2f − 1) t，K 由 y(t_C) = y_C 给出，其中 t_C ≈ t_B 由 (8.42) 给出 t_B = ln[z_A (1 − q) / (z_B − q)]。由此可解出 y(t)，再由 z ≈ y(1 − y) / (2f) 还原 z(t)；DA 段时间同样可忽略。

回到真实 BZ 反应：z 与 y 分别是 Ce⁴⁺ 与 Br⁻ 的无量纲浓度。从图 8.4(a) 的 A 出发，极限环轨迹经 BCDA 回到 A，对应图 8.1 实验振荡的相应字母：[Br⁻] 从 A 指数下降到 B（达到临界 Br⁻ 浓度即第 8.1 节描述的过程 I→II 切换阈值），然后快速跌至 C；C 至 D 是 O(1) 时间的 [Br⁻] 增大，由 (8.42)、(8.43) 刻画；接着快速升回 A 处的 [Br⁻]。弛豫振子近似 (8.26) 模拟了真实实验振荡 (8.1)。与实测量的具体对比：Field 与 Noyes (1974) 实验给出 f = 0.5；在 f → 0.5 极限下 T_CD = 2 ln 2 − 1，由于 q 很小周期的主要贡献来自 T_AB。代入 (8.4) 的参数并用 (8.31) 的极限值换算 [Br⁻]，表 8.1 列出弛豫振子近似所得与 Field 等 (1972) 实测值的对比。Period 近似为 183–228 s、实测约 110 s；[Br⁻]_B 计算 1.7 × 10⁻⁵ [BrO₃⁻]、实测 2 × 10⁻⁵ [BrO₃⁻]；[Br⁻]_C 计算 0.3 [Br⁻]_crit、实测 0.3 [Br⁻]_crit；[Br⁻]_A 计算 1.6 × 10⁻³ [BrO₃⁻] = 90 [Br⁻]_crit、实测 3 [Br⁻]_crit。考虑到反应本身的复杂度与所引入的多层近似，理论与实验吻合程度相当好。

当 f 落在图 8.4(b) 与 8.4(c) 范围时，机制 (8.26) 不产生振荡（见第 7 章第 7.3 节），但可能出现阈值行为（图 7.5(a)(b)）。若允许 FKN 基础机制 (8.2) 的反应可逆，则可能出现三个正稳态其中两个稳定（图 7.4(d)），对应第 7.3 节讨论的生物开关行为（图 7.5(c)）；可逆模型的数值研究还显示 bursting 行为与混沌。若干 BZ 反应的可能模型预测了多种反常现象，许多后来都被实验证实（Tyson 1985、Field 与 Burger 1985、Scott 1991）：Barkley 等 (1987) 证明存在包括周期 bursting、滞回、周期-混沌序列的复杂行为；Györgyi 与 Field (1991) 给出 BZ 反应中确定性混沌的模型。本章只讨论了均匀或良好搅拌系统；后续章节讨论耦合生物振子与无搅拌系统（须含扩散）时会出现一系列新的惊人现象；BZ 反应又一次成为验证理论预测的关键实验体系。

本章个人批注

读 8.1 时最冲击我的是 Belousov 工作的"被拒稿"经历——一个重要发现因为违反当时教条而几乎被埋没。Murray 在这里不厌其详地叙述这段史话（Field 与 Burger 1985 编辑的书中也有翻译原文），是有用意的：科学共同体对新现象的开放度是一个变量，振荡反应从"不可能"到"原型反应"用了 30 年。Winfree 1984 是另一个重要历史记述，值得找来看看。

8.2 节关于无量纲化的讨论是一处隐藏的精华。Murray 明确说"通常有几种无量纲化方式"，并采用 Tyson (1985) 的标度——这个选择的目的是让小参数（ε, δ, q）显现出来，从而为后续奇异摄动分析铺路。读 ODE 教材时这一点常被忽略：不同的标度选择会暴露或掩盖不同的数学结构。这让我想到自己做力学建模时也常常在无量纲化这一步反复权衡。

8.3 节的矩形盒约束集是本章最优雅的"工程"技巧。Hastings 与 Murray (1975) 给出的不是抽象的拓扑论证，而是显式构造——这给读者一个可以具体验证的对象。y 面的处理尤其精彩：通过在 (x, z) 的整个矩形上对 2f z / (q + x) 取极值，得出 y₁ = 2f q / (q + 1)、y₂ = f / q。这种"在最坏情形下也成立"的边界估计是常微分方程分析中的标准技巧，但在具体三阶系统上手工完成并不常见。

8.4–8.5 节是全章的高潮。Murray 把 Van der Pol 振子作为教学例子先把弛豫振子概念讲清，再把同样的思路套到 BZ 上；这种"先简单后复杂"的写法让我很受用。T_AB 与 T_CD 的积分细节相当繁琐，但 Murray 反复强调"all we want here is a reasonable approximation"——他不假装能给出精确的周期，而是给出量级正确、参数依赖合理的解析估计。表 8.1 的对比把这种"参数依赖合理"实证地呈现出来，理论与实验相差在 2–30 倍量级，但 30 年后我们仍能引用它作为 FKN 模型在弛豫极限下有效性的证据。

最后我注意到一个连贯性线索：第 7 章的极限环分析工具（Hopf 分岔、约束集）在第 8 章被原封不动地搬到了三阶的 FKN 系统上。Murray 的写法让两章之间有清晰的"概念复用"关系——这是教科书最好的样子。

与上下章的衔接（一段话）

第 7 章是"原型振子概念引入"，从反馈控制、Goodwin 模型、神经/激素的振荡讲到二维相平面上的极限环、阈值与开关；它为读者装备了"分岔 + 极限环 + 约束集"这一套分析工具，并明确指出 BZ 反应是研究反应动力学振子的关键实验体系，但那时还未深入到具体反应机制。第 8 章紧随其后，把这套分析工具搬到了具体的 BZ 反应：8.1 节引入 Belousov–Zhabotinskii 反应的历史与定性图像，并建立 FKN 三物种模型；8.2 节做标准的线性稳定性分析，找到分岔参数域 (f, δ)；8.3 节用矩形盒给出约束集并引用 Hastings–Murray 的极限环存在性证明；8.4–8.5 节进一步用 ε ≪ 1 的奇异摄动把模型降阶为二阶弛豫振子，进而解析估计周期、极值并与 Field 等的实验对比。读完后获得的不仅是一个具体反应的分析范例，更是"如何用前章工具拆解一个真实系统"的完整示范。第 9 章将把视线从单一振子（化学反应器）转到耦合振子与扩散系统——再次回到"空间结构"这个贯穿全书的母题。