跳转至

第6章：行波在兴奋系统中的传播

作者

本章作者与 Ch1-Ch5 相同——James P. Keener 与 James Sneyd。本章是 I 卷"细胞生理学"中从单细胞时空动力学到组织水平的关键章，角色是把 Ch5 的兴奋性单细胞 ODE 与 Ch4 的电缆方程结合，揭示动作电位沿轴突传播的物理机制。本章是数学物理学经典——行波、奇异性摄动、eikonal-curvature 方程、螺旋波——这些不仅是生理学工具，也是应用数学的不朽成果。具体内容：（1）前缘（traveling front）与脉冲（traveling pulse）的分类与求解（§6.1-§6.3）；（2）双稳态方程的精确解（§6.2）——奇异性摄动 + shooting 方法；（3）髓鞘化神经的离散模型（§6.2.2-§6.2.3）——跳跃式传播与多发性硬化；（4）HH 与 FHN 的行波（§6.3.1-§6.3.2）——理论预测与实验的精确吻合；（5）色散曲线与运动学理论（§6.4）——周期波列的传播速度-周期关系；（6）高维行波、eikonal-curvature 方程、螺旋波、卷波（§6.5）——心肌纤颤、视网膜幻觉的物理基础。

内容概述

本章围绕"兴奋系统如何产生行波"这一核心问题展开。

第一部分（§6.1-§6.2）双稳态方程与前缘。两类行波：traveling front（"拉链"式——把介质从静息态"拉"到激发态）和 traveling pulse（"脉冲"式——激发后回到静息）。双稳态方程（6.1）\(\partial V/\partial t = \partial^2 V/\partial x^2 + f(V)\)，其中 \(f\) 有三个零点 \(0, \alpha, 1\)（\(0 < \alpha < 1\)，\(V = 0, 1\) 稳定）。行波 ODE（6.5）\(U_{\xi\xi} - c U_\xi + f(U) = 0\)，转化为 2 维 \((U, W)\) 系统（6.6）–(6.7)。核心定理（§6.2.1）：shooting 方法证明唯一速度 \(c\) 存在（介于 \(c=0\) 过小和 \(c\) 过大之间），速度符号由 (6.8) 的面积 \(\int_0^1 f(u) du\) 决定——\(f\) 的"非对称性"决定前缘方向。显式解（6.17）–(6.19)：cubic 情形下解为 \(\tanh\) 曲线；piecewise-linear 情形（6.12）给出 \(c = (1-2\alpha)/\sqrt{\alpha - \alpha^2}\) 显式速度。

§6.2.1 末段"还原为物理参数"：(6.21) 给出轴突上的物理速度 \(s = (c/2)\sqrt{d/(R_m R_c)}\)——与实验 \(s = 21.2\) mm/ms 吻合（鱿鱼 \(d = 500\) μm）。比例律（6.22）\(s = \sqrt{d/10^{-6} m}\) mm/ms——直觉上"轴突越粗传导越快"的精确版本。

§6.2.2 髓鞘化：髓鞘把膜电阻提高 100 倍、膜电容降低 100 倍，节点间跳跃传导（saltatory conduction），速度提升因子 \(\sqrt{100 d/\mu} \approx 10\)（节点间距 \(L \approx 100 d\)，节点宽度 \(\mu \approx 1\) μm）。疾病学意义：多发性硬化的脱髓鞘使传导失败。

§6.2.3 离散双稳态方程：(6.27) 离散化。关键定理（§6.2.3 末段）：存在耦合阈值 \(D^* = \alpha(1-\alpha)/(2\alpha-1)^2\)——\(D \le D^*\) 时传播失败。这是连续方程没有的特性。

