4.6 Article

AN ADAPTIVE PARAREAL ALGORITHM: APPLICATION TO THE SIMULATION OF MOLECULAR DYNAMICS TRAJECTORIES

期刊

SIAM JOURNAL ON SCIENTIFIC COMPUTING
卷 44, 期 1, 页码 B146-B176

出版社

SIAM PUBLICATIONS
DOI: 10.1137/21M1412979

关键词

parallel-in-time simulation; molecular dynamics; adaptive algorithm

资金

  1. ANR [ANR-15-CE23-0019-06]
  2. Alexander von Humboldt foundation
  3. DFG [EXC-2046/1-390685689]
  4. European Research Council (ERC) under the European Union [810367]
  5. European High-Performance Computing Joint Undertaking (JU) [955701]
  6. European Union's Horizon 2020 research and innovation programme

向作者/读者索取更多资源

The aim of this article is to design parareal algorithms for thermostated molecular dynamics simulations. The traditional parareal algorithm is not suitable for molecular dynamics due to its limitations. This article proposes a modified version of the parareal algorithm that is better suited for molecular dynamics simulations. However, the modified algorithm still has some limitations, including intermediate trajectory blow-up, encounters with undefined values, and no computational advantage in long time horizons. Through numerical experiments, this article demonstrates that the adaptive algorithm overcomes the limitations of the standard algorithm and achieves significant improvements.
The aim of this article is to design parareal algorithms in the context of thermostated molecular dynamics. In its original setup, the fine and coarse propagators used in the parareal algorithm solve the same dynamics with different time-steps, with the goal of achieving accuracy in the limit of small time-step of the integrators involved. This is typically not useful in molecular dynamics, where one is interested in extremely long trajectories and where the time-step of the fine propagator is in practice chosen as large as possible, that is, close to the limit of stability of the numerical scheme. In this article, we consider a version of the parareal algorithm which is better suited to molecular dynamics simulations, and wherein the propagators involved use the same time step while employing different potential energy landscapes to drive the dynamics. Although the parareal algorithm always converges, it suffers from various limitations in this context: intermediate blow-up of the trajectory (which makes it impossible to postprocess) may be observed; in certain cases the trajectory encounters undefined values before converging (the way the algorithm handles them might depend on the computer architecture); the algorithm does not provide any computational gain in terms of wall-clock time (compared to a standard sequential integration) in the limit of increasingly long time horizons. We highlight these issues with numerical experiments and provide some elements of theoretical analysis. We then present a modified version of the parareal algorithm wherein the algorithm adaptively divides the entire time horizon into smaller time-slabs where the aforementioned issues are circumvented. Using numerical experiments on toy examples, we show that the adaptive algorithm overcomes the various limitations of the standard parareal algorithm, thereby allowing for significantly improved gains.

作者

我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。

评论

主要评分

4.6
评分不足

次要评分

新颖性
-
重要性
-
科学严谨性
-
评价这篇论文

推荐

暂无数据
暂无数据