Maybe I should weigh in on this, since I've done a number of these calcs (yay, I'm even cited in Neese's review).
As cgk says, the only way to get the value theoretically is by explicit calculation, I'm afraid, and DFT would be the most likely candidate. While DFT calcs on, say, simple organic molecules and reactions can be done today without a lot of in-depth knowledge, calculating open-shell metalloorganic systems take a bit of knowing what you're doing.
To begin with there's the simple question of what functional to use. As cgk mentions, the relative energies of your spin states depend on the exchange energy, and so they depend almost entirely on which functional you choose. Especially with hybrid functionals where you might choose different values for the HF-exchange component, and the improvement on one aspect doesn't necessarily translate to other ones. There's no consensus on what's best. (According to Per Siegbahn, B3LYP* with 15% HF-exchange gives better energies for transition-metal atoms. According to Martin Kaupp, BHLYP, with 50% HF-exchange gives better spin-densities. But it also increases the risk of spin contamination, which is always a problem.) On top of all that, there's a purely practical issue in that most programs have convergence problems with that kind of systems, and usually require explicit starting guesses, vshifting etc to 'find' the correct wave function. It can be determined experimentally though, using EPR spectroscopy.