We compute energies of intermediates using B3LYP hybrid density functional theory to predict the ground spin state of various intermediates along the catalytic cycle. We then vary the amount of Hartree-Fock exchange (aHF) in the functional to understand the dependence of energetics on the functional. For each intermediate, we compute the difference of energies between various spin states, and thereby identify the ground spin state.