Physics > Chemical Physics
[Submitted on 3 Nov 2019 (this version), latest version 15 Mar 2020 (v3)]
Title:Dimension-free path-integral molecular dynamics without preconditioning
View PDFAbstract:Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics (PIMD) - including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods - is that for a given MD timestep, the overlap between the exact ring-polymer Boltzmann distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte-Carlo/MD sampling schemes, and it also causes the well-known divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these problems can be avoided through the introduction of "dimension-free" numerical integration schemes for which the sampled ring-polymer distribution has non-zero overlap with the exact distribution in the infinite-bead limit. Moreover, we show that this can be achieved by using a (previously introduced) strongly stable method for the free-ring-polymer evolution in combination with a (newly introduced) mollification of the forces from the external physical potential. The resulting dimension-free numerical integration schemes yield finite error bounds for a given MD timestep, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic cases. Importantly, dimension-free RPMD achieves these benefits while preserving strong stability, symplecticity, time reversibility, and global second-order accuracy; and it remains a simple, black-box method by avoiding computational costs, tunable parameters, or system-specific implementations.
Submission history
From: Roman Korol [view email][v1] Sun, 3 Nov 2019 17:14:02 UTC (81 KB)
[v2] Fri, 17 Jan 2020 19:49:55 UTC (1,136 KB)
[v3] Sun, 15 Mar 2020 04:54:35 UTC (1,118 KB)
Current browse context:
physics.chem-ph
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.