Premium
Determining the Local Shear Viscosity of a Lipid Bilayer System by Reverse Non‐Equilibrium Molecular Dynamics Simulations
Author(s) -
Müller Thomas J.,
MüllerPlathe Florian
Publication year - 2009
Publication title -
chemphyschem
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.016
H-Index - 140
eISSN - 1439-7641
pISSN - 1439-4235
DOI - 10.1002/cphc.200900156
Subject(s) - bilayer , dipalmitoylphosphatidylcholine , molecular dynamics , lipid bilayer , chemistry , lipid bilayer mechanics , chemical physics , viscosity , thermodynamics , slip (aerodynamics) , shear (geology) , lipid bilayer phase behavior , materials science , computational chemistry , physics , membrane , phospholipid , biochemistry , phosphatidylcholine , composite material
The parallel shear viscosity of a dipalmitoylphosphatidylcholine (DPPC) bilayer system is studied by reverse non‐equilibrium molecular dynamics simulations (RNEMD) with two different united‐atom force fields. The results are related to diffusion coefficients and structural distributions obtained by equilibrium molecular simulations. We investigate technical issues of the algorithm in the bilayer setup, namely, the dependence of the velocity profiles on the imposed flux and the influence of the thermostat on the calculated shear viscosity. We introduce the concept of local shear viscosity and investigate its dependence on the slip velocity of the monolayers and the particle density at the headgroup–water interface and the tail–tail interface. With this we demonstrate that the lipid bilayer is more viscous than the surrounding water phase, and that slip takes place near the headgroup region and in the centre of the bilayer where the alkyl tails meet. We also quantify the apparent increase in viscosity of the water molecules entangled at the water–headgroup interface.