We present a computational study on the spectroscopic properties of UV-Vis absorbing dyes in water solution. We model the solvation environment by using both continuum and discrete models, with and without polarization, to establish how the physical and chemical properties of the solute-solvent interaction may affect the spectroscopic response of aqueous systems. Seven different compounds were chosen, representing different classes of organic molecules. The classical atomistic description of the solvent molecules was enriched with polarization effects treated by means of the fluctuating charges (FQ) model, propagated to the first-order response function of the quantum-mechanical (QM) solute to include its effects withing the modeling of the electronic excitations of the systems. Results obtained with the QM/FQ model were compared with those from continuum solvation models as well as nonpolarizable atomistic models, and then confronted with the experimental values to determine the accuracy that can be expected with each level of theory. Moreover, a thorough structural analysis using molecular dynamics simulations is provided for each system. K E Y W O R D S excitation energies, QM/FQ, QM/MM, QM/PCM, solvent effects, TD-DFT 1 | I N TR O DU C TI O N One-photon absorption spectroscopy within the UV-Visible range is often the most direct and inexpensive analytical tool that can be used to study the electronic properties of a system. Most commonly, such measurements are carried out on solvated samples, with water being a ubiquitous choice.With the gradual increase in the complexity of the systems under investigation, the correct interpretation of experimental data is increasingly reliant on their calculated ab initio counterparts. Many theoretical models based on quantum mechanics (QM), accompanied by their computational implementations, have been presented over the years offering different levels of compromise between the computational cost and the accuracy of the results. [1][2][3] At present, methods based on density functional theory (DFT) and its time-dependent counterpart (TD-DFT) have become the most popular choice for the simulation of absorption spectra of medium-large organic molecular systems thanks to their versatility stemming from the freedom of choice of density functional and basis set, as well as the favorable scaling with system size which allows their application to increasingly large systems. [1,[4][5][6] Many benchmarks studies have been presented elaborating on the merits and limitations of TD-DFT for the simulation of UV-Vis spectroscopy, as well as on the most appropriate choice of functional and basis set combination for different types of system. [6][7][8][9][10][11][12][13][14][15][16] And though many computational studies are carried out on isolated systems, solvent effects should not be neglected for the presence of the solvation environment can significantly alter the electronic absorption properties of a system, both qualitatively and quantitatively. [17][18][19][20][21][22][23][24][25][2...