We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQEs) in molecular simulations: the adaptive quantum thermal bath (adQTB). In this method, zero-point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation−dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to that of classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.

Quantum thermal bath (QTB) simulations reproduce statistical nuclear quantum effects via a Langevin equation with a colored random force. Although this approach has proven efficient for a variety of chemical and condensed-matter problems, the QTB, as many other semiclassical methods, suffers from zero-point energy leakage (ZPEL). The absence of a reliable criterion to quantify the ZPEL without resorting to demanding comparisons with path integral based calculations has so far hindered the use of the QTB for the simulation of real systems. In this work, we establish a quantitative connection between ZPEL in the QTB framework and deviations from the quantum fluctuation-dissipation theorem (FDT) that can be monitored along the simulation. This provides a rigorous general criterion to detect and quantify the ZPEL without any a priori knowledge of the system under study. We then use this criterion to build an adaptive QTB method that strictly enforces the quantum FDT at all frequencies via an on-the-fly, spectrally resolved fine tuning of the system-bath coupling coefficients. The validity of the adaptive approach is first demonstrated on a simple two-oscillator model. It is then applied to two more realistic problems: the description of the vibrational properties of a model aluminium crystal at low temperature and the simulation of the liquid-solid phase transition in a 13-atom neon cluster. In both systems, the standard QTB results are strongly altered by the ZPEL, which can be essentially eliminated using the adaptive approach.

We introduce a new parametrization of the AMOEBA polarizable force field for water denoted Q-AMOEBA, for use in simulations that explicitly account for nuclear quantum effects (NQEs). This study is made possible thanks to the recently introduced adaptive Quantum Thermal Bath (adQTB) simulation technique which computational cost is comparable to classical molecular dynamics. The flexible Q-AMOEBA model conserves the initial AMOEBA functional form, with an intermolecular potential including an atomic multipole description of electrostatic interactions (up to quadrupole), a polarization contribution based on the Thole interaction model and a buffered 14–7 potential to model van der Waals interactions. It has been obtained by using a ForceBalance fitting strategy including high-level quantum chemistry reference energies and selected condensed-phase properties targets. The final Q-AMOEBA model is shown to accurately reproduce both gas-phase and condensed-phase properties, notably improving the original AMOEBA water model. This development allows the fine study of NQEs on water liquid phase properties such as the average H–O–H angle compared to its gas-phase equilibrium value, isotope effects, and so on. Q-AMOEBA also provides improved infrared spectroscopy prediction capabilities compared to AMOEBA03. Overall, we show that the impact of NQEs depends on the underlying model functional form and on the associated strength of hydrogen bonds. Since adQTB simulations can be performed at near classical computational cost using the Tinker-HP package, Q-AMOEBA can be extended to organic molecules, proteins, and nucleic acids opening the possibility for the large-scale study of the importance of NQEs in biophysics.