A systematic theoretical study of the work function behavior for Si1−xGexheterostructures over the whole composition range, from Si (x = 0) to Ge (x = 1), is presented. Our results, obtained through Density Functional Theory calculations and in good agreement with experimental evidences, show that increasing the Ge content lowers the work function value. We find that in order to exactly reproduce this behaviour in relation to the work function of pure Ge and Si systems and their concentrations, a deviation from the linear Vegard's rule is necessary. However, the calculated bowing parameter is very small, thus making the simple linear interpolation a valid approximation to obtain the work function of complex SiGe alloys.