[Paper Review] Variational Problems in Elastic Theory of Biomembranes, Smectic-a Liquid Crystals, and Carbon Related Structures

Z. C. Tu, Z. C. Ou-Yang|ArXiv.org|Jun 21, 2005

Lipid Membrane Structure and Behavior17 references4 citations

TL;DR

This paper formulates variational problems in elastic theories of biomembranes, smectic-A liquid crystals, and carbon nanostructures using exterior differential forms and the moving frame method. It derives Euler-Lagrange equations for equilibrium shapes and computes second-order variations to analyze mechanical stability, showing that the membrane skeleton significantly enhances the critical pressure for spherical cell membranes, with a predicted critical pressure of ~2 Pa under physiological conditions.

ABSTRACT

After a brief introduction to several variational problems in the study of shapes of thin thickness structures, we deal with variational problems on 2-dimensional surface in 3-dimensional Euclidian space by using exterior differential forms. The morphological problems of lipid bilayers and stabilities of cell membranes are also discussed. The key point is that the first and the second order variations of the free energy determine equilibrium shapes and mechanical stabilities of structures.

Motivation & Objective

• To develop a unified variational framework for thin elastic structures including biomembranes, smectic-A liquid crystals, and carbon-based nanostructures.
• To analyze equilibrium shapes of lipid bilayers and cell membranes by minimizing free energy functionals incorporating curvature and surface energy.
• To investigate mechanical stability of spherical cell membranes under internal pressure, particularly the role of the membrane skeleton.
• To derive and solve boundary conditions for open lipid bilayers with free edges using differential forms.
• To explain the formation of focal conic domains in smectic-A liquid crystals via curvature elastic energy minimization.

Proposed method

• Uses exterior differential forms and the moving frame method to derive Euler-Lagrange equations on 2D surfaces embedded in R³.
• Applies the variational principle to free energy functionals involving mean curvature (H), Gaussian curvature (K), and thickness (t) for different systems.
• Derives shape equations and boundary conditions for open lipid bilayers, including line tension and geodesic curvature effects.
• Introduces a modified free energy functional for cell membranes that includes in-plane strain contributions from the membrane skeleton.
• Employs Hodge decomposition to express perturbations in terms of scalar functions and computes second-order variations of the free energy.
• Uses spherical harmonic expansions to analyze the stability of spherical membranes and determine the critical pressure for buckling.

Experimental results

Research questions

• RQ1 What are the equilibrium shapes of lipid bilayers and cell membranes under pressure and curvature elasticity?
• RQ2 How does the membrane skeleton influence the mechanical stability of spherical cell membranes?
• RQ3 Why do focal conic structures form in smectic-A liquid crystals despite the preference for flat layers?
• RQ4 What are the shape equations and boundary conditions for open lipid bilayers with free edges?
• RQ5 How do curvature elastic energy and surface energy balance to determine the morphology of smectic-A domains?

Key findings

• The membrane skeleton increases the critical pressure for mechanical instability of spherical cell membranes from ~0.2 Pa (without skeleton) to ~2 Pa (with skeleton).
• The critical pressure is determined by minimizing over l-mode instabilities, yielding p_c ≈ 2 Pa for typical values: k_c ≈ 20k_B T, k_d ≈ 6×10⁻⁴ k_B T/nm², R ≈ 1 μm.
• The second-order variation of the free energy, δ²F = G₁ + G₂, is positive definite when p < p_l, ensuring mechanical stability against shape fluctuations.
• The stability condition depends on the competition between curvature elasticity (k_c), membrane skeleton stiffness (k_d), and pressure (p), with distinct regimes depending on k_d R² relative to 121k_c.
• For open lipid bilayers, the derived boundary conditions (Eqs. 4–6) include contributions from normal curvature, geodesic curvature, and torsion, ensuring mechanical equilibrium at the edge.
• The model explains focal conic domain formation in smectic-A LCs as a balance between curvature elastic energy and the Gibbs free energy difference between isotropic and SmA phases.