Photocatalytic water splitting has attracted considerable attention for renewable energy production. Since the first reported photocatalytic water splitting by titanium dioxide, this material remains one of the most promising photocatalysts, due to its suitable band gap and band-edge positions. However, predicting both of these properties is a challenging task for existing computational methods. Here we show how Koopmans spectral functionals can accurately predict the band structure and level alignment of rutile, anatase, and brookite TiO₂ using a computationally efficient workflow that only requires (a) a DFT calculation of the photocatalyst/vacuum interface and (b) a Koopmans spectral functional calculation of the bulk photocatalyst. The success of this approach for TiO₂ suggests that this strategy could be deployed for assessing the suitability of novel photocatalyst candidates. This record contains for each polymorph: (a) a calculation of a slab; (b) KI and pKIPZ band structure calculation of the bulk system.