Molecular mechanics 2 force field accuracy range

Equation (4-5) can be directly utilized in statistical mechanical Monte Carlo and molecular dynamics simulations by choosing an appropriate QM model, balancing computational efficiency and accuracy, and MM force fields for biomacromolecules and the solvent water. Our group has extensively explored various QM/MM methods using different quantum models, ranging from semiempirical methods to ab initio molecular orbital and valence bond theories to density functional theory, applied to a wide range of applications in chemistry and biology. Some of these studies have been discussed before and they are not emphasized in this article. We focus on developments that have not been often discussed. 

Twenty years ago Car and Parrinello introduced an efficient method to perform Molecular Dynamics simulation for classical nuclei with forces computed on the fly by a Density Functional Theory (DFT) based electronic calculation [1], Because the method allowed study of the statistical mechanics of classical nuclei with many-body electronic interactions, it opened the way for the use of simulation methods for realistic systems with an accuracy well beyond the limits of available effective force fields. In the last twenty years, the number of applications of the Car-Parrinello ab-initio molecular d3mam-ics has ranged from simple covalent bonded solids, to high pressure physics, material science and biological systems. There have also been extensions of the original algorithm to simulate systems at constant temperature and constant pressure [2], finite temperature effects for the electrons [3], and quantum nuclei [4]. 

---