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

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.

arXiv commentpublished in PCCP

References in corpus (29)

• Van der Waals Density Functional for General Geometries
• DeePMD-kit: A deep learning package for many-body potential energy representation and molecular dynamics
• Perspective: How good is DFT for water?
• An exchange functional that tests the robustness of the plasmon description of the van der Waals density functional
• Competing quantum effects in the dynamics of a flexible water model
• Ab initio theory and modeling of water
• Efficient stochastic thermostatting of path integral molecular dynamics
• How to remove the spurious resonances from ring polymer molecular dynamics
• Density, structure and dynamics of water: the effect of Van der Waals interactions
• Anomalous Nuclear Quantum Effects in Ice
• Unraveling quantum mechanical effects in water using isotopic fractionation
• The Many-Body Expansion Combined with Neural Networks
• Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances
• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• Effect of quantum nuclear motion on hydrogen bonding
• Water: one molecule, two surfaces, one mistake
• Direct Measurement of Competing Quantum Effects on the Kinetic Energy of Heavy Water upon Melting
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Classical-path integral adaptive resolution in molecular simulation: towards a smooth quantum-classical coupling
• Room temperature compressibility and diffusivity of liquid water from first principles
• Quantum Ring-Polymer Contraction Method: Including nuclear quantum effects at no additional computational cost in comparison to ab-initio molecular dynamics
• Simulation of quantum zero-point effects in water using a frequency-dependent thermostat
• Nuclear quantum effects in liquid water from path-integral simulations using an ab initio force matching approach
• Efficient mixed-force first-principles molecular dynamics
• Multiple Time Step Integrators in Ab Initio Molecular Dynamics
• Optimization of an exchange-correlation density functional for water
• Electronic and Nuclear Quantum Effects on the Ice XI/Ice Ih Phase Transition
• Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface
• From Classical to Quantum and Back: Hamiltonian Adaptive Resolution Path Integral, Ring Polymer, and Centroid Molecular Dynamics