This enables using the Chebyshev expansion for the exponential. The span is computed for the discrete Hamiltonian with factors coming from the FD discretization of the gradient.