The lipid bilayer at the mesoscale : a physical continuum model

Abstract

We study a continuum model of the lipid bilayer based on minimizing the free energy of a mixture of water and lipid molecules. This paper extends previous work by Blom and Peletier (2004) in the following ways. (a) It formulates a more physical model of the hydrophobic effect to facilitate connections with microscale simulations. (b) It clarifies the meaning of the model parameters. (c) It outlines a method for determining parameter values so that physically-realistic bilayer density profiles can be obtained, for example for use in macroscale simulations. Points (a)-(c) suggest that the model has potential to robustly connect some micro- and macroscale levels of multiscale blood fl.ow simulations. The mathematical modelling in point (a) is based upon a consideration of the underlying physics of intermolecular forces. The governing equations thus obtained are minimized by gradient flows via a novel numerical approach; this enables point (b). The numerical results are shown to behave physically in terms of the effect of background concentration, in contrast to the earlier model which is shown here to not display the expected behaviour. A "short-tail" approximation of the lipid molecules also gives an analytical tool which yields critical values of some parameters under certain conditions. Point (c) involves the first quantitative comparison of the numerical data with physical experimental results.