  An iterative approach is introduced, which allows an efficient solution of
the hierarchical equations of motion (HEOM) for the steady state of open
quantum systems. The approach combines the method of matrix equations with an
efficient preconditioning technique to reduce the numerical effort for solving
the HEOM. Illustrative applications to simulate nonequilibrium charge transport
in single-molecule junctions demonstrate the performance of the method.
