The hierarchical equations of motion are comprised of a number of terms, each of which operates on an order of auxiliary density operator (ADO). These auxiliaries are evolved through a series of sub-equations and define contributions to the full density matrix. They are structured as an infinite hierarchy and their positions within that hierarchy dictates their contribution to the density matrix. The HEOM take the following form,