Logo image
Back
Numerical Methods for Biomembranes Based on Piecewise Linear Surfaces
Book chapter   Peer reviewed

Numerical Methods for Biomembranes Based on Piecewise Linear Surfaces

John P. Brogan, Yilin Yang, Thomas P. -Y. Yu and Yonghong F Yang
Numerical Mathematics and Advanced Applications ENUMATH 2017, pp 809-817
05 Jan 2019
DOI: https://doi.org/10.1007/978-3-319-96415-7_76

Abstract

The shapes of phospholipid bilayer biomembranes are modeled by the celebrated Canham-Evans-Helfrich model as constrained Willmore minimizers. Several numerical treatments of the model have been proposed in the literature, one of which was used extensively by biophysicists over two decades ago to study real lipid bilayer membranes. While the key ingredients of this algorithm are implemented in Brakke’s well-known surface evolver software, some of its glory details were never explained by either the geometers who invented it or the biophysicists who used it. As such, most of the computational results claimed in the biophysics literature are difficult to reproduce. In this note, we give an exposition of this method, connect it with some related ideas in the literature, and propose a modification of the original method based on replacing mesh smoothing with harmonic energy regularization. We present a theoretical finding and related computational observations explaining why such a smoothing or regularization step is indispensable for the success of the algorithm. A software package called WMINCON is available for reproducing the experiments in this and related articles.

Metrics

17 Record Views
1 citations in Scopus

Details

Logo image