As part of a multi-scale approach for modeling hydrogen embrittlement in hardened steels we have investigated, employing density functional theory methods, the stability and concentrations of the point defect clusters present in metastable Fe-C-H alloys with vacancies. The defect spectrum is found to be strongly dependent on the local vacancy concentration, and for low hydrogen levels sharp hig...