Towards an accurate modeling of optical rotation for solvated systems: anharmonic vibrational contributions coupled to continuum solvation models