Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: From Force Fields to Machine Learning Potentials

Abstract

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, and 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 thermostatting. 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 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 help to shed new light on the quantum nature of important interactions in biological matter.

Deciphering the Influence of Ground-State Distributions on the Calculation of Photolysis Observables

“…QT can properly treat both high and low-frequency modes and performs well even for (moderately) anharmonic systems. , We note that a different implementation of a similar idea was proposed by Dammak et al and termed quantum thermal bath (QTB). QTB was applied to studying vibrational spectra or the structure of liquid water. , A recent extension of QTB was devised to tackle the ZPE leakage issue. , Our group has recently compared the impact of QT and Wigner sampling on the prediction of photoabsorption cross-sections of several exemplary VOCs within the NEA . QT was found superior whenever low-frequency anharmonic modes play a role in the photochemistry/photophysics of a molecule.…”

Scalable hybrid deep neural networks/polarizable potentials biomolecular simulations including long-range effects

“…Such performances should continue to improve thanks to further Deep-HP optimizations, TorchANI updates and GPUs hardware evolutions. Deep-HP will enable the implementation of the next generation of improved MLPs [84][85][86] and has been designed to be a place for their further development. It will include direct neural networks coupling with physics-driven contributions going beyond multipolar electrostatics and polarization through inclusion of many-body dispersion models.…”

Force-field-enhanced neural network interactions: from local equivariant embedding to atom-in-molecule properties and long-range effects

“…A promising strategy that could provide a solid middle-ground is the ne-tuning approach where a general model is slightly re-trained with limited data to better represent a specic target system. Further work will focus on the development and use of more general datasets (such as the SPICE dataset 97 ), the extension of the FENNIX approach and its inclusion in the optimized Deep-HP platform 98,99 present within the Tinker-HP package. 12,100 Appendix: dissociations across chemical space…”