Bayesian uncertainty quantification and propagation in molecular dynamics simulations: a high performance computing framework.
We present a Bayesian probabilistic framework for quantifying and propagating the uncertainties in the parameters of force fields employed in molecular dynamics (MD) simulations. We propose a highly parallel implementation of the transitional Markov chain Monte Carlo for populating the posterior probability distribution of the MD force-field parameters. Efficient scheduling algorithms are proposed to handle the MD model runs and to distribute the computations in clusters with heterogeneous architectures. Furthermore, adaptive surrogate models are proposed in order to reduce the computational cost associated with the large number of MD model runs. The effectiveness and computational efficiency of the proposed Bayesian framework is demonstrated in MD simulations of liquid and gaseous argon.
Computational Science and Engineering Laboratory, ETH Zürich, CH-8092 Zurich, Switzerland. panagiotis.angelikopoulos@mavt.ethz.ch
Panagiotis
P
Angelikopoulos
Costas
C
Papadimitriou
Petros
P
Koumoutsakos
D001499:Bayes Theorem; D002138:Calibration; D008390:Markov Chains; D056004:Molecular Dynamics Simulation; D009010:Monte Carlo Method; D035501:Uncertainty