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

L2 · 第 2 章

弯一片膜要付出多少能量

膜抵抗的是偏离偏爱曲率的弯曲。把每一小块偏离的平方加起来,就得到 Helfrich 弯曲能。它解释了一个初看奇怪的结果:没有自发曲率时,大球和小球的弯曲能相同。

从每一小块开始相加

先把膜视为厚度可忽略的光滑曲面。局部偏爱的总曲率记为 c0c_0,单位是长度倒数。膜偏离它时,能量密度按 (2H−c0)2(2H-c_0)^2 增长。前面的弯曲刚度 κ\kappa 以能量计,越大就越难弯。

积分符号表示把整个表面划成越来越小的面积块 dA\mathrm dA,计算每块的能量后相加。高斯项还有自己的系数 κˉ\bar\kappa。两种系数都属于曲率弹性;直接拉长膜的代价还需要另一项弹性能。

Fb=∫A[κ2(2H−c0)2+κˉK] dAF_b=\int_A\left[\frac{\kappa}{2}(2H-c_0)^2+\bar\kappa K\right]\,\mathrm dA
局域曲率平方乘以面积后没有量纲,故两个刚度系数都具有能量单位。这里假定连续、各向同性、流体双层,忽略膜厚尺度上的细节。

来源:W. Helfrich

球越大,为什么能量未必越大

球上每一点的 2H=2/R2H=2/R 相同,可以直接从积分中提出。面积 4πR24\pi R^2 随半径平方增加,曲率平方却随半径平方减小,两者恰好抵消。先取 c0=0c_0=0 并只计算平均曲率项,就得到 8πκ8\pi\kappa。

允许自发曲率后,球面结果改为 2πκ(2−c0R)22\pi\kappa(2-c_0R)^2。此时半径重新进入能量,膜存在偏爱的弯曲尺度。翻转曲率符号约定时要同时翻转 c0c_0,物理预测才保持相同。

对于无边界且拓扑固定的闭合曲面,Gauss–Bonnet 定理给出 ∫K dA=4π(1−g)\int K\,\mathrm dA=4\pi(1-g),gg 是孔洞数。球拓扑的高斯项是常数;比较同一拓扑下的形状时它不影响排序。融合或穿孔会改变这个前提。

Fb,sphere(H)=2πκ(2−c0R)2,Fb,sphere(H)∣c0=0=8πκF_{b,\rm sphere}^{(H)}=2\pi\kappa(2-c_0R)^2,\qquad F_{b,\rm sphere}^{(H)}\big|_{c_0=0}=8\pi\kappa
上标只标识平均曲率部分。球拓扑的高斯贡献还需另加一个与形状无关的常数。

来源:U. Seifert

张力和拉伸刚度扮演不同角色

张力是增加少量面积需要付出的能量除以面积,也可理解为沿切口每单位长度的拉力。若把总面积严格固定,张力可以作为保证面积不变的拉格朗日乘子:一个在求极值时使用的辅助量,其数值由约束决定。

面积压缩模量 KAK_A 描述张力随相对面积伸长变化的快慢。设轻松状态的面积为 A0A_0,小应变为 α=(A−A0)/A0\alpha=(A-A_0)/A_0,二次拉伸能给出下面两式。张力与 KAK_A 都以能量除以面积计,但前者是当前的应力,后者是应力对应变的斜率。

囊泡常先用展开热皱褶的方式增加可见面积,随后才明显拉开脂质。实验中用哪一种面积定义应变,会影响解释;解释脂质分子之间的距离变化时,需要先分开皱褶展开和真实面积拉伸。

Fs=KAA02α2,σ=∂Fs∂A=KAαF_s=\frac{K_AA_0}{2}\alpha^2,\qquad \sigma=\frac{\partial F_s}{\partial A}=K_A\alpha
求导就是在同一材料状态附近比较一小段面积变化和能量变化。此式限于小应变的弹性拉伸模型。

完整例题:把能量换成热能

取教学参数 κ=20kBT\kappa=20k_BT,没有自发曲率。球面平均曲率能为 160πkBT≃503kBT160\pi k_BT\simeq503k_BT。在 T=298 KT=298\,\mathrm K 时,kBT≃4.11 pN nmk_BT\simeq4.11\,\mathrm{pN\,nm},总能约为 2.07×103 pN nm2.07\times10^3\,\mathrm{pN\,nm}。

把球半径从五微米增到十微米,上面这个结果保持不变。若膜从来是平的,要封成球仍须付出这份弯曲代价。球在已闭合、面积与体积固定的其他形状中是否最低,则是另一个比较问题。

再设 KA=0.20 N m−1K_A=0.20\,\mathrm{N\,m^{-1}}、α=0.01\alpha=0.01,小应变模型给出 σ=2.0 mN m−1\sigma=2.0\,\mathrm{mN\,m^{-1}}。这些数值只用来演算量纲,并未标定某种脂质。

把推理带进实验

把球面作为能量标尺

  1. 把约化体积调到一、自发曲率与极区曲率设为零,预测能量除以弯曲刚度后的数值。
  2. 比较数值积分与球面解析值,再降低约化体积,查看两个允许形状的能量。
  3. 保持几何约束,只改变自发曲率,解释局部偏好怎样改变形状排序。

练习与解释性答案

  1. 若球满足 c0R=2c_0R=2,它的平均曲率弯曲能是多少?总能一定为零吗?

    展开解释

    平均曲率项为零,因为整张球面恰好达到偏爱曲率。总能仍可能含高斯项、拉伸项和外力功;必须先说明采用的能量与参考零点。

  2. 没有自发曲率时,把同一形状整体缩小十倍,平均曲率弯曲能怎样变?

    展开解释

    保持不变。每处曲率增大十倍、平方增大一百倍,面积减到百分之一。这个缩放结论要求刚度不随位置和尺度改变,且连续膜模型仍适用。

  3. 为什么比较固定面积的两种形状时,不能靠张力乘以面积来选胜者?

    展开解释

    两者面积相同,这一项相同。要比较的是弯曲能和其他确实随形状改变的项。求解时的面积乘子仍决定驻点条件与力学响应。

来源与复算范围

  1. W. Helfrich, Elastic Properties of Lipid Bilayers: Theory and Possible Experiments, Z. Naturforsch. C 28, 693–703 (1973)

    曲率弹性原论文,693—703 页。课程采用局域二次曲率能,并另行声明球面为正的曲率约定;球面、膜管和小斜率例题在这里推导,没有复制原图或实测材料参数。

  2. U. Seifert, K. Berndl and R. Lipowsky, Shape transformations of vesicles, Physical Review A 44, 1182–1202 (1991)

    已发表论文:自发曲率模型、双层耦合模型与轴对称形状分支。配套实验在有限径向变分族中固定面积和体积寻优;完整相图、极窄颈极限和口形囊泡不属于这次复算。

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

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