Abstract:Collective variable-based enhanced sampling methods are routinely used on systems with metastable states, where high free energy barriers impede proper sampling of the free energy landscapes when using conventional molecular dynamics simulations. One such method is variationally enhanced sampling (VES), which is based on a variational principle where a bias potential in the space of some chosen slow degrees of freedom, or collective variables, is constructed by minimizing a convex functional. In practice, the bias potential is taken as a linear expansion in some basis function set. So far, primarily basis functions delocalized in the collective variable space, like plane waves, Chebyshev, or Legendre polynomials, have been used. However, there has not been an extensive study of how the convergence behavior is affected by the choice of the basis functions. In particular, it remains an open question if localized basis functions might perform better. In this work, we implement, tune, and validate Daubechies wavelets as basis functions for VES\@. The wavelets construct orthogonal and localized bases that exhibit an attractive multiresolution property. We evaluate the performance of wavelet and other basis functions on various systems, going from model potentials to the calcium carbonate association process in water. We observe that wavelets exhibit excellent performance and much more robust convergence behavior than all other basis functions, as well as better performance than metadynamics. In particular, using wavelet bases yields far smaller fluctuations of the bias potential within individual runs and smaller differences between independent runs. Based on our overall results, we can recommend wavelets as basis functions for VES.

On the usefulness of Meyer wavelets for deconvolution and density estimation

Uniform almost sure convergence and asymptotic distribution of the wavelet-based estimators of partial derivatives of multivariate density function under weak dependence

Nonlinear Wavelet Density Estimation with Data Missing at Random when Covariates Are Present

Wavelet methods in statistics: Some recent developments and their applications

Convergence Rate of Wavelet Density Estimator with Data Missing Randomly when Covariables Are Present

Minimax adaptive wavelet estimator for the simultaneous blind deconvolution with fractional Gaussian noise

Bayesian nonparametric regression using complex wavelets

Debiasing Welch's Method for Spectral Density Estimation

A practical guide to the recovery of wavelet coefficients from Fourier measurements

On the empirical spectral distribution of large wavelet random matrices based on mixed-Gaussian fractional measurements in moderately high dimensions

On high-dimensional wavelet eigenanalysis

Short-time homomorphic wavelet estimation

Multiplicative deconvolution under unknown error distribution

Nonparametric Density Estimation Using Wavelet Transformation And Scale-Space Zero-Crossing Reconstruction

On posterior distribution of Bayesian wavelet thresholding

Statistical inference for wavelet curve estimators of symmetric positive definite matrices

A Note on Block-Thresholded Wavelet Estimators with Correlated Noise

Shannon wavelet approximations of linear differential operators

Intrinsic wavelet regression for surfaces of Hermitian positive definite matrices

Improving the Efficiency of Variationally Enhanced Sampling with Wavelet-Based Bias Potentials

Adapting to Unknown Smoothness by Aggregation of Thresholded Wavelet Estimators