AN EXAMINATION OF A DENSITY-FUNCTIONAL MOLECULAR MECHANICAL COUPLED POTENTIAL

in this article we describe the coupling of a density functional (DF) Hamiltonian with the molecular mechanics (MM) potential function AMBER. We examine a series of test cases in which we compare the binding energies and geometries of the complexes predicted by this coupled potential with those predicted by other theoretical methods and experiment to establish the relative accuracy of the DF/MM coupled potential. We find that the DF/MM coupled potential performs well in most cases studied and, in general, outperforms the semiempirical/MM approach. The interaction energies and structures obtained using this method appear to be insensitive to the use of nonlocal (NL) corrections to the DF method. The is fortuitous because the NL treatment is significantly more computationally expensive than the local treatment. However, NL corrections may be required to predict accurately the shape of potential energy surfaces that involve bond breaking and formation. The DF/MM method has also been applied to the determination of the solvation free energy for a series of ions using free-energy perturbation methods. The results obtained are good and can be improved by a simple scaling of the Lennard-Jones parameters for the ion in question. (C) 1995 by John Wiley and Sons, Inc.