DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Reproducing Quantum Probability Distributions at the Speed of Classical Dynamics: A New Approach for Developing Force-Field Functors

Abstract

Modeling nuclear quantum effects is required for accurate molecular dynamics (MD) simulations of molecules. The community has paid special attention to water and other biomolecules that show hydrogen bonding. Standard methods of modeling nuclear quantum effects like Ring Polymer Molecular Dynamics (RPMD) are computationally costlier than running classical trajectories. A force-field functor (FFF) is an alternative method that computes an effective force field that replicates quantum properties of the original force field. In this work, we propose an efficient method of computing FFF using the Wigner–Kirkwood expansion. As a test case, we calculate a range of thermodynamic properties of Neon, obtaining the same level of accuracy as RPMD, but with the shorter runtime of classical simulations. Lastly, by modifying existing MD programs, the proposed method could be used in the future to increase the efficiency and accuracy of MD simulations involving water and proteins.

Authors:
 [1];  [1]; ORCiD logo [2]
  1. Harvard Univ., Cambridge, MA (United States). Dept. of Chemistry and Chemical Biology
  2. Harvard Univ., Cambridge, MA (United States). Dept. of Chemistry and Chemical Biology; Canadian Inst. for Advanced Research, Toronto, ON (Canada)
Publication Date:
Research Org.:
Energy Frontier Research Centers (EFRC) (United States). Center for Excitonics (CE)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1470507
Grant/Contract Number:  
SC0001088
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Physical Chemistry Letters
Additional Journal Information:
Journal Volume: 9; Journal Issue: 7; Related Information: CE partners with Massachusetts Institute of Technology (lead); Brookhaven National Laboratory; Harvard University; Journal ID: ISSN 1948-7185
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; solar (photovoltaic); solid state lighting; photosynthesis (natural and artificial); charge transport; optics; synthesis (novel materials); synthesis (self-assembly); synthesis (scalable processing)

Citation Formats

