BH-DFTB/DFT calculations for iron clusters