第二部分（§6.3）行波脉冲。核心区别：front 是 heteroclinic 轨道（连接两个不同稳态），pulse 是 homoclinic 轨道（从同一稳态出发再返回）。FHN 6.48）–(6.49) 的脉冲解：奇异性摄动分析（Keener 1980a）——\(\epsilon \ll 1\) 时，\(v\) 是快变量，\(w\) 是慢变量。外方程（6.67）–(6.68）\(dw/dt = G^\pm(w)\)——\(v\) 沿 cubic 零线的稳定分支滑动；内方程（6.71）–(6.72）——过渡层（wave front / wave back）。几何（Fig. 6.7）：脉冲由一个 upjump（v 从 \(V^-\) 跳到 \(V^+\)）+ 沿 \(V^+\) 分支的慢动力学 + 一个 downjump 组成。速度条件（6.74）–(6.75)：upjump 速度 \(c(w^+) > 0\) + downjump 速度 \(-c(w^-)\)，要求 \(c(w^+) = -c(w^-)\)。激发态持续时间（6.76）\(T_e = \int_{w^+}^{w^-} dw/G^+(w)\)。

§6.3.2 HH 行波脉冲（6.80）–(6.81）：shooting 方法从静息态的 4 维不稳定流形出发。HH 数值解：\(c = 21\) mm/ms（vs 实验 21.2 mm/ms——0.5% 误差）。这是 Hodgkin-Huxley 1952 最著名的"理论预测"。

第三部分（§6.4）周期波列。色散曲线（dispersion curve）：传播速度 vs 周期的关系。典型形状（Fig. 6.8-6.9）：快分支（典型动作电位，stable）+ 慢分支（小振幅，unstable），两支在 kneepoint 相遇。生物意义：神经传导是快分支上的稳定波列。运动学理论（§6.4.3）：不追踪每个动作电位的细节，只跟踪"第 \(n\) 个动作电位到达位置 \(x\) 的时刻 \(t_n(x)\)"，由 \(dt_n/dx = 1/c(t_n - t_{n-1})\) 描述（6.90）。核心守恒律（6.95）：\(t'_r(c_{n+1}) dc_{n+1}/dx = 1/c_{n+1} - 1/c_n - A'(c_n) dc_n/dx\)——动作电位持续时间 (APD) 与前速度的功能关系是核心。

第四部分（§6.5）高维行波。eikonal-curvature 方程（6.112）\(\tau R_t \cdot n = c_0 \sqrt{1/2} - \sqrt{2}\kappa\)——曲率 \(\kappa\) 减慢波前速度。物理直觉：弯曲的波前要把能量"散开"到更大弧长 → 速度较低。实验验证（Foerster et al. 1988）：BZ 反应螺旋波的曲率-速度关系验证了 eikonal-curvature 公式。

§6.5.2 螺旋波：(6.113)-(6.125) 给出螺旋波的 eikonal + 奇异性摄动解。关键关系（6.128）\(\omega = c^2 m^*/\epsilon\)，\(m^* = 0.330958\)——频率与速度的平方成正比。这是数值 + 解析的一致性检验。生物医学意义（§6.5.2 开头）：螺旋波在心脏 = 心律失常 / 纤颤；在皮层 = 癫痫发作；在视网膜 = 幻觉。

§6.5.2 末段 3D 推广：卷波（scroll wave，Winfree 1973）——螺旋波在 3D 的推广，包括闭合卷环、纽结卷、链状卷对。心律失常的 3D 机制是现代心脏电生理的核心。

主要数学结果： 1. 双稳态方程的 shooting 方法（6.5）-（6.10）+ 唯一性 2. 显式行波解（6.17）cubic 情形和（6.12）piecewise-linear 情形 3. 物理速度公式（6.21）+ 经验律 \(s = \sqrt{d/10^{-6}\text{m}}\) mm/ms 4. 离散双稳态的耦合阈值 \(D^* = \alpha(1-\alpha)/(2\alpha-1)^2\) 5. FHN 脉冲的奇异性摄动分析：外/内方程 + 速度匹配 6. HH 数值解 \(c = 21\) mm/ms（0.5% 误差） 7. 色散曲线：快分支 + 慢分支 + kneepoint 8. 运动学理论（6.95）：APD 与前速度的功能关系 9. eikonal-curvature 方程（6.112）+ 螺旋波频率-速度关系 \(\omega = c^2 m^*/\epsilon\) 10. 卷波（3D 螺旋波）的拓扑分类

前置知识：ODE（相平面、稳定流形、Hopf 分岔）、PDE（双稳态方程、奇异性摄动）、Ch5 兴奋性 + Ch4 电缆方程。

核心方程与概念

1. 双稳态方程与前缘

双稳态方程（6.1）：\(\partial V/\partial t = \partial^2 V/\partial x^2 + f(V)\)，\(f\) 有 3 零点 \(0, \alpha, 1\)（\(0 < \alpha < 1\)），\(f(0) = f(1) = 0\) 稳定。

行波 ODE（6.5）：\(U_{\xi\xi} - c U_\xi + f(U) = 0\)，\(\xi = x + ct\)。这是 2 维动力系统（6.6）–(6.7）。

面积定理（§6.2.1）：(6.8) \(c \cdot \int W^2 d\xi = \int_0^1 f(u) du\)——\(c\) 的符号由 \(f\) 在 \([0,1]\) 上面积决定。

Shooting 唯一性（§6.2.1）：\(c=0\) 时轨迹过冲不到 \(U=1\)；\(c\) 大时轨迹过冲过头；连续性 + 唯一性保证某个 \(c_0\) 让 heteroclinic 连接存在。

显式解（6.17）：

\[ U(\xi) = \frac{u_0 + u_2}{2} + \frac{u_2 - u_0}{2} \tanh\left(\frac{A}{\sqrt{2}} \frac{u_2 - u_0}{2} \xi\right) \]

物理速度（6.21）：

\[ s = \frac{c \lambda_m}{\tau_m} = \frac{c}{2} \sqrt{\frac{d}{R_m R_c}} \]

经验律（6.22）\(s = \sqrt{d / 10^{-6} \text{ m}}\) mm/ms——轴突直径决定速度。

2. 髓鞘化与离散传播

髓鞘化效果（§6.2.2）：\(R_m\) 提高 100 倍，\(C_m\) 降低 100 倍——\(\tau_m\) 降低 10⁴ 倍，\(\lambda_m\) 提高 10 倍。跳跃传导（saltatory）：动作电位只在节点发生，节点间是纯电阻。

离散方程（6.27）：

\[ \frac{dV_n}{d\tau} = f(V_n) + D(V_{n+1} - 2V_n + V_{n-1}) \]

耦合阈值（6.47）：\(D^* = \alpha(1-\alpha)/(2\alpha-1)^2\)——\(D \le D^*\) 传播失败。这是连续方程没有的特性，对多发性硬化的脱髓鞘导致传导阻滞至关重要。

3. 行波脉冲的奇异性摄动分析

FHN 方程（6.48）–(6.49）：

\[ \epsilon \partial_t v = \epsilon^2 \partial_x^2 v + f(v, w), \quad \partial_t w = g(v, w) \]

\(\epsilon\) 物理意义：快（v，膜电位）vs 慢（w，恢复变量）的时间尺度比。

外方程（\(\epsilon = 0\)）（6.67）–(6.68）：

\[ f(v, w) = 0 \Rightarrow v = V^\pm(w), \quad \frac{\partial w}{\partial t} = G^\pm(w) = g(V^\pm(w), w) \]

内方程（边界层）（6.71）–(6.72）：

\[ v_{\xi\xi} + y'(\tau) v_\xi + f(v, w) = 0, \quad y'(\tau) w_\xi = 0 \]

关键：(6.71) 是双稳态方程——对固定 \(w\)，有 \(c(w)\) 速度。脉冲由 upjump + 沿 \(V^+\) 分支 + downjump 组成。

速度条件（§6.3.1 末段）：\(c(w^+) = c(w^-)\)，其中 \(w^+\) 是 upjump 处的恢复变量值，\(w^-\) 是 downjump 处的值。

激发态持续时间（6.76）：\(T_e = \int_{w^+}^{w^-} dw/G^+(w)\)。

4. HH 行波脉冲

4 维 ODE 系统（6.80）–(6.81)：\(\xi = x/c + t\) 坐标系。Shooting 方法从 4 维不稳定流形出发。

关键结果（§6.3.2 末段）：\(c = 3.24 \lambda_m\) ms⁻¹ = 21 mm/ms（vs 实验 21.2 mm/ms）——HH 1952 论文最著名的理论预测。

通道密度的作用（Table 6.1）：节点 of Ranvier 2100 channels/μm² 是无髓鞘的 10 倍——跳跃传导的几何补偿。

5. 周期波列与色散曲线

色散关系 \(c = C(T)\)：传播速度 vs 周期。典型形状（Fig. 6.8-6.9）： - 快分支：典型动作电位，stable - 慢分支：小振幅，unstable - Kneepoint：绝对不应期

运动学理论（6.90）：

\[ \frac{dt_n}{dx} = \frac{1}{c(t_n(x) - t_{n-1}(x))} \]

——不需要追踪每个动作电位的细节。

守恒律（6.95）：

\[ t'_r(c_{n+1}) \frac{dc_{n+1}}{dx} = \frac{1}{c_{n+1}} - \frac{1}{c_n} - A'(c_n) \frac{dc_n}{dx} \]

——APD 与前速度的功能关系是核心。这是 Ch12 心肌组织动作电位传播的工程化基础。

6. 高维行波：eikonal-curvature 方程

eikonal-curvature 方程（6.112）：

\[ \tau R_t \cdot n = c_0 \sqrt{\frac{1}{2}} - \sqrt{2} \kappa \]

物理意义： - \(c_0 \sqrt{1/2}\)：平面波速度 - \(-\sqrt{2} \kappa\)：曲率减慢（负号因为弯曲的波前扩散能量）

实验验证（Foerster et al. 1988）：BZ 反应螺旋波的曲率-速度关系精确吻合 (6.112)。

关键应用（§6.5.2 开头 + Ch12）：心律失常/纤颤的物理基础——螺旋波在心肌组织中旋转。比正常心率快 10-20 倍的螺旋波就是室颤。

7. 螺旋波与卷波

螺旋波（Fig. 6.12-6.13）：旋转的自持续波。关键关系（6.128）：

\[ \omega = \frac{c^2 m^*}{\epsilon}, \quad m^* = 0.330958 \]

——频率与速度的平方成正比。这是数值与解析的桥梁。

卷波（3D 螺旋波，Winfree 1973）：包括闭合卷环、纽结卷、链状卷对。生物医学意义：心律失常的 3D 机制（Ch12 的核心主题）。

关键结论

  1. 行波 = 兴奋性 + 空间耦合（§6.1）：当兴奋系统的局部动力学与扩散/电缆方程结合，自然涌现行波。这是神经传导、心律、钙波、BZ 反应的共同数学基础。

  2. 行波速度由"面积"决定（(6.8) 面积定理）：\(c \propto \int_0^1 f(u) du\)——\(f\) 的非对称性编码了传播方向。这给出了为什么兴奋系统能定向传播的简洁解释。

  3. HH 1952 预测的 \(c = 21\) mm/ms 是 0.5% 误差（§6.3.2 末段）：从电压钳数据 + 4 维 ODE + 行波 shooting → 精确传导速度预测。这是生物数学建模的金标准。

  4. 髓鞘化使速度提升 10 倍（§6.2.2）：\(R_m\) 提高 100 + \(C_m\) 降低 100 → 跳跃传导。这是为什么脊椎动物有髓鞘神经的物理基础。

  5. 离散系统有传播失败阈值（(6.47) 耦合阈值 \(D^*\)）：连续方程没有的特性。这是脱髓鞘疾病（多发性硬化）导致神经传导阻滞的数学基础。

  6. 运动学理论 (6.95) 把行波简化为 1 维 ODE（§6.4.3）：不追踪每个动作电位细节，只跟踪当前 APD 与前速度。这是心律组织模拟的工程基础（Ch12）。

  7. eikonal-curvature 方程 (6.112) 把曲率纳入速度（§6.5.1）：弯曲的波前速度减慢。这是螺旋波的核心物理——曲率让"打转"成为可能。

  8. 螺旋波是心律失常的元凶（§6.5.2 开头）：心肌组织的螺旋波导致室颤。现代除颤理论（defibrillation theory, Keener 1990s-2000s）就是用强电场"破坏"螺旋波——这是临床电生理学的数学基础。

  9. 3D 卷波是螺旋波的拓扑推广（§6.5.2 末段）：闭合卷环、纽结卷、链状卷对。这是 Ch12 心律失常 3D 分析的前沿。

  10. \(\omega = c^2 m^*/\epsilon\) 关系（6.128）：螺旋波频率与速度的平方成正比。这是螺旋波"自持续"的几何必然——相速度必须与曲率匹配。

挑战和开放性问题

  1. FHN 速度匹配 \(c(w^+) = -c(w^-)\) 的存在性（§6.3.1 末段）：当 \(w\) 区间内 \(c(w)\) 单调时，方程 \(c(w^-) = -c(w^+)\) 总有解。当 \(c(w)\) 非单调时（\(w\) 区间内有 \(c(w) = 0\)），可能出现phase wave——下行跳变由外动力学决定。这一边界的严格判据是开放问题。

  2. 离散系统传播失败的严格判据（(6.47) 仅对 piecewise-linear 给出）：对一般 \(f\)，\(D^*\) 存在但显式公式依赖 \(f\)。这是 Keener 1987 给出的定性存在性——精确公式需要数值。

  3. 多维螺旋波严格存在性（§6.5.2 末段）：奇异性摄动解是渐近有效的但不是严格证明。Pelce-Sun 1991、Keener 1992、Kessler-Kupferman 1996 部分解决了单/双扩散变量情形，但完全严格仍是开放问题。

  4. 螺旋波稳定性（§6.5.2 末段）：哪些螺旋波是稳定的？Winfree 1991、Jahnke-Winfree 1991、Barkley 1994、Karma 1993/1994 都给出部分答案——但完整稳定性理论没有。

  5. 高维各向异性（Exercise 22）：\(D\) 是矩阵时 eikonal-curvature 方程修正为 \(S_t = c_0 \sqrt{\nabla S \cdot D \nabla S}/(Dk) + \nabla \cdot [D\nabla S / \sqrt{\nabla S \cdot D\nabla S}]\)——心肌纤维的各向异性（纤维方向 vs 横向）下螺旋波形态如何变？这是心律失常建模的核心。

  6. 除颤理论（§6.5.2 + 隐含）：强电场施加 → 螺旋波"破坏"。Pumir-Plonsey 1990s、Keener 2000s 用 eikonal-curvature 方程研究除颤阈值——临床实践 vs 理论的定量吻合仍是开放。

  7. 运动学理论的"反问题"（§6.4.3 末段）：给定 \(c = C(T)\) 色散曲线 → 推断 \(A(c), t_r(c)\) 函数。真实心脏组织的 \(C(T)\) 由实验测量（电生理标测）——如何反推组织异质性？这是现代心律失常诊断的核心。

  8. 卷波的拓扑与动力学（§6.5.2 末段）：闭合卷环、纽结卷、链状卷对的几何分类（Winfree-Strogatz 1983-1984）很完善，但这些拓扑结构的稳定性和它们与心律失常的对应关系仍不清楚。

  9. 细胞异质性的影响（§6.5 开头 + Ch12）：真实组织中 \(D, k, f, g\) 都是空间变化的。异质性对行波稳定性的影响——简单回答是"破坏螺旋波"，但定量没有。

  10. 二维 vs 三维（§6.5.2 末段）：2D 螺旋波 + 3D 卷波的对应关系。真实心脏是 3D 的——2D 模型只是近似。卷波环的自持续 vs 螺旋波的自持续是否有根本差异？

个人反思与批判性分析

建模哲学层面：

Ch6 是 I 卷"建模哲学的应用章"——Ch5 建立的兴奋性单细胞模型，与 Ch4 的电缆方程结合，涌现出行波这一新层次的现象。这是"局部 → 全局"涌现的经典例子——单细胞时间尺度分离（Ch5）→ 行波空间结构（Ch6）。

值得专门提出的是：Shooting 方法（§6.2.1, §6.3.2）是"局部 → 全局"涌现的工程技巧——从静息态的 4 维不稳定流形出发（Ch5 的局部分析），用 ODE 积分找到全空间的行波解（Ch6 的全局分析）。Hastings 1975、Carpenter 1977 用此方法严格证明 HH 行波的存在性——这是生物数学中少见的严格存在性证明。

eikonal-curvature 方程（§6.5.1）是几何光学的物理类比——在曲面上的波前传播用 Hamilton-Jacobi 方程描述，曲率项是几何约束的"代价"。这把行波物理与经典力学（最小作用原理）联系起来——这是物理学统一性的范例。

与后续章节的衔接：

  • Ch7（钙动力学）：钙火花、钙波是行波物理在亚细胞尺度的应用——eikonal-curvature 方程在 Xenopus 卵母细胞中的验证（Lechleiter 1991, Sneyd-Atri 1993）。
  • Ch8（细胞间通讯）：gap junction 介导的细胞间耦合——异质扩散 (6.96) 的应用。
  • Ch9（神经内分泌）：簇发放的"内行波"——单个细胞的钙/电压联合动力学。
  • Ch10（细胞功能调控）：基因调控网络的双稳态 + 行波（如胚胎发育中的形态发生波）。
  • Ch12（心脏电生理）：心肌组织中的螺旋波 + 卷波 + 除颤理论——Ch6 的直接应用。

值得重新推导的方程：

  • (6.5) 行波 ODE：从 (6.1) 出发，代入 \(\xi = x + ct\) → 一阶 + 二阶 ODE 转化为 2 维动力系统 (6.6)-(6.7)。
  • (6.8) 面积定理：从 (6.5) 出发乘以 \(U_\xi\) 积分，利用 \(U_{\xi\xi} U_\xi = d(W^2/2)/d\xi\) 的技巧。
  • (6.16) 显式速度：从 (6.13) cubic 出发猜测 \(W = -B(U-u_0)(U-u_2)\) 形式，代入 (6.5) 比较 \(U\) 的系数。
  • (6.21) 物理速度：从 (6.20) 无量纲到有量纲的速度变换——\(\lambda_m, \tau_m\) 的几何意义。
  • (6.71) 内方程：从 (6.69) 出发做空间重标度 \(\xi = (x - y(t))/\epsilon\)——\(\epsilon \to 0\) 时退化为双稳态方程。
  • (6.76) 激发态持续时间：从 (6.67) 外方程沿 \(V^+\) 分支积分。
  • (6.112) eikonal-curvature 方程：从 (6.107) 出发识别为主 + 比较平面波 (6.111)。
  • (6.128) 螺旋波 \(\omega = c^2 m^*/\epsilon\) 关系：从 (6.125) 出发做尺度分析——\(r \to \alpha r\) 等价于 \(c \to c/\alpha, \omega \to \omega/\alpha^2\)。

对作者建模选择的评价：

  • 优点：（1）Shooting 方法（§6.2.1）的几何直观 + 数学严格性是教学典范；（2）奇异性摄动（§6.3.1, §6.4.2）把"快-慢"分离用得淋漓尽致；（3）HH 行波（§6.3.2）的"0.5% 误差"是建模哲学的高光；（4）运动学理论（§6.4.3）把"细节 ODE"化简为"1 维动力学"——这是工程思想的胜利；（5）eikonal-curvature（§6.5.1）跨越应用数学与生物医学——是 Ch6 最有深远影响的章节。
  • 不足：（1）螺旋波的严格存在性（§6.5.2 末段）仍是开放问题——作者承认"奇异性摄动解是渐近的，不是严格的"；（2）高维各向异性（Exercise 22）只在习题中提及——真实心肌组织是各向异性的；（3）除颤理论（临床应用）没有专门讨论——这是 Ch12 的主题但应该在 Ch6 有所铺垫；（4）随机行波（Ch3 随机 Markov + Ch5 HH + Ch6 行波 = 随机行波）没有讨论——真实神经元的有限通道数让传导速度有涨落。

若能与作者对话，我会问：

  1. 螺旋波严格存在性：奇异性摄动解的"渐近有效性"是否在所有参数下都成立？是否有"反例"（如解不存在或不稳定）？Kessler-Kupferman 1996 的"完全严格"是否已经覆盖了 FHN 的所有参数空间？
  2. 3D 卷波与心律失常：真实心脏的 3D 卷波与 2D 螺旋波有何本质差异？除颤理论是否需要 3D 模型？卷波环的自持续需要什么条件？
  3. 离散系统传播失败的"现实"：多发性硬化的脱髓鞘 → 节点间 \(D\) 降低 → 接近 \(D^*\) 阈值。临床上这个阈值如何测量？MRI 的弥散张量成像 (DTI) 能否预测传导失败？
  4. 随机行波：真实神经元有 100-1000 个 Na⁺ 通道 / μm²（Ch6 Table 6.1）——单通道 Markov 的涨落（Ch3）如何影响行波速度？噪声能否"帮助"或"破坏"行波？
  5. 高维各向异性：心肌纤维是高度各向异性的（纤维方向 vs 横向 \(D\) 差 10 倍）。eikonal-curvature 方程 (Exercise 22) 的各向异性修正如何影响螺旋波形态？临床除颤电场的"方向性"是否利用了这个各向异性？

重要参考文献

[X1] Hodgkin AL, Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol 1952;117:500–544. （Ch5 已引——§6.3.2 的 HH 行波预测。）

[X2] Hodgkin AL. The Conduction of the Nervous Impulse. Liverpool University Press, 1964. （动作电位传播的早期综述。）

[X3] Hodgkin AL, Rushton WAH. The electrical constants of a crustacean nerve fibre. Proc R Soc Lond B 1946;133:444–479. （Ch4 已引——电缆方程在神经中的应用。）

[X4] FitzHugh R. Impulses and physiological states in theoretical models of nerve membrane. Biophys J 1961;1:445–466. （Ch5 已引——§6.3 FHN 简化的源头。）

[X5] Rinzel J, Keller JB. Traveling wave solutions of a nerve conduction equation. Biophys J 1973;13:1313–1337. （§6.3 分段线性 FHN 行波解。）

[X6] Casten RG, Cohen H, Lagerstrom PA. Perturbation analysis of an approximation to Hodgkin-Huxley theory. Q Appl Math 1975;32:365–402. （§6.3.1 奇异性摄动分析的早期工作。）

[X7] Keener JP. Waves in excitable media. SIAM J Appl Math 1980a;39:528–548. （§6.3.1 FHN 奇异性摄动分析的完整论文。）

[X8] Rauch J, Smoller J. Qualitative theory of the FitzHugh-Nagumo equations. Adv Math 1978;27:12–44.

[X9] Hastings SP. The existence of homoclinic solutions for the FitzHugh-Nagumo equations. Q J Mech Appl Math 1975;28:387–394. （HH 行波严格存在性的经典证明。）

[X10] Carpenter B. A geometric approach to singular perturbation problems of the nerve equation. J Differ Equ 1977;23:335–363.

[X11] Jones CKR. Stability of the traveling wave solution of the FitzHugh-Nagumo equations. Trans Am Math Soc 1984;286:431–469. （快脉冲稳定性。）

[X12] Maginu K. Stability of travelling wave solutions of the FitzHugh-Nagumo equations. Jpn J Appl Math 1985;2:19–37. （慢脉冲不稳定性。）

[X13] Keener JP. Propagation and its failure in coupled systems of discrete excitable cells. SIAM J Appl Math 1987;47:556–572. （§6.2.3 离散传播失败的严格分析。）

[X14] Rinzel J, Maginu K. Kinematic analysis of wave pattern formation in excitable media. In: Nonlinear Diffusion, Fife PC, Smoller J (eds.), Springer, 1984. （§6.4.3 运动学理论。）

[X15] Dockery JD, Keener JP. Diffusive effects on dispersion in excitable media. SIAM J Appl Math 1989;49:539–566.

[X16] Tyson JJ, Keener JP. Singular perturbation theory of traveling waves in excitable media (a review). Physica D 1988;32:327–361. （§6.3.1 综述。）

[X17] Murray JD. Mathematical Biology (3rd ed., 2 vols). Springer, 2002/2003. （Ch1 已引——行波传播的标准教材。）

[X18] Britton NF. Reaction-Diffusion Equations and Their Applications to Biology. Academic Press, 1986. （反应-扩散方程的应用数学综述。）

[X19] Grindrod P. Patterns and Waves: The Theory and Applications of Reaction-Diffusion Equations. Oxford University Press, 1991.

[X20] Fife PC. Mathematical Aspects of Reacting and Diffusing Systems. Lecture Notes in Biomathematics 28, Springer, 1979. （行波理论的经典。）

[X21] Keener JP, Tyson JJ. The dynamics of scroll waves. Physica D 1988;31:307–324. （§6.5.2 卷波理论。）

[X22] Winfree AT. The Geometry of Biological Time (2nd ed.). Springer, 2001. （生物时间几何——螺旋波的科普 + 数学。）

[X23] Winfree AT. Varieties of spiral wave behavior: An experimentalist's approach to the theory of excitable media. Chaos 1991;1:303–334.

[X24] Barkley D. Spiral wave dynamics. In: Chemical Waves and Patterns, Kapral R, Showalter K (eds.), Kluwer, 1995, pp. 163–189.

[X25] Karma A. Spiral wave annhilation by paced stimuli: a resonance effect. Phys Rev Lett 1994;72:3813–3816. （除颤理论。）

[X26] Foerster P, Müller SC, Hess B. Curvature and propagation velocity of chemical waves. Science 1988;241:685–687. （eikonal-curvature 方程的实验验证。）

[X27] Tyson JJ, Fife PC. Target patterns in a realistic model of the Belousov-Zhabotinskii reaction. J Chem Phys 1980;73:2224–2237.

[X28] Jafri MS, Keizer J. On the roles of Ca²⁺ diffusion, Ca²⁺ buffers, and the endoplasmic reticulum in IP3-induced Ca²⁺ waves. Biophys J 1995;69:2139–2153. （钙波中的 eikonal 方程。）

[X29] Courtemanche M, Glass L, Keener JP. A model of a two-dimensional reentrant cardiac arrhythmia. Phys Rev Lett 1993;69:3031–3034. （Beeler-Reuter 心脏模型的螺旋波。）

[X30] Lechleiter JD, Clapham DE. Molecular mechanisms of intracellular calcium excitability in Xenopus oocytes. Cell 1992;69:283–294. （Xenopus 卵母细胞钙波。）

[X31] Cytrynbaum E, Keener JP. Kinematic theory of reaction-diffusion wave fronts. SIAM J Appl Math 2002;62:1064–1080. （§6.4.3 运动学理论的进一步发展。）

[X32] Pelce P, Sun J. Wave front interaction in steadily rotating spirals. Physica D 1991;48:353–366.

[X33] Keener JP. The core of the spiral wave: a geometrical study. SIAM J Appl Math 1992;52:1370–1390.

[X34] Kessler DA, Kupferman R. Spirals in excitable media: a kinematic approach. Physica D 1996;97:509–516.

[X35] Zinner B. Existence of traveling wave solutions for the bistable differential-difference equation. J Differ Equ 1992;96:203–215.

[X36] Burton WK, Cabrera N, Frank FC. The growth of crystals and the equilibrium structure of their surfaces. Philos Trans R Soc Lond A 1951;243:299–358. （eikonal 方程在晶体生长中的应用。）

[X37] Frankel D, Sivashinsky GI. On effects due to thermal expansion and non-uniform heating in wave propagation in reacting media. Combust Sci Technol 1987;55:331–340.

[X38] Greenberg JM. A note on the Nagumo equation. Q J Mech Appl Math 1972;25:225–230.