def expected(self, observable, beta=1e5):
        """Wrapper to the expected_value function to fix the eigenbasis"""
        return expected_value(observable,
                         self.eig_energies,
                         self.eig_states,
                         beta)