Jump to main content
Jump to site search
Access to RSC content Close the message box

Continue to access RSC content when you are not at your institution. Follow our step-by-step guide.


Issue 1, 2019
Previous Article Next Article

Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost

Author affiliations

Abstract

It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge–Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call “monomer PIMD”, correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.

Graphical abstract: Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost

Back to tab navigation

Article information


Submitted
27 Sep 2018
Accepted
30 Nov 2018
First published
03 Dec 2018

Phys. Chem. Chem. Phys., 2019,21, 409-417
Article type
Paper
Author version available

Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost

D. C. Elton, M. Fritz and M. Fernández-Serra, Phys. Chem. Chem. Phys., 2019, 21, 409
DOI: 10.1039/C8CP06077K

Social activity

Search articles by author

Spotlight

Advertisements