  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.
