2019/11/04 by Giberti, Federico, Cheng, Bingqing, Tribello, Gareth Aneurin +1 · 1 citation
#Computational Physics (physics.comp-ph) #FOS: Physical sciences
paper · doi:10.48550/arxiv.1911.01140
Atomistic modelling of phase transitions, chemical reactions, or other rare events that involve overcoming high free energy barriers usually entails prohibitively long simulation times. Introducing a bias potential as a function of an appropriately-chosen set of collective variables can significantly accelerate the exploration of phase space, albeit at the price of distorting the distribution of microstates. Efficient re-weighting to recover the unbiased distribution can be nontrivial when employing adaptive sampling techniques such as Metadynamics, Variationally Enhanced Sampling or Parallel Bias Metadynamics, in which the system evolves in a quasi-equilibrium manner under a time-dependent bias. We introduce an iterative unbiasing scheme that makes efficient use of all the trajectory data, and that does not require the distribution to be evaluated on a grid. The method can thus be used even when the bias has a high dimensionality. We benchmark this approach against some of the existing schemes, on models systems with different complexities and dimensionalities.