Non-unitary time evolution via the Chebyshev expansion method

The Chebyshev expansion method is a well-established technique for computing the time evolution of quantum states, particularly in Hermitian systems with a bounded spectrum. Here, we show that the applicability of the Chebyshev expansion method extends well beyond this constraint: It remains valid across the entire complex plane and is thus suitable for arbitrary non-Hermitian matrices. We identify numerical rounding errors as the primary source of errors encountered when applying the method outside the conventional spectral bounds, and they are not caused by fundamental limitations. By carefully selecting the spectral radius and the time step, we show how these errors can be effectively suppressed, enabling accurate time evolution calculations in non-Hermitian systems. We derive an analytic upper bound for the rounding error, which serves as a practical guideline for selecting time steps in numerical simulations. As an application, we illustrate the performance of the method by computing the time evolution of wave packets in the Hatano-Nelson model.