Computing accurate absorption spectra of riboflavin (RBF) has proven a difficult task for computational chemistry. Time-dependent density functional theory have herein been employed using a wide range of recent range-separated and hybrid meta functionals to investigate ultraviolet and visible spectra of RBF to determine if any progress has been made through recent developments. It is concluded that B3LYP and PBE0 perform the best throughout the entire test set. However, since all methods deviate from experimental results by at least 40 nm when computing the spectra in vacuum, two approaches to describe aqueous solution are implemented together with the MPWB1K, B3LYP, and PBE0 functionals: implicitly using integral equation formulation of the polarized continuum model (minor improvement) and explicitly through molecular dynamics (MD) simulations of the molecule embedded in a water cluster whereafter snapshots of RBF-water clusters are extracted and time-dependent density functional theory calculations performed. The resulting averaged spectra from the MD-simulated clusters show a constant blue-shift for all peaks by similar to 20 nm compared to experimental data at the TD-B3LYP/6-31+G(d,p) level.