We present benchmark calculations on excitation energies based on time-dependent density functional theory (TDDFT) as well as orbital relaxed self-consistent and constricted variational DFT (RSCF-CV-DFT) with and without use of the Tamm-Dancoff approximation. The compilation contains results for the 3d complexes MnO₄⁻, CrO₄²⁻, and VO₄³⁻, as well as the 4d congeners RuO₄, TcO₄⁻, and MoO₄²⁻, and 5d homologues OsO₄, ReO₄⁻, and WO₄²⁻. Considerations have been given to the local density approximation (LDA) and the functionals BP86 and PBE based on the generalized gradient approximation (GGA), as well as the hybrids B3LYP, BHLYP, and PBE0 and the length corrected functional LCBP86. We find for the 3d complexes that RSCF-CV-DFT fares better than TDDFT. Thus, in the case of RSCF-CV-DFT, the average root-mean-square deviations (RMSDs) are 0.25-0.3 eV for GGAs, 0.1 eV for B3LYP, and 0.45 eV for BHLYP. TDDFT affords RMSDs that on average range from 0.3 eV for local functionals to 0.7 eV for BHLYP with the largest fraction of Hartree-Fock (HF) exchange. TDDFT is seen to fare better among the heavier tetraoxo systems. For the 4d and 5d systems, the three functionals B3LYP, PBE0 with an intermediate fraction of HF exchange, and LCBP86 have the lowest RMSD of 0.2 eV, whereas the local functionals (LDA, BP86, BPE) and BHLYP with the highest HF fraction and LCBP86* have a somewhat larger RMSD of 0.3 eV. Nearly the same performance is observed for RSCF-CV-DFT with respect to the different functionals in the case of the 4d and 5d systems. Thus, for the heavier tetraoxo systems, the two DFT schemes are comparable in accuracy.