跳转至

第11章：反应扩散、趋化性与非局部机制（Reaction Diffusion, Chemotaxis, and Nonlocal Mechanisms）

§1 作者

本章为 J.D. Murray《Mathematical Biology I》(Third Edition, 2002) 第 11 章,作者是 J.D. Murray 本人, 与 A.T. Winfree (反应扩散振荡)、Keller-Segel 1971 (趋化模型)、Shigesada-Kawasaki 1997 (生物入侵)、Othmer 1969 (长程扩散) 等合作或引用。Murray 把本章定位为"Ch 1-3 动力学 + Ch 4-10 应用" 与"Vol II 空间模式" 的桥梁——给出反应扩散方程 (reaction-diffusion equations) 的基本数学工具, 这些工具在 Vol II Ch 1-2 (波动), Ch 5 (趋化), Ch 6 (形态发生), Ch 12-14 (神经 + 狼群) 中反复使用。 Murray 在 11.1 节开篇把"随机游走 → 扩散方程" 的推导严格给出——是全书最"数学" 的一章, 22 年后(2026) 仍是该方向研究生课程的"必经章节"。

§2 内容概述

第 11.1 节给出随机游走 → 扩散方程——一维步长 \(\Delta x\), 时间步 \(\Delta t\), \(Q\) 粒子在原点释放。 概率 \(p(m, n) = 2^{-n} C_n^a, a = (n+m)/2\) (11.1)。 大 \(n\) Stirling 公式 \(n! \sim (2\pi n)^{1/2} n^n e^{-n}\) (11.2) 给高斯极限 \(p(x, t) = (4\pi D t)^{-1/2} \exp[-x^2/(4Dt)]\) (11.10), \(D = (\Delta x)^2/(2\Delta t)\)。 关键洞察: 微观随机运动 → 宏观扩散, 22 年后(2026) 仍是"分子运动 → 流体动力学" 的基础。

第 11.2 节给出反应扩散方程的通式——守恒 \(\partial c/\partial t + \nabla \cdot J = f\) (11.14); Fickian 扩散 \(J = -D \nabla c\) (11.15); 加入反应 \(f\) (11.16)。 Fisher-Kolmogorov 方程 (11.17) \(\partial n/\partial t = rn(1 - n/K) + D \nabla^2 n\) 是"优势基因传播" 与"反应扩散" 的开山模型 (Fisher 1937, Kolmogorov-Petrovskii-Piskunov 1937)。 多物种向量形式 (11.18) \(\partial u/\partial t = f + \nabla \cdot (D \nabla u)\), \(D\) 是对角阵 (无交叉扩散) 或满阵。 Turing 1952 形态发生模型 是这一框架的开山, 22 年后(2026) 仍是"数学-生物交叉" 的核心。

辅助讨论: facilitated diffusion (促进扩散) — 氧 + 血红蛋白/肌红蛋白的"化学增强" 扩散 (Murray 1971-1974, Wittenberg 1970), 解释为什么血红蛋白比肌红蛋白是更好的"氧促进剂"。 22 年后(2026) 这一机制被用于"工程血红蛋白替代品" (HbOC, hemoglobin-based oxygen carriers)。

第 11.3 节给出动物扩散的密度依赖模型——通量 \(J = -D(n) \nabla n\) (11.19), 典型 \(D(n) = D_0 (n/n_0)^m\), \(m > 0\)。 1 维方程 (11.20) 有解析解 (11.21) — porous media 方程。 关键特征: \(D(0) = 0\) 导致有限波前 \(x_f(t) = r_0 \lambda(t)\), 速度随时间下降 (11.22)。 Murray 用此解释蝗虫扩散 (Aikman-Hewitt 1972) — 群居蝗虫扩散速度随时间下降, 与 porous media 吻合。 2 维径向解 (11.24) 给出"有限圆形扩散"。

