Adiabatic terminator for fermionic hierarchical equations of motion^†

The hierarchical equation of motion method has become one of the most popular numerical methods for describing the dissipative dynamics of open quantum systems linearly coupled to environment. However, its applications to systems with strong electron correlation are largely restrained by the computational cost, which is mainly caused by the high truncation tier L required to accurately characterize the strong correlation effect. In this work, we develop an adiabatic terminator by decoupling the principal dissipation mode with the fastest dissipation rate from the slower ones. The adiabatic terminator leads to substantially enhanced convergence with respect to L as demonstrated by the numerical tests carried out on a single impurity Anderson model. Moreover, the adiabatic terminator alleviates the numerical instability problems in the long-time dissipative dynamics.