Sundar, Vikram, Gelbwaser-Klimovsky, David, and Aspuru-Guzik, Alán. Reproducing Quantum Probability Distributions at the Speed of Classical Dynamics: A New Approach for Developing Force-Field Functors. United States: N. p., 2018. Web. doi:10.1021/acs.jpclett.7b03254.
Sundar, Vikram, Gelbwaser-Klimovsky, David, & Aspuru-Guzik, Alán. Reproducing Quantum Probability Distributions at the Speed of Classical Dynamics: A New Approach for Developing Force-Field Functors. United States. https://doi.org/10.1021/acs.jpclett.7b03254
Sundar, Vikram, Gelbwaser-Klimovsky, David, and Aspuru-Guzik, Alán. Mon . "Reproducing Quantum Probability Distributions at the Speed of Classical Dynamics: A New Approach for Developing Force-Field Functors". United States. https://doi.org/10.1021/acs.jpclett.7b03254. https://www.osti.gov/servlets/purl/1470507.
@article{osti_1470507,
title = {Reproducing Quantum Probability Distributions at the Speed of Classical Dynamics: A New Approach for Developing Force-Field Functors},
author = {Sundar, Vikram and Gelbwaser-Klimovsky, David and Aspuru-Guzik, Alán},
abstractNote = {Modeling nuclear quantum effects is required for accurate molecular dynamics (MD) simulations of molecules. The community has paid special attention to water and other biomolecules that show hydrogen bonding. Standard methods of modeling nuclear quantum effects like Ring Polymer Molecular Dynamics (RPMD) are computationally costlier than running classical trajectories. A force-field functor (FFF) is an alternative method that computes an effective force field that replicates quantum properties of the original force field. In this work, we propose an efficient method of computing FFF using the Wigner–Kirkwood expansion. As a test case, we calculate a range of thermodynamic properties of Neon, obtaining the same level of accuracy as RPMD, but with the shorter runtime of classical simulations. Lastly, by modifying existing MD programs, the proposed method could be used in the future to increase the efficiency and accuracy of MD simulations involving water and proteins.},
doi = {10.1021/acs.jpclett.7b03254},
journal = {Journal of Physical Chemistry Letters},
number = 7,
volume = 9,
place = {United States},
year = {2018},
month = {3}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 1 work
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Force-field functor theory: classical force-fields which reproduce equilibrium quantum distributions
journal, January 2013


Neutron diffraction of liquid neon and xenon along the coexistence line
journal, March 1992

  • Bellissent-Funel, M. C.; Buontempo, U.; Filabozzi, A.
  • Physical Review B, Vol. 45, Issue 9
  • DOI: 10.1103/physrevb.45.4605

Ring-Polymer Molecular Dynamics: Quantum Effects in Chemical Dynamics from Classical Trajectories in an Extended Phase Space
journal, April 2013


Asymptotic expansions and effective potentials for almost classicalN-body systems
journal, October 1989


An introduction to biomolecular simulations and docking
journal, August 2014


Amber 2017, University of California, San Francisco
dataset, January 2017


Nuclear Quantum Effects and Enzyme Dynamics in Dihydrofolate Reductase Catalysis
journal, March 2002

  • Agarwal, Pratul K.; Billeter, Salomon R.; Hammes-Schiffer, Sharon
  • The Journal of Physical Chemistry B, Vol. 106, Issue 12
  • DOI: 10.1021/jp020190v

Wigner-Kirkwood expansion: Calculation of ‘‘almost classical’’ static properties of a Lennard-Jones many-body system
journal, September 1987


Improving the convergence of closed and open path integral molecular dynamics via higher order Trotter factorization schemes
journal, August 2011

  • Pérez, Alejandro; Tuckerman, Mark E.
  • The Journal of Chemical Physics, Vol. 135, Issue 6
  • DOI: 10.1063/1.3609120

Biomolecular Simulation: A Computational Microscope for Molecular Biology
journal, June 2012


An introduction to biomolecular simulations and docking
journal, August 2014


Force fields and molecular dynamics simulations
journal, January 2011

  • González, M. A.
  • École thématique de la Société Française de la Neutronique, Vol. 12
  • DOI: 10.1051/sfn/201112009

Polarizable force field development for lipids and their efficient applications in membrane proteins
journal, April 2017

  • Chu, Huiying; Cao, Liaoran; Peng, Xiangda
  • WIREs Computational Molecular Science, Vol. 7, Issue 5
  • DOI: 10.1002/wcms.1312

Nuclear Quantum Effects in Water and Aqueous Systems: Experiment, Theory, and Current Challenges
journal, April 2016


Nuclear Quantum Effects and Hydrogen Bonding in Liquids
journal, July 2003

  • Raugei, Simone; Klein, Michael L.
  • Journal of the American Chemical Society, Vol. 125, Issue 30
  • DOI: 10.1021/ja0351995

Nuclear Quantum Effects and Enzyme Dynamics in Dihydrofolate Reductase Catalysis
journal, March 2002

  • Agarwal, Pratul K.; Billeter, Salomon R.; Hammes-Schiffer, Sharon
  • The Journal of Physical Chemistry B, Vol. 106, Issue 12
  • DOI: 10.1021/jp020190v

Nuclear quantum effects on an enzyme-catalyzed reaction with reaction path potential: Proton transfer in triosephosphate isomerase
journal, March 2006

  • Wang, Mingliang; Lu, Zhenyu; Yang, Weitao
  • The Journal of Chemical Physics, Vol. 124, Issue 12
  • DOI: 10.1063/1.2181145

Exploiting the isomorphism between quantum theory and classical statistical mechanics of polyatomic fluids
journal, April 1981

  • Chandler, David; Wolynes, Peter G.
  • The Journal of Chemical Physics, Vol. 74, Issue 7
  • DOI: 10.1063/1.441588

Path-Integral Calculation of the Two-Particle Slater Sum for He 4
journal, March 1966


A fast path integral method for polarizable force fields
journal, January 2009

  • Fanourgakis, George S.; Markland, Thomas E.; Manolopoulos, David E.
  • The Journal of Chemical Physics, Vol. 131, Issue 9
  • DOI: 10.1063/1.3216520

Ring-Polymer Molecular Dynamics: Quantum Effects in Chemical Dynamics from Classical Trajectories in an Extended Phase Space
journal, April 2013


Centroid molecular dynamics: A quantum dynamics method suitable for the parallel computer
journal, July 2000


Accelerating the convergence of path integral dynamics with a generalized Langevin equation
journal, February 2011

  • Ceriotti, Michele; Manolopoulos, David E.; Parrinello, Michele
  • The Journal of Chemical Physics, Vol. 134, Issue 8
  • DOI: 10.1063/1.3556661

Applications of higher order composite factorization schemes in imaginary time path integral simulations
journal, November 2001

  • Jang, Seogjoo; Jang, Soonmin; Voth, Gregory A.
  • The Journal of Chemical Physics, Vol. 115, Issue 17
  • DOI: 10.1063/1.1410117

Improving the convergence of closed and open path integral molecular dynamics via higher order Trotter factorization schemes
journal, August 2011

  • Pérez, Alejandro; Tuckerman, Mark E.
  • The Journal of Chemical Physics, Vol. 135, Issue 6
  • DOI: 10.1063/1.3609120

Forward–backward semiclassical dynamics for quantum fluids using pair propagators: Application to liquid para -hydrogen
journal, October 2003

  • Nakayama, Akira; Makri, Nancy
  • The Journal of Chemical Physics, Vol. 119, Issue 16
  • DOI: 10.1063/1.1611473

Practical evaluation of condensed phase quantum correlation functions: A Feynman–Kleinert variational linearized path integral method
journal, December 2003

  • Poulsen, Jens Aage; Nyman, Gunnar; Rossky, Peter J.
  • The Journal of Chemical Physics, Vol. 119, Issue 23
  • DOI: 10.1063/1.1626631

An accurate and simple quantum model for liquid water
journal, November 2006

  • Paesani, Francesco; Zhang, Wei; Case, David A.
  • The Journal of Chemical Physics, Vol. 125, Issue 18
  • DOI: 10.1063/1.2386157

Quantum dynamical effects in liquid water: A semiclassical study on the diffusion and the infrared absorption spectrum
journal, October 2009

  • Liu, Jian; Miller, William H.; Paesani, Francesco
  • The Journal of Chemical Physics, Vol. 131, Issue 16
  • DOI: 10.1063/1.3254372

An efficient ring polymer contraction scheme for imaginary time path integral simulations
journal, July 2008

  • Markland, Thomas E.; Manolopoulos, David E.
  • The Journal of Chemical Physics, Vol. 129, Issue 2
  • DOI: 10.1063/1.2953308

Force-field functor theory: classical force-fields which reproduce equilibrium quantum distributions
journal, January 2013


On the Quantum Correction For Thermodynamic Equilibrium
journal, June 1932


Quantum Statistics of Almost Classical Assemblies
journal, July 1933


Wigner function approach to the quantum Brownian motion of a particle in a potential
journal, January 2007

  • Coffey, W. T.; Kalmykov, Yu. P.; Titov, S. V.
  • Physical Chemistry Chemical Physics, Vol. 9, Issue 26
  • DOI: 10.1039/b614554j

Theory and application of the quantum phase-space distribution functions
journal, August 1995


Time dependence of a general class of quantum distribution functions
journal, November 1984


CHARMM: The biomolecular simulation program
journal, July 2009

  • Brooks, B. R.; Brooks, C. L.; Mackerell, A. D.
  • Journal of Computational Chemistry, Vol. 30, Issue 10
  • DOI: 10.1002/jcc.21287

Nuclear quantum effects in liquid water from path-integral simulations using an ab initio force-matching approach
journal, November 2014


Liquid argon: Monte carlo and molecular dynamics calculations
journal, January 1971


Molecular-dynamics calculation of the quantum corrections to the pair-distribution function for a Lennard-Jones fluid
journal, March 1984


Molecular dynamics study of the structural and thermodynamic properties of MgO crystal with quantum correction
journal, July 1989

  • Matsui, Masanori
  • The Journal of Chemical Physics, Vol. 91, Issue 1
  • DOI: 10.1063/1.457484

Quantum-corrected pair distribution function of liquid neon
journal, June 1985


Many-body and quantum effects in the radial distribution function of liquid neon and argon
journal, November 1995


Structure Study of Liquid Neon by Neutron Diffraction
journal, November 1971

  • de Graaf, L. A.; Mozer, B.
  • The Journal of Chemical Physics, Vol. 55, Issue 10
  • DOI: 10.1063/1.1675607

The equation of state of neon between 27 and 70 K
journal, August 1969


OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation
journal, November 2012

  • Eastman, Peter; Friedrichs, Mark S.; Chodera, John D.
  • Journal of Chemical Theory and Computation, Vol. 9, Issue 1
  • DOI: 10.1021/ct300857j

Direct determination of phase coexistence properties of fluids by Monte Carlo simulation in a new ensemble
journal, July 1987


Phase equilibria by simulation in the Gibbs ensemble: Alternative derivation, generalization and application to mixture and membrane equilibria
journal, March 1988


Influence of simulation protocols on the efficiency of Gibbs ensemble Monte Carlo simulations
journal, December 2013

  • Cortés Morales, Angel D.; Economou, Ioannis G.; Peters, Cornelis J.
  • Molecular Simulation, Vol. 39, Issue 14-15
  • DOI: 10.1080/08927022.2013.828209

MCCCS Towhee: a tool for Monte Carlo molecular simulation
journal, December 2013


Accurate and precise determination of critical properties from Gibbs ensemble Monte Carlo simulations
journal, September 2015

  • Dinpajooh, Mohammadhasan; Bai, Peng; Allan, Douglas A.
  • The Journal of Chemical Physics, Vol. 143, Issue 11
  • DOI: 10.1063/1.4930848