If n is greater than 0, one can interpret the corresponding P +n as the density matrix for the (2n + 1) -anion of the two chromophoresolvent complex, with molecular orbitals taken from DFT calculations for the singly anionic complex. The inclusion of low-lying virtual orbitals of the complex in the construction of the generalized density matrix P +n allows us to not only identify bath orbitals that entangle with the chromophore occupied orbitals but also more precisely capture those entangled with the low-lying virtual orbitals of the chromophore, both of which are expected to be essential for an accurate description of the excitation/de-excitation of a solvated chromophore.