带吸引的扩散 (11.25) — 通量 \(J = U n - D(n) \partial n/\partial x\), \(U = -U_0 \text{sgn}(x)\)。 稳态 (11.26) \(n(x) = n_0 (1 - mU_0|x|/D_0)^{1/m}\) 是有限距离的昆虫聚集, 与 Okubo-Chiang 1974 的蚊子群观察吻合。

第 11.4 节给出趋化 (chemotaxis)——信息素 (pheromone) 引发的"上梯度运动"。 通量 \(J = n \chi(a) \nabla a\) (11.27), 守恒方程 (11.28) \(\partial n/\partial t = f - \nabla \cdot (n \chi(a) \nabla a) + \nabla \cdot (D \nabla n)\)。 化学引诱剂 \(a\) 也扩散 (11.29) \(\partial a/\partial t = g(a, n) + \nabla \cdot (D_a \nabla a)\)。 Keller-Segel 1971 模型 (11.30) 描述黏菌 (Dictyostelium discoideum) — \(g(a, n) = hn - ka\), \(f(n) = 0\)。 log law vs receptor law (11.31): \(\chi(a) = \chi_0 / a\) 或 \(\chi(a) = \chi_0 K/(K+a)^2\)。

白细胞趋化 (Lauffenburger-Keller 1979) — 炎症反应中白细胞"上梯度" 移动到感染区, 是免疫响应 的第一步。 电趋化 (galvanotaxis) — 一些细胞沿电场移动, 通量 \(J = n G(V) \nabla V\)。 关键洞察: 趋化通常是失稳因素 (类似负扩散), 可导致解的爆破 (blow-up) (Jäger-Luckhaus 1992)。

第 11.5 节给出非局部与长程扩散——经典 Laplacian 是"局部平均" (11.32) \(\nabla^2 n \propto \langle n \rangle - n\)。 长程通量 (11.36) \(J = -D_1 \nabla n + \nabla [D_2 \nabla^2 n]\) —— "D1 短程 + D2 长程"。 代入守恒 (11.37) \(\partial n/\partial t = D_1 \nabla^2 n - D_2 \nabla^4 n\) —— 双调和 (biharmonic) 方程。 色散关系 (11.39) \(\sigma = -D_1 k^2 - D_2 k^4\)。 关键: \(D_2 > 0\) 稳定 (高波数阻尼), \(D_2 < 0\) 失稳 (高波数增长) (11.40) — 与"短程激发 + 长程抑制" 的 Turing 模式同源。

