Many-body potentials, in either the Embedded Atom or Finnis-Sinclair form, have gained wide popularity recently. The major difficulty in implementing the method concerns the derivation of suitable forms for the pair potential, electron density, and embedding, function which reproduce a range of empirically observable parameters such as elastic constants, defect formation energies, and defect Green's functions. This is a particularly difficult problem for niobium, which shows a variety of anomalous features in its phonon dispersion. Embedding functions which match only elastic constants may do a poor job of reproducing short wavelength behavior, and hence provide poor defect modeling. A straightforward method of deriving embedding functions for homonuclear BCC and
FCC metals will be presented which provides excellent agreement with experimental phonon dispersion curves and elastic constants, as well as Griineisen coefficients, vacancy formation energies, lattice constants and heats of sublimation. The results of the application of a set of many-body potentials derived in this fashion to nitrogen irradiated niobium will be presented.