We use density functional calculations to model the FeFe-cofactor of the iron-only nitrogenase. We determine the ground state geometical and magnetic structure and show that our results are in agreement with experimental observations. We model the reactivity of the FeFeco towards hydrogen and