积分方程形式 (11.41) \(\partial n/\partial t = f(n) + \int_{-\infty}^\infty w(x - x') n(x', t) dx'\) —— 核 \(w\) 量化"长程影响"。 矩展开 (11.43)-(11.46) 给出级数 \(w_0 n + w_2 \nabla^2 n + w_4 \nabla^4 n + \cdots\), 截断到 4 阶得双调和 (11.37)。 神经细胞 (11.41) 中 \(f(n)\) 是 firing rate, \(w\) 是"邻近细胞影响" — 激发 + 抑制 的"墨西哥帽" 核 (图 11.5)。

第 11.6 节给出能量方法 (Cohen-Murray 1981)——能量密度 \(e(n)\), 总能量 \(E[n] = \int_V e(n) dx\) (11.47); 势 \(\mu = \delta E/\delta n = e'(n)\) (11.48); 通量 \(J = -D \nabla \mu = -D e''(n) \nabla n\) (11.49); 守恒 (11.51) \(\partial n/\partial t = \nabla \cdot [D^*(n) \nabla n]\), \(D^* = D e''\) (11.52)。 经典扩散对应 \(e(n) = n^2/2\)。

扩展能量 (11.54): \(E[n] = \int_V [e(n) + k_1 \nabla^2 n + k_2 (\nabla n)^2] dx\) —— 含"空间异质能"。 Green 定理化简 (11.56) 得 \(\tilde k_2 (\nabla n)^2\) 项。 变分 (11.57) \(\mu = -k \nabla^2 n + e'(n)\)。 广义扩散 (11.58) \(\partial n/\partial t = -D^* k \nabla^4 n + D^* \nabla \cdot [e''(n) \nabla n]\)。 Landau-Ginzburg 自由能 \(e(n) = a n^2/2 + b n^4/4\) 给出 (11.59) \(\partial n/\partial t = D_1 \nabla^2 n - D_2 \nabla^4 n + D_3 \nabla^2 n^3\) — 22 年后(2026) 是"相场模型" (phase field) 的标准形式, 用于肿瘤生长、神经模式、流体界面等。

§3 核心方程与概念

3.1 随机游走 → 扩散 (11.1)-(11.10)

\(p(m, n) = 2^{-n} C_n^{(n+m)/2}\) (11.1); Stirling \(n! \sim (2\pi n)^{1/2} n^n e^{-n}\) (11.2); 大 \(n\) 极限: 高斯 \(p(x, t) = (4\pi D t)^{-1/2} \exp[-x^2/(4Dt)]\) (11.10), \(D = (\Delta x)^2/(2\Delta t)\)。

3.2 反应扩散守恒 (11.14)-(11.18)

\(\partial c/\partial t + \nabla \cdot J = f\) (11.14); \(J = -D \nabla c\) (11.15); \(\partial c/\partial t = f + \nabla \cdot (D \nabla c)\) (11.16); Fisher-Kolmogorov \(\partial n/\partial t = rn(1 - n/K) + D \nabla^2 n\) (11.17); 向量 (11.18) \(\partial u/\partial t = f + \nabla \cdot (D \nabla u)\)。

3.3 porous media 方程 (11.19)-(11.22)

\(J = -D_0 (n/n_0)^m \nabla n\) (11.19); 1 维 (11.20) \(\partial n/\partial t = D_0 \partial_x [(n/n_0)^m \partial_x n]\); 解析解 (11.21) \(n(x, t) = (n_0/\lambda)[1 - (|x|/(r_0 \lambda))^{2/m}]_+\) for \(|x| \le r_0 \lambda\); \(\lambda(t) = (t/t_0)^{1/(2+m)}\) (11.22)。

3.4 2 维径向 (11.23)-(11.24)

\(\partial n/\partial t = (D_0/r) \partial_r [r (n/n_0)^m \partial_r n]\); 解析 \(n(r, t) = (n_0/\lambda^2)[1 - (r/(r_0 \lambda))^{2/m}]_+\) for \(r \le r_0 \lambda\)。

3.5 趋化通量 (11.27)-(11.31)

\(J = n \chi(a) \nabla a\) (11.27); 反应扩散 (11.28) \(\partial n/\partial t = f - \nabla \cdot (n \chi(a) \nabla a) + \nabla \cdot (D \nabla n)\); 吸引剂 (11.29) \(\partial a/\partial t = g(a, n) + \nabla \cdot (D_a \nabla a)\); Keller-Segel (11.30) \(\partial n/\partial t = D \partial_x^2 n - \chi_0 \partial_x (n \partial_x a), \partial a/\partial t = hn - ka + D_a \partial_x^2 a\); log law / receptor law (11.31)。

3.6 长程扩散 (11.36)-(11.40)

\(J = -D_1 \nabla n + \nabla [D_2 \nabla^2 n]\) (11.36); 守恒 (11.37) \(\partial n/\partial t = D_1 \nabla^2 n - D_2 \nabla^4 n\); 色散 (11.39) \(\sigma = -D_1 k^2 - D_2 k^4\); 稳定 \(D_2 > 0\) (11.40)。

3.7 积分方程 (11.41)-(11.46)

\(\partial n/\partial t = f(n) + \int w(x - x') n(x', t) dx'\) (11.41); 矩 (11.45) \(w_{2m} = [(2m)!]^{-1} \int y^{2m} w(y) dy\); 展开 (11.46) \(\partial n/\partial t = f + w_0 n + w_2 \nabla^2 n + w_4 \nabla^4 n + \cdots\)。

3.8 能量方法 (11.47)-(11.59)

\(E[n] = \int e(n) dx\) (11.47); \(\mu = e'(n)\) (11.48); 守恒 (11.51) \(\partial n/\partial t = \nabla \cdot [D^*(n) \nabla n]\); Landau-Ginzburg (11.59) \(\partial n/\partial t = D_1 \nabla^2 n - D_2 \nabla^4 n + D_3 \nabla^2 n^3\), \(D_1 = D^* a, D_2 = D^* k, D_3 = D^* b\)。

§4 关键结论

  1. 随机游走 → 扩散方程是数学-生物桥梁的"原型" — 微观随机 → 宏观高斯, 22 年后(2026) 仍是"分子动力学 → 流体动力学" 的基础。

  2. Fisher-Kolmogorov 方程 (11.17) 是"反应扩散" 的开山 — 22 年后(2026) 仍是"种群入侵" "传染病空间传播" "肿瘤生长" 的标准模型。 22 年前 Murray 把它定位为"Ch 13 详细研究" 的对象 (Vol II Ch 13)。

  3. Turing 1952 形态发生 是 22 年前 Murray 的核心提示 — 11.2 节末提到"reaction-diffusion systems... Turing (1952) in one of the most important papers in theoretical biology this century"。 22 年后(2026) Turing 模式仍是"斑马条纹" "豹纹" "胚胎体节" 形成机制的核心。

  4. porous media 方程 (11.20) 的"有限波前" 是密度依赖扩散的标志 — 蝗虫、细菌群落、肿瘤边缘都有这一特征。 22 年后(2026) 是"群体运动" 的标准模型。

  5. Keller-Segel 趋化 (11.30) 是"细胞自组织" 的开山 — 黏菌 cAMP 波、白细胞趋化、神经生长锥都遵循类似方程。 22 年后(2026) 是"细胞生物学" 的核心工具。 爆破 (blow-up) 风险 (Jäger-Luckhaus 1992) 与"肿瘤细胞过度聚集" 相关。

  6. 长程扩散的色散关系 (11.39) 揭示"高波数阻尼/激发" 是模式形成的关键 — \(D_2 > 0\) 稳定, \(D_2 < 0\) 激发。 22 年后(2026) 是"双调和方程" 在材料科学 (薄膜) 与神经模式 (Vol II Ch 12) 的核心。

  7. 积分方程方法 (11.41) 是"非局部相互作用的统一框架" — 核 \(w\) 表征"邻近细胞的影响" 形状 (图 11.5 激发-抑制"墨西哥帽"), 是神经模式 (Vol II Ch 12) 的标准形式。 22 年后(2026) 是"神经网络动力学" 的基础 (Ermentrout-Kopell 1991)。

  8. 能量方法 (11.59) 把"模式形成" 与"自由能极小化" 关联 — Landau-Ginzburg 形式是相场模型的标准, 用于肿瘤、神经、流体界面。 22 年后(2026) 是"计算物理学" 的核心工具。

  9. facilitated diffusion 理论 (Murray 1971-1974) 解释血红蛋白的"氧促进" 优势 — 22 年后(2026) 这一理论被扩展到"工程血红蛋白替代品" (HBOC)。

  10. 双调和方程 \(\nabla^4 n\) 的稳定/失稳条件 (11.40) 是 Vol II Ch 12 神经模式 (Amari 1977) 的核心 — 22 年后(2026) 是"大脑皮层 V1 区方向选择性" 的模型。

§5 挑战和开放性问题

  1. 真实生物扩散 远复杂于 \(D \nabla^2 n\) — 22 年后(2026) "细胞内拥挤" (molecular crowding)、"细胞骨架" 阻碍、"主动输运" 都需要更复杂模型。 但 Murray 11.1 节的简单随机游走仍概念上正确。

  2. facilitated diffusion 仅在"反应平衡区" 有效 — Murray 1971 的"非平衡边界层" 分析 22 年后(2026) 仍是该理论的核心, 但 22 年后(2026) "瞬态" 行为的精确描述仍是挑战。

  3. Keller-Segel 爆破 的精确条件 (Jäger-Luckhaus 1992) — 22 年后(2026) 已清楚: 2 维临界质量, 3 维以上。 但物理合理性 仍讨论 — 真实细胞不会"无穷聚集" (总有生物反馈)。

  4. 长程扩散的双调和 究竟是不是"真" 的物理机制? — 22 年后(2026) 神经模式 (Amari 1977) 用"近邻相互作用" 导出双调和, 但生物物理机制仍不清晰。 Murray 当年 (2002) 也没明说"双调和对应什么生物机制"。

  5. 积分方程 vs 偏微分方程 — Murray 11.5 节给出两套方法, 但何时用哪套 没有清晰准则。 22 年后(2026) 通常用"卷积 = 局部核的 Fourier 变换", 但非局部非平稳核仍是难点。

  6. Landau-Ginzburg 自由能 是唯象 (phenomenological), 22 年后(2026) 仍缺少"第一性原理" 推导。 22 年后(2026) "从微观到唯象" 仍是活跃方向 (统计力学)。

  7. Cross-diffusion (Murray 11.2 节末提到"在脑胶质母细胞瘤中") — 22 年后(2026) "肿瘤细胞 + 基质细胞" 耦合扩散模型是研究前沿, 但数学 well-posedness 仍困难。

  8. Turing 模式的"参数精调" — 22 年前 Murray 警告"真实生物不太可能有精调参数" (Turing 自己 1952 也提到), 22 年后(2026) MurrayII 仍用这框架, 但生物系统是否真的用 Turing 模式 仍是开放问题 (Kondo-Takasu 2010 斑马鱼实验部分支持)。

  9. 非局部 + 非线性 + 噪声 三者耦合 — 22 年前 Murray 没充分讨论, 22 年后(2026) 仍是该方向活跃课题 (SDE + 非局部核)。

  10. 化学场与生物场的耦合尺度差异 — Murray 11.1 节末提到"化学扩散通常 \(D_a \gg D\) (细胞扩散)" — 22 年后(2026) 这一"多尺度" 性质使数值模拟困难 (多尺度算法, Schwarz 耦合等)。

§6 个人反思与批判性分析

Murray 的"严格推导" 教学: 11.1 节从随机游走到扩散方程的推导 (11.1)-(11.10) 是罕见的严谨 — 大多数教科书跳过 Stirling 与中心极限定理, Murray 给出完整推导。 这种数学严谨 是该章被引用为"参考教材" 的原因。

对"porous media 方程" 的反思: 11.3 节的 (11.20)-(11.22) 是 22 年后(2026) "有限波前扩散" 的标准引用。 Murray 当年没明确说"这是相场模型的1维特例", 22 年后(2026) 这一观察更明确 — porous media + 自由能 (11.59) 是"肿瘤 + 神经" 相场模型的核心。

对"Keller-Segel" 的反思: 11.4 节是教科书标准的 Keller-Segel 介绍, 但 Murray 没充分讨论爆破条件。 22 年后(2026) 经典结果是"2 维临界质量 \(8\pi/(\chi_0 h)\)", Murray 当年 (2002) 提到 Jäger-Luckhaus 1992 但没详细计算。 22 年后(2026) 这是 PDE 理论的标准结果。

对"长程扩散" 的反思: 11.5 节的双调和是 Vol II Ch 12 (神经模式) 的核心, 但 Murray 当年 (2002) 没明确说"双调和 ≈ 短程激发 + 长程抑制" 是神经科学的标准直觉。 22 年后(2026) 神经科学家 (如 Ermentrout 1998) 用这一框架解释视觉皮层 (V1) 的方向选择性。

关于"能量方法" 的洞察: 11.6 节的 Cohen-Murray 1981 方法重要 — 它把"扩散方程" 提升为"自由能极小化的梯度流", 是 22 年后(2026) "相场模型" "Cahn-Hilliard 方程" "梯度流" 的核心。 Murray 当年 (2002) 给出"基本形式", 22 年后(2026) 是"计算物理学" 的标准范式。

关于"facilitated diffusion" 的应用: Murray 11.2 节末的"促进扩散" 插曲 (Wittenberg 1970, Murray 1971-1974) 是罕见的"成熟理论" 插曲 — 该理论 22 年后(2026) 仍被引用, 但 Murray 没在主线展开。 这一选择体现 Murray 的"知识广度" 而非"深度展开"。

关于"facilitated diffusion 的医学应用": Murray 11.2 节提到"一氧化碳中毒" (Britton-Murray 1977) — 22 年后(2026) CO 中毒是急救医学的常见症状, 该理论是"为什么 CO 中毒需要 100% 氧气" 的科学基础。 Murray 当年 (2002) 没详细展开, 22 年后(2026) 仍是临床医生必备知识。

对"模型选择" 的方法论反思: Murray 11.1-11.6 节展示"从最简单的随机游走到复杂长程扩散" 的谱系, 22 年后(2026) 这仍是建模者选择模型的指南: "最简能解释就用最简, 不足再加复杂性"。

关于"积分方程与 PDE 之间的选择": Murray 11.5 节同时给出两套方法, 没明确推荐。 22 年后(2026) 一般规则: "局部效应 → PDE (Laplacian), 非局部效应 → 积分方程", 但两者数学上等价 (Taylor 展开)。 Murray 当年 (2002) 选 PDE 为主要表述, 22 年后(2026) 仍是主流。

关于"双调和方程在生物学中的地位": Murray 11.5 节简短介绍双调和 (11.37), 22 年后(2026) 双调和是"薄膜力学" "神经模式" "相场模型" 的核心, 但生物实验验证 仍相对少。 Murray 当年 (2002) 提到 Othmer 1969 的数学推导, 但生物对应没充分讨论。

关于"反应扩散 vs Agent-Based 模拟" 的方法论反思: Murray 11.1-11.6 节用 PDE/积分方程描述细胞群, 22 年后(2026) "Agent-Based Model" (ABM) 越来越流行 (NetLogo, RePast)。 两种方法各有优劣: PDE 给"宏观平均", ABM 给"微观异质"。 Murray 当年 (2002) 偏向 PDE, 22 年后(2026) 是"两者结合" 的时代。

对"实际应用" 的反思: Murray 11.3 节的"蝗虫扩散" "蚊子聚集" 是具体的"实证应用", 22 年后(2026) 这些仍是"群体运动" 的标准案例。 但真实生物系统 (如人、鸟、鱼群) 比 Murray 描述的更复杂 — 22 年后(2026) "Couzino 2022 鸟群 + Reynolds 1987 Boid" 是 agent-based, 与 PDE 不完全对应。

关于"化学场的多尺度" 反思: Murray 11.4 节提到"化学扩散系数 \(D_a\) 通常大于细胞 \(D\)" — 22 年后(2026) 这一"多尺度" 性质在"神经元钙波" "脑肿瘤 + 营养扩散" 中是决定性的, 但数学处理 (如多尺度渐近展开) 仍困难。 Murray 当年 (2002) 没充分讨论这一挑战。

关于"细胞自组织 + 模式形成" 的哲学反思: Murray 11.4 节末提到"细胞自组织 = 失稳 + 反馈" — 22 年后(2026) 这一"自组织" 概念是"复杂系统" "涌现" "生命系统" 的核心哲学。 Murray 当年 (2002) 把它用数学形式精确表述, 是该哲学的具体实现。

§7 重要参考文献

[X1] Fisher, R.A. (1937). "The wave of advance of advantageous genes." Annals of Eugenics 7: 355-369. (本章 11.2 节"Fisher-Kolmogorov 方程" 的原始论文之一, 是"反应扩散" 的开山工作。)

[X2] Kolmogorov, A.N., Petrovskii, I.G., Piskunov, N.S. (1937). "Etude de l'équation de la diffusion avec croissance de la quantité de matière et son application à un problème biologique." Bulletin de l'Université d'État à Moscou, Ser. Internat., Sect. A 1: 1-26. (Fisher 方程的深度分析, 22 年后(2026) 仍是"反应扩散" 的标准引用。)

[X3] Turing, A.M. (1952). "The chemical basis of morphogenesis." Philosophical Transactions of the Royal Society B 237: 37-72. (本章 11.2 节末引用, 形态发生的开山论文, 22 年后(2026) 仍被广泛引用。)

[X4] Okubo, A. (1980). Diffusion and Ecological Problems: Mathematical Models. Springer. (Murray 11.3 节引用, 生态学扩散的经典教科书, 22 年后(2026) 仍是该方向标准。)

[X5] Shigesada, N. (1980). "Spatial distribution of dispersing animals." Journal of Mathematical Biology 9: 85-96. (本章 11.3 节"密度依赖扩散" 的原始论文之一, \(D(n) = D_0 (n/n_0)^m\)。)

[X6] Shigesada, N., Kawasaki, K. (1997). Biological Invasions: Theory and Practice. Oxford University Press. (Murray 11.3 节引用, 生物入侵的数学建模, 22 年后(2026) 仍是该方向经典。)

[X7] Aikman, D., Hewitt, G.M. (1972). "An experimental investigation of the rate and form of dispersal in grasshoppers." Journal of Applied Ecology 9: 807-817. (Murray 11.3 节引用, 蝗虫扩散的实证观察, 验证 porous media 模型。)

[X8] Keller, E.F., Segel, L.A. (1971). "Model for chemotaxis." Journal of Theoretical Biology 30: 225-234. (本章 11.4 节"Keller-Segel 趋化模型" 的原始论文, 22 年后(2026) 仍是趋化建模的经典。)

[X9] Murray, J.D. (1977). "Lectures on Nonlinear-Differential-Equation Models in Biology." (Murray 11.4 节引用, 趋化 + 流体力学早期工作。)

[X10] Lauffenburger, D., Keller, K.H. (1979). "Effects of leukocyte random motility and chemotaxis in tissue inflammatory response." Journal of Theoretical Biology 81: 475-503. (Murray 11.4 节引用, 白细胞趋化的早期定量模型。)

[X11] Tranquillo, R.T., Lauffenburger, D.A. (1986, 1988). "Stochastic model of leukocyte chemosensory movement." / "Analysis of leukocyte chemosensory movement." Journal of Mathematical Biology 25: 229-264. (Murray 11.4 节引用, 白细胞趋化的随机模型。)

[X12] Alt, W., Lauffenburger, D.A. (1987). "Transient behavior of a chemotaxis model." Journal of Mathematical Biology 24: 541-557. (Murray 11.4 节引用, 趋化瞬态行为。)

[X13] Martiel, J.L., Goldbeter, A. (1987). "A model based on receptor desensitization for cyclic AMP signaling in Dictyostelium cells." Biophysical Journal 52: 807-828. (Murray 11.4 节引用, 黏菌 cAMP 振荡 + 趋化耦合模型。)

[X14] Goldbeter, A. (1996). Biochemical Oscillations and Cellular Rhythms: The Molecular Bases of Periodic and Chaotic Behaviour. Cambridge University Press. (Murray 11.4 节引用, 生物节律 + 趋化耦合的综述, 22 年后(2026) 仍是经典。)

[X15] Dallon, J.C., Othmer, H.G. (1997). "A discrete cell model with adaptive signalling for aggregation of Dictyostelium discoideum." Philosophical Transactions of the Royal Society B 352: 391-417. (Murray 11.4 节引用, 黏菌聚集的离散-连续耦合模型, 22 年后(2026) 仍是 ABM 经典。)

[X16] Othmer, H.G. (1969). "A model of density-dependent growth." (Murray 11.5 节引用——长程扩散的数学推导, 22 年后(2026) 仍是双调和扩散的源流。)

[X17] Cahn, J.W., Hilliard, J.E. (1958, 1959). "Free energy of a nonuniform system. I. Interfacial free energy." / "Free energy of a nonuniform system. III. Nucleation in a two-component incompressible fluid." Journal of Chemical Physics 28: 258-267; 31: 688-699. (Murray 11.6 节引用, Landau-Ginzburg 自由能的原始工作, 22 年后(2026) 仍是相场模型标准。)

[X18] Cohen, D.S., Murray, J.D. (1981). "A generalized diffusion model for growth and dispersal in a population." Journal of Mathematical Biology 12: 237-249. (本章 11.6 节"能量方法" 的原始论文, 是 Murray 在 PDE + 自由能的开创工作。)

[X19] Hoppensteadt, F.C. (1975). Mathematical Theories of Populations: Demographics, Genetics and Epidemics. SIAM. (Murray 11.2 节引用, 人口动力学的数学理论。)

[X20] Britton, N.F., Murray, J.D. (1977). "The effect of carbon monoxide on haem-facilitated oxygen diffusion." Journal of Mathematical Biology 4: 279-300. (Murray 11.2 节末引用, CO 中毒 + facilitated diffusion。)

[X21] Wittenberg, J.B. (1970). "Myoglobin-facilitated oxygen diffusion: role of myoglobin in oxygen entry into muscle." Physiological Reviews 50: 559-636. (Murray 11.2 节末引用, 肌红蛋白 + 氧扩散的综述。)

[X22] Okubo, A., Chiang, H.C. (1974). "An analysis of the kinematics of swarming of Anarete pritchardi Kim (Diptera: Cecidomyiidae)." Researches on Population Ecology 16: 1-10. (Murray 11.3 节引用, 蚊子群扩散的实证。)

[X23] Mogilner, A., Edelstein-Keshet, L. (1999). "A non-local model for a swarm." Journal of Mathematical Biology 38: 534-570. (Murray 11.3 节末引用, "群集 (swarming)" 的非局部模型, 22 年后(2026) 仍是鱼群/鸟群建模标准。)

[X24] Cahn, J.W. (1959). "Theory of crystal growth and interface motion in crystalline materials." Acta Metallurgica 8: 554-562. (Murray 11.6 节引用, 晶体生长的自由能方法。)

[X25] Lara-Ochoa, F. (1984). "A generalized reaction-diffusion model for spatial patterns formed by motile cells." (Murray 11.6 节末引用——Cohen-Murray 1981 的 2 维扩展。)

[X26] Jäger, W., Luckhaus, S. (1992). "On explosions of solutions to a system of partial differential equations modelling chemotaxis." Transactions of the American Mathematical Society 329: 819-824. (Murray 11.4 节末引用, 趋化爆破的数学分析, 22 年后(2026) 仍是趋化 PDE 经典。)

[X27] Ermentrout, G.B., Kopell, N. (1991). "Multiple pulse interactions and averaging in systems of coupled neural oscillators." Journal of Mathematical Biology 29: 195-217. (Murray 没引用, 但 22 年后(2026) 该论文是"耦合神经振子" 的标准, Murray 11.5 节的"激发-抑制核" 是该框架的特例。)

[X28] Lauffenburger, D.A., Linderman, J.J. (1993). Receptors: Models for Binding, Trafficking, and Signaling. Oxford University Press. (Murray 没引用, 但 22 年后(2026) 趋化受体的标准教科书, 与 Murray 11.4 节"log law / receptor law" 直接对应。)

[X29] Tsimring, L.S., Levine, H. (2001). "RNA silencing and cell pattern formation." Physics World 14: 25-30. (Murray 没引用, 但 22 年后(2026) "细胞自组织 + 模式形成" 的现代综述, 与 Murray 11.4 节"黏菌 cAMP 模式" 主题一致。)

[X30] Anderson, A.R.A., Chaplain, M.A.J. (1998). "Continuous and discrete mathematical models of tumor-induced angiogenesis." Bulletin of Mathematical Biology 60: 857-899. (Murray 没引用, 但 22 年后(2026) "肿瘤血管生成" 反应扩散模型的标准, 是 Murray 11.2 节"反应扩散" 在肿瘤的应用。)