While E solute is computed as in vacuum, the wavefunction parameters are determined in the presence of the solvent's reaction field. Half of E int accounts for the work required to polarize the dielectric medium when a polarized solute molecule is inserted into the cavity, making the total energy a free energy. Thus, the free energy is obtained with consistent densities for both the solute and solvent.