分子动力学
统计物理学
加速度
采样(信号处理)
马尔可夫链
玻尔兹曼分布
亚稳态
蒙特卡罗方法
计算机科学
Atom(片上系统)
物理
经典力学
数学
量子力学
机器学习
统计
滤波器(信号处理)
计算机视觉
嵌入式系统
作者
Leon Klein,Andrew Y. K. Foong,Tor Erlend Fjelde,Bruno Mlodozeniec,Marc Brockschmidt,Sebastian Nowozin,Frank Noé,Ryota Tomioka
出处
期刊:Cornell University - arXiv
日期:2023-02-03
被引量:20
标识
DOI:10.48550/arxiv.2302.01170
摘要
Molecular dynamics (MD) simulation is a widely used technique to simulate molecular systems, most commonly at the all-atom resolution where equations of motion are integrated with timesteps on the order of femtoseconds ($1\textrm{fs}=10^{-15}\textrm{s}$). MD is often used to compute equilibrium properties, which requires sampling from an equilibrium distribution such as the Boltzmann distribution. However, many important processes, such as binding and folding, occur over timescales of milliseconds or beyond, and cannot be efficiently sampled with conventional MD. Furthermore, new MD simulations need to be performed for each molecular system studied. We present Timewarp, an enhanced sampling method which uses a normalising flow as a proposal distribution in a Markov chain Monte Carlo method targeting the Boltzmann distribution. The flow is trained offline on MD trajectories and learns to make large steps in time, simulating the molecular dynamics of $10^{5} - 10^{6}\:\textrm{fs}$. Crucially, Timewarp is transferable between molecular systems: once trained, we show that it generalises to unseen small peptides (2-4 amino acids) at all-atom resolution, exploring their metastable states and providing wall-clock acceleration of sampling compared to standard MD. Our method constitutes an important step towards general, transferable algorithms for accelerating MD.
科研通智能强力驱动
Strongly Powered by AbleSci AI