Ab initio atom-atom potentials using CamCASP: Many-body potentials for the pyridine dimer
In Part I of this two-part investigation we described a methodology for the development of robust, analytic, many-body atom-atom potentials for small organic molecules from first principles and demonstrated how the CamCASP program can be used to derive the damped, distributed multipole models for pyridine. Here we demonstrate how the theoretical ideas for the short-range models described in Part I, which are implemented in the CamCASP suite of programs, can be used to develop a series of many-body potentials for the pyridine system. Even the simplest of these potentials exhibit r.m.s. errors of only about 0.6kJ mol-1 for the low-energy pyridine dimers, significantly surpassing the best empirical potentials. Our best model is shown to support eight stable minima, four of which have not been reported in the literature before. Further, the functional form can be made systematically more elaborate so as to improve the accuracy without a significant increase in the human-time spent in their generation. We investigate the effects of anisotropy, rank of multipoles, and choice of polarizability and dispersion models.