跳到主要内容
Sandbox Physics
形态发生与反应扩散

FR21 · L1–L4

形态发生与反应扩散

均匀场可以长出斑点,局部刺激可以卷成螺旋,收缩的皮层可以把材料拉成密度峰。沿着每种过程的方程往下走,再问一个更难的问题:观察什么,才能把这些机制分开?

沿着一个问题往下走

先看一小块区域里的反应和邻格交换,再学会预测、量出波纹的增长率。第四章转向激发与恢复;第五、六章加入组织生长和主动收缩。最后用一份合成干预记录比较机制,并亲手重算图灵原文中一个范围明确的结果。

每章都列出先修。可以先看现象,再回到相应推导补上方法。共同数学基础可沿教材计划按需补读。 教材计划

  1. L1 · 01

    扩散为什么有时会帮忙长出花纹

    扩散会抹平浓度差,反应却可能放大它。从相邻小格之间的交换,理解两者怎样一起选出空间图案。

  2. L2 · 02

    算出哪些波纹会长大

    从 Schnakenberg 反应的均匀稳态出发,逐步推导 Jacobian、增长谱与图灵条件,并完成一个可手算的例子。

  3. L2 · 03

    从模拟记录量出增长率

    有限区域只容得下某些波长。用保存的早期浓度场提取 Fourier 振幅,区分真实增长、网格误差与非线性饱和。

  4. L2 · 04

    一阵刺激怎样变成绕行的螺旋

    从激发阈值和恢复过程理解可激发波:为什么波能传播、相撞会消失,断开的波前又怎样绕成螺旋。

  5. L3 · 05

    组织长大时,浓度为什么会变淡

    从物质量守恒推导生长方域上的反应扩散。用随材料移动的坐标,分清扩散尺度、几何稀释与总量变化。

  6. L3 · 06

    收缩怎样把更多材料拉到一起

    用一维活性黏性皮层推导密度、应力与流速的反馈,再用局部抑制比较机械输运和扩散。

  7. L4 · 07

    怎样用一次干预分辨图案机制

    用带噪声的合成干预记录比较主动皮层与被动扩散:先校准,再预测未参与拟合的条件,检查噪声怎样削弱分辨力。

  8. L4 · 08

    重算图灵论文里的两细胞例子

    打开图灵 1952 年论文第 42–43 页,重算两细胞的反应、交换与指数增长,再用完整四变量数值记录检查结果。

四个可以比较的实验

还可以走进已有的免疫监视实验,看分子相遇怎样连到细胞和组织尺度。它研究搜索与捕获,适合作为尺度背景;这里则逐式研究反应、扩散与力学反馈。 多尺度免疫监视

课程来源与计算范围

推导、示意图与数值例由本课程编写。下面分别注明采用的方程、教学扩展和实际重算的已发表例子。交互场与带噪观测均来自所列模型。

  1. A. M. Turing, The Chemical Basis of Morphogenesis (1952), Philosophical Transactions B 237, 37–72

    核读第 4 节第 42–43 页的两细胞数值例与第 6 节的模态分析。本课程复算第 4 节早期线性增长;没有复现原文图,也没有把 Schnakenberg 实验称为原文模型。

  2. J. Schnakenberg, Simple Chemical Reaction Systems with Limit Cycle Behaviour (1979), Journal of Theoretical Biology 81, 389–400

    原始反应模型的文献入口。这里明确写出采用的无量纲、不可逆教学形式,稳态、Jacobian 与数值参数均在课程中重新推导;没有声称使用生物标定参数。

  3. D. Barkley, A Model for Fast Computer Simulation of Waves in Excitable Media (1991), Physica D 49, 61–70

    核读第 2 节的两变量方程与文末模型适用性说明。本实验使用所列连续方程及自己的离散求解器,初态和预设由本课程设计,没有复现原文图像或转动频率。

  4. E. J. Crampin, E. A. Gaffney and P. K. Maini, Reaction and Diffusion on Growing Domains (1999), Bulletin of Mathematical Biology 61, 1093–1120

    核读第 3 节方程 9–20 的守恒、输运与域增长推导。原文主要讨论一维生长;本课程另推导二维均匀膨胀方域的稀释项,作为教学扩展,不复现原文倍频图。

  5. M. Nishikawa, S. R. Naganathan, F. Jülicher and S. W. Grill, Controlling Contractile Instabilities in the Actomyosin Cortex (2017), eLife 6, e19595

    核读 Appendix 方程 4–9 的密度通量、活性黏性应力、饱和函数与力平衡。这里把原文两组分调控简化成单密度与线性周转,另加局部应力抑制;没有拟合真实胚胎数据或复现原文图 2。原论文采用 CC BY 许可。

  6. A. L. Krause, M. A. Ellis and R. A. Van Gorder, Influence of Curvature, Growth, and Anisotropy on the Evolution of Turing Patterns on Growing Manifolds (online 2018; issue 2019)

    核读方程 4–6 的移动域守恒、14–15 的材料坐标及 44–45 的 Schnakenberg 反应与 Jacobian。这里的周期方域和参数是教学选择,没有复现原文椭圆几何或图 1–4。

  7. L. Giomi, Geometry and Topology of Turbulence in Active Nematics, arXiv:1409.1555v2 (2015); Physical Review X 5, 031003

    只作前沿分支阅读:二维不可压缩活性向列流、涡旋统计与拓扑缺陷。它需要取向场和流体方程,归入 FR22;这里的一维标量皮层模型没有复现这些结果。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策