Abstract
Nonequilibrium simulation protocols based on shear deformations are applied to determine the surface viscosity and interleaflet friction of lipid bilayers. At high shear rates, a non-Newtonian shear thinning regime is observed, but lower shear rates yield a Newtonian plateau and results that are consistent with equilibrium measurements based on fluctuation-dissipation theorems. Application to all-atom bilayers modeled with the CHARMM36 parameter set yields values for the surface viscosity that are consistent with microscopic measurements based on membrane protein diffusion but are approximately 10 times lower than more macroscopic experimental measurements. The interleaflet friction is about 10 times lower than experimental measurements. Trends across different lipids, temperatures, and ternary liquid-disordered phase mixtures produce results that are consistent with experimental diffusion constants. Application of the protocol to the liquid-ordered phase fails to yield a Newtonian plateau, suggesting more complex rheology.