High-productivity, high-performance workflow for virus-scale electrostatic simulations with Bempp-Exafmm

03/01/2021
by   Tingyu Wang, et al.
0

Biomolecular electrostatics is key in protein function and the chemical processes affecting it.Implicit-solvent models expressed by the Poisson-Boltzmann (PB) equation can provide insights with less computational power than full atomistic models, making large-system studies – at the scale of viruses, for example – accessible to more researchers. This paper presents a high-productivity and high-performance computational workflow combining Exafmm, a fast multipole method (FMM) library, and Bempp, a Galerkin boundary element method (BEM) package. It integrates an easy-to-use Python interface with well-optimized computational kernels that are written in compiled languages. Researchers can run PB simulations interactively via Jupyter notebooks, enabling faster prototyping and analyzing. We provide results that showcase the capability of the software, confirm correctness, and evaluate its performance with problem sizes between 8,000 and 2 million boundary elements. A study comparing two variants of the boundary integral formulation in regards to algebraic conditioning showcases the power of this interactive computing platform to give useful answers with just a few lines of code. As a form of solution verification, mesh refinement studies with a spherical geometry as well as with a real biological structure (5PTI) confirm convergence at the expected 1/N rate, for N boundary elements. Performance results include timings, breakdowns, and computational complexity. Exafmm offers evaluation speeds of just a few seconds for tens of millions of points, and 𝒪(N) scaling. This allowed computing the solvation free energy of a Zika virus, represented by 1.6 million atoms and 10 million boundary elements, at 80-min runtime on a single compute node (dual 20-core Intel Xeon Gold 6148). All results in the paper are presented with utmost care for reproducibility.

READ FULL TEXT

page 10

page 11

page 12

research
08/24/2021

Towards optimal boundary integral formulations of the Poisson-Boltzmann equation for molecular electrostatics

The Poisson-Boltzmann equation offers an efficient way to study electros...
research
09/20/2020

Efficient mesh refinement for the Poisson-Boltzmann equation with boundary elements

The Poisson-Boltzmann equation is a widely used model to study the elect...
research
08/31/2021

An Efficient Finite Element Iterative Method for Solving a Nonuniform Size Modified Poisson-Boltzmann Ion Channel Model

In this paper, a nonuniform size modified Poisson-Boltzmann ion channel ...
research
09/19/2023

Julia as a unifying end-to-end workflow language on the Frontier exascale system

We evaluate using Julia as a single language and ecosystem paradigm powe...
research
05/03/2021

[Re] Three-dimensional wake topology and propulsive performance of low-aspect-ratio pitching-rolling plates

This article reports on a full replication study in computational fluid ...
research
04/27/2021

Including monopoles to a fully desingularized boundary element method for acoustics

The inclusion of domain (point) sources into a three dimensional boundar...

Please sign up or login with your details

Forgot password? Click here to reset