A Wang-Landau study of a lattice model for lipid bilayer self-assembly

The Wang-Landau (WL) Monte Carlo method has been applied to simulate the self-assembly of a lipid bilayer on a 3D lattice. The WL method differs from conventional Monte Carlo methods in that a complete density of states is obtained directly for the system, from which properties, such as the free energy, can be derived. Furthermore, from a single WL simulation, continuous curves of the average energy and heat capacity can be determined, which provide a complete picture of the phase behavior. The lipid model studied consists of 3 or 5 coarse-grained segments on lattices of varying sizes, with the empty lattice sites representing water. A bilayer structure is found to form at low temperatures, with phase transitions to clusters as temperature increases. For 3-segment chains, varying lattice sizes were studied, with the observation that the ratio of chain number to lattice area (i.e., area per lipid) affects the phase transition temperature. At small ratios, only one phase transition occurs between the bilayer and cluster phases, while at high lipid ratios the phase transition occurs in a two-step process with a stable intermediate phase. This second phase transition was not observed in conventional Metropolis Monte Carlo simulations on the same model, demonstrating the advantage of being able to perform a complete scan of the whole temperature range with the WL method. For longer 5-segment chains similar phase transitions are also observed with changes in temperature. In the WL method, due to the extensive nature of the energy, the number of energy bins required to represent the density of states increases as the system size increases and so limits its practical application to larger systems. To improve this, an extension of the WL algorithm, the statistical-temperature Monte Carlo method that allows simulations with larger energy bin sizes, has recently been proposed and is implemented in this work for the 3-segment lattice model. The results obtained are in good agreement with the original WL method and appear to be independent of the energy bin size used.