On the Consistency of Approximate Quantum Dynamics Simulation Methods for Vibrational Spectra in the Condensed Phase

Abstract

Including quantum mechanical effects on the dynamics of nuclei in the condensed phase is challenging, because the complexity of exact methods grows exponentially with the number of quantum degrees of freedom. Efforts to circumvent these limitations can be traced down to two approaches: methods that treat a small subset of the degrees of freedom with rigorous quantum mechanics, considering the rest of the system as a static or classical environment, and methods that treat the whole system quantum mechanically, but using approximate dynamics. Here we perform a systematic comparison between these two philosophies for the description of quantum effects in vibrational spectroscopy, taking the Embedded Local Monomer (LMon) model and a mixed quantum-classical (MQC) model as representatives of the first family of methods, and centroid molecular dynamics (CMD) and thermostatted ring polymer molecular dynamics (TRPMD) as examples of the latter. We use as benchmarks D${}_2$O doped with HOD and pure H${}_2$O at three distinct thermodynamic state points (ice Ih at 150K, and the liquid at 300K and 600K), modeled with the simple q-TIP4P/F potential energy and dipole moment surfaces. With few exceptions the different techniques yield IR absorption frequencies that are consistent with one another within a few tens of cm${}^{-1}$. Comparison with classical molecular dynamics demonstrates the importance of nuclear quantum effects up to the highest temperature, and a detailed discussion of the discrepancies between the various methods let us draw some (circumstantial) conclusions about the impact of the very different approximations that underlie them. Such cross validation between radically different approaches could indicate a way forward to further improve the state of the art in simulations of condensed-phase quantum dynamics.

Bundled together with the Supporting Materials

References in corpus (6)

• Canonical sampling through velocity-rescaling
• Competing quantum effects in the dynamics of a flexible water model
• Efficient stochastic thermostatting of path integral molecular dynamics
• How to remove the spurious resonances from ring polymer molecular dynamics
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• Derivation of a true (t->0+) quantum transition-state theory. Part I: Uniqueness and equivalence to ring-polymer molecular dynamics transition-state-theory

Cited by in corpus (26)

• Nuclear quantum effects enter the mainstream
• Raman Spectrum and Polarizability of Liquid Water from Deep Neural Networks
• Fine Tuning Classical and Quantum Molecular Dynamics using a Generalized Langevin Equation
• 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
• High Order Path Integrals Made Easy
• Stability of Complex Biomolecular Structures: Vander Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects
• Path-integral dynamics of water using curvilinear centroids
• An analysis of nonadiabatic ring-polymer molecular dynamics and its application to vibronic spectra
• Inexpensive modelling of quantum dynamics using path integral generalized Langevin equation thermostats
• All-Electron, Real-Space Perturbation Theory for Homogeneous Electric Fields: Theory, Implementation, and Application within DFT
• Decisive role of nuclear quantum effects on surface mediated water dissociation at finite temperature
• Quantum Ring-Polymer Contraction Method: Including nuclear quantum effects at no additional computational cost in comparison to ab-initio molecular dynamics
• Quantum dynamics using path integral coarse-graining
• Fully First-Principles Surface Spectroscopy with Machine Learning
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• Mean-field Matsubara dynamics: Analysis of path-integral curvature effects in rovibrational spectra
• Should Thermostatted Ring Polymer Molecular Dynamics be used to calculate thermal reaction rates?
• How to obtain thermostatted ring polymer molecular dynamics from exact quantum dynamics and when to use it
• Discretized hierarchical equations of motion in mixed Liouville--Wigner space for two-dimensional vibrational spectroscopies of liquid water
• Fast quasi-centroid molecular dynamics for water and ice
• Improved torque estimator for condensed-phase quasicentroid molecular dynamics
• Ab initio study of water dissociation on a charged Pd(111) surface
• Dynamical Heterogeneity in Supercooled Water and its Spectroscopic Fingerprints
• Representing spherical tensors with scalar-based machine-learning models
• Nuclear Quantum Effects in Multi-Step Condensed Matter Chemistry: A Path Integral Molecular Dynamics Study of Thermal Decomposition
• Quantitative Molecular Simulations