Abstract:We introduce a novel approach for decomposing and learning every scale of a given multiscale objective function in $\mathbb{R}^d$, where $d\ge 1$. This approach leverages a recently demonstrated implicit bias of the optimization method of gradient descent by Kong and Tao, which enables the automatic generation of data that nearly follow Gibbs distribution with an effective potential at any desired scale. One application of this automated effective potential modeling is to construct reduced-order models. For instance, a deterministic surrogate Hamiltonian model can be developed to substantially soften the stiffness that bottlenecks the simulation, while maintaining the accuracy of phase portraits at the scale of interest. Similarly, a stochastic surrogate model can be constructed at a desired scale, such that both its equilibrium and out-of-equilibrium behaviors (characterized by auto-correlation function and mean path) align with those of a damped mechanical system with the original multiscale function being its potential. The robustness and efficiency of our proposed approach in multi-dimensional scenarios have been demonstrated through a series of numerical experiments. A by-product of our development is a method for anisotropic noise estimation and calibration. More precisely, Langevin model of stochastic mechanical systems may not have isotropic noise in practice, and we provide a systematic algorithm to quantify its covariance matrix without directly measuring the noise. In this case, the system may not admit closed form expression of its invariant distribution either, but with this tool, we can design friction matrix appropriately to calibrate the system so that its invariant distribution has a closed form expression of Gibbs.

Hyper Molecular Dynamics with A Local Bias Potential

Local bias potential in hyper molecular dynamics method

Ridge-based Bias Potentials to Accelerate Molecular Dynamics.

Molecular hyperdynamics coupled with the nonorthogonal tight-binding approach: Implementation and validation

Application of the Accelerated Molecular Dynamics Simulations to the Folding of A Small Protein

Hyperdynamics Simulations with Ab Initio Forces

Free energy calculations: An efficient adaptive biasing potential method

Enhanced Sampling of Multidimensional Free-Energy Landscapes Using Adaptive Biasing Forces

A Random-Batch Monte Carlo Method for Many-Body Systems with Singular Kernels.

Refining Potential Energy Surface through Dynamical Properties via Differentiable Molecular Simulation

Generalization of Biasing Particle Insertion and Deletion for Grand Canonical Monte Carlo Simulation

Efficient molecular dynamics simulations with many-body potentials on graphics processing units

A Method for Molecular Dynamics on Curved Surfaces

Constrained Broyden Dimer Method with Bias Potential for Exploring Potential Energy Surface of Multistep Reaction Process.

A local resampling trick for focused molecular dynamics

Accelerated Bayesian Inference for Molecular Simulations using Local Gaussian Process Surrogate Models

A fast molecular dynamics method with short-range interaction

Open-boundary Hamiltonian adaptive resolution. From grand canonical to non-equilibrium molecular dynamics simulations

Automated construction of effective potential via algorithmic implicit bias

Estimating Free-Energy Surfaces and Their Convergence from Multiple, Independent Static and History-Dependent Biased Molecular-Dynamics Simulations with Mean Force Integration