Search

We report the implementation of a multi-CPU and multi-GPU massively parallel platform dedicated to the explicit inclusion of nuclear quantum effects (NQEs) in the Tinker-HP molecular dynamics (MD) package. The platform, denoted Quantum-HP, exploits two simulation strategies: the Ring-Polymer Molecular Dynamics (RPMD) that provides exact structural properties at the cost of a MD simulation in an extended space of multiple replicas, and the adaptive Quantum Thermal Bath (adQTB) that imposes the quantum distribution of energy on a classical system via a generalized Langevin thermostat and provides computationally affordable and accurate (though approximate) NQEs. We discuss some implementation details, efficient numerical schemes, parallelization strategies and quickly review the GPU acceleration of our code. Our implementation allows an efficient inclusion of NQEs in MD simulations for very large systems, as demonstrated by scaling tests on water boxes with more than 200,000 atoms (simulated using the AMOEBA polarizable force field). We test the compatibility of the approach with Tinker-HP's recently introduced Deep-HP machine learning potentials module by computing water properties using the DeePMD potential with adQTB thermostating. Finally, we show that the platform is also compatible with the alchemical free energy estimation capabilities of Tinker-HP and fast enough to perform simulations. Therefore, we study how the NQEs affect the hydration free energy of small molecules solvated with the recently developed Q-AMOEBA water force field. Overall, the Quantum-HP platform allows users to perform routine quantum MD simulations of large condensed-phase systems and will participate to shed a new light on the quantum nature of important interactions in biological matter.

International audience; We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQE) 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 we 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 classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.

Increasing interest has risen for nuclear quantum effects (NQE) in the recent past. Indeed, NQE such as proton tunneling and zero point energy often play a crucial role in the properties of hydrogen-containing materials. The standard methods to simulate NQE are based on path integrals. An alternative to these methods is the Quantum Thermal Bath (QTB): it is based on a Langevin equation where the classical degrees of freedom are coupled to an ensemble of quantum harmonic oscillators. In the classical Langevin equation, the random force is a white noise and fulfills the classical fluctuation-dissipation theorem, while within the QTB formalism, it fulfills the quantum fluctuation-dissipation theorem. We investigate through simple models the reliability and the limits of the QTB and show that the QTB enables realistic simulations including NQE of condensed-phase systems, generating static and dynamic information such as pair correlation functions and vibrational spectra which can be confronted with experimental results. We show that the QTB is particularly successful in the study of the symmetrization of hydrogen bonds in several systems. Indeed, the difficulty lies in the identification of a precise transition pressure since this phase transition is often blurred by quantum or thermal fluctuations. In high-pressure ice, it depends on the oxygen-oxygen distance but it can be affected by ionic impurities and by the asymmetric environment of hydrogen bonds as in the delta phase of AlOOH. Moreover, by comparing results from QTB and standard ab initio simulations, we are able to disentangle the respective roles of NQE and thermal fluctuations in these phase transitions.; L’étude des effets quantiques nucléaires (NQE) suscite de plus en plus d’intérêt. En effet, les effets quantiques comme l’effet tunnel ou l’énergie de point zéro, peuvent profondément modifier les propriétés de matériaux constitués d'atomes légers comme l'hydrogène. Les méthodes standards de simulation des NQE sont basées sur les intégrales de chemin. Le bain thermique quantique (QTB) constitue une alternative à ces méthodes: le principe est que les degrés de liberté classiques du système obéissent à une équation de Langevin et sont couplés à des oscillateurs harmoniques quantiques. Dans l’équation de Langevin classique, la force aléatoire est un bruit blanc et le théorème de fluctuation-dissipation classique est vérifié; avec le QTB, le théorème de fluctuation-dissipation quantique est vérifié. Nous étudierons à travers des modèles simples la validité et les limites du QTB et montrerons qu'il permet de simuler des systèmes de la matière condensée en incluant les NQE en générant leurs propriétés structurales et dynamiques. Nous montrerons que le QTB est particulièrement adapté à l’étude de la symétrisation de liaisons hydrogènes et permet d'identifier précisément une pression de transition. Celle-ci dépend de la distance entre deux oxygènes voisins comme dans la glace sous haute pression, mais est modifiée par la présence d'impuretés ioniques ou par l'environnement atomique des liaisons hydrogènes comme dans la phase delta de AlOOH. De plus, en comparant des simulations classiques à des simulations QTB, nous pouvons identifier les rôles respectifs des effets quantiques et thermiques dans ces transitions de phase.

International audience; We derive the stochastic master equations which describe the evolution of open quantum systems in contact with a heat bath and undergoing indirect measurements. These equations are obtained as a limit of a quantum repeated measurement model where we consider a small system in contact with an infinite chain at positive temperature. At zero temperature it is well-known that one obtains stochastic differential equations of jump-diffusion type. At strictly positive temperature, we show that only pure diffusion type equations are relevant.