Periodic Reporting for period 2 - BeStMo (Beyond Static Molecules: Modeling Quantum Fluctuations in Complex Molecular Environments)
Reporting period: 2018-09-01 to 2020-02-29
Numerous fascinating and useful dynamical phenomena in complex molecular systems stem from an interplay between quantum-mechanical electronic and nuclear fluctuations. Many of the anomalies of water – the most ubiquitous liquid on Earth – are the result of a subtle coupling between its peculiar hydrogen bonded network and quantized nuclear fluctuations. Another particularly striking example is the aspirin molecular crystal, in which the prevalent form I polymorph is stabilized by entropy stemming from an unexpected dynamical coupling between non-covalent van der Waals interactions and quantized lattice vibrations. Additional prominent examples of coaction of electronic and nuclear quantum fluctuations include the peculiar ability of protons to permeate atomically-thin membranes at room temperature, the pervasive observation of delocalized boson peaks in disordered molecular materials, and non-trivial wavelike behavior of large molecules in electric and optical fields. Our ability to atomistically model and understand all of these important phenomena requires substantial breakthroughs.
Our final goal is to bridge the accuracy of quantum mechanics with the efficiency of force fields, enabling large-scale predictive quantum molecular dynamics simulations for complex systems containing 1000s of atoms, and leading to novel conceptual insights into quantum-mechanical fluctuations in large molecular systems. The project goes well beyond the presently possible applications and once successful will pave the road towards having a suite of first-principles-based modeling tools for a wide range of realistic materials, such as biomolecules, nanostructures, disordered solids, and organic/inorganic interfaces.
In WP1, we have developed a coarse-grained multiscale approach to van der Waals interactions in arbitrary external fields including electronic and nuclear fluctuations. Using this model, we have obtained surprising results on modification of van der Waals interactions in systems subjected to external environments. These results have been reported in multidisciplinary journals Science Advances, Physical Review Letters, and Nature Communications (3 publications). In the next period, we will release two open source code package and two comprehensive publications that will allow others to do calculations on complex systems.
In WP2, we have extended the previously published MBD method to treat systems with more than 10,000 atoms. We applied these developments to study the quantum mechanics of fully solvated proteins. The corresponding manuscript has been published in Science Advances. At least three more drafts are being prepared for publication in high-impact journals. At the moment we are also implementing a coarse-grained version of the MBD model that will push its applicability to systems with 100,000 atoms.
In WP3, we have developed a suite of machine learning models for essentially exact molecular dynamics of small molecules. The three ensuing papers have been published in Science Advances (1 paper) and Nature Communications (2 papers). The community has received this work very enthusiastically and I would consider these developments as true breakthroughs. At the moment, we are able to run molecular dynamics simulations with fully quantized electrons and nuclei for molecules with 50-100 atoms for the first time. The remaining challenge is to treat long-range non-covalent interactions with physical approaches such as MBD and derive a coupling potential between machine-learned force field and physics-based Hamiltonians. We already obtained preliminary results indicating that this approach is promising.
In WP4, in collaboration with the group of Michele Ceriotti at EPFL, we have developed a method for carrying out efficient high-order path integral simulations of molecular systems. This work has been published in the Journal of Chemical Theory and Computation and implemented in publically available i-PI code.
Overall, I feel that the original plan of the proposal has been fully accomplished and, indeed, exceeded in WP2 and WP3. During the next period, we will concentrate on the mentioned developments mainly in WP2 and WP3.