We apply the frozen density embedding method, using a full relaxation of embedded densities through a freeze-and-thaw procedure, to study the electronic structure of several benchmark ground-state charge-transfer complexes, in order to assess the merits and limitations of the approach for this class of systems. The calculations are performed using both semilocal and hybrid exchange-correlation (XC) functionals. The results show that embedding calculations using semilocal XC functionals yield rather large deviations with respect to the corresponding supermolecular calculations. Due to a large error cancellation effect, however, they can often provide a relatively good description of the electronic structure of charge-transfer complexes, in contrast to supermolecular calculations performed at the same level of theory. On the contrary, when hybrid XC functionals are employed, both embedding and supermolecular …