We present a quantum-mechanical theory to study excitation energy transfers between molecular systems in solution. The model is developed within the time-dependent (TD) densityfunctional theory and the solvent effects are introduced in terms of the polarizable continuum model (PCM). Unique characteristic of this model is that both "reaction field" and screening effects are included in a coherent and self-consistent way. This is obtained by introducing proper solvent-specific operators in the Kohn-Sham equations and in the corresponding TD scheme. The solvation model exploits the integral equation formalism (IEF) version of PCM and it defines the solvent operators on a molecular cavity modeled on the real three-dimensional (3D) structure of the solute systems. Applications to EET in dimers of ethylene and naphtalene are presented and discussed.