Gaussian accelerated molecular dynamics (GaMD) is an enhanced sampling method for MD simulations that can efficiently sample thermodynamic properties such as the free energy landscape of the system. When the system potential V (r), where r denotes the position vector of an N -atom system, is lower than a threshold energy E, GaMD fills the energy wells by adding a harmonic boost potential ∆V (r), i.e.,