Towards an accurate modeling of optical rotation for solvated systems: anharmonic vibrational contributions coupled to the polarizable continuum model