Obtaining an accurate theoretical model for the activation of dinitrogen by three-coordinate molybdenum amide complexes (e.g. Mo(NH2)3) is difficult due to the interaction of various high- and low-spin open-shell complexes along the reaction coordinate which must be treated with comparable levels of accuracy in order to obtain reasonable potential energy surfaces. Density functional theory with present-day functionals is a popular choice in this situation; however, the dinitrogen activation reaction energetics vary substantially with the choice of functional. An assessment of the reaction using specialized wave function based methods indicates that although current density functionals in general agree qualitatively on the mechanistic details of the reaction, a variety of high-level electron correlation methods (including CCSD(T), OD(T), CCSD(2), KS-CCSD(T), and spin-flip CCSD) provide a consistent but slightly different representation of the system.