跳到主要内容
Sandbox Physics
膜与囊泡力学

L2 · 第 4 章

膜的抖动怎样显出弯曲刚度

膜越容易弯,热运动就越能把它推离平均形状。测量不同尺度的抖动,可以把弯曲刚度与张力分开:短波起伏更怕弯曲,长波起伏更容易受到张力影响。温度提供了上一章的静态形状所缺少的能量标尺。

先把一张起伏的膜拆成波纹

取半径 RR 附近的准球形膜,把径向位移除以半径,得到无量纲小形变 u(θ,ϕ)u(\theta,\phi)。球谐函数是一套在球面上振荡的标准图案,类似圆周上的正弦波。阶数 ℓ\ell 越高,表面波纹越细;同一阶的 mm 区分绕轴方向。

零阶主要改变平均半径,一阶对应整体平移。把中心移到囊泡中心后,研究从 ℓ=2\ell=2 开始的形变。以下二次理论取零自发曲率、小幅形变和均匀张力;涨落模型使用给定张力,而不逐帧强制固定过剩面积。

r(θ,ϕ)=R[1+u(θ,ϕ)],u=∑ℓ,muℓmYℓmr(\theta,\phi)=R[1+u(\theta,\phi)],\qquad u=\sum_{\ell,m}u_{\ell m}Y_{\ell m}
球谐函数按球面面积角归一化。平均半径的二阶修正可用于保持体积;大幅形变需要回到完整曲面能量。

来源:B. Loubet

一个二次自由度怎样分到热能

先看一个实数模振幅 qq。若能量为 Kqq2/2K_q q^2/2,平衡概率与 exp⁡[−Kqq2/(2kBT)]\exp[-K_q q^2/(2k_BT)] 成正比。这是以零为中心的高斯分布,比较指数便得方差 kBT/Kqk_BT/K_q。方差是振幅平方的平均值,描述长期抖动有多大。

把 Helfrich 能在球面附近展开,球谐函数的正交性会消掉不同模之间的交叉项,每个实自由度都变成这样一个二次项。曲率给高阶模更大的惩罚,张力则增加所有表面起伏的代价。

Eℓ=(ℓ−1)(ℓ+2)[ℓ(ℓ+1)+σˉ],σˉ=σR2κE_\ell=(\ell-1)(\ell+2)[\ell(\ell+1)+\bar\sigma],\qquad \bar\sigma=\frac{\sigma R^2}{\kappa}
这里的张力加横线表示无量纲;它与以力除以长度计的物理张力不同。
F(2)=κ2∑jEℓ(j)qj2,⟨qj2⟩=kBTκEℓ(j)F^{(2)}=\frac{\kappa}{2}\sum_j E_{\ell(j)}q_j^2,\qquad \langle q_j^2\rangle=\frac{k_BT}{\kappa E_{\ell(j)}}
下标逐一枚举归一化的实球谐自由度。要求各保留模的二次系数为正,热平衡方差才存在。

来源:B. Loubet

完整例题:同一个模,在两种张力下

取 κ=20kBT\kappa=20k_BT。最低形变阶 ℓ=2\ell=2 在零无量纲张力时有 E2=24E_2=24,单个归一化实模的方差为 1/480≃0.002081/480\simeq0.00208,均方根振幅约 0.04560.0456。这是球谐系数,局部径向位移还要乘半径和该处球谐函数的值。

把无量纲张力升到 σˉ=6\bar\sigma=6,得到 E2=48E_2=48,方差减半。对 ℓ=10\ell=10,张力从零升到六只把括号从 110110 改成 116116,影响较小。测到多个尺度,才容易识别哪一部分抑制来自张力。

如果只看很窄一段模数,刚度与张力可能互相补偿。增加帧数能降低随机误差;分辨两种参数还需要足够宽的尺度范围。

赤道上的一条线包含许多球面模

显微镜常跟踪赤道轮廓,相当于在 θ=π/2\theta=\pi/2 取截线,再沿方位角做 Fourier 分解。一次圆周谐波可能由多个球面阶数共同贡献。把它误认成单独一个球谐模,会改变拟合曲线。

若归一化的赤道余弦振幅写成 Am=∑ℓ≥mwℓmqℓmA_m=\sum_{\ell\geq m}w_{\ell m}q_{\ell m},其中 wℓm=2NℓmPℓm(0)w_{\ell m}=\sqrt2N_{\ell m}P_\ell^m(0),独立模的方差相加便得到下面的谱。角向正弦分量有同样方差。有限阶数模型只计算所保留的项。

真实物镜还会收集有限焦深中的荧光或轮廓信息,观测不一定等于数学上的零厚度截线。这里先把赤道切片、角向像素平均和曝光分开,下一章再跟着相机走完整条测量链。

⟨Am2⟩=∑ℓ≥max⁡(2,m)wℓm2kBTκEℓ\langle A_m^2\rangle=\sum_{\ell\geq\max(2,m)}w_{\ell m}^2\frac{k_BT}{\kappa E_\ell}
权重来自球谐函数在赤道上的值;奇偶性使其中一些权重恰好为零。局部几何、球面阶数和轮廓模数各自承担不同角色。

来源:S. A. Rautu et al. · J. Pécréaux et al.

把推理带进实验

从不同尺度的抖动辨认刚度

  1. 先预测增大刚度会怎样改变整体抖动,增大张力又最先压低哪段谱。
  2. 一次只改变一个参数,对照理想轮廓和保存后得到的谱。
  3. 导出记录,作为进阶练习离线比较只保留低模数与保留全谱的拟合。解释为什么尺度范围窄时,参数更容易互相补偿。

练习与解释性答案

  1. 温度不变、刚度加倍且无量纲张力保持不变,单模方差和均方根振幅怎样变?

    展开解释

    方差减半,均方根振幅除以根号二。固定无量纲张力意味着物理张力也要随刚度调整,不能把两种固定条件混为一谈。

  2. 一阶球谐系数随时间漂移,是否说明膜有特别软的长波形变?

    展开解释

    一阶主要对应整体平移。应先跟踪中心并校正漂移,再解释形变谱。中心定位误差也可能把信号混入低模数。

  3. 若高模数的谱突然变平,应立刻推断膜变软吗?

    展开解释

    先检查轮廓定位噪声、像素与焦深。噪声底可把真正下降的谱托平;下一章会把这些相机效应放进观测模型。

来源与复算范围

  1. J. Pécréaux et al., Refined contour analysis of giant unilamellar vesicles, European Physical Journal E 13, 277–290 (2004)

    轮廓涨落谱、相机积分与观察平面。课程自行生成带种子的有限模记录,并推导曝光积分;本课程的数值拟合针对合成记录,不是对原论文脂质测量的重新拟合。

  2. S. A. Rautu et al., The role of optical projection in the analysis of membrane fluctuations, Soft Matter 13, 3480–3483 (2017)

    式 (1)—(3)、图 1 与有限焦深的测量问题。本课程模型采用零焦深赤道切片,另行计算角向平均与曝光;完整光学投影算子作为后续方法讨论。

  3. B. Loubet, U. Seifert and M. A. Lomholt, Effective tension and fluctuations in active membranes, arXiv:1112.2919v2; Physical Review E 85, 031913 (2012)

    2011 年 12 月 14 日 v2:准球形二次模与具有指数时间相关的直接主动力。本课程的同方差响应比较采用明确给出的有限线性模型;不复算论文中自洽过剩面积张力曲线,也不覆盖渗流力和曲率力分支。

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

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