UCL Discovery
UCL home » Library Services » Electronic resources » UCL Discovery

Large-scale parallelised boundary element method electrostatics for biomolecular simulation

Fallaize, D.R.; (2012) Large-scale parallelised boundary element method electrostatics for biomolecular simulation. Doctoral thesis , UCL (University College London). Green open access


Download (73MB)


Large-scale biomolecular simulations require a model of particle interactions capable of incorporating the behaviour of large numbers of particles over relatively long timescales. If water is modelled as a continuous medium then the most important intermolecular forces between biomolecules can be modelled as long-range electrostatics governed by the Poisson- Boltzmann Equation (PBE). We present a linearised PBE solver called the "Boundary Element Electrostatics Program"(BEEP). BEEP is based on the Boundary Element Method (BEM), in combination with a recently developed O(N) Fast Multipole Method (FMM) algorithm which approximates the far-�field integrals within the BEM, yielding a method which scales linearly with the number of particles. BEEP improves on existing methods by parallelising the underlying algorithms for use on modern cluster architectures, as well as taking advantage of recent progress in the �field of GPGPU (General Purpose GPU) Programming, to exploit the highly parallel nature of graphics cards. We found the stability and numerical accuracy of the BEM/FMM method to be highly dependent on the choice of surface representation and integration method. For real proteins we demonstrate the critical level of surface detail required to produce converged electrostatic solvation energies, and introduce a curved surface representation based on Point-Normal G1-continuous triangles which we �find generally improves numerical stability compared to a simpler surface constructed from planar triangles. Despite our improvements upon existing BEM methods, we �find that it is not possible to directly integrate BEM surface solutions to obtain intermolecular electrostatic forces. It is, however, practicable to use the total electrostatic solvation energy calculated by BEEP to drive a Monte-Carlo simulation.

Type: Thesis (Doctoral)
Title: Large-scale parallelised boundary element method electrostatics for biomolecular simulation
Open access status: An open access version is available from UCL Discovery
Language: English
UCL classification: UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Engineering Science > Dept of Computer Science > CoMPLEX: Mat&Phys in Life Sci and Exp Bio
URI: https://discovery.ucl.ac.uk/id/eprint/1343954
Downloads since deposit
Download activity - last month
Download activity - last 12 months
Downloads by country - last 12 months

Archive Staff Only

View Item View Item