The ongoing debate with regards to the origin of peak splittings in the lineshape of the molecule, as well as the strong dependence of the splitting on solvent interactions, make CV an ideal target for a detailed study using theoretical/computational spectroscopy predictions. Firstprinciples computational predictions of optical properties of solvated dyes can serve as an excellent tool to provide insight and interpretation of experimental data. However, to obtain truly predictive simulated spectra in real condensed phase systems, it is often necessary to explicitly describe the interaction of the chromophore with its complex environment, which may demand large-scale electronic structure calculations. Very accurate calculations of a complicated system are resource-intensive, thus a good trade-off between the accuracy and computational cost is essential. To accurately model absorption lineshapes of molecules embedded in complex environments, vibronic couplings, direct solute-solvent interactions, temperature effects and environmental polarization effects all have to be accounted for. In recent years, significant efforts have been dedicated to developing robust and computationally affordable approaches to capture these effects from first principles. The most widely-used method, the Franck-Condon (FC) approach, captures vibronic effects by approximating the ground and excited state potential energy surfaces (PESs) of the chromophore as harmonic around their respective minima. In addition to the approximations inherent in this method, the FC method often neglects direct solvent interactions, since the calculations of ground-and excited state vibrational frequencies require optimized molecular geometries that are only straightforwardly accessible when representing the solvent environment through polarizable continuum models. The ensemble method, although capable of capturing anharmonicity and solvent effects due to a molecular dynamics (MD) sampling of solute-solvent configurations, treats nuclear degrees of freedom purely classically and thus cannot describe vibronic effects. A number of hybrid approaches combining features of both methods have been introduced recently, but they often require an arbitrary timescale separation between slow solvent and fast chromophore motions that might not be applicable in systems like CV that undergo collective chromophore-environment motion.