Theory Reference — GUV Fluctuation Suite

01Helfrich Free Energy

The elastic free energy of a lipid bilayer is described by the Helfrich functional. The local energy density depends on curvature and tension:

e = \frac{1}{2}\kappa(2H - C_0)^2 + \sigma

Integrating over the vesicle surface gives:

F_\text{Hel} = \int_S \frac{\kappa}{2}(2H - C_0)^2 \, dA + \sigma A

κ	Bending rigidity	units: kBT ≈ 10⁻²¹ J
σ	Membrane surface tension	units: kBT/µm²
C₀	Spontaneous curvature	units: µm⁻¹
H	Mean curvature = (c₁+c₂)/2	c₁,c₂ = principal curvatures

        Physical meaning: κ penalises bending. σ penalises area change.
        C₀ ≠ 0 means the membrane has an intrinsic preferred curvature (asymmetry between leaflets).
      

02Quasi-Spherical Vesicles

GUVs are nearly spherical. We describe small deviations from a perfect sphere of radius R₀:

r(\theta, \phi) = R_0 \left[ 1 + u(\theta, \phi) \right], \quad |u| \ll 1

The fluctuation field u(θ,φ) can be expanded in spherical harmonics Y_l^m(θ,φ), which are the natural eigenmodes of shape deformations on a sphere:

u(\theta, \phi) = \sum_{l=0}^{\infty} \sum_{m=-l}^{l} u_{lm} \, Y_l^m(\theta, \phi)

        Physical modes:
        l=0 → volume change (suppressed by incompressibility) ·
        l=1 → rigid translation (zero energy) ·
        l≥2 → genuine shape fluctuations (physical, analysed here)
      

03Fluctuation Spectrum & Equipartition

Inserting the spherical harmonic expansion into the Helfrich energy and retaining quadratic terms yields a sum of independent harmonic oscillators:

F = \sum_{l,m} \frac{1}{2} \alpha_l |u_{lm}|^2

where the mode stiffness is:

\alpha_l = \kappa(l-1)(l+2)\left[l(l+1) + \Sigma\right] \] \[ \Sigma = \frac{\sigma R^2}{\kappa} + \frac{R^2 C_0^2}{2} + 2RC_0

The equipartition theorem assigns ½kBT of energy per quadratic degree of freedom:

\frac{1}{2} \alpha_l \langle|u_{lm}|^2\rangle = \frac{1}{2} k_B T \] \[ \Rightarrow \quad \langle|u_{lm}|^2\rangle = \frac{k_B T}{\alpha_l} = \frac{k_B T}{\kappa(l-1)(l+2)\left[l(l+1)+\Sigma\right]}

        Key insight: Stiffer modes (larger α_l) have smaller amplitudes.
        Low-l modes probe surface tension Σ; high-l modes probe bending rigidity κ.
      

04Angular Auto-Correlation Function

In experiments we image the vesicle equator. The equatorial contour gives r_eq(φ,t) = R[1 + u(π/2, φ, t)]. The instantaneous angular average is:

r_\text{avg}(t) = \frac{1}{2\pi}\int_0^{2\pi} r_{eq}(\phi, t)\, d\phi

The equatorial angular auto-correlation function quantifies how the contour shape at angle φ correlates with the shape at φ+γ:

\xi_{eq}(\gamma, t) = \frac{1}{2\pi r_\text{avg}^2} \int_0^{2\pi} \left[ r_{eq}(\phi+\gamma, t)\, r_{eq}(\phi, t) - r_\text{avg}^2 \right] d\phi

For small fluctuations (|u| ≪ 1), using Taylor expansion of the normalisation factor:

\xi_{eq}(\gamma, t) \approx \frac{1}{2\pi} \int_0^{2\pi} u(\phi+\gamma, t)\, u(\phi, t)\, d\phi

        Interpretation:
        ξ(0) > 0 always (self-correlation) ·
        ξ(π/2) < 0 indicates anti-correlation (quadrupolar mode dominates) ·
        Double-bell shape ∝ P₂(cos γ) is the signature of l=2 mode dominance.
      

05Legendre Expansion

The time-averaged theoretical autocorrelation expands naturally as:

\langle \xi^{th}(\gamma) \rangle = \sum_{l \geq 2} b_l \, P_l(\cos\gamma)

where the Legendre coefficients are:

b_l = \frac{2l+1}{4\pi} \langle|u_{lm}|^2\rangle = \frac{2l+1}{4\pi} \cdot \frac{k_B T}{\kappa(l-1)(l+2)[l(l+1)+\Sigma]}

The experimental ξ^exp(γ) is fitted to this expression. The fit parameters are: j = kBT/κ and s = Σ.

Double-Bell Structure (l = 2)

The l=2 (quadrupolar) mode gives P₂(cos γ) = ½(3cos²γ − 1). This has:

\text{maxima at } \gamma = 0, \pi \quad (+1) \qquad \text{minima at } \gamma = \tfrac{\pi}{2}, \tfrac{3\pi}{2} \quad (-\tfrac{1}{2})

This characteristic double-bell shape in ξ(γ) is the signature of an ellipsoidal (quadrupolar) deformation mode. Points at γ=0 and γ=π are along the elongation axis (positively correlated), while points at γ=π/2 are perpendicular (anti-correlated).

06Parameter Extraction (χ² Method)

The theoretical mean Legendre coefficient is:

\bar{B}_l(j, s) = \frac{j}{p_l + s\, q_l} \] \[ p_l = \frac{4\pi(l-1)(l+2)(l+1)l}{2l+1} \qquad q_l = \frac{4\pi(l-1)(l+2)}{2l+1} \] \[ j = \frac{k_B T}{\kappa} \qquad s = \Sigma

Define the goodness-of-fit as:

\chi^2(j, s) = \sum_{l=l_{min}}^{l_{max}} \left[ \frac{B_l - \bar{B}_l(j,s)}{\sigma_l} \right]^2

Setting ∂χ²/∂j = 0 gives j as an explicit function of s:

j(s) = \frac{\displaystyle\sum_l \frac{p_l B_l}{(p_l \sigma_l)^2} \left(1 + s\frac{q_l}{p_l}\right)} {\displaystyle\sum_l \frac{1}{(p_l \sigma_l)^2} \left(1 + s\frac{q_l}{p_l}\right)^2}

Substituting j(s) back into χ²(j(s), s) gives a 1D function of s only. Its minimum is found numerically (golden-section search). The optimal s = Σ is then inserted into j(s) to get j, and hence κ = kBT/j.

07Error Estimation

Near the χ² minimum, a Gaussian approximation gives the covariance matrix:

H' = \begin{pmatrix} \partial^2\chi^2/\partial j^2 & \partial^2\chi^2/\partial j \partial s \\ \partial^2\chi^2/\partial s\partial j & \partial^2\chi^2/\partial s^2 \end{pmatrix}^{-1}_\text{min}

The variances of the fitted parameters are:

\sigma^2(j) = H'_{11} \qquad \sigma^2(s) = H'_{22} \] \[ \sigma^2(\kappa) = H'_{11} \cdot \frac{\kappa^2}{j^2} \cdot k_B^2 T^2

08References

[1] Brochard & Lennon, J. de Physique 36, 1035–1047 (1975) — Frequency spectrum of flicker phenomenon in erythrocytes.

[2] Helfrich, Z. Naturforsch. C 28, 693–703 (1973) — Elastic properties of lipid bilayers.

[3] Kumar et al., Chem. Phys. Lipids 259, 105374 (2024) — Bottom-up approach to explore α-amylase assisted membrane remodelling.

[4] Dimova & Marques, The Giant Vesicle Book, CRC Press (2019).

[5] Faucon et al., J. de Physique 50, 2389–2414 (1989) — Bending elasticity and thermal fluctuations of lipid membranes.

From Flicker to Physics