--- title: "Reading bending rigidity from thermal flicker" description: "Derive quasi-spherical mode variance from equipartition, separate rigidity and tension, and project spherical modes onto an equatorial contour." canonical_url: "https://sandboxphysics.com/stories/membrane-mechanics/thermal-modes/" language: "en" page_type: "story" alternate_languages: zh-CN: "https://sandboxphysics.com/zh/stories/membrane-mechanics/thermal-modes/" last_modified: "2026-10-28" --- [Membrane mechanics](https://sandboxphysics.com/stories/membrane-mechanics/) L2 · Chapter 4 # Reading bending rigidity from thermal flicker Sandbox Physics · 2026-10-28 A membrane that bends easily fluctuates more under thermal motion. Measuring several length scales can separate bending rigidity from tension: short wavelengths pay strongly for bending, and long wavelengths are more sensitive to tension. Temperature supplies the energy scale missing from a static shape. ## Before you begin Curvature and bending energy from Chapters 1–2. Variance, independent modes and equipartition are introduced here. Think of spherical harmonics initially as distinguishable ripple patterns on a sphere. - Derive the variance of an independent quadratic mode and explain how scale separates tension from rigidity. - Explain why a spherical mode and a Fourier mode of the equatorial contour are different objects. ## Decompose a fluctuating surface into ripples For a nearly spherical membrane of radius $R$, divide radial displacement by radius to obtain a small dimensionless deformation $u(\theta,\phi)$. Spherical harmonics are standard ripple patterns, analogous to sines on a circle. Increasing $\ell$ gives finer surface structure; $m$ distinguishes variation around the axis. The zero mode mainly changes mean radius and first-order modes translate the vesicle. After centering, shape fluctuations begin at $\ell=2$. The quadratic theory below assumes zero spontaneous curvature, small deformation and uniform tension. Our fluctuation model specifies tension instead of enforcing exactly fixed excess area in every frame. $$ r(\theta,\phi)=R[1+u(\theta,\phi)],\qquad u=\sum_{\ell,m}u_{\ell m}Y_{\ell m} $$ Harmonics are normalized over solid angle. A second-order mean-radius correction can maintain volume; large deformations require the full surface energy. Source: B. Loubet ## Thermal energy in one quadratic degree of freedom Take one real mode amplitude $q$ with energy $K_q q^2/2$. Equilibrium probability is proportional to $\exp[-K_q q^2/(2k_BT)]$. Comparing this with a centered Gaussian gives variance $k_BT/K_q$. Variance is the long-time mean squared amplitude. Expand the Helfrich energy around a sphere. Orthogonality removes cross terms between different harmonics, leaving independent quadratic degrees of freedom. Curvature penalizes high-order ripples especially strongly. Tension adds a cost to surface undulations. $$ E_\ell=(\ell-1)(\ell+2)[\ell(\ell+1)+\bar\sigma],\qquad \bar\sigma=\frac{\sigma R^2}{\kappa} $$ The bar marks dimensionless tension. Physical tension still has units of force per length. $$ 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)}} $$ The index enumerates normalized real harmonic degrees of freedom. Each retained quadratic coefficient must be positive for a finite equilibrium variance. Source: B. Loubet ## Worked example at two tensions Set $\kappa=20k_BT$. For $\ell=2$ and zero dimensionless tension, $E_2=24$. One normalized real coefficient has variance $1/480\simeq0.00208$ and root-mean-square amplitude $0.0456$. This is a harmonic coefficient; local radial displacement also includes the radius and the local harmonic value. Raise dimensionless tension to $\bar\sigma=6$. Then $E_2=48$ and the variance halves. For $\ell=10$, the bracket increases only from $110$ to $116$. Measuring several scales helps identify which suppression comes from tension. Over a narrow mode range, rigidity and tension can compensate. More frames reduce random uncertainty; separating the parameters also requires a sufficiently broad scale range. ## One equatorial wave contains several spherical modes Tracking an equatorial contour means taking $\theta=\pi/2$ and Fourier-decomposing along azimuth. A single circular harmonic can receive contributions from several spherical orders. Treating it as one spherical harmonic changes the fitted model. Write its normalized cosine amplitude as $A_m=\sum_{\ell\geq m}w_{\ell m}q_{\ell m}$, where $w_{\ell m}=\sqrt2N_{\ell m}P_\ell^m(0)$. Independent mode variances add to give the spectrum below. Sine amplitudes have the same variance. A finite-order model sums only its retained terms. An objective also collects information over a finite focal depth, so real optical observation need not equal an infinitesimally thin equatorial section. We separate this issue from angular pixel averaging and exposure. The next chapter follows the measurement from membrane to image to estimate. $$ \langle A_m^2\rangle=\sum_{\ell\geq\max(2,m)}w_{\ell m}^2\frac{k_BT}{\kappa E_\ell} $$ Weights are the harmonic values on the equator. Parity makes some vanish. Surface order, contour mode number and observation geometry play distinct roles. Source: S. A. Rautu et al. · J. Pécréaux et al. ## Try the reasoning ### [Use several fluctuation scales](https://sandboxphysics.com/labs/membrane-fluctuations/) 1. Predict how rigidity changes the full flicker amplitude and which part of the spectrum responds most to tension. 2. Change one parameter at a time. Compare the ideal contour with the spectrum derived from saved observations. 3. As an advanced offline exercise, export the record and compare a low-mode-only fit with the full spectrum. Explain why a narrow scale range couples the parameters more strongly. ## Exercises and explanations 1. At fixed temperature and dimensionless tension, double rigidity. What happens to variance and root-mean-square amplitude? Read the explanation Variance halves and root-mean-square amplitude falls by the square root of two. Keeping dimensionless tension fixed requires adjusting physical tension with rigidity; the two fixed conditions differ. 2. A first-order harmonic drifts. Does this indicate an exceptionally soft shape mode? Read the explanation First-order harmonics mainly translate the vesicle. Track and correct the center before interpreting shape fluctuations. Center-localization error can also contaminate low contour modes. 3. A high-mode spectrum flattens. Must the membrane have become softer? Read the explanation First examine localization noise, pixels and focal depth. A noise floor can flatten a descending physical spectrum. The next chapter includes camera effects in the observation model. [Previous: Several shapes for the same membrane](https://sandboxphysics.com/stories/membrane-mechanics/shape-minima/) [Next: Why a camera can make a membrane look stiffer](https://sandboxphysics.com/stories/membrane-mechanics/camera-and-inference/) ## Sources and calculated scope 1. [J. Pécréaux et al., Refined contour analysis of giant unilamellar vesicles, European Physical Journal E 13, 277–290 (2004)](https://doi.org/10.1140/epje/i2004-10001-9) Contour fluctuation spectroscopy, camera integration and the observation plane. The course constructs its own seeded finite-mode records and derives its own exposure integral. Numerical fits here concern those synthetic records, not the paper’s lipid measurements. 2. [S. A. Rautu et al., The role of optical projection in the analysis of membrane fluctuations, Soft Matter 13, 3480–3483 (2017)](https://doi.org/10.1039/C7SM00108H) Equations (1)–(3), Figure 1 and the finite focal-depth measurement problem. The ideal equatorial slice in our model has zero focal depth. Finite angular averaging and exposure are modeled separately; the full optical projection operator is discussed as an extension. 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)](https://arxiv.org/html/1112.2919v2) Version 2, 14 December 2011: quadratic quasi-spherical modes and a direct active force with exponential temporal correlations. Our equal-variance response comparison is an explicitly specified finite linear model. We do not reproduce the paper’s self-consistent excess-area tension curves or permeation and curvature-force cases.