We have conducted a test of running NRLMOL with different basis sets to comapre against Gaussian for a series of atoms, differences wher found, here are the energies using basis STO-3G.
We looked at Vanadium to account for the difference.
We tried refining the mesh, it did not work.
We compared the two basis sets used (the were the same).
So the difference in the calculation lies elsewhere.