EC version now reproduces non-interacting results! Horray!
TODO:
    Optimize - no part should scale as NFCI*environment - factorize them
    Complex to add broadening to get spectra
    More perturbations - non neutral excitations
    -omega term?
    IC version working
    2D version
    What to do with the correlation potential for full hubbard model?