Jump to main content
Jump to site search

Issue 2, 2016
Previous Article Next Article

Modeling quantum nuclei with perturbed path integral molecular dynamics

Author affiliations

Abstract

The quantum nature of nuclear motions plays a vital role in the structure, stability, and thermodynamics of molecules and materials. The standard approach to model nuclear quantum fluctuations in chemical and biological systems is to use path-integral molecular dynamics. Unfortunately, conventional path-integral simulations can have an exceedingly large computational cost due to the need to employ an excessive number of coupled classical subsystems (beads) for quantitative accuracy. Here, we combine perturbation theory with the Feynman–Kac imaginary-time path integral approach to quantum mechanics and derive an improved non-empirical partition function and estimators to calculate converged quantum observables. Our perturbed path-integral (PPI) method requires the same ingredients as the conventional approach, but increases the accuracy and efficiency of path integral simulations by an order of magnitude. Results are presented for the thermodynamics of fundamental model systems, an empirical water model containing 256 water molecules within periodic boundary conditions, and ab initio simulations of nitrogen and benzene molecules. For all of these examples, PPI simulations with 4 to 8 classical beads recover the nuclear quantum contribution to the total energy and heat capacity at room temperature within a 3% accuracy, paving the way toward seamless modeling of nuclear quantum effects in realistic molecules and materials.

Graphical abstract: Modeling quantum nuclei with perturbed path integral molecular dynamics

Back to tab navigation

Supplementary files

Publication details

The article was received on 12 Sep 2015, accepted on 28 Oct 2015 and first published on 30 Oct 2015


Article type: Edge Article
DOI: 10.1039/C5SC03443D
Author version
available:
Download author version (PDF)
Citation: Chem. Sci., 2016,7, 1368-1372
  • Open access: Creative Commons BY-NC license
  •   Request permissions

    Modeling quantum nuclei with perturbed path integral molecular dynamics

    I. Poltavsky and A. Tkatchenko, Chem. Sci., 2016, 7, 1368
    DOI: 10.1039/C5SC03443D

    This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. Material from this article can be used in other publications provided that the correct acknowledgement is given with the reproduced material and it is not used for commercial purposes.

    Reproduced material should be attributed as follows:

    • For reproduction of material from NJC:
      [Original citation] - Published by The Royal Society of Chemistry (RSC) on behalf of the Centre National de la Recherche Scientifique (CNRS) and the RSC.
    • For reproduction of material from PCCP:
      [Original citation] - Published by the PCCP Owner Societies.
    • For reproduction of material from PPS:
      [Original citation] - Published by The Royal Society of Chemistry (RSC) on behalf of the European Society for Photobiology, the European Photochemistry Association, and RSC.
    • For reproduction of material from all other RSC journals:
      [Original citation] - Published by The Royal Society of Chemistry.

    Information about reproducing material from RSC articles with different licences is available on our Permission Requests page.

Search articles by author

Spotlight

Advertisements