Numerical analysis of viscoelasticity of two-dimensional fluid membranes under oscillatory loadings

TL;DR

This study models the viscoelastic behavior of two-dimensional biomembranes under oscillatory loads, revealing frequency-dependent transitions from elastic to viscous dominance, which enhances understanding of membrane mechanics.

Contribution

It introduces a numerical analysis of membrane viscoelasticity using hydrodynamic equations derived from Onsager's principle, focusing on frequency effects and state transitions.

Findings

Membrane behavior shifts from elastic to viscous with increasing frequency.

Transitions occur between 10 to 100 Hz, largely independent of surface tension.

Eigenvalue analysis explains the relaxation mechanisms behind state changes.

Abstract

Biomembranes consisting of two opposing phospholipid monolayers, which comprise the so-called lipid bilayer, are largely responsible for the dual solid-fluid behavior of individual cells and viruses. Quantifying the mechanical characteristics of biomembrane, including the dynamics of their in-plane fluidity, can provide insight not only into active or passive cell behaviors but also into vesicle design for drug delivery systems. Despite numerous studies on the mechanics of biomembranes, their dynamical viscoelastic properties have not yet been fully described. We thus quantify their viscoelasticity based on a two-dimensional (2D) fluid membrane model, and investigate this viscoelasticity under small amplitude oscillatory loadings in micron-scale membrane area. We use hydrodynamic equations of bilayer membranes, obtained by Onsager's variational principle, wherein the fluid membrane is assumed to be an almost planar bilayer membrane. Simulations are performed for a wide range of oscillatory frequencies $f$ and membrane tensions. Our numerical results show that as frequencies increase, membrane characteristics shift from an elastic-dominant to viscous-dominant state. Furthermore, such state transitions obtained with a 1-$\mu$m-wide loading profile appear with frequencies between $O(f) = 10^1-10^2$ Hz, and almost independently of surface tensions. We discuss the formation mechanism of the viscous- or elastic-dominant transition based on relaxation rates that correspond to the eigenvalues of the dynamical matrix in the governing equations.