Bayesian Uncertainty Quantification for Particle-Based Simulation of Lipid Bilayer Membranes

Abstract

A number of problems of interest in applied mathematics and biology involve the quantification of uncertainty in computational and real-world models. A recent approach to Bayesian uncertainty quantification using transitional Markov chain Monte Carlo (TMCMC) is extremely parallelizable and has opened the door to a variety of applications which were previously too computationally intensive to be practical. In this chapter, we first explore the machinery required to understand and implement Bayesian uncertainty quantification using TMCMC. We then describe dissipative particle dynamics, a computational particle simulation method which is suitable for modeling biological structures on the subcellular level, and develop an example simulation of a lipid membrane in fluid. Finally, we apply the algorithm to a basic model of uncertainty in our lipid simulation, effectively recovering a target set of parameters (along with distributions corresponding to the uncertainty) and demonstrating the practicality of Bayesian uncertainty quantification for complex particle simulations.