extend Hubbard model implementation