Force-field parametrization based on radial and energy distribution functions

Implementation and Validation of an OpenMM Plugin for the Deep Potential Representation of Potential Energy

Machine learning potentials, particularly the deep potential (DP) model, have revolutionized molecular dynamics (MD) simulations, striking a balance between accuracy and computational efficiency. To facilitate the DP model's integration with the popular MD engine OpenMM, we have developed a versatile OpenMM plugin. This plugin supports a range of applications, from conventional MD simulations to alchemical free energy calculations and hybrid DP/MM simulations. Our extensive validation tests encompassed energy conservation in microcanonical ensemble simulations, fidelity in canonical ensemble generation, and the evaluation of the structural, transport, and thermodynamic properties of bulk water. The introduction of this plugin is expected to significantly expand the application scope of DP models within the MD simulation community, representing a major advancement in the field.

"Dividing and Conquering" and "Caching" in Molecular Modeling

Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, most important methodological advancements in more than half century of molecular modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but results are not transferable. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science. This framework is based on a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.

Force-field parametrization based on radial and energy distribution functions

We propose a novel force‐field‐parametrization procedure that fits the parameters of potential functions in a manner that the pair distribution function (DF) of molecules derived from candidate parameters can reproduce the given target DF. Conventionally, approaches to minimize the difference between the candidate and target DFs employ radial DFs (RDF). RDF itself has been reported to be insufficient for uniquely identifying the parameters of a molecule. To overcome the weakness, we introduce energy DF (EDF) as a target DF, which describes the distribution of the pairwise energy of molecules. We found that the EDF responds more sensitively to a small perturbation in the pairwise potential parameters and provides better fitting accuracy compared to that of RDF. These findings provide valuable insights into a wide range of coarse graining methods, which determine parameters using information obtained from a higher‐level calculation than that of the developed force field. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.