In this contribution to the proceedings of the 11th Mathematical Society of Japan (MSJ) Seasonal Institute (July 2018) we give an overview of some recent work on a mathematical model for small deformations of a spherical membrane. The idea is to consider perturbations to minimisers of a surface geometric energy. The model is obtained from consideration of second order approximations to a perturbed energy. In particular, the considered problems involve particle constraints and surface phase field energies.
We derive and analyse an energy to model lipid raft formation on biological membranes involving a coupling between the local mean curvature and the local composition. We apply a perturbation method recently introduced by Fritz, Hobbs and the first author to describe the geometry of the surface as a graph over an undeformed Helfrich energy minimising surface. The result is a surface Cahn-Hilliard functional coupled with a small deformation energy We show that suitable minimisers of this energy exist and consider a gradient flow with conserved Allen-Cahn dynamics, for which existence and uniqueness results are proven. Finally, numerical simulations show that for the long time behaviour raft-like structures can emerge and stablise, and their parameter dependence is further explored.
In this paper we introduce a mathematical model for small deformations induced by external forces of closed surfaces that are minimisers of Helfrich-type energies. Our model is suitable for the study of deformations of cell membranes induced by the cytoskeleton. We describe the deformation of the surface as a graph over the undeformed surface. A new Lagrangian and the associated Euler-Lagrange equations for the height function of the graph are derived. This is the natural generalisation of the well known linearisation in the Monge gauge for initially flat surfaces. We discuss energy perturbations of point constraints and point forces acting on the surface. We establish existence and uniqueness results for weak solutions on spheres and on tori. Algorithms for the computation of numerical solutions in the general setting are provided. We present numerical examples which highlight the behaviour of the surface deformations in different settings at the end of the paper.
We consider a model of a biomembrane with attached proteins. The membrane is represented by a near spherical continuous surface and attached proteins are described as discrete rigid structures which attach to the membrane at a finite number of points. The resulting surface minimises a quadratic elastic energy (obtained by a perturbation of the Canham-Helfrich energy) subject to the point constraints which are imposed by the attachment of the proteins. We calculate the derivative of the energy with respect to protein configurations. The proteins are constrained to move tangentially by translation and by rotation in the axis normal to a reference point. Previous studies have typically restricted themselves to a nearly flat membrane and circular inclusions. A numerically accessible representation of this derivative is derived and employed in some numerical experiments.
We consider sharp interface asymptotics for a phase field model of two phase near spherical biomembranes involving a coupling between the local mean curvature and the local composition proposed by the first and second authors. The model is motivated by lipid raft formation. We introduce a reduced diffuse interface energy depending only on the membrane composition and derive the $Gamma-$limit. We demonstrate that the Euler-Lagrange equations for the limiting functional and the sharp interface energy coincide. Finally, we consider a system of gradient flow equations with conserved Allen-Cahn dynamics for the phase field model. Performing a formal asymptotic analysis we obtain a system of gradient flow equations for the sharp interface energy coupling geodesic curvature flow for the phase interface to a fourth order PDE free boundary problem for the surface deformation.
A Rayleigh-Schrodinger type of perturbation scheme is employed to study weak self-interacting scalar potential perturbations occurring in scalar field models describing 1D domain kinks and 3D domain walls. The solutions for the unperturbed defects are modified by the perturbing potentials. An illustration is provided by adding a cubic potential to the familiar quartic kink potential and solving for the first order correction to the kink solution, using a slab approximation. A result is the appearance of an asymmetric scalar potential with different, nondegenerate, vacuum values and the subsequent formation of vacuum bubbles.