Abstract
Exceptional points (EPs) are singularities in the spectra of non-Hermitian operators where eigenvalues and eigenvectors coalesce. Open quantum systems have recently been explored as EP testbeds due to their non-Hermitian nature. However, most studies focus on the Markovian limit, leaving a gap in understanding EPs in the non-Markovian regime. This work addresses this gap by proposing a general framework based on two numerically exact descriptions of non-Markovian dynamics: the pseudomode equation of motion (PMEOM) and the hierarchical equations of motion (HEOM). The PMEOM is particularly useful due to its Lindblad-type structure, aligning with previous studies in the Markovian regime while offering deeper insights into EP identification. This framework incorporates non-Markovian effects through auxiliary degrees of freedom, enabling the discovery of additional or higher-order EPs that are inaccessible in the Markovian regime. We demonstrate the utility of this approach using the spin-boson model and linear bosonic systems.