---
title: "蛋白和脂质怎样一起改变膜形状"
description: "把组分作为可移动自由度，推导曲率耦合后的线性稳定性，复算有限波数失稳，并区分预设极区与相分离。"
canonical_url: "https://sandboxphysics.com/zh/stories/membrane-mechanics/curvature-and-composition/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/membrane-mechanics/curvature-and-composition/"
last_modified: "2026-10-28"
---

[膜与囊泡力学](https://sandboxphysics.com/zh/stories/membrane-mechanics/)

L3 · 第 7 章

# 蛋白和脂质怎样一起改变膜形状

Sandbox Physics · 2026-10-28

能让膜弯曲的分子若可以移动，就会寻找适合自己的弯曲区域；它们聚集后又改变膜偏爱的曲率。形状和组分因此互相影响。预先指定一块曲率区域能研究其中一半，描述相分离还需要让组分也参与寻优或运动。

## 开始前需要什么

第二章的局域弯曲能、第三章的约束寻优、第四章的模分解。下面引入小斜率曲面与二次稳定性；矩阵只用来整理两个耦合振幅。

- 建立形状与组分的耦合能量，说明预设分布和自由迁移的区别。
- 消去一个高度模，计算组分有效刚度并判断均匀态何时失稳。

## 把分子的分布也当成未知量

设 $\phi$ 表示某类分子的局部比例相对平均值的偏离，取无量纲。分子既可能偏爱某种曲率，也可能偏爱同类邻居。前者改变局域自发曲率，后者写入混合自由能。组分变化太急还要付出梯度能，它与两种组分交界处的代价有关。

一块预设极区只规定分子在哪里；真正的相分离则问它们为什么聚在那里。若总分子数不变，组分演化还必须保持总量。Jülicher 与 Lipowsky 研究了域边界和弯曲竞争导致的出芽；下面采用另一套明确给出的连续小形变教学模型来分析均匀态的起始失稳。

来源：F. Jülicher and R. Lipowsky · W. Helfrich

## 在近乎平的膜上写出最少的耦合

用平面坐标 $\mathbf x$ 和高度 $h(\mathbf x)$ 描述一小片膜，斜率很小。向上的法向与球面为正的约定给出 $2H\simeq-\nabla^2h$。定义有符号耦合 $c$，局域自发曲率取为 $c_0=-c\phi$。于是弯曲项可写成高度 Laplace 算子与组分的差的平方。

系数 $a$ 是均匀混合态附近的局域二阶刚度；它为负时混合本身偏爱分开。$b>0$ 抑制太细的组分起伏，$g>0$ 的四次项限制振幅无限长大。$\kappa$ 以能量计，$\sigma$ 与 $a,g$ 以能量除以面积计，$b$ 以能量计。

$$
F=\int\!\mathrm d^2x\left[\frac{\kappa}{2}(\nabla^2h-c\phi)^2+\frac{\sigma}{2}|\nabla h|^2+\frac a2\phi^2+\frac b2|\nabla\phi|^2+\frac g4\phi^4\right]
$$

这是课程自定的小斜率连续组分模型。讨论初始线性稳定性时只保留二次项；大幅出芽和最终域形状需要非线性曲面模型。

## 让高度先松弛，看看组分还剩多少代价

取同一波数 $q$ 的高度和组分余弦波。二次能量可以整理成一个对称矩阵。高度单独改变时的刚度来自弯曲与张力；交叉项告诉我们组分波能够用高度变化降低多少能量。

固定组分振幅，让高度取最低能量。把高度方程解出再代回，得到有效组分刚度 $S(q)$。若它为负，一个很小的该波数组分起伏就能降低总能，均匀态便在线性意义下不稳定。平均组分固定时，零波数本身不参与变化，但长波极限仍可作比较。

$$
H_{hh}=\kappa q^4+\sigma q^2,\quad H_{h\phi}=\kappa cq^2,\quad H_{\phi\phi}=a+bq^2+\kappa c^2
$$

模振幅按同一面积归一化。高度与组分的交叉项符号来自前述曲率与耦合约定。

$$
S(q)=H_{\phi\phi}-\frac{H_{h\phi}^2}{H_{hh}}=a+bq^2+\frac{\kappa c^2\sigma}{\kappa q^2+\sigma}\quad(q>0)
$$

减去的项是允许高度松弛释放的能量。膜若被钉死成平面，就不能释放这一部分。

## 完整例题：长波稳定，中间尺度先降低能量

选定长度与能量单位，把系数写为无量纲教学参数 $\kappa=20$、$\sigma=0.2$、$c=0.4$、$a=-1$、$b=0.1$。长波极限为 $S(0)=a+\kappa c^2=2.20$，为正。很高波数又会被正的梯度项抬高。最低点可能落在中间。

对 $q^2$ 求导后得到下式，代入可得 $q_*\simeq0.745$、$S(q_*)\simeq-0.888$，对应波长 $2\pi/q_*\simeq8.43$ 个所选长度单位。允许高度弯曲后，这段组分起伏能够降低能量。

这组数值可用二维二次能量直接对高度做数值最小化，再与消元公式比较。它检验的是线性失稳的起点；最终域的大小、形状和生长速度还取决于非线性项、守恒动力学与流体运动。增长最快的波数还取决于动力学系数和守恒条件。

$$
q_*^2=|c|\sqrt{\frac{\sigma}{b}}-\frac{\sigma}{\kappa}
$$

右端为正时给出内部极小值；否则最低点转到长波边界。数值和生成脚本保存在本路线的组分基准中。

## 从固定极区走向可移动的分子

形状实验中的 $c_0R_A=c_{\rm global}+c_{\rm polar}\cos^8\theta$ 是指定在两极的偏好。它能比较固定异质性怎样改变形状，却没有求解 $\phi$ 的迁移。要把它推进到相分离，需同时更新形状和组分，保持分子总量，并检验网格与界面宽度的收敛。

真实蛋白还可能有方向、占据面积和结合动力学。若这些量对问题关键，就应把它们作为独立模型变量；曲率参数描述的是已指定的局部弯曲偏好。

## 把推理带进实验

### [看固定极区怎样改变形状](https://sandboxphysics.com/zh/labs/vesicle-shapes/)

1. 先预测极区偏好曲率增加后，两个分支的能量会怎样变化。
2. 只改变极区曲率，比较轮廓和能量，保持面积、体积相同。
3. 解释这一步对应固定分布的哪部分物理，再为可迁移组分提出一个额外观测：例如局部浓度与曲率的时间相关。

## 练习与解释性答案

1. 把 $c$ 反号，线性有效刚度改变吗？局部高度关系改变吗？

   展开解释

   有效刚度只含耦合的平方，因此不变；最优高度与组分之间的符号关系会反转。能量稳定性相同，不代表膜向同一侧弯。

2. 为什么不能把最低刚度波数直接称为最终域大小？

   展开解释

   二次分析只判断小扰动能否降低能量。域长大后四次项、界面能、守恒动力学和曲面几何都参与；最终尺度需要求解并验证后续过程。

3. 若把膜完全钉平，上述教学参数还会在线性意义下失稳吗？

   展开解释

   钉平后组分刚度为 $a+bq^2+\kappa c^2=2.2+0.1q^2$，始终为正。这个例子中，允许高度松弛是有限波数失稳的关键。

[上一章：拉出一根膜管要多大的力](https://sandboxphysics.com/zh/stories/membrane-mechanics/membrane-tethers/) [下一章：同样抖得厉害，怎样看出主动驱动](https://sandboxphysics.com/zh/stories/membrane-mechanics/active-fluctuations/)

## 来源与复算范围

1. [W. Helfrich, Elastic Properties of Lipid Bilayers: Theory and Possible Experiments, Z. Naturforsch. C 28, 693–703 (1973)](https://doi.org/10.1515/znc-1973-11-1209)

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

2. [F. Jülicher and R. Lipowsky, Shape transformations of vesicles with intramembrane domains, Physical Review E 53, 2670–2683 (1996)](https://doi.org/10.1103/PhysRevE.53.2670)

   引言与域能量：弯曲和边界线能之间的竞争。本课程连续组分—高度矩阵及有限波数例题为自行推导的教学扩展，不是对论文域形状相图的复现。