Refractory transition metal nitrides exhibit a plethora of polymorphic expressions and chemical stoichiometries. To afford a better understanding of how defects may play a role in the structural and thermodynamics of these nitrides, using density-functional theory calculations, we investigate the influence of point and pair defects in bulk metastable γ-MoN and its (001) surface. We report favorable formation of Schottky defect pairs of neighboring Mo and N vacancies in bulk γ-MoN and apply this as a defect-mediated energy correction term to the surface energy of γ-MoN(001) within the ab initio atomistic thermodynamics approach. We also inspect the structural distortions in both bulk and surfaces of γ-MoN by using the partial radial distribution function, g(r), of Mo-N bond lengths. Large atomic displacements are found in both cases, leading to a broad spread of Mo-N bond length values when compared to their idealized bulk values. We propose that these structural and thermodynamic analyses may provide some insight into a better understanding of metastable materials and their surfaces.