Infrared photodissociation is used to record the vibrational spectrum of FeO2+(He)2–4 which shows three bands at 1035, 980, and 506 cm−1. Quantum chemical multi-reference configuration interaction calculations (MRCISD) of structures and harmonic frequencies show that these bands are due to two different isomers, an inserted dioxo complex with Fe in the +V oxidation state and a side-on superoxo complex with Fe in the +II oxidation state. These two are separated by a substantial barrier, 53 kJ/mol, whereas the third isomer, an end-on complex between Fe+ and an O2 molecule, is easily converted into the side-on complex. For all three isomers, states of different spin multiplicity have been considered. Our best energies are computed at the MRCISD+Q level, including corrections for complete active space and basis set extension, core-valence correlation, relativistic effects, and zero-point vibrational energy. The average coupled pair functional (ACPF) yields very similar energies. Density functional theory (DFT) differs significantly from our best estimates for this system, with the TPSS functional yielding the best results. The other functionals tested are BP86, PBE, B3LYP, TPSSh, and B2PLYP. Complete active space second order perturbation theory, (CASPT2) performs better than DFT, but less good than ACPF.