Simulation of multiscale hydrophobic lipid dynamics via efficient integral equation methods
Document Type
Article
Publication Date
1-1-2020
Abstract
In this paper, a mathematical model for long-range, hydrophobic attraction between amphiphilic particles is developed to quantify the macroscopic assembly and mechanics of a lipid bilayer membrane in solvents. The nonlocal interactions between amphiphilic particles are obtained from the first domain variation of a hydrophobicity functional, giving rise to forces and torques (between particles) that dictate the motion of both particles and the surrounding solvent. The functional minimizer (that accounts for hydrophobicity at molecular-aqueous interfaces) is a solution to a boundary value problem of the screened Laplace equation. We reformulate the boundary value problem as a second-kind integral equation (SKIE), discretize the SKIE using a Nyström discretization and Quadrature by Expansion (QBX), and solve the resulting linear system iteratively using GMRES. We evaluate the required layer potentials using the GIGAQBX fast algorithm, a variant of the Fast Multipole Method (FMM), yielding the required particle interactions with asymptotically optimal cost. A mobility problem formulation supplies the motion for the rigid particles in a viscous fluid. The simulated fluid-particle systems exhibit a variety of multiscale behaviors over both time and length. Over short time scales, the numerical results show self-assembly for model lipid particles. For large system simulations, the particles form realistic configurations like micelles and bilayers. Over long time scales, the bilayer shapes emerging from the simulation appear to minimize a form of bending energy.
Identifier
85083900065 (Scopus)
Publication Title
Multiscale Modeling and Simulation
External Full Text Location
https://doi.org/10.1137/18M1219503
e-ISSN
15403467
ISSN
15403459
First Page
79
Last Page
103
Issue
1
Volume
18
Grant
DMS-1418918
Fund Ref
National Science Foundation
Recommended Citation
Fu, Szu Pei P.; Ryham, Rolf; Klöckner, Andreas; Wala, Matt; Jiang, Shidong; and Young, Yuan Nan, "Simulation of multiscale hydrophobic lipid dynamics via efficient integral equation methods" (2020). Faculty Publications. 5761.
https://digitalcommons.njit.edu/fac_pubs/5761
