Frozen-density embedding theory (FDET) provides the formal framework for multilevel numerical simulations, such that a selected subsystem is described at the quantum mechanical level, whereas its environment is described by means of the electron density (frozen density; ρB( r →) ) The frozen density ρB( r →) is usually obtained from some lower-level quantum mechanical methods applied to the environment, but FDET is not limited to such choices for ρB( r →). The present work concerns the application of FDET, in which ρB( r →) is the statistically averaged electron density of the solvent <ρB( r →)> . The specific solute–solvent interactions are represented in a statistical manner in <ρB( r →)>. A full self-consistent treatment of solvated chromophore, thus involves a single geometry of the chromophore in a given state and the corresponding <ρB( r →)>. We show that the coupling between the two descriptors might be made in an approximate manner that is applicable for both absorption and emission. The proposed protocol leads to accurate (error in the range of 0.05 eV) descriptions of the solvatochromic shifts in both absorption and emission.