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

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

设 $\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$ 以能量计。

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

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

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

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

选定长度与能量单位，把系数写为无量纲教学参数 $\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$ 个所选长度单位。允许高度弯曲后，这段组分起伏能够降低能量。

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

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

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

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