On the practical truncation tier of fermionic hierarchical equations of motion.

The fermionic hierarchical equations of motion (HEOM) approach has found wide application in the exploration of open quantum systems, and extensive efforts have been committed to improving its efficiency and accuracy in practical calculations. In this work, by scrutinizing the stationary-state and dynamic properties of Kondo-correlated quantum impurity systems, we show that the strength of Kondo correlation induced by the system-environment entanglement primarily determines the converged hierarchical truncation tier of the HEOM method. This complements the rule of thumb regarding the positive correlation between the height of hierarchy and system-environment coupling strength. These insights will provide useful guidelines for developing a more sophisticated fermionic HEOM method for the investigation of many-body open quantum systems.