Fast quasi-centroid molecular dynamics for water and ice
arXiv:2310.12246 · doi:10.1021/acs.jpcb.3c05028
Abstract
We describe how the fast quasi-centroid molecular dynamics (f-QCMD) method can be applied to condensed phase systems by approximating the quasi-centroid potential of mean force as a sum of inter- and intra-molecular corrections to the classical interaction potential. The corrections are found by using a regularised iterative Boltzmann inversion procedure to recover the inter- and intra-molecular quasi-centroid distribution functions obtained from a path integral molecular dynamics simulation. The resulting methodology is found to give good agreement with a previously published QCMD dipole absorption spectrum for liquid water, and satisfactory agreement for ice. It also gives good agreement with spectra from a recent implementation of CMD that uses a pre-computed elevated temperature potential of mean force. Modern centroid molecular dynamics methods therefore appear to be reaching a consensus regarding the impact of nuclear quantum effects on the vibrational spectra of water and ice.
10 pages, 4 figures, plus supplementary information and toc figure
References in corpus (7)
- Competing quantum effects in the dynamics of a flexible water model
- How to remove the spurious resonances from ring polymer molecular dynamics
- Boltzmann-conserving classical dynamics in quantum time-correlation functions: Matsubara dynamics
- On the Consistency of Approximate Quantum Dynamics Simulation Methods for Vibrational Spectra in the Condensed Phase
- Quantum dynamics using path integral coarse-graining
- Well-posedness of the Iterative Boltzmann Inversion
- Improved torque estimator for condensed-phase quasicentroid molecular